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.     

Interpreting NMR data for beta-peptides using molecular dynamics simulations

NMR is one of the most used techniques to resolve structure of proteins and peptides in solution. However, inconsistencies may occur due to the fact that a polypeptide may adopt more than one conformation. Since the NOE distance bounds and (3)J-values used in such structure determination represent a nonlinear average over the total ensemble of conformers, imposition of NOE or (3)J-value restraints to obtain one unique conformation is not an appropriate procedure in such cases. Here, we show that unrestrained MD simulation of a solute in solution using a high-quality force field yields a conformational ensemble that is largely compatible with the experimental NMR data on the solute. Four 100 ns MD simulations of two forms of a nine-residue beta-peptide in methanol at two temperatures produced conformational ensembles that were used to interpret the NMR data on this molecule and resolve inconsistencies between the experimental NOEs. The protected and unprotected forms of the beta-peptide adopt predominantly a 12/10-helix in agreement with the qualitative interpretation of the NMR data. However, a particular NOE was not compatible with this helix indicating the presence of other conformations. The simulations showed that 3(14)()-helical structures were present in the ensemble of the unprotected form and that their presence correlates with the fulfillment of the particular NOE. Additionally, all inter-hydrogen distances were calculated to compare NOEs predicted by the simulations to the ones observed experimentally. The MD conformational ensembles allowed for a detailed and consistent interpretation of the experimental data and showed the small but specific conformational differences between the protected and unprotected forms of the peptide.     

Suppression of phase and amplitude J(HH) modulations in HSQC experiments

The amplitude and the phase of cross peaks in conventional 2D HSQC experiments are modulated by both proton–proton, J(HH), and proton–carbon, ¹J(CH), coupling constants. It is shown by spectral simulation and experimentally that J(HH) interferences are suppressed in a novel perfect‐HSQC pulse scheme that incorporates perfect‐echo INEPT periods. The improved 2D spectra afford pure in‐phase cross peaks with respect to ¹J(CH) and J(HH), irrespective of the experiment delay optimization. In addition, peak volumes are not attenuated by the influence of J(HH), rendering practical issues such as phase correction, multiplet analysis, and signal integration more appropriate. Copyright © 2014 John Wiley & Sons, Ltd.     

A conformational dynamics study of alpha-l-Rhap-(1-->2)[alpha-l-Rhap-(1-->3)]-alpha-l-Rhap-OMe in solution by NMR experiments and molecular simulations

The conformational preference of alpha-l-Rhap-(1-->2)[alpha-l-Rhap-(1-->3)]-alpha-l-Rhap-OMe in solution has been studied by NMR spectroscopy using one-dimensional (1)H,(1)H T-ROESY experiments and measurement of trans-glycosidic (3)J(C,H) coupling constants. Molecular dynamics (MD) simulations with a CHARMM22 type of force field modified for carbohydrates were performed with water as the explicit solvent. The homonuclear cross-relaxation rates, interpreted as effective proton-proton distances, were compared to those obtained from simulation. Via a Karplus torsional relationship, (3)J(C,H) values were calculated from simulation and compared to experimental data. Good agreement was observed between experimental data and the MD simulation, except for one inter-residue T-ROE between protons in the terminal sugar residues. The results show that the trisaccharide exhibits substantial conformational flexibility, in particular along the psi glycosidic torsion angles. Notably, for these torsions, a high degree of correlation (77%) was observed in the MD simulation revealing either psi(2)(+) psi(3)(+) or psi(2)(-)psi(3)(-) states. The simulations also showed that non-exoanomeric conformations were present at the phi torsion angles, but to a limited extent, with the phi(3) state populated to a larger extent than the phi(2) state. Further NMR analysis of the trisaccharide by translational diffusion measurements and (13)C T(1) relaxation experiments quantified global reorientation using an anisotropic model together with interpretation of the internal dynamics via the "model-free" approach. Fitting of the dynamically averaged states to experimental data showed that the psi(2)(+)psi(3)(+) state is present to approximately 49%, psi(2)(-) psi(3)(-) to approximately 39%, and phi(3) (non-exo) to approximately 12%. Finally, using a dynamic and population-averaged model, (1)H,(1)H T-ROE buildup curves were calculated using a full relaxation matrix approach and were found to be in excellent agreement with experimental data, in particular for the above inter-residue proton-proton interaction between the terminal residues.     

