Conformations of Oxidized Glutathione in Aqueous Urea Solution by All-Atom Molecular Dynamic Simulations and 2D-NOESY Spectrum

All-atom molecular simulations and two-dimensional nuclear overhauser effect spectra have been used to study the conformations and interactions of oxidized glutathione (GSSG) in aqueous urea solution. The simulations were characterized by intramolecular distance, radius of gyration, solvent-accessible surface area, and root-mean-square deviation. Interestingly, the two chains connected by the GSSG disulfide linkage exhibited different flexibilities in the aqueous urea solution. GSSG can convert from “extended” to “folded” states in the simulations. Its preferred conformation in aqueous urea solutions is “extended”, which was confirmed by the 2D nuclear magnetic resonance (NMR) experiment. The two different types of amide hydrogen atoms in cysteine and glycine also showed different capacities to form N–H?O hydrogen bonds. The results were confirmed by temperature-dependent NMR experiment.     

Molecular Dynamics Simulations and NMR Experimental Study of Oxidized Glutathione in Aqueous Solution

All-atom molecular simulations and NMR experiments have been used to study the conformations and interactions of oxidized glutathione (GSSG) in aqueous solution. The simulations are characterized by the radius of gyration, intramolecular distance, root-mean-square deviation and solvent-accessible surface area. The variations in these properties show time dependences. Interestingly, the two chains connected by the disulfide linkage in GSSG show different flexibilities in aqueous solution. The conformations of GSSG can convert from ??extended?? to ??folded?? states. Also, the two different kinds of amide hydrogen atoms in cysteine (Cys) and glycin (Gly) show different capabilities in forming N?CH?O hydrogen bonds with water molecules. Temperature-dependent NMR results show agreements with the MD simulations.     

Multiexponential electronic spin relaxation and Redfield's limit in Gd(III) complexes in solution: consequences for 17O/1H NMR and EPR simultaneous analysis

Multiple experiments (17O NMR, 1H NMR, and EPR) have been performed in the past to understand the microscopic parameters that control the magnetic relaxation rate enhancement induced by paramagnetic molecules on neighboring water protons, the so-called relaxivity. The generally accepted theories of the electron spin relaxation of S = 7/2 ions such as Gd3+ (Solomon-Bloembergen-Morgan or simplified Hudson-Lewis) are unsatisfactory for a simultaneous analysis. Recently, an improved theory, where the electron spin relaxation is due to the combination of a static (thus explicitly linked to the molecular structure) and a dynamic zero field splitting, has been developed and tested on experimental EPR data. The model has also been extended beyond the electronic Redfield limit using Monte Carlo simulations. Using the aqua ion [Gd(H2O)8]3+ as a test case, we present here the first simultaneous analysis of 17O NMR, 1H NMR, and EPR relaxation data using this rigorous approach of the electron spin relaxation. We discuss the physical meaning of the calculated parameters. The consequences on future experiments are also considered, especially regarding the analysis of nuclear magnetic relaxation dispersion (NMRD) profiles in the study of Gd3+ complexes.     

NMR solvent shifts of acetonitrile from frozen density embedding calculations

We present a density functional theory (DFT) study of solvent effects on nuclear magnetic shielding parameters. As a test example we have focused on the sensitive nitrogen shift of acetonitrile immersed in a selected set of solvents, namely water, chloroform, and cyclohexane. To include the effect of the solvent environment in an accurate and efficient manner, we employed the frozen-density embedding (FDE) scheme. We have included up to 500 solvent molecules in the NMR computations and obtained the cluster geometries from a large set of conformations generated with molecular dynamics. For small solute-solvent clusters comparison of the FDE results with conventional supermolecular DFT calculations shows close agreement. For the large solute-solvent clusters the solvent shift values are compared with experimental data and with values obtained using continuum solvent models. For the water --> cyclohexane shift the obtained value is in very good agreement with experiments. For the water --> chloroform NMR solvent shift the classical force field used in the molecular dynamics simulations is found to introduce an error. This error can be largely avoided by using geometries taken from Car-Parrinello molecular dynamics simulations.     

Use of NMR spectroscopy in the analysis of carnosine and free amino acids in fermented sausages during ripening

