All-atom Molecular Dynamics Simulationsand NMR Spectroscopy Study on Interactions and Structures in N-Glycylglycine Aqueous Solution

All-atom molecular dynamics (MD) simulation and the NMR spectra are used to investi-gate the interactions in N-glycylglycine aqueous solution. Different types of atoms exhibit different capability in forming hydrogen bonds by the radial distribution function analysis. Some typical dominant aggregates are found in different types of hydrogen bonds by the statistical hydrogen-bonding network. Moreover, temperature-dependent NMR are used to compare with the results of the MD simulations. The chemical shifts of the three hydrogen atoms all decrease with the temperature increasing which reveals that the hydrogen bonds are dominant in the glycylglycine aqueous solution. And the NMR results show agreement with the MD simulations. All-atom MD simulations and NMR spectra are successful in revealing the structures and interactions in the N-glycylglycine-water mixtures.     

Mode of binding of camptothecins to double helix oligonucleotides

We report an NMR study on the interaction of topotecan (Tpt) and other camptothecins (Cpts) with several double helix and single strand oligonucleotides. The results obtained by (31)P NMR spectroscopy, nuclear Overhauser experiments (NOE) and molecular dynamics (MD) simulations show that Cpt drugs do not intercalate into the double helix, as suggested by many authors. Phosphorus NMR spectra indicated that no deformation occurs at any level of the phosphodiester backbone, while 2D NOESY experiments allowed the detection of several contacts between the aromatic protons of Cpts and those of the double helix. Models of the drug/oligonucleotide complexes, built on the basis of NOE data, show that the drug is located at the end of the double helix, by stacking the A and B rings with the guanine or cytidine of the terminal CG base pairs, with a preference for the 3[prime or minute]-terminal end sites. Cpts interact with double strand, as well as with single strand oligomers, as can be seen from the NMR shift variation observed on the drug protons; but this shielding effect cannot be an evidence of intercalation, as it is largely due to external non-specific interactions of the positively charged drug with the negatively charged ionic surface of the oligonucleotide. The molecular weight of one of the complexes was obtained from the correlation time value. The conformational behaviour of the DNA fragment d(CGTACG)(2) was studied by MD simulations on a ns time scale in the presence of water molecules and Na(+) ions. Different models were examined and the deformations induced on the phosphodiester backbone by molecules that are known to intercalate, were monitored by MD simulations.     

Signatures of beta-peptide unfolding in two-dimensional vibrational echo spectroscopy: a simulation study.

An ensemble of exciton Hamiltonians for the amide-I band of the folded and unfolded states of a helical beta-heptapeptide is generated using a molecular dynamics (MD) simulation. The correlated fluctuations of its parameters and their signatures in two-dimensional (2D) vibrational echo spectroscopy are computed. This technique uses infrared pulse sequences to provide ultrafast snapshots of molecular structural fluctuations, in analogy with multidimensional NMR. The present study demonstrates that, by combining a method of calculating the vibrational Hamiltonian from MD snapshots and the nonlinear exciton equations (NEE), it may be possible to simulate realistic multidimensional IR spectra of chemically and biologically interesting systems.     

Ni(II).Arg-Gly-His-DNA interactions: investigation into the basis for minor-groove binding and recognition

A study of the minor-groove recognition of A/T-rich DNA sites by Ni(II).L-Arg-Gly-His and Ni(II).D-Arg-Gly-His was carried out with a fluorescence-based binding assay, one- and two-dimensional (1D and 2D) NMR methodologies, and molecular simulations. Fluorescence displacement titrations revealed that Ni(II).L-Arg-Gly-His binds to A/T-rich sequences better than the D-Arg diastereomer, while NMR investigations revealed that both metallopeptides bind to the minor groove of an AATT core region as evidenced by an intermolecular nuclear Overhauser effect (NOE) between each metallopeptide His imidazole C4 proton and the C2 proton of adenine. Results from molecular dynamics simulations of these systems were consistent with the experimental data and indicated that the His imidazole N-H, the N-terminal peptide amine, and Arg side chains of each metallopeptide are major determinants of minor-groove recognition by functioning as H-bond donors to the O2 of thymine residues or N3 of adenine residues.     

