Effects of osmolytes on the helical conformation of model peptide: molecular dynamics simulation

Co-solvents such as glycerol and sorbitol are small organic molecules solvated in the cellular solutions that can have profound effects on the protein structures. Here, the molecular dynamics simulations and comparative structural analysis of magainin, as a peptide model, in pure water, 2,2,2-trifluoroethanol∕water, glycerol∕water, and sorbitol∕water are reported. Our results show that the peptide NMR structure is largely maintained its native structure in osmolytes-water mixtures. The simulation data indicates that the stabilizing effect of glycerol and sorbitol is induced by preferential accumulation of glycerol and sorbitol molecules around the nonpolar and aromatic residues. Thus, the presence of glycerol and sorbitol molecules decreases the interactions of water molecules with the hydrophobic residues of the peptide, and the alpha helical structure is stabilized.     

Interactions between hydrophobic and ionic solutes in aqueous guanidinium chloride and urea solutions: lessons for protein denaturation mechanism

In order to clarify the mechanism of denaturant-induced unfolding of proteins we have calculated the interactions between hydrophobic and ionic species in aqueous guanidinium chloride and urea solutions using molecular dynamics simulations. Hydrophobic association is not significantly changed in urea or guanidinium chloride solutions. The strength of interaction between ion pairs is greatly diminished by the guanidinium ion. Although the changes in electrostatic interactions in urea are small, examination of structures, using appropriate pair functions, of urea and water around the solutes show strong hydrogen bonding between urea's carbonyl oxygen and the positively charged solute. Our results strongly suggest protein denaturation occurs by the direct interaction model according to which the most commonly used denaturants unfold proteins by altering electrostatic interactions either by solvating the charged residues or by engaging in hydrogen bonds with the protein backbone. To further validate the direct interaction model we show that, in urea and guanidinium chloride solutions, unfolding of an unusually stable helix (H1) from mouse PrPC (residues 144-153) occurs by hydrogen bonding of denaturants to charged side chains and backbone carbonyl groups.     

Examination of the folding of a short alanine-based helical peptide with salt bridges using molecular dynamics simulation

A molecular dynamics simulation of the folding of a short alanine-based helical peptide of 17 residues with three Glu...Lys (i, i + 4) salt bridge pairs, referred to as the AEK17 peptide, was carried out. The simulation gave an estimated simulation folding time of 2.5 ns, shorter than 12 ns for an alanine-based peptide of 16 residues with three Lys residues only, referred to as the AK16 peptide, simulated previously. After folded, the AEK17 peptide had a helical content of 77%, in excellent agreement with the experimentally determined value of 80%. An examination of the folding pathways of AEK17 indicated that the peptide proceeded via three-turn helix conformations more than the helix-turn-helix conformation in the folding pathways. An analysis of interactions indicated that the formation of hydrogen bonds between Lys residue side chains and backbone carbonyls is a major factor in the abundant conformation of the three-turn helix intermediate. The substitution of three Ala with Glu residues reduces the extent of hydrophobic interaction in alanine-based AK peptides with the result that the breaking of the interactions of Lys epsilon-NH3+(side chain)...C=O(backbone) is a major activation action for the AEK17 to achieve a complete fold, in contrast to the AK16 peptide, in which breaking non-native hydrophobic interaction is the rate-determining step.     

Discriminating the helical forms of peptides by NMR and molecular dynamics simulation

The HNCO NMR pulse sequence was applied to three selectively labeled (15)N and (13)C isotopic homologues of the peptide Ac-WAAAH(AAARA)(3)A-NH(2) to probe directly for hydrogen bonds between residues 8 and 11 (characteristic of a 3(10)-helix), 8 and 12 (alpha-helix), and 8 and 13 (pi-helix). The experiments demonstrate conclusively, and in agreement with circular dichroism studies, that the center of the peptide is alpha-helical; there is no discernible 3(10)- or pi-helix at these specific positions. Molecular dynamics simulations of the preceding peptide and Ac-(AAAAK)(3)A-NH(2) in water using the potential energy parameter set CHARMM22/CMAP correctly yield an alpha-helix, in contrast to simulations with the set CHARMM22, which result in a pi-helix.     

