Dynamic exchange between stabilized conformations predicted for hyaluronan tetrasaccharides: comparison of molecular dynamics simulations with available NMR data

Studies of the hyaluronan (HA) tetrasaccharides are important for understanding hydrogen-bonding in the HA polymer, as they are probably the smallest oligomers in which characteristics of the constituent monosaccharides and the polymer are simultaneously exhibited. Here we present extensive molecular dynamics simulations of the two tetrasaccharides of HA in dilute aqueous solution. These simulations have confirmed the existence of intramolecular hydrogen-bonds between the neighboring sugar residues of HA in solution, as proposed by Scott (1989). However, our simulations predict that these intramolecular hydrogen-bonds are not static as previously proposed, but are in constant dynamic exchange on the sub-nanosecond time-scale. This process results in discrete internal motion of the HA tetrasaccharides where they rapidly move between low energy conformations. Specific interactions between water and intramolecular hydrogen-bonds involving the hydroxymethyl group were found to result in differing conformations and dynamics for the two alternative tetrasaccharides of HA. This new observation suggests that this residue may play a key role in the entropy and stability of HA in solution, allowing it to stay soluble up to high concentration. The vicinal coupling constants3 J NHCH of the acetamido groups have been calculated from our aqueous simulations of HA. We found that high values of 3J NHCH approximately 8 Hz, as experimentally measured for HA, are consistent with mixtures of both trans and cis conformations, and thus3 J NHCH cannot be used to imply a purely trans conformation of the acetamido. The rapid exchange of intramolecular hydrogen-bonds indicates that although the structure is at any moment stabilized by these hydrogen-bonds, no one hydrogen-bond exists for an extended period of time. This could explain why NMR often fails to provide evidence for intramolecular hydrogen-bonds in HA and other aqueous carbohydrate structures.      

On the dependence of molecular conformation on the type of solvent environment: a molecular dynamics study of cyclosporin A

The dependence of the conformation of cyclosporin A (CPA), a cyclic undecapeptide with potent immunosuppressive activity, on the type of solvent environment is examined using the computer simulation method of molecular dynamics (MD). Conformational and dynamic properties of CPA in aqueous solution are obtained from MD simulations of a CPA molecule dissolved in a box with water molecules. Corresponding properties of CPA in apolar solution are obtained from MD simulations of CPA in a box with carbontetrachloride. The results of these simulations in H2O and in CCl4 are compared to each other and to those of previous simulations of crystalline CPA and of an isolated CPA molecule. The conformation of the backbone of the cyclic polypeptide is basically independent of the type of solvent. In aqueous solution the beta-pleated sheet is slightly weaker and the gamma-turn is a bit less pronounced than in apolar solution. Side chains may adopt different conformations in different solvents. In apolar solution the hydrophobic side chain of the MeBmt residue is in an extended conformation with its hydroxyl group hydrogen bonded to the backbone carbonyl group. In aqueous solution this hydrophobic side chain folds over the core of the molecule and the mentioned hydrogen bond is broken in favor of hydrogen bonding to water molecules. The conformation obtained from the MD simulation in CCl4 nicely agrees with experimental atom-atom distance data as obtained from nmr experiments in chloroform. In aqueous solution the relaxation of atomic motion tends to be slower than in apolar solution.     

The crystal and molecular structures of cellulose I and II

The paper describes molecular dynamics (MD) simulations on the crystal structures of the Iβ and II phases of cellulose. Structural proposals for each of these were made in the 1970s on the basis of X-ray diffraction data. However, due to the limited resolution of these data some controversies remained and details on hydrogen bonding could not be directly obtained. In contrast to structure factor amplitudes in X-ray diffraction, energies, as obtained from MD simulations, are very sensitive to the positions of the hydroxyl hydrogen atoms. Therefore the latter technique is very suitable for obtaining such structural details. MD simulations of the Iβ phase clearly shows preference for one of the two possible models in which the chains are packed in a parallel orientation. Only the parallel-down mode (in the definition of Gardner and Blackwell (1974) J Biopolym 13: 1975-2001) presents a stable structure. The hydrogen bonding consists of two intramolecular hydrogen bonds parallel to the glycosidic linkage for both chains, and two intralayer hydrogen bonds. The layers are packed hydrophobically. All hydroxymethyl group are positioned in the tg conformation. For the cellulose II form it was found that, in contrast to what seemed to emerge from the X-ray fibre diffraction data, both independent chains had the gt conformation. This idea already existed because of elastic moduli calculations and 13C-solid state NMR data. Recently, the structure of cellotetraose was determined. There appear to be a striking similarity between the structure obtained from the MD simulations and this cellotetraose structure in terms of packing of the two independent molecules, the hydrogen bonding network and the conformations of the hydroxymethyl group, which were also gt for both molecules. The structure forms a 3D hydrogen bonded network, and the contribution from electrostatics to the packing is more pronounced than in case of the Iβ structure. In contrast to what is expected, in view of the irreversible transition of the cellulose I to II form, the energies of the Iβ form is found to be lower than that of II by 1 kcal mol-1 per cellobiose. This revised version was published online in November 2006 with corrections to the Cover Date.     