Conformational flexibility of the pentasaccharide LNF-2 deduced from NMR spectroscopy and molecular dynamics simulations

Human milk oligosaccharides (HMOs) are important as prebiotics since they stimulate the growth of beneficial bacteria in the intestine and act as receptor analogues that can inhibit the binding of pathogens. The conformation and dynamics of the HMO Lacto-N-fucopentaose 2 (LNF-2), α-L-Fucp-(1 → 4)[β-D-Galp-(1 → 3)]-β-D-GlcpNAc-(1 → 3)-β-D-Galp-(1 → 4)-D-Glcp, having a Lewis A epitope, has been investigated employing NMR spectroscopy and molecular dynamics (MD) computer simulations. 1D (1)H,(1)H-NOESY experiments were used to obtain proton-proton cross-relaxation rates from which effective distances were deduced and 2D J-HMBC and 1D long-range experiments were utilized to measure trans-glycosidic (3)J(CH) coupling constants. The MD simulations using the PARM22/SU01 force field for carbohydrates were carried out for 600 ns with explicit water as solvent which resulted in excellent sampling for flexible glycosidic torsion angles. In addition, in vacuo MD simulations were performed using an MM3-2000 force field, but the agreement was less satisfactory based on an analysis of heteronuclear trans-glycosidic coupling constants. LNF-2 has a conformationally well-defined region consisting of the terminal branched part of the pentasaccharide, i.e., the Lewis A epitope, and a flexible β-D-GlcpNAc-(1 → 3)-β-D-Galp-linkage towards the lactose unit, which is situated at the reducing end. For this β-(1 → 3)-linkage a negative ψ torsion angle is favored, when experimental NMR data is combined with the MD simulation in the analysis. In addition, flexibility on a similar time scale, i.e., on the order of the global overall molecular reorientation, may also be present for the ? torsion angle of the β-D-Galp-(1 → 4)-D-Glcp-linkage as suggested by the simulation. It was further observed from a temperature variation study that some (1)H NMR chemical shifts of LNF-2 were highly sensitive and this study indicates that Δδ/ΔT may be an additional tool for revealing conformational dynamics of oligosaccharides.     

Modelling a nematic liquid crystal

The bulk phase liquid crystalline behaviour of a cyclic siloxane with a pentamethylcyclosiloxane core and biphenyl-4-allyloxybenzoate mesogens (BCS) was studied using molecular dynamics (MD) and wide angle X-ray analysis. This material exhibits partial crystallinity at room temperature and liquid crystalline behaviour above 120° C. For the MD simulations an ensemble of 27 molecules with 135 mesogenic units was simulated and a molecular mechanics force field was used to model the structural anisotropy of the siloxane molecules. Simulations were carried out both at room temperature and at an elevated temperature (425 K). Room temperature simulations showed that, contrary to our initial assumptions, the low energy molecular conformations were not cylindrical but splayed in shape. During the simulation a smectic-like, tilted layer structure was found to evolve for the cluster when full atom potentials were used, while no such development was observed when electrostatic interactions were neglected. The presence of a tilted layered structure was also suggested by the X-ray data. These results indicate that long range electrostatic interactions are significant for the molecular system under study. In order to calculate the orientational order parameter, the orientation of the molecular axis had to be determined. This was achieved by describing the mesogen shapes to be ellipsoidal and defining the principal axis of the ellipsoids to be the molecular directors. By sampling over 200 ps of simulation at 425 K, the time averaged order parameter (S) was calculated. The calculated S of 0.36 was comparable to the value of 0.4-0.45 found from the experimental data. Apart from providing insight into the relative importance of the various competing forces in the formation of the liquid crystalline phase, these simulations are also expected to be useful in predicting the mesophase behaviour of liquid crystalline systems.     