Thermophoresis in liquids: a molecular dynamics simulation study

Thermophoresis in liquids is studied by molecular dynamics simulation (MD). A theory is developed that divides the problem in the way consistent with the characteristic scales. MD is then conducted to obtain the solution of each problem, which is to be all combined for macroscopic predictions. It is shown that when the temperature gradient is applied to the nonconducting liquid bath that contains neutral particles, there occurs a pressure gradient tangential to the particle surface at the particle-liquid interface. This may induce the flow in the interfacial region and eventually the particle to move. This applies to the material system that interacts through van der Waals forces and may be a general source of the thermophoresis phenomenon in liquids. The particle velocity is linearly proportional to the temperature gradient. And, in a large part of the given temperature range, the particle motion is in the direction toward the cold end and decreases with respect to the temperature. It is also shown that the particle velocity decreases or even reverses its sign in the lowest limit of the temperature range or with a particle of relatively weak molecular interactions with the liquid. The characteristics of the phenomenon are analyzed in molecular details.     

Molecular dynamics simulation of folding of a short helical peptide with many charged residues

A molecular dynamics simulation of the folding of conantokin-T (con-T), a short helical peptide with 5 helical turns of 21 amino acids with 10 charged residues, was carried out to examine folding pathways for this peptide and to predict the folding rate. In the 18 trajectories run at 300 K, 16 trajectories folded, with an averaged folding time of approximately 50 ns. Two trajectories did not fold in up to 200 ns simulation. The folded structure in folded trajectories is in good agreement with experimental structure. An analysis of the trajectories showed that, at the beginning of a few nanoseconds, helix formation started from residues 5-9 with assistance of a hydrophobic clustering involving Tyr5, Met8, and Leu9. The peptide formed a U-shape mainly due to charge-charge interactions between charged residues at the N- and C-terminus segments. In the next approximately 10 ns, several nonnative charge-charge interactions were broken and nonnative Gla10-Lys18 (this denotes a salt bridge between Gal10 and Lys18) and/or Gla10-Lys19 interactions appeared more frequently in this folding step and the peptide became a fishhook J-shape. From this structure, the peptide folded to the folded state in 7 of all 16 folded trajectories in approximately 15 ns. Alternatively, in approximately 30 ns, the con-T went to a conformation in an L-shape with 4 helical turns and a kink at the Arg13 and Gla14 segment in the other 9 trajectories. Con-T in the L-shape then required another approximately 15 ns to fold into the folded state. In addition, in overall folding times, the former 7 trajectories folded faster with the total folding times all shorter than 45 ns, while the latter 9 trajectories folded at a time longer than 45 ns, resulting in an average folding time of approximately 50 ns. Two major folding intermediates found in 2 nonfolded trajectories are stabilized by charge clusters of 5 and 6 charged residues, respectively. With inclusion of friction and solvent-solvent interactions, which were ignored in the present GB/SA solvation model, the folding time obtained above should be multiplied by a factor of 1.25-1.7 according to a previous, similar simulation study. This results in a folding time of 65-105 ns, slightly shorter than the folding time of 127 ns for an alanine-based peptide of the same length. This suggests that the energy barrier of folding for this type of peptides with many charged residues is slightly lower than alanine-based helical peptides by less than 1 kcal/mol.     

Pore-spanning lipid membrane under indentation by a probe tip: a molecular dynamics simulation study

