Studies on the conformational behaviour of GlcNAc-Man3-GlcNAc2 oligosaccharides using molecular dynamics simulations

Three-dimensional structures of the natural substrate unit for the enzyme N-acetylglucosamine-transferase II, GlcNAc-Man3-GlcNAc2, were investigated by molecular modelling methods. Molecular dynamics (MD) and molecular mechanics calculations on two hexasaccharides, namely GlcNAc-Man3-GlcNAc2-Asn and GlcNAc-Man3-GlcNAc2-OMe were performed by the Biosym/MSI software using the CVFF and CFF95 force fields in vacuum. The MD simulations were calculated for 3 ns at different simulation temperatures and for two values of dielectric constant, = 1 and = 4. From each 3 ns trajectory, 3050 structures have been optimized. The local minima obtained have been clustered into families exhibiting similar values of glycosidic torsional angles phi, psi, and omega. The influence of the simulation conditions and force fields used on the conformational behaviour and structure of the title oligosaccharides is discussed.     

A systematic molecular dynamics simulation study of temperature dependent bilayer structural properties

Although lipid force fields (FFs) used in molecular dynamics (MD) simulations have proved to be accurate, there has not been a systematic study on their accuracy over a range of temperatures. Motivated by the X-ray and neutron scattering measurements of common phosphatidylcholine (PC) bilayers (Ku?erka et al. BBA. 1808: 2761, 2011), the CHARMM36 (C36) FF accuracy is tested in this work with MD simulations of six common PC lipid bilayers over a wide range of temperatures. The calculated scattering form factors and deuterium order parameters from the C36 MD simulations agree well with the X-ray, neutron, and NMR experimental data. There is excellent agreement between MD simulations and experimental estimates for the surface area per lipid, bilayer thickness (D_B), hydrophobic thickness (D_C), and lipid volume (V_L). The only minor discrepancy between simulation and experiment is a measure of (D_B − D_HH) / 2 where D_HH is the distance between the maxima in the electron density profile along the bilayer normal. Additional MD simulations with pure water and heptane over a range of temperatures provide explanations of possible reasons causing the minor deviation. Overall, the C36 FF is accurate for use with liquid crystalline PC bilayers of varying chain types and over biologically relevant temperatures.     

Anisotropic atomic motions in high-resolution protein crystallography molecular dynamics simulations

Molecular dynamics (MD) simulations using empirical force fields are popular for the study of proteins. In this work, we compare anisotropic atomic fluctuations in nanosecond-timescale MD simulations with those observed in an ultra-high-resolution crystal structure of crambin. In order to make our comparisons, we have developed a compact graphical technique for assessing agreement between spatial atomic distributions determined by MD simulations and observed anisotropic temperature factors.     

MD simulation of subtilisin BPN' in a crystal environment.

In this paper we present a molecular dynamics (MD) simulation of subtilisin BPN' in a crystalline environment containing four protein molecules and solvent. Conformational and dynamic properties of the molecules are compared with each other and with respect to the X-ray structure to test the validity of the force field. The agreement between simulated and experimental structure using the GROMOS force field is better than that obtained in the literature using other force fields for protein crystals. The overall shape of the molecule is well preserved, as is the conformation of alpha-helices and beta-strands. Structural differences are mainly found in loop regions. Solvent networks found in the X-ray structure were reproduced by the simulation, which was unbiased with respect to the crystalline hydration structure. These networks seem to play an important role in the stability of the protein; evidence of this is found in the structure of the active site. The weak ion binding site in the X-ray structure of subtilisin BPN' is occupied by a monovalent ion. When a calcium ion is placed in the initial structure, three peptide ligands are replaced by 5 water ligands, whereas a potassium ion retains (in part) its original ligands. Existing force fields yield a reliable method to probe local structure and short-time dynamics of proteins, providing an accuracy of about 0.1 nm.     

Molecular dynamics simulation and steered molecular dynamics simulation on irisin dimers