Solution structure of a peptide nucleic acid duplex from NMR data: features and limitations

This paper describes the results of a 1D and 2D NMR spectroscopy study of a palindromic 8-base pair PNA duplex GGCATGCC in H2O and H2O-D2O solutions. The (1)H NMR peaks have been assigned for most of the protons of the six central base pairs, as well as for several amide protons of the backbone. The resulting 36 interbase and base-backbone distance restraints were used together with Watson-Crick restraints to generate the PNA duplex structure in the course of 10 independent simulated annealing runs followed by restrained molecular dynamics (MD) simulations in explicit water. The resulting PNA structures correspond to a P-type helix with helical parameters close to those observed in the crystal structures of PNA. Based on the current limited number of restraints obtained from NMR spectra, alternative structures obtained by MD from starting PNA models based on DNA cannot be ruled out and are also discussed.     

Characterization of Doubly Ionic Hydrogen Bonds in Protic Ionic Liquids by NMR Deuteron Quadrupole Coupling Constants: Differences to H‐bonds in Amides,Peptides, and Proteins

We present the first deuteron quadrupole coupling constants (DQCCs) for selected protic ionic liquids (PILs) measured by solid‐state NMR spectroscopy. The experimental data are supported by dispersion‐corrected density functional theory (DFT‐D3) calculations and molecular dynamics (MD) simulations. The DQCCs of the N−D bond in the triethylammonium cations are the lowest reported for deuterons in PILs, indicating strong hydrogen bonds between ions. The NMR coupling parameters are compared to those in amides, peptides, and proteins. The DQCCs show characteristic behavior with increasing interaction strength of the counterion and variation of the H‐bond motifs. We report the similar presence of the quadrupolar splitting pattern and the narrow liquid line in the NMR spectra over large temperature ranges, indicating the heterogeneous nature of PILs.     

Studies on the Structures and Hydrogen-bonding Interactions in Aqueous Glycine Solutions Using All-atom Molecular Dynamics Simulations and NMR Spectroscopy

All-atom molecular dynamics （MD） simulations and chemical shifts were used to study interactions and structures in the glycine-water system. Radial distribution functions and the hydrogen-bond network were applied in MD simulations. Aggregates in the aqueous glycine solution could be classified into different regions by analysis of the hydrogen-bonding network. Temperature-dependent NMR spectra and the viscosity of glycine in aqueous solutions were measured to compare with the results of MD simulations. The variation tendencies of the hydrogen atom chemical shifts and viscosity with concentration of glycine agree with the statistical results of hydrogen bonds in the MD simulations.     

Hyperbranched polyesters: synthesis, characterization, and molecular simulations

New types of hyperbranched polyesters were synthesized by the reaction of 2,2-bis(hydroxymethyl) propionic acid as an AB2-type monomer with pentaerythritol, trimethylol propane, or glycerol as the core moiety. The obtained globular networks were characterized by NMR and MALDI-TOF spectroscopic techniques. Molecular weights determined by MALDI-TOF were confirmed by gel permeation chromatography. Fourier transform infrared (FTIR) spectroscopy was used for the quantitative evaluation of hydrogen bonding as well as to study the structure-property relationship. To investigate the changes and types of intermolecular H-bonding interactions in hyperbranched polyesters with a variation in molecular structure, the deconvolution of FTIR spectra was carried out using Origin 6.0 software through the Gaussian curve-fitting method. Molecular simulations were performed through molecular mechanics and molecular dynamics (MD) calculations using the DISCOVER module. Cohesive energy density, solubility parameters, and surface properties of the hyperbranched polyesters were calculated. Further, vibrational analysis was computed using MD simulations for all the hyperbranched polyesters developed in this work.     

All-atom molecular dynamic simulations and relative NMR spectra study of weak C-H...O contacts in amide-water systems

