Molecular dynamics studies of side chain effect on the beta-1,3-D-glucan triple helix in aqueous solution

beta-1,3-D-glucans have been isolated from fungi as right-handed 6(1) triple helices. They are categorized by the side chains bound to the main triple helix through beta-(1-->6)-D-glycosyl linkage. Indeed, since a glucose-based side chain is water soluble, the presence and frequency of glucose-based side chains give rise to significant variation in the physical properties of the glucan family. Curdlan has no side chains and self-assembles to form an water-insoluble triple helical structure, while schizophyllan, which has a 1,6-D-glucose side chain on every third glucose unit along the main chain, is completely water soluble. A thermal fluctuation in the optical rotatory dispersion is observed for the side chain, indicating probable co-operative interaction between the side chains and water molecules. This paper documents molecular dynamics simulations in aqueous solution for three models of the beta-1,3-D-glucan series: curdlan (no side chain), schizophyllan (a beta-(1-->6)-D-glycosyl side-chain at every third position), and a hypothetical triple helix with a side chain at every sixth main-chain glucose unit. A decrease was observed in the helical pitch as the population of the side chain increased. Two types of hydrogen bonding via water molecules, the side chain/main chain and the side chain/side chain hydrogen bonding, play an important role in determination of the triple helix conformation. The formation of a one-dimensional cavity of diameter about 3.5 A was observed in the schizophyllan triple helix, while curdlan showed no such cavity. The side chain/side chain hydrogen bonding in schizophyllan and the hypothetical beta-1,3-D-glucan triple helix could cause the tilt of the main-chain glucose residues to the helix.     

Dielectric study of the cooperative order–disorder transition in aqueous solutions of schizophyllan,a triple-helical polysaccharide

Schizophyllan exists in aqueous solution as a triple helix, which is intact at room temperature. Its aqueous solution forms some ordered structure at low temperatures but undergoes a sharp transition to a disordered structure as the temperature is raised. The transition temperature T_c is about 7 and 18°C for H₂O and D₂O solutions, respectively. This transition was followed by time-domain reflectometry to investigate dynamic aspects of the transition. In addition to a major peak around 10 GHz, the dielectric dispersion curve of a 20 wt % schizophyllan in D₂O exhibited a small peak around 100 MHz below T_c and around 10 MHz above T_c. The major peak is due to bulk water, whereas the 100 MH_z peak is assigned to “bound” or “structured” water, and that around 10 MHz to side-chain glucose residues. However, unlike usual bound water reported for biopolymer solutions, this “structured” water disappears abruptly when the temperature becomes close to T_c without accompanying a conformational transition of the main chain. The above assignment is consistent with the structure of the ordered phase derived from previous static data that it consists of side-chain glucose residues along with nearby water molecules surrounding the helix core that are interacting with each other loosely through hydrogen bonds, and spreads radially only a layer of one or two water molecules but a long distance along the helix axis. © 1995 John Wiley & Sons, Inc.     

Molecular dynamics simulation of the stability of a 22-residue α-helix in water and 30% trifluoroethanol

A molecular dynamics (MD) simulation was performed on the α-helix H8-HC5, the C-terminal part of myoglobin (residue 132–153), under periodic boundary conditions in two different solutions, water and water with 30% (v/v) 2,2,2-trifluoroethanol (TFE), at 300 K to investigate the stability of the helix. In both simulations, the initial configuration was a canonical right-handed α-helix. In the course of the MD trajectory in water (200 ps), the helix clearly destabilized and began to unfold after 100 ps. In the TFE solution, two stable parts of helical regions were observed after 70 ps of a 200-ps MD simulation, supporting the notion that TFE acts as a structure-forming solvent. © 1993 John Wiley & Sons, Inc.     

The Role of Cross-Chain Ionic Interactions for the Stability of Collagen Model Peptides