Irisin is found closely associated with promoting the browning of beige fat cells in white adipose tissue. The crystal structure reveals that irisin forms a continuous inter-subunit β-sheet dimer. Here, molecular dynamics (MD) simulation and steered molecular dynamics (SMD) simulation were performed to investigate the dissociation process and the intricate interactions between the two irisin monomers. In the process of MD, the interactions between the monomers were roughly analyzed through the average numbers of both hydrophobic contacts and H-bonds. Then, SMD was performed to investigate the accurate interaction energy between the monomers. By the analysis of dissociation energy, the van der Waals (vdW) force was identified as the major energy to maintain the dimer structure, which also verified the results of MD simulation. Meanwhile, 11 essential residues were discovered by the magnitude of rupture force during dissociation. Among them, residues Arg75, Glu79, Ile77, Ala88, and Trp90 were reported in a previous study using the method of mutagenesis and size exclusion chromatography, and several new important residues (Arg72, Leu74, Phe76, Gln78, Val80, and Asp91) were also identified. Interestingly, the new important residues that we discovered and the important residues that were reported are located in the opposite side of the β-sheet of the dimer.     

Crystal Structure Determination of New Antimitotic Agent Bis(p-fluorobenzyl)trisulfide

The purpose of this research was to investigate the physical characteristics and crystalline structure of bis(p-fluorobenzyl)trisulfide, a new anti-tumor agent. Methods used included X-ray single crystal diffraction, X-ray powder diffraction (XRPD), Fourier-transform infrared (FT-IR) spectroscopy, differential scanning calorimetric (DSC) and thermogravimetric (TG) analyses. The findings obtained with X-ray single crystal diffraction showed that a monoclinic unit cell was a = 12.266(1) A, b = 4.7757(4) A, c = 25.510(1) A, beta = 104.25(1) degrees ; cell volume = 1,448.4(2) A(3), Z = 4, and space group C2/c. The XRPD studies of the four crystalline samples, obtained by recrystallization from four different solvents, indicated that they had the same diffraction patterns. The diffraction pattern stimulated from the crystal structure data is in excellent agreement with the experimental results. In addition, the identical FT-IR spectra of the four crystalline samples revealed absorption bands corresponding to S-S and C-S stretching as well as the characteristic aromatic substitution. Five percent weight loss at 163.3 degrees C was observed when TG was used to study the decomposition process in the temperature range of 20-200 degrees C. DSC also allowed for the determination of onset temperatures at 60.4(1)-60.7(3) degrees C and peak temperatures at 62.1(3)-62.4(3) degrees C for the four crystalline samples studied. The results verified that the single crystal structure shared the same crystal form with the four crystalline samples investigated.     

Molecular dynamics simulations of a beta-hairpin fragment of protein G: balance between side-chain and backbone forces

How is the native structure encoded in the amino acid sequence? For the traditional backbone centric view, the dominant forces are hydrogen bonds (backbone) and phi-psi propensity. The role of hydrophobicity is non-specific. For the side-chain centric view, the dominant force of protein folding is hydrophobicity. In order to understand the balance between backbone and side-chain forces, we have studied the contributions of three components of a beta-hairpin peptide: turn, backbone hydrogen bonding and side-chain interactions, of a 16-residue fragment of protein G. The peptide folds rapidly and cooperatively to a conformation with a defined secondary structure and a packed hydrophobic cluster of aromatic side-chains. Our strategy is to observe the structural stability of the beta-hairpin under systematic perturbations of the turn region, backbone hydrogen bonds and the hydrophobic core formed by the side-chains, respectively. In our molecular dynamics simulations, the peptides are solvated. with explicit water molecules, and an all-atom force field (CFF91) is used. Starting from the original peptide (G41EWTYDDATKTFTVTE56), we carried out the following MD simulations. (1) unfolding at 350 K; (2) forcing the distance between the C(alpha) atoms of ASP47 and LYS50 to be 8 A; (3) deleting two turn residues (Ala48 and Thr49) to form a beta-sheet complex of two short peptides, GEWTYDD and KTFTVTE; (4) four hydrophobic residues (W43, Y45, F52 and T53) are replaced by a glycine residue step-by-step; and (5) most importantly, four amide hydrogen atoms (T44, D46, T53, and T55, which are crucial for backbone hydrogen bonding), are substituted by fluorine atoms. The fluorination not only makes it impossible to form attractive hydrogen bonding between the two beta-hairpin strands, but also introduces a repulsive force between the two strands due to the negative charges on the fluorine and oxygen atoms. Throughout all simulations, we observe that backbone hydrogen bonds are very sensitive to the perturbations and are easily broken. In contrast, the hydrophobic core survives most perturbations. In the decisive test of fluorination, the fluorinated peptide remains folded under our simulation conditions (5 ns, 278 K). Hydrophobic interactions keep the peptide folded, even with a repulsive force between the beta-strands. Thus, our results strongly support a side-chain centric view for protein folding.     

