Team:Freiburg/Modelling

• < home >
• < team >
  • → members
  • → attributions
  • → freiburg
• < project >
  • → description
  • → results
  • → modeling
  • → notebook
  • → sample data
• < human practices >
  • → oath
  • → ethics
  • → safety
• < gallery >
• < sponsors >
• < to the forums >

Modelling: Rational protein design

The Idea First there was the idea to replace the chemically produced Ni-NTA columns by some biological, reproducible material. It seemed unreasonable to produce an ion precipitating substance in a complicated, energy consuming way – completely artificial, while so many proteins in nature are capable of the same task – and much more complicated things. But how to harness that? The second milestone in the generation of the idea came up in a lecture of Prof. Martin, about Nickel allergy. He found a mutation in the TLR-4 receptor that introduces Histidines in the LRR motif of the protein, which are capable of binding nickel and thus forming complexes of receptors on the cell surface of immune cells, triggering the inflammatory reaction.

caption

This was the necessary clue that gave answer to the question. In the following work, a lot of investigation about the perfect Nickel precipitating protein was done. It seemed that there is no optimal natural protein which would serves our needs. Furthermore, many of the natural LRR proteins – of course – have their own biological function, we we were afraid could interfere in the cellular system we wanted to express them in. Luckily there are over a thousand already known structures of LRR proteins, what gave us a broad choice of motifs to choose and compare. The Construction of a new protein After a long and detailed search, we found out that it is most reasonable to use bacterial LRR motifs, since they seemed very well conserved in sequence, and they are the shortest – with only 21 amino acids per LRR repeat. (Wei 2008, Kajava 1998). We wanted to have the protein as simple and as predictable in behavior and structure as possible. Several search inquiries led to the most conserved and shortest of all bacterial LRR (PDB: NO), unluckily this protein is a toxin derived from Yersinia pestis. We of course did not want to mess around with toxins,

Caption

especially not with one of the black plague. We also did not know what amino acid pattern on the LRR motif would cause the toxic effect. There probably would have been no harm from this protein alone, even if we did clone and express it, but why should we do that when we have a huge choice for free at hand anyway. Finally we settled on the bacterial ligase PDB:3CVR. A ligase would be definitively harmless, especially since we planned to restructure it completely.

To do this we needed to understand how the structure of the Ligase was made. It was that our desired protein LRR motif was only a part of a bigger protein “factory” which included several domains that seemed to serve distinct functions – which we did not want. So we took only the LRR domain. It looked to us, as if this domain had mostly a structural role in this “factory”, but we could not tell that exactly. To avoid any unspecific biological function it was necessary to rearrange the amino acid composition. We fed the protein sequence in a free online consensus sequence generator called weblogo from the Berkeley University.

Caption

By analyzing the logo it was obvious which positions of the LRR motif were conserved and which not – AND: which of the non conserved amino acids appeared in what kind of patterns. Were there positions in the protein that required a polar amino acid? Or non polar, hydrophobic /-philic, charged, non-charged?. We compared the consensus sequence with the 3D structure, using PYMOL, to extract as much information as possible and then came up with this ideal consensus sequence:

Caption

The next consideration we had to do, was: how many Nickel do we need on the surface of our ideal Nickel binding protein, in what pattern, with what distances between, and at what angles towards each other to allow proper ion complexion? Nickel can be complexed by imidazole structures from 4 planar orthogonal directions, as well as to axial positions. It can however only take four ligands at once, preferably in a planar orientation. Cobalt has a bipyramidal setup for ligand-binding, too, but can take up to six ligands. The distances from a N-atom in the Imidazole ring to the ion had to be between 3 and 6 Angström. We found a crystal structure of a different protein (PDB:) which was resolved with three Histidines complexing a Nickel ion for a comparison, as well as some old publications that analyzed peptide ion bonds that taught us how the complex should look like. (Jordan 1974) We wanted our protein to be short. First, to cause as little unspecific binding as possible, second, to allow an easy expression and third to keep costs los for gene synthesis. For a proper purification of His-tagged proteins in Ni-NTA columns we knew that there are up to three Nickel involved in the interaction between the His tag (which consists of 6 or 7 Histidines) and the column. These have to be in close proximity to one another. There are always two Histidines of the tag binding to one ion.

This was what we had to implement into our LRR backbone.

To get a first impression how Histidines are positioned in a bacterial Ligase, we replaced all non conserved amino acids by Histidines and sent the sequence for a 3D sequence prediction to I-TASSER, an online structure prediction software tool.

Caption

Caption

We received this pdb file with a C-score= 0.53. Obviously this was not a realistic model, since the folding would surely be impaired by so many Histidines everywhere. But is was useful to find several good locations, were the Nickel could fit in in respect to the requirements it has for binding ligands. We manually fit in a Nickel ion, as it was crystallized in the mentioned PDB file, and measured the distances and evaluated the 3dimensional orientation of the Histidines towards the Nickel. Histidines can coordinate ligands with their free electron pair pointing planar away from the imidazole ring. What we further realized from the structure file was, that the end of the LRR segments were “open”, that means, the hydrophobic core of the protein was exposed and as it is visible in the prediction, curled in on one end into a sort of helix. This means the protein folding is not reliable and the structure needs some caps on both ends to stabilize the LRR core motif.

Caption

A solution to this was shown by Schmidt et al 2010, who crystallized the TLR-4 receptor. In this very nice piece of work he dissected the TLR4(PDB: 3FXI) into 3 parts, since it was to large and unhandy to be crystallized at once

Caption

To overcome this problem of an exposed hydrophobic core, he used the N- and C-terminal protein fragments of a LRR protein derived from hagfish. They tried a variety of different versions of how to glue together his fragments and these N- and C-terminal caps until he found a working one.

