The structure of aqueous L-proline amino acid has been the subject of much debate centering on the validity of various proposed models, differing widely in the extent to which local and long-range correlations are present. Here, aqueous proline is investigated by atomistic, replica exchange molecular dynamics simulations, and the results are compared to neutron diffraction and small angle neutron scattering (SANS) data, which have been reported recently (McLain, S.; Soper, A.; Terry, A.; Watts, A. J. Phys. Chem. B 2007, 111, 4568). Comparisons between neutron experiments and simulation are made via the static structure factor S(Q) which is measured and computed from several systems with different H/D isotopic compositions at a concentration of 1:20 molar ratio. Several different empirical water models (TIP3P, TIP4P, and SPC/E) in conjunction with the CHARMM22 force field are investigated. Agreement between experiment and simulation is reasonably good across the entire Q range although there are significant model-dependent variations in some cases. In general, agreement is improved slightly upon application of approximate quantum corrections obtained from gas-phase path integral simulations. Dimers and short oligomeric chains formed by hydrogen bonds (frequently bifurcated) coexist with apolar (hydrophobic) contacts. These emerge as the dominant local motifs in the mixture. Evidence for long-range association is more equivocal: No long-range structures form spontaneously in the MD simulations, and no obvious low-Q signature is seen in the SANS data. Moreover, associations introduced artificially to replicate a long-standing proposed mesoscale structure for proline correlations as an initial condition are annealed out by parallel tempering MD simulations. However, some small residual aggregates do remain, implying a greater degree of long-range order than is apparent in the SANS data.

The structure of a short fragment of the human HIV-1 membrane glycoprotein gp41 has been examined using a combination of parallel tempering molecular dynamics (PTMD) and far UV circular dichroism spectroscopy. The aim is to resolve conflicting reports on the solution state conformational bias in this membrane proximal domain spanning the epitope for the 2F5 monoclonal antibody. We conclude that gp41(659-671) exhibits conformational plasticity in which competing folding propensities are present and can be influenced by local microenvironment. Contrary to previous reports, the 3(10) helix does not emerge as a dominant motif from either simulation or experiment, and this peptide is therefore not a model system for this fold type. Other fold groups such as turn motifs are identifiable at elevated temperatures in the PTMD trajectories and are potentially relevant in antibody binding. Helical populations in pure water are significantly overestimated according to the CHARMM parametrization. However, circular dichroism (CD) data show that helices are promoted in membrane mimetic solvents. As this is a membrane proximal peptide, the helical motif may well have physiological significance.