Conformational study of cyclolinopeptide A. A distance geometry and molecular dynamics approach

The conformation of cyclolinopeptide A, c(Pro-Pro-Phe-Phe-Leu-Ile-Ile-Leu-Val), a naturally occurring peptide with remarkable cytoprotective activity, has been investigated by means of distance geometry calculations and molecular dynamics simulations. The starting points for all the calculations were an X-ray structure and other structures obtained from distance geometry calculations based on NMR data. Restrained and unrestrained molecular dynamics simulations are reported in vacuo and in CCl4. Structural and dynamic properties are investigated and compared with those experimentally determined. The conformation obtained from the MD simulations which best reproduces the NMR parameters is at the same time one of the most stable ones and is also fairly similar to the crystal structure. An explanation for the occurrence of multiple conformations in solution at room temperature is given.     

Structure and dynamics of Mycobacterium tuberculosis truncated hemoglobin N: insights from NMR spectroscopy and molecular dynamics simulations

The potent nitric oxide dioxygenase (NOD) activity (trHbN-Fe2?-O? + (?)NO → trHbN-Fe3?-OH? + NO??) of Mycobacterium tuberculosis truncated hemoglobin N (trHbN) protects aerobic respiration from inhibition by (?)NO. The high activity of trHbN has been attributed in part to the presence of numerous short-lived hydrophobic cavities that allow partition and diffusion of the gaseous substrates (?)NO and O? to the active site. We investigated the relation between these cavities and the dynamics of the protein using solution NMR spectroscopy and molecular dynamics (MD). Results from both approaches indicate that the protein is mainly rigid with very limited motions of the backbone N-H bond vectors on the picoseconds-nanoseconds time scale, indicating that substrate diffusion and partition within trHbN may be controlled by side-chains movements. Model-free analysis also revealed the presence of slow motions (microseconds-milliseconds), not observed in MD simulations, for many residues located in helices B and G including the distal heme pocket Tyr33(B10). All currently known crystal structures and molecular dynamics data of truncated hemoglobins with the so-called pre-A N-terminal extension suggest a stable α-helical conformation that extends in solution. Moreover, a recent study attributed a crucial role to the pre-A helix for NOD activity. However, solution NMR data clearly show that in near-physiological conditions these residues do not adopt an α-helical conformation and are significantly disordered and that the helical conformation seen in crystal structures is likely induced by crystal contacts. Although this lack of order for the pre-A does not disagree with an important functional role for these residues, our data show that one should not assume an helical conformation for these residues in any functional interpretation. Moreover, future molecular dynamics simulations should not use an initial α-helical conformation for these residues in order to avoid a bias based on an erroneous initial structure for the N-termini residues. This work constitutes the first study of a truncated hemoglobin dynamics performed by solution heteronuclear relaxation NMR spectroscopy.     

The conformational free-energy map for solvated neocarrabiose

A Ramachandran map of the conformational potential of mean force (pmf) for neocarrabiose in water was obtained using molecular dynamics (MD) simulations with umbrella sampling. The potential energy map calculated in a previous study for this molecule in vacuum exhibited a global minimum located at (phi = 81 degrees, psi = -141 degrees). However, the global minimum on the new pmf map in aqueous solution is located in an area centered around (phi = 175 degrees, psi = 180 degrees), indicating a considerable solvent shift. This new global minimum-energy solution conformation was found to correspond to the experimental value obtained from NMR-NOE measurements, and is also consistent with the experimental crystal structure for neocarrabiose and the fiber diffraction conformation for iota-carrageenan. The global minimum of the solution pmf and its local topology were found to be approximately reproduced by quick vacuum conformational energy mapping using several approximations that mimic solvation effects by de-emphasizing intramolecular hydrogen bonding.     

Biological properties of hyaluronan in aqueous solution are controlled and sequestered by reversible tertiary structures,defined by NMR spectroscopy

