Classical and quantum mechanical/molecular mechanical molecular dynamics simulations of alanine dipeptide in water: comparisons with IR and vibrational circular dichroism spectra

We have implemented the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulations of alanine dipeptide in water along with the polarizable and nonpolarizable classical MD simulations with different models of water. For the QM/MM MD simulation, the alanine dipeptide is treated with the AM1 or PM3 approximations and the fluctuating solute dipole moment is calculated by the Mulliken population analysis. For the classical MD simulations, the solute is treated with the polarizable or nonpolarizable AMBER and polarizable CHARMM force fields and water is treated with the TIP3P, TIP4P, or TIP5P model. It is found that the relative populations of right-handed alpha-helix and extended beta and P(II) conformations in the simulation trajectory strongly depend on the simulation method. For the QM/MM MD simulations, the PM3/MM shows that the P(II) conformation is dominant, whereas the AM1/MM predicts that the dominant conformation is alpha(R). Polarizable CHARMM force field gives almost exclusively P(II) conformation and other force fields predict that both alpha-helical and extended (beta and P(II)) conformations are populated with varying extents. Solvation environment around the dipeptide is investigated by examining the radial distribution functions and numbers and lifetimes of hydrogen bonds. Comparing the simulated IR and vibrational circular dichroism spectra with experimental results, we concluded that the dipeptide adopts the P(II) conformation and PM3/MM, AMBER03 with TIP4P water, and AMBER polarizable force fields are acceptable for structure determination of the dipeptide considered in this paper.     

Global energy minimization of alanine dipeptide via barrier function methods

This paper presents an interior point method to determine the minimum energy conformation of alanine dipeptide. The CHARMM energy function is minimized over the internal coordinates of the atoms involved. A barrier function algorithm to determine the minimum energy conformation of peptides is proposed. Lennard-Jones 6-12 potential which is used to model the van der Waals interactions in the CHARMM energy equation is used as the barrier function for this algorithm. The results of applying the algorithm for the alanine dipeptide structure as a function of varying number of dihedral angles are reported, and they are compared with that obtained from genetic algorithm approach. In addition, the results for polyalanine structures are also reported.     

Alpha-gamma hybrid peptides that contain the conformationally constrained gabapentin residue: characterization of mimetics of chain reversals

The crystal structures of four dipeptides that contain the stereochemically constrained gamma-amino acid residue gabapentin (1-(aminomethyl)cyclohexaneacetic acid Gpn) are described. The molecular conformation of Piv-Pro-Gpn-OH (1), reveals a beta-turn mimetic conformation, stabilized by a ten atom C[bond]H...O hydrogen bond between the Piv CO group and the pro S hydrogen of the Gpn CH(2)[bond]CO group. The peptides Boc-Gly-Gpn-OH (2), Boc-Aib-Gpn-OH (3), and Boc-Aib-Gpn-OMe (4) form compact, folded structures, in which a distinct reversal of polypeptide chain direction is observed. In all cases, the Gpn residue adopts a gauche,gauche (g,g) conformation about the C(gamma)[bond]C(beta) (theta(1)) and C(beta)[bond]C(alpha) (theta(2)) bonds. Two distinct Gpn conformational families are observed. In peptides 1 and 3, the average backbone torsion angle values for the Gpn residue are phi=98 degrees, theta(1)=-62 degrees, theta(2)=-73 degrees, and psi=79 degrees, while in peptide 2 and 4 the average values are phi=-103 degrees, theta(1)=-46 degrees, theta(2)=-49 degrees, and psi=-92 degrees. In the case of 1 and 3, an intramolecular nine-membered O[bond]H...O hydrogen bond is formed between the C[double bond]O of the preceding residue and the terminal carboxylic acid OH group. All four alpha-gamma dipeptide sequences yield compact folded backbone conformations; this suggests that the Gpn residue may be employed successfully in the design of novel folded structures.     

Extending the treatment of backbone energetics in protein force fields: limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations

Computational studies of proteins based on empirical force fields represent a powerful tool to obtain structure-function relationships at an atomic level, and are central in current efforts to solve the protein folding problem. The results from studies applying these tools are, however, dependent on the quality of the force fields used. In particular, accurate treatment of the peptide backbone is crucial to achieve representative conformational distributions in simulation studies. To improve the treatment of the peptide backbone, quantum mechanical (QM) and molecular mechanical (MM) calculations were undertaken on the alanine, glycine, and proline dipeptides, and the results from these calculations were combined with molecular dynamics (MD) simulations of proteins in crystal and aqueous environments. QM potential energy maps of the alanine and glycine dipeptides at the LMP2/cc-pVxZ//MP2/6-31G* levels, where x = D, T, and Q, were determined, and are compared to available QM studies on these molecules. The LMP2/cc-pVQZ//MP2/6-31G* energy surfaces for all three dipeptides were then used to improve the MM treatment of the dipeptides. These improvements included additional parameter optimization via Monte Carlo simulated annealing and extension of the potential energy function to contain peptide backbone phi, psi dihedral crossterms or a phi, psi grid-based energy correction term. Simultaneously, MD simulations of up to seven proteins in their crystalline environments were used to validate the force field enhancements. Comparison with QM and crystallographic data showed that an additional optimization of the phi, psi dihedral parameters along with the grid-based energy correction were required to yield significant improvements over the CHARMM22 force field. However, systematic deviations in the treatment of phi and psi in the helical and sheet regions were evident. Accordingly, empirical adjustments were made to the grid-based energy correction for alanine and glycine to account for these systematic differences. These adjustments lead to greater deviations from QM data for the two dipeptides but also yielded improved agreement with experimental crystallographic data. These improvements enhance the quality of the CHARMM force field in treating proteins. This extension of the potential energy function is anticipated to facilitate improved treatment of biological macromolecules via MM approaches in general.     

The strength with which a peptide group can form a hydrogen bond varies with the internal conformation of the polypeptide chain

The strength of the H-bond formed between a dipeptide and a proton acceptor molecule is assessed by correlated ab initio quantum calculations for a broad range of different conformations of the dipeptide. The H-bond energy is very sensitive to the internal (phi,psi) angles, even when the geometry of the H-bond does not vary significantly from one conformation to another. This result indicates that the peptide NH is a much less potent proton donor in certain conformations than in others. In particular, extended conformations of a polypeptide are capable of only weak H-bonds. Thus, the interstrand NH...O H-bonds in parallel and antiparallel beta-sheets are expected to be significantly weaker than those found in other conformations, such as helices, ribbons, and beta-bends, even if the H-bond geometries are similar.     

An algorithm to find minimum free-energy paths using umbrella integration

The calculation of free-energy barriers by umbrella sampling and many other methods is hampered by the necessity for an a priori choice of the reaction coordinate along which to sample. We avoid this problem by providing a method to search for saddle points on the free-energy surface in many coordinates. The necessary gradients and Hessians of the free energy are obtained by multidimensional umbrella integration. We construct the minimum free-energy path by following the gradient down to minima on the free-energy surface. The change of free energy along the path is obtained by integrating out all coordinates orthogonal to the path. While we expect the method to be applicable to large systems, we test it on the alanine dipeptide in vacuum. The minima, transition states, and free-energy barriers agree well with those obtained previously with other methods.     

Improved treatment of the protein backbone in empirical force fields

Empirical force field-based calculations of proteins, including protein-folding studies, have improved our understanding of the relationship of their structure to their biological function. However, limitations in the accuracy of empirical force fields in the treatment of the peptide backbone exist. Presented is a grid correction approach to improve the treatment of the peptide backbone phi/psi conformational energies. Inclusion of this correction with the CHARMM22 all-atom protein force field is shown to lead to significant improvement in the treatment of the conformational energies of both the peptide model compound, the alanine dipeptide, and of proteins in their crystal environment. The developed approach is suggested to lead to significant improvements in the accuracy of empirical force fields to treat peptides and proteins.     

A theoretical study of conformational properties of N-methyl azapeptide derivatives