Low-temperature molecular dynamics simulations of horse heart cytochrome c and comparison with inelastic neutron scattering data

Molecular dynamics (MD) simulation combined with inelastic neutron scattering can provide information about the thermal dynamics of proteins, especially the low-frequency vibrational modes responsible for large movement of some parts of protein molecules. We performed several 30-ns MD simulations of cytochrome c (Cyt c) in a water box for temperatures ranging from 110 to 300 K and compared the results with those from experimental inelastic neutron scattering. The low-frequency vibrational modes were obtained via dynamic structure factors, S(Q, ω), obtained both from inelastic neutron scattering experiments and calculated from MD simulations for Cyt c in the same range of temperatures. The well known thermal transition in structural movements of Cyt c is clearly seen in MD simulations; it is, however, confined to unstructured fragments of loops Ω₁ and Ω₂; movement of structured loop Ω₃ and both helical ends of the protein is resistant to thermal disturbance. Calculated and experimental S(Q, ω) plots are in qualitative agreement for low temperatures whereas above 200 K a boson peak vanishes from the calculated plots. This may be a result of loss of crystal structure by the protein–water system compared with the protein crystal.     

Studies of short-wavelength collective molecular motions in lipid bilayers using high resolution inelastic X-ray scattering

We summarize a series of experimental results made with the newly developed high resolution X-ray scattering (IXS) instrument on two pure lipid bilayers, including dimyristoylphosphatidylcholine (DMPC) and dilauroylphosphatidylcholine (DLPC) in both gel and liquid crystal phases, and lipid bilayers containing cholesterol. By analyzing the IXS data based on the generalized three effective eigenmode model (GTEE), we obtain dispersion relations of the high frequency density oscillations (phonons) of lipid molecules in these bilayers. We then compare the dispersion relations of pure lipid bilayers of different chain lengths among themselves and the dispersion relations of pure lipid bilayers with those of the cholesterol containing bilayers. We also compare our experimental results with collective dynamics data generated by computer molecular dynamics (MD) simulations for dipalmitoylphosphatidylcholine (DPPC) in gel phase and DMPC in liquid crystal phase.     

The effect of hydrogen bonding in two crystal modifications of aquabis(N,N-dimethylglycinato-κN,κO)copper(II): Experimental and theoretical study

From the crystals of trans aquabis(N,N-dimethylglycinato-κN,κO)copper(II) dihydrate (compound 1, space group P2₁2₁2₁) novel crystal structure of trans aquabis(N,N-dimethylglycinato-κN,κO)copper(II) (compound 2, space group Pbca) was obtained and analysed by X-ray diffraction. In the crystal structure 1, the O-H?O hydrogen bonds form three-dimensional network. In the crystal structure 2, two-dimensional layers stacking to each other are formed, with non-polar N,N-dimethyl groups placed on the opposite sides of the layers, and with the polar part in the middle forming CO?O-H and C-H?O hydrogen bonds. Different hydrogen bonding patterns in 1 and 2 do not pronouncedly affect molecular geometry of the title compound. Molecular mechanics force field suited for studying the properties of bis(amino acidato)copper(II) complexes in the solid state can follow the differences between the experimental molecular structures in the two diverse crystalline surroundings. To make possible direct comparison between crystal lattices, the force field was applied to predict unit cell packing of supposed anhydrous bis(N,N-dimethylglycinato)copper(II) in space group Pbca. Relative intermolecular energies of hypothetic anhydrous crystal and simulated 1 and 2 crystals are discussed. On the basis of experimental and theoretical results we conclude that the main effect of two water molecules of crystallisation in 1 is to stabilise the crystal packing via hydrogen bonding, whilst similar pyramidal copper(II) coordination geometry in 1 and 2 is due to axially coordinated water molecule and its intermolecular interactions.     

Ensemble MD simulations restrained via crystallographic data: Accurate structure leads to accurate dynamics