Hyaluronan (HA) is a ubiquitous polysaccharide of (predominantly) animal origin that has important medical applications in joint, skin, and eye conditions. Biological activities shown by HA fragments in angiogenesis, inflammation, etc. are absent from highly polymerized HA. We propose that HA physiological properties are controlled by molecular mass dependent transitions between tertiary structures (e.g., beta sheets) and 2-fold helices, - reversible "denaturation", which is central to HA solution behavior. We demonstrate this phenomenon by 13C NMR. Four different acetamido C=O resonances, assigned to secondary, tertiary, and disordered HA structures, monitored "denaturation" by (a) warming, (b) alkalinizing to pH >12.0, (c) hyaluronidase digestion, and (d) methylation of carboxylates. (a) and (b) acted reversibly but (c) and (d) are irreversible. 1H NMR implicated H-bonded acetamido NH in (b). Temperature dependencies of other 13C chemical shifts were small and unspecific. Arrhenius plots indicate that hyaluronan tertiary structures are on the edge of instability under physiological conditions. The results help to explain the appearance of biological activities on "denaturation" or degradation of HA.     

Protein conformational dynamics of crambin in crystal, solution and in the trajectories of molecular dynamics simulations

Atomic displacement parameters — B factors of the eight crambin crystal structures obtained at 0.54–1.5 Å resolution and temperatures of 100–293 K have been analyzed. The comparable contributions to the B factor values are the intramolecular motions which are modeled by the harmonic vibration calculations and derived from the molecular dynamics simulation (MD) as well as rigid body changes in the position of a protein molecule as a whole. In solution for the average NMR structure of crambin the amplitudes of the backbone atomic fluctuations of the most residues of the segments with the regular backbone conformations are close to the amplitudes of the small scale harmonic vibrations. For the same residues the probability of the medium scale fluctuations fixed by the hydrogen exchange method is very low. The restricted conformational mobility of those segments is coupled with the depressed amplitudes of the fluctuation changes of the tertiary structure registered by the residue accessibility changes in an ensemble of NMR structures that forms the average NMR structure of crambin. The amplitudes of temperature fluctuations of backbone atoms and the tertiary structure raise in the segment with the irregular conformations, turn and loops. In the same segments the amplitudes of the calculated harmonic vibrations also increase, but to a lesser extent and especially in the interhelical loop with the most strong and complicated fluctuation changes of the backbone conformation. In solution for the NMR structure in this loop the conformational transitions occur between the conformational substates separated by the energy barriers, but they are not observed even in the long 100 ns trajectories from the MD simulation of crambin. These strong local fluctuation changes of the structure may play a key role in the protein functioning and modern performance improvements in the MD simulation techniques are oriented to increase the probability of protein appearance in the trajectories from the MD simulations.     

On the Occurrence of Three-Center Hydrogen Bonds in Cyclodextrins in Crystalline Form and in Aqueous Solution: Comparison of Neutron Diffraction and Molecular Dynamics Results

Abstract

Three-center (bifurcated) hydrogen bonds may play a role by serving as an intermediate state between different dynamically changing hydrogen bonding patterns. Hydrogen bonding configurations can be studied experimentally by neutron diffraction and theoretically by computer simulation techniques. Here, both methods are used to analyse the occurrence of three-center hydrogen bonds in crystals of cyclodextrins.

Almost all experimentally observed three-center hydrogen bonds in the crystal are reproduced in the molecular dynamics (MD) simulations, even as far as the detailed asymmetric geometry is concerned. On the basis of this result a MD simulation of cyclodextrin in aqueous solution is searched for the occurrence of three-center hydrogen bonds. Significant differences are found. In solution more different three-center hydrogen bonds per α-cyclodextrin molecule are observed than in the crystal but the population (existence as percent of the simulation period) of each three-center hydrogen bond is lower in solution than in crystal. These may indeed serve as intermediate states in the process of changing one hydrogen bonding pattern into another.     

Solution structure of the glycosylphosphatidylinositol membrane anchor glycan of Trypanosoma brucei variant surface glycoprotein

The average solution conformation of the glycosylphosphatidylinositol (GPI) membrane anchor of Trypanosoma brucei variant surface glycoprotein (VSG) has been determined by using a combination of two-dimensional 1H-1H NMR methods together with molecular orbital calculations and restrained molecular dynamics simulations. This allows the generation of a model to describe the orientation of the glycan with respect to the membrane. This shows that the glycan exists in an extended configuration along the plane of the membrane and spans an area of 600 A2, which is similar to the cross-sectional area of a monomeric N-terminal VSG domain. Taken together, these observations suggest a possible space-filling role for the GPI anchor that may maintain the integrity of the VSG coat. The potential importance of the GPI glycan as a chemotherapeutic target is discussed in light of these observations.     