We study the indentation of a free-standing lipid membrane suspended over a nanopore on a hydrophobic substrate by means of molecular dynamics simulations. We find that in the course of indentation the membrane bends at the point of contact and the fringes of the membrane glide downward intermittently along the pore edges and stop gliding when the fringes reach the edge bottoms. The bending continues afterward, and the large strain eventually induces a phase transition in the membrane, transformed from a bilayered structure to an interdigitated structure. The membrane is finally ruptured when the indentation goes deep enough. Several local physical quantities in the pore regions are calculated, which include the tilt angle of lipid molecules, the nematic order, the included angle, and the distance between neighboring lipids. The variations of these quantities reveal many detailed, not-yet-specified local structural transitions of lipid molecules under indentation. The force-indentation curve is also studied and discussed. The results make a connection between the microscopic structure and the macroscopic properties and provide deep insight into the understanding of the stability of a lipid membrane spanning over nanopore.     

Characterizing hydrophobicity at the nanoscale: a molecular dynamics simulation study

We use molecular dynamics (MD) simulations of water near nanoscopic surfaces to characterize hydrophobic solute-water interfaces. By using nanoscopic paraffin like plates as model solutes, MD simulations in isothermal-isobaric ensemble have been employed to identify characteristic features of such an interface. Enhanced water correlation, density fluctuations, and position dependent compressibility apart from surface specific hydrogen bond distribution and molecular orientations have been identified as characteristic features of such interfaces. Tetrahedral order parameter that quantifies the degree of tetrahedrality in the water structure and an orientational order parameter, which quantifies the orientational preferences of the second solvation shell water around a central water molecule, have also been calculated as a function of distance from the plate surface. In the vicinity of the surface these two order parameters too show considerable sensitivity to the surface hydrophobicity. The potential of mean force (PMF) between water and the surface as a function of the distance from the surface has also been analyzed in terms of direct interaction and induced contribution, which shows unusual effect of plate hydrophobicity on the solvent induced PMF. In order to investigate hydrophobic nature of these plates, we have also investigated interplate dewetting when two such plates are immersed in water.     

Structure and dynamics of surfactant and hydrocarbon aggregates on graphite: a molecular dynamics simulation study

We have examined the structure and dynamics of sodium dodecyl sulfate (SDS) and dodecane (C12) molecular aggregates at varying surface coverages on the basal plane of graphite via classical molecular dynamics simulations. Our results suggest that graphite-hydrocarbon chain interactions favor specific molecular orientations at the single-molecule level via alignment of the tail along the crystallographic directions. This orientational bias is reduced greatly upon increasing the surface coverage for both molecules due to intermolecular interactions, leading to very weak bias at intermediate surface coverages. Interestingly, for complete monolayers, we find a re-emergent orientational bias. Furthermore, by comparing the SDS behavior with C12, we demonstrate that the charged head group plays a key role in the aggregate structures: SDS molecules display a tendency to form linear file-like aggregates while C12 forms tightly bound planar ones. The observed orientational bias for SDS molecules is in agreement with experimental observations of hemimicelle orientation and provides support for the belief that an initial oriented layer governs the orientation of hemimicellar aggregates.     

Thermal denaturation of polyalanine peptide in water by molecular dynamics simulations and theoretical prediction of infrared spectra: helix-coil transition kinetics

Perspectives in the helix-coil transition kinetics of secondary structures are examined by temperature-jump molecular dynamics (T-jump MD) simulations and theoretically calculated infrared (IR) spectra. Homopolymeric polyalanine, Ac-(A)(21)-NHMe (A21), is unfolded in water by T-jumps from 273 to 300 K approximately 450 K using AMBER ff99 and ff03 force fields. MD simulation results provide in silico evidence that 3(10)-helix and type I beta-turn motifs are highly probable in both ff99 and ff03 results. Temperature-dependent difference IR spectra of A21 do not possess an isosbestic point in both results, and isotope-labeled difference IR spectra in ff03 results predict characteristic profiles observed in experiments. Unfolding rates obtained from simulated time-resoled IR spectra are in good agreement with those estimated by helical contents, but they are still an order of magnitude smaller than experimental values. We demonstrate that the conventional criteria such as single-exponential fit of transient amide I absorbance, sigmoidal fit of temperature-dependent amide I absorbance, and Arrhenius plot of relaxation rates cannot guarantee the validity of assuming a two-state mechanism. We suggest a way of determining T(m) by the temperature dependence of center frequency and full width at half-maximum of amide I band. Overall, both ff99 and ff03 force fields give consistent results in reproducing key aspects concerned experimentally, but are not predominantly satisfactory in quantitative aspects.     