Currently, the best existing molecular dynamics (MD) force fields cannot accurately reproduce the global free‐energy minimum which realizes the experimental protein structure. As a result, long MD trajectories tend to drift away from the starting coordinates (e.g., crystallographic structures). To address this problem, we have devised a new simulation strategy aimed at protein crystals. An MD simulation of protein crystal is essentially an ensemble simulation involving multiple protein molecules in a crystal unit cell (or a block of unit cells). To ensure that average protein coordinates remain correct during the simulation, we introduced crystallography‐based restraints into the MD protocol. Because these restraints are aimed at the ensemble‐average structure, they have only minimal impact on conformational dynamics of the individual protein molecules. So long as the average structure remains reasonable, the proteins move in a native‐like fashion as dictated by the original force field. To validate this approach, we have used the data from solid‐state NMR spectroscopy, which is the orthogonal experimental technique uniquely sensitive to protein local dynamics. The new method has been tested on the well‐established model protein, ubiquitin. The ensemble‐restrained MD simulations produced lower crystallographic R factors than conventional simulations; they also led to more accurate predictions for crystallographic temperature factors, solid‐state chemical shifts, and backbone order parameters. The predictions for ¹⁵N relaxation rates are at least as accurate as those obtained from conventional simulations. Taken together, these results suggest that the presented trajectories may be among the most realistic protein MD simulations ever reported. In this context, the ensemble restraints based on high‐resolution crystallographic data can be viewed as protein‐specific empirical corrections to the standard force fields.     

Measuring gluconeogenesis using a low dose of 2H2O: advantage of isotope fractionation during gas chromatography

The contribution of gluconeogenesis to glucose production can be measured by enriching body water with (2)H(2)O to approximately 0.5% (2)H and determining the ratio of (2)H that is bound to carbon-5 vs. carbon-2 of blood glucose. This labeling ratio can be measured using gas chromatography-mass spectrometry after the corresponding glucose carbons are converted to formaldehyde and then to hexamethylenetetramine (HMT). We present a technique for integrating ion chromatograms that allows one to use only 0.05% (2)H in body water (i.e., 10 times less than the current dose). This technique takes advantage of the difference in gas chromatographic retention times of naturally labeled HMT and [(2)H]HMT. We discuss the advantage(s) of using a low dose of (2)H(2)O to quantify the contribution of gluconeogenesis.     

How the headpiece hinge angle is opened: New insights into the dynamics of integrin activation

How the integrin head transitions to the high-affinity conformation is debated. Although experiments link activation with the opening of the hinge angle between the betaA and hybrid domains in the ligand-binding headpiece, this hinge is closed in the liganded alpha(v)beta3 integrin crystal structure. We replaced the RGD peptide ligand of this structure with the 10th type III fibronectin module (FnIII10) and discovered through molecular dynamics (MD) equilibrations that when the conformational constraints of the leg domains are lifted, the betaA/hybrid hinge opens spontaneously. Together with additional equilibrations on the same nanosecond timescale in which small structural variations impeded hinge-angle opening, these simulations allowed us to identify the allosteric pathway along which ligand-induced strain propagates via elastic distortions of the alpha1 helix to the betaA/hybrid domain hinge. Finally, we show with steered MD how force accelerates hinge-angle opening along the same allosteric pathway. Together with available experimental data, these predictions provide a novel framework for understanding integrin activation.     

Microscopic structure and properties changes of cassava stillage residue pretreated by mechanical activation

This study has focused on the pretreatment of cassava stillage residue (CSR) by mechanical activation (MA) using a self-designed stirring ball mill. The changes in surface morphology, functional groups and crystalline structure of pretreated CSR were examined by using scanning electron microscopy (SEM), Fourier transform infrared (FTIR) spectroscopy and X-ray diffraction (XRD) under reasonable conditions. The results showed that MA could significantly damage the crystal structure of CSR, resulting in the variation of surface morphology, the increase of amorphous region ratio and hydrogen bond energy, and the decrease in crystallinity and crystalline size. But no new functional groups generated during milling, and the crystal type of cellulose in CSR still belonged to cellulose I after MA.     

A test of improved force field parameters for urea: molecular-dynamics simulations of urea crystals