All-atom molecular dynamics simulations using orientational constraints from anisotropic NMR samples

Orientational constraints obtained from solid state NMR experiments on anisotropic samples are used here in molecular dynamics (MD) simulations for determining the structure and dynamics of several different membrane-bound molecules. The new MD technique is based on the inclusion of orientation dependent pseudo-forces in the COSMOS-NMR force field. These forces drive molecular rotations and re-orientations in the simulation, such that the motional time-averages of the tensorial NMR properties approach the experimentally measured parameters. The orientational-constraint-driven MD simulations are universally applicable to all NMR interaction tensors, such as chemical shifts, dipolar couplings and quadrupolar interactions. The strategy does not depend on the initial choice of coordinates, and is in principle suitable for any flexible molecule. To test the method on three systems of increasing complexity, we used as constraints some deuterium quadrupolar couplings from the literature on pyrene, cholesterol and an antimicrobial peptide embedded in oriented lipid bilayers. The MD simulations were able to reproduce the NMR parameters within experimental error. The alignment of the three membrane-bound molecules and some aspects of their conformation were thus derived from the NMR data, in good agreement with previous analyses. Furthermore, the new approach yielded for the first time the distribution of segmental orientations with respect to the membrane and the order parameter tensors of all three systems.     

Conformation and dynamics of the three-helix bundle UBA domain of p62 from experiment and simulation

The ubiquitin associated domain of p62 is a small three-helix bundle of approximately 50 residues that mediates the recognition of polyubiquitin chains and ubiquitylated substrates. The solution structure of a 52 residue construct containing this domain has been characterized using heteronuclear nuclear magnetic resonance (NMR) methods. The resulting ensemble of NMR-derived structures was used in molecular dynamics (MD) simulations to investigate the equilibrium conformation and dynamics of this domain. NOE and (15)N relaxation data have been used to validate the structural ensemble produced by the MD simulations and show a good correlation for residues in regions of secondary structure. A similar approach was taken using an ensemble of structures from the MD simulations to calculate electronic circular dichroism (CD) and IR spectra from first principles with an encouraging correlation with the experimental CD and IR data.     

Conformation of a trimannoside bound to mannose-binding protein by nuclear magnetic resonance and molecular dynamics simulations

A model of the carbohydrate recognition domain of the serum form of mannose-binding protein (MBP) from rat complexed with methyl 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside is presented. Allowed conformations for the bound sugar were derived from simulated annealing protocols incorporating distance restraints computed from transferred NOESY spectra. The resulting sugar conformations were then modeled into the MBP binding site, and these models of the complex were refined using molecular dynamics (MD) simulations in the presence of solvent water. These studies indicate that only one of the two major conformations of the alpha(1-->6) linkage found in solution is significantly populated in the bound state (omega = 60 degrees ), whereas the alpha(1-->3) linkage samples at least two states, similar to its behavior in free solution. The bound conformation allows direct hydrogen bonds to form between the sugar and K182 of MBP, in addition to other water-mediated hydrogen bonds. Estimates of binding constants of candidate complexes based on changes in solvent-accessible surface areas upon binding support the NMR and MD results. These estimates further suggest that the enthalpic gains of the additional sugar-MBP interactions in a trisaccharide as opposed to a monosaccharide are offset by entropic penalties, offering an explanation for previous binding data.     

Effect of conformation on the conversion of cyclo-(1,7)-Gly-Arg-Gly-Asp-Ser-Pro-Asp-Gly-OH to its cyclic imide degradation product.

The objective of this study was to explain the increased propensity for the conversion of cyclo-(1,7)-Gly-Arg-Gly-Asp-Ser-Pro-Asp-Gly-OH (1), a vitronectin-selective inhibitor, to its cyclic imide counterpart cyclo-(1,7)-Gly-Arg-Gly-Asu-Ser-Pro-Asp-Gly-OH (2). Therefore, we present the conformational analysis of peptides 1 and 2 by NMR and molecular dynamic simulations (MD). Several different NMR experiments, including COSY, COSY-Relay, HOHAHA, NOESY, ROESY, DQF-COSY and HMQC, were used to: (a) identify each proton in the peptides; (b) determine the sequential assignments; (c) determine the cis-trans isomerization of X-Pro peptide bond; and (d) measure the NH-HCalpha coupling constants. NOE- or ROE-constraints were used in the MD simulations and energy minimizations to determine the preferred conformations of cyclic peptides 1 and 2. Both cyclic peptides 1 and 2 have a stable solution conformation; MD simulations suggest that cyclic peptide 1 has a distorted type I beta-turn at Arg2-Gly3-Asp4-Ser5 and cyclic peptide 2 has a pseudo-type I beta-turn at Ser5-Pro6-Asp7-Gly1. A shift in position of the type I beta-turn at Arg2-Gly3-Asp4-Ser5 in peptide 1 to Ser5-Pro6-Asp7-Gly1 in peptide 2 occurs upon formation of the cyclic imide at the Asp4 residue. Although the secondary structure of cyclic peptide 1 is not conducive to succinimide formation, the reaction proceeds via neighbouring group catalysis by the Ser5 side chain. This mechanism is also supported by the intramolecular hydrogen bond network between the hydroxyl side chain and the backbone nitrogen of Ser5. Based on these results, the stability of Asp-containing peptides cannot be predicted by conformational analysis alone; the influence of anchimeric assistance by surrounding residues must also be considered.     