The conformational properties of azapeptide derivatives, Ac-azaGly-NHMe (1), Ac-azaAla-NHMe (2), Ac-NMe-azaGly-NHMe (3), Ac-NMe-azaAla-NHMe (4), Ac-azaGly-NMe(2) (5), Ac-azaAla-NMe(2) (6), Ac-NMe-azaGly-NMe(2) (7), and Ac-NMe-azaAla-NMe(2) (8), were systematically examined by using ab initio MO and DFT methods. Structural perturbations in azapeptides resulting from cyclic substitution of a methyl group at three N-positions of an azaamino acid were studied on the basis of the structure of the simplest model azapeptide, 1. Potential energy surfaces were generated at the HF/6-31G level for 1-4 and at the HF/6-31G//HF/3-21G level for 5-8 by rotating two key dihedral angles (phi, psi) in increments of 30 degrees. The backbone (phi, psi) angles of the minima for 1-4 are observed at the i + 2 position to form the betaI(I')-, betaII(II')-, betaVI-turns or the polyproline II structure according to the orientation of the acetyl group and the positions of the N-methyl groups. Compounds 5-8 coupled to a secondary amine were found to preferentially adopt polyproline II, betaI(III)-turn, or alpha-helical structure or even extended conformations depending on the orientation of the acetyl group and the positions of the N-methyl groups. Furthermore, N-methyl groups, depending on their positions, were found to affect the orientation of the amide group in the lowest energy conformations, the pyramidality of the N2 atom, and the bond length in azapeptide derivatives. These unique theoretical conformations of N-methyl azapeptide derivatives could be utilized in the definite design of secondary structure for peptides and proteins, and in the development of new drugs and molecular machines.     

The relationship between peptide plane rotation (PPR) and similar conformations

Currently, several energy functions and conformational search methods have been developed that are based on the observed distribution of phi and psi angles in protein structures. The definition of phi and psi angles is directly related to the orientation of the peptide plane (CA CO NH CA). Starting from one conformation and rotating a single peptide plane, the angles psi for one residue and phi for the consecutive residue that are linked by the peptide plane, display a continuous range of values within one global conformation. When peptide plane rotation is analyzed in several different conformations generated from a restricted conformation database, a large number of these conformations are related. Based on these observations, a new simplified all-atom representation for protein folding simulations is presented where only one torsion angle variable is required for each residue. The underlying theme of this article is that conformational search methods using phi and psi torsion space, search through many redundant conformations. These conformations are related by anticorrelated torsion changes of peptide plane rotations. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 947–955, 1999     

Dihedral angles of trialanine in D2O determined by combining FTIR and polarized visible Raman spectroscopy

We have measured the polarized visible Raman and FTIR spectra of trialanine and triglycine in D(2)O at acid, neutral, and alkaline pD. From the Raman spectra we obtained the isotropic and the anisotropic scattering. A self-consistent spectral analysis of the region between 1550 and 1800 cm(-1) was carried out to obtain the intensities, frequencies, and halfwidths of the respective amide I bands. A model was developed by means of which the intensity ratios of the amide I bands in all spectra and the respective frequency differences were utilized to determine the orientational angle theta between the peptide groups and the strength of excitonic coupling between the corresponding amide I modes. By exploiting results from a recent ab initio study on triglycine (Torii, H; Tasumi, M. J. Raman Spectrosc. 1998, 29, 81), we used these parameters to determine the dihedral angles phi and psi between the peptide groups. Our results show that trialanine adopts a 3(1)-helical structure in D(2)O for all of its three protonation states. The structure is insensitive to the carboxylate protonation and changes only slightly with N-terminal protonation. Triglycine is structurally more heterogeneous in the zwitterionic and the cationic state. Our spectral analysis suggests that 3(1)-helices coexist with right-handed alpha-helical and/or with beta-turn conformations. The N-terminal protonation stabilizes the 3(1)-structure. Our study provides compelling evidence that tripeptides adopt stable conformations in aqueous solution and that they are suitable model systems to investigate the initiation of secondary structure formation.     

Estimating kinetic rates from accelerated molecular dynamics simulations: alanine dipeptide in explicit solvent as a case study

Molecular dynamics (MD) simulation is the standard computational technique used to obtain information on the time evolution of the conformations of proteins and many other molecular systems. However, for most biological systems of interest, the time scale for slow conformational transitions is still inaccessible to standard MD simulations. Several sampling methods have been proposed to address this issue, including the accelerated molecular dynamics method. In this work, we study the extent of sampling of the phi/psi space of alanine dipeptide in explicit water using accelerated molecular dynamics and present a framework to recover the correct kinetic rate constant for the helix to beta-strand transition. We show that the accelerated MD can drastically enhance the sampling of the phi/psi conformational phase space when compared to normal MD. In addition, the free energy density plots of the phi/psi space show that all minima regions are accurately sampled and the canonical distribution is recovered. Moreover, the kinetic rate constant for the helix to beta-strand transition is accurately estimated from these simulations by relating the diffusion coefficient to the local energetic roughness of the energy landscape. Surprisingly, even for such a low barrier transition, it is difficult to obtain enough transitions to accurately estimate the rate constant when one uses normal MD.     