Molecular-dynamics (MD) simulations of urea crystals of different shapes (cubic, rectangular prismatic, and sheet) have been performed using our previously published force field for urea. This force field has been validated by calculating values for the cohesive energy, sublimation temperature, and melting point from the MD data. The cohesive energies computed from simulations of cubic and rectangular prismatic urea crystals in vacuo at 300 K agreed very well with the experimental sublimation enthalpies reported at 298 K. We also found very good agreement between the melting points as observed experimentally and from simulations. Annealing the crystals just below the melting point leads to reconstruction to form crystal faces that are consistent with experimental observations. The simulations reveal a melting mechanism that involves surface (corner/edge) melting well below the melting point, and rotational disordering of the urea molecules in the corner/edge regions of the crystal, which then facilitates the translational motion of these molecules.     

Diffraction-based density restraints for membrane and membrane-peptide molecular dynamics simulations

We have recently shown that current molecular dynamics (MD) atomic force fields are not yet able to produce lipid bilayer structures that agree with experimentally-determined structures within experimental errors. Because of the many advantages offered by experimentally validated simulations, we have developed a novel restraint method for membrane MD simulations that uses experimental diffraction data. The restraints, introduced into the MD force field, act upon specified groups of atoms to restrain their mean positions and widths to values determined experimentally. The method was first tested using a simple liquid argon system, and then applied to a neat dioleoylphosphatidylcholine (DOPC) bilayer at 66% relative humidity and to the same bilayer containing the peptide melittin. Application of experiment-based restraints to the transbilayer double-bond and water distributions of neat DOPC bilayers led to distributions that agreed with the experimental values. Based upon the experimental structure, the restraints improved the simulated structure in some regions while introducing larger differences in others, as might be expected from imperfect force fields. For the DOPC-melittin system, the experimental transbilayer distribution of melittin was used as a restraint. The addition of the peptide caused perturbations of the simulated bilayer structure, but which were larger than observed experimentally. The melittin distribution of the simulation could be fit accurately to a Gaussian with parameters close to the observed ones, indicating that the restraints can be used to produce an ensemble of membrane-bound peptide conformations that are consistent with experiments. Such ensembles pave the way for understanding peptide-bilayer interactions at the atomic level.     

Calculations of the φ -ψ Conformational Contour Maps for N-Acetyl Alanine N'-Methyl Amide and of the Characteristic Ratios of Poly-L-alanine Using Various Molecular Mechanics Forcefields

Abstract

We generated φ -ψ conformational energy contour maps for the of N-acetyl alanine N'-methyl amide using the molecular mechanics forcefields AMBER, AMBER3, BI085, CFF91, CVFF, MM2, MM3, MM+, and SYBYL. With MM2, MM3, and MM+, we used a dielectric constant of ? = 1.5, the default effective value for these forcefields. With the other forcefields we used ? = 1 and 4, except with SYBYL, which, in Spartan 3.1, has no electrostatic term. All forcefields yielded the C^eq ₇ conformation as a low-energy minimum or the global minimum. Most of the forcefields also yielded a minimum-energy conformation in the C₅,α_R, and α_t. regions of the φ -ψ contour map. Fewer of the forcefields yielded a minimum in the C^ax ₇ region; however, adiabatic relaxation frequently lowered the relative energy of this region. Based on the appearance of the φ -ψ maps, the following pairs of forcefields were broadly similar (but not identical) to each other but dissimilar to the other pairs: AMBER3 and AMBER, BI085 and CHARMM, MM+ and MM2, SYBYL and ECEPP, and CFF91 and MM3. We used the data from the φ -ψ contour maps to compute the characteristic ratio of poly-L-alanine. Most of the computed values deviated significantly from the experimental value. Only the computed characteristic ratio of CFF91 without adiabatic relaxation at ? = 4 and MM3 without adiabatic relaxation at ? = 1.5 agreed with the experimental value.     

Co-crystallization between a thymine and a metal complex connected by triple hydrogen bonds

The synthesis and crystal structure of co-crystal between bis(dithiobiureto)platinum(II) ([Pt(dtb)₂]) with thymine is reported. The crystalline structure of [Pt(dtb)₂](thymine)₂ shows two-dimensional arrays created by hydrogen-bonding interactions. One mononuclear complex and two thymine molecules form a building unit connected by triple hydrogen bond employing ADA-DAD arrangements (A, hydrogen bond acceptor; D, hydrogen bond donor). The building unit is linked to adjacent units via additional hydrogen bonds to form planar sheet.