Orthogonal coiled-coils.html
From 2009.igem.org
Line 12: | Line 12: | ||
<br><br> | <br><br> | ||
Based on the published rules for the prediction of stability of coiled-coils (Hageman et al., 2008) we prepared a program to predict the stability of each tested pair. | Based on the published rules for the prediction of stability of coiled-coils (Hageman et al., 2008) we prepared a program to predict the stability of each tested pair. | ||
- | More about software (download and abstract) follow link | + | More about software (download and abstract) follow link [https://2009.igem.org/Team:Slovenia/Software HERE] |
<br><br> | <br><br> | ||
In order to restrict the number of possible combinations to manageable we limited the variability at defined positions to a limited set of residues that have been previously shown to stabilize the interactions. We started by choosing the most appropriate amino acids for the every position in sequence, similar as reported by (Bromley et al., 2009). Generally we followed the rule that hydrophobic residues at positions “a” or “d” and opposite charge at positions “e” and “g” on the equivalent positions between the two helices stabilize, whereas burial of polar residues (Asn) at positions “a” or “d” and the same charge at positions “e” and “g” destabilize the structure (see [https://2009.igem.org/Modular_protein_nanostructures.html Introduction]). In this simplified model we considered a limited number of variable residues only at positions that significantly affect the stability (“a”, “d”, “e”, “g”) and neglecting the effect on other positions (positions “b”, “c”, “f”). It has to be noted that the parameters to predict stability have been extracted from parallel dimers of basic leucine zipper (bZIP) peptides, therefore the predictions for this type of interaction are most reliable. Although we can probably estimate also contributions to the stability of antiparallel orientation, we tried to use experimentally parallel orientations whenever it was possible. | In order to restrict the number of possible combinations to manageable we limited the variability at defined positions to a limited set of residues that have been previously shown to stabilize the interactions. We started by choosing the most appropriate amino acids for the every position in sequence, similar as reported by (Bromley et al., 2009). Generally we followed the rule that hydrophobic residues at positions “a” or “d” and opposite charge at positions “e” and “g” on the equivalent positions between the two helices stabilize, whereas burial of polar residues (Asn) at positions “a” or “d” and the same charge at positions “e” and “g” destabilize the structure (see [https://2009.igem.org/Modular_protein_nanostructures.html Introduction]). In this simplified model we considered a limited number of variable residues only at positions that significantly affect the stability (“a”, “d”, “e”, “g”) and neglecting the effect on other positions (positions “b”, “c”, “f”). It has to be noted that the parameters to predict stability have been extracted from parallel dimers of basic leucine zipper (bZIP) peptides, therefore the predictions for this type of interaction are most reliable. Although we can probably estimate also contributions to the stability of antiparallel orientation, we tried to use experimentally parallel orientations whenever it was possible. | ||
Line 31: | Line 31: | ||
</html> | </html> | ||
<br> | <br> | ||
- | We extended this algorithm to pairs with different lenghts and register as well as antiparallel peptides - LINK | + | We extended this algorithm to pairs with different lenghts and register as well as antiparallel peptides - [https://2009.igem.org/Team:Slovenia/Software LINK] |
<br><br> | <br><br> | ||
- | Complete algorithm (in C, easily applicable to C++/Matlab/Mathematica): LINK | + | Complete algorithm (in C, easily applicable to C++/Matlab/Mathematica): [https://2009.igem.org/Team:Slovenia/Software LINK] |
- | and in Excel (just for pairs of with set register): LINK | + | and in Excel (just for pairs of with set register): [https://2009.igem.org/Team:Slovenia/Software LINK] |
<br><br> | <br><br> | ||
===Identification of a set of orthogonal coiled-coil pairs=== | ===Identification of a set of orthogonal coiled-coil pairs=== | ||
Line 48: | Line 48: | ||
First we calculated the interaction energy between all possible pairs of 8 different heptads. The interaction of coiled-coil segment is additive, therefore the interaction between larger coiled-coil segments, composed of several heptads could be deduced from the interacting contributions of all heptads. Next we generated a table where we calculated the interaction energy between all possible combinations of coiled-coil-forming segments composed of four heptads in both parallel and antiparallel orientation. | First we calculated the interaction energy between all possible pairs of 8 different heptads. The interaction of coiled-coil segment is additive, therefore the interaction between larger coiled-coil segments, composed of several heptads could be deduced from the interacting contributions of all heptads. Next we generated a table where we calculated the interaction energy between all possible combinations of coiled-coil-forming segments composed of four heptads in both parallel and antiparallel orientation. | ||
<br><br> | <br><br> | ||
- | Then we prepared another algorithm (download and abstract follow link | + | Then we prepared another algorithm (download and abstract follow link [https://2009.igem.org/Team:Slovenia/Software HERE]) to pick the best set of orthogonal pairs from the generated table which was selected to maximize the difference between the least stable desired pair and most stable undesired pair in the selected set. Predicted temperature difference between the least stable desired pair and most stable undesired pair was predicted to be more than 70 ºC. |
<br><br> | <br><br> | ||
As a result of this approach we found a system of 8 orthogonal designed peptides: P1, P2, P3, P4, P5, P6, P7 and P8, which form four pairs: P1+P2, P3+P4, P5+P6 and P7+P8. | As a result of this approach we found a system of 8 orthogonal designed peptides: P1, P2, P3, P4, P5, P6, P7 and P8, which form four pairs: P1+P2, P3+P4, P5+P6 and P7+P8. |
Revision as of 21:13, 21 October 2009
|
De novo design of orthogonal coiled-coil forming segments
In this way selected system of coiled-coil forming segments is called orthogonal. At the N-terminal segment of each coiled-coil segment we introduced a Ser-Pro-Glu-Asp tetrapeptide, which is composed of contributions of the preferred residues at the N-termini of coiled-coil helices (Sun et al., 2000; Cochran et al., 2001; Penel et al.,l 1999). Side chain of those residues can form hydrogen bond with polypeptide backbone of the residues in the fist turn, interact with a dipole moment of helix and have conformational properties compatible with the beginning of a helix and at the same time serve to prevent formation of an extended helix, that might form between the consecutive coiled-coil forming segments.
To predict coiled-coil interactions we examined two algorithms (see Hagemann et al., 2009; Fong et al., 2006) and in the end chose to follow Hageman et al.:
There are many factors that affect stability: hydrophobic burial, propensity, solubility, electrostatic interaction of flanking residues and others. This algorithm considers three factors: core, electrostatic and propensity which seem to have the key role in CC formation.
The calculation is shown in following scheme (shown on a fictional protein pair)(Figure 1):
Figure 1 We extended this algorithm to pairs with different lenghts and register as well as antiparallel peptides - LINK Complete algorithm (in C, easily applicable to C++/Matlab/Mathematica): LINK and in Excel (just for pairs of with set register): LINK Identification of a set of orthogonal coiled-coil pairs
Figure 2 Figure 2 above shows our approach: First we calculated the interaction energy between all possible pairs of 8 different heptads. The interaction of coiled-coil segment is additive, therefore the interaction between larger coiled-coil segments, composed of several heptads could be deduced from the interacting contributions of all heptads. Next we generated a table where we calculated the interaction energy between all possible combinations of coiled-coil-forming segments composed of four heptads in both parallel and antiparallel orientation. Then we prepared another algorithm (download and abstract follow link HERE) to pick the best set of orthogonal pairs from the generated table which was selected to maximize the difference between the least stable desired pair and most stable undesired pair in the selected set. Predicted temperature difference between the least stable desired pair and most stable undesired pair was predicted to be more than 70 ºC. As a result of this approach we found a system of 8 orthogonal designed peptides: P1, P2, P3, P4, P5, P6, P7 and P8, which form four pairs: P1+P2, P3+P4, P5+P6 and P7+P8. In the following chart we can see theoretically predicted melting temperatures (higher the temperature, the more stable the coiled-coil) Figure 3 |