A conformational study of alpha-L-Rhap-(1----2)-alpha-L-Rhap-(1----OMe) by NMR nuclear Overhauser effect spectroscopy (NOESY) and molecular dynamics calculations.

The conformational preference of the disaccharide alpha-L-Rhap-(1----2)-alpha-L-Rhap-(1----OMe) (1) about the glycosidic torsion angles, phi and psi, was studied by NMR NOESY spectroscopy and molecular mechanics calculations. The NOE data were consistent with either of two distinct conformations close to minima on a calculated phi/psi potential energy surface. Starting from the lowest energy conformation, a 1-ns molecular dynamics (MD) trajectory was computed in vacuo, from which the NOE curves were simulated and compared to the experimentally observed NOESY data.     

Conformational studies of human milk oligosaccharides using (1)H-(13)C one-bond NMR residual dipolar couplings

1H-(13)C one-bond dipolar coupling values were measured for natural abundance samples of the human milk oligosaccharides "lacto-N-fucopentaose" (LNF-1 LNF-2, and LNF-3), "lacto-N-difucohexaose" (LND-1), "lacto-N-tetraose" (LNT), and "lacto-N-neo-tetraose" (LNnT), four of which have Lewis blood group epitopes. Each oligosaccharide was dissolved in a 7.5% solution of 1, 2-dimyristoyl-sn-glycero-3-phosphocholine/1, 2-dihexanoyl-sn-glycero-3-phosphocholine (DMPC/DHPC) bicelle liquid crystals oriented in the NMR magnetic field. The dipolar coupling data and NOE were fitted to conformational models with calculations of an optimum orientation tensor which best represents the dipolar coupling values for a fragment hypothesized to adopt a single conformation. In the case of LNF-1, LNF-2, LNF-3, and LND-1, the models confirm previous conformational models for the Lewis epitopes based on NOE and molecular dynamics simulations. Extensions of the model provided new structural data for the remaining residues. In all cases, upper limits for the errors in the glycosidic angles of the models were estimated. Since residual dipolar coupling provides information on long-range order, it is a valuable complement to other types of NMR data such as NOE and scalar coupling for exploring conformations of complex oligosaccharides.     

Sucrose in aqueous solution revisited, Part 2: adaptively biased molecular dynamics simulations and computational analysis of NMR relaxation

We report 100 ns molecular dynamics simulations, at various temperatures, of sucrose in water (with concentrations of sucrose ranging from 0.02 to 4M), and in a 7:3 water‐DMSO mixture. Convergence of the resulting conformational ensembles was checked using adaptive‐biased simulations along the glycosidic Φ and ψ torsion angles. NMR relaxation parameters, including longitudinal (R₁) and transverse (R₂) relaxation rates, nuclear Overhauser enhancements (NOE), and generalized order parameter (S²) were computed from the resulting time‐correlation functions. The amplitude and time scales of molecular motions change with temperature and concentration in ways that track closely with experimental results, and are consistent with a model in which sucrose conformational fluctuations are limited (with 80–90% of the conformations having ??ψ values within 20° of an average conformation), but with some important differences in conformation between pure water and DMSO‐water mixtures. © 2011 Wiley Periodicals, Inc. Biopolymers 97: 289–302, 2012.     

Conformational analysis and molecular dynamics simulations of maltose

Constrained conformational energy minimizations have been used to calculate an adiabatic (Φ, ψ) potential energy surface for the disaccharide β-maltose. The inclusion of molecular flexibility in the conformational energy analysis of the disaccharide was found to significantly lower the barriers to conformational transitions, as has been observed previously for other systems. Several low energy wells were identified on the adiabatic surface which differ in energy by small amounts and with low absolute barriers separating them, indicating the possibility of a non-negligible equilibrium population distribution in each well. If such a distribution of conformations existed in the physical system, the conformation observed by NMR NOE measurements would thus be a “virtual” conformation. Molecular dynamics simulations of the motions of this molecule in vacuum were also conducted and indicate that the rate of relaxation of the molecule to the adiabatic surface may be slower than the typical timescale of conformational fluctuations. This effect is apparently due to an unphysical persistence of hydrogen bond patterns in vacuum which does not occur in aqueous solution. Trajectories undergoing transitions between wells were calculated and the effects of such conformational transitions upon the ensemble mean structure, such as might be observed in an NMR experiment, were demonstrated.     

