Dependence of the L-alanyl-L-alanine conformation on molecular charge determined from ab initio computations and NMR spectra

The l-alanyl-l-alanine (AA) molecule behaves differently in acidic, neutral, and basic environments. Because of its molecular flexibility and strong interaction with the aqueous environment, its behavior has to be deduced from the NMR spectra indirectly, using statistical methods and comparison with ab initio predictions of geometric and spectral parameters. In this study, chemical shifts and indirect spin-spin coupling constants of the AA cation, anion, and zwitterion were measured and compared to values obtained by density functional computations for various conformers of the dipeptide. The accuracy and sensitivity of the quantum methods to the molecular charge was also tested on the (mono)-alanine molecule. Probable AA conformers could be identified at two-dimensional potential energy surfaces and verified by the comparison of the computed parameters with measured NMR data. The results indicate that, whereas the main-chain peptide conformations of the cationic (AA+) and zwitterionic (AAZW) forms are similar, the anion (AA-) adopts also another, approximately equally populated conformer in the aqueous solution. Additionally, the NH2 group can rotate in the two main chain conformations of the anionic form AA-. According to a vibrational quantum analysis of the two-dimensional energy surfaces, higher-energy conformers might exist for all three charged AA forms but cannot be detected directly by NMR spectroscopy because of their small populations and short lifetimes. In accord with previous studies, the NMR parameters, particularly the indirect nuclear spin-spin coupling constants, often provided an excellent probe of a local conformation. Generalization to peptides and proteins, however, has to take into account the environment, molecular charge, and flexibility of the peptide chain.     

Theoretical and Experimental Conformational Studies of Oligoglucosides Structurally Related to Fragments of Fungal Cell Wall β-(1→3)-D-Glucan

The conformational behavior of linear oligo-β-(1→3)-D-glucosides was studied using NMR experiments and molecular modeling. The explicit solvent model in calculations yielded the best coincidence between experimental and theoretical values of NOE and spin-spin coupling constants to evidence the strong influence of solvation upon the conformations of the oligoglucosides. Long-range coupling constants calculated for di- and trimeric clusters of the studied glucosides fit the experimental data much better than the single-molecule approach. It was shown that conformational properties of disaccharide fragments in studied oligoglucosides depended on neither their position in the chain nor the length of the chain.     

Calculation of structural behavior of indirect NMR spin-spin couplings in the backbone of nucleic acids

Calculated indirect NMR spin-spin coupling constants (J-couplings) between (31)P, (13)C, and (1)H nuclei were related to the backbone torsion angles of nucleic acids (NAs), and it was shown that J-couplings can facilitate accurate and reliable structural interpretation of NMR measurements and help to discriminate between their distinct conformational classes. A proposed stepwise procedure suggests assignment of the J-couplings to torsion angles from the sugar part to the phosphodiester link. Some J-couplings show multidimensional dependence on torsion angles, the most prominent of which is the effect of the sugar pucker. J-couplings were calculated in 16 distinct nucleic acid conformations, two principal double-helical DNAs, B- and A-, the main RNA form, A-RNA, as well as in 13 other RNA conformations. High-level quantum mechanics calculations used a baseless dinucleoside phosphate as a molecular model, and the effect of solvent was included. The predicted J-couplings correlate reliably with available experimental data from the literature.     

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.     

Calculating nuclear magnetic resonance chemical shifts in solvated systems

The nuclear magnetic resonance (NMR) chemical shift is extremely sensitive to molecular geometry, hydrogen bonding, solvent, temperature, pH, and concentration. Calculated magnetic shielding constants, converted to chemical shifts, can be valuable aids in NMR peak assignment and can also give detailed information about molecular geometry and intermolecular effects. Calculating chemical shifts in solution is complicated by the need to include solvent effects and conformational averaging. Here, we review the current state of NMR chemical shift calculations in solution, beginning with an introduction to the theory of calculating magnetic shielding in general, then covering methods for inclusion of solvent effects and conformational averaging, and finally discussing examples of applications using calculated chemical shifts to gain detailed structural information.     

Predicting the NMR spectra of nucleotides by DFT calculations: cyclic uridine monophosphate