The contribution of ionic interactions to the stability of the collagen triple helix was studied using molecular dynamics (MD) simulations and biophysical methods. To this end, we examined the stability of a host-guest collagen model peptide, Ac-GPOGPOGPYGXOGPOGPO-NH₂, substituting KGE, KGD, EGK, and DGK for the YGX sequence. All-atom, implicit solvent MD simulations show that the fraction of cross-chain ionic interactions formed is different, with the most pronounced in the KGE and KGD sequences, and the least in the DGK sequence. To test whether the fraction of cross-chain ionic interactions correlates with the stability, experimental measurements of thermostability were done using differential scanning calorimetry and circular dichroism spectroscopy. It was found that the melting temperature is very similar for KGE and KGD peptides, whereas the EGK peptide has lower thermostability and the DGK peptide is the least thermostable. A novel, to our knowledge, computational protocol termed temperature-scan MD was applied to estimate the relative stabilities of the peptides from MD simulations. We found an excellent correlation between transition temperatures obtained from temperature-scan MD and those measured experimentally. These results suggest the importance of cross-chain ionic interactions for the stability of collagen triple helix and confirm the utility of MD simulations in predicting interactions and stability in this system.     

The Role of Cross-Chain Ionic Interactions for the Stability of Collagen Model Peptides

The contribution of ionic interactions to the stability of the collagen triple helix was studied using molecular dynamics (MD) simulations and biophysical methods. To this end, we examined the stability of a host-guest collagen model peptide, Ac-GPOGPOGPYGXOGPOGPO-NH₂, substituting KGE, KGD, EGK, and DGK for the YGX sequence. All-atom, implicit solvent MD simulations show that the fraction of cross-chain ionic interactions formed is different, with the most pronounced in the KGE and KGD sequences, and the least in the DGK sequence. To test whether the fraction of cross-chain ionic interactions correlates with the stability, experimental measurements of thermostability were done using differential scanning calorimetry and circular dichroism spectroscopy. It was found that the melting temperature is very similar for KGE and KGD peptides, whereas the EGK peptide has lower thermostability and the DGK peptide is the least thermostable. A novel, to our knowledge, computational protocol termed temperature-scan MD was applied to estimate the relative stabilities of the peptides from MD simulations. We found an excellent correlation between transition temperatures obtained from temperature-scan MD and those measured experimentally. These results suggest the importance of cross-chain ionic interactions for the stability of collagen triple helix and confirm the utility of MD simulations in predicting interactions and stability in this system.     

Protein flexibility: multiple molecular dynamics simulations of insulin chain B

Multiple molecular dynamics simulations totaling more than 100 ns were performed on chain B of insulin in explicit solvent at 300 K and 400 K. Despite some individual variations, a comparison of the protein dynamics of each simulation showed similar trends and most structures were consistent with NMR experimental values, even at the elevated temperature. The importance of packing interactions in determining the conformational transitions of the protein was observed, sometimes resulting in conformations induced by localized hydrophobic interactions. The high temperature simulation generated a more diverse range of structures with similar elements of secondary structure and populated conformations to the simulations at room temperature. A broad sampling of the conformational space of insulin chain B illustrated a wide range of conformational states with many transitions at room temperature in addition to the conformational states observed experimentally. The T-state conformation associated with insulin activity was consistently present and a possible mechanism of behavior was suggested.     

Molecular dynamics simulation of proteins under high pressure: Structure,function and thermodynamics

BackgroundMolecular dynamics (MD) simulation is well-recognized as a powerful tool to investigate protein structure, function, and thermodynamics. MD simulation is also used to investigate high pressure effects on proteins. For conducting better MD simulation under high pressure, the main issues to be addressed are: (i) protein force fields and water models were originally developed to reproduce experimental properties obtained at ambient pressure; and (ii) the timescale to observe the pressure effect is often much longer than that of conventional MD simulations.Scope of reviewFirst, we describe recent developments in MD simulation methodologies for studying the high-pressure structure and dynamics of protein molecules. These developments include force fields for proteins and water molecules, and enhanced simulation techniques. Then, we summarize recent studies of MD simulations of proteins in water under high pressure.Major conclusionsRecent MD simulations of proteins in solution under pressure have reproduced various phenomena identified by experiments using high pressure, such as hydration, water penetration, conformational change, helix stabilization, and molecular stiffening.General significanceMD simulations demonstrate differences in the properties of proteins and water molecules between ambient and high-pressure conditions. Comparing the results obtained by MD calculations with those obtained experimentally could reveal the mechanism by which biological molecular machines work well in collaboration with water molecules.     