Amide-water mixtures are studied by all-atom molecular dynamics (MD) simulations and the relative temperature-dependent NMR experiment. The weak C-H...O contacts are found in the amide-water systems theoretically and experimentally. The statistical results of the average numbers of hydrogen bonds indicate that the methyl groups in amide molecules represent different capabilities in forming the weak C-H...O contacts. The statistics also imply that the C-H...O contacts are more obvious in the amide-rich region than those in the water-rich region. The temperature-dependent NMR spectra are also adopted to investigate the weak C-H...O contacts in the amide-water systems. The relative chemical shifts of the methyl groups are in good agreement with the MD simulations.     

Study of the interaction of GFG tripeptide with cesium perfluorooctanoate micelles by means of NMR spectroscopy and MD simulations

The interaction of glycyl-phenylalanyl-glycine (GFG) with bilayers formed by cesium perfluorooctanoate (CsPFO) in water was investigated in the isotropic phase by means of 1H NMR and molecular dynamics (MD) simulations. Details on the preferential location of the different residues of GFG were obtained from selective variations of chemical shift with peptide concentration and of line width in the presence of the paramagnetic ion Mn2+. The analysis of 1H NMR spectra recorded at different concentrations and temperatures allowed the association constant and the enthalpy change upon binding to be evaluated. MD simulations highlighted the hydrogen bonds formed between the different GFG functional groups and the micelle. Both NMR and MD gave indications of high affinity of GFG with the micelle, with the N-terminal residue anchoring on the surface via hydrogen bonds with the micelle COO(-) groups.     

Comprehensive Characterization of the Self-Folding Cavitand Dynamics

The conformational equilibria and guest exchange process of a resorcin[4]arene derived self-folding cavitand receptor have been characterized in detail by molecular dynamics simulations (MD) and ¹H EXSY NMR experiments. A multi-timescale strategy for exploring the fluxional behaviour of this system has been constructed, exploiting conventional MD and accelerated MD (aMD) techniques. The use of aMD allows the reconstruction of the folding/unfolding process of the receptor by sampling high-energy barrier processes unattainable by conventional MD simulations. We obtained MD trajectories sampling events occurring at different timescales from ns to s: 1) rearrangement of the directional hydrogen bond seam stabilizing the receptor, 2) folding/unfolding of the structure transiting partially open intermediates, and 3) guest departure from different folding stages. Most remarkably, reweighing of the biased aMD simulations provided kinetic barriers that are in very good agreement with those determined experimentally by ¹H NMR. These results constitute the first comprehensive characterization of the complex dynamic features of cavitand receptors. Our approach emerges as a valuable rational design tool for synthetic host-guest systems     

Two-dimensional Raman spectra of atomic solids and liquids

We calculate third- and fifth-order Raman spectra of simple atoms interacting through a soft-core potential by means of molecular-dynamics (MD) simulations. The total polarizability of molecules is treated by the dipole-induced dipole model. Two- and three-body correlation functions of the polarizability at various temperatures are evaluated from equilibrium MD simulations based on a stability matrix formulation. To analyze the processes involved in the spectroscopic measurements, we divide the fifth-order response functions into symmetric and antisymmetric integrated response functions; the symmetric one is written as a simple three-body correlation function, while the antisymmetric one depends on a stability matrix. This analysis leads to a better understanding of the time scales and molecular motions that govern the two-dimensional (2D) signal. The 2D Raman spectra show novel differences between the solid and liquid phases, which are associated with the decay rates of coherent motions. On the other hand, these differences are not observed in the linear Raman spectra.     

Deuterium NMR studies of inositol liquid crystals [1]