We present an experimental and quantum chemical NMR study of the mononucleotide cyclic uridine monophosphate in water. Spectral parameters ((1)H and (13)C chemical shifts and (1)H--(1)H, (13)C--(1)H, (31)P--(13)C and (31)P--(1)H spin-spin coupling constants) have been carefully obtained experimentally and calculated using DFT methods including the solvent effect and the conformational flexibility of the solute. This study confirms that the (1)H and (13)C spectra of polar, flexible molecules in aqueous solution can be predicted with a high level of accuracy, comparable to that obtained for less complex systems.     

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.     

Non-empirical calculations of the nuclear magnetic resonance spectra of imidazole and hydrated imidazole

The magnetic shielding constants of the ¹H, ¹³C and ¹⁵N nuclei of imidazole are calculated for the isolated and hydrated molecules. The results show that the hydrogen bonds produce not only large variations of the chemical shifts for the nitrogen nuclei and the NH proton which are directly involved in the intermolecular bonding, but also measurable shifts for the carbon nuclei. The calculated shielding constants and their variation with hydration are discussed in relation to experimental results concerning imidazole, the 5-membered ring of the purine bases and the imidazole ring of histidine. The calculated values of the spin-spin coupling constants confirm that it is possible to study the tautomeric equilibrium of the imidazole ring from the measurement of these coupling constants and that spin-spin coupling constants are not very sensitive to solvent effects.     

Conformational study of a guaiacyl beta-O-4 lignin model compound by NMR. Examination of intramolecular hydrogen bonding interactions and conformational flexibility in solution

Intramolecular H-bonding interactions were investigated in solution for the threo and erythro diastereomeric forms of a guaiacyl beta-O-4 lignin model compound by using the NMR data obtained from hydroxyl protons. Temperature coefficients of the chemical shifts (ddelta/dT) and coupling constants (3J(HCOH)) were measured in aprotic and protic solutions: DMSO-d6, acetone-d6 and acetone-d6-water. The NMR parameters do not support the existence of strong and persistent intramolecular H-bonds that could participate in the stabilization of the guaiacyl beta-O-4 structure in solution, but instead indicate that intermolecular H-bonds to solvent predominate. 1D NOE experiments nevertheless revealed the presence of a direct chemical exchange between the hydroxyl protons, suggesting the possible existence of weak and transient intramolecular H-bonding interactions. The conformational flexibility of the threo structure was also investigated in acetone solution from the measurement of long-range 1H, 1H and 1H, 13C coupling constants and from NOESY experiments. The NMR data are not consistent with any single conformation, indicating that different conformers co-exist in solution. The experimental results support the conformational flexibility predicted by molecular dynamics simulations performed in a previous study. Finally, both experimental and theoretical approaches indicate that weak intramolecular H-bonds can exist transiently in solution, breaking and reforming as the beta-O-4 molecule undergoes conformational interconversion, but cannot be invoked as possible means of conferring rigidity to the beta-O-4 structure.     

Interpreting protein structural dynamics from NMR chemical shifts

In this investigation, semiempirical NMR chemical shift prediction methods are used to evaluate the dynamically averaged values of backbone chemical shifts obtained from unbiased molecular dynamics (MD) simulations of proteins. MD-averaged chemical shift predictions generally improve agreement with experimental values when compared to predictions made from static X-ray structures. Improved chemical shift predictions result from population-weighted sampling of multiple conformational states and from sampling smaller fluctuations within conformational basins. Improved chemical shift predictions also result from discrete changes to conformations observed in X-ray structures, which may result from crystal contacts, and are not always reflective of conformational dynamics in solution. Chemical shifts are sensitive reporters of fluctuations in backbone and side chain torsional angles, and averaged (1)H chemical shifts are particularly sensitive reporters of fluctuations in aromatic ring positions and geometries of hydrogen bonds. In addition, poor predictions of MD-averaged chemical shifts can identify spurious conformations and motions observed in MD simulations that may result from force field deficiencies or insufficient sampling and can also suggest subsets of conformational space that are more consistent with experimental data. These results suggest that the analysis of dynamically averaged NMR chemical shifts from MD simulations can serve as a powerful approach for characterizing protein motions in atomistic detail.     