Caption

We used this knowledge for our purpose and took the same sequences and attached them to our protein sequence. These caps partially still show the typical LRR consensus sequences (which luckily is highly conserved in all kingdoms of life!), which made it possible to fit them onto our stack of LRR loops in the right position.

To the outer end the protein has a helix on the N-terminal side shielding off the core and the C-terminus the LRR disappears slowly turn by turn with more and more hydrophilic amino acids replacing the LRR consensus sequence.

After combining these caps with our structural analysis for possible Nickel binding sites and filling the empty positions of the consensus sequence with the ideal sequence derived from the ligase sequence(see above), we were still unsure whether this would work out. To increase our chances of success we designed six different versions with distinct Histidine patterns on the LRR core.

After evaluating the C-Scores of I-TASSER we saw that our favorite versions 1,2 and 4 gave acceptable scores. We then reverse translated them with a codon usage optimized for E. coli, using http://www.bioinformatics.org/sms2/rev_trans.html The sequence was handed over to our sponsor ATG:biosynthetics, who ran their own analysis on it, to optimize the sequence for RNA trafficking and secondary structures. In pictures below the alignment of our reverse translated sequence and the returned sequence is shown: ours above, ATG:biosynthetics sequence below. Underlined in red are the mutations the company introduced to optimize the expression. They also were so kind to synthesize the three genes for us as a gift.

Rational 1

ATGtgcccgagccgttgcagctgtagcggcaccgaaattcgctgcaacagcaaaggcctgaccagcgtgccgaccggcattccgagcagcgcgacccgcctggaactggaaag caacaaactgcaaagcctgccggatggcgtgtttgataaactgaccgcgctgcataaaagcaacaaccagctgaccagcctgccggataacctgccggcgagcctggaacatct ggcggtgagcaacaaccagctgaccagcctgccggataacctgccggcgagcctggaagcgctgcatgtgagcaacaaccagctgaccagcctgccggataacctgccggcgag cctggaacatctggcggtgagcaacaaccagctgaccagcctgccggataacctgccggcgagcctggaagcgctgcatgtgagcaacaaccagctgaccagcctgccggataa cctgccggcgagcctggaacatctggcggtgagcaacaaccagctgaccagcctgccggataacctgccggcgagcctgaaagcgctgcatctggataccaaccagctgaaaag cgtgccggatggcatttttgatcgcctgaccagcctgcaaaaaatttggctgcaaaccaacccgtgggattgcagctgcccgcgcattgattatctgagccgctggctgaacaa aaacagccagaaagaacagggcagcgcgaaatgcagcggcagcggcaaaccggtgcgcagcattatttgcccgTAGTAA

Rational 2

ATGtgcccgagccgttgcagctgtagcggcaccgaaattcgctgcaacagcaaaggcctgaccagcgtgccgaccggcattccgagcagcgcgacccgcctggaactggaaag caacaaactgcaaagcctgccggatggcgtgtttgataaactgacccagctgaccaaaagcaacaaccatctgcatagcctgccggataacctgccggcgagcctggaagtgctg gatgtgagcaacaaccatctgcatagcctgccggataacctgccggcgagcctggaagtgctggatgtgagcaacaaccatctgcatagcctgccggataacctgccggcgagcc tggaagtgctggatgtgagcaacaaccatctgcatagcctgccggataacctgccggcgagcctggaagtgctggatgtgagcaacaaccatctgcatagcctgccggataacctg ccggcgagcctggaagtgctggatgtgagcaacaaccatctgcatagcctgccggataacctgccggcgagcctgaaagaactggcgctggataccaaccatctgcatagcgtgcc ggatggcatttttgatcgcctgaccagcctgcaaaaaatttggctgcaaaccaacccgtgggattgcagctgcccgcgcattgattatctgagccgctggctgaacaaaaacagcc agaaagaacagggcagcgcgaaatgcagcggcagcggcaaaccggtgcgcagcattatttgcccgTAGTAA

Rational 4

ATGtgcccgagccgttgcagctgtagcggcaccgaaattcgctgcaacagcaaaggcctgaccagcgtgccgaccggcattccgagcagcgcgacccgcctggaactggaaagc aacaaactgcaaagcctgccggatggcgtgtttgataaactgacccagctgaccaaaagcaacaaccagctgaccagcctgccggatcatctgccgcatagcctggaagtgctgg atgtgagcaacaaccagctgaccagcctgccggatcatctgccgcatagcctggaagtgctggatgtgagcaacaaccagctgaccagcctgccggatcatctgccgcatagcctg gaagtgctggatgtgagcaacaaccagctgaccagcctgccggatcatctgccgcatagcctggaagtgctggatgtgagcaacaaccagctgaccagcctgccggatcatctgc cgcatagcctggaagtgctggatgtgagcaacaaccagctgaccagcctgccggatcatctgccgcatagcctgaaagaactggcgctggataccaaccagctgaaaagcgtgc cggatcatatttttcatcgcctgaccagcctgcaaaaaatttggctgcaaaccaacccgtgggattgcagctgcccgcgcattgattatctgagccgctggctgaacaaaaacag ccagaaagaacagggcagcgcgaaatgcagcggcagcggcaaaccggtgcgcagcattatttgcccgTAGTAA

Step 1

The Precipitator protein (P) binds the plastic column (C) by its plastic binding domain. In reality the parameter “column” is not a soluble reagent, but the surface of the serological pipette that offers binding sites for the precipitator. We omitted the calculation for a surface reaction, since the simplified model is sufficient for our needs.

Steady state:

700px