Atomistic Brownian dynamics simulation of peptide phosphorylation

We report the implementation of an all-atom Brownian dynamics simulation model of peptides using the constraint algorithm LINCS. The algorithm has been added as a part of UHBD. It uses adaptive time steps to achieve a balance between computational speed and stability. The algorithm was applied to study the effect of phosphorylation on the conformational preference of the peptide Gly-Ser-Ser-Ser. We find that the middle serine residue experiences considerable conformational change from the C(7eq) to the alpha(R) structure upon phosphorylation. NMR (3)J coupling constants were also computed from the Brownian trajectories using the Karplus equation. The calculated (3)J results agree reasonably well with experimental data for phosphorylated peptide but less so for doubly charged phosphorylated one.     

Solvent effects on one-bond B-Li coupling constants in boryllithium compounds

EOM-CCSD 11B-7Li coupling constants and B chemical shifts have been computed for Li-diazaborole and its complexes with one H2O or FLi molecule. B-Li coupling constants for a model compound H(2)BLi and its complexes with up to 4 H2O or FLi molecules have also been obtained in an attempt to resolve discrepancies between the computed values of these properties for isolated Li-diazaborole and experimentally determined values for boryllithium in a THF solution. The presence of solvent molecules increases the ion-pair character of the B-Li bond, with the result that 1J(B-Li) decreases systematically as the basicity and the number of solvent molecules increases. In the presence of even a single solvent molecule, the boron chemical shift for Li-diazaborole increases, and approaches the experimental value. The computed results emphasize the role of the solvent in determining these NMR properties.     

Structure and NMR spectra of some [2.2]paracyclophanes. The dilemma of [2.2]paracyclophane symmetry

Density functional theory (DFT) quantum chemical calculations of the structure and NMR parameters for highly strained hydrocarbon [2.2]paracyclophane 1 and its three derivatives are presented. The calculated NMR parameters are compared with the experimental ones. By least-squares fitting of the (1)H spectra, almost all J(HH) coupling constants could be obtained with high accuracy. Theoretical vicinal J(HH) couplings in the aliphatic bridges, calculated using different basis sets (6-311G(d,p), and Huz-IV) reproduce the experimental values with essentially the same root-mean-square (rms) error of about 1.3 Hz, regardless of the basis set used. These discrepancies could be in part due to a considerable impact of rovibrational effects on the observed J(HH) couplings, since the latter show a measurable dependence on temperature. Because of the lasting literature controversies concerning the symmetry of parent compound 1, D(2h) versus D(2), a critical analysis of the relevant literature data is carried out. The symmetry issue is prone to confusion because, according to some literature claims, the two hypothetical enantiomeric D(2) structures of 1 could be separated by a very low energy barrier that would explain the occurrence of rovibrational effects on the observed vicinal J(HH) couplings. However, the D(2h) symmetry of 1 with a flat energy minimum could also account for these effects.     

Experimental (NMR) and Theoretical (MD Simulations) Studies on the Conformational Preference of Three Cycloalkanols when Included in β-Cyclodextrin

The inclusion complexes between three cycloalkanols (cyclopentanol, cyclohexanol and cycloheptanol) and β-cyclodextrin (β-CD) have been studied by NMR experiments, and by molecular dynamics (MD) simulations. Complexes present medium to small association constants. All experimental data show the equatorial conformer as the most stable after complexation because no changes were detected in the coupling constants of the H1 protons. Intermolecular ROE experiments suggest that while cyclopentanol is deeply included into the β-CD cavity, cyclohexanol and cycloheptanol occupy mainly the wider entrance. The MD simulations agree with the experimental data (equatorial conformers are always the most stable), and average geometries coincide with those deduced from the ROE experiments.     

