Molecular dynamics simulation of diffusion of hydrogen in binary hydrogen–tetrahydrofuran hydrate

The binary structure II hydrogen–tetrahydrofuran (THF) hydrate was studied with molecular dynamics simulation. The simulations were carried out at 300, 310 K and 10.1 MPa, and with various contents of hydrogen and THF. The migrations of hydrogen molecules from cage to cage were observed. The migration process of hydrogen was also analysed, and the diffusion coefficients of hydrogen in the hydrate were calculated. The calculated diffusion coefficients qualitatively agreed with the experimental data. Double and quintet occupancies of hydrogen molecules were observed in the small and large cages, respectively, without changing the hydrate structure.     

Quantum chemical molecular dynamics simulation of carbon nanotube–graphene fusion

Abstract

Here we report a quantum mechanical molecular dynamics (QM/MD) study of a fusion process of an open-ended carbon nanotube on a graphene hole, which results in the formation of a so-called pillared graphene structure – a three-dimensional nanomaterial consisting entirely of sp²-carbons. The self-consistent-charge density-functional tight-binding potential was adopted in this study. Two different sizes of graphene holes with 12 or 24 central carbon atoms removed from a graphene flake, and a (6,6) carbon nanotube with a compatible diameter were adopted. Formations of 6–7–6/5–8–5 defect structures were found on the fusion border between tube and graphene hole. The 6–7–6 structure was found to bear less curvature-induced strain energy and therefore to be more stable and much easier to form than the 5–8–5 structure.     

Molecular dynamics simulation of Axillaridine–A: A potent natural cholinesterase inhibitor

Molecular Dynamics (MD) simulations were carried out for human acetylcholinesterase (hAChE) and its complex with Axillaridine–A, in order to dynamically explore the active site of the protein and the behaviour of the ligand at the peripheral binding site. Simulation of the enzyme alone showed that the active site of AChE is located at the bottom of a deep and narrow cavity whose surface is lined with rings of aromatic residues while Tyr72 is almost perpendicular to the Trp286, which is responsible for stable π -π interactions. The complexation of AChE with Axillaridine-A, results in the reduction of gorge size due to interaction between the ligand and the active site residues. The gorge size was determined by the distance between the center of mass of Glu81 and Trp286. As far as the geometry of the active site is concerned, the presence of ligand in the active site alters its specific conformation, as revealed by stable hydrogen bondings established between amino acids. With the increasing interaction between ligand and the active amino acids, size of the active site of the complex decreases with respect to time. Axillaridine-A, forms stable π -π interactions with the aromatic ring of Tyr124 that results in inhibition of catalytic activity of the enzyme. This π -π interaction keeps the substrate stable at the edge of the catalytic gorge by inhibiting its catalytic activity. The MD results clearly provide an explanation for the binding pattern of bulky steroidal alkaloids at the active site of AChE.     

Molecular dynamics simulation of the graphene–water interface: comparing water models

In this work, different water models (i.e. SPC/E, TIP3P, TIP4P/2005, TIP5P, SPC/Fw, TIP4P/2005f and SWM4_DP) are implemented to simulate water on neutral, negatively charged and positively charged graphene. In all cases ambient conditions are considered. Structural and dynamical properties for water are calculated to quantify the differences among various water models. The results show that SPC/E, TIP4P/2005, SPC/Fw, TIP4P/2005f and SWM4_DP water models yield a similar structure for interfacial water on graphene, whether it is neutral, negatively charged or positively charged. TIP5P is the model whose predictions for the structure of the interface deviate the most from those of the other models. Although qualitatively the results are for the most part similar, a large quantitative variation is observed among the dynamical properties predicted when various water models are implemented. Although experimental data are not available to discriminate the most/least accurate of the model predictions, our results could be useful for comparing results for interfacial water obtained implementing different models. Such critical comparison will benefit practical applications such as the development of energy-storage and water-desalination devices (e.g. electric double-layer capacitors), among others.     

What is the role of SNARE proteins in membrane fusion?

Soluble N-ethylmaleimide-sensitive factor activating protein receptor (SNARE) proteins have been at the fore-front of research on biological membrane fusion for some time. The subcellular localization of SNAREs and their ability to form the so-called SNARE complex may be integral to determining the specificity of intracellular fusion (the SNARE hypothesis) and/or serving as the minimal fusion machinery. Both the SNARE hypothesis and the idea of the minimal fusion machinery have been challenged by a number of experimental observations in various model systems, suggesting that SNAREs may have other functions. Considering recent advances in the SNARE literature, it appears that SNAREs may actually function as part of a complex fusion "machine." Their role in the machinery could be any one or a combination of roles, including establishing tight membrane contact, formation of a scaffolding on which to build the machine, binding of lipid surfaces, and many others. It is also possible that complexations other than the classic SNARE complex participate in membrane fusion.     

Fluorescent membrane probes’ behavior in lipid bilayers: insights from molecular dynamics simulations

Fluorescence spectroscopy and microscopy have been used as tools to study membrane biophysics for decades now. Because phospholipids are non-fluorescent, the use of extrinsic membrane probes in this context is commonplace. Two major points of concern arise regarding this matter, namely the incomplete understanding of the probe behavior inside the bilayer and the perturbation of the latter resulting from probe incorporation. To this effect, molecular dynamics (MD) simulations, by providing detailed atomic-scale information, represent a valuable way to characterize the location and dynamics of bilayer-inserted membrane probes, as well as the magnitude of perturbation they induce on the host lipid structure, and several important classes of reporter molecules have been studied in recent years. This article reviews the state of the art of MD simulations of bilayer-inserted fluorescent probes, focusing on the information that has been obtained from previous studies and hinting at future perspectives in this rapidly emerging field.     

Molecular dynamics simulation of β₂-microglobulin in denaturing and stabilizing conditions

β₂‐Microglobulin has been a model system for the study of fibril formation for 20 years. The experimental study of β₂‐microglobulin structure, dynamics, and thermodynamics in solution, at atomic detail, along the pathway leading to fibril formation is difficult because the onset of disorder and aggregation prevents signal resolution in Nuclear Magnetic Resonance experiments. Moreover, it is difficult to characterize conformers in exchange equilibrium. To gain insight (at atomic level) on processes for which experimental information is available at molecular or supramolecular level, molecular dynamics simulations have been widely used in the last decade. Here, we use molecular dynamics to address three key aspects of β₂‐microglobulin, which are known to be relevant to amyloid formation: (1) 60 ns molecular dynamics simulations of β₂‐microglobulin in trifluoroethanol and in conditions mimicking low pH are used to study the behavior of the protein in environmental conditions that are able to trigger amyloid formation; (2) adaptive biasing force molecular dynamics simulation is used to force cis‐trans isomerization at Proline 32 and to calculate the relative free energy in the folded and unfolded state. The native‐like trans‐conformer (known as intermediate 2 and determining the slow phase of refolding), is simulated for 10 ns, detailing the possible link between cis‐trans isomerization and conformational disorder; (3) molecular dynamics simulation of highly concentrated doxycycline (a molecule able to suppress fibril formation) in the presence of β₂‐microglobulin provides details of the binding modes of the drug and a rationale for its effect. Proteins 2011. © 2010 Wiley‐Liss, Inc.     

Molecular dynamics simulation of the melting processes of core–shell and pure nanoparticles

Core–shell nanoparticles are nanosized particles that consist of a core and a shell, constructed from different metallic elements. Core–shell nanoparticles have received extensive attention, owing to their various potential applications such as paints, optical films and catalysts. Herein, we investigate the melting behaviours of different core–shell nanoparticles under continuous heating using molecular dynamics simulation. Different metallic elements were examined as core–shell and pure nanoparticles. Five different processes were observed during the melting of core–shell nanoparticles. In contrast, only one process was identified during the melting of pure nanoparticles. These processes were influenced by the nanoparticle size, shell thickness and differences between the lattice constants and melting point temperatures of the metallic elements. Our simulation provides microscopic insights into the melting behaviours of existing and proposed core–shell nanoparticles that would be highly beneficial towards the fabrication of materials with different chemical coatings.     

Molecular dynamics simulation of the gGAPDH–NAD+ complex from Trypanosoma cruzi

A 50-ns molecular dynamics simulation has been used to study the homotetramer of the enzyme glycosomal glyceraldehyde 3-phosphate dehydrogenase (gGAPDH) complexes, from Trypanosoma cruzi, with nicotinamide adenine dinucleotide (NAD⁺) cofactors in aqueous solution. The root mean square deviation indicates that the overall structure of the homotetramer does not undergo significant change. The largest structural change observed was in the NAD⁺ binding domain of subunit (chain) D; as a consequence, the NAD⁺ cofactor was dislocated from its initial position. However, the other subunits were not affected, suggesting that the gGAPDH enzyme exhibits non-cooperative behaviour. Our simulation estimates that the NAD⁺ binding domain rotates about 4.8° relative to the catalytic domain in the apo–holo form transition. The hydrogen bond analysis reveals that the residues R12, I13, D38 and M39 are essential for gGAPDH–NAD⁺ interaction. Furthermore, two promising cavities to be explored in drug design were found: one formed by residues I13, R12, T197, T199, E336 and Y339, and the other by residues C166, H194, R249, I13, R12, T197, T199, E336 and Y339. The results presented in this paper offer new insight into the search for inhibitors of the gGAPDH enzyme of T. cruzi protozoan.     

Molecular dynamics simulation of helium–argon gas mixture under various wall conditions

As computational capabilities increase, molecular dynamics (MD) simulations become important tools of simulating reality. These simulations are especially useful for compressible gas mixture problems. In this study, binary diffusion of helium and argon was examined using a hard-sphere MD simulation method. For the sake of computational speed, low spacing ratios were chosen. Binary mass diffusion of gases in two equally sized halves of a box was simulated for identical initial kinetic energies and number densities. It has been noted that a purely mass diffusion mechanism of different gases is not physically possible. The resultant gas mixtures of several diffusion simulations were used as initial conditions for combined heat transfer – Couette flow, and heating and cooling experiments. The results showed the interesting behaviour of the mixture, which was subjected to various wall conditions. Energy of heavier molecules is found to be more sensitive to the wall velocities and less sensitive to the wall temperatures than lighter molecules. Diffusion, heat transfer, viscosity and heat capacity coefficients are deduced as well.     

Molecular dynamics simulation of α-melanocyte stimulating hormone in a water-membrane model interface

The conformation of the tridecapeptide α-melanocyte stimulating hormone in the presence of a double water-membrane interface was studied by molecular dynamics simulation, using the computational package THOR. In this program the solvent is represented by a continuous medium with dielectric constant ɛ, and the interface between different media is simulated by a surface of discontinuity of the dielectric constant. The electrostatic image method was used to write down the terms, added to the force field, that describe the polarisation effects induced in the interface by the atomic charges. The program was further improved by the introduction of a second surface, parallel to the first one, to mimic the membrane. A conformational search using the software Prelude was employed to find an initial geometry for the peptide in water. The molecular dynamics simulation performed during 10 ns showed that the peptide structure is flexible in water, without stabilisation of any preferential conformation. In the presence of the model membrane, the peptide moved to the medium representing the interior of the membrane. Inside the low dielectric constant medium, the structure of the peptide showed a turn in the central sequence of amino acids and a packed conformation remained stabilised during more than 7.0 ns of simulation. Received: 27 November 1998 / Revised version: 11 March 1999 / Accepted: 8 April 1999     

Molecular dynamics simulation of diffusion behaviour of gas molecules within oil–paper insulation system

The diffusion behaviour of hydrogen, carbon monoxide, carbon dioxide, methane, acetylene, ethylene and ethane in oil and paper medium was examined using molecular dynamics to reveal the diffusion mechanism of gas molecules in transformer oil–paper insulation system at the microscopic level. These compounds are commonly used in the dissolved gas analysis of power transformers and produced during the ageing process of oil–paper composite insulating material. Two groups of models were constructed using molecular dynamics simulation software to simulate the diffusion behaviour of the aforementioned seven types of small gas molecules in oil and paper. The diffusion coefficients, displacement features, free volume characteristics and interaction energies of the gas molecules were investigated. In particular, the diffusion micro-mechanism of the gas molecules was observed. The differences in diffusion features among the gas molecules were discussed, and the factors influencing the diffusion of the gas molecules were compared. Simulation results indicate that the diffusion coefficients of gas molecules in cellulose is an order of magnitude lower than that in oil, and the diffusion coefficients of these gas molecules in the two types of insulation media have different orders. Free volume of gas molecules is the main factor that influences the diffusion behaviour in oil, whereas intermolecular interaction is the main influencing factor of diffusion behaviour in cellulose.     

Oxidation of phosphatidylserine: a mechanism for plasma membrane phospholipid scrambling during apoptosis?

Selective oxidation of phosphatidylserine (PS) during apoptosis precedes its externalization in plasma membrane and is essential for the engulfment of apoptotic cells. To experimentally test whether PS oxidation stimulates its externalization via its effects on aminophospholipid translocase (APT) or by enhanced PS scrambling, action of oxidized PS (PSox) was studied using leukemia HL-60 cells and lymphoma Raji cells. Both PS and PSox were equally well recognized by APT. PSox did not inhibit APT. Rate of transmembrane PS diffusion was fourfold higher in cells with integrated PSox than with PS. Thus, PSox acts as a "non-enzymatic scramblase" likely contributing to PS externalization.     

Molecular dynamics simulation approach for the prediction of transmembrane helix–helix heterodimers assembly

Computational methods are useful to identify favorable structures of transmembrane (TM) helix oligomers when experimental data are not available or when they cannot help to interpret helix-helix association. We report here a global search method using molecular dynamics (MD) simulations to predict the structures of transmembrane homo and heterodimers. The present approach is based only on sequence information without any experimental data and is first applied to glycophorin A to validate the protocol and to the HER2-HER3 heterodimer receptor. The method successfully reproduces the experimental structures of the TM domain of glycophorin A (GpA(TM)) with a root mean square deviation of 1.5 A. The search protocol identifies three energetically stable models of the TM domain of HER2-HER3 receptor with favorable helix-helix arrangement, including right-handed and left-handed coiled-coils. The predicted TM structures exhibit the GxxxG-like motif at the dimer interface which is presumed to drive receptor oligomerization. We demonstrate that native structures of TM domain can be predicted without quantitative experimental data. This search protocol could help to predict structures of the TM domain of HER heterodimer family.     

Molecular dynamics simulation of the reciprocal fused LiF–KBr mixture: local structure and self-diffusion coefficients

This paper describes the molecular dynamics simulation of the reciprocal fused LiF–KBr mixture, which is located above the critical mixing point, in the temperature range 1280–1450 K. The first coordination sphere is found to form as follows: a smaller ion is formed around a smaller counter-ion, and a larger ion is formed around a larger counter-ion. The calculated concentration dependence of the self-diffusion coefficients and the radial distribution functions of all ion pairs indicate that the degree of association of the Li–F pair increases as the lithium fluoride fraction in the mixture decreases.     

Molecular mechanism of action of β-hairpin antimicrobial peptide arenicin: oligomeric structure in dodecylphosphocholine micelles and pore formation in planar lipid bilayers

The membrane-active, cationic, β-hairpin peptide, arenicin, isolated from marine polychaeta Arenicola marina exhibits a broad spectrum of antimicrobial activity. The peptide in aqueous solution adopts the significantly twisted β-hairpin conformation without pronounced amphipathicity. To assess the mechanism of arenicin action, the spatial structure and backbone dynamics of the peptide in membrane-mimicking media and its pore-forming activity in planar lipid bilayers were studied. The spatial structure of the asymmetric arenicin dimer stabilized by parallel association of N-terminal strands of two β-hairpins was determined using triple-resonance nuclear magnetic resonance (NMR) spectroscopy in dodecylphosphocholine (DPC) micelles. Interaction of arenicin with micelles and its oligomerization significantly decreased the right-handed twist of the β-hairpin, increased its amphipathicity, and led to stabilization of the peptide backbone on a picosecond to nanosecond time scale. Relaxation enhancement induced by water-soluble (Mn(2+)) and lipid-soluble (16-doxylstearate) paramagnetic probes pointed to the dimer transmembrane arrangement. Qualitative NMR and circular dichroism study of arenicin-2 in mixed DPC/1,2-dioleoyl-sn-glycero-3-phosphoglycerol bicelles, sodium dodecyl sulfate micelles, and lipid vesicles confirmed that a similar dimeric assembly of the peptide was retained in membrane-mimicking systems containing negatively charged lipids and detergents. Arenicin-induced conductance was dependent on the lipid composition of the membrane. Arenicin low-conductivity pores were detected in the phosphatidylethanolamine-containing lipid mixture, whereas the high-conductivity pores were observed in an exclusively anionic lipid system. The measured conductivity levels agreed with the model in which arenicin antimicrobial activity was mediated by the formation of toroidal pores assembled of two, three, or four β-structural peptide dimers and lipid molecules. The structural transitions involved in arenicin membrane-disruptive action are discussed.     

Molecular basis for the glycosphingolipid-binding specificity of α-synuclein: key role of tyrosine 39 in membrane insertion

Cell surface glycosphingolipids (GSLs) including gangliosides play a key role in the regulation of the conformation, oligomerization, and fibrillation of amyloidogenic proteins. Correspondingly, most amyloidogenic proteins possess a functional GSL-binding motif (GBM). Sequence alignments of GSL-binding proteins against the GBM of α-synuclein allowed the establishment of a consensus GBM sequence defined as K/H/R/-X(1-4)-Y/F-X(4-5)-K/H/R, where at least one of the X(1-4) residues is glycine. The GBMs of α-synuclein (34-KEGVLYVGSKTK-45) and Alzheimer's disease β-amyloid peptide (Aβ) (5-RHDSGYEVHHQK-16) consist of a structurally related loop centered on tyrosine (Y39 for α-synuclein, Y10 for Aβ). Surface pressure measurements of GSL monolayers at the air-water interface allowed us to determine the following order for α-synuclein-GSL interactions: GM3 > Gb3 > GalCer-NFA > GM1 > sulfatide > GalCer-HFA > LacCer > GM4 > GM2 > asialo-GM1 > GD3, indicating a marked preference for GSLs with one, three, or five sugar units. The insertion of α-synuclein into sphingomyelin-containing monolayers was strongly stimulated by the presence of GM3. This effect was not observed with phosphatidylcholine monolayers, suggesting that the ganglioside facilitated the insertion of α-synuclein into raft-like membrane domains. Molecular dynamics simulations suggested that the side chain of Y39 was deeply inserted between GM3 head groups. Monolayer experiments with mutant GBM peptides showed that Y39, K34, and K45 were important for GM3 binding, whereas only Y39 appeared critical for GM1 recognition. The interaction of Aβ 5-16 with GM1 involved R5, H13, H14, and K16, but not Y10. These data indicate that subtle amino acid variations in the consensus GBM of α-synuclein and Aβ conferred distinct GSL-binding properties.     

The role of membrane thickness in charged protein–lipid interactions

Charged amino acids are known to be important in controlling the actions of integral and peripheral membrane proteins and cell disrupting peptides. Atomistic molecular dynamics studies have shed much light on the mechanisms of membrane binding and translocation of charged protein groups, yet the impact of the full diversity of membrane physico-chemical properties and topologies has yet to be explored. Here we have performed a systematic study of an arginine (Arg) side chain analog moving across saturated phosphatidylcholine (PC) bilayers of variable hydrocarbon tail length from 10 to 18 carbons. For all bilayers we observe similar ion-induced defects, where Arg draws water molecules and lipid head groups into the bilayers to avoid large dehydration energy costs. The free energy profiles all exhibit sharp climbs with increasing penetration into the hydrocarbon core, with predictable shifts between bilayers of different thickness, leading to barrier reduction from 26 kcal/mol for 18 carbons to 6 kcal/mol for 10 carbons. For lipids of 10 and 12 carbons we observe narrow transmembrane pores and corresponding plateaus in the free energy profiles. Allowing for movements of the protein and side chain snorkeling, we argue that the energetic cost for burying Arg inside a thin bilayer will be small, consistent with recent experiments, also leading to a dramatic reduction in pK_a shifts for Arg. We provide evidence that Arg translocation occurs via an ion-induced defect mechanism, except in thick bilayers (of at least 18 carbons) where solubility-diffusion becomes energetically favored. Our findings shed light on the mechanisms of ion movement through membranes of varying composition, with implications for a range of charged protein–lipid interactions and the actions of cell-perturbing peptides. This article is part of a Special Issue entitled: Membrane protein structure and function.     

Molecular dynamics studies of simple membrane — Water interfaces: Structure and functions in the beginnings of cellular life

Molecular dynamics computer simulations of the structure and functions of a simple membrane are performed in order to examine whether membranes provide an environment capable of promoting protobiological evolution. Our model membrane is composed of glycerol 1-monooleate. It is found that the bilayer surface fluctuates in time and space, occasionally creating thinning defects in the membrane. These defects are essential for passive transport of simple ions across membranes because they reduce the Bom barrier to this process by approximately 40%. Negative ions are transferred across the bilayer more readily than positive ions due to favorable interactions with the electric field at the membrane-water interface. Passive transport of neutral molecules is, in general, more complex than predicted by the solubility-diffusion model. In particular, molecules which exhibit sufficient hydrophilicity and lipophilicity concentrate near membrane surfaces and experience interfacial resistance to transport. The membrane-water interface forms an environment suitable for heterogeneous catalysis. Several possible mechanisms leading to an increase of reaction rates at the interface are discussed. We conclude that vesicles have many properties that make them very good candidates for earliest protocells. Some potentially fruitful directions of experimental and theoretical research on this subject are proposed.     

IB–LBM simulation of the haemocyte dynamics in a stenotic capillary

To study the behaviour of a haemocyte when crossing a stenotic capillary, the immersed boundary–lattice Boltzmann method was used to establish a quantitative analysis model. The haemocyte was assumed to be spherical and to have an elastic cell membrane, which can be driven by blood flow to adopt a highly deformable character. In the stenotic capillary, the spherical blood cell was stressed both by the flow and the wall dimension, and the cell shape was forced to be stretched to cross the stenosis. Our simulation investigated the haemocyte crossing process in detail. The velocity and pressure were anatomised to obtain information on how blood flows through a capillary and to estimate the degree of cell damage caused by excessive pressure. Quantitative velocity analysis results demonstrated that a large haemocyte crossing a small stenosis would have a noticeable effect on blood flow, while quantitative pressure distribution analysis results indicated that the crossing process would produce a special pressure distribution in the cell interior and to some extent a sudden change between the cell interior and the surrounding plasma.