Structure, dynamics and energetics of initiation sites in protein folding: I. Analysis of a 1 ns molecular dynamics trajectory of an early folding unit in water: the helix I/loop I-fragment of barnase

Molecular dynamics simulations of the unfolding of barnase in water and 8 M aqueous urea

Molecular dynamics simulations of barnase have been conducted both in water and in 8 M urea solution for 500 ps at 25 degrees C and for 2000 ps at 85 degrees C. The final structure of the aqueous simulation at room temperature matches closely the structure obtained by NMR and the experimentally observed protections from isotopic exchange. The comparison of the structures generated by the aqueous simulation at 85 degrees C reveals a trajectory composed of groups of geometrically related structures separated by narrow regions of rapid change in structure. The first of these regions displays changes in backbone rmsd to the crystal structure and solvent-accessible area suggestive of a transition state, while the properties observed during the final 300 ps of the simulation are consistent with a stable intermediate. These assignments were confirmed by calculation of the "progress along the reaction coordinate" phi-values using an empirical equation based on a linear response method. The pathway of unfolding defined in this fashion agrees well with the experimental results of site-directed mutagenesis in terms of secondary structure content of the transition state and the intermediate and reproduces the relative stability of the different elements of secondary structure. The results of the simulations in urea suggest a mechanism at the molecular level for its well-known enhancement of the denaturation of proteins. The analysis of radial distribution functions shows that the first solvation shell of the protein is enriched in urea relative to the bulk solvent. The displacement of water molecules allows greater exposure of hydrophobic side chains, as witnessed particularly in the analysis of solvent-accessible surface areas at the higher temperature. Almost all urea molecules in the first shell form at least one hydrogen bond with the protein. They provide a more favorable environment for accommodation of the remaining water molecules, and they facilitate the separation of secondary structure elements by acting as a bridge between groups previously forming intraprotein hydrogen bonds.

Comparison of database potentials and molecular mechanics force fields

The advantages and disadvantages of database and molecular mechanics force fields for the study of macromolecules are compared, with emphasis on the ability to distinguish between correct and incorrect structures. Molecular mechanics force fields have the advantage of resting on a clear theoretical basis, permitting an in-depth analysis of different contributions. On the other hand, large simplifications are necessary for tractable computing, and there has so far been little effective testing at the macromolecular level. Database potentials allow greater freedom of functional form and have been shown to be effective at discriminating between correct and incorrect complete structures. The principal negative is a controversial relationship to free energy. More testing and comparison of both sorts of potential are needed.

Molecular dynamics simulations of N-terminal peptides from a nucleotide binding protein

Molecular dynamics (MD) simulations of N-terminal peptides from lactate dehydrogenase (LDH) with increasing length and individual secondary structure elements were used to study their stability in relation to folding. Ten simulations of 1-2 ns of different peptides in water starting from the coordinates of the crystal structure were performed. The stability of the peptides was compared qualitatively by analyzing the root mean square deviation (RMSD) from the crystal structure, radius of gyration, secondary and tertiary structure, and solvent accessible surface area. In agreement with earlier MD studies, relatively short (< 15 amino acids) peptides containing individual secondary structure elements were generally found to be unstable; the hydrophobic alpha 1-helix of the nucleotide binding fold displayed a significantly higher stability, however. Our simulations further showed that the first beta alpha beta supersecondary unit of the characteristic dinucleotide binding fold (Rossmann fold) of LDH is somewhat more stable than other units of similar length and that the alpha 2-helix, which unfolds by itself, is stabilized by binding to this unit. This finding suggests that the first beta alpha beta unit could function as an N-terminal folding nucleus, upon which the remainder of the polypeptide chain can be assembled. Indeed, simulations with longer units (beta-alpha-beta-alpha and beta-alpha-beta-alpha beta-beta) showed that all structural elements of these units are rather stable. The outcome of our studies is in line with suggestions that folding of the N-terminal portion of LDH in vivo can be a cotranslational process that takes place during the ribosomal peptide synthesis.

Structure of a cyclic peptide with a catalytic triad, determined by computer simulation and NMR spectroscopy

We report the design of a cyclic, eight-residue peptide that possesses the catalytic triad residues of the serine proteases. A manually built model has been relaxed by 0.3 ns of molecular dynamics simulation at room temperature, during which no major changes occurred in the peptide. The molecule has been synthesised and purified. Two-dimensional NMR spectroscopy provided 35 distance and 7 torsion angle constraints, which were used to determine the three-dimensional structure. The experimental conformation agrees with the predicted one at the beta-turn, but deviates in the arrangement of the disulphide bridge that closes the backbone to a ring. A 1.2 ns simulation at 600 K provided extended sampling of conformation space. The disulphide bridge reoriented into the experimental arrangement, producing a minimum backbone rmsd from the experimental conformation of 0.8 A. At a later stage in the simulation, a transition at Ser3 produced more pronounced high-temperature behaviour. The peptide hydrolyses p-nitrophenyl acetate about nine times faster than free histidine.