Structure determination of [d(ATATATAUAT)]2 via two-dimensional NOE spectroscopy and molecular dynamics calculations

Proton homonuclear two-dimensional (2D) NOE spectra were obtained for the decamer [d(ATATATAUAT)]2 as a function of mixing time, and proton resonance assignments were made. Quantitative assessment of the 2D NOE cross-peak intensities was used in conjunction with the program MARDIGRAS, which entails a complete relaxation matrix analysis of the 2D NOE peak intensities, to obtain a set of upper and lower bound interproton distance constraints. The analysis with MARDIGRAS was carried out using three initial models: A-DNA, B-DNA and Z-DNA. The distance constraints determined were essentially the same regardless of initial structure. These experimental structural constraints were used with restrained molecular dynamics calculations to determine the solution structure of the decamer. The molecular dynamics program AMBER was run using A-DNA or B-DNA as starting model. The root-mean-square (rms) difference between these two starting models is 0.504 nm. The two starting models were subjected to 22.5 ps of restrained molecular dynamics calculations. The coordinates of the last 10.5 ps of the molecular dynamics runs were averaged to give two final structures. MDA and MDB. The rms difference between these two structures is 0.09 nm, implying convergence of the two molecular dynamics runs. The 2D NOE spectral intensities calculated for the derived structures are in good agreement with experimental spectra, based on sixth-root residual index analysis of intensities. A detailed examination of the structural features suggests that while the decamer is in the B-family of DNA structures, many torsion angle and helical parameters alternate from purine to pyrimidine, with kinks occurring at the U-A steps.     

1H NMR assignments of sidechain conformations in proteins using a high-dimensional potential in the simulated annealing calculations

A high-dimensional potential representing distance constraints for stereospecifically assignable diastereotopic proton or methyl pairs was incorporated into the dynamical simulated annealing protocol to calculate structure with stereospecifically determined sidechain conformations. The protocol is tested on nuclear magnetic resonance cross-relaxation data of a trypsin inhibitor from squash seeds, CMTI-I, and compared with two other methods of stereospecific assignment, the floating chirality and coupling constant methods. There is good agreement between the three methods in predicting the same stereospecific assignments. Because the high-dimensional potential uses more relaxed absolute distance constraints and also takes into account the relative distance constraint patterns, it avoids possible overinterpretation of the NOE data.     

Statistical analysis of double NOE transfer pathways in proteins as measured in 3D NOE-NOE spectroscopy

Summary The recent development of three-dimensional NMR spectroscopy has alleviated the problem of overlap of resonances. However, also for the 3D experiments resonance assignment strategies have usually relied upon knowledge about spin systems, combined with information about short (sequential) distances. For doubly (¹⁵N/¹³C)-labelled molecules, a novel assignment strategy has been developed. In this paper we address the possibilities of an assignment strategy for proteins, based solely upon the use of NOE data. For this, the 3D NOE-NOE experiment seems most suitable. Therefore, we have made a theoretical evaluation of double NOE transfer pathways in 28 protein crystal structures. We identify 95 connectivities which are most likely to be observed as cross peaks in a 3D NOE-NOE spectrum of a protein. Given the occurrence of one of these 95 connectivities, we evaluate the chances of occurrence for the others. Analysis of these conditional probabilities allowed the construction of five patterns of related, highly correlated cross peaks which resemble the conventional idea of spin systems to some extent and may provide a basis for assignment and secondary structure analysis from 3D NOE-NOE data alone.Dedicated to the memory of Professor V.F. Bystrov     

Solution structure of [d(GTATATAC)]2 via restrained molecular dynamics simulations with nuclear magnetic resonance constraints derived from relaxation matrix analysis of two-dimensional nuclear Overhauser effect experiments