First-principle computational study on the full conformational space of L-threonine diamide, the energetic stability of cis and trans isomers

First-principle computations were carried out on the conformational space of trans and cis peptide bond isomers of HCO-Thr-NH2. Using the concept of multidimensional conformational analysis (MDCA), geometry optimizations were performed at the B3LYP/6-31G(d) level of theory, and single-point energies as well as thermodynamic functions were calculated at the G3MP2B3 level of theory for the corresponding optimized structures. Two backbone Ramachandran-type potential energy surfaces (PESs) were computed, one each for the cis and trans isomers, keeping the side chain at the fully extended orientation (chi1=chi2=anti). Similarly, two side chain PESs for the cis and trans isomers were generated for the (phi=psi=anti) orientation corresponding to approximately the betaL backbone conformation. Besides correlating the relative Gibbs free energy of the various stable conformations with the number of stabilizing hydrogen bonds, the process of trans-->cis isomerization is discussed in terms of intrinsic stabilities as measured by the computed thermodynamic functions.     

Free energy surfaces for the alpha(1 --> 4)-glycosidic linkage: implications for polysaccharide solution structure and dynamics

We present a potential of mean force surface for rotation about phi and psi dihedral angles of the alpha(1 --> 4)-glycosidic linkage in the maltose disaccharide (4-O-alpha-d-glucopyranosyl-d-glucopyranose) in aqueous solution. Comparison of the vacuum and solution free energy surfaces for maltose shows the principal effects of water to be an increase in the rotational freedom of the alpha(1 --> 4) linkage brought about by lowering the energy barrier for syn to anti conformational changes as well as expansion of the range of low-energy phi,psi conformations. This free energy analysis thus provides a thermodynamic and conformational rationale for the effects of water on alpha(1 --> 4)-linked polysaccharides and carbohydrate glasses.     

Conformational preferences and cis-trans isomerization of L-lactic acid residue

The conformational study on N-acetyl- N'-methylamide of l-lactic acid (Ac-Lac-NHMe, the Lac dipeptide) is carried out using ab initio HF and density functional methods with the self-consistent reaction field method to explore its backbone conformational preferences and cis-trans isomerization for the depsipeptide with an ester bond in the gas phase and in solution. In the gas phase and in chloroform, the conformation tB with a trans depsipeptide bond is most preferred for the Lac dipeptide, whose backbone torsion angles are phi approximately -150 degrees and psi approximately -5 degrees , juxtaposed to those of the 3 10-helical structure. The larger shift in phi is brought to reduce the repulsion between the two carbonyl carbons of the acetyl and NHMe groups. However, the polyproline II-like tF conformation becomes more populated and the relative stability of conformation tB decreases significantly as the solvent polarity increases. This may be ascribed to weakening a C(5) hydrogen bond between the depsipeptidyl oxygen and the carboxyl amide hydrogen that plays a role in stabilizing the conformation tB in the gas phase and in chloroform. The cis populations about the depsipeptide bond are nearly negligible in the gas phase and in solution. The rotational barriers to the cis-trans isomerization of the depsipeptide bond for the Lac dipeptide are calculated to be about 11 kcal/mol, which is about half of those for the Ala dipeptide, although they increase somewhat with the increase of solvent polarity. The cis-trans isomerization of the depsipeptide bond proceeds through either clockwise or anticlockwise rotations with torsion angles of about +90 degrees or -90 degrees , respectively, in the gas phase and in solution, whereas it has been known that the isomerization proceeds through only the clockwise rotation for alanyl and prolyl peptide bonds. The pertinent distances between the depsipeptidyl oxygen and the carboxyl amide hydrogen can describe the role of this hydrogen bond in stabilizing the transition state structures in the gas phase and in solution.     

Flexibility of "polyunsaturated fatty acid chains" and peptide backbones: A comparative ab initio study

