A combined 2D-NMR and molecular dynamics analysis of the structure of the actinomycin D: d(ATGCAT)2 complex

We present a comparative analysis of an NMR experiment and molecular and harmonic dynamics simulations of an actinomycin D: d(ATGCAT)2 complex. A comparison of NOE measurements and 1/R6 weighted proton-proton distances confirm the general correctness of the Actinomycin D-DNA model proposed by Sobell. There are, however, some substantial differences between the proton-proton distances inferred from the NOE results and the molecular and harmonic dynamics simulations. The remaining discrepancies could either come from contributions of other conformations to the average properties of the complex or from uncertainties in the NMR distance analysis. An analysis of the molecular dynamics helix properties, sugar puckers, hydrogen bonding, rms fluctuations and torsional properties are qualitatively consistent with those from previous simulations, but the presence of an intercalated drug leads to some new structural and dynamical features.     

Structural Studies of 5-Ethyl-2′-Deoxyuridine by Selective Pulse 1H DPFGSE NOE Spectroscopy and PM3 Calculations

Abstract

The solution conformation of 5-ethyl-2′-deoxyuridine (EDU) has been calculated from the vicinyl proton-proton NMR coupling constants and nuclear Overhauser (NOE) distances using excitation sculpting of selective pulses (Double Pulsed Field Gradient Spin Echo NOE) at 500 MHz and molecular modelling (PM3) studies.     

Conformational investigation of a cyclic enterobacterial common antigen employing NMR spectroscopy and molecular dynamics simulations

The three-dimensional structure of a cyclic enterobacterial common antigen (ECA) having four trisaccharide repeating units has been investigated by NMR spectroscopy and molecular dynamics simulations. Three different NMR parameters were determined: (a) (1)H,(1)H cross-relaxation rates from NOE experiments were used for determination of proton-proton distances; (b) trans-glycosidic (3)J(C,H) scalar coupling constants analyzed via a Karplus-type relationship provided information on torsion angles; and (c) (1)H,(13)C one-bond dipolar couplings obtained in a dilute liquid-crystalline medium were interpreted in terms of the orientational order and molecular conformations. The molecular dynamics simulations of the dodecasaccharide were performed with explicit water and counterions, which are important factors that strongly influence molecular conformation. Subsequently, the results from computer simulation were used to generate a three-dimensional structure of the cyclic ECA which is consistent with the experimental NMR parameters.     

A two-dimensional nuclear Overhauser enhancement (2D NOE) experiment for the elucidation of complete proton-proton cross-relaxation networks in biological macromolecules

The recently developed technique of two-dimensional (2D) cross-relaxation spectroscopy is utilized for systematic measurements of selective nuclear Overhauser enhancements (NOE) in the high resolution ¹H nuclear magnetic resonance (NMR) spectra of biological macromolecules in solution. Compared to conventional one-dimensional NOE studies, the 2D NOE experiment has the principal advantage that it avoids detrimental effects arising from the limited selectivity of preirradiation in crowded spectral regions. Furthermore, it yields with a single instrument setting a complete network of NOE's between all the protons in the macromolecule. The resulting information on intramolecular proton-proton distances provides a new avenue for studies of the spatial structures of biopolymers.     

Structural studies of 5-ethyl-2'-deoxyuridine by selective pulse 1H DPFGSE NOE spectroscopy and PM3 calculations.

The solution conformation of 5-ethyl-2'-deoxyuridine (EDU) has been calculated from the vicinyl proton-proton NMR coupling constants and nuclear Overhauser (NOE) distances using excitation sculpting of selective pulses (Double Pulsed Field Gradient Spin Echo NOE) at 500 MHz and molecular modelling (PM3) studies.     

Validation of the GROMOS force-field parameter set 45A3 against nuclear magnetic resonance data of hen egg lysozyme

