Free energy based populations of interconverting microstates of a cyclic peptide lead to the experimental NMR data

Analysis of nuclear Overhauser enhancement (NOE) intensities data of interconverting microstates of a peptide is a difficult problem in nmr. A new statistical mechanics methodology has been proposed recently, consisting of several steps: (1) potential energy wells on the energy surface of the molecule are identified (the corresponding regions are called wide microstates); (2) each wide microstate is then spanned by a Monte Carlo (MC) or molecular dynamics simulation starting from a representative structure, and the corresponding relative populations are obtained from the free energy calculated with the local states method; and (3) the overall NOEs and 3J coupling constants are obtained as averages over the corresponding contributions of the samples, weighted by the populations. Extending this methodology to cyclic peptides, we are treating here the hexapeptide cyclo(D-Pro1-Phe2-Ala3-Ser4-Phe5-Phe6) in DMSO, which was studied by Kessler et al. using nmr (Journal of the American Chemical Society, 1992, Vol. 114, pp. 4805-4818). They found that at least two structures are required to explain their NOE data, a conclusion also corroborated by our analysis (Journal of the American Chemical Society, 1998, Vol. 120, pp. 800-812) and led to a novel derivation of atomic solvation parameters (ASPs) for DMSO. Thus, the overall interactions within the peptide-solvent system are described approximately by Etot = EGRO + summation operator sigmaiAi, where EGRO is the energy of the GROMOS force field, Ai is the solvent-accessible surface area of atom i, and sigmai is the ASP. In the present paper the validity of these ASPs within the framework of the entire methodology is verified. This requires taking into account 23 microstates. A very good agreement is obtained between experimental and calculated NOEs and 3J coupling constants. The free energy based populations lead to the best results, which means that entropic effects should not be ignored. We have also studied the behavior of the internal angular fluctuations of the proton-proton vectors and discovered that they have a negligible effect on the calculated NOEs; this is due to the relatively concentrated wide microstates spanned by the MC simulations. The applicability of our ASPs to other cyclic peptides in DMSO is being studied in another work and preliminary results are discussed.