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 1,063 records · Page 59Linked to original sources
PubMed2005

Comparison of morphology and mechanical properties of surfactant aggregates at water-silica and water-graphite interfaces from molecular dynamics simulations.

Cationic surfactants are important for a wide range of applications, including controlled drug delivery systems, emulsifiers, and chemical mechanical polishing. It is therefore important to better understand surfactant structure and properties at the solid-liquid interface. Here, classical molecular dynamics simulations with empirical potentials are used to compare the structures and mechanical properties of cationic surfactant micelles at hydrophobic (graphite) and hydrophilic (silica) surface-water interfaces. In particular, the morphology of monolayers and bilayers of C12TAB (n-dodecyltrimethylammoniumbromide) at these interfaces, and their responses to atomic force microscopy indentation, are examined. The simulations predict that surfactant monolayers and bilayers on silica evolve into a spherical micelle structure, in agreement with theoretical models of surfactant morphology. In contrast, surfactant monolayers on graphite evolve into a hemi-cylindrical structure, in agreement with experimental findings. In the simulated indentation of the micelle/silica system, the spherical micelle breaks apart and forms a surfactant monolayer. The indentation force curve has a maximum value of 2.25 nN. On the other hand, the simulated indentation of the micelle/graphite system causes the hemi-cylindrical micelle structure to break apart and the surfactant tails to wrap around the graphite indenter. The indentation force curve has a maximum value of 13 nN.

Journal Article↗
PubMed2005

Molecular dynamics simulation of titanium dioxide nanoparticle sintering.

Nanoparticles have been an area of active research in recent years due to their properties, which can be greatly different from the bulk. In this work, we study the sintering of TiO2 nanoparticles using molecular dynamics simulations. Such sintering occurs in flame reactors where nanotitania is prepared via the chloride process. Decrease in free energy due to reduction in surface area is the main driving force for sintering of particles. Simulations, at various starting temperatures and orientations, indicate that the process of sintering is strongly affected by temperature and initial orientation. Extremely high diffusion of ions in the neck region of sintering nanoparticles supports the idea that solid-state diffusion is significant in metal-oxide nanoparticle sintering. It is found that the dipole-dipole interaction between sintering nanoparticles plays a very important role at temperatures away from the melting point. The duration of the simulation is not enough to observe the complete sintering process, but important initial stages are well studied.

Journal Article↗
PubMed2006

Molecular dynamics simulations of multilayer polyelectrolyte films: effect of electrostatic and short-range interactions.

The effect of the strength of electrostatic and short-range interactions on the multilayer assembly of oppositely charged polyelectrolytes at a charged substrate was studied by molecular dynamics simulations. The multilayer buildup was achieved through sequential adsorption of charged polymers in a layer-by-layer fashion from dilute polyelectrolyte solutions. The strong electrostatic attraction between oppositely charged polyelectrolytes at each deposition step is a driving force behind the multilayer growth. Our simulations have shown that a charge reversal after each deposition step is critical for steady multilayer growth and that there is a linear increase in polymer surface coverage after the first few deposition steps. Furthermore, there is substantial intermixing between chains adsorbed during different deposition steps. We show that the polymer surface coverage and multilayer structure are each strongly influenced by the strength of electrostatic and short-range interactions.

Journal Article↗
PubMed2004

Heat capacity effects associated with the hydrophobic hydration and interaction of simple solutes: a detailed structural and energetical analysis based on molecular dynamics simulations.

