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.     

Molecular dynamics study of methane in water: diffusion and structure

In this study, theoretical analysis on the geometries and electronic properties of various conjugated oligomers based on thiophene (Th) or bicyclic non-classical Th units is reported. The dihedral angle, bond length, bond-length alternation, bond critical point (BCP) properties, nucleus-independent chemical shift (NICS) and Wiberg bond index (WBIs) are analysed and correlated with conduction properties. The changes of bond length, BCP properties, NICS and WBIs all show that the conjugational degree is increased systematically with main chain extension. As a result, the highest occupied molecular orbital–lowest unoccupied molecular orbital energy separation (E _g) is decreased upon chain elongation. The E _g of oligomers based on bicyclic non-classical Th unit is much lower than that of Th-based oligomers due to the narrower E _g of bicyclic non-classical Ths compared to Th, which indicates that the narrow E _g of the bicyclic non-classical Ths can be carried over to their polymers by using them as building blocks for the polymers. The band structures and density of states analysis show that the four polymers all have small band gap ( < 0.9 eV), wide highest occupied bandwidth and relatively small effective mass of hole, which indicate that those proposed polymers may be potential conductors.     

Condensed State Molecular Dynamics in Sorbitol and Maltitol: Mobility Gradients and Conformation Transitions

Abstract

Molecular mobility in sorbitol and maltitol is studied in order to understand their differences near the junction between the α and β relaxations. The molecular dynamics simulations performed on the polyols in their bulk state give support to the ¹³C NMR results and imply that the mobility of a carbon atom located at the extremity of the chain is higher than that of any other carbon. Moreover, the difference in carbon atoms mobility is greater within the sorbitol moiety of maltitol than in sorbitol and seems intimately related to the junction temperature of the α and β relaxation processes. The reorientation of the C–H vectors as probed by NMR is shown to be mainly the effect of conformation transitions in the case of a carbon atom located at the end of the chain.     

Conformation, dynamics, solvation and relative stabilities of selected beta-hexopyranoses in water: a molecular dynamics study with the GROMOS 45A4 force field

The present article reports long timescale (200 ns) simulations of four beta-D-hexopyranoses (beta-D-glucose, beta-D-mannose, beta-D-galactose and beta-D-talose) using explicit-solvent (water) molecular dynamics and vacuum stochastic dynamics simulations together with the GROMOS 45A4 force field. Free-energy and solvation free-energy differences between the four compounds are also calculated using thermodynamic integration. Along with previous experimental findings, the present results suggest that the formation of intramolecular hydrogen-bonds in water is an 'opportunistic' consequence of the close proximity of hydrogen-bonding groups, rather than a major conformational driving force promoting this proximity. In particular, the conformational preferences of the hydroxymethyl group in aqueous environment appear to be dominated by 1,3-syn-diaxial repulsion, with gauche and solvation effects being secondary, and intramolecular hydrogen-bonding essentially negligible. The rotational dynamics of the exocyclic hydroxyl groups, which cannot be probed experimentally, is found to be rapid (10-100 ps timescale) and correlated (flip-flop hydrogen-bonds interconverting preferentially through an asynchronous disrotatory pathway). Structured solvent environments are observed between the ring and lactol oxygen atoms, as well as between the 4-OH and hydroxymethyl groups. The calculated stability differences between the four compounds are dominated by intramolecular effects, while the corresponding differences in solvation free energies are small. An inversion of the stereochemistry at either C(2) or C(4) from equatorial to axial is associated with a raise in free energy. Finally, the particularly low hydrophilicity of beta-D-talose appears to be caused by the formation of a high-occurrence hydrogen-bonded bridge between the 1,3-syn-diaxial 2-OH and 4-OH groups. Overall, good agreement is found with available experimental and theoretical data on the structural, dynamical, solvation and energetic properties of these compounds. However, this detailed comparison also reveals some discrepancies, suggesting the need (and providing a solid basis) for further refinement.     

Molecular dynamics simulation studies of mechanical properties of different carbon nanotube systems