Simultaneous measurement of J(HH) and two different nJ(CH) coupling constants from a single multiply edited 2D cross‐peak

Three different J‐editing methods (IPAP, E.COSY and J‐resolved) are implemented in a single NMR experiment to provide spin‐state‐edited 2D cross‐peaks from which a simultaneous measurement of different homonuclear and heteronuclear coupling constants can be performed. A new J‐selHSQMBC‐IPAP experiment is proposed for the independent measurement of two different ⁿJ(CH) coupling constants along the F2 and F1 dimensions of the same 2D cross‐peak. In addition, the E.COSY pattern provides additional information about the magnitude and relative sign between J(HH) and ⁿJ(CH) coupling constants. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

Molecular Dynamics Simulations of Small Peptides: Can One Derive Conformational Preferences from ROESY Spectra?

Folding properties of beta-peptides were investigated by means of NMR experiments and MD simulations of beta-dipeptides, which serve as small test systems to study the influence of stereocenters and side chains on hydrogen-bond and consequently on secondary-structure formation. Two stereoisomers, SR and SS, of a Val-Phe dipeptide, and of the corresponding Ala-Ala dipeptide, and a Gly-Gly dipeptide were simulated in methanol for 40 ns. In agreement with experiment, the isomers of the Val-Phe dipeptide adopt quite different conformers at 298 K, the differences being reduced at 340 K. Interestingly, the SR isomer shows enhanced hydrogen bonding at the higher temperature. The adopted conformations are primarily determined by the R or S side chain substitution, and less by the type of side chain. Back-calculation of (1)H ROESY spectra and (3)J coupling constants from the MD simulations and comparison with the experimental data for the Val-Phe dipeptides shows good agreement between simulation and experiment, and reveals possible problems and pitfalls, when deriving structural properties of a small and extremely flexible molecule from NMR data only. Inclusion of all aspects of internal dynamics is essential to the correct prediction of the NMR spectra of these small molecules. Cross comparison of calculated with experimental spectra for both isomers shows that only a few out of many ROESY peaks reflect the sizeable conformational differences between the isomers at 298 K.     

RNA unrestrained molecular dynamics ensemble improves agreement with experimental NMR data compared to single static structure: a test case

Nuclear magnetic resonance (NMR) provides structural and dynamic information reflecting an average, often non-linear, of multiple solution-state conformations. Therefore, a single optimized structure derived from NMR refinement may be misleading if the NMR data actually result from averaging of distinct conformers. It is hypothesized that a conformational ensemble generated by a valid molecular dynamics (MD) simulation should be able to improve agreement with the NMR data set compared with the single optimized starting structure. Using a model system consisting of two sequence-related self-complementary ribonucleotide octamers for which NMR data was available, 0.3 ns particle mesh Ewald MD simulations were performed in the AMBER force field in the presence of explicit water and counterions. Agreement of the averaged properties of the molecular dynamics ensembles with NMR data such as homonuclear proton nuclear Overhauser effect (NOE)-based distance constraints, homonuclear proton and heteronuclear ¹H–³¹P coupling constant (J) data, and qualitative NMR information on hydrogen bond occupancy, was systematically assessed. Despite the short length of the simulation, the ensemble generated from it agreed with the NMR experimental constraints more completely than the single optimized NMR structure. This suggests that short unrestrained MD simulations may be of utility in interpreting NMR results. As expected, a 0.5 ns simulation utilizing a distance dependent dielectric did not improve agreement with the NMR data, consistent with its inferior exploration of conformational space as assessed by 2-D RMSD plots. Thus, ability to rapidly improve agreement with NMR constraints may be a sensitive diagnostic of the MD methods themselves.     

Convergence of replica exchange molecular dynamics

