Molecular dynamics study of water penetration in staphylococcal nuclease

The ionization properties of Lys and Glu residues buried in the hydrophobic core of staphylococcal nuclease (SN) suggest that the interior of this protein behaves as a highly polarizable medium with an apparent dielectric constant near 10. This has been rationalized previously in terms of localized conformational relaxation concomitant with the ionization of the internal residue, and with contributions by internal water molecules. Paradoxically, the crystal structure of the SN V66E variant shows internal water molecules and the structure of the V66K variant does not. To assess the structural and dynamical character of interior water molecules in SN, a series of 10-ns-long molecular dynamics (MD) simulations was performed with wild-type SN, and with the V66E and V66K variants with Glu66 and Lys66 in the neutral form. Internal water molecules were identified based on their coordination state and characterized in terms of their residence times, average location, dipole moment fluctuations, hydrogen bonding interactions, and interaction energies. The locations of the water molecules that have residence times of several nanoseconds and display small mean-square displacements agree well with the locations of crystallographically observed water molecules. Additional, relatively disordered water molecules that are not observed crystallographically were found in internal hydrophobic locations. All of the interior water molecules that were analyzed in detail displayed a distribution of interaction energies with higher mean value and narrower width than a bulk water molecule. This underscores the importance of protein dynamics for hydration of the protein interior. Further analysis of the MD trajectories revealed that the fluctuations in the protein structure (especially the loop elements) can strongly influence protein hydration by changing the patterns or strengths of hydrogen bonding interactions between water molecules and the protein. To investigate the dynamical response of the protein to burial of charged groups in the protein interior, MD simulations were performed with Glu66 and Lys66 in the charged state. Overall, the MD simulations suggest that a conformational change rather than internal water molecules is the dominant determinant of the high apparent polarizability of the protein interior.     

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.     

Molecular dynamics study of aquaporin-1 water channel in a lipid bilayer

The aquaporin-1 water channel was modeled in a palmitoyl-oleoyl-phosphatidyl-choline lipid bilayer, by means of molecular dynamics simulations. Interaction of the protein with the membrane and inter-monomer interactions were analyzed. Structural features of the channel important for its biological function, including the Asn-Pro-Ala (NPA) motifs, and the diffusion of water molecules into the channels, were investigated. Simulations revealed the formation of single file water inside the channels for certain relative positions of the NPA motifs.     

Molecular dynamics study of interaction of dimyristoyl phosphotidyl choline with water

We present here results of molecular dynamics (MD) simulations on hydrated bilayers of 40 molecules of 1-2-dimyristoyl-sn-glycero-3-phosphatidyl choline (DMPC) in liquid crystalline (Lα) phase using two different models (i) with same (A) conformation for all DMPC molecules, (ii) with alternate rows having different (A and B reported in crystallographic studies on DMPC) conformations. The bilayers were hydrated using 776 and 1064 water molecules. Simulations have been carried out at 310K with AMBER 4.0 program, using united atom force field for 200 pico seconds (ps) after equilibration. During heating and equilibration constant pressure temperature (PT) conditions were maintained while in simulation of equillibrated bilayers constant volume temperature (VT) conditions were used. Subaveraged atomic coordinates were used to calculate geometric parameters of lipid molecules and lipid water interaction. Our results show larger flexibility of polar head group and glycerol region in Lα phase compared to gel or non-hydrated bilayers. Chain disorder was more towards end. Sn-2 chains were more disordered. Use of two types of starting conformations increased disorder. Trans fraction of chain torsional angle was higher in non-hydrated bilayer. However it was more disordered due to ‘swing’ movement of chains because of distortion in torsional angles α2 and 03 due to absence of water molecules. Trans fraction of the chains, order parameter and water penetration showed general agreement with the available experimental results. On the whole MD technique was found to be quite useful for depicting microscopic behaviour of liquid crystalline system and correlating the same with macroscopic changes observed experimentally.     

Molecular dynamics study of onset of water gelation around the collagen triple helix

The onset of water gelation around a collagen-like triple helix peptide was studied at ambient temperature and pressure by performing Molecular Dynamics simulations. The radial distribution functions of the oxygen and hydrogen atoms of water are distorted below 4 A from the peptide. The distortion is accompanied by the breakdown of the tetrahedral coordination of the hydrogen-bonded network of water molecules. The water shell around the peptide consists of alternating regions of higher and lower density. In agreement with experiments we find that the first hydration shell is kinetically labile, with a residence time in the order of picoseconds for a water molecule. From the computed diffusion coefficient, a key measure of the collective dynamics, we estimate the average diffusion speed decreases by a factor of 1.5 close to the peptide compared to the liquid. Our results give new insight in gel formation and structure on a molecular level.     

Molecular dynamics of water in the neighborhood of aquaporins

Present knowledge obtained by molecular dynamics (MD) simulation studies regarding the dynamics of water, both in the vicinity of biological membranes and within the proteinaceous water channels, also known as aquaporins (AQPs), is reviewed. A brief general summary of the water models most extensively employed in MD simulations (SPC, SPC/E, TIP3P, TIP4P), indicating their most relevant pros and cons, is likewise provided. Structural considerations of water are also discussed, based on different order parameters, which can be extracted from MD simulations as well as from experiments. Secondly, the behaviour of water in the neighbourhood of membranes by means of molecular dynamics simulations is addressed. Consequently, the comparison with previous experimental evidence is pointed out. In living cells, water is transported across the plasma membrane through the lipid bilayer and the aforementioned AQPs, which motivates this review to focus mostly on MD simulation studies of water within AQPs. Relevant contributions explaining peculiar properties of these channels are discussed, such as selectivity and gating. Water models used in these studies are also summarised. Finally, based on the information presented here, further MD studies are encouraged.     

Molecular dynamics study of peptide-bilayer adsorption

Two 6-ns simulations of the somatostatin analog sandostatin and a 1-palmityl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer are presented. In the first simulation, the peptide was placed in a region of bulk water density and allowed to spontaneously move toward and bind to the bilayer surface. An attractive force between the peptide and bilayer drove the binding process, which was opposed by a significant frictional force caused by the solvent (water). During the approach of the peptide toward the bilayer the area of the interacting surface between the species was inversely proportional to the distance between them, supporting the application of such a relationship in continuum calculations of peptide-bilayer binding free energies. In the second simulation, the N-terminus of the surface-bound peptide was deprotonated. Consistent with experiment, this strengthened interactions between the peptide and the bilayer. Details of both peptide-bilayer complexes, including the orientation, percent buried surface area, and orientation of the lipid headgroups are in good agreement with those obtained from experiment. The location of the different side chains in the bilayer is in direct correlation with an interfacial hydrophobicity scale developed using model peptides. The aromatic side chains of the Phe and Trp residues all lie flat with respect to the bilayer surface in both complexes. Changes in lipid and water ordering due to peptide binding suggest a possible domination of lipophobic over hydrophobic effects, as proposed by other workers. Where appropriate, peptide and lipid properties in the bound states are compared with separate simulations of sandostatin and the bilayer in water, respectively, so as to monitor the response of the system to the binding process.     

Molecular dynamics study of talin-vinculin binding

Cells can sense mechanical force in regulating focal adhesion assembly. One vivid example is the force-induced recruitment of vinculin to reinforce initial contacts between a cell and the extracellular matrix. Crystal structures of the unbound proteins and bound complex between the vinculin head subdomain (Vh1) and the talin vinculin binding site 1 (VBS1) indicate that vinculin undergoes a conformational change upon binding to talin. However, the molecular basis for this event and the precise nature of the binding pathway remain elusive. In this article, molecular dynamics is used to investigate the binding mechanism of Vh1 and VBS1 under minimal constraints to facilitate binding. One simulation demonstrates binding of the two molecules in the complete absence of external force. VBS1 makes early hydrophobic contact with Vh1 by positioning the critical hydrophobic residues (L608, L615, and L622) in the groove formed by helices 1 and 2 of Vh1. The solvent-exposed hydrophobic residues (V619 and L623) then gradually penetrate the hydrophobic core of Vh1, thus further separating helix 1 from helix 2. These critical residues are highly conserved as large hydrophobic side groups in other vinculin binding sites; studies also have demonstrated that these residues are essential in Vh1-VBS1 binding. Similar binding mechanisms are also demonstrated in separate molecular dynamics simulations of Vh1 binding to other vinculin binding sites both in talin and α-actinin.     

Molecular dynamics simulations of Leu-enkephalin in water and DMSO.

The structure of Leu-enkephalin (L-Enk) and Met-enkephalin (M-Enk) have frequently been studied, in particular by nuclear magnetic resonance spectroscopy. After more than 20 years of research, it was concluded that enkephalins have no preferred structure in aqueous solution, but that they may have in other solvents. We have performed molecular dynamics simulations of zwitterionic L-Enk in water, and zwitterionic as well as neutral L-Enk dimethyl sulfoxide (DMSO). In water the peptide is very flexible, although there seems to be a preference for compact conformations. In DMSO, the peptide forms a clear salt bridge in the zwitterionic form, but has no preferred conformation in the neutral form. This difference in conformation may provide an explanation for measurements in DMSO in which multiple conformations were found to exist. In this paper we introduce a new formulation for a dihedral angle autocorrelation function, and apply it to study side-chain dynamics in L-Enk. We find that the side-chain dynamics of the large Tyr and Phe residues cannot be adequately sampled in 2.0-ns simulations, while this does seem to be possible for the smaller Leu side chain.     

Molecular dynamics study of iduronate ring conformation

Molecular dynamics (MD) simulations of the conformation of the iduronate ring in a methyl glycoside and as the central residue in a trisaccharide have been carried out. Separate simulations were carried out with initial ¹C₄, ²S₀, and ⁴C₁ iduronate ring conformations. Simulations were followed by observing the time development of the Cremer–Pople ring puckering parameters θ,?₂. Starting with chair geometries gave trajectories showing only ring oscillations close to the initial geometry. Simulations were performed with a ²S₀ starting geometry using explicit water and in vacuum with dielectric constants (ε) of 1 and 80, as well as with distance-dependent dielectric functions of 2r and 4r. In both the explicit water simulation and the vacuum (ε = 80) simulations, extensive pseudorotational motion was observed in which boat and twist-boat ring conformers interconvert. The overall range of ?₂2 variation in the trisaccharide was about half of that observed in the methyl glycoside. The Haasnoot procedure for calculating H-H coupling constants in saccharides was applied to structures obtained from MD trajectories. Using MD time averaged couplings along with experimental data allowed the relative fractions of chair and boat/twist-boat forms to be derived. © 1993 John Wiley & Sons, Inc.     

Molecular dynamics simulation of deoxy and carboxy murine neuroglobin in water

Globins are respiratory proteins that reversibly bind dioxygen and other small ligands at the iron of a heme prosthetic group. Hemoglobin and myoglobin are the most prominent members of this protein family. Unexpectedly a few years ago a new member was discovered and called neuroglobin (Ngb), being predominantly expressed in the brain. Ngb is a single polypeptide of 151 amino acids and despite the small sequence similarity with other globins, it displays the typical globin fold. Oxygen, nitric oxide, or carbon monoxide can displace the distal histidine which, in ferrous Ngb as well as in ferric Ngb, is bound to the iron, yielding a reversible adduct. Recent crystallographic data on carboxy Ngb show that binding of an exogenous ligand is associated to structural changes involving heme sliding and a topological reorganization of the internal cavities; in particular, the huge internal tunnel that connects the bulk with the active site, peculiar to Ngb, is heavily reorganized. We report the results of extended (90 ns) molecular dynamics simulations in water of ferrous deoxy and carboxy murine neuroglobin, which are both coordinated on the distal site, in the latter case by CO and in the former one by the distal His(64)(E7). The long timescale of the simulations allowed us to characterize the equilibrated protein dynamics and to compare protein structure and dynamical behavior coupled to the binding of an exogenous ligand. We have characterized the heme sliding motion, the topological reorganization of the internal cavities, the dynamics of the distal histidine, and particularly the conformational change of the CD loop, whose flexibility depends ligand binding.     

Molecular dynamics simulations of interaction between sub-bituminous coal and water

The high moisture content of sub-bituminous coal is associated with the interactions between coal and water. Because of complex composition and structure, the graphite surface modified by hydroxyl, carboxyl and carbonyl groups was used to represent the surface model of sub-bituminous coal according to XPS results. Density profiles for oxygen atoms and hydrogen atoms indicate that the coal surface properties affect the structural and dynamic characteristics of the interfacial water molecules. The interfacial water exhibits much more ordering than bulk water. The results of radial distribution functions, mean square displacement and local self-diffusion coefficient for water molecule related to three oxygen moieties confirmed that the water molecules prefer to absorb with carboxylic groups, and adsorption of water molecules at the hydroxy and carbonyl is similar.     

Molecular dynamics of solid-state lysozyme as affected by glycerol and water: a neutron scattering study

Glycerol has been shown to lower the heat denaturation temperature (T(m)) of dehydrated lysozyme while elevating the T(m) of hydrated lysozyme (. J. Pharm. Sci. 84:707-712). Here, we report an in situ elastic neutron scattering study of the effect of glycerol and hydration on the internal dynamics of lysozyme powder. Anharmonic motions associated with structural relaxation processes were not detected for dehydrated lysozyme in the temperature range of 40 to 450K. Dehydrated lysozyme was found to have the highest T(m) by. Upon the addition of glycerol or water, anharmonicity was recovered above a dynamic transition temperature (T(d)), which may contribute to the reduction of T(m) values for dehydrated lysozyme in the presence of glycerol. The greatest degree of anharmonicity, as well as the lowest T(d), was observed for lysozyme solvated with water. Hydrated lysozyme was also found to have the lowest T(m) by. In the regime above T(d), larger amounts of glycerol lead to a higher rate of change in anharmonic motions as a function of temperature, rendering the material more heat labile. Below T(d), where harmonic motions dominate, the addition of glycerol resulted in a lower amplitude of motions, correlating with a stabilizing effect of glycerol on the protein.     

Molecular dynamics simulation of water confined in a nanopore of amorphous silica

Molecular dynamics simulations are performed to study the transport and structural properties of water confined in a cylindrical silica nanopore. The pore wall is amorphous and mimics a typical mesoporous silica material. The diameters of silica pores studied are 4.75, 9.51, 20 and 25 Å. The self-diffusion of water calculated decreases with pore size and indicates much slower transport compared to the bulk phase. Strong adsorption of water to the silica wall is observed in the density profiles, indicating the hydrophilic nature of the wall. The hydrogen-bonding network is strongly affected by water–silica wall interaction. The average number of hydrogen bonds per water decreased with decreasing pore diameter.     

Molecular dynamics simulation of water permeation through the alpha-hemolysin channel

The alpha-hemolysin (AHL) nanochannel is a non-selective channel that allows for uncontrolled transport of small molecules across membranes leading to cell death. Although it is a bacterial toxin, it has promising applications, ranging from drug delivery systems to nano-sensing devices. This study focuses on the transport of water molecules through an AHL nanochannel using molecular dynamics (MD) simulations. Our results show that AHL can quickly transport water across membranes. The first-passage time approach was used to estimate the diffusion coefficient and the mean exit time. To study the energetics of transport, the potential of mean force (PMF) of a water molecule along the AHL nanochannel was calculated. The results show that the energy barriers of water permeation across a nanopore are always positive along the channel and the values are close to thermal energy (k_BT). These findings suggest that the observed quick permeation of water is due to small energy barriers and a hydrophobic inner channel surface resulting in smaller friction. We speculate that these physical mechanisms are important in how AHL causes cell death.     

Molecular dynamics study of bipolar tetraether lipid membranes

Membranes composed of bipolar tetraether lipids have been studied by a series of 25-ns molecular dynamics simulations to understand the microscopic structure and dynamics as well as membrane area elasticity. By comparing macrocyclic and acyclic tetraether and diether archaeal lipids, the effect of tail linkage of the two phytanyl-chained lipids on the membrane properties is elucidated. Tetraether lipids show smaller molecular area and lateral mobility. For the latter, calculated diffusion coefficients are indeed one order-of-magnitude smaller than that of the diether lipid. These two tetraether membranes are alike in many physical properties except for membrane area elasticity. The macrocyclic tetraether membrane shows a higher elastic area expansion modulus than its acyclic counterpart by a factor of three. Free energy profiles of a water molecule crossing the membranes show no major difference in barrier height; however, a significant difference is observed near the membrane center due to the lack of the slip-plane in tetraether membranes.     

Molecular dynamics study of the KcsA potassium channel

The structural, dynamical, and thermodynamic properties of a model potassium channel are studied using molecular dynamics simulations. We use the recently unveiled protein structure for the KcsA potassium channel from Streptomyces lividans. Total and free energy profiles of potassium and sodium ions reveal a considerable preference for the larger potassium ions. The selectivity of the channel arises from its ability to completely solvate the potassium ions, but not the smaller sodium ions. Self-diffusion of water within the narrow selectivity filter is found to be reduced by an order of magnitude from bulk levels, whereas the wider hydrophobic section of the pore maintains near-bulk self-diffusion. Simulations examining multiple ion configurations suggest a two-ion channel. Ion diffusion is found to be reduced to approximately (1)/(3) of bulk diffusion within the selectivity filter. The reduced ion mobility does not hinder the passage of ions, as permeation appears to be driven by Coulomb repulsion within this multiple ion channel.