Effect of the solvent on the 13C NMR spectra of pyrimidine nucleosides

The ¹³C NMR spectra of uridine and cytidine in d₆-DMSO, d₇-DMF, d₉-trimethyl phosphate, d₄-methanol, d₅-pyridine, and D₂O were investigated. A linear correlation between the C₁' chemical shifts and the J₁' vicinal spin-spin coupling constants of the protons was established. From the experimental data it may be assumed that the chief reason for the effect of the solvent on the C₁'₂' chemical shift is the different contribution of the effect of the base as a consequence of a change in the conformational equilibrium of the ribose ring. Deviations from the correlation in aqueous solutions and solutions of cytidine in pyridine are observed as a result of a change in the electron density in the base. The effect of the nature of the solvent on the position of the conformational equilibrium of the base relative to the glycoside bond was demonstrated.Translated from Khimiya Geterotsiklicheskikh Soedinenii, No. 11, pp. 1550–1552, November, 1980.     

Fluorine NMR spectra of pentafluorophenyl derivatives

The fluorine chemical shifts and spin-spin coupling constants of 65 pentafluorophenyl derivatives with widely varying organic substituents were examined. Useful correlations of the three meta coupling constants with the chemical shifts of the para fluorine were found. It is suggested that these relationships be extended to all compounds of the type considered in order to determine the signs and approximate values of meta coupling constants. Equations for correlation of the fluorine chemical shifts with the Taft constants are presented. The possibility of calculating the Taft constants from the ¹⁹F NMR spectra of pentafluorophenyl compounds is being discussed.     

Between Two Chairs: Combination of Theory and Experiment for the Determination of the Conformational Dynamics of Xylosides

The conformational properties of monosaccharides constitute fundamental features of oligosaccharides. While the energy landscape of monosaccharides can be altered by a specific biochemical environment or by chemical modifications, the analysis of resulting dynamic conformational equilibria is not feasible by experimental means alone. In this work, a series of β-d -xylopyranosides is used to outline how a combination of experimental NMR parameters and computed molecular properties can be used to determine conformers and quantify the composition of conformational equilibria. We demonstrate that identifying the most stable conformers using energy calculations is challenging and computing of NMR shieldings is typically not sensitive enough. On the other hand, computed spin-spin coupling constants for the xyloside ring can be used to unambiguously assign experimental NMR data of dynamic conformational equilibria and quantify the ratio of different conformers in the mixture. As a proof of principle, this procedure allowed to analyze a hitherto unknown dynamic equilibrium of a diamino-xyloside as a precursor of a molecular switch.     

Aromatic C-H...O interactions in a series of bindone analogues. NMR and quantum mechanical study

The aromatic C-H...O hydrogen bonding within the series of the structurally relative indenone derivatives has been studied. The presence of the hydrogen bonds is corroborated by the large low-field chemical shifts of the protons involved in the hydrogen bond observed experimentally and reproduced by quantum mechanical calculations. Further confirmation is provided by analysis of the orbital overlap coefficients, (13)C NMR chemical shifts, and one-bond spin-spin coupling constants J((13)C-(1)H). The relationship between molecular geometry and (1)H NMR chemical shifts of involved protons has a complex nature, but the C-H...O distance is the principal factor.     

Intramolecular hydrogen bonding in disubstituted ethanes: general considerations and methodology in quantum mechanical calculations of the conformational equilibria of succinamate monoanion

The importance of intramolecular hydrogen bonding between the carboxylate oxygen and amide proton of succinamate anion has been investigated by quantum mechanical simulations as a function of solvent for comparison with conformational equilibria estimated by NMR spectroscopy. The focus is on those methodological considerations of general interest to the conformational equilibrium problem, which are also particularly relevant to the quantum calculations. The roughly planar symmetry of the amide and carboxylate substituents of succinamate anion and the possibility of intramolecular hydrogen-bond formation together suggest that the orientational degrees of freedom of the substituents could play an important role in the equilibrium of the CH2-CH2 torsion. To test this hypothesis, two-dimensional potential-energy surfaces (PESs) were mapped out from the quantum mechanical calculations, with coordinates corresponding to the CH2-CH2 torsional and amide group rotational degrees of freedom. The Boltzmann populations obtained from two-dimensional PESs and those obtained from a one-dimensional adiabatic surface for the CH2-CH2 torsion were compared with the experimental results. In these comparisons, the agreement of calculated gauche fractions with corresponding experimental values was checked, as well as the agreement between predicted coupling constants and those determined from experimental spectra. In polar protic and aprotic solvents, where highly polar trans conformations can be stabilized by dipole-dipole and hydrogen-bonding interactions with the solvent, the orientational degree of freedom of the amide substituent appears to play a sufficiently important role in the CH2-CH2 torsional equilibrium that it cannot be safely ignored in the simulations.     

