Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Molecular Dynamic simulation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37Linked to original sources
PubMed2006

Perfect wetting along a three-phase line: theory and molecular dynamics simulations.

Wetting behavior along a three-phase equilibrium has been obtained by density gradient theory (DGT) and molecular dynamics simulations for a type-II equal size Lennard-Jones mixture. In order to perform a consistent comparison between both methodologies, the molecular parameters of this type of mixture were defined from the global phase diagram of equal size Lennard-Jones mixtures. We have found excellent agreement between predictions from the DGT (coupled to a Lennard-Jones equation for the bulk phases) and simulations results for both the phase and interface behavior, in the whole temperature, pressure, and concentration ranges. For all conditions explored in this work, this type-II mixture shows a three-phase equilibrium composed by a bulk immiscible liquid phase (L1) and a bulk gas phase (G) separated by a second immiscible liquid phase (L2). A similar phase distribution is obtained from the interfacial concentration profile in the whole range of conditions used in this work. This type of structure is a clear evidence that L2 completely wets the GL1 interface. The wetting behavior is also confirmed by the values and evolution of the interfacial tensions. In summary, this kind of type-II mixture does not show wetting transitions and exhibits a permanent perfect wetting in all the thermodynamic conditions explored here.

Journal Article↗
PubMed2004

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.

Carbohydrate Conformation↗
PubMed2003

Binding interactions between the core central domain of 16S rRNA and the ribosomal protein S15 determined by molecular dynamics simulations.

The goal of the current study is to utilize molecular dynamic (MD) simulations to investigate the dynamic behavior of 16S rRNA in the presence and absence of S15 and to identify the binding interactions between these two molecules. The simulations show that: (i) 16S rRNA remains in a highly folded structure when it is bound to S15; (ii) in the absence of S15, 16S rRNA significantly alters its conformation and transiently forms conformations that are similar to the bound structure that make it available for binding with S15; (iii) the unbound rRNA spends the majority of its time in extended conformations. The formation of the extended conformations is a result of the molecule reaching a lower electrostatic energy and the formation of the highly folded, crystal-like conformation is a result of achieving a lower solvation energy. In addition, our MD simulations show that 16S rRNA and S15 bind across the major groove of helix 22 (H22) via electrostatic interactions. The negatively charged phosphate groups of G658, U740, G741 and G742 bind to the positively charged S15 residues Lys7, Arg34 and Arg37. The current study provides a dynamic view of the binding of 16S rRNA with S15.

Base Sequence↗
PubMed1995

Molecular dynamics simulations of chlorambucil/DNA adducts. A structural basis for the 5'-GNC interstrand DNA crosslink formed by nitrogen mustards.

The alkylation of DNA by chlorambucil has been studied using a computational approach. Molecular dynamics simulations were performed on the fully solvated non-covalent complex, two monoadducts and a crosslinked diadduct of chlorambucil with the d(CGG3G2CGC).-d(GCG1CCCG) duplex, in which the N7 atoms of G1, G2 and G3 are potential alkylation sites. The results provide a structural basis for the preference of nitrogen mustards to crosslink DNA duplexes at a 5'-GNC site (a 1,3 crosslink, G1-G3) rather than at a 5'-GC sites (a 1,2 crosslink, G1-G2). In the non-covalent complex simulation the drug reoriented from a non-interstrand crosslinking location to a position favorable for G1-G3 diadduct formation. It proved possible to construct a G1-G3 diadduct from a structure from the non-covalent simulation, and continue the molecular dynamics calculation without further disruption of the DNA structure. A crosslinked diadduct developed with four BII conformations on the 3' side of each alkylated guanine and of their respective complementary cytosine. In the first monoadduct simulation the starting point was the same DNA conformation used in the crosslinked diadduct simulation with alkylation at G1. In this simulation the DNA deformation was reduced, with the helix returning to a more canonical form. A second monoadduct simulation was started from a canonical DNA conformation alkylated at G3. Here, no significant motion towards a potential crosslinking conformation occurred. Collectively, the results suggest that crosslink formation is dependent upon the drug orientation prior to alkylation and the required deformation of the DNA to permit 1,3 crosslinking can largely be achieved in the non-covalent complex.