Various mechanical properties of single-walled carbon nanotubes (SWCNT) and double-walled carbon nanotubes (DWCNT) are evaluated using molecular dynamics (MD) simulations. A tensioning process was first performed on a SWCNT whose interaction is based on the Brenner’s ‘second generation’ potential under varying length–diameter ratios and strain rates, in order to understand the SWCNT’s behaviour under axial tension. The results showed an increase in the SWCNT’s ultimate tensile strength and a decrease in critical strain given the conditions of increasing strain rate and a decreasing length–diameter ratio. Comparison was done with previous studies on axial tensioning of SWCNT to validate the results obtained from the set-up, based on the general stress–strain relationship and key mechanical properties such as the strain at failure and the Young’s modulus. A DWCNT was then constructed, and Lennard-Jones ‘12-6’ potential was used to describe the energy present between the nanotube layers. Extraction of the inner tube in a DWCNT was performed using two inner wall tubings of different diameters to draw comparison to the energies needed to separate fully the outer and inner tubing. Finally, a bending test was performed on two DWCNTs with different intertube separations. Insights into the entire bending process were obtained through analyses of the variations in the strain energy characteristic of the surface atoms near the bending site, as the DWCNT is gradually bent until failure.     

Molecular dynamics simulation of liquid water under the influence of an external electric field

Molecular dynamics simulations of liquid water were performed at 258K and a density of 1.0?g/cm³ under various applied external electric field, ranging 0～10¹⁰?V/m. The influence of external field on structural and dynamical properties of water was investigated. The simple point charge (SPC) model is used for water molecules. An enhancement of the water hydrogen bond structure with increasing strength of the electric field has been deduced from the radial distribution functions and the analysis of hydrogen bonds structure. With increasing field strength, water system has a more perfect structure, which is similar to ice structure. However, the electrofreezing phenomenon of liquid water has not been detected since the self-diffusion coefficient was very large. The self-diffusion coefficient decreases remarkably with increasing strength of electric field and the self-diffusion coefficient is anisotropic.     

Molecular dynamics simulation of silicate glasses

Summary Although the structure of glasses is not really accessible by experimental methods, molecular dynamics is a very useful alternative, as we have tried to demonstrate in this chapter. The simulations reproduce the broad macroscopic features found in these glasses, both structural and transport-related, providing a basis for the more detailed atomic scale features found in the simulated structures. An understanding of important aspects of alkali ion transport, such as the mixed alkali effect and anomalous behaviour in some alumino-silicates, can thus be approached from the atomistic pictures of the glasses produced by the simulations. Although there is room for improvements to the potential models available, it should be clear that the further application of computer simulation methods, such as molecular dynamics, promises to provide much needed advances in glass science and engineering.     

Molecular dynamics simulations of polyampholytes inside a slit

Alternating and diblock polyampholytes confined in a slit with and without an electric field have been simulated by the molecular dynamics method with a Langevin thermostat. It is shown that the slit has a strong effect on the properties of the polyampholyte. The effect is stronger when the electric field is weak, or the temperature is not too high. When a polyampholyte chain moves close to the slit wall, its radius of gyration perpendicular to the wall becomes smaller but that parallel to the wall becomes larger. Owing to the confinement of the slit, the polyampholyte chain closer to the slit wall tends to lie on the wall and becomes more flat. The width of the slit has only a little influence on the properties of solutions near the slit wall, values of several physical statistics are very close with different widths. However, when the electric field strength is strong enough in a narrow slit, the obtained properties obviously differ.     

Molecular dynamics simulations on hydrogen adsorption in finite single walled carbon nanotube bundles

Molecular dynamics simulations of the adsorption of hydrogen molecules in finite single-walled carbon nanotube bundles are presented using a curvature dependent force field. The heat of formation and the effective adsorption capacity are expressed as a function of H₂ distance from adsorbent. The heat of adsorption decreases rapidly with the distance and increasing H₂ loading results in weakening adsorption strength. The effects of nanotube packing and bundle thickness on hydrogen adsorption strength were investigated and the results show that the heat of adsorption can be improved slightly if hydrogen molecules are placed in thicker and inhomogeneously packed nanotube bundles. Only very small diameter nanotube bundles were found to hold promise for significant hydrogen storage for onboard applications.     

1H NMR studies on the hydrogen-bonding network in mono-altro-beta-cyclodextrin and its complex with adamantane-1-carboxylic acid

The hydrogen-bond network in mono-altro-beta-cyclodextrin and in its inclusion complex with adamantane-1-carboxylic acid were investigated by (1)H NMR spectroscopy using the chemical shifts, temperature coefficients and vicinal coupling constants of the exchangeable hydroxy protons. The chemical shifts of the 3-OH signals indicated that the hydrogen-bond network between the 2-OH and 3-OH groups was disturbed not only on each side of the altrose residue, but also along the whole dextrin chain. Upon addition of adamantane-1-carboxylic acid, altrose underwent a conformational change from the (1)C(4) to the (O)S(2) form, allowing a more continuous belt of hydrogen bonding, as evidenced by the downfield shift experienced by the 3-OH proton signals of the glucose residues.     