The quality of molecular dynamics (MD) simulations of proteins depends critically on the biomolecular force field that is used. Such force fields are defined by force-field parameter sets, which are generally determined and improved through calibration of properties of small molecules against experimental or theoretical data. By application to large molecules such as proteins, a new force-field parameter set can be validated. We report two 3.5 ns molecular dynamics simulations of hen egg white lysozyme in water applying the widely used GROMOS force-field parameter set 43A1 and a new set 45A3. The two MD ensembles are evaluated against NMR spectroscopic data NOE atom–atom distance bounds, ³J_NH and ³J coupling constants, and ¹5N relaxation data. It is shown that the two sets reproduce structural properties about equally well. The 45A3 ensemble fulfills the atom–atom distance bounds derived from NMR spectroscopy slightly less well than the 43A1 ensemble, with most of the NOE distance violations in both ensembles involving residues located in loops or flexible regions of the protein. Convergence patterns are very similar in both simulations atom-positional root-mean-square differences (RMSD) with respect to the X-ray and NMR model structures and NOE inter-proton distances converge within 1.0–1.5 ns while backbone ³J_HN-coupling constants and ¹H– ¹5N order parameters take slightly longer, 1.0–2.0 ns. As expected, side-chain ³J-coupling constants and ¹H– ¹5N order parameters do not reach full convergence for all residues in the time period simulated. This is particularly noticeable for side chains which display rare structural transitions. When comparing each simulation trajectory with an older and a newer set of experimental NOE data on lysozyme, it is found that the newer, larger, set of experimental data agrees as well with each of the simulations. In other words, the experimental data converged towards the theoretical result.     

Structure refinement of a cyclic peptide from two-dimensional NMR data and molecular modeling

The conformational and dynamic properties of a cyclic peptide designed to inhibit human renin have been examined by using NMR and molecular modeling. From a quantitative analysis of a series of two-dimensional NOE data sets, proton-proton distances were calculated. Several different methods were explored and compared to incorporate these distance constraints as well as those derived from vicinal spin-spin coupling constants into computer-generated three-dimensional structures. These methods included interactive manual manipulation of the structures to fit the NMR-determined distance constraints, distance geometry, constrained energy minimizations, and constrained molecular dynamics. The advantages and disadvantages of the methods are discussed. In addition, to gain insight into the conformations accessible to the cyclic peptide and the relative flexibility of the different parts of the molecule, molecular dynamics calculations were performed at three different temperatures. Average interproton distances and dihedral angles were obtained from the structures generated in the dynamics trajectories and compared to those obtained from the NMR experiments. Despite the four methylene groups and ether linkage contained in the cyclic portion of the peptide, our NMR results indicated a preferred conformation for the macrocyclic ring of the peptide and supported the presence of a cis Phe-Ala peptide bond. In contrast, both the molecular dynamics and NMR data indicated a considerable amount of flexibility for the remaining noncyclic portion of the molecule. These results are used to propose an explanation for the cyclic peptide's inability to inhibit human renin.     

Methods of NMR structure refinement: molecular dynamics simulations improve the agreement with measured NMR data of a C-terminal peptide of GCN4-p1

The C-terminal trigger sequence is essential in the coiled-coil formation of GCN4-p1; its conformational properties are thus of importance for understanding this process at the atomic level. A solution NMR model structure of a peptide, GCN4p16–31, encompassing the GCN4-p1 trigger sequence was proposed a few years ago. Derived using a standard single-structure refinement protocol based on 172 nuclear Overhauser effect (NOE) distance restraints, 14 hydrogen-bond and 11 ϕ torsional-angle restraints, the resulting set of 20 NMR model structures exhibits regular α-helical structure. However, the set slightly violates some measured NOE bounds and does not reproduce all 15 measured ³J(H_N-H_Cα)-coupling constants, indicating that different conformers of GCN4p16–31 might be present in solution. With the aim to resolve structures compatible with all NOE upper distance bounds and ³J-coupling constants, we executed several structure refinement protocols employing unrestrained and restrained molecular dynamics (MD) simulations with two force fields. We find that only configurational ensembles obtained by applying simultaneously time-averaged NOE distance and ³J-coupling constant restraining with either force field reproduce all the experimental data. Additionally, analyses of the simulated ensembles show that the conformational variability of GCN4p16–31 in solution admitted by the available set of 187 measured NMR data is larger than represented by the set of the NMR model structures. The conformations of GCN4p16–31 in solution differ in the orientation not only of the side-chains but also of the backbone. The inconsistencies between the NMR model structures and the measured NMR data are due to the neglect of averaging effects and the inclusion of hydrogen-bond and torsional-angle restraints that have little basis in the primary, i.e. measured NMR data.     

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.     