NMR structure of duplex DNA containing the α‐OH‐PdG·dA base pair: A mutagenic intermediate of acrolein

Acrolein, a cell metabolic product and main component of cigarette smoke, reacts with DNA generating α‐OH‐PdG lesions, which have the ability to pair with dATP during replication thereby causing G to T transversions. We describe the solution structure of an 11‐mer DNA duplex containing the mutagenic α‐OH‐PdG·dA base pair intermediate, as determined by solution nuclear magnetic resonance (NMR) spectroscopy and retrained molecular dynamics (MD) simulations. The NMR data support a mostly regular right‐handed helix that is only perturbed at its center by the presence of the lesion. Undamaged residues of the duplex are in anti orientation, forming standard Watson‐Crick base pairs alignments. Duplication of proton signals at and near the damaged base pair reveals the presence of two enantiomeric duplexes, thus establishing the exocyclic nature of the lesion. The α‐OH‐PdG adduct assumes a syn conformation pairing to its partner dA base that is protonated at pH 6.6. The three‐dimensional structure obtained by restrained molecular dynamics simulations show hydrogen bond interactions that stabilize α‐OH‐PdG in a syn conformation and across the lesion containing base pair. We discuss the implications of the structures for the mutagenic bypass of acrolein lesions. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 391–401, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Right-handed 14-helix in beta 3-peptides from L-aspartic acid monomers

beta-Peptides made from L-aspartic acid monomers form a new class of beta 3-peptides. Here we report the first three-dimensional NMR solution structure of a beta 3-hexapeptide (1) from L-aspartic acid monomers in 2,2,2-trifluoroethanol (TFE). We show that 1 forms a right-handed 14-helical structure in TFE. alpha-peptides from naturally occurring L-amino acids adopt a right-handed alpha-helix whereas beta 3-peptides formed from beta 3-amino acids derived from naturally occurring L-amino acids form left-handed 14-helices. The right-handed 14-helical conformation of 1 is a better mimic of alpha-peptide conformations. Using the NMR structure of 1 in TFE, we further study the conformation of 1 in water, as well as two similar beta 3-peptides (2 and 3) in water and TFE by molecular dynamics (MD) simulations. NMR and MD results suggest loss of secondary structure of 1 in water and show that it forms a fully extended structure. 2 and 3 contain residues with oppositely charged side chains that engage in salt-bridge interactions and dramatically stabilize the 14-helical conformation in aqueous media.     

Exploring the energy landscape of protein folding using replica-exchange and conventional molecular dynamics simulations

Two independent replica-exchange molecular dynamics (REMD) simulations with an explicit water model were performed of the Trp-cage mini-protein. In the first REMD simulation, the replicas started from the native conformation, while in the second they started from a nonnative conformation. Initially, the first simulation yielded results qualitatively similar to those of two previously published REMD simulations: the protein appeared to be over-stabilized, with the predicted melting temperature 50-150K higher than the experimental value of 315K. However, as the first REMD simulation progressed, the protein unfolded at all temperatures. In our second REMD simulation, which starts from a nonnative conformation, there was no evidence of significant folding. Transitions from the unfolded to the folded state did not occur on the timescale of these simulations, despite the expected improvement in sampling of REMD over conventional molecular dynamics (MD) simulations. The combined 1.42 micros of simulation time was insufficient for REMD simulations with different starting structures to converge. Conventional MD simulations at a range of temperatures were also performed. In contrast to REMD, the conventional MD simulations provide an estimate of Tm in good agreement with experiment. Furthermore, the conventional MD is a fraction of the cost of REMD and continuous, realistic pathways of the unfolding process at atomic resolution are obtained.