Quantitative changes of carnosine and free amino acids in high-fat (43–50 mass %) fermented sausages during ripening were analysed using a 600 MHz VNMRS NMR spectrometer. Seven free amino acids were identified in the samples and a relatively high content of carnosine was observed in the final stage of ripening. The NMR method for the determination of free amino acids and carnosine content applied in this work has been used for the first time and it has proven to be suitable for the analysis of fermented sausages.     

Radio frequency (RF-) saturation spectroscopy as detected via optical nuclear polarization,mechanisms due to ESR and NMR transitions in excited triplet states

The rate equation treatment of Optical Nuclear Polarization (ONP) by an optical excitation cycle passing through an excited molecular triplet state has been extended to include rf-saturation of ESR and NMR transitions in the excited state. Specifically, the ONP changes due to allowed and forbidden ESR transitions are described and compared with each other. NMR transitions are considered in connection with double rf-saturation when nuclear polarization is generated first by saturation of ESR frequencies. The relevant predictions of the model calculations are summarized and used to explain some recent experimental results.     

Conformational polymorphism of the PrP106-126 peptide in different environments: a molecular dynamics study

Extensive molecular dynamic simulations (approximately 240 ns) have been used to investigate the conformational behavior of PrP106-126 prion peptide in four different environments (water, dimethyl sulfoxide, hexane, and trifluoroethanol) and under both neutral and acidic conditions. The conformational polymorphism of PrP106-126 in solution observed in the simulations supports the role of this fragment in the structural transition of the native to the abnormal form of prion protein in response to changes in the local environmental conditions. The peptide in solution is primarily unstructured. The simulations show an increased presence of helical structure in an apolar solvent, in agreement with the results from circular dichroism spectroscopy. In water solution, beta-sheet elements were observed between residues 108-112 and either residues 115-121 or 121-126. An alpha-beta transition was observed under neutral conditions. In DMSO, the peptide adopted an extended conformation, in agreement with nuclear magnetic resonance experiments.     

NMR chemical shifts as a tool to analyze first principles molecular dynamics simulations in condensed phases: the case of liquid water

We present (1)H NMR chemical shift calculations of liquid water based on first principles molecular dynamics simulations under periodic boundary conditions. We focus on the impact of computational parameters on the structural and spectroscopic data, which is an important question for understanding how sensitive the computed (1)H NMR resonances are upon variation of the simulation setup. In particular, we discuss the influence of the exchange-correlation functional and the size of the basis set, the choice for the fictitious electronic mass and the use of pseudopotentials for the nuclear magnetic resonance (NMR) calculation on one hand and the underlying Car-Parrinello-type molecular dynamics simulations on the other hand. Our findings show that the direct effect of these parameters on (1)H shifts is not big, whereas the indirect dependence via the structural data is more important. The (1)H NMR chemical shifts clearly reflect the induced structural changes, illustrating once again the sensitivity of (1)H NMR observables on small changes in the local chemical structure of complex hydrogen-bonded liquids.     

Studies on the Conformations and Hydrogen Bonding of ACE Inhibitory Tripeptide VEF by All-atom Molecular Dynamics Simulations and Molecular Docking

Molecular dynamic simulations and molecular docking are performed to study the conformations and hydrogen bonding interactions of ACE inhibitory tripeptide VEF. Intramolecular distance, radius of gyration, solvent-accessible surface, and root-mean-square deviations are used to characterize the properties of VEF in aqueous solution. The VEF molecule is highly flexible in water and conformations can shift between the extended and folded states. The VEF molecule exists in extended state mostly in aqueous solution and the conformations bonded with ACE are also the extended ones. The findings indicate that MD simulations have a good agreement with the molecular docking analysis.     

Computing the 7Li NMR chemical shielding of hydrated Li+ using cluster calculations and time-averaged configurations from ab initio molecular dynamics simulations

Ab initio molecular dynamics (AIMD) simulations have been used to predict the time-averaged Li NMR chemical shielding for a Li(+) solution. These results are compared to NMR shielding calculations on smaller Li(+)(H(2)O)(n) clusters optimized in either the gas phase or with a polarizable continuum model (PCM) solvent. The trends introduced by the PCM solvent are described and compared to the time-averaged chemical shielding observed in the AIMD simulations where large explicit water clusters hydrating the Li(+) are employed. Different inner- and outer-coordination sphere contributions to the Li NMR shielding are evaluated and discussed. It is demonstrated an implicit PCM solvent is not sufficient to correctly model the Li shielding, and that explicit inner hydration sphere waters are required during the NMR calculations. It is also shown that for hydrated Li(+), the time averaged chemical shielding cannot be simply described by the population-weighted average of coordination environments containing different number of waters.     