Explicit-solvent molecular dynamics simulations of the polysaccharide schizophyllan in water

Schizophyllan is a beta(1-->3)-D-glucan polysaccharide with beta(1-->6)-branched lateral glucose residues that presents a very stiff triple-helical structure under most experimental conditions. Despite the remarkable stability of this structure (which persists up to 120 degrees C in aqueous solution), schizophyllan undergoes a major change of state around 7 degrees C in water that has been hypothesized to result from an order-disorder transition in the lateral residues. This hypothesis is only supported by indirect experimental evidence and detailed knowledge (at the atomic level) concerning hydrogen-bonding networks, interactions with the solvent molecules, orientational freedom of the lateral residues, and orientational correlations among them is still lacking. In this study explicit-solvent molecular dynamics simulations of a schizophyllan fragment (complemented by simulations of its tetrasaccharide monomer) are performed at three different temperatures (273 K, 350 K, and 450 K) and with two different types of boundary conditions (finite nonperiodic or infinite periodic fragment) as an attempt to provide detailed structural and dynamical information about the triple-helical conformation in solution and the mechanism of the low-temperature transition. These simulations suggest that three important driving forces for the high stability of the triple helix are i), the limited conformational work involved in its formation; ii), the formation of a dense hydrogen-bonding network at its center; and iii), the formation of interchain hydrogen bonds between main-chain and lateral glucose residues. However, these simulations evidence a moderate and continuous variation of the simulated observables upon increasing the temperature, rather than a sharp transition between the two lowest temperatures (that could be associated with the state transition). Although water-mediated hydrogen-bonded association of neighboring lateral residues is observed, this interaction is not strong enough to promote the formation of an ordered state (correlated motions of the lateral residues), even at the lowest temperature considered.     

The investigation of the secondary structure propensities and free-energy landscapes of peptide ligands by replica exchange molecular dynamics simulations

The conformational states of two peptide sequences that bind to staphylococcal enterotoxin B are sampled by replica exchange molecular dynamic (REMD) simulations in explicit water. REMD simulations were treated with 52 replicas in the range of 280–501 K for both peptides. The conformational ensembles of both peptides are dominated by random coil, bend and turn structures with a small amount of helical structures for each temperature. In addition, while an insignificant presence of β-bridge structures were observed for both peptides, the β-sheet structure was observed only for peptide 3. The results obtained from simulations at 300 K are consistent with the experimental results obtained from circular dichroism spectroscopy. From the analysis of REMD results, we also calculated hydrophobic and hydrophilic solvent accessible surface areas for both peptides, and it was observed that the hydrophobic segments of the peptides tend to form bend or turn structures. Moreover, the free-energy landscapes of both peptides were obtained by principal component analysis to understand how the secondary structural properties change according to their complex space. From the free-energy analysis, we have found several minima for both peptides at decreased temperature. For these obvious minima of both peptides, it was observed that the random coil, bend and turn structures are still dominant and the helix, β-bridge or β-sheet structures can appear or disappear with respect to minima. On the other hand, when we compare the results of REMD with conventional MD simulations for these peptides, the configurations of peptide 3 might be trapped in energy minima during the conventional MD simulations. Hence, it can be said that the REMD simulations have provided a sufficiently high sampling efficiency.     

Competing interactions contributing to alpha-helical stability in aqueous solution.