Structural stability of myoglobin and glycomyoglobin: a comparative molecular dynamics simulation study

Glycoproteins are formed as the result of enzymatic glycosylation or chemical glycation in the body, and produced in vitro in industrial processes. The covalently attached carbohydrate molecule(s) confer new properties to the protein, including modified stability. In the present study, the structural stability of a glycoprotein form of myoglobin, bearing a glucose unit in the N-terminus, has been compared with its native form by the use of molecular dynamics simulation. Both structures were subjected to temperatures of 300 and 500 K in an aqueous environment for 10 ns. Changes in secondary structures and RMSD were then assessed. An overall higher stability was detected for glycomyoglobin, for which the most stable segments/residues were highlighted and compared with the native form. The simple addition of a covalently bound glucose is suggested to exert its stabilizing effect via increased contacts with surrounding water molecules, as well as a different pattern of interactions with neighbor residues.

Electronic supplementary material

The online version of this article (doi:10.1007/s10867-015-9383-2) contains supplementary material, which is available to authorized users.     

Molecular aspects of the interaction between amphotericin B and a phospholipid bilayer: molecular dynamics studies

Amphotericin B (AmB) is a polyene macrolide antibiotic used to treat systemic fungal infections. The molecular mechanism of AmB action is still only partly characterized. AmB interacts with cell-membrane components and forms membrane channels that eventually lead to cell death. The interaction between AmB and the membrane surface can be regarded as the first (presumably crucial) step on the way to channel formation. In this study molecular dynamics simulations were performed for an AmB–lipid bilayer model in order to characterize the molecular aspects of AmB–membrane interactions. The system studied contained a box of 200 dimyristoylphosphatidylcholine (DMPC) molecules, a single AmB molecule placed on the surface of the lipid bilayer and 8,065 water molecules. Two molecular dynamics simulations (NVT ensemble), each lasting 1 ns, were performed for the model studied. Two different programs, CHARMM and NAMD2, were used in order to test simulation conditions. The analysis of MD trajectories brought interesting information concerning interactions between polar groups of AmB and both DMPC and water molecules. Our studies show that AmB preferentially took a vertical position, perpendicular to the membrane surface, with no propensity to enter the membrane. Our finding may suggest that a single AmB molecule entering the membrane is very unlikely.Figure The figure presents the whole structure of the system simulated—starting point. AmB is presented as a space-filling model, DMPC molecules—green sticks, water molecules—red sticks     

Molecular dynamics studies of the transmembrane domain of gp41 from HIV-1

Helix-helix interactions in the putative three-helix bundle formation of the gp41 transmembrane (TM) domain may contribute to the process of virus-cell membrane fusion in HIV-1 infection. In this study, molecular dynamics is used to analyze and compare the conformations of monomeric and trimeric forms of the TM domain in various solvent systems over the course of 4 to 23-ns simulations. The trimeric bundles of the TM domain were stable as helices and remained associated in a hydrated POPE lipid bilayer for the duration of the 23-ns simulation. Several stable inter-chain hydrogen bonds, mostly among the three deprotonated arginine residues located at the center of each of the three TM domains, formed in a right-handed bundle embedded in the lipid bilayer. No such bonds were observed when the bundle was left-handed or when the central arginine residue in each of the three TM helices was replaced with isoleucine (R_I mutant), suggesting that the central arginine residues may play an essential role in maintaining the integrity of the three-helix bundle. These observations suggest that formation of the three-helix bundle of the TM domain may play a role in the trimerization of gp41, thought to occur during the virus-cell membrane fusion process.     

Comparative studies on the dynamics of crosslinked insulin