Different Behaviors of Glutathione in Aqueous and DMSO Solutions: Molecular Dynamics Simulation and NMR Experimental Study

An all-atom molecular simulation and NMR experiments have been carried out to investigate the interactions and conformations of glutathione (GSH) in aqueous and DMSO solutions. The simulations started, from different initial conformations, are characterized by intramolecular distance, radius of gyration, root-mean-square deviation, and solvent-accessible surface. Interestingly, different behaviors are found in the two different solutions. GSH is highly flexible in an aqueous solution with transitions to the extended, semifolded, and folded states. However, once GSH reaches the folded state in DMSO, it remains there and becomes difficult to break down. The NMR results show agreement with the MD simulations. The water molecule is small. It is also a good proton donor and a good proton acceptor. Water molecules can easily break down the “folded” conformation. In DMSO solution, the stronger hydrogen bonds and the hydrophobic interactions are more important, which can make the GSH in the folded state stable. Variations in the distribution of conformations and the hydrogen-bonding network may play an important role in its function under physiological conditions.     

Structure of human telomeric DNA in crowded solution

G-quadruplex structures formed by DNA at the human telomeres are attractive anticancer targets. Human telomeric sequences can adopt a diverse range of intramolecular G-quadruplex conformations: a parallel-stranded conformation was observed in the crystalline state, while at least four other forms were seen in K(+) solution, raising the question of which conformation is favored in crowded cellular environment. Here, we report the first NMR structure of a human telomeric G-quadruplex in crowded solution. We show that four different G-quadruplex conformations are converted to a propeller-type parallel-stranded G-quadruplex in K(+)-containing crowded solution due to water depletion. This study also reveals the formation of a new higher-order G-quadruplex structure under molecular crowding conditions. Our molecular dynamics simulations of solvent distribution provide insights at molecular level on the formation of parallel-stranded G-quadruplex in environment depleted of water. These results regarding human telomeric DNA can be extended to oncogenic promoters and other genomic G-rich sequences.     

Anisotropic NMR interaction tensors in the decamethylaluminocenium cation

Solid-state NMR experiments, analytical and numerical simulations of solid NMR powder patterns, ab initio self-consistent field and hybrid density functional theory calculations, and single-crystal X-ray diffraction are used to characterize the molecular structure and anisotropic NMR interaction tensors in the bis(pentamethylcyclopentadienyl)aluminum cation, [Cp(2)Al](+). This highly symmetric main group metallocene has a structure analogous to that of ferrocene and the cobaltocenium cation. The single-crystal X-ray diffraction structure is reported for [Cp(2)Al][AlCl(4)]. Solid-state (27)Al[(1)H] magic-angle spinning and static NMR experiments are used to study the aluminum chemical shielding and electric field gradient tensors, revealing axial symmetry in both cases with a large chemical shielding span of Omega = 83(3) ppm and a small nuclear quadrupole coupling constant, C(Q)((27)Al) = 0.86(10) MHz. Carbon-13 CPMAS NMR experiments in combination with ab initio calculations and simulations of the effects of chemical exchange on (13)C static powder patterns reveal dynamic rotation of rings and suggest a low internal rotational barrier for this process. Theoretical computations of interaction tensors using the Gaussian 98 and Amsterdam Density Functional software packages are in good agreement with experiment and lend insight into the molecular origin of these NMR interactions. Orientations of the NMR tensors determined from theory, the large chemical shielding span, and the very small value of C(Q)((27)Al) can all be rationalized in terms of the high molecular symmetry.     

Modelling zwitterions in solution: 3-fluoro-γ-aminobutyric acid (3F-GABA)

