The trigger sequence in the GCN4 leucine zipper: α-helical propensity and multistate dynamics of folding and dimerization

Modulation of Coiled-Coil Dimer Stability through Surface Residues while Preserving Pairing Specificity

The coiled-coil dimer is a widespread protein structural motif and, due to its designability, represents an attractive building block for assembling modular nanostructures. The specificity of coiled-coil dimer pairing is mainly based on hydrophobic and electrostatic interactions between residues at positions a, d, e, and g of the heptad repeat. Binding affinity, on the other hand, can also be affected by surface residues that face away from the dimerization interface. Here we show how design of the local helical propensity of interacting peptides can be used to tune the stabilities of coiled-coil dimers over a wide range. By designing intramolecular charge pairs, regions of high local helical propensity can be engineered to form trigger sequences, and dimer stability is adjusted without changing the peptide length or any of the directly interacting residues. This general principle is demonstrated by a change in thermal stability by more than 30 °C as a result of only two mutations outside the binding interface. The same approach was successfully used to modulate the stabilities in an orthogonal set of coiled-coils without affecting their binding preferences. The stability effects of local helical propensity and peptide charge are well described by a simple linear model, which should help improve current coiled-coil stability prediction algorithms. Our findings enable tuning the stabilities of coiled-coil-based building modules match a diverse range of applications in synthetic biology and nanomaterials.

Structural propensities and entropy effects in peptide helix-coil transitions

The helix-coil transition in peptides is a critical structural transition leading to functioning proteins. Peptide chains have a large number of possible configurations that must be accounted for in statistical mechanical investigations. Using hydrogen bond and local helix propensity interaction terms, we develop a method for obtaining and incorporating the degeneracy factor that allows the exact calculation of the partition function for a peptide as a function of chain length. The partition function is used in calculations for engineered peptide chains of various lengths that allow comparison with a variety of different types of experimentally measured quantities, such as fraction of helicity as a function of both temperature and chain length, heat capacity, and denaturation studies. When experimental sensitivity in helicity measurements is properly accounted for in the calculations, the calculated curves fit well with the experimental curves. We determine values of interaction energies for comparison with known biochemical interactions, as well as quantify the difference in the number of configurations available to an amino acid in a random coil configuration compared to a helical configuration.

Free-energy landscapes and thermodynamic parameters of complex molecules from nonequilibrium simulation trajectories

Thermodynamic parameters such as free energies and heat capacities are important quantities for understanding processes involving structural transitions in complex molecules such as proteins. Computational investigations provide simulated data that can be used for calculating thermodynamic parameters. However, calculations give accurate results only if the simulations sample all of configuration space with the appropriate temperature-dependent Boltzmann equilibrium probabilities. For many systems, truly comprehensive sampling of configuration space is not computationally feasible. We present an approximation technique for the calculations that will give accurate values for thermodynamic parameters when the data is incomplete. Our work is applicable to systems in which there are two distinct, important regions of configuration space that must be sampled. Importantly, the results are also valid when the system is more complex than two-state systems. Transition pathways that involve intermediate configurations between two stable regions are allowed in this treatment, and therefore the results are valid for multistate systems.

Dry and wet interfaces: influence of solvent particles on molecular recognition

We present a coarse-grained lattice model to study the influence of water on the recognition process of two rigid proteins. The basic model is formulated in terms of the hydrophobic effect. We then investigate several modifications of our basic model showing that the selectivity of the recognition process can be enhanced by considering the explicit influence of single solvent particles. When the number of cavities at the interface of a protein-protein complex is fixed as an intrinsic geometric constraint, there typically exists a characteristic fraction that should be filled with water molecules such that the selectivity exhibits a maximum. In addition, the optimum fraction depends on the hydrophobicity of the interface so that one has to distinguish between dry and wet interfaces.