Molecular dynamics simulation of oleyl oleate swollen micelles system

Problems with transdermal drug delivery were directly associated with the skin barrier which is the lipid bilayer at the stratum corneum. Chemical penetration enhancers such as swollen micelles that formed from the solubilisation of the surfactants in the nano-emulsion system could provide an effective solution. However, the structural properties of swollen micelles from nano-emulsions of palm-oil esters, whose behaviour is related to colloidal systems, have not been studied in great detail. In this paper, we report on the use of molecular dynamics (MD) simulations to investigate the structural properties of swollen micelles of oleyl oleate (OE). Five series of 10 ns MD simulations were performed at different micelle compositions to determine the structural evolution of OE/Span20 (S20) swollen micelles. We also carried out four MD simulations on the structure of S20, OE/S20, Tween80 (T80) and OE/T80 micelles to study the effect of different surfactants and the addition of OE into the systems. The shapes of the swollen micelles were observed to vary by the difference in the micelle composition, the surfactants used and the addition of OE. The results were correlated with published theory, and consistent with experimental results on the phase behaviour of the nano-emulsion system.     

Molecular dynamics simulation for shape change of water-in-oil droplets

We performed molecular dynamics (MD) simulations of water-in-oil droplet shape transformations induced by the addition of polymer chains. In a prior experiment, transformations of spherical droplets to rod-like, worm-like and network-like droplets were observed. In our previous study, we reproduced rod-like droplets via coarse-grained MD simulations, and the mechanism for the droplet shape change was elucidated by considering the contact area between the chains and the surfactant head groups. However, in that simulation model, we could not reproduce the worm-like and network-like droplets. In this study, we improved the simulation model. For a small number of chains, several spherical droplets were obtained. As the number of chains increased, the spherical droplets were transformed to rod-like, worm-like and network-like shapes by coalescence of the droplets. The calculated and experimental results agreed well, and we verified that the mechanism for the droplet shape transformations observed in the present simulations could be explained by the mechanism suggested in the previous study.     

Molecular dynamics simulations from putative transition states of alpha-spectrin SH3 domain

A series of molecular dynamics simulations in explicit solvent were started from nine structural models of the transition state of the SH3 domain of alpha-spectrin, which were generated by Lindorff-Larsen et al. (Nat Struct Mol Biol 2004;11:443-449) using molecular dynamics simulations in which experimental Phi - values were incorporated as restraints. Two of the nine models were simulated 10 times for 200 ns and the remaining models simulated two times for 200 ns. Complete folding was observed in one case, while in the other simulations partial folding or unfolding events were observed, which were characterized by a regularization of elements of secondary structure. These results are consistent with recent experimental evidence that the folding of SH3 domains involves low populated intermediate states.     

Molecular dynamics simulation of content of bound water in ethylene glycol and glycerol solution

In this paper, the content of bound water was studied to evaluate the cryoprotective properties of ethylene glycol and glycerol solution. Molecular dynamic models for the solution were built, the classification principle and statistical methods of water molecules in solutions were presented, respectively. The content of bound water with various hydroxyl molarity at different temperatures was obtained through molecular dynamic simulation. The results reveal that the content of bound water increases with increasing hydroxyl molarity, but decreases with increasing temperature. It was found that, the content of bound water in ethylene glycol solution is always slightly more than that in glycerol solution, regardless of whether the temperature increases or not.     

Molecular dynamics simulation of intrinsically disordered proteins

Intrinsically disordered proteins are biomolecules that do not have a definite 3D structure; therefore, their dynamical simulation cannot start from a known list of atomistic positions, such as a Protein Data Bank file. We describe a method to start a computer simulation of these proteins. The first step of the procedure is the creation of a multi-rod configuration of the molecule, derived from its primary sequence. This structure is dynamically evolved in vacuo until its gyration radius reaches the experimental average value; at this point solvent molecules, in explicit or implicit implementation, are added to the protein and a regular molecular dynamics simulation follows. We have applied this procedure to the simulation of tau, one of the largest totally disordered proteins.     

Molecular dynamics simulation of hydration in myoglobin