Interfacial fluctuations of block copolymers: a coarse-grain molecular dynamics simulation study

The lamellar and cylindrical phases of block copolymers have a number of technological applications, particularly when they occur in supported thin films. One such application is block copolymer lithography, the use of these materials to subdivide or enhance submicrometer patterns defined by optical or electron beam methods. A key parameter of all lithographic methods is the line edge roughness (LER), because the electronic or optical activities of interest are sensitive to small pattern variations. While mean-field models provide a partial picture of the LER and interfacial width expected for the block interface in a diblock copolymer, these models lack chemical detail. To complement mean-field approaches, we have carried out coarse-grain molecular dynamics simulations on model poly(ethyleneoxide)-poly(ethylethylene) (PEO-PEE) lamellae, exploring the influence of chain length and hypothetical chemical modifications on the observed line edge roughness. As expected, our simulations show that increasing chi (the Flory-Huggins parameter) is the most direct route to decreased roughness, although the addition of strong specific interactions at the block interface can also produce smoother patterns.     

Thermal decomposition of a honeycomb-network sheet: A molecular dynamics simulation study

The thermal degradation of a graphene-like two-dimensional honeycomb membrane with bonds undergoing temperature-induced scission is studied by means of Molecular Dynamics simulation using Langevin thermostat. We demonstrate that at lower temperature the probability distribution of breaking bonds is highly peaked at the rim of the membrane sheet whereas at higher temperature bonds break at random everywhere in the hexagonal flake. The mean breakage time τ is found to decrease with the total number of network nodes N by a power law τ ∝ N(-0.5) and reveals an Arrhenian dependence on temperature T. Scission times are themselves exponentially distributed. The fragmentation kinetics of the average number of clusters can be described by first-order chemical reactions between network nodes n(i) of different coordination. The distribution of fragments sizes evolves with time elapsed from initially a δ-function through a bimodal one into a single-peaked again at late times. Our simulation results are complemented by a set of 1st-order kinetic differential equations for n(i) which can be solved exactly and compared to data derived from the computer experiment, providing deeper insight into the thermolysis mechanism.     

Contact angle hysteresis at the nanoscale: a molecular dynamics simulation study

In this paper, the contact angle hysteresis (CAH) of nanodroplets on both rigid and flexible substrates with different wettabilities was studied using molecular dynamics (MD) simulations. The critical shear stress (CSS) that determines the motion of the contact line (CL) was investigated. A theoretical correlation between CAH and CSS was proposed. Both CAH and CSS reflect the energy dissipation at the CL of the droplet in response to the exerted force. MD results of CAH are qualitatively consistent with the theoretical model. Simulation results also show that, for the same liquid–solid interactions, CAH on the flexible substrate is larger than that on the rigid substrate. These findings aim to enhance our understanding of the mechanism of the CAH at the nanoscale.     

IBIsCO: a molecular dynamics simulation package for coarse-grained simulation