The stability of a 15-residue peptide has been investigated using CD spectroscopy and molecular simulation techniques. The sequence of the peptide was designed to include key features that are known to stabilize alpha-helices, including ion pairs, helix dipole capping, peptide bond capping, and aromatic interactions. The degree of helicity has been determined experimentally by CD in three solvents (aqueous buffer, methanol, and trifluoroethanol) and at two temperatures. Simulations of the peptide in the aqueous system have been performed over 500 ps at the same two temperatures using a fully explicit solvent model. Consistent with the CD data, the degree of helicity is decreased at the higher temperature. Our analysis of the simulation results has focused on competition between different side-chain/side-chain and side-chain/main-chain interactions, which can, in principle, stabilize the helix. The unfolding in aqueous solution occurs at the amino terminus because the side-chain interactions are insufficient to stabilize both the helix dipole and the peptide hydrogen bonds. Loss of capping of the peptide backbone leads to water insertion within the first peptide hydrogen bond and hence unfolding. In contrast, the carboxy terminus of the alpha-helix is stable in both simulations because the C-terminal lysine residue stabilizes the helix dipole, but at the expense of an ion pair.     

Characterization of the denaturation of human alpha-lactalbumin in urea by molecular dynamics simulations

Molecular dynamics (MD) simulations were used to characterize the non-cooperative denaturation of the molten globule A-state of human alpha-lactalbumin by urea. A solvent of explicit urea and water molecules was used, corresponding to a urea concentration of approximately 6M. Three simulations were performed at temperatures of 293K, 360K and 400K, with lengths of 2 ns, 8 ns and 8 ns respectively. The results of the simulations were compared with experimental data from NMR studies of human alpha-lactalbumin and related peptides. During the simulations, hydrogen bonds were formed from the protein to both urea and water molecules as intra-protein hydrogen bonds were lost. Urea was shown to compete efficiently with water as both a hydrogen bond donor and acceptor. Radial distribution functions of water and urea around hydrophobic side chain atoms showed a significant increase in urea molecules in the solvation shell as the side chains became exposed during denaturation. A considerable portion of the native-like secondary structure persisted throughout the simulations. However, in the simulations at 360K and 400K, there were substantial changes in the packing of aromatic and other hydrophobic side chains in the protein, and many native contacts were lost. The results suggest that during the non-cooperative denaturation of the molten globule, secondary structure elements are stabilized by non-specific, non-native interactions.     

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.     

Spermine: an "invisible" component in the crystals of B-DNA. A grand canonical Monte Carlo and molecular dynamics simulation study

The association of spermine(4+) (Spm(4+)), Mg(2+) and monovalent (M(+)) ions with DNA in crystal form, have been studied using grand canonical Monte Carlo (GCMC) and molecular dynamics (MD) computer simulations. GCMC calculations were used to calculate the distribution of Spm(4+), Mg(2+), and M(+) between the equilibrating solvent and the DNA crystal under conditions mimicking the crystal-growing protocols reported in a number of recent X-ray diffraction studies of DNA oligomers. The GCMC simulations show that the composition of ions neutralizing the negative charge of DNA can vary in a broad range. The GCMC simulations were used to provide appropriate conditions for subsequent 6 ns constant pressure and temperature MD simulations of DNA in a typical crystalline environment consisting of three DNA double helix decamers in a periodic hexagonal cell, containing 1200 water molecules, eight Spm(4+), 32 Na(+) and four Cl(-) ions. Based on the simulation results, it seems possible to give an explanation why spermine molecules are usually not detected in X-ray studies in spite of their high concentration in the preparatory samples used as the crystallizing agent. It appears that this flexible polyamine molecule has several binding modes, interacting in fairly irregular manner with different sites on DNA and showing no regular ordering in the DNA crystals. Ions of Na(+) and Spm(4+) compete with each other and with water molecules in binding to bases in the minor groove and they influence the structure of the DNA hydration shell in different ways.     

Simulations of the solvent structure for macromolecules. I. Solvation of B-DNA double helix at T = 300 K

Monte Carlo simulations are reported for a system of 447 water molecules enclosing a B-DNA double-helix fragment with 12 base pairs and the corresponding sugar and phosphate units. From a detailed analysis on the interaction energies and probability distributions (at a simulated temperature of 300 K), the water molecules can be partitioned into clusters strongly interacting with (1) the phosphates, (2) the sugars, (3) the sugars and the bases, and (4) the base pairs. In addition, transgroove and interphosphate filament of hydrogen-bonded water molecules have been detected. From simulations performed with variable numbers of water molecules, a theoretical isotherm has been obtained, with the characteristic sigmoidal shape, known from absorption–desorption experiments on related systems. The expected main features for the structure of water molecules solvating B-DNA with Na⁺ counterions are briefly discussed at the end of the paper.     