This study was carried out to evaluate the stability of the 89 bound water molecules that were observed in the neutron diffraction study of CO myoglobin. The myoglobin structure derived from the neutron analysis was used as the starting point in the molecular dynamics simulation using the software package CHARMM. After solvation of the protein, energy minimization and equilibration of the system, 50 ps of Newtonian dynamics was performed. This data showed that only 4 water molecules are continously bound during the length of this simulation while the other solvent molecules exhibit considerable mobility and are breaking and reforming hydrogen bonds with the protein. At any instant during the simulation, 73 of the hydration sites observed in the neutron structure are occupied by water. © 1995 Wiley-Liss, Inc.     

Wetting transition of nanodroplets of water on textured surfaces: a molecular dynamics study

The crossover behaviour of water droplet's state from the Wenzel state to the Cassie state with varying pillar height and surface fraction is examined critically using molecular dynamics. We report the effect of the system size on the wetting behaviour of water droplets by examining the contact angle for both regimes. We observe that when the droplet size is comparable to the pillar dimension, the contact angle of droplets fluctuates with increasing droplet size because of the contact line pinning, which is more pronounced in the Wenzel regime. We further demonstrate the phantom-wall method to evaluate free energy of intermediate wetting states.     

Molecular dynamics simulation of interactions in glycolytic enzymes

Two glycolytic enzymes, phosphoglycerate mutase (PGM) and enolase from Saccharomyces cerevisiae, have been chosen to detect complex formation and possible channeling, using molecular dynamics simulation. The enzymes were separated by 10 angstroms distance and placed in a water-filled box of size 173 x 173 x 173 angstroms. Three different orientations have been investigated. The two initial 3-phosphoglycerate substrate molecules near the active centers of the initial structure of PGM have been replaced with final product (2-phosphoglycerate) molecules, and 150 mM NaCl together with three Mg2+ ions have been added to the system to observe post-catalytic activity under near-physiological conditions. Analysis of interaction energies and conformation changes for 3 nsec simulation indicates that PGM and enolase do show binding affinity between their near active regions, which is necessary for channeling to occur. Interaction of the C-terminal residues Ala239 and Val240 of PGM (which partially "cap" the 2-phosphoglycerate) with enolase also favors the existence of channeling.     

Molecular dynamics simulation of self-assembled nanocubes in methanol

Molecular dynamics simulations were performed for the hexameric nanocubes of methylated (1₆) and demethylated (2₆) gear-shaped amphiphiles in pure methanol to reveal the difference in structural fluctuation between 1₆ and 2₆. Within our simulation time of 2.0 ns, the cubic structure of 1₆ in methanol is maintained, whereas that of 2₆ is collapsed. We found that the triple π-stacking moieties consisting of the three 3-pyridyl groups in 2₆ are more fluctuated than those in 1₆. This suggests that methyl groups serve to reduce structural fluctuation for nanocubes. We also found that the existence of the solvent molecules near the nanocube is an important factor for the collapse of the 2₆ structure.     

Molecular dynamics simulation of a decasaccharide fragment of heparin in aqueous solution

Molecular dynamics (MD) simulations on heparin-water-sodium systems were carried out in order to establish a simulation protocol able to represent heparin solution conformation under physiological conditions. Atomic charges suitable for heparin oligosaccharides were obtained from ab initio quantum-mechanical computations, at the 6-31G(**) level. The GROMACS forcefield, the SPC, and SPC/E water models were employed. Also heparin was simulated with IdoA residues in 1C(4) or 2S(0) conformational states. The results of the performed MD simulations are in agreement with the available experimental data, suggesting that this approach can be applied for the study of heparin interactions with its target proteins and thus play a role in the development of new antithrombotic agents.     

Molecular dynamics simulation of formation and control of nanodroplets in piezoelectric nanoejection systems

In this paper, the formation of nanodroplets in piezoelectric nanoejection processes is investigated by non-equilibrium molecular dynamics simulation. By compressing liquid propane molecules with various specific pushing periods of oscillation, the phenomena of liquid thread breakup and droplet formation are simulated. The simulation results revealed that various features aid the piezoelectric nanoejection system. Two breakup shapes including double-cone and long tail structures were found in this process. To analyse the ejection process in detail, 2D contour plots and thermal properties for various pushing periods are shown and discussed in this paper. The results show that the sizes of nanodroplets are linear depending on the pushing periods. The findings show a new control factor and mechanism for nanodroplet formation through piezoelectric nanoejection processes.     