Alkylation↗
PubMed2002

Large-scale molecular dynamics simulations of general anesthetic effects on the ion channel in the fully hydrated membrane: the implication of molecular mechanisms of general anesthesia.

Interactions of volatile anesthetics with the central nervous system are characterized by low yet specific binding affinities. Although neurotransmitter-gated ion channels are considered the primary anesthetic targets, the mechanism of action at the molecular level remains elusive. We consider here the theoretical implications of channel dynamics on anesthetic action in a simplified membrane-channel system. Large-scale 2.2-ns all-atom molecular dynamics simulations were performed to study the effects of halothane, a clinical anesthetic, on a gramicidin A (gA) channel in a fully hydrated dimyristoyl phosphatidylcholine membrane. In agreement with experimental results, anesthetics preferentially target the anchoring residues at the channel-lipid-water interface. Although the anesthetic effect on channel structure is minimal, the presence of halothane profoundly affects channel dynamics. For 2.2-ns simulation, the rms fluctuation of gA backbone in the lipid core increases from approximately equal 1 A in the absence of anesthetics to approximately equal 1.5 A in the presence of halothane. Autocorrelation analysis reveals that halothane (i) has no effect on the subpicosecond librational motion, (ii) prolongs the backbone autocorrelation time in the 10- to 100-ps time scale, and (iii) significantly decreases the asymptotic values of generalized order parameter and correlation time of nanosecond motions for the inner but not the outer residues. The simulation results discount the viewpoint of a structure-function paradigm that overrates the importance of structural fitting between general anesthetics and yet-unidentified hydrophobic protein pockets. Instead, the results underscore the global, as opposed to local, effects of anesthetics on protein dynamics as the underlying mechanisms for the action of general anesthetics and possibly of other low-affinity drugs.

Anesthesia, General↗
PubMed1986

Molecular dynamics simulations of native and substrate-bound lysozyme. A study of the average structures and atomic fluctuations.

Molecular dynamics simulations of hen egg-white lysozyme in the free and substrate-bound states are reported and the nature of the average structures and atomic fluctuations are analyzed. Crystallographic water molecules of structural importance, as determined by hydrogen-bonding, were included in the simulations. Comparisons are made between the dynamics and the X-ray results for the atomic positions, the main-chain and side-chain dihedral angles, and the hydrogen-bonding geometry. Improvements over earlier simulations in the potential energy function and methodology resulted in stable trajectories with the C alpha co-ordinates within 1.5 A of the starting X-ray structure. Structural features analyzed in the simulations agreed well with the X-ray results except for some surface residues. The Asx chi 2 dihedral distribution and the geometry of hydrogen bonding at reverse turns show differences; possible causes are discussed. The relation between the magnitudes and time-scales of the residue fluctuations and secondary structural features, such as helices beta-sheets and coiled loops, is examined. Significant differences in the residue mobilities between the simulations of the free and substrate-bound states were found in a region of the enzyme that is in direct contact with the substrate and in a region that is distant from the active-site cleft. The dynamic behavior of the structural water molecules is analyzed by examining the correlation between the fluctuations of the water oxygens and the lysozyme heavy-atoms to which they are hydrogen-bonded.

Animals↗
PubMed2005

Effect of nanostructure on the properties of water at the water-hydrophobic interface: a molecular dynamics simulation.

The local structure of water near hydrophobic surfaces of different surface topographies has been analyzed by molecular dynamics simulation. An alkane crystal has been taken as the parent model for a hydrophobic surface. Surface structures were created by placing pits into it, which were half a nanometer deep and several nanometers wide. Around all structures, the water has a lower density, less orientational ordering, fewer water-water hydrogen bonds, and fewer surface contacts than for a flat unstructured surface. This indicates that the structured surfaces are more hydrophobic than the flat surface. Of the structures investigated, pits with a diameter of approximately 2.5 nm were effective in increasing the hydrophobic character of the surface.

