Modeling microdomains: the surface area of globin helices.

The accessible surface areas of 53 high-resolution globin helices are correlated with molecular weight. The linear fit is assessed for statistical accuracy using a boot-strap analysis, and by comparison to the areas of 13 ideal polyalanine alpha-helices. The accessible area of the unfolded helices is compared with the folded values before helix-helix packing. An analytical physical model is presented to explain the correlation, and to provide an analytical value for the surface area parameter in the diffusion-collision model of protein folding.     

Assembly of amyloid protofibrils via critical oligomers--a novel pathway of amyloid formation

The amyloid formation of phosphoglycerate kinase (PGK) was investigated by static and dynamic light-scattering. The time-course of the scattering intensity and the hydrodynamic radius scale with initial monomer concentration in a linear fashion over a range of about 50 in concentration. This sets limits on theories for aggregation kinetics that can be used, and points towards irreversible, cascade type models. In addition, circular dichroism (CD) was used to monitor the transition between a predominantly alpha-helical spectrum to a beta-sheet enriched one. The time-course of the CD also proves to scale linearly with initial monomer concentration. Electron microscopy shows that small oligomers as well as protofibrils are present during aggregation. The found coupling between growth of intermediates and acquisition of beta-sheet structure is interpreted in terms of a generalized diffusion-collision model, where stabilization of beta-strands takes place by intermolecular interactions.     

The early folding kinetics of apomyoglobin.

The folding pathway of apomyoglobin has been experimentally shown to have early kinetic intermediates involving the A, B, G, and H helices. The earliest detected kinetic events occur on a ns to micros time scale. We show that the early folding kinetics of apomyoglobin may be understood as the association of nascent helices through a network of diffusion-collision-coalescence steps G + H <--> GH + A <--> AGH + B <--> ABGH obtained by solving the diffusion-collision model in a chemical kinetics approximation. Our reproduction of the experimental results indicates that the model is a useful way to analyze folding data. One prediction from our fit is that the nascent A and H helices should be relatively more helix-like before coalescence than the other apomyoglobin helices.     

Quantifying the structural requirements of the folding transition state of protein A and other systems

The B-domain of protein A is a small three-helix bundle that has been the subject of considerable experimental and theoretical investigation. Nevertheless, a unified view of the structure of the transition-state ensemble (TSE) is still lacking. To characterize the TSE of this surprisingly challenging protein, we apply a combination of psi analysis (which probes the role of specific side-chain to side-chain contacts) and kinetic H/D amide isotope effects (which measures hydrogen-bond content), building upon previous studies using mutational phi analysis (which probes the energetic influence of side-chain substitutions). The second helix is folded in the TSE, while helix formation appears just at the carboxy and amino termini of the first and third helices, respectively. The experimental data suggest a homogenous yet plastic TS with a native-like topology. This study generalizes our earlier conclusion, based on two larger alpha/beta proteins, that the TSEs of most small proteins achieve approximately 70% of their native state's relative contact order. This high percentage limits the degree of possible TS heterogeneity and requires a reevaluation of the structural content of the TSE of other proteins, especially when they are characterized as small or polarized.     

Diffusing and colliding: the atomic level folding/unfolding pathway of a small helical protein

Proteins with ultra-fast folding/unfolding kinetics are excellent candidates for study by molecular dynamics. Here, we describe such simulations of a three helix bundle protein, the engrailed homeodomain (En-HD), which folds via the diffusion-collision model. The unfolding pathway of En-HD was characterized by seven simulations of the protein and 12 simulations of its helical fragments yielding over 1.1 micros of simulation time in water. Various conformational states along the unfolding pathway were identified. There is the compact native-like transition state, a U-shaped helical intermediate and an unfolded state with dynamic helical segments. Each of these states is in good agreement with experimental data. Examining these states as well as the transitions between them, we find the role of long-range tertiary contacts, specifically salt-bridges, important in the folding/unfolding pathway. In the folding direction, charged residues form long-range tertiary contacts before the hydrophobic core is formed. The formation of HII is assisted by a specific salt-bridge and by non-specific (fluctuating) tertiary contacts, which we call contact-assisted helix formation. Salt-bridges persist as the protein approaches the transition state, stabilizing HII until the hydrophobic core is formed. To complement this information, simulations of fragments of En-HD illustrate the helical propensities of the individual segments. By thermal denaturation, HII proved to be the least stable helix, unfolding in less than 450 ps at high temperature. We observed the low helical propensity of C-terminal residues from HIII in fragment simulations which, when compared to En-HD unfolding simulations, link the unraveling of HIII to the initial event that drives the unfolding of En-HD.