The conformational properties of omega-3 type of polyunsaturated fatty acid (PUFA) chains and their fragments were studied using Hartree-Fock (RHF/3-21G) and DFT (B3LYP/6-31G(d)) methods. Comparisons between a unit (U) fragment of the PUFA chain and a mono N-Ac-glycine-NHMe residue show that both structures have the same sequence of sp2-sp3-sp2 atoms. The flexibility of PUFA originates in the internal rotation about the above pairs of sigma bonds. Therefore, potential energy surfaces (PESs) were generated by a scan around the terminal dihedral angles (phi t1 and phi t2) as well as the phi 1 and psi 1 dihedrals of both 1U congeners (Me-CHCH-CH2-CHCHMe and MeCONH-CH2-CONHMe) at the RHF/3-21G level of theory. An interesting similarity was found in the flexibility between the cis allylic structure and the trans peptide models. A flat landscape can be seen in the cis 1U (hepta-2,5-diene) surface, implying that several conformations are expected to be found in this (PES). An exhaustive search carried out on the 1U and 2U models revealed that straight chain structures such as trans and cis beta (phi 1 approximately psi 1 approximately 120 degrees; phi 2 approximately psi 2 approximately -120 degrees) or trans and cis extended (phi 1 approximately psi 1 approximately phi 2 approximately psi 2 approximately 120 degrees) can be formed at the lowest energy of both isomers. However, forming helical structures, such as trans helix (phi 1 approximately -120 degrees, psi 1 approximately 12 degrees; phi 2 approximately -120 degrees, psi 2 approximately 12 degrees) or cis helix (phi 1 approximately -130 degrees, psi 1 approximately 90 degrees; phi 2 approximately -145 degrees, psi 2 approximately 90 degrees) will require more energy. These six conformations, found in 2U, were selected to construct longer chains such as 3U, 4U, 5U, and 6U to obtain the thermochemistry of secondary structures. The variation in the extension or compression of the chain length turned out to be a factor of 2 between the helical and nonhelical structures. The inside diameter of the "tube" of cis helix turned out to be 3.5 A after discounting the internal H atoms. Thermodynamic functions were computed at the B3LYP/6-311+G(2d,p)//B3LYP/6-31G(d). The cis-trans isomerization energy of 1.7 +/- 0.2 kcal mol(-1) unit(-1) for all structure pairs indicates that the conformer selection was consistent. A folding energy of 0.5 +/- 0.1 kcal mol(-1) unit(-1) has been extracted from the energy comparison of the helices and most extended nonhelical structures. The entropy change associated with the folding (Delta S(folding)) is decreases faster with the degree of polymerization (n) for the cis than for the trans isomer. As a consequence, the linear relationships between (Delta G(folding)) and n for the cis and trans isomer crossed at about n = 3. This suggested that the naturally occurring cis isomer less ready to fold than the trans isomer since a greater degree of organization is exhibited by the cis isomer during the folding process. The result of this work leads to the question within the group additivity rule: could the method applied in our study of the folding of polyallylic hydrocarbons be useful in investigating the thermochemistry of protein folding?     

Role of solvent in determining conformational preferences of alanine dipeptide in water