Journal Article↗
PubMed2007

Molecular dynamics simulation study on controlling the adsorption behavior of polyethylene by fine tuning the surface nanodecoration of graphite.

Molecular dynamics simulations are applied to study the adsorption of polyethylene with different chain lengths on patterned graphite surfaces that contain nanoscale protrusions. The influence of the nanostructure on the strong attractive interaction inherently in the hydrophobic polyethylene and hydrophobic graphite system is investigated by modifying the top surface area and the height and the shape of the protrusions. The results are analyzed in terms of the chain configuration, the adsorption energy, the global orientational order parameter, and the normalized surface-chain contacting pair number in the first adsorption layer. When the size of the protrusion increases, the adsorption energy, the order parameter, and the normalized surface-chain contacting pair number decrease at a fixed chain length. When the size of the protrusion is fixed, the average adsorption energy per monomer and the order parameter decrease with increasing chain length because of the stronger intramolecular interactions between the monomers. Changing the protrusion shape in a suitable way will effectively reduce the strong surface-chain interaction.

Journal Article↗
PubMed2006

Molecular dynamics simulations of hydration, dissolution and nucleation processes at the alpha-quartz (0001) surface in liquid water.

Molecular dynamics simulations have been employed to investigate the hydration and dissolution of alpha-quartz (0001) surfaces in a liquid water environment. Our study indicates that the structure of the water layers near the surfaces is affected by the nature of the substrate surface and by temperature. Ordered mono-layers of interfacial water molecules form in the region of the substrate where the surface is highly charged and built up of Si-O-Si bridges. As the temperature is increased this ordered mono-layer structure is gradually lost. When the surface is terminated by silanol groups, the water retains liquid-like properties even at low temperature and the molecules are distributed in a random manner, without the formation of distinct ordered mono-layers of water molecules near the surface. Taking into account the entropy of the system, the calculated energies of stepwise dissolution of a silicon species from the surface suggest that on thermodynamic grounds the complete dissolution of silicon atoms from the quartz surfaces in a liquid water environment is an endothermic process, but that the formation of a -Si(OH)3 species at the surface would be possible. In addition, if the Si(OH)(4) species were to be dissolved, it would remain near the surface, and re-deposition at the defect-free surface is thermodynamically preferred, although there is an activation enthalpy to the first step in the process of nucleation of Si(OH)4 at the perfect surface.

Journal Article↗
PubMed2004

Short-time Fourier transform analysis of ab initio molecular dynamics simulation: collision reaction between NH+4 (NH3)2 and NH3.

An analyzing technique of the ab initio molecular dynamics simulation is proposed with the use of short-time Fourier transform (ST-FT). The ST-FT analysis demonstrates the dynamical change of the vibrational states in the simulated system. Numerical assessments are preformed for the collision reaction of the ammonia cluster ion NH+(4) (NH3)(2) with the ammonia monomer NH3. Spectrogram obtained by the ST-FT method, which corresponds to the time evolution of vibrational power spectra, clarifies the relationship between the vibrational states and the reaction channels such as nonreactive collision, substitution, and incorporation.

Journal Article↗
PubMed2005

Molecular dynamics simulations of side chain liquid crystal polymer molecules in isotropic and liquid-crystalline melts.