Replica exchange molecular dynamics (REMD) method is one of the generalized-ensemble algorithms which performs random walk in energy space and helps a system to escape from local energy traps. In this work, we studied the accuracy and efficiency of REMD by examining its ability to reproduce the results of multiple extended conventional molecular dynamics (MD) simulations and to enhance conformational sampling. Two sets of REMD simulations with different initial configurations, one from the fully extended and the other from fully helical conformations, were conducted on a fast-folding 21-amino-acid peptide with a continuum solvent model. Remarkably, the two REMD simulation sets started to converge even within 1.0 ns, despite their dramatically different starting conformations. In contrast, the conventional MD within the same time and with identical starting conformations did not show obvious signs of convergence. Excellent convergence between the REMD sets for T>300 K was observed after 14.0 ns REMD simulations as measured by the average helicity and free-energy profiles. We also conducted a set of 45 MD simulations at nine different temperatures with each trajectory simulated to 100.0 and 200.0 ns. An excellent agreement between the REMD and the extended MD simulation results was observed for T>300 K, showing that REMD can accurately reproduce long-time MD results with high efficiency. The autocorrelation times of the calculated helicity demonstrate that REMD can significantly enhance the sampling efficiency by 14.3+/-6.4, 35.1+/-0.2, and 71.5+/-20.4 times at, respectively, approximately 360, approximately 300, and approximately 275 K in comparison to the regular MD. Convergence was less satisfactory at low temperatures (T<300 K) and a slow oscillatory behavior suggests that longer simulation time was needed to reach equilibrium. Other technical issues, including choice of exchange frequency, were also examined.     

First-principles molecular dynamics evaluation of thermal effects on the NMR (1)J(Li,C) spin-spin coupling

Car-Parrinello (CP) molecular dynamics were applied to sample conformations of various models of organolithium aggregates which are chosen to estimate (1)J(Li,C) NMR coupling constants. The results show that the deviations from the values computed using static (optimized) geometries are small provided no large-amplitude motions occur within the timescale of the simulations. In the case of the vinyllithium dimer, for which rotation of the vinyl chain is observed, this approach allows analysis of the various contributions to the experimentally measured constants. For the trisolvated methyllithium monomer, partial decoordination of solvating dimethyl ether is observed and results in a significant shift of (1)J(Li,C). All these results highlight that a varied physicochemical machinery is hidden behind general empirical formulas, such as the Bauer-Winchester-Schleyer rule used experimentally.     

Theoretical studies of nuclear magnetic resonance parameters for the proton-exchange pathways in porphyrin and porphycene

The nuclear magnetic resonance (NMR) parameters in porphyrin and porphycene have been calculated to investigate their changes during the process of proton exchange, using density-functional theory (DFT) for both the spin-spin coupling constants and the shielding constants. In addition, in calculations on the smaller 1,3-bis(arylimino)isoindoline molecule, we have tested the performance of our computational approach against experimental data. The calculated nuclear spin-spin coupling constants and shielding constants have been analyzed as functions of the progress of the proton transfer between two nitrogen atoms. The one-bond couplings between proton and nitrogen, dominated by the Fermi-contact term, decay steeply as the internuclear distance increases. The small changes in the intramolecular J(HH) coupling between two inner protons are mainly determined by the sum of relatively large spin-orbit terms. The isotropic shielding constant shows a strong deshielding of the nitrogen nuclei as the proton migrates away. Both the isotropic shielding of the exchanged protons and the shielding anisotropy exhibit a minimum close to the transition states.     

High-resolution NMR spectra under inhomogeneous fields via intermolecular double-quantum coherences

High-resolution NMR spectroscopy is a powerful tool for analyzing molecular structures and compositions. Line-widths of conventional liquid NMR signals are directly proportional to the overall magnetic field inhomogeneity the sample experiences. In many circumstances, spatial and temporal homogeneity of the magnetic field is degraded. In this paper, a modified pulse sequence based on intermolecular double-quantum coherences (iDQCs) was proposed to obtain 1D high-resolution NMR spectra under inhomogeneous fields using 2D acquisition. Analytical expressions were derived from the intermolecular multiple-quantum coherence (iMQC) treatments. Both experimental and simulated spectra provide high-resolution 1D projection spectra similar to conventional 1D high-resolution spectra. Moreover, the apparent J coupling constants are threefold magnified, which allows a more accurate measurement of small J coupling constants.     