Deuterium NMR investigations are presented on members of two new mesogenic series derived from the naturally occurring stereoisomers myo- and scW/o-inositol. Tetraethers of these two series exhibit thermotropic columnar phases in which the columns are apparently formed by stacked hydrogen bonded dimers of these molecules which chemically are vicinal diols. Deuterium NMR measurements were performed on the tetraoctyl homologues 2e (a cis diol) and 3e (a trans diol) of these series. We have investigated mixtures of these diols with small amounts of benzene-d₆ as probe molecules as well as samples of the neat diol compounds deuteriated at their hydroxyl groups. The results obtained show that the mesophases of both compounds are uniaxial and align partially in a magnetic field upon slow cooling from their isotropic liquids. The alignment is with the director parallel to the field direction indicating that the anisotropic magnetic susceptibility of this mesophase is positive. The deuterium quadrupole splitting of the benzene-d₆ probe in both systems is temperature dependent and in the trans diol 3e it even changes sign. This is interpreted in terms of a model in which the benzene-d₆ probe equilibrates rapidly between two (or more) solvation sites with quadrupole splittings of opposite signs The deuterium spectra of the neat deuterium labelled cis diol 2e exhibit two different signals due to the two deuterons which are located at the axial and equatorial hydroxyl groups. This indicates that there is no fast intra- or intermolecular exchange of the hydroxyl hydrogens. The overall quadrupole splittings of the hydroxyl deuterons in this compound are highly reduced compared to their static values and this is interpreted in terms of motional modes involving both reorientation of the hydroxyl deuterons about their C-O axis and overall reorientation of the molecules (or pairs of molecules) around the columnar axes. The corresponding spectra of the neat deuteriated trans diol 3e exhibit a single spectrum indicating that both hydroxyl deuterons in this compound are equivalent, or very nearly so. Within the mesophase region the spectrum undergoes gradual changes due to the increase in the molecular mobility, but the overall motional narrowing is less than in the cis isomer 2e. Apparently due to stronger hydrogen bonding in the trans isomer 3e the precession of the hydroxyl groups is hindered and a fast molecular reorientation is only possible at high temperatures.     

Force-field development and molecular dynamics simulations of ferrocene-peptide conjugates as a scaffold for hydrogenase mimics

The increasing importance of hydrogenase enzymes in the new energy research field has led us to examine the structure and dynamics of potential hydrogenase mimics, based on a ferrocene-peptide scaffold, using molecular dynamics (MD) simulations. To enable this MD study, a molecular mechanics force field for ferrocene-bearing peptides was developed and implemented in the CHARMM simulation package, thus extending the usefulness of the package into peptide-bioorganometallic chemistry. Using the automated frequency-matching method (AFMM), optimized intramolecular force-field parameters were generated through quantum chemical reference normal modes. The partial charges for ferrocene were derived by fitting point charges to quantum-chemically computed electrostatic potentials. The force field was tested against experimental X-ray crystal structures of dipeptide derivatives of ferrocene-1,1'-dicarboxylic acid. The calculations reproduce accurately the molecular geometries, including the characteristic C2-symmetrical intramolecular hydrogen-bonding pattern, that were stable over 0.1 micros MD simulations. The crystal packing properties of ferrocene-1-(D)alanine-(D)proline-1'-(D)alanine-(D)proline were also accurately reproduced. The lattice parameters of this crystal were conserved during a 0.1 micros MD simulation and match the experimental values almost exactly. Simulations of the peptides in dichloromethane are also in good agreement with experimental NMR and circular dichroism (CD) data in solution. The developed force field was used to perform MD simulations on novel, as yet unsynthesized peptide fragments that surround the active site of [Ni-Fe] hydrogenase. The results of this simulation lead us to propose an improved design for synthetic peptide-based hydrogenase models. The presented MD simulation results of metallocenes thereby provide a convincing validation of our proposal to use ferrocene-peptides as minimal enzyme mimics.     

Study of conformational properties of a biologically active peptide of fibronectin by circular dichroism, NMR and molecular dynamics simulation