Exploring the dynamic information content of a protein NMR structure: Comparison of a molecular dynamics simulation with the NMR and X‐ray structures of Escherichia coli ribonuclease HI

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7‐ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high‐resolution (1.48 Å) crystal structure than with either the lower‐resolution (2.05 Å) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side‐chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn‐pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean‐square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non‐rotameric/syn‐pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X‐ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed. Proteins 1999;36:87–110. © 1999 Wiley‐Liss, Inc.     

The determination of the conformational properties of nucleic acids in solution from NMR data

A program, NUCFIT, has been written for simulating the effects of conformational averaging on nuclear Overhauser enhancement (NOE) intensities for the spin systems found in nucleic acids. Arbitrary structures can be generated, and the NOE time courses can be calculated for truncated one-dimensional NOEs, two-dimensional NOE and rotating frame NOE spectroscopy (NOESY and ROESY) experiments. Both isotropic and anisotropic molecular rotation can be treated, using Woessner's formalism (J. Chem. Phys. (1962) 37, 647-654). The effects of slow conformational averaging are simulated by taking population-weighted means of the conformations present. Rapid motions are allowed for by using order parameters which can be supplied by the user, or calculated for specific motional models using the formalism of Tropp (J. Chem. Phys. (1980) 72, 6035-6043). NOE time courses have been simulated for a wide variety of conformations and used to determine the quality of structure determinations using NMR data for nucleic acids. The program also allows grid-searching with least-squares fitting of structures to experimental data, including the effects of spin-diffusion, conformational averaging and rapid internal motions. The effects of variation of intra and internucleotide conformational parameters on NOE intensities has been systematically explored. It is found that (i) the conformation of nucleotides is well determined by realistic NOE data sets, (ii) some of the helical parameters, particularly the base pair roll, are poorly determined even for extensive, noise-free data sets, (iii) conformational averaging of the sugars by pseudorotation has at most second-order influence on the determination of other parameters and (iv) averaging about the glycosidic torsion bond also has, in most cases, an insignificant effect on the determination of the conformation of nucleotides.     

Modeling conformational ensembles of slow functional motions in Pin1-WW

Protein-protein interactions are often mediated by flexible loops that experience conformational dynamics on the microsecond to millisecond time scales. NMR relaxation studies can map these dynamics. However, defining the network of inter-converting conformers that underlie the relaxation data remains generally challenging. Here, we combine NMR relaxation experiments with simulation to visualize networks of inter-converting conformers. We demonstrate our approach with the apo Pin1-WW domain, for which NMR has revealed conformational dynamics of a flexible loop in the millisecond range. We sample and cluster the free energy landscape using Markov State Models (MSM) with major and minor exchange states with high correlation with the NMR relaxation data and low NOE violations. These MSM are hierarchical ensembles of slowly interconverting, metastable macrostates and rapidly interconverting microstates. We found a low population state that consists primarily of holo-like conformations and is a "hub" visited by most pathways between macrostates. These results suggest that conformational equilibria between holo-like and alternative conformers pre-exist in the intrinsic dynamics of apo Pin1-WW. Analysis using MutInf, a mutual information method for quantifying correlated motions, reveals that WW dynamics not only play a role in substrate recognition, but also may help couple the substrate binding site on the WW domain to the one on the catalytic domain. Our work represents an important step towards building networks of inter-converting conformational states and is generally applicable.     

The conformational behaviour of the cardiac glycoside digoxin as indicated by NMR spectroscopy and molecular dynamics calculations.

The 1H- and 13C-NMR spectra of digoxin in solution in Me2SO-d6 have been assigned completely. Measurement of the 3JC,H values has enabled estimation of the torsional angles involving the bonds linking the digitoxose residues, between the inner digitoxose and the genin unit, and for the unsaturated gamma-lactone ring. These values have been supplemented by 1H-1H NOE data. In general, there is good agreement between the conformations in solution (NMR data) and the solid state (X-ray data), and that derived from theoretical modelling which shows evidence of conformational flexibility. The major difference occurs for the torsion between the genin and the innermost digitoxose residue where molecular dynamics predict the presence of two conformations, one similar to that seen by NMR and the other similar to the X-ray structure.     

Multiple helical conformations of the helix‐turn‐helix region revealed by NOE‐restrained MD simulations of tryptophan aporepressor,TrpR