We examine the SPCE [H. J. C. Berendsen et al., J. Chem. Phys. 91, 6269 (1987)] and TIP5P [M. W. Mahoney and W. L. Jorgensen, J. Chem. Phys 112, 8910 (2000)] water models using a temperature series of molecular-dynamics simulations in order to study heat-capacity effects associated with the hydrophobic hydration and interaction of xenon particles. The temperature interval between 275 and 375 K along the 0.1-MPa isobar is studied. For all investigated models and state points we calculate the excess chemical potential for xenon employing the Widom particle insertion technique. The solvation enthalpy and excess heat capacity is obtained from the temperature dependence of the chemical potentials and, alternatively, directly by Ewald summation, as well as a reaction field based method. All three methods provide consistent results. In addition, the reaction field technique allows a separation of the solvation enthalpy into solute/solvent and solvent/solvent parts. We find that the solvent/solvent contribution to the excess heat capacity is dominating, being about one order of magnitude larger than the solute/solvent part. This observation is attributed to the enlarged heat capacity of the water molecules in the hydration shell. A detailed spatial analysis of the heat capacity of the water molecules around a pair of xenon particles at different separations reveals that even more enhanced heat capacity of the water located in the bisector plane between two adjacent xenon atoms is responsible for the maximum of the heat capacity found for the desolvation barrier distance, recently reported by Shimizu and Chan [J. Am. Chem. Soc. 123, 2083 (2001)]. The about 60% enlarged heat capacity of water in the concave part of the joint xenon-xenon hydration shell is the result of a counterplay of strengthened hydrogen bonds and an enhanced breaking of hydrogen bonds with increasing temperature. Differences between the two models with respect to the heat capacity in the xenon-xenon contact state are attributed to the different water model bulk heat capacities, and to the different spatial extension of the structure effect introduced by the hydrophobic particles. Similarities between the different states of water in the joint xenon-xenon hydration shell and the properties of stretched water are discussed.

Journal Article↗
PubMed1995

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 continuously 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.

Computer Simulation↗
PubMed2004

Base pairing at the stem-loop junction in the SL1 kissing complex of HIV-1 RNA: a thermodynamic study probed by molecular dynamics simulation.

The thermodynamics of the opening/closure process of a GC base pair located at the stem-loop junction of the SL1 sequence from HIV-1(Lai) genomic RNA was investigated in the context of a loop-loop homodimer (or kissing complex) using molecular dynamics simulation. The free energy, enthalpy and entropy changes for the closing reaction are 0 kcal x mol(-1), -11 kcal x mol(-1) and -0.037 kcal x mol(-1) x K(-1) at 300 degrees K respectively. Furthermore it is found that the free energy change is the same for the formation of a 11 nucleotide loop closed with UG and for the formation of a 9 nucleotide loop closed with GC. Our study evidences the high flexibility of the nucleotides at the stem-loop junction explaining the weak stability of this structure.

Base Pairing↗
PubMed2005

Molecular dynamics simulation of sucrose- and trehalose-coated carboxy-myoglobin.

We performed a room temperature molecular dynamics (MD) simulation on a system containing 1 carboxy-myoglobin (MbCO) molecule in a sucrose-water matrix of identical composition (89% [sucrose/(sucrose + water)] w/w) as for a previous trehalose-water-MbCO simulation (Cottone et al., Biophys J 2001;80:931-938). Results show that, as for trehalose, the amplitude of protein atomic mean-square fluctuations, on the nanosecond timescale, is reduced with respect to aqueous solutions also in sucrose. A detailed comparison as a function of residue number evidences mobility differences along the protein backbone, which can be related to a different efficacy in bioprotection. Different heme pocket structures are observed in the 2 systems. The joint distribution of the magnitude of the electric field at the CO oxygen atom and of the angle between the field and the CO unit vector shows a secondary maximum in sucrose, absent in trehalose. This can explain the CO stretching band profile (A substates distribution) differences evidenced by infrared spectroscopy in sucrose- and trehalose-coated MbCO (Giuffrida et al., J Phys Chem B 2004;108:15415-15421), and in particular the appearance of a further substate in sucrose. Analysis of hydrogen bonds at the protein-solvent interface shows that the fraction of water molecules shared between the protein and the sugar is lower in sucrose than in trehalose, in spite of a larger number of water molecules bound to the protein in the former system, thus indicating a lower protein-matrix coupling, as recently observed by Fourier transform infrared (FTIR) experiments (Giuffrida et al., J Phys Chem B 2004;108:15415-15421).

Binding Sites↗
PubMed2006

Structural and dynamical properties of different protonated states of mutant HIV-1 protease complexed with the saquinavir inhibitor studied by molecular dynamics simulations.