A detailed molecular dynamics simulation study is described for a polysiloxane side chain liquid crystal polymer (SCLCP). The simulations use a coarse-grained model composed of a combination of isotropic and anisotropic interaction sites. On cooling from a fully isotropic polymer melt, we see spontaneous microphase separation into polymer-rich and mesogen-rich regions. Upon application of a small aligning potential during cooling, the structures that form on microphase separation anneal to produce a smectic-A phase in which the polymer backbone is largely confined between the smectic layers. Several independent quenches from the melt are described that vary in the strength of the aligning potential and the degree of cooling. In each quench, defects were found where the backbone chains hop from one backbone-rich region to the next by tunneling through the mesogenic layers. As expected, the number of such defects is found to depend strongly on the rate of cooling. In the vicinity of such a defect, the smectic-A structure of the mesogen-rich layers is disrupted to give nematiclike ordering. Additionally, several extensive annealing runs of approximately 40 ns duration have been carried out at the point of microphase separation. During annealing the polymer backbone is seen to be slowly excluded from the mesogenic layers and lie perpendicular to the smectic-A director. These observations agree with previous assumptions about the structure of a SCLCP and with interpretations of x-ray diffraction and small angle neutron scattering data. The flexible alkyl spacers, which link the backbone to the mesogens, are found to form sublayers around the backbone layer.

Journal Article↗
PubMed2002

The conformation of dehydroalanine in short homopeptides: molecular dynamics simulations of a 6-residue chain.

A molecular dynamics study about the conformational preferences in a chloroform solution of a homo-oligomer constituted by six residues of dehydroalanine is presented. For this purpose, two sets of force-field parameters and explicit solvent molecules have been used. Furthermore, ab initio calculations have been performed in order to estimate 1[H]-NMR chemical shifts. Results have been compared with experimental data.

Alanine↗
PubMed2004

Molecular dynamics simulations of discoidal bilayers assembled from truncated human lipoproteins.

Human apolipoprotein A-1 (apo A-1) is the major protein component of high-density lipoproteins. The apo A-1 lipid-binding domain was used as a template for the synthesis of amphipathic helical proteins termed membrane scaffold proteins, employed to self-assemble soluble monodisperse discoidal particles called Nanodiscs. In these particles, membrane scaffold proteins surround a lipid bilayer in a belt-like fashion forming bilayer disks of discrete size and composition. Here we investigate the structure of Nanodiscs through molecular dynamics simulations in which Nanodiscs were built from scaffold proteins of various lengths. The simulations showed planar or deformed Nanodiscs depending on optimal length and alignment of the scaffold proteins. Based on mean surface area per lipid calculations, comparison of small-angle x-ray scattering curves, and the relatively planar shape of Nanodiscs made from truncated scaffold proteins, one can conclude that the first 17 to 18 residues of the 200-residue apo A-1 lipid-binding domain are not involved in formation of the protein "belts" surrounding the lipid bilayer. To determine whether the addition of an integral membrane protein has an effect on the overall structure of a Nanodisc, bacteriorhodopsin was embedded into a Nanodisc and simulated using molecular dynamics, revealing a planar disk with a slightly rectangular shape.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed2006

Interaction between K+ channel gate modifier hanatoxin and lipid bilayer membranes analyzed by molecular dynamics simulation.

Hanatoxin (HaTx) is an ellipsoidal-shaped peptide that binds to the voltage sensor of voltage-dependent channels. Of physicochemical interest, HaTx has a "ring" of charged residues around its periphery and a hydrophobic protrusion. It has previously been postulated that HaTx binds to and functions on the surface of membranes, but a recent fluorescent-quenching study has implied a fairly deep positioning of HaTx in the lipid bilayer membrane. We carried out numerous molecular dynamic simulations of HaTx1, a well-studied variant of HaTx, in fully hydrated phospholipid bilayers. The system reproduced the surface-binding mode of HaTx1, in which HaTx1 resided in the extracellular side (outer) of the water/membrane interface with the hydrophobic patch of HaTx1 facing the membrane interior. On the other hand, analyses with various parameter settings suggested that the surface-binding mode was unstable because of the substantial attractive electrostatic force between HaTx1 and the lipid head groups of the inner (opposite) leaflet. Compared with this electrostatic force, the energetic cost for membrane deformation involving meniscus formation appeared to be small. In an attempt to interpret the quenching data, we consider the possibility of dimpling (meniscus formation) that brings HaTx1 inward (only ~0.7-0.8 nm above the bilayer center), while accounting for the flexibility of both leaflets of the membrane and the long-range interaction between positively charged residues of the membrane-bound peptide and the polar head groups of the opposite leaflet of the membrane. It is suggested that molecular dynamics simulations taking into account the flexibility of the membrane surface is potentially useful in interpreting the fluorescence-quenching data.