Two-dimensional nuclear Overhauser effect (2D NOE) spectra have been used as the experimental basis for determining the solution structure of the duplex [d(GTATATAC)]2 employing restrained molecular dynamics (rMD) simulations. The MARDIGRAS algorithm has been employed to construct a set of 233 interproton distance constraints via iterative complete relaxation matrix analysis utilizing the peak intensities from the 2D NOE spectra obtained for different mixing times and model structures. The upper and lower bounds for each of the constraints, defining size of a flat-well potential function term used in the rMD simulations, were conservatively chosen as the largest or smallest value calculated by MARDIGRAS. Three different starting models were utilized in several rMD calculations: energy-minimized A-DNA, B-DNA, and a structure containing wrinkled D-DNA in the interior. Considerable effort was made to define the appropriate force constants to be employed with the NOE terms in the AMBER force field, using as criteria the average constraints deviation, the constraints violation energy and the total energy. Of the 233 constraints, one was generated indirectly, but proved to be crucial in defining the structure: the cross-strand A5-H2 A5-H2 distance. As those two protons resonate isochronously for the self-complementary duplex, the distance cannot be determined directly. However, the general pattern of 2D NOE peak intensities, spin-lattice relaxation time (T1) values, and 31P nuclear magnetic resonance spectra lead to use of the A3-H2 A7-H2 distance for A5-H2 A5-H2 as well. Five rMD runs, with different random number seeds, were made for each of the three starting structures with the full distance constraint set. The average structure from all 15 runs and the five-structure averages from each starting structure were all quite similar. Two rMD runs for each starting structure were made with the A5-H2 A5-H2 constraint missing. The average of these six rMD runs revealed differences in structure, compared to that with the full set of constraints, primarily for the middle two base-pairs involving the missing cross-strand constraint but global deviations also were found. Conformational analysis of the resulting structures revealed that the inner four to six base-pairs differed in structure from the termini. Furthermore, an alternating structure was suggested with features alternating for the A-T and T-A steps.     

Relaxation matrix refinement of the solution structure of squash trypsin inhibitor

The structure of the small squash trypsin inhibitor CMTI-I is refined by directly minimizing the difference between the observed two-dimensional nuclear Overhauser enhancement (NOE) intensities and those calculated by the full relaxation matrix approach. To achieve this, a term proportional to this difference was added to the potential energy function of the molecular dynamics program X-PLOR. Derivatives with respect to atomic co-ordinates are calculated analytically. Spin diffusion effects are thus accounted for fully during the refinement. Initial structures for the refinement were those determined recently by solution nuclear magnetic resonance using the isolated two-spin approximation to derive distance range estimates. The fits to the nuclear magnetic resonance data improve significantly with only small shifts in the refined structures during a few cycles of conjugate gradient minimization. However, larger changes (approximately 1 A) in the conformation occur during simulated annealing, which is accompanied by a further reduction of the difference between experimental and calculated two-dimensional NOE intensities. The refined structures are closer to the X-ray structure of the inhibitor complexed with trypsin than the initial structures. The root-mean-square difference for backbone atoms between the initial structures and the X-ray structure is 0.96 A, and that between the refined structures and the X-ray structure 0.61 A.     

Solution structure of a DNA octamer containing the Pribnow box via restrained molecular dynamics simulation with distance and torsion angle constraints derived from two-dimensional nuclear magnetic resonance spectral fitting.