Evidence from a variety of spectroscopic probes indicates that (phi, psi) values corresponding to the left-handed polyproline II helix (P(II)) are preferred for short alanine-based peptides in water. On the basis of results from theoretical studies, it is believed that the observed preference is dictated by favorable peptide-solvent interactions, which are realized through formation of optimal hydrogen-bonding water bridges between peptide donor and acceptor groups. In the present study, we address this issue explicitly by analyzing the hydration structure and thermodynamics of 16 low-energy conformers of the alanine dipeptide (N-acetylalanine-N'-methylamide) in liquid water. Monte Carlo simulations in the canonical ensemble were performed under ambient conditions with all-atom OPLS parameters for the alanine dipeptide and the TIP5P model for water. We find that the number of hydrogen-bonded water molecules connecting the peptide group donor and acceptor atoms has no effect on the solvation thermodynamics. Instead, the latter are determined by the work done to fully hydrate the peptide. This work is minimal for conformations that are characterized by a minimal overlap of the primary hydration shells around the peptide donor and acceptor atoms. As a result, peptide-solvent interactions favor "compact" conformations that do not include P(II)-like geometries. Our main conclusion is that the experimentally observed preference for P(II) does not arise due to favorable direct interactions between the peptide and water molecules. Instead, the latter act to unmask underlying conformational preferences that are a consequence of minimizing intrapeptide steric conflicts.     

A conformational dynamics study of alpha-l-Rhap-(1-->2)[alpha-l-Rhap-(1-->3)]-alpha-l-Rhap-OMe in solution by NMR experiments and molecular simulations

The conformational preference of alpha-l-Rhap-(1-->2)[alpha-l-Rhap-(1-->3)]-alpha-l-Rhap-OMe in solution has been studied by NMR spectroscopy using one-dimensional (1)H,(1)H T-ROESY experiments and measurement of trans-glycosidic (3)J(C,H) coupling constants. Molecular dynamics (MD) simulations with a CHARMM22 type of force field modified for carbohydrates were performed with water as the explicit solvent. The homonuclear cross-relaxation rates, interpreted as effective proton-proton distances, were compared to those obtained from simulation. Via a Karplus torsional relationship, (3)J(C,H) values were calculated from simulation and compared to experimental data. Good agreement was observed between experimental data and the MD simulation, except for one inter-residue T-ROE between protons in the terminal sugar residues. The results show that the trisaccharide exhibits substantial conformational flexibility, in particular along the psi glycosidic torsion angles. Notably, for these torsions, a high degree of correlation (77%) was observed in the MD simulation revealing either psi(2)(+) psi(3)(+) or psi(2)(-)psi(3)(-) states. The simulations also showed that non-exoanomeric conformations were present at the phi torsion angles, but to a limited extent, with the phi(3) state populated to a larger extent than the phi(2) state. Further NMR analysis of the trisaccharide by translational diffusion measurements and (13)C T(1) relaxation experiments quantified global reorientation using an anisotropic model together with interpretation of the internal dynamics via the "model-free" approach. Fitting of the dynamically averaged states to experimental data showed that the psi(2)(+)psi(3)(+) state is present to approximately 49%, psi(2)(-) psi(3)(-) to approximately 39%, and phi(3) (non-exo) to approximately 12%. Finally, using a dynamic and population-averaged model, (1)H,(1)H T-ROE buildup curves were calculated using a full relaxation matrix approach and were found to be in excellent agreement with experimental data, in particular for the above inter-residue proton-proton interaction between the terminal residues.     

The conformation of tetraalanine in water determined by polarized Raman, FT-IR, and VCD spectroscopy

The present article reports the conformation of cationic tetraalanine in aqueous solution. The determination of the dihedral angles of the two central amino acid residues was achieved by analyzing the amide I' band profile in the respective polarized visible Raman, Fourier transform-IR, and vibrational circular dichroism (VCD) spectra by means of a novel algorithm which utilizes the excitonic coupling between the amide I modes of nearest neighbor and second nearest peptide groups. It is an extension of a recently developed theory (Schweitzer-Stenner, R. Biophys. J., 2002, 83, 523-532). UV electronic circular dichroism (ECD) spectra of the peptides were used to validate the results of the structure analysis. The analyses yielded the dihedral angles (phi(12), psi(12)) = (-70 degrees, 155 degrees ) and (phi(23), psi(23)) = (-80 degrees, 145 degrees ). The obtained values are very close to the Ramachandran coordinates of the polyproline II helix (PPII). The data suggest that this is the conformation predominantly adopted by the peptide at room temperature. This notion was corroborated by the corresponding electronic circular dichroism spectrum. Tetraalanine exhibits a higher propensity for PPII than trialanine for which a 50:50 mixture of polyproline II and an extended beta-strand-like conformation was obtained from recent spectroscopic studies (Eker et al., J. Am. Chem. Soc. 2002, 124, 14330-14341). The temperature dependence of the CD spectra rule out that any cooperativity is involved in the strand if PPII transition. This led to the conclusion that solvent-peptide interactions give rise to the observed PPII stability. Our result can be utilized to understand why the denaturation of helix-forming peptides generally yields a PPII rather than a heterogeneous random conformation.     

Predictive in silico all-atom folding of a four-helix protein with a free-energy model

We report the predictive all-atom folding of the 60 amino acid four-helix bacterial ribosomal protein (BRP) L20 with a stochastic evolutionary optimization method in a free-energy force field. The energetically best, as well as six of the 10 lowest conformations, converge to a near-native structure. All of the 10 best energy conformations share the secondary structure elements of the native conformation, but differ in their tertiary alignment. The best conformation has a backbone root-mean-square deviation of 4.6 A to the native conformation and reproduces most distance constraints of the NMR experiment to 1.5 A resolution. Starting from random initial conditions, the native content of the simulated population increases more than 60-fold in the course of the simulation. These data demonstrate the feasibility of predictive unbiased all-atom protein folding with present day computational resources for the BRP L20.