IBIsCO is a parallel molecular dynamics simulation package developed specially for coarse-grained simulations with numerical potentials derived by the iterative Boltzmann inversion (IBI) method (Reith et al., J Comput Chem 2003, 24, 1624). In addition to common features of molecular dynamics programs, the techniques of dissipative particle dynamics (Groot and Warren, J Chem Phys 1997, 107, 4423) and Lowe-Andersen dynamics (Lowe, Europhys Lett 1999, 47, 145) are implemented, which can be used both as thermostats and as sources of friction to compensate the loss of degrees of freedom by coarse-graining. The reverse nonequilibrium molecular dynamics simulation method (Müller-Plathe, Phys Rev E 1999, 59, 4894) for the calculation of viscosities is also implemented. Details of the algorithms, functionalities, implementation, user interfaces, and file formats are described. The code is parallelized using PE_MPI on PowerPC architecture. The execution time scales satisfactorily with the number of processors.     

A molecular dynamics study of chloride binding by the cryptand SC24

The capture of chloride from water by the tetraprotonated form of the spherical macrotricyclic molecule SC24 was studied using molecular dynamics simulation methods. This model ionophore represents a broad class of molecules which remove ions from water. Two binding sites for the chloride were found, one inside and one outside the ligand. These sites are separated by a potential energy barrier of approximately 20 kcal mol⁻¹. The major contribution to this barrier comes from dehydration of the chloride. The large, unfavorable dehydration effect is compensated for by an increase in electrostatic attraction between the oppositely charged chloride and cryptand, and by energetically favorable rearrangements of water structure. Additional assistance in crossing the barrier and completing the dehydration of the ion is provided by the shift of three positively charged hydrogen atoms of the cryptand towards the chloride. This structural flexibility of the cryptand, leads to a decrease in the energy barrier, whereas, its structural rigidity is partially responsible for its selectivity.     

Crystallization of alkane melts induced by carbon nanotubes and graphene nanosheets: a molecular dynamics simulation study

The crystallization of alkane melts on carbon nanotubes (CNT) and the surface of graphene nanosheets (GNS) is investigated using molecular dynamics (MD) simulations. The crystallization process of the alkane melts is analyzed in terms of the bond-orientational order parameter, atomic radial distribution for the CNT/alkane, atomic longitudinal distribution for the GNS/alkane, and diffusion properties. The dimensional effects of the different carbon-based nanostructures on the crystallization of alkane melts are shown. It is found that one-dimensional CNT has a stronger ability to induce the crystallization of the polymer than that of two-dimensional GNS, which provides a support at molecular level for the experimental observation [Li et al., J. Am. Chem. Soc., 2006, 128, 1692 and Xu et al., Macromolecules, 2010, 43, 5000]. From the MD simulations, we also find that the crystallization of alkane molecules has been completed with the highly cooperative processes of adsorption and orientation.     

Molecular inclusion of PCB126 by beta-cyclodextrin: a combined molecular dynamics simulation and quantum chemical study

The effective enrichment and identification of lowly concentrated polychlorinated biphenyls (PCBs) in the environment is attracting much research attention due to human health concerns raised from their emissions. Cyclodextrins (CDs) are known to be capable to form inclusion complexes with a variety of organic molecules. The purpose of this study is to provide theoretical evidences whether CDs can form energetically stable inclusion complexes with PCBs through a host–guest interaction, and if so, whether infrared and Raman techniques are suitable for the detection of CD-modified PCBs. Focusing on a representative PCB molecule, 3,3′,4,4′,5-pentachlorobiphenyl (PCB126), we studied its molecular inclusion by β-CD (BCD) by performing molecular dynamics simulations and density functional theory calculations. Calculated results show that PCB126 and BCD preferentially form the stable 1:1 inclusion complex. The calculated IR spectra of the 1:1 inclusion complexes mainly present the spectra features of BCD and give only a slight indication for bands of the guest molecule. In contrast, the characteristic vibration modes of the guest molecule are remarkably prominent in the Raman spectra of the inclusion complexes. Based on the present results, we propose that BCD can potentially serve as a candidate for including PCB126 to form the stable 1:1 host–guest complex, and that Raman spectroscopy technology is expected to be suitable for the identification of the CD-modified PCBs, whereas IR spectroscopy is not feasible for such an application.