The DNA octamer [d(GTATAATG].[(CATATTAC)], containing the prokaryotic upstream consensus recognition sequence, has been examined via proton homonuclear two-dimensional nuclear Overhauser effect (2D NOE) and double-quantum-filtered correlation (2QF-COSY) spectra. All proton resonances, except those of H5' and H5" protons, were assigned. A temperature dependence study of one-dimensional nuclear magnetic resonance (NMR) spectra, rotating frame 2D NOE spectroscopy (ROESY), and T1 rho measurements revealed an exchange process that apparently is global in scope. Work at lower temperatures enabled a determination of structural constraints that could be employed in determination of a time-averaged structure. Simulations of the 2QF-COSY cross-peaks were compared with experimental data, establishing scalar coupling constant ranges of the individual sugar ring protons and hence pucker parameters for individual deoxyribose rings. The rings exhibit a dynamic equilibrium of N and S-type conformers with 80 to 100% populations of the latter. A program for iterative complete relaxation matrix analysis of 2D NOE spectral intensities, MARDIGRAS, was employed to give interproton distances for each mixing time. According to the accuracy of the distance determination, upper and lower distance bounds were chosen. The distance bounds define the size of a flat-well potential function term, incorporated into the AMBER force-field, which was employed for restrained molecular dynamics calculations. Torsion angle constraints in the form of a flat-well potential were also constructed from the analysis of the sugar pucker data. Several restrained molecular dynamics runs of 25 picoseconds were performed, utilizing 184 experimental distance constraints and 80 torsion angle constraints; three different starting structures were used: energy minimized A-DNA, B-DNA, and wrinkled D-DNA, another member of the B-DNA family. Convergence to similar structures obtained with root-mean-square deviations between resulting structures of 0.37 to 0.92 A for the central hexamer of the octamer. The average structure from the nine different molecular dynamics runs was subjected to final restrained energy minimization. The resulting final structure was in good agreement with the structures derived from different molecular dynamics runs and exhibited a substantial improvement in the 2D NOE sixth-root residual index in comparison with the starting structures. An approximation of the structure in the terminal base-pairs, which displayed experimental evidence of fraying, was made by maintaining the structure of the inner four base-pairs and performing molecular dynamics simulations with the experimental structural constraints observed for the termini.(ABSTRACT TRUNCATED AT 400 WORDS)     

Metropolis Monte Carlo calculations of DNA structure using internal coordinates and NMR distance restraints: An alternative method for generating a high-resolution solution structure

Summary A new method, a restrained Monte Carlo (rMC) calculation, is demonstrated for generating high-resolution structures of DNA oligonucleotides in solution from interproton distance restraints and bounds derived from complete relaxation matrix analysis of two-dimensional nuclear Overhauser effect (NOE) spectral peak intensities. As in the case of restrained molecular dynamics (rMD) refinement of structures, the experimental distance restraints and bounds are incorporated as a pseudo-energy term (or penalty function) into the mathematical expression for the molecular energy. However, the use of generalized helical parameters, rather than Cartesian coordinates, to define DNA conformation increases efficiency by decreasing by an order of magnitude the number of parameters needed to describe a conformation and by simplifying the potential energy profile. The Metropolis Monte Carlo method is employed to simulate an annealing process. The rMC method was applied to experimental 2D NOE data from the octamer duplex d(GTA-TAATG)·d(CATTATAC). Using starting structures from different locations in conformational space (e.g. A-DNA and B-DNA), the rMC calculations readily converged, with a root-mean-square deviation (RMSD) of <0.3 Å between structures generated using different protocols and starting structures. Theoretical 2D NOE peak intensities were calculated for the rMC-generated structures using the complete relaxation matrix program CORMA, enabling a comparison with experimental intensities via residual indices. Simulation of the vicinal proton coupling constants was carried out for the structures generated, enabling a comparison with the experimental deoxyribose ring coupling constants, which were not utilized in the structure determination in the case of the rMC simulations. Agreement with experimental 2D NOE and scalar coupling data was good in all cases. The rMC structures are quite similar to that refined by a traditional restrained MD approach (RMSD<0.5 Å) despite the different force fields used and despite the fact that MD refinement was conducted with additional restraints imposed on the endocyclic torsion angles of deoxyriboses. The computational time required for the rMC and rMD calculations is about the same. A comparison of structural parameters is made and some limitations of both methods are discussed with regard to the average nature of the experimental restraints used in the refinement.Abbreviations MC Monte Carlo - rMC restrained Monte Carlo - MD molecular dynamics - rMD restrained molecular dynamics - DG distance geometry - EM energy minimization - 2D NOE two-dimensional nuclear Overhauser effect - DQF-COSY double-quantum-filtered correlation spectroscopy - RMSD root-mean-square deviation To whom correspondence should be addressed.     

Solution structure of the EcoRI DNA octamer containing 5-fluorouracil via restrained molecular dynamics using distance and torsion angle constraints extracted from NMR spectral simulations.

The self-complementary DNA octamer [d(GGAATUFCC)]2, containing the EcoRI recognition sequence with one of the thymines replaced by 5-fluorouracil (UF), was synthesized. Proton homonuclear two-dimensional nuclear Overhauser effect (2D NOE) and double-quantum-filtered correlation (2QF-COSY) spectra, as well as one-dimensional spectra at different temperatures, were recorded for the octamer. Consequently, all proton resonances were assigned. The thermally induced transition from the duplex to single strands has been followed, demonstrating the stability of the duplex containing 5-fluorouracil. Simulations of the 2QF-COSY cross-peaks by means of the programs SPHINX and LINSHA were compared with experimental data, establishing scalar coupling constants for the sugar ring protons and hence sugar pucker parameters. The deoxyribose rings exhibit a dynamic equilibrium of N- and S-type conformers with 75-95% populations of the latter. Two programs used for complete relaxation matrix analysis 2D NOE spectra, CORMA and MARDIGRAS, were modified to account for the influence of the fluorines on dipolar interactions in the proton system. Quantitative assessment of the 2D NOE cross-peak intensities for different mixing times, in conjunction with the program MARDIGRAS, gave a set of interproton distances for each mixing time. The largest and smallest values of each of the interproton distances were chosen as the upper and lower bounds for each distance constraint. The distance bounds define the size of a flat-well potential function term, incorporated into the AMBER force field, which was employed for restrained molecular dynamics calculations. Torsion angle constraints in the form of a flat-well potential were also constructed from the analysis of the sugar pucker data. Several restrained molecular dynamics runs of 35 ps were performed, utilizing 284 experimental distance and torsion angle constraints and two different starting structures, energy-minimized A- and B-DNA. Convergence to similar structures with a root-mean-square deviation of 1.2 A was achieved for the central hexamer of the octamer, starting from A- and B-DNA. The average structure from six different molecular dynamics runs was subjected to final restrained energy minimization. The resulting final structure was in good agreement with the structures derived from different molecular dynamics runs and showed a substantial improvement of the 2D NOE sixth-root residual index in comparison with classical and energy-minimized B-DNA. A detailed analysis of the conformation of the final structure and comparison with structures of similar sequences, obtained by different methods, were performed.     

Refinement of the NMR structures for acyl carrier protein with scalar coupling data

Structure determination of small proteins using NMR data is most commonly pursued by combining NOE derived distance constraints with inherent constraints based on chemical bonding. Ideally, one would make use of a variety of experimental observations, not just distance constraints. Here, coupling constant constraints have been added to molecular mechanics and molecular dynamics protocols for structure determination in the form of a psuedoenergy function that is minimized in a search for an optimum molecular conformation. Application is made to refinement of a structure for a 77 amino acid protein involved in fatty acid synthesis, Escherichia coli acyl carrier protein (ACP). 54 3JHN alpha coupling constants, 12 coupling constants for stereospecifically assigned side chain protons, and 450 NOE distance constraints were used to calculate the 3-D structure of ACP. A three-step protocol for a molecular dynamics calculation is described, in analogy to the protocol previously used in molecular mechanics calculations. The structures calculated with the molecular mechanics approach and the molecular dynamics approach using a rigid model for the protein show similar molecular energies and similar agreement with experimental distance and coupling constant constraints. The molecular dynamics approach shows some advantage in overcoming local minimum problems, but only when a two-state averaging model for the protein was used, did molecular energies drop significantly.     

Determination of cross relaxation rates by 3D NOE-NOE spectroscopy

Summary A method for quantitative determination of cross-relaxation rates of macromolecules in solution is developed. The method is based on the analysis of the intensities of cross peaks in 3D NOE-NOE spectra. The linear combination of the intensities of 3D peaks (spin-diffusion peaks, back-transfer peaks) results in an expression directly proportional to the cross-relaxation rate. The proposed approach allows to determine interproton distances in macromolecules more accurately.     

Solution structure of the TnAn DNA duplex GCCGTTAACGCG containing the HpaI restriction site.

The solution structure of the self-complementary DNA duplex [d(GCCGTTAACGGC)]2, which contains the HpaI restriction site GTTAAC, has been elucidated by two-dimensional NMR, distance geometry (DG), and NOE back-calculation methods. Initial distance constraints were determined by polynomial fitting the two-spin initial NOE rates; backbone constraints from NOE and J-coupling observations (Kim et al., 1992) were included. RMSDs between initial-distance-refined structures derived from random-embedded DG, A-DNA, and B-DNA starting structures were all in the range 0.5-1.0 A, indicating good convergence properties of the algorithm, regardless of the starting structure. A semiautomatic back-calculation refinement procedure was developed and used to generate more refined structures for which the BKCALC-simulated NOE volumes matched the experimental data. The six final structures refined from various starting structures exhibit very good agreement with the experimental data (R values = 0.18) and converge well to within 0.8-A RMSD differences for the central 8 base pairs. The torsion and pseudorotation phase angles were found to be well determined by the data, and the local helical parameters for each base step converged quite well. The final structures show that the central T6-A7 step is somewhat underwound (twist angle ca. 29 degrees), with a large negative cup and a normal (wide) minor groove width, while the T5-T6 and A7-A8 steps have a partially narrowed minor groove.     

“Ensemble” iterative relaxation matrix approach: A new NMR refinement protocol applied to the solution structure of crambin

The structure in solution of crambin, a small protein of 46 residues, has been determined from 2D NMR data using an iterative relaxation matrix approach (IRMA) together with distance geometry, distance bound driven dynamics, molecular dynamics, and energy minimization. A new protocol based on an “ensemble” approach is proposed and compared to the more standard initial rate analysis approach and a “single structure” relaxation matrix approach. The effects of fast local motions are included and R-factor calculations are performed on NOE build-ups to describe the quality of agreement between theory and experiment. A new method for stereospecific assignment of prochiral groups, based on a comparison of theoretical and experimental NOE intensities, has been applied. The solution structure of crambin could be determined with a precision (rmsd from the average structure) of 0.7 Å on backbone atoms and 1.1 Å on all heavy atoms and is largely similar to the crystal structure with a small difference observed in the position of the side chain of Tyr-29 which is determined in solution by both J-coupling and NOE data. Regions of higher structural variability (suggesting higher mobility) are found hi the solution structure, in particular for the loop between the two helices (Gly-20 to Pro-22). © 1993 Wiley-Liss, Inc.     

Combined procedure of distance geometry and restrained molecular dynamics techniques for protein structure determination from nuclear magnetic resonance data: application to the DNA binding domain of lac repressor from Escherichia coli

The technique of two-dimensional nuclear magnetic resonance (2D-NMR) has recently assumed an active role in obtaining information on structures of polypeptides, small proteins, sugars, and DNA fragments in solution. In order to generate spatial structures from the atom-atom distance information obtained by the NMR method, different procedures have been developed. Here we introduce a combined procedure of distance geometry (DG) and molecular dynamics (MD) calculations for generating 3D structures that are consistent with the NMR data set and have reasonable internal energies. We report the application of the combined procedure on the lac repressor DNA binding domain (headpiece) using a set of 169 NOE and 17 "hydrogen bond" distance constraints. Eight of ten structures generated by the distance geometry algorithm were refined within 10 ps MD simulation time to structures with low internal energies that satisfied the distance constraints. Although the combination of DG and MD was designed to combine the good sampling properties of the DG algorithm with an efficient method of lowering the internal energy of the molecule, we found that the MD algorithm contributes significantly to the sampling as well.     

The solution structure of a [3Fe-4S] ferredoxin: oxidised ferredoxin II from Desulfovibrio gigas

The use of standard 2D NMR experiments in combination with 1D NOE experiments allowed the assignment of 51 of the 58 spin systems of oxidised [3Fe-4S] ferredoxin isolated from Desulfovibrio gigas. The NMR solution structure was determined using data from 1D NOE and 2D NOESY spectra, as distance constraints, and information from the X-ray structure for the spin systems not detected by NMR in torsion angle dynamics calculations to produce a family of 15 low target function structures. The quality of the NMR family, as judged by the backbone r.m.s.d. values, was good (0.80?Å), with the majority of φ/ψ angles falling within the allowed region of the Ramachandran plot. A comparison with the X-ray structure indicated that the overall global fold is very similar in solution and in the solid state. The determination of the solution structure of ferredoxin II (FdII) in the oxidised state (FdII_ox) opens the way for the determination of the solution structure of the redox intermediate state of FdII (FdII_int), for which no X-ray structure is available.     

Systematic study of nuclear Overhauser effects vis-à-vis local helical parameters, sugar puckers, and glycosidic torsions in B DNA: insensitivity of NOE to local transitions in B DNA oligonucleotides due to internal structural compensations.

A method has been developed to solve structures of DNA oligomers in solution from the experimental NOE data. The method is a combination of two approaches: (1) full matrix NOESY simulations and (2) conformational calculations of DNA double helix based on generalized helical parameters. The process of the refinement of a solution structure does not involve NMR-derived interproton distance constraints; rather it consists of a direct fitting of a structure to the experimental NOE data, a weighted sum of energy, and R factor being under minimization. A helical parameters-based generation of DNA forms makes it possible to organize the search for the optimal structure more effectively, systematically varying starting conformations. The method has been used to calculate a structure for the self-complementary DNA hexamer GGATCC, which is consistent with the available experimental data. The structure belongs to the B family of forms, although the local structural heterogeneity is very strong. Sugar puckers vary from O4'-exo to C3'-exo; helical steps are open with different magnitudes toward the minor groove. Next, we have addressed the question of how uniquely the structure is defined by the existing NMR data. Different structural parameters have been systematically varied, and their effect on individual NOE's and the R factor has been studied. Two energetically conjugated parameters, sugar puckers and glycosidic angles, can be determined very reliably, because of the strong dependences of the intraresidue H6/H8 to H2'/H2'/H3' NOE's. In contrast, the local helical conformation of DNA and the geometry of base pairs proved to be underdetermined by the existing NOE information, because the effect of any helical parameter on interproton distances can be compensated by the concerted changes in other parameters.     

Structural homology of spinach acyl carrier protein and Escherichia coli acyl carrier protein based on NMR data

Acyl carrier proteins (ACPs) from spinach and from Escherichia coli have been used to demonstrate the utility of proton NMR for comparison of homologous structures. The structure of E. coli ACP had been previously determined and modeled as a rapid equilibrium among multiple conformational forms (Kim and Prestegard, Biochemistry 28:8792–8797, 1989). Spinach ACP showed two slowly exchanging forms and could be manipulated into one form for structural study. Here we compare this single form to postulated multiple forms of E. coli ACP using the limited amount of NOE data available for the spinach protein. A number of long-range NOE contacts were present between homologous residues in both spinach and E. coli ACP, suggesting tertiary structural homology. To allow a more definitive structural comparison, a method was developed to use spinach ACP NOE constraints to search for regions of structural divergence from two postulated forms of E. coli ACP. The homologous regions of the two protein sequences were aligned, additional distance constraints were extracted from the E. coli structure, and these were mapped onto the spinach sequence. These distance constraints were combined with experimental NOE constraints and a distance geometry simulated annealing protocol was used to test for compatibility of the constraints. All of the experimental spinach NOE constraints could be successfully combined with the E. coli data, confirming the general hypothesis of structural homology. A better fit was obtained with one form, suggesting a preferential stabilization of that form in the spinach case. Proteins 27:131–143 © 1997 Wiley-Liss, Inc.     

Intermolecular relaxation has little effect on intra-peptide exchange-transferred NOE intensities

Exchange-transferred nuclear Overhauser enhancement (etNOE) provides a useful method for determining the 3-dimensional structure of a ligand bound to a high-molecular-weight complex. Some concern about the accuracy of such structures has arisen because indirect relaxation can occur between the ligand and macromolecule. Such indirect relaxation, or spin diffusion, would lead to errors in interproton distances used as restraints in structure determination. We address this concern by assessing the extent of intermolecular spin diffusion in nineteen peptide-protein complexes of known structure. Transferred NOE intensities were simulated with the program CORONA (Calculated OR Observed NOESY Analysis) using the rate-matrix approach to include contributions from indirect relaxation between protein-ligand and intraligand proton pairs. Intermolecular spin diffusion contributions were determined by comparing intensities calculated with protonated protein to those calculated with fully deuterated protein. The differences were found to be insignificant overall, and to diminish at short mixing times and high mole ratios of peptide to protein. Spin diffusion between the peptide ligand and the protein contributes less to the etNOE intensities and alters fewer cross peaks than the well-studied intramolecular spin diffusion effects. Errors in intraligand interproton distances due to intermolecular relaxation effects were small on average and can be accounted for with the restraint functions commonly used in NMR structure determination methods. In addition, a rate-matrix approach to calculate distances from etNOESY intensities using a volume matrix comprising only intraligand intensities was found to give accurate values. Based on these results, we conclude that structures determined from etNOESY data are no less accurate due to spin diffusion than structures determined from conventional NOE intensities.     

Complete relaxation matrix refinement of NMR structures of proteins using analytically calculated dihedral angle derivatives of NOE intensities

Summary A new method for refining three-dimensional (3D) NMR structures of proteins is described, which takes account of the complete relaxation pathways. Derivatives of the NOE intensities with respect to the dihedral angles are analytically calculated, and efficiently evaluated with the use of a filter technique for identifying the dominant terms of these derivatives. This new method was implemented in the distance geometry program DIANA. As an initial test, we refined 30 rigid distorted helical structures, using a simulated data set of NOE distance constraints for a rigid standard -helix. The final root-mean-square deviations of the refined structures relative to the standard helix were less than 0.1 Å, and the R-factors dropped from values between 7% and 32% to values of less than 0.5% in all cases, which compares favorably with the results from distance geometry calculations. In particular, because spin diffusion was not explicitly considered in the evaluation of exact¹H–¹H distances corresponding to the simulated NOE intensities, a group of nearly identical distance geometry structures was obtained which had about 0.5 Å root-mean-square deviation from the standard -helix. Further test calculations using an experimental NOE data set recorded for the protein trypsin inhibitor K showed that the complete relaxation matrix refinement procedure in the DIANA program is functional also with systems of practical interest.Abbreviations RMSD root-mean-square deviation - NOE nuclear Overhauser enhancement - NOESY 2-dimensional nuclear Overhauser enhancement spectroscopy - CPU central processing unit     

Solution conformations of GM3 gangliosides containing different sialic acid residues as revealed by NOE-based distance mapping, molecular mechanics, and molecular dynamics calculations.

The conformation of the GM3 ganglioside, Neu5Ac alpha 2-3Gal beta 1-4Glc beta 1-1 Cer, and its analogs containing the Neu5Gc or Neu4Ac5Gc residues (Gc = glycolyl, CH2OHCO) was investigated in Me2SO-d6 solution with the aid of a distance-mapping procedure based on rotating-frame NOE contacts, with hydroxyl protons being used as long-range sensors defining the distance constraints. A pronounced flexibility found for both the Neu-Gal and Gal-Glc linkages was confirmed by 1000-ps molecular dynamics simulations. Similar results, although based on a smaller number of NOE constraints, were obtained for GM3 gangliosides anchored in mixed D2O/dodecylphosphocholine-d38 micelles and for the Neu5Ac-, Neu5Gc-, and Neu5,9Ac2-sialyllactoses dissolved in D2O. No noteworthy differences in conformational behavior of the glycan chains of the three gangliosides or sialyllactoses were observed in either of the media.