Penetratin-membrane association: W48/R52/W56 shield the peptide from the aqueous phase

Using molecular dynamics simulations, we studied the mode of association of the cell-penetrating peptide penetratin with both a neutral and a charged bilayer. The results show that the initial peptide-lipid association is a fast process driven by electrostatic interactions. The homogeneous distribution of positively charged residues along the axis of the helical peptide, and especially residues K46, R53, and K57, contribute to the association of the peptide with lipids. The bilayer enhances the stability of the penetratin helix. Oriented parallel to the lipid-water interface, the subsequent insertion of the peptide through the bilayer headgroups is significantly slower. The presence of negatively charged lipids considerably enhances peptide binding. Lateral side-chain motion creates an opening for the helix into the hydrophobic core of the membrane. The peptide aromatic residues form a pi-stacking cluster through W48/R52/W56 and F49/R53, protecting the peptide from the water phase. Interaction with the penetratin peptide has only limited effect on the overall membrane structure, as it affects mainly the conformation of the lipids which interact directly with the peptide. Charge matching locally increases the concentration of negatively charged lipids, lateral lipid diffusion locally decreases. Lipid disorder increases, through decreased order parameters of the lipids interacting with the penetratin side chains. Penetratin molecules at the membrane surface do not seem to aggregate.     

The structural properties of magainin in water, TFE/water, and aqueous urea solutions: molecular dynamics simulations

Here, the MD simulations and comparative structural analysis of Magainin in water, TFE/water, and 2M, 4M, and BM urea solutions is reported. For MAG-TFE/water and MAG-2M urea the largely alpha helical conformation of the peptide is maintained throughout the 9-ns simulation. While in water, 4M urea, and 8M urea, the helix length decreases and at the same time helix radius increases. This suggests a more destabilized magainin secondary structure. Our simulation data reveals that the stabilizing effect of TFE is induced by preferential accumulation of TFE molecules around the alpha helical peptide. These results indicate that an aqueous urea solution solvates the surface of polypeptide chain more favorably than pure water. Urea molecules interact more favorably with nonpolar groups of the peptide in comparison with water, and the presence of urea improves the interactions of water molecules with the hydrophilic groups of the peptide. At 8M urea, there are more direct interactions between the urea and solute, and the helix is destabilized. At 2M urea, the interaction of urea molecules and nonpolar residues are weak, therefore, the presence of urea molecules decreases the interactions of water molecules with hydrophilic groups. Urea could not deteriorate the peptide secondary structure with time from an initial helix structure.     

Molecular dynamics simulations of a winter flounder “antifreeze” polypeptide in aqueous solution

A winter flounder antifreeze polypeptide (HPLC-6) has been studied in vacuo and in aqueous solution using molecular dynamics computer simulation techniques. The helical conformation of this polypeptide was found to be stable both in vacuum and in solution. The major stabilizing interactions were found to be the main-chain hydrogen bonds, a salt-bridge interaction, and solute–solvent hydrogen bonds. A significant bending in the middle of the polypeptide chain was observed both in vacuo and in solvent at 300 K. Possible causes of the bending are discussed. From simulations of mutant polypeptide molecules in vacuo, it is concluded that the bend in the native polypeptide was caused by side chain to backbone hydrogen bond competition involving the Thr 24 side chain and facilitated by strains on the helix resulting from the Lys 18-Glu 22 salt bridge. © 1993 John Wiley & Sons, Inc.     

The Three-dimensional Structure of Two Mutants of the Signal Transduction Protein CheY Suggest its Molecular Activation Mechanism