Solution structure by nuclear magnetic resonance of the two lantibiotics 97518 and NAI‐107

Lantibiotics 97518 and NAI‐107, produced by the related genera Planomonospora and Microbispora respectively, are members of a family of nisin‐related compounds. They represent promising compounds to treat infections caused by multiresistant Gram‐positive pathogens. Despite their similar structure and a similar antibacterial spectrum, the two lantibiotics exhibit significant differences in their potency. To gain an insight into the structure–activity relationships, their conformational properties in solution are determined by NMR. After carrying out an NOE analysis of 2D ¹H NMR spectra, high‐resolution 3D structures are determined using molecular dynamics simulations. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Molecular Dynamics and NMR Spectroscopy Studies of E.?coli Lipopolysaccharide Structure and Dynamics

Lipopolysaccharide (LPS), a component of Gram-negative bacterial outer membranes, comprises three regions: lipid A, core oligosaccharide, and O-antigen polysaccharide. Using the CHARMM36 lipid and carbohydrate force fields, we have constructed a model of an Escherichia coli R1 (core) O6 (antigen) LPS molecule. Several all-atom bilayers are built and simulated with lipid A only (LIPA) and varying lengths of 0 (LPS0), 5 (LPS5), and 10 (LPS10) O6 antigen repeating units; a single unit of O6 antigen contains five sugar residues. From ¹H,¹H-NOESY experiments, cross-relaxation rates are obtained from an O-antigen polysaccharide sample. Although some experimental deviations are due to spin-diffusion, the remaining effective proton-proton distances show generally very good agreement between NMR experiments and molecular dynamics simulations. The simulation results show that increasing the LPS molecular length has an impact on LPS structure and dynamics and also on LPS bilayer properties. Terminal residues in a LPS bilayer are more flexible and extended along the membrane normal. As the core and O-antigen are added, per-lipid area increases and lipid bilayer order decreases. In addition, results from mixed LPS0/5 and LPS0/10 bilayer simulations show that the LPS O-antigen conformations at a higher concentration of LPS5 and LPS10 are more orthogonal to the membrane and less flexible. The O-antigen concentration of mixed LPS bilayers does not have a significant effect on per-lipid area and hydrophobic thickness. Analysis of ion and water penetration shows that water molecules can penetrate inside the inner core region, and hydration is critical to maintain the integrity of the bilayer structure.     

Molecular Dynamics and NMR Spectroscopy Studies of E. coli Lipopolysaccharide Structure and Dynamics

Lipopolysaccharide (LPS), a component of Gram-negative bacterial outer membranes, comprises three regions: lipid A, core oligosaccharide, and O-antigen polysaccharide. Using the CHARMM36 lipid and carbohydrate force fields, we have constructed a model of an Escherichia coli R1 (core) O6 (antigen) LPS molecule. Several all-atom bilayers are built and simulated with lipid A only (LIPA) and varying lengths of 0 (LPS0), 5 (LPS5), and 10 (LPS10) O6 antigen repeating units; a single unit of O6 antigen contains five sugar residues. From ¹H,¹H-NOESY experiments, cross-relaxation rates are obtained from an O-antigen polysaccharide sample. Although some experimental deviations are due to spin-diffusion, the remaining effective proton-proton distances show generally very good agreement between NMR experiments and molecular dynamics simulations. The simulation results show that increasing the LPS molecular length has an impact on LPS structure and dynamics and also on LPS bilayer properties. Terminal residues in a LPS bilayer are more flexible and extended along the membrane normal. As the core and O-antigen are added, per-lipid area increases and lipid bilayer order decreases. In addition, results from mixed LPS0/5 and LPS0/10 bilayer simulations show that the LPS O-antigen conformations at a higher concentration of LPS5 and LPS10 are more orthogonal to the membrane and less flexible. The O-antigen concentration of mixed LPS bilayers does not have a significant effect on per-lipid area and hydrophobic thickness. Analysis of ion and water penetration shows that water molecules can penetrate inside the inner core region, and hydration is critical to maintain the integrity of the bilayer structure.     

Structure determination of a pseudotripeptide zinc complex with the COSMOS-NMR force field and DFT methods