To understand the basis of drug resistance, particularly of the HIV-1 PR, three molecular dynamics (MD) simulations of HIV-1 PR mutant species, G48V, complexed with saquinavir (SQV) in explicit aqueous solution with three protonation states, diprotonation on Asp25 and Asp25' (Di-pro) and monoprotonation on each Asp residue (Mono-25 and Mono-25'). For all three states, H-bonds between saquinavir and HIV-1 PR were formed only in the two regions, flap and active site. It was found that conformation of P2 subsite of SQV in the Mono-25 state differs substantially from the other two states. The rotation about 177 degrees from the optimal structure of the wild type was observed, the hydrogen bond between P2 and the flap residue (Val48) was broken and indirect hydrogen bonds with the three residues (Asp29, Gly27, and Asp30) were found instead. In terms of complexation energies, interaction energy of -37.3 kcal/mol for the Mono-25 state is significantly lower than those of -30.7 and -10.7kcal/mol for the Mono-25' and Di-pro states, respectively. It was found also that protonation at the Asp25 leads to a better arrangement in the catalytic dyad, i.e., the Asp25-Asp25' interaction energy of -8.8 kcal/mol of the Mono-25 is significantly lower than that of -2.6kcal/mol for the Mono-25' state. The above data suggest us to conclude that interaction in the catalytic area should be used as criteria to enhance capability in drug designing and drug screening instead of using the total inhibitor/enzyme interaction.

Binding Sites↗
PubMed2002

Molecular dynamics simulations of the acyl-enzyme and the tetrahedral intermediate in the deacylation step of serine proteases.

Despite the availability of many experimental data and some modeling studies, questions remain as to the precise mechanism of the serine proteases. Here we report molecular dynamics simulations on the acyl-enzyme complex and the tetrahedral intermediate during the deacylation step in elastase catalyzed hydrolysis of a simple peptide. The models are based on recent crystallographic data for an acyl-enzyme intermediate at pH 5 and a time-resolved study on the deacylation step. Simulations were carried out on the acyl enzyme complex with His-57 in protonated (as for the pH 5 crystallographic work) and deprotonated forms. In both cases, a water molecule that could provide the nucleophilic hydroxide ion to attack the ester carbonyl was located between the imidazole ring of His-57 and the carbonyl carbon, close to the hydrolytic position assigned in the crystal structure. In the "neutral pH" simulations of the acyl-enzyme complex, the hydrolytic water oxygen was hydrogen bonded to the imidazole ring and the side chain of Arg-61. Alternative stable locations for water in the active site were also observed. Movement of the His-57 side-chain from that observed in the crystal structure allowed more solvent waters to enter the active site, suggesting that an alternative hydrolytic process directly involving two water molecules may be possible. At the acyl-enzyme stage, the ester carbonyl was found to flip easily in and out of the oxyanion hole. In contrast, simulations on the tetrahedral intermediate showed no significant movement of His-57 and the ester carbonyl was constantly located in the oxyanion hole. A comparison between the simulated tetrahedral intermediate and a time-resolved crystallographic structure assigned as predominantly reflecting the tetrahedral intermediate suggests that the experimental structure may not precisely represent an optimal arrangement for catalysis in solution. Movement of loop residues 216-223 and P3 residue, seen both in the tetrahedral simulation and the experimental analysis, could be related to product release. Furthermore, an analysis of the geometric data obtained from the simulations and the pH 5 crystal structure of the acyl-enzyme suggests that since His-57 is protonated, in some aspects, this crystal structure resembles the tetrahedral intermediate.

Acylation↗
PubMed2004

Molecular dynamics simulation of the C-terminal sterile alpha-motif domain of human p73alpha: evidence of a dynamical relationship between helices 3 and 5.

We used molecular dynamics simulation to evaluate the association properties of C-terminal sterile alpha-motif (SAM) domain of human p73alpha. To test the dimerization propensity of this structure we carried out four simulations: EphB2 X-ray dimer, p73 modeled dimer, p73 NMR monomer, and p73 modeled monomer with an elongated helix 5. The results show a direct interaction between helix 5 and helix 3 since a conformational collapse of helix 3 is observed when dimer contact and/or an elongation of helix 5 is introduced by modeling in p73 SAM domain. On the basis of these results we suggest that the recognition properties of the SAM domains may be modulated by the conformational state of helix 5.

Amino Acid Motifs↗
PubMed2006