The three-dimensional crystal structures of the single mutant M17G and the triple mutant F14G-S15G-M17G of the response regulator protein CheY have been determined to 2.3 and 1.9 Å, respectively. Both mutants bind the essential Mg²⁺cation as determined by the changes in stability, but binding does not cause the intrinsic fluorescence quenching of W58 observed in the wild-type protein. The loop β4-α4 appears to be very flexible in both mutants and helix α4, which starts at N94 in the native Mg²⁺-CheY and at K91 in the native apo-CheY, starts in both mutants at residue K92. The side-chain of K109 appears to be more mobile because of the space freed by the M17G mutation. In the triple mutant the main chain of K109 and adjacent residues (loop β5-α5) is displaced almost by 2Å affecting the main chain at residues T87 to E89 (C terminus of β4). The triple mutant structure has a Mg²⁺bound at the active site, but although the Mg²⁺coordination is similar to that of the native Mg²⁺-CheY, the structural consequences of the metal binding are quite different. It seems that the mutations have disrupted the mechanism of movement transmission observed in the native protein. We suggest that the side-chain of K109, packed between V86, A88 and M17 in the native protein, slides forwards and backwards upon activation and deactivation dragging the main chain at the loop β5-α5 and triggering larger movements at the functional surface of the protein.     

The effect of structural water molecules on the normal mode spectrum of dTn . dAn x dTn DNA

In this paper we present a theoretical treatment of triplex B type DNA hydration using normal mode calculation techniques. Discrete solvent is added as spines of hydration in the Watson-Crick and Crick-Hoogsteen grooves as well as water bridges between the Phosphate groups. The effect of binding the discrete structural waters on the normal mode of vibration of the system was studied by introducing a parameter, Xw, that is proportional to the degree of water binding and inversely proportional to the relative humidity (RH) of the system. We examined the variation of the dipole moments of characteristic modes with Xw. The results show that there is a direct relationship between the degree of binding of the water molecules to the atoms in the triple helix, the relative humidity of the system and the conformation and stability of the triple helix. At high RH and Xw = 0:0 the triple helix has mostly B type conformation characteristics, with C'2 -endo sugars. The emergence of normal modes of vibration characteristic to the A type conformation (C'3 - endo sugars) at Xw = 0:4 and 60% RH indicates a conformational shift towards A-type for some of the sugars between Xw = 0.2 (80% RH) and Xw = 0.4 (60% RH). These results are in agreement with the "economy of hydration hypothesis" of Saenger (Saenger et al., 1986) which maintains that the main difference in the hydration of A- and B- forms of DNA is the presence of water bridges between adjacent Phosphate groups in the low-hydration A-form but not in the B- form. Free energy calculations for the triplex DNA with structural waters show that there is a minimum of the free energy at Xw = 0.2 and the free energy increases with Xw and becomes larger than the free energy of the B conformation without structural waters for Xw equal to and larger than 0.4. This result indicates that the B conformation is more stable with bound structural water molecules (for degrees of water binding that are not over 20% higher than the degree of binding between bulk water molecules). The structural water molecules are bound much tighter in the A conformation than in the B conformation. The model predicts that the B to A transition occurs at higher relative humidities in D2O than in H2O. Part of these results (Dadarlat, 1997) have been subsequently confirmed by the experimental work and MD simulations of Ouali (Ouali et al., 1997). The experimental results showed that the N-type sugars corresponding to the A conformation are clearly detected below 75% RH.     

Hydration of the dTn.dAn x dTn parallel triple helix: a Fourier transform infrared and gravimetric study correlated with molecular dynamics simulations.

We present a comparative analysis of the water organization around the dTn.dAn x dTn triple helix and the Watson-Crick double helix dTn.dAn respectively by means of gravimetric measurements, infrared spectroscopy and molecular dynamics simulations. The hydration per nucleotide determined by gravimetric and spectroscopic methods correlated with the molecular dynamics simulations shows that at high relative humidity (98% RH) the triple helix is less solvated than the duplex (17 +/- 2 water molecules per nucleotide instead of 21 +/-1). The experimental desorption curves are different for both structures and indicate that below 81% RH the triplex becomes more hydrated than the duplex. At this RH the FTIR spectra show the emergence of N-type sugars in the adenosine strand of the triplex. When the third strand is bound in the major groove of the Watson-Crick duplex molecular dynamics simulations show the formation of a spine of water molecules between the two thymidine strands.