The nature of flexibility in the helix‐turn‐helix region of E. coli trp aporepressor has been unexplained for many years. The original ensemble of nuclear magnetic resonance (NMR structures showed apparent disorder, but chemical shift and relaxation measurements indicated a helical region. Nuclear Overhauser effect (NOE) data for a temperature‐sensitive mutant showed more helical character in its helix‐turn‐helix region, but nevertheless also led to an apparently disordered ensemble. However, conventional NMR structure determination methods require all structures in the ensemble to be consistent with every NOE simultaneously. This work uses an alternative approach in which some structures of the ensemble are allowed to violate some NOEs to permit modeling of multiple conformational states that are in dynamic equilibrium. Newly measured NOE data for wild‐type aporepressor are used as time‐averaged distance restraints in molecular dynamics simulations to generate an ensemble of helical conformations that is more consistent with the observed NMR data than the apparent disorder in the previously reported NMR structures. The results indicate the presence of alternating helical conformations that provide a better explanation for the flexibility of the helix‐turn‐helix region of trp aporepressor. Structures representing these conformations have been deposited with PDB ID: 5TM0. Proteins 2017; 85:731–740. © 2016 Wiley Periodicals, Inc.     

Conformational properties of a cyclic oligosaccharide: cyclotrikis-(1→6)-[α-D-glucopyranosyl-(1→4)-β-D-glucopyranosyl]

The title compound is a cyclic oligosaccharide having six glucopyranose residues linked alternatively by -(14) and -(16) glycosidic linkages. Like cyclodextrin analogues it is expected to exhibit an internal cavity and to form inclusion complexes with other species. In order to investigate its conformational preferences, an extensive conformational search was carried out using a combination of Metropolis Monte-Carlo (MMC) procedure in the glycosidic torsion angle space and molecular mechanics procedures. To this end a specific program (METROCYCLIX) was developed. To reduce the MMC search, conformational maps of parent disaccharides were considered as starting entries. Fully minimized conformations were gathered into families using a clustering technique based on RMS fitting over the glycosidic torsion angle values. A wide range of local energy minima were identified in spite of ring closure conditions that constrained the structure of the oligosaccharide. Low energy conformers were stabilized by intramolecular interactions between distant residues. From the Bolzmann population of the best structures derived from the clustering results, various average properties were calculated and compared with experimental data obtained by high resolution NMR. Interpretation of these experimental values (heteronuclear coupling constants, rotating frame nuclear Overhauser effects, relaxation times) relies on the use of Karplus like equations (coupling constants) and analysis of the full relaxation rate matrix treatment (ROE). The quality of the molecular modelling strategy used is assessed by the agreement obtained between calculated and measured observables.     

New strategy for the conformational analysis of carbohydrates based on NOE and 13C NMR coupling constants. Application to the flexible polysaccharide of Streptococcus mitis J22.

For complex oligosaccharides, which are relatively rigid with modest excursions from a single minimum energy conformation, it is straightforward to build conformational models from NOE data. Other oligosaccharides are more flexible with transitions between distinct minima separated by substantial energy barriers. We show that modeling based on scalar coupling data is superior to NOE-based modeling for the latter case. Long range 13C-13C and 13C-1H coupling constants measured for the heptasaccharide repeating subunit of the cell wall polysaccharide from Streptococcus mitis J22 are correlated with individual glycosidic dihedral angles, effectively uncoupling the degrees of freedom of the oligosaccharide and allowing a search for combinations of dihedral angles which are energetically reasonable, i.e., with no bad van der Waals contacts, and which can be combined to satisfy all the measured J values. Allowed values of the individual angles can then be combined to search for overall oligosaccharide conformations which contribute to the ensemble. We show that while the polysaccharide from S. mitis J22 is flexible, requiring multiple conformations, most of the flexibility is localized to a few bonds and only a rather small number of conformations is required to reproduce the experimental NOE and scalar coupling data.     

The confromational behaviour of the cardiac glycoside digoxin as indicated by NMR spectroscopy and molecular dynamics calculations

The ¹H- and ¹³C-NMR spectra of dogoxin in solution in Me₂So-d₆ have been assigned completely. Measurement of the ³J_C,H values has enabled estimation of the torsional angles involving the bonds linking the digitoxose residues, between the inner digitoxose and the genin unit, and for the unsaturated γ-lactone ring. These values have been supplemented by ¹H---¹H NOE data. In general, there is good agreement between the conformations in solution (NMR data) and the solid state (X-ray data), and that derived from theological modelling which shows evidence of conformational flexibility. The major difference occurs for the torsion between the genin and the innermost digitoxose residue where molecular dynamics predict the presence of two conformations, one similar to that seen by NMR and the other similar to the X-ray structure.     

Interaction of berenil with the EcoRI dodecamer d(CGCGAATTCGCG)2 in solution studied by NMR

The conformation of the EcoRI dodecamer d(CGCGAATTCGCG)2 has been examined in solution by 1H and 31P NMR. Spin-spin coupling constants and nuclear Overhauser (NOE) enhancement spectroscopy show that all deoxyriboses lie in the south domain, with a small admixture of the north conformation (0-20%). The time dependence of the nuclear Overhauser enhancements also reveals a relatively uniform conformation at the glycosidic bonds (average angle, chi = -114 degrees). The average helical twist is 36.5 degrees (9.8 base pairs per turn). Tilt angles are small (in the range 0 to -10 degrees), and roll angles are poorly determined. Unlike single-crystal X-ray studies of the same sequence, there is no evidence for asymmetry in the structure. Both the NOE intensities and 31P relaxation data imply conformational anomalies at the C3-G4/C9-G10 and the A5-A6/T7-T8 steps. Berenil binds in 1:1 stoichiometry to the dodecamer with high affinity (Kd = 1 microM at 298 K) and causes substantial changes in chemical shifts of the sugar protons of nucleotides Ado 5-Cyt 9 and of the H2 resonances of the two Ado residues. No significant asymmetry appears to be induced in the DNA conformation on binding, and there is no evidence for intercalation, although the binding site is not centrosymmetric. NOEs are observed between the aromatic protons of berenil and the H1' of both Thy 7 and Thy 8, as well as to Ado 5 and Ado 6 H2. These results firmly establish that berenil binds via the minor groove and closely approaches the nucleotides Ado 6, Thy 7, and Thy 8. On the basis of quantitative NOE spectroscopy and measurements of spin-spin coupling constants, changes in the conformations of the nucleotides are found to be small. Using the observed NOEs between the ligand and the DNA together with the derived glycosidic torsion angles, we have built models that satisfy all of the available solution data. The berenil molecule binds at the 5'-AAT (identical to 5'-ATT on the complementary strand) site such that (i) favorable hydrogen bonds are formed between the charged amidinium groups and the N3 atoms of Ado 6 and Ado 18 and (ii) the ligand is closely isohelical with the floor of the minor groove.     

NMR and molecular modeling studies of the interaction between wheat germ agglutinin and the beta-D-GlcpNAc-(1-->6)-alpha-D-Manp epitope present in glycoproteins of tumor cells

The beta-D-GlcpNAc-(1-->6)-alpha-D-Manp disaccharide is a constituent of highly branched cell-surface glycoconjugates that are malignancy markers. The conformational preference of the disaccharide beta-D-GlcpNAc-(1-->6)-alpha-D-Manp-OMe in solution has been studied by molecular modeling and NMR spectroscopy including 1D (1)H,(1)H T-ROESY experiments and analysis of (3)J(H,H) of the hydroxymethyl group being part of the glycosidic linkage of the disaccharide, which revealed the relative populations of the omega torsion angle as gt = 0.60, gg = 0.35, and tg = 0.05. Good agreement was obtained between the effective proton-proton distances from the experiment and those obtained by molecular modeling when the flexibility at the omega torsion angle was taken into account. Molecular modeling of the disaccharide in the binding sites of the lectin wheat germ agglutinin indicates that several conformations could be adopted in the bound state. (1)H NMR and transfer NOESY experiments confirmed that binding took place, and trans-glycosidic proton-proton interactions indicated that a conformational preference was present in the bound state, as observed by the relative change of the NOEs from H1' to H6(pro-R) and H6(pro-S). STD NMR experiments showed that binding occurred in the region of the N-acetyl group of the terminal sugar residue. In addition, the O-methyl group received saturation transfer because of the proximity to the protein. (1)H,(1)H NOEs indicated that the two methyl groups were close in space, as observed in only one of the predicted bound conformations. Experimental and theoretical data therefore agree that one conformation with a gt conformation of the hydroxymethyl group and a negative sign for the psi torsion angle is indeed selected by the lectin upon binding.     

Conformational dynamics of sialyl Lewisx in aqueous solution and its interaction with selectinE. A study by molecular dynamics

Three dimensional structures of sialyl Lewis(x) (SLe(x)) in aqueous solution and bound to selectinE are described based on an exhaustive conformational analysis and several long molecular dynamics simulations using different glycosidic regions as starting conformations. It appears from this study that when the oligosaccharide is free in solution the NeuNAcalpha(2-3)Gal segment favors glycosidic conformation in three different regions in the (Phi,Psi) plane with propensity of populations in the ratio 1:8:1. Each one of these structures is characteristically stabilized by specific hydrogen bonding interaction between NeuNAc and Gal. On the other hand, the Gal-GlcNAc-Fuc segment can exist in four different conformational states. Based on the topology of SLe(x) we are able to predict that out of all the allowed conformations in solution only two of these structures possess a geometry that would fit without steric clashes into the binding location of selectinE. In both of these binding modes, segment Gal-GlcNAc-Fuc adopts a unique conformation. The only difference between the two SLe(x) conformers that can successfully bind to selectinE is given by two possible regions in glycosidic space in the fragment NeuNAcalpha(2-3)Gal. A large conformational departure from the crystallographic data is observed for two lysine residues at the binding site of selectinE. These two residues play an important role when SLe(x) binds selectinE in aqueous solution. These findings help reconcile the X-ray data, in which these residues appear to be 1 nm away from SLe(x), with recent liquid NMR data reporting couplings between these protein residues and the sugar.     

The use of NMR residual dipolar couplings in aqueous dilute liquid crystalline medium for conformational studies of complex oligosaccharides

C-H dipolar coupling values were measured for a natural-abundance sample of the pentasaccharide beta-D-Galp-(1-->3)-[alpha-L-Fucp-(1-->4)]-beta-D-GlcNAcp-(1 -->3)-beta-D- Galp-(1-->4)-beta-D-Glcp ('lacto-N-fucopentaose 2') (LNF-2), in a 7.5% solution of dimyristoyl phosphatidylcholine-dihexanoyl phosphatidylcholine bicelle liquid crystals oriented in the NMR magnetic field. Interpretation of the dipolar coupling data and NOE confirms the conformational model for the Lewis(a) trisaccharide epitope based on NOE, molecular dynamics simulations, and scalar coupling data and provided new structural information for the remaining residues of the pentasaccharide. Since residual dipolar coupling provides information on long-range order, it is a valuable complement to other types of NMR data such as NOE and scalar coupling for exploring conformations of complex oligosaccharides.     

Reevaluation of rotamer populations for 1,6 linkages: reconciliation with potential energy calculations

Applications of ensemble averaging to the solution conformation of model compounds for N-linked glycans are further investigated. Specifically, the interpretation usually applied to observed values of J5,6', a parameter reflecting the rotameric distribution about C5-C6 bonds (torsion angle omega in 1,6 glycosidic linkages) in 6-O-substituted hexopyranosides, was found to be inconsistent with populations derived from potential energy calculations. However, agreement between observed and calculated, ensemble-averaged values of J5,6' was obtained and the distribution of omega rotamers reinterpreted. Values of J5,6' that were previously interpreted as indicative of equipartition between two rotamers in fact reflect a marked preference for one of them. Additional potential energy terms, previously absent from energy calculations, are introduced and shown to be without effect on interpretations of omega rotamer distributions. From comparisons with both NMR relaxation and scalar coupling constant data, it is concluded that a simple empirical algorithm, HSEA, calculating van der Waals, exo-anomeric, and (as appropriate) hydrogen-bonding terms, is best suited for describing the population distributions in solution for oligosaccharides and N-linked glycans.     

Structural study of Ac-Phe-[Orn-Pro-dCha-Trp-Arg], a potent C5a receptor antagonist, by NMR

The cyclic hexapeptide Ac(0)-Phe(1)-[Orn(2)-Pro(3)-dCha(4)-Trp(5)-Arg(6)] (the square brackets denote cyclization) is a potent antagonist against C5a (the a-fragment of complement protein C5) binding to C5a receptor (C5aR) and an excellent candidate to become a therapeutic agent against diseases that involve unregulated activation of the complement system. We present the solution structure determination of this cyclic C5aR peptide antagonist (cC5aR-pa), using nuclear magnetic resonance (NMR) data and restrained molecular dynamics-based simulated annealing in torsion angle space with NMR-derived distance and torsion angle restraints. The calculated NMR ensemble of structures demonstrates the presence of a predominant conformation of a distorted type II' beta-turn in the segment Pro(3)-dCha(4)-Trp(5)-Arg(6). We critically examine the calculated structure with measured NMR parameters, such as nuclear Overhauser enhancement (NOE) connectivity patterns and intensities characteristic of specific structures, (3)J(H(N)-H(alpha)) scalar coupling constants, temperature coefficients for NH groups, and differences between observed chemical shifts and their random coil values. The raw NMR data are consistent with the presence of the type II' beta-turn, but also indicate the presence of conformational inter-conversion. The calculated three-dimensional coordinates for cC5aR-pa will form the basis for further computational studies and for the development of pharmacophore models.