Advantages of a Lowe-Andersen thermostat in molecular dynamics simulations.

The Lowe-Andersen thermostat is a momentum conserving and Galilean invariant analog of the Andersen thermostat. Like the Andersen thermostat it has the advantage of being local. We show that by using a minimal thermostat interaction radius in a molecular dynamics simulation, it perturbs the system dynamics to a far lesser extent than the Andersen method. This alleviates a well known drawback of the Andersen thermostat by allowing high thermostatting rates without the penalty of significantly suppressed diffusion in the system.

Journal Article↗
PubMed1992

A model of the molten globule state from molecular dynamics simulations.

It is generally accepted that a protein's primary sequence determines its three-dimensional structure. It has proved difficult, however, to obtain detailed structural information about the actual protein folding process and intermediate states. We present the results of molecular dynamics simulations of the unfolding of reduced bovine pancreatic trypsin inhibitor. The resulting partially "denatured" state was compact but expanded relative to the native state (11-25%); the expansion was not caused by an influx of water molecules. The structures were mobile, with overall secondary structure contents comparable to those of the native protein. The protein experienced relatively local unfolding, with the largest changes in the structure occurring in the loop regions. A hydrophobic core was maintained although packing of the side chains was compromised. The properties displayed in the simulation are consistent with unfolding to a molten globule state. Our simulations provide an in-depth view of this state and details of water-protein interactions that cannot yet be obtained experimentally.

Animals↗
PubMed2004

Shear-induced undulation of smectic-A: Molecular Dynamics simulations vs. analytical theory.

Experiments on a variety of systems have shown that layered liquids are unstable under shear even if the liquid layers are planes of constant velocity. We investigate the stability of smectic- A like liquids under shear using Molecular Dynamics simulations and a macroscopic hydrodynamic theory (including the layer normal and the director as independent variables). Both methods show an instability of the layers, which sets in above a critical shear rate. We find a remarkable qualitative and reasonable quantitative agreement between both methods for the spatial homogeneous state and the onset of the instability.

Computer Simulation↗
PubMed2005

Limited flexibility of lactose detected from residual dipolar couplings using molecular dynamics simulations and steric alignment methods.

The conformational flexibility of lactose in solution has been investigated by residual dipolar couplings (RDCs). One-bond carbon-proton and proton-proton coupling constants have been measured in two oriented media and interpreted in combination with molecular dynamics simulations (MD). Two different approaches, known as PALES (Zweckstetter et al., J. Am. Chem. Soc. 2000, 122, 3791-3792) and TRAMITE (Azurmendi et al., J. Am. Chem. Soc. 2002, 124, 2426-2427), have been used to determine the alignment tensor from a shape-induced alignment model with the oriented medium. The steric alignment of the structures from several MD trajectories has provided ensemble averaged RDCs that have been compared with the experimental ones. The obtained results reveal the almost exclusive presence of a major low energy region defined as syn-phi/syn-psi (> 97%), for which sampling occurs in a dynamic manner. This result satisfactorily agrees with that determined by standard NOE-based methods.

Carbohydrate Conformation↗
PubMed2003

Atomistic molecular dynamics simulations of chemical force microscopy.

Chemical force microscopy and related force measurement techniques have emerged as powerful tools for studying fundamental interactions central to understanding adhesion and tribology at the molecular scale. However, detailed interpretation of these interactions requires knowledge of chemical and physical processes occurring in the region of the tip-sample junction that experiments cannot provide, such as atomic-scale motions and distribution of forces. In an effort to address some of these open issues, atomistic molecular dynamics simulations were performed modeling a chemical force microscope stylus covered with a planar C12 alkylthiolate self-assembled monolayer (SAM) interacting with a solid wall. A complete loading-unloading sequence was simulated under conditions of near-constant equilibrium, approximating the case of infinitely slow tip motion. In the absence of the solid wall, the stylus film existed in a fluid state with structural and dynamic properties similar to those of the analogous planar SAM at an elevated temperature. When the wall was brought into contact with the stylus and pressed against it, a series of reversible changes occurred culminating with solidification of the SAM film at the largest compressive force. During loading, the chemical composition of the contact changed, as much of the film's interior was exposed to the wall. At all tip heights, the distribution of forces within the contact zone was uneven and subject to large local fluctuations. Analysis using the Johnson-Kendall-Roberts, Derjaguin-Muller-Toporov, and Hertz contacts mechanics models revealed significant deviations from the simulation results, with the JKR model providing best overall agreement. Some of the discrepancies found would be overlooked in an actual experiment, where, unlike the simulations, contact area is not separately known, possibly producing a misleading or incorrect interpretation of experimental results. These shortcomings may be improved upon by using a model that correctly accounts for the finite thickness of the compliant components and nonlinear elastic effects.