A theoretical study of the large Hg-Hg spin-spin coupling constants in Hg(2)(2+), Hg(3)(2+), and Hg(2)(2+)-crown ether complexes

Nuclear spin-spin coupling constants (1)J(Hg-Hg) in the systems Hg(2)(2+) and Hg(3)(2+) represent the largest coupling constants so far observed in NMR experiments. We have performed a computational study on these ions, on Hg(2)(2+) complexes with 18-crown-6 and 15-crown-5, and on Hg(3)(2+) with solvent molecules and counterions. The results obtained with our recently developed program for the density functional computation of heavy nucleus spin-spin coupling constants are in good agreement with experiments. The data reveal that the bare ions Hg(2)(2+) and Hg(3)(2+) would afford much larger coupling constants than those experimentally observed, with an upper limit of approximately 0.9 MHz for Hg(2)(2+). This limit is much larger than that previously estimated by Hückel theory. It is demonstrated that in solution or due to complexation the experimentally determined values are much smaller than the free ion's coupling constants. With the help of intuitive MO arguments, it is illustrated how the environment strongly reduces the coupling constants in Hg(2)(2+) and Hg(3)(2+). The two-bond coupling constant (2)J(Hg-Hg) in Hg(3)(2+) is also examined.     

Structures in solutions from joint experimental-computational analysis: applications to cyclic molecules and studies of noncovalent interactions

The potential of an approach combining nuclear magnetic resonance (NMR) spectroscopy, molecular dynamics (MD) simulations, and quantum mechanical (QM) calculations for full structural characterizations in solution is assessed using cyclic organic compounds, namely, benzazocinone derivatives 1-3 with fused five- and eight-membered aliphatic rings, camphoric anhydride 4, and bullvalene 5. Various MD simulations were considered, using force field and semiempirical QM treatments, implicit and explicit solvation, and high-temperature MD calculations for selecting plausible molecular geometries for subsequent QM geometry optimizations using mainly B3LYP, M062X, and MP2 methods. The QM-predicted values of NMR parameters were compared to their experimental values for verification of the final structures derived from the MD/QM analysis. From these comparisons, initial estimates of quality thresholds (calculated as rms deviations) were 0.7-0.9 Hz for (3)J(HH) couplings, 0.07-0.11 ? for interproton distances, 0.05-0.08 ppm for (1)H chemical shifts, and 1.0-2.1 ppm for (13)C chemical shifts. The obtained results suggest that the accuracy of the MD analysis in predicting geometries and relative conformational energies is not critical and that the final geometry refinements of the structures selected from the MD simulations using QM methods are sufficient for correcting for the expected inaccuracy of the MD analysis. A unique example of C(sp(3))-H···N(sp(3)) intramolecular noncovalent interaction is also identified using the NMR/MD/QM and the natural bond orbital analyses. As the NMR/MD/QM approach relies on the final QM geometry optimization, comparisons of geometric characteristics predicted by different QM methods and those from X-ray and neutron diffraction measurements were undertaken using rigid and flexible cyclic systems. The joint analysis shows that intermolecular noncovalent interactions present in the solid state alter molecular geometries significantly compared to the geometries of isolated molecules from QM calculations.     

Probing the proton-transfer coordinate of complexes with F--H...P hydrogen bonds using one- and two-bond spin-spin coupling constants

Scalar coupling constants have been computed using the EOM-CCSD method for equilibrium structures of complexes stabilized by F--H...P hydrogen bonds, as well as structures along the proton-transfer coordinates of these complexes. Variations in the signs and absolute values of (1)J(F--H), (1h)J(H--P) and (2h)J(F--P) have been analyzed and interpreted in terms of changing hydrogen bond type. Of the three phosphorus bases (phosphine, trimethylphosphine and phosphinine) investigated in this study, trimethylphosphine forms the strongest complex with FH, and has the largest two-bond F--P coupling constant. Among the relatively simple phosphorus bases, it would appear to be a leading candidate for experimental NMR study. Similarities and differences are noted between the corresponding coupling constants (J) and the reduced coupling constants (K) across F--H...P and F--H...N hydrogen bonds.