Computer Simulation↗
PubMed2005

Molecular dynamics simulations of peptide-surface interactions.

Proteins, which are bioactive molecules, adsorb on implants placed in the body through complex and poorly understood mechanisms and directly influence biocompatibility. Molecular dynamics modeling using empirical force fields provides one of the most direct methods of theoretically analyzing the behavior of complex molecular systems and is well-suited for the simulation of protein adsorption behavior. To accurately simulate protein adsorption behavior, a force field must correctly represent the thermodynamic driving forces that govern peptide residue-surface interactions. However, since existing force fields were developed without specific consideration of protein-surface interactions, they may not accurately represent this type of molecular behavior. To address this concern, we developed a host-guest peptide adsorption model in the form of a G(4)-X-G(4) peptide (G is glycine, X is a variable residue) to enable determination of the contributions to adsorption free energy of different X residues when adsorbed to functionalized Au-alkanethiol self-assembled monolayers (SAMs). We have previously reported experimental results using surface plasmon resonance (SPR) spectroscopy to measure the free energy of peptide adsorption for this peptide model with X = G and K (lysine) on OH and COOH functionalized SAMs. The objectives of the present research were the development and assessment of methods to calculate adsorption free energy using molecular dynamics simulations with the GROMACS force field for these same peptide adsorption systems, with an oligoethylene oxide (OEG) functionalized SAM surface also being considered. By comparing simulation results to the experimental results, the accuracy of the selected force field to represent the behavior of these molecular systems can be evaluated. From our simulations, the G(4)-G-G(4) and G(4)-K-G(4) peptides showed minimal to no adsorption to the OH SAM surfaces and the G(4)-K-G(4) showed strong adsorption to the COOH SAM surface, which is in agreement with our SPR experiments. Contrary to our experimental results, however, the simulations predicted a relatively strong adsorption of G(4)-G-G(4) peptide to the COOH SAM surface. In addition, both peptides were unexpectedly predicted to adsorb to the OEG surface. These findings demonstrate the need for GROMACS force field parameters to be rebalanced for the simulation of peptide adsorption behavior on SAM surfaces. The developed methods provide a direct means of assessing, modifying, and validating force field performance for the simulation of peptide and protein adsorption to surfaces, without which little confidence can be placed in the simulation results that are generated with these types of systems.

Adsorption↗
PubMed1994

Molecular dynamics simulation of a phospholipid membrane.

We present the results of molecular dynamics (MD) simulations of a phospholipid membrane in water, including full atomic detail. The goal of the simulations was twofold: first we wanted to set up a simulation system which is able to reproduce experimental results and can serve as a model membrane in future simulations. This goal being reached it is then further possible to gain insight in to those properties that are experimentally more difficult to access. The system studied is dipalmitoylphosphatidylcholine/water, consisting of 5408 atoms. Using original force field parameters the membrane turned out to approach a gel-like state. With slight changes of the parameters, the system adopted a liquid-crystalline state. Separate 80 ps runs were performed on both the gel and liquid-crystalline systems. Comparison of MD results with reliable experimental data (bilayer repeat distance, surface area per lipid, tail order parameters, atom distributions) showed that our simulations, especially the one in the liquid-crystalline phase, can serve as a realistic model for a phospholipid membrane. Further analysis of the trajectories revealed valuable information on various properties. In the liquid-crystalline phase, the interface turns out to be quite diffuse, with water molecules penetrating into the bilayer to the position of the carbonyl groups. The 10-90% width of the interface turns out to be 1.3 nm and the width of the hydrocarbon interior 3.0 nm. The headgroup dipoles are oriented at a small angle with respect to the bilayer plane. The resulting charge distribution is almost completely cancelled by the water molecules. The electron density distribution shows a large dip in the middle of the membrane. In this part the tails are more flexible. The mean life time between dihedral transitions is 20 ps. The average number of gauche angles per tail is 3.5. The occurrence of kinks is not a significant feature.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed2005