Circular dichroism (CD), and NMR spectra have been recorded and molecular dynamics (MD) simulations have been performed in water and water-trifluoroethanol (TFE) mixed solvent for a synthetic biologically active 13-amino-acid fragment of human fibronectin and two related peptides. The CD results are interpreted on the basis of statistical analyses of MD trajectories and of ensuing calculations of CD spectra based on Schellman's matrix method. It is observed that the peptide conformation is quite variable in water and loses its mobility with the addition of TFE. (1)H-NOE data were found to be consistent with the most abundant calculated conformation.     

Structure elucidation and 3D solution conformation of the antibiotic enduracidin determined by NMR spectroscopy and molecular dynamics

Enduracidin and ramoplanin belong to the large family of cyclodepsipeptide antibiotics, highly effective against Gram-positive bacteria. The primary and 3D solution structure of ramoplanin is already well known, and the primary structure of enduracidin has been determined by a combination of chemical and NMR spectroscopic methods. Both antibiotics share a similar peptide core of 17 amino acids and differ mainly in the length of the acyl chain and the presence of two D-mannose moieties in ramoplanin. Based on the high sequence homology with ramoplanin, the structure in solution of enduracidin is modeled as a cyclic peptide. The tertiary structure thus obtained was refined through molecular dynamics (MD) simulation, in which the interatomic NOE-derived distance restraints were imposed. MD simulations yielded a family of representative 3D structures (RMSD = 0.89), which highlighted a backbone geometry similar to that of ramoplanin in its beta-hairpin arrangement. In contrast, enduracidin displays a different arrangement of the side-chain and of the residues forming the hydrophobic core.     

Motion of aromatic hydrocarbons in the microporous aluminum methylphosphonates AlMePO-alpha and AlMePO-beta

(2)H wide-line NMR has been used, in conjunction with molecular dynamics simulations where appropriate, to follow the reorientation of the monoaromatic compounds benzene, toluene, and p-xylene within the one-dimensional channels of the alpha- and beta-polymorphs of aluminum methylphosphonate, Al(2)(CH(3)PO(3))(3). Variable-temperature, static, (2)H NMR spectra of adsorbed d(6)-benzene, d(3)-, d(5)-, and d(8)-toluenes, and d(3),d(3)-p-xylene were matched by line shape simulation. The motion of p-xylene in both polymorphs is approximated by the long axis of the molecule describing a cone within the channels, the half-angle of which is greater for the slightly wider channels in AlMePO-beta (27-30 degrees cf. 18-19 degrees). The (2)H NMR of d(3)-toluene is simulated using a similar model, whereas the signal from aromatic deuterons in d(5)- and d(8)-toluenes is simulated by a ring undergoing 2pi/3 flips around the para axis. The reorientation of benzene shows the largest differences between the two pore structures. In AlMePO-beta it tumbles with little restriction, although at low temperatures the spectral details are better matched by allowing the molecule to spend a greater proportion of its time closer to the wall. In AlMePO-alpha the much broader line shape arises from constrained motion within the strongly triangular channels. Molecular dynamics simulations of benzene in the two structures confirm the differences. They support a model for benzene in AlMePO-alpha where its motion is restricted to rotations about its 6-fold axis and 2pi/3 jumps between symmetry-related sites in the pores, so that the plane of the aromatic ring remains approximately parallel to the c-axis.     

Conformations of Carnosine in Aqueous Solutions by All-Atom Molecular Dynamics Simulations and 2D-NOSEY Spectrum

All-atom molecular simulations and two-dimensional nuclear overhauser effect spectrum have been used to study the conformations of carnosine in aqueous solution. Intramolecular distances, root-mean-square deviation, radius of gyration, and solvent-accessible surface are used to characterize the properties of the carnosine. Carnosine can shift between extended and folded states, but exists mostly in extended state in water. Its preference for extension in pure water has been proven by the 2D nuclear magnetic resonance (NMR) experiment. The NMR experimental results are consistent with the molecular dynamics simulations.