The conformations and relative stabilities of folded and extended 3-fluoro-γ-aminobutyric acid (3F-GABA) conformers were studied using explicit solvation models. Geometry optimisations in the gas phase with one or two explicit water molecules favour folded and neutral structures containing intramolecular NH···O-C hydrogen bonds. With three or five explicit water molecules zwitterionic minima are obtained, with folded structures being preferred over extended conformers. The stability of folded versus extended zwitterionic conformers increases on going from a PCM continuum solvation model to the microsolvated complexes, though extended structures become less disfavoured with the inclusion of more water molecules. Full explicit solvation was studied with a hybrid quantum-mechanical/molecular-mechanical (QM/MM) scheme and molecular dynamics simulations, including more than 6000 TIP3P water molecules. According to free energies obtained from thermodynamic integration at the PM3/MM level and corrected for B3LYP/MM total energies, the fully extended conformer is more stable than folded ones by about -4.5 kJ mol(-1). B3LYP-computed (3)J(F,H) NMR spin-spin coupling constants, averaged over PM3/MM-MD trajectories, agree best with experiment for this fully extended form, in accordance with the original NMR analysis. The seeming discrepancy between static PCM calculations and experiment noted previously is now resolved. That the inexpensive semiempirical PM3 method performs so well for this archetypical zwitterion is encouraging for further QM/MM studies of biomolecular systems.     

Prediction of the 1H and 13C NMR spectra of alpha-D-glucose in water by DFT methods and MD simulations

We have applied computational protocols based on DFT and molecular dynamics simulations to the prediction of the alkyl 1H and 13C chemical shifts of alpha-d-glucose in water. Computed data have been compared with accurate experimental chemical shifts obtained in our laboratory. 13C chemical shifts do not show a marked solvent effect. In contrast, the results for 1H chemical shifts provided by structures optimized in the gas phase are only fair and point out that it is necessary to take into account both the flexibility of the glucose structure and the strong effect exerted by solvent water thereupon. Thus, molecular dynamics simulations were carried out to model both the internal geometry as well as the influence of solvent molecules on the conformational distribution of the solute. Snapshots from the simulation were used as input to DFT NMR calculations with varying degrees of sophistication. The most important factor that affects the accuracy of computed 1H chemical shifts is the solute geometry; the effect of the solvent on the shielding constants can be reasonably accounted for by self-consistent reaction field models without the need of explicitly including solvent molecules in the NMR property calculation.     

Dynamical behavior of one-dimensional water molecule chains in zeolites: nanosecond time-scale molecular dynamics simulations of bikitaite

Nanosecond scale molecular dynamics simulations of the behavior of the one-dimensional water molecule chains adsorbed in the parallel nanochannels of bikitaite, a rare lithium containing zeolite, were performed at different temperatures and for the fully and partially hydrated material. New empirical potential functions have been developed for representing lithium-water interactions. The structure and the vibrational spectrum of bikitaite were in agreement both with experimental data and Car-Parrinello molecular dynamics results. Classical molecular dynamics simulations were extended to the nanosecond time scale in order to study the flip motion of water molecules around the hydrogen bonds connecting adjacent molecules in the chains, which has been observed by NMR experiments, and the dehydration mechanism at high temperature. Computed relaxation times of the flip motion follow the Arrhenius behavior found experimentally, but the activation energy of the simulated system is slightly underestimated. Based on the results of the simulations, it may be suggested that the dehydration proceeds by a defect-driven stepwise diffusion. The diffusive mechanism appears as a single-file motion: the molecules never pass one another, even at temperatures as high as about 1000 K, nor can they switch between different channels. However, the mean square displacement (MSD) of the molecules, computed with respect to the center of mass of the simulated system, shows an irregular trend from which the single-file diffusion cannot be clearly evidenced. If the MSDs are evaluated with respect to the center of mass of the molecules hosted in each channel, the expected dependence on the square root of time finally appears.     

The low-temperature dynamic crossover phenomenon in protein hydration water: simulations vs experiments

A super-Arrhenius-to-Arrhenius dynamic crossover phenomenon has been observed in the translational alpha-relaxation time and in the inverse of the self-diffusion constant both experimentally and by simulations for lysozyme hydration water in the temperature range of TL = 223 +/- 2 K. MD simulations are based on a realistic hydrated powder model, which uses the TIP4P-Ew rigid molecular model for the hydration water. The convergence of neutron scattering, nuclear magnetic resonance and molecular dynamics simulations supports the interpretation that this crossover is a result of the gradual evolution of the structure of hydration water from a high-density liquid to a low-density liquid form upon crossing of the Widom line above the possible liquid-liquid critical point of water.     

Molecular dynamics simulations of cyclosporin A: The crystal structure and dynamic modelling of a structure in apolar solution based on NMR data