Molecular dynamics simulation of cocaine binding with human butyrylcholinesterase and its mutants.

Molecular dynamics (MD) simulations were carried out to study cocaine binding with wild-type human butyrylcholinesterase (BChE) and its mutants based on a recently reported X-ray crystal structure of human BChE. For each BChE-cocaine system, we simulated both the nonprereactive and prereactive complexes in water. Despite the significant difference found at the acyl binding pocket, the simulated structures confirm the fundamental structural and mechanistic insights obtained from earlier computational studies of wild-type BChE with cocaine based on a homology model, e.g. the rate-determining step for BChE-catalyzed hydrolysis of biologically active (-)-cocaine is the (-)-cocaine rotation in the active site from the nonprereactive BChE-(-)-cocaine complex to the prereactive complex. It has been demonstrated that the MD simulations on both the nonprereactive and prereactive BChE-cocaine complexes can clearly reveal whether specific mutations produce the desired BChE-(-)-cocaine binding structures in which the (-)-cocaine rotation is less hindered while the required prereactive BChE-(-)-cocaine binding is maintained. Based on the MD simulations, both A328W/Y332A and A328W/Y332G BChE's are expected to have catalytic activity for (-)-cocaine hydrolysis higher than that of wild-type BChE and the activity of A328W/Y332G BChE should be slightly higher than that of A328W/Y332A BChE due to the less-hindered (-)-cocaine rotation in the mutant BChE's. However, the less-hindered (-)-cocaine rotation is only a necessary condition for a higher activity mutant BChE. The (-)-cocaine rotation is also less hindered in A328W/Y332A/Y419S BChE, but (-)-cocaine binds with A328W/Y332A/Y419S BChE in a way that is not suitable for the catalysis. Thus, A328W/Y332A/Y419S BChE is expected to lose the catalytic activity. The computational predictions were confirmed by our experimental kinetic data, demonstrating that the MD simulation-based computational protocol used in this study is reliable in prediction of the catalytic activity of BChE mutants for (-)-cocaine hydrolysis.

Butyrylcholinesterase↗
PubMed2006

Vibrational dynamics of DNA. III. Molecular dynamics simulations of DNA in water and theoretical calculations of the two-dimensional vibrational spectra.

A theoretical description of the vibrational excitons in DNA is presented by using the vibrational basis mode theory developed in Papers I and II. The parameters obtained from the density functional theory calculations, such as vibrational coupling constants and basis mode frequencies, are used to numerically simulate two-dimensional (2D) IR spectra of dG(n):dC(n) and dA(n):dT(n) double helices with n varying from 1 to 10. From the molecular dynamics simulations of dG(5)C(5) and dA(5)T(5) double helices in D(2)O solution, it is found that the thermally driven internal motions of these systems in an aqueous solution do not induce strong fluctuations of basis mode frequencies nor vibrational couplings. In order to construct the two-exciton Hamiltonian, the vibrational anharmonicities of eight basis modes are obtained by carrying out B3LYP6-31G(*) calculations for the nine basis modes. The simulated 2D IR spectra of dG(n):dC(n) double helix in D(2)O solution are directly compared with closely related experimental results. The 2D IR spectra of dG(n):dC(n) and dA(n):dT(n) are found to be weakly dependent on the number of base pairs. The present work demonstrates that the computational procedure combining quantum chemistry calculation and molecular dynamics simulation methods can be of use to predict 2D IR spectra of nucleic acids in solutions.

Base Pairing↗