Folding thermodynamics of a model three-helix-bundle protein

The calculated folding thermodynamics of a simple off-lattice three-helix-bundle protein model under equilibrium conditions shows the experimentally observed protein transitions: a collapse transition, a disordered-to-ordered globule transition, a globule to native-state transition, and the transition from the active native state to a frozen inactive state. The cooperativity and physical origin of the various transitions are explored with a single “optimization” parameter and characterized with the Lindemann criterion for liquid versus solid-state dynamics. Below the folding temperature, the model has a simple free energy surface with a single basin near the native state; the surface is similar to that calculated from a simulation of the same three-helix-bundle protein with an all-atom representation [Boczko, E. M. & Brooks III, C. L. (1995) Science 269, 393–396].

The folding mechanism of larger model proteins: role of native structure.

The folding mechanism of a 125-bead heteropolymer model for proteins is investigated with Monte Carlo simulations on a cubic lattice. Sequences that do and do not fold in a reasonable time are compared. The overall folding behavior is found to be more complex than that of models for smaller proteins. Folding begins with a rapid collapse followed by a slow search through the semi-compact globule for a sequence-dependent stable core with about 30 out of 176 native contacts which serves as the transition state for folding to a near-native structure. Efficient search for the core is dependent on structural features of the native state. Sequences that fold have large amounts of stable, cooperative structure that is accessible through short-range initiation sites, such as those in anti-parallel sheets connected by turns. Before folding is completed, the system can encounter a second bottleneck, involving the condensation and rearrangement of surface residues. Overly stable local structure of the surface residues slows this stage of the folding process. The relation of the results from the 125-mer model studies to the folding of real proteins is discussed.

Role of hydrophobic interactions and desolvation in determining the structural properties of a model alpha beta peptide.

Model AB, a 20-amino acid peptide that was designed to adopt an alpha beta tertiary structure stabilized by hydrophobic interactions between residues in adjacent helical and extended segments, exhibited large pKa shifts of several ionizable groups and slow hydrogen/deuterium exchange rates of nearly all the peptide amide groups [Butcher, D. J., Bruch, M. D. & Moe, G. T. (1995) Biopolymers 36, 109-120]. These properties, which depend on structure and hydration, are commonly observed in larger proteins but are quite unusual for small peptides. To identify which of several possible features of the peptide design are most important in determining these properties, several closely related analogs of Model AB were characterized by CD and NMR spectroscopy. The results show that hydrophobic interactions between adjacent helical and extended segments are structure-determining and have the additional effect of altering water-peptide interactions over much of the peptide surface. These results may have important implications for understanding mechanisms of protein folding and for the design of independently folding peptides.

A constant radius of curvature model for the organization of DNA in toroidal condensates.

Toroidal DNA condensates have received considerable attention for their possible relationship to the packaging of DNA in viruses and in general as a model of ordered DNA condensation. A spool-like model has primarily been supported for DNA organization within toroids. However, our observations suggest that the actual organization may be considerably different. We present an alternate model in which DNA for a given toroid is organized within a series of equally sized contiguous loops that precess about the toroid axis. A related model for the toroid formation process is also presented. This kinetic model predicts a distribution of toroid sizes for DNA condensed from solution that is in good agreement with experimental data.