A His-X-His pseudotripeptide zinc complex (X is a N-alkyl glycine derivative) similar to the catalytic center of the carbonic anhydrase was computer designed and experimentally synthesized. Using 2D-NMR techniques, all proton, carbon chemical shifts and nuclear overhauser effect signals were assigned. The three-dimensional structure of the complex was determined with the COSMOS (computer simulation of molecular structures) force field by applying ¹³C bond polarization theory chemical shift pseudo forces and restrictions for NOE distances. From molecular dynamics, simulated annealing simulations and geometry optimizations, the three best force field structures were taken for a final investigation by density functional theory calculations.     

Simulating structural and thermodynamic properties of carcinogen-damaged DNA

A pair of stereoisomeric covalent adducts to guanine in double-stranded DNA, derived from the reaction of mutagenic and tumorigenic metabolites of benzo[a]pyrene, have been well characterized structurally and thermodynamically. Both high-resolution NMR solution structures and an array of thermodynamic data are available for these 10S (+)- and 10R (-)-trans-anti -[BP]-N(2)-dG adducts in double-stranded deoxyoligonucleotides. The availability of experimentally well-characterized duplexes containing these two stereoisomeric guanine adducts provides an opportunity for evaluating the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method for computing thermodynamic properties from molecular dynamics ensembles. We have carried out 3-ns molecular dynamics simulations, using NMR solution structures as the starting models for the 10S (+)- and 10R (-)-trans-anti-dG adducts in a DNA duplex 11-mer using AMBER 6.0. We employed the MM-PBSA method to compute the free energies, enthalpies, and entropies of the two adducts. Our complete thermodynamic analysis agrees quite well with the full experimental thermodynamic characterization of these adducts, showing essentially equal stabilities of the two adducts. We also calculated the nuclear Overhauser effect (NOE) distances from the molecular dynamics trajectories, and compared them against the experimental NMR-derived NOE distances. Our results showed that the simulated structures are in good agreement with the NMR experimental NOE data. Furthermore, the molecular dynamics simulations provided new structural and biological insights. Specifically, the puzzling observation that the BP aromatic ring system in the 10S (+)-trans-anti-dG adduct is more exposed to the aqueous solvent than the 10R (-)-trans-anti-dG adduct, is rationalized in terms of the adduct structures. The structural and thermodynamic features of these stereoisomeric adducts are also discussed in relation to their reported low susceptibilities to nucleotide excision repair.     

NMR studies of the conformation of the natural sweetener rebaudioside A

Rebaudioside A is a natural sweetener from Stevia rebaudiana in which four β-d-glucopyranose units are attached to the aglycone steviol. Its ¹H and ¹³C NMR spectra in pyridine-d₅ were assigned using 1D and 2D methods. Constrained molecular dynamics of solvated rebaudioside using NMR constraints derived from ROESY cross peaks yielded the orientation of the β-d-glucopyranose units. Hydrogen bonding was examined using the temperature coefficients of the hydroxyl chemical shifts, ROESY and long-range COSY spectra, and proton-proton coupling constants.     

On the similarity of properties in solution or in the crystalline state: A molecular dynamics study of hen lysozyme

As protein crystals generally possess a high water content, it is assumed that the behaviour of a protein in solution and in crystal environment is very similar. This assumption can be investigated by molecular dynamics (MD) simulation of proteins in the different environments. Two 2ns simulations of hen egg white lysozyme (HEWL) in crystal and solution environment are compared to one another and to experimental data derived from both X-ray and NMR experiments, such as crystallographic B-factors, NOE atom–atom distance bounds, ³J_{H N}-coupling constants, and ¹H-¹⁵N bond vector order parameters. Both MD simulations give very similar results. The crystal simulation reproduces X-ray and NMR data slightly better than the solution simulation.     

Characterization of the internal dynamics and conformational space of zinc-bound amyloid β peptides by replica-exchange molecular dynamics simulations