Structural and thermodynamic analysis of heteroassociation of daunomycin and flavin mononucleotide molecules in water by <Superscript>1</Superscript>H NMR spectroscopy

Heteroassociation of the antitumor antibiotic daunomycin (DAU) with flavin mononucleotide (FMN) has been investigated by one-and two-dimensional ¹H NMR spectroscopy (500 MHz) in a water solution to determine the molecular mechanism of the combined action of the antibiotic and vitamin in the FMN-DAU system. The equilibrium constants of the reactions, induced proton chemical shifts, and thermodynamic parameters (ΔH, ΔS) of heteroassociation were determined from the concentration and temperature dependences of the proton chemical shifts in the interacting aromatic molecules. Analysis of the results indicate that heterocomplexes of riboflavin mononucleotide and daunomycin are formed due to stacking interactions between aromatic chromophores. The most probable spatial structure of the 1:1 DAU-FMN heterocomplex was determined by the molecular dynamics method using the X-PLOR program and the results of the analysis of the induced proton chemical shifts in molecules. Calculation of the relative content of self-and hetero-complexes of daunomycin for different values of the ratio (r) between the concentrations of flavin mononucleotide and daunomycin demonstrated that for r > 3, the contribution of DAU-FMN heterocomplexes to the equilibrium distribution of associates in aqueous solution is dominant. It is concluded that the aromatic molecules of vitamins, in particular, riboflavin, can form energetically strong heteroassociates with antitumor antibiotics in water solution and can thereby affect their medical and biological activity.     

Complex Formation in the Ternary System Tl(III)-CN(-)-Cl(-) in Aqueous Solution. A (205)Tl NMR Study

The existence of mixed complexes of the general formula Tl(CN)(m)()Cl(n)()(3)(-)(m)()(-)(n)() (m + n 相似文献   

1H and 13C NMR chemical shifts and spin-spin coupling constants in trans- and cis-decalins

The NMR parameters characterizing the spectra of trans- and cis-decalins were determined from theoretical calculations and experimental spectra. The calculated values of the shielding constants are in good agreement with the measured chemical shifts, with a small but noticeable difference in accuracy for the bridgehead atoms. Of all the spin-spin coupling constants, only most of (1)J(C,C) and (1)J(C,H) values could be extracted from the spectra, and the corresponding computed values are in good agreement with experiment. It appears that the applied density functional theory (DFT) approach overestimates slightly the J(C,C) coupling and underestimates the differences between one-bond (1)J(C,H) coupling constants. For all these constants [J(C,C), J(C,H) and J(H,H)] through one to three bonds, which could not be obtained experimentally, the predicted values are in good agreement with the general rules relating spin-spin coupling to the number and spatial arrangement of the intervening bonds.     

Car-Parrinello molecular dynamics simulations and EPR property calculations on aqueous ubisemiquinone radical anion

Car-Parrinello molecular dynamics (CPMD) simulations have been performed on ubisemiquinone radical anion in aqueous solution. The different types of hydrogen bonding formed between the semiquinone and the solvent were studied in terms of frequency and directionality, in comparison with the parent benzosemiquinone radical anion. The EPR parameters (g-tensors and hyperfine coupling constants) were obtained from quantum chemical property calculations performed on snapshots along the MD trajectory, and the contributions of different solvation effects to the EPR parameters have been evaluated. The influence of the anion’s conformational behaviour was examined, including the orientation-dependent effects of hyperconjugation on side-chain hyperfine coupling.