Journal Article↗
PubMed2003

Molecular dynamics simulations of pentapeptides at interfaces: salt bridge and cation-pi interactions.

Peptide-membrane interactions are important for understanding the binding, partitioning, and folding of membrane proteins; the activity of antimicrobial and fusion peptides; and a number of other processes. We describe molecular dynamics simulations (10-25 ns) of two pentapeptides Ace-WLXLL (with X = Arg or Lys side chain) (White, S. H., and Wimley, W.C. (1996) Nat. Struct. Biol. 3, 842-848) in water and three different membrane mimetic systems: (i) a water/cyclohexane interface, (ii) water-saturated octanol, and (iii) a solvated dioleoylphosphatidylcholine bilayer. A salt bridge is found between the protonated Arg or Lys side chains with the carboxyl terminus at the three interfaces. In water/cyclohexane, the salt bridge is most exposed to the water phase and least stable. In water/octanol and the lipid bilayer systems, the salt bridge once formed persists throughout the simulations. In the lipid bilayer, the salt bridge is more stable when the peptide penetrates deeper into the bilayer. In one of two peptides, a cation-pi interaction between the Arg and the Trp side chains is stable in the lipid bilayer for about 15 ns before breaking. In all cases, the conformations of the peptides are restricted by their presence at the interface and can be assigned to a few major conformational clusters. Side chains facing the water phase are most mobile. In the lipid bilayer, the peptides remain in the interface area, where they overlap with the carbonyl area of the lipid bilayer and perturb the local density profile of the bilayer. The tryptophan side chain remains in the water-lipid interface, where it interacts with the lipid choline group and forms hydrogen bonds with the ester carbonyl of the lipid and with water in the interface.

Arginine↗
PubMed2005

Characterization of the molten globule state of retinol-binding protein using a molecular dynamics simulation approach.

Retinol-binding protein transports retinol, and circulates in the plasma as a macromolecular complex with the protein transthyretin. Under acidic conditions retinol-binding protein undergoes a transition to the molten globule state, and releases the bound retinol ligand. A biased molecular dynamics simulation method has been used to generate models for the ensemble of conformers populated within this molten globule state. Simulation conformers, with a radius of gyration at least 1.1 A greater than that of the native state, contain on average 37%beta-sheet secondary structure. In these conformers the central regions of the two orthogonal beta-sheets that make up the beta-barrel in the native protein are highly persistent. However, there are sizable fluctuations for residues in the outer regions of the beta-sheets, and large variations in side chain packing even in the protein core. Significant conformational changes are seen in the simulation conformers for residues 85-104 (beta-strands E and F and the E-F loop). These changes give an opening of the retinol-binding site. Comparisons with experimental data suggest that the unfolding in this region may provide a mechanism by which the complex of retinol-binding protein and transthyretin dissociates, and retinol is released at the cell surface.

Computer Simulation↗
PubMed2005

Correlated dynamics determining x-ray diffuse scattering from a crystalline protein revealed by molecular dynamics simulation.

The dynamical origin of the x-ray diffuse scattering by crystals of a protein, Staphylococcal nuclease, is determined using molecular dynamics simulation. A smooth, nearly isotropic scattering shell at originates from equal contributions from correlations in nearest-neighbor water molecule dynamics and from internal protein motions, the latter consisting of -helix pitch and inter--strand fluctuations. Superposed on the shell are intense, three-dimensional scattering features that originate from a very small number of slowly varying (>10 ns) collective motions. The individual three-dimensional features are assigned to specific collective motions in the protein, and some of these explicitly involve potentially functional active-site deformations.

Computer Simulation↗