Amyloid β (Aβ) peptides and metal ions have been associated with the pathogenesis of Alzheimer’s disease. The conformational space of Aβ fragments of different length with and without binding of metal ions has been extensively investigated by replica-exchange molecular dynamics (REMD) simulation. However, only trajectories extracted at relatively low temperatures have been used for this analysis. The capability of REMD simulations to characterize the internal dynamics of such intrinsically disordered proteins (IDPs) as Aβ has been overlooked. In this work, we use an approach recently developed by Xue and Skrynnikov (J Am Chem Soc 133:14614–14628, 2011) to calculate NMR observables, including ¹⁵N relaxation rates and ¹⁵N–¹H nuclear Overhauser enhancement (NOE), from the high-temperature trajectory of REMD simulations for zinc-bound Aβ peptides. The time axis of the trajectory was rescaled to correct for the effect of the high temperature (408 K) compared with the experimental temperature (278 K). Near-quantitative agreement between simulated values and experimental results was obtained. When the structural properties and free-energy surfaces of zinc-bound Aβ_(1–40) and Aβ_(1–42) were compared at the physiological temperature 310 K it was found that zinc-bound Aβ_(1–42) was more rigid than Aβ_(1–40) at the C terminus, and its conformational transitions were also more preferred. The self-consistent results derived from trajectories at high and low temperatures demonstrate the capability of REMD simulations to capture the internal dynamics of IDPs.     

Structural Basis for the Binding Affinity of a Homologous Series of Synthetic Phenoxazone Drugs with DNA: NMR and Molecular Mechanics Analysis

Abstract

The molecular basis of the marked structure-activity relationship for a homologous series of DNA-binding phenoxazone drugs (ActII-ActIV) has been investigated by NMR spectroscopy and molecular mechanics. The spatial structures of the complexes between the drugs and a model deoxytetranucleotide, 5′-d(TpGpCpA), have been determined by molecular mechanics methods using homonuclear ¹H-¹H 2D-NOESY and heteronuclear ¹H-³¹P (HMBC) NMR spectroscopic data. Observed intermolecular NOE contacts and equilibrium binding studies confirm that the binding affinity of the synthetic phenoxazone derivatives with d(TGCA) decreases with an increase in the number of CH₂ groups in the dimethylami- noalkyl side chains, i.e., ActII > ActIII > ActIV, in agreement with the observed biological activity of these compounds. Molecular mechanics calculations of the spatial structures of the intercalated complexes of ActII-ActIV with d(TGCA) indicate that the different binding constants of the phenoxazone derivatives with the DNA oligomer are due to the different degrees of intercalation of the chromophore and the different steric arrangements of aminoalkyl side chains in the minor groove of the tetramer duplex; this results in different distances between the negatively-charged phosphates of the DNA duplex and the terminal positively-charged N(CH₃)₂ groups of the side chains.     

Salt induced deprotonation of initially salt-free transfer RNA.

¹H NMR double resonance studies of valinomycin in (CD₃)₂SO were conducted at 90 MHz (FT-mode) and 250 MHz (correlation mode) to determine the mechanism of intramolecular nuclear Overhauser effects (NOE). These studies set specific constraints on any model for the conformation of valinomycin in (CD₃)₂SO and illustrate that NOE experiments of this type hold considerable promise for conformational studies of peptides, proteins and other biomolecules. The NOE's are positive at the lower frequency and negative at the higher frequency. Consideration of the theoretical dependence of the NOE on the proton-proton internuclear correlation time and on the resonance frequency indicates that these results are explained by a predominantly dipolar relaxation mechanism.     

Determination of preferred conformations of brassinosteroids by means of NMR investigations and Boltzmann statistical analysis of simulated annealing calculations

Brassinosteroids are a class of steroidal phytohormones with high growth-promoting properties. The preferred side-chain conformations of 10 brassinosteroids were determined by means of detailed NMR investigations and molecular-modeling studies. Vacuum conformations obtained by simulated annealing calculations and Boltzmann statistical analysis were compared with solution conformations derived from NOE experiments and molecular dynamic simulations, and with X-ray structures. In general, results from simulated annealing calculations and NMR-supported molecular dynamics simulations are in good agreement. For some of the compounds investigated the conformation was less well-defined at the end of the side-chain. It could be shown that the energetically most favorable and most probable conformations also include the conformations obtained by NMR supported molecular-dynamics calculations and by X-ray analysis. For the most bioactive compound brassinolide (1) the majority of conformations show a side-chain bent towards the -face of the steroid skeleton, whereas for the less bioactive brassinosteroids, conformations with straight side-chains or side-chains bent towards the -face are more frequent.Electronic Supplementary Material available.