Summary The conformation of the immunosuppressive drug cyclosporin A (CPA), both in apolar solution and in crystalline state, has been studied by computer simulation techniques. Three molecular dynamics (MD) simulations have been performed: one modelling the crystal structure and two modelling the structure in apolar solution, using a restrained MD approach in which data from nuclear magnetic resonance (NMR) and infrared (IR) spectroscopy are taken into account. The simulation of the crystalline state (MDC) concerns a system of 4 unit cells containing 16 cyclosporin A molecules and 22 water molecules, which is simulated using crystalline periodic boundary conditions. The simulations modelling the apolar solvent conformation (MDS) concern one isolated cyclosporin A molecule. In these simulations an extra term in the interatomic potential function is used, which forces the molecule to satisfy a set of 57 atom-atom distance constraints originating from nuclear Overhauser effects (NOEs) obtained from NMR spectroscopy and one distance constraint deduced from IR spectroscopy.From a comparison of the results of the crystal simulation to those of the X-ray experiment in terms of structure, atomic fluctuations, hydrogen bond pattern, etc., it is concluded that the force field that is used yields an adequate representation of crystalline cyclosporin A. Secondly, it is shown that the dynamic modelling technique that is used to obtain a structure in a polar solution from NMR distance information works well. Starting from initial conformations which have a root mean square difference of 0.14 nm both distance restrained MD simulations converge to the same final solution structure. A comparison of the crystal structure of cyclosporin A and the one in apolar solution shows that there are significant differences. The overall difference in atomic positions is 0.09 nm for the C_x atoms and 0.17 nm for all atoms. In apolar solution, the molecule is slightly more bent and the side chains of 1 MeBmt and 10 MeLeu adopt a different conformation.Abbreviations MeBmt (4R)-4[(E)-2-butenyl]-4-methyl-l-Threonine - MD Molecular dynamics - EM Energy minimization - MDC Molecular dynamics simulation of the crystal - MDS1 Restrained molecular dynamics simulation to obtain the structure in solution starting from the crystal structure - MDS2 Like MDS1, but starting from the SMS structure - SMS Proposed structure in solution, obtained by model building - XRAY An X-ray structure - CPA Cyclosporin A - NMR Nuclear magnetic resonance spectroscopy - NOE Nuclear Overhauser enhancement - MDS1 Mean simulated structure obtained by averaging over the time period 20–40 ps of the MDS1 simulation - MDS2 Mean simulated structure obtained by averaging over the time period 10–30 ps of the MDS2 simulation - Mean simulated structure obtained by averaging over the time period 7–15 ps and over the 16 asymmetric units in the computational box of the MDC simulation.     

Towards Accurate Predictions of Proton NMR Spectroscopic Parameters in Molecular Solids

The factors contributing to the accuracy of quantum-chemical calculations for the prediction of proton NMR chemical shifts in molecular solids are systematically investigated. Proton chemical shifts of six solid amino acids with hydrogen atoms in various bonding environments (CH, CH₂, CH₃, OH, SH and NH₃) were determined experimentally using ultra-fast magic-angle spinning and proton-detected 2D NMR experiments. The standard DFT method commonly used for the calculations of NMR parameters of solids is shown to provide chemical shifts that deviate from experiment by up to 1.5 ppm. The effects of the computational level (hybrid DFT functional, coupled-cluster calculation, inclusion of relativistic spin-orbit coupling) are thoroughly discussed. The effect of molecular dynamics and nuclear quantum effects are investigated using path-integral molecular dynamics (PIMD) simulations. It is demonstrated that the accuracy of the calculated proton chemical shifts is significantly better when these effects are included in the calculations.     

Refinement of the NMR structures of alpha-conotoxin MI using molecular dynamics simulation with explicit solvent water and a full molecular force field

Three NMR structures of alpha-conotoxin MI, a potent antagonist of the nicotinic acetylcholine receptor, have been refined using molecular dynamics (MD) simulation with explicit water. Although the convergence of the NMR structures of alpha-conotoxin MI was not sufficient to provide detailed structural features, the average structures obtained from MD simulations converged to one conformation, providing structural characteristics. The resulting structure was also found to be consistent with the results of amide proton-exchange experiments. These results demonstrate that MD simulation with explicit solvent water is very useful in refining NMR structures.