Molecular dynamics simulations were carried out on an insulin crosslinked between the N-terminal A chain and the C-terminal B chain to form a so-called mini-proinsulin: N -^A1-N -^B29-diaminosuberoyl insulin (DASI). To investigate the influence of crosslinking on the dynamics of the insulin moiety, the bridge was removed from a transient DASI structure and simulation was carried on independently with the then unlinked (ULKI) as well as with the crosslinked species. The effects of crystal packing and quaternary interactions were checked by simulating both types of monomers and dimers known from the hexamer structure. All simulations were compared to previous ones of native insulin. DASI shows general similarity to the native simulations in most parts of the structure. Deviations are visible in the segments to which the bridge is directly connected, i.e. their flexibility is reduced. Upon removal of the bridge the ULKI simulations reapproach those of native insulin. The influence of the bridge spreads over the whole molecule, but all of its main structural features remain intact. The simulations suggest that the displacement of the C-terminal B chain of native insulin, considered important for receptor interaction, is prevented by the bridge, which also partially shields some binding residues. This is in accordance with the poor biological potency of A1-B29-crosslinked insulins.Abbreviations DASI-insulin(DASI) bovineN -^A1-N -^B29-di-aminosuberoyl insulin - ULK-insulin (ULKI) Native beef insulin with the bridge of DASI removed     

Investigations of peptide hydration using NMR and molecular dynamics simulations: A study of effects of water on the conformation and dynamics of antamanide

Summary The influence of water binding on the conformational dynamics of the cyclic decapeptide antamanide dissolved in the model lipophilic environment chloroform is investigated by NMR relaxation measurements. The water-peptide complex has a lifetime of 35 s at 250 K, which is longer than typical lifetimes of water-peptide complexes reported in aqueous solution. In addition, there is a rapid intracomplex mobility that probably involves librational motions of the bound water or water molecules hopping between different binding sites. Water binding restricts the flexibility of antamanide. The experimental findings are compared with GROMOS molecular dynamics simulations of antamanide with up to eight bound water molecules. Within the simulation time of 600 ps, no water molecule leaves the complex. Additionally, the simulations show a reduced flexibility for the complex in comparison with uncomplexed antamanide. Thus, there is a qualitative agreement between the experimental NMR results and the computer simulations.     

Conformational properties of glucose-based disaccharides investigated using molecular dynamics simulations with local elevation umbrella sampling

Explicit-solvent molecular dynamics (MD) simulations of the 11 glucose-based disaccharides in water at 300 K and 1 bar are reported. The simulations were carried out with the GROMOS 45A4 force-field and the sampling along the glycosidic dihedral angles ? and ψ was artificially enhanced using the local elevation umbrella sampling (LEUS) method. The trajectories are analyzed in terms of free-energy maps, stable and metastable conformational states (relative free energies and estimated transition timescales), intramolecular H-bonds, single molecule configurational entropies, and agreement with experimental data. All disaccharides considered are found to be characterized either by a single stable (overwhelmingly populated) state ((1→n)-linked disaccharides with n = 1, 2, 3, or 4) or by two stable (comparably populated and differing in the third glycosidic dihedral angle ; gg or gt) states with a low interconversion barrier ((1→6)-linked disaccharides). Metastable (anti-? or anti-ψ) states are also identified with relative free energies in the range of 8-22 kJ mol⁻¹. The 11 compounds can be classified into four families: (i) the α(1→1)α-linked disaccharide trehalose (axial-axial linkage) presents no metastable state, the lowest configurational entropy, and no intramolecular H-bonds; (ii) the four α(1→n)-linked disaccharides (n = 1, 2, 3, or 4; axial-equatorial linkage) present one metastable (anti-ψ) state, an intermediate configurational entropy, and two alternative intramolecular H-bonds; (iii) the four β(1→n)-linked disaccharides (n = 1, 2, 3, or 4; equatorial-equatorial linkage) present two metastable (anti-? and anti-ψ) states, an intermediate configurational entropy, and one intramolecular H-bond; (iv) the two (1→6)-linked disaccharides (additional glycosidic dihedral angle) present no (isomaltose) or a pair of (gentiobiose) metastable (anti-?) states, the highest configurational entropy, and no intramolecular H-bonds. The observed conformational preferences appear to be dictated by four main driving forces (ring conformational preferences, exo-anomeric effect, steric constraints, and possible presence of a third glycosidic dihedral angle), leaving a secondary role to intramolecular H-bonding and specific solvation effects. In spite of the weak conformational driving force attributed to solvent-exposed H-bonds in water (highly polar protic solvent), intramolecular H-bonds may still have a significant influence on the physico-chemical properties of the disaccharide by decreasing its hydrophilicity. Along with previous work, the results also complete the suggestion of a spectrum of approximate transition timescales for carbohydrates up to the disaccharide level, namely: ∼30 ps (hydroxyl groups), ∼1 ns (free lactol group, free hydroxymethyl groups, glycosidic dihedral angle in (1→6)-linked disaccharides), ∼10 ns to 2 μs (ring conformation, glycosidic dihedral angles ? and ψ). The calculated average values of the glycosidic torsional angles agree well with the available experimental data, providing validation for the force-field and simulation methodology employed.     

分子动力学模拟与分子生物力学

生物大分子的微观结构动力学决定其生物学功能,其力学-化学耦合规律是分子生物力学的重点关注方向。分子动力学模拟是耦合生物大分子力学-化学性质微观结构动力学基础的有效手段,其结果可用于预测结构-功能关系、指导实验设计和诠释实验结果。本文简要介绍了分子动力学模拟的方法学特点、基本工作原理及其在分子生物力学中的应用,并展望了未来可能的发展方向和应用前景。     

Molecular modeling and dynamics studies of purine nucleoside phosphorylase from Bacteroides fragilis

Purine nucleoside phosphorylase (PNP) catalyzes the reversible phosphorolysis of N-ribosidic bonds of purine nucleosides and deoxynucleosides, except adenosine, to generate ribose 1-phosphate and the purine base. This work describes for the first time a structural model of PNP from Bacteroides fragilis (Bf). We modeled the complexes of BfPNP with six different ligands in order to determine the structural basis for specificity of these ligands against BfPNP. Comparative analysis of the model of BfPNP and the structure of HsPNP allowed identification of structural features responsible for differences in the computationally determined ligand affinities. The molecular dynamics (MD) simulation was assessed to evaluate the overall stability of the BfPNP model. The superposition of the final onto the initial minimized structure shows that there are no major conformational changes from the initial model, which is consistent with the relatively low root mean square deviation (RMSD). The results indicate that the structure of the model was stable during MD, and does not exhibit loosely structured loop regions or domain terminals.     

Molecular dynamics, crystallography and mutagenesis studies on the substrate gating mechanism of prolyl oligopeptidase

Altered prolyl oligopeptidase (PREP) activity is found in many common neurological and other genetic disorders, and in some cases PREP inhibition may be a promising treatment. The active site of PREP resides in an internal cavity; in addition to the direct interaction between active site and substrate or inhibitor, the pathway to reach the active site (the gating mechanism) must be understood for more rational inhibitor design and understanding PREP function. The gating mechanism of PREP has been investigated through molecular dynamics (MD) simulation combined with crystallographic and mutagenesis studies. The MD results indicate the inter-domain loop structure, comprised of 3 loops at residues, 189-209 (loop A), 577-608 (loop B), and 636-646 (loop C) (porcine PREP numbering), are important components of the gating mechanism. The results from enzyme kinetics of PREP variants also support this hypothesis: When loop A is (1) locked to loop B through a disulphide bridge, all enzyme activity is halted, (2) nicked, enzyme activity is increased, and (3) removed, enzyme activity is only reduced. Limited proteolysis study also supports the hypothesis of a loop A driven gating mechanism. The MD results show a stable network of H-bonds that hold the two protein domains together. Crystallographic study indicates that a set of known PREP inhibitors inhabit a common binding conformation, and this H-bond network is not significantly altered. Thus the domain separation, seen to occur in lower taxa, is not involved in the gating mechanism for mammalian PREP. In two of the MD simulations we observed a conformational change that involved the breaking of the H-bond network holding loops A and B together. We also found that this network was more stable when the active site was occupied, thus decreasing the likelihood of this transition.     

Molecular dynamics study of deformation mechanisms of poly(vinyl alcohol) hydrogel

Molecular dynamics (MD) simulations have been performed on the physically crosslinking poly(vinyl alcohol) (PVA) hydrogel to study the deformation mechanisms under uniaxial tensile conditions. The distributions of hydroxyl oxygens and dihedral angle and the number of hydrogen bonds have been analysed to study the structure of the hydrogel. The water content and temperature dependency of mechanical properties have been investigated. The energy contributions from the partially united atom potential have been calculated as a function of strain. It is found that the stress–strain curve comprises toe region, linear region and yield and failure region which is close to most biomaterials. In the toe and yield region, all the contributions to the internal energy change a little. However, in the linear region, the bond stretching and angle bending energy increase rapidly and mainly dominate the region, and the energy increases more rapidly with the increasing water content but the decreasing temperature. The degree of crosslinking decreases with the increasing deformation. The polymer chains occur significant torsional activity in the toe region. Hydrogen bonds are stable in the toe and yield region, but the hydrogen bonds between hydroxyl groups and waters decrease in the linear region.