Molecular dynamics simulation of thread break-up and formation of droplets in nanoejection system

This study investigated nanojet processes by a non-equilibrium molecular dynamics simulation. The phenomena of liquid thread break-up and droplet formation were simulated by compressing liquid propane molecules with various compressing velocities. Properties' distributions show that, at the nanoscale, density and pressure were neither uniform nor continuous during the ejection process. Shear heating phenomena were found in the contact area of the nozzle channel. A linear relationship between the length of liquid threads and the compressing velocity was also found in this study. The results from different trials using various compressing velocities show that higher compressing velocities in nanojet processes result in longer liquid thread lengths and liquid molecules with higher energy levels. Therefore, the ejection process is more unstable, resulting in an increase in the number of evaporating molecules and satellite droplets. Results that illustrate various features are presented to aid in the comprehension of the nanojet processes.     

Molecular dynamics simulation of a human thiopurine S-methyltransferase complexed with 6-mercaptopurine model

Human thiopurine S-methyltransferase (TPMT) is an essential protein in 6-mercaptopurine (6MP) drug metabolism. To understand the pharmacogenetics of TPMT and 6MP, X-ray co-crystal structures of TPMT complexes with S-adenosyl-L-methionine (AdoMet) and 6MP are required. However, the co-crystal structure of this complex has not been reported because 6MP is poorly water soluble. We used molecular dynamics (MD) simulation to predict the structure of the complex of human TPMT-AdoHcy(CH₂)6MP, where the sulfur atoms of AdoHcy and 6MP were linked by a CH₂ group. After 1300 picoseconds of MD simulation, the trajectory showed that 6MP was stabilized in the TPMT active site by formation of non-bonded interactions between 6MP and Phe40, Pro196 and Arg226 side chains of TPMT. The intersulfur distance between AdoHcy and 6MP as well as the binding modes and the interactions of our TPMT-AdoHcy model are consistent with those observed in the X-ray crystal structure of murine TPMT-AdoHcy-6MP complex. The predicted binding modes of AdoHcy and 6MP in our model are consistent with those observed in murine TPMT X-ray crystal structures, which provides structural insights into the interactions of TPMT, AdoHcy, and 6MP at the atomic level and may be used as a starting point for further study of thiopurine drug pharmacogenetics.     

Molecular dynamics simulation of a phosphatidylglycerol membrane

Although molecular dynamics simulations are an important tool for studying membrane systems, relatively few simulations have used anionic lipids. This paper reports the first simulation of a pure phosphatidylglycerol (PG) bilayer. The properties of this equilibrated palmitoyloleoylphosphatidylglycerol membrane agree with experimental observations of PG membranes and with previous simulations of monolayers and mixed bilayers containing PG lipids. These simulations also provide interesting insights into hydrogen bonding interactions in PG membranes. This equilibrated membrane will be a useful starting point for simulations of membrane proteins interacting with PG lipids.     

Molecular dynamics simulation of methane hydrate dissociation by depressurisation

Methane (CH₄) hydrate dissociation and the mechanism by depressurisation are investigated by molecular dynamics (MD) simulation. The hydrate decomposition processes are studied by the ‘vacuum removal method’ and the normal method. It is found that the hydrate decomposition is promoted by depressurisation. The quasi-liquid layer is formed in the hydrate surface layer. The driving force of dissociation is found to be controlled by the concentration gradient between the H₂O molecules of the hydrate surface layer and the H₂O molecules of the hydrate inner layer. The clathrates collapse gradually, and the hydrate decomposes layer by layer. Relative to our previous MD simulation results, this study shows that the rate of the hydrate dissociation by depressurisation is slower than that by the thermal stimulation and the inhibitor injection. This study illustrated that MD simulation can play a significant role in investigating the hydrate decomposition mechanisms.     

Molecular dynamics study on the conformational flexibility and energetics in aqueous solution of methylated beta-cyclodextrins

All possible methylated beta-cyclodextrins (CDs) with C7-symmetry have been studied by molecular dynamics simulations, in gas phase and in water solution. Energetic and structural information were obtained from the trajectory analysis. CD flexibility increases with degree of methylation, very likely due to the concomitant reduction of the intramolecular hydrogen bonds. Solvation-free energy was computed for each of the studied CDs using the MM/GBSA method. An analysis of radial distribution functions was used to determine distribution of solvent molecules around the O2, O3, and O6. The number of solvent molecules around these oxygens decreases with an increase in the degree of methylation. The DeltaS contribution from solvent thus becomes more positive when the degree of methylation increases and, consequently, the overall DeltaG in water diminishes.     

Investigating the disorder–order transition of calmodulin binding domain upon binding calmodulin using molecular dynamics simulation

We have studied the conformational transition of the calmodulin binding domains (CBD) in several calmodulin‐binding kinases, in which CBD changes from the disordered state to the ordered state when binding with calmodulin (CaM). Targeted molecular dynamics simulation was used to investigate the binding process of CaM and CBD of CaM‐dependent kinase I (CaMKI–CBD). The results show that CaMKI–CBD began to form an α‐helix and the interaction free energy between CaM and CaMKI–CBD increased once CaM fully encompassed CaMKI–CBD. Two series of CaM/CBD complex systems, including the complexes of CaM with the initially disordered and the final ordered CBD, were constructed to study the interaction using molecular dynamics simulations. Our analyses suggest that the VDW interaction plays a dominant role in CaM/CBD binding and is a key factor in the disorder–order transition of CBD. Additionally, the entropy effect is not in favor of the formation of the CaM/CBD complex, which is consistent with the experimental evidence. Based on the results, it appears that the CBD conformational change from a non‐compact extended structure to compact α‐helix is critical in gaining a favorable VDW interaction and interaction free energy. Copyright © 2009 John Wiley & Sons, Ltd.     

Extending the capabilities of targeted molecular dynamics: simulation of a large conformational transition in plasminogen activator inhibitor 1

Plasminogen activator inhibitor type 1 (PAI-1) is an inhibitor of plasminogen activators such as tissue-type plasminogen activator or urokinase-type plasminogen activator. For this molecule, different conformations are known. The inhibiting form that interacts with the proteinases is called the active form. The noninhibitory, noncleavable form is called the latent form. X-ray and modeling studies have revealed a large change in position of the reactive center loop (RCL), responsible for the interaction with the proteinases, that is inserted into a beta-sheet (s4A) in the latent form. The mechanism underlying this spontaneous conformational change (half-life = 2 h at 37 degrees C) is not known in detail. This investigation attempts to predict a transition path from the active to the latent structure at the atomic level, by using simulation techniques. Together with targeted molecular dynamics (TMD), a plausible assumption on a rigid body movement of the RCL was applied to define an initial guess for an intermediate. Different pathways were simulated, from the active to the intermediate, from the intermediate to the latent structure and vice versa under different conditions. Equilibrium simulations at different steps of the path also were performed. The results show that a continuous pathway from the active to the latent structure can be modeled. This study also shows that this approach may be applied in general to model large conformational changes in any kind of protein for which the initial and final three-dimensional structure is known.     

Molecular dynamics simulations of Factor Xa: insight into conformational transition of its binding subsites

Protein flexibility and conformational diversity is well known to be a key characteristic of the function of many proteins. Human blood coagulation proteins have multiple substrates, and various protein-protein interactions are required for the smooth functioning of the coagulation cascade to maintain blood hemostasis. To address how a protein may cope with multiple interactions with its structurally diverse substrates and the accompanied structural changes that may drive these changes, we studied human Factor X. We employed 20 ns of molecular dynamics (MD) and steered molecular dynamics (SMD) simulations on two different conformational forms of Factor X, open and closed, and observed an interchangeable conformational transition from one to another. This work also demonstrates the roles of various aromatic residues involved in aromatic-aromatic interactions, which make this dynamic transition possible.