Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “molecular 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 73 records · Page 4Linked to original sources
PubMed2001

Cholesterol effects on the phosphatidylcholine bilayer nonpolar region: a molecular simulation study.

A 15-ns molecular dynamics (MD) simulation of the fully hydrated dimyristoylphosphatidylcholine-cholesterol (DMPC-Chol) bilayer in the liquid-crystalline state was carried out to investigate the effect of Chol on the hydrocarbon chain region of the bilayer. The last 8-ns fragment of the generated trajectory was used for analyses. As a reference system, a pure DMPC bilayer (M. Pasenkiewicz-Gierula, Y. Takaoka, H. Miyagawa, K. Kitamura, and A. Kusumi, 1999, Biophys. J. 76:1228-1240) simulated for 14 ns was used. The study shows that a Chol-induced increase of the bulk molecular order parameter along both beta- and gamma-chain is mainly caused by a decrease of the average tilt of the chains, because the bulk average number of gauche rotamers/myristoyl chain is not significantly changed by Chol. Nevertheless, for DMPCs located near Chol molecules both the number of gauche rotamers/chain and the chain tilt are decreased. The magnitude of the Chol effect on the PC alkyl chains depends, in addition to the PC-Chol distance, on the side of the Chol molecule (alpha- or beta-face) that the chains are in contact with. This study provides some new insight into the properties of the coexistence region of the partial phase diagram for DMPC-Chol bilayers.

Cholesterol↗
PubMed2003

Molecular simulation study of phospholipid bilayers and insights of the interactions with disaccharides.

Molecular simulations of hydrated dipalmitoylphosphatidylcholine lipid bilayers have been performed for temperatures in the range of 250-450 K. The area per headgroup increases with temperature from 58 to 77 A(2). Other properties such as hydration number, alkyl tail order parameter, diffusion coefficients, and radial distribution functions exhibit a clear dependence on temperature. Simulations of bilayers have also been performed in the presence of two disaccharides, namely trehalose and sucrose, at concentrations of up to 18 wt % (lipid-free basis). The simulated area per headgroup of the bilayer is not affected by the presence of the disaccharides, suggesting that the overall structure of the bilayer remains undisturbed. The results of simulations reveal that the interaction of disaccharide molecules with the bilayer occurs at the surface of the bilayer, and it is governed by the formation of multiple hydrogen bonds to specific groups of the lipid. Disaccharide molecules are observed to adopt specific conformations to fit onto the surface topology of the bilayer, often interacting with up to three different lipids simultaneously. At high disaccharide concentrations, the results of simulations indicate that disaccharides can serve as an effective replacement for water under anhydrous conditions, which helps explain their effectiveness as lyophilization agents for liposomes and cells.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed2005

Separation of CO2 and N2 by adsorption in C168 schwarzite: a combination of quantum mechanics and molecular simulation study.

Using a hierarchical multiscale approach combining quantum mechanics and molecular simulation, we have investigated the adsorption of pure CO(2) and N(2) and their mixture at room temperature in C(168) schwarzite, as a model for nanoporous carbons. First, the adsorbate-adsorbent interaction potential is determined using ab initio quantum mechanics computations, and then the adsorption is predicted using full atomistic Monte Carlo simulations. The extents of adsorption, adsorption energies, and isosteric heats of pure CO(2) and N(2) simulated with the ab initio potential are found to be higher than those with the empirical Steele potential that had been developed from gas adsorption on planar graphite. The inclusion of the electric quadrupole moment of adsorbate in simulation has no discernible effect on N(2) adsorption but results in a larger extent of CO(2) adsorption at high coverages. The selectivity of CO(2) over N(2) in the C(168) schwarzite from a model flue gas is predicted to be significantly larger with the ab initio potential than with the Steele potential. This illustrates the importance of an accurate adsorbate-adsorbent interaction potential in determining gas adsorption and suggests that nanoporous carbons might be useful for the separation of flue gases. As a comparison, the adsorption and selectivity of CO(2) and N(2) in ZSM-5 zeolites are also simulated with the experimentally validated potential parameters. The selectivity in the C(168) schwarzite predicted with the ab initio or Steele potential is found to be larger than the selectivity in all-silica ZSM-5, but less than that in Na-exchanged ZSM-5 zeolites.

Journal Article↗
PubMed2002

Reaction mechanism of soluble epoxide hydrolase: insights from molecular dynamics simulations.

Molecular dynamics simulations have been performed to gain insights into the catalytic mechanism of the hydrolysis of epoxides to vicinal diols by soluble epoxide hydrolase (sEH). The binding of a substrate, 1S,2S-trans-methylstyrene oxide, was studied in two conformations in the active site of the enzyme. It was found that only one is likely to be found in the active enzyme. In the preferred conformation the phenyl group of the substrate is pi-sandwiched between two aromatic residues, Tyr381 and His523, whereas the other conformation is pi-stacked with only one aromatic residue, Trp334. Two simulations were carried out to 1 ns for each conformation to evaluate the protonation state of active site residue His523. It was found that a protonated histidine is essential for keeping the active site from being disrupted. Long time scale, 4 ns, molecular dynamics simulation was done for the structure with the most likely combination of binding conformation and protonation state of His523. Near Attack Conformers (NACs) are present 5.3% of the time and nucleophilic attack on either epoxide carbon atom, approximately 75% on C(1) and approximately 25% on C(2), is found. A maximum of one hydrogen bond between the epoxide oxygen and either of the active site tyrosines, Tyr465 and Tyr381, is present, in agreement with experimental mutagenesis results that reveal a slight loss in activity if one tyrosine is mutated and essential loss of all activity upon double mutation of the two tyrosines in question. It was found that a hydrogen bond from Tyr465 to the substrate oxygen is essential for controlling the regioselectivity of the reaction. Furthermore, a relationship between the presence of this hydrogen bond and the separation of reactants was found. Two groups of amino acid segments were identified each as moving collectively. Furthermore, an overall anti-correlation was found between the movements of these two individually collectively moving groups, made up by parts of the cap-region, including the two tyrosines, and the site of the catalytic triad, respectively. This overall anti-correlated collective domain motion is, perhaps, involved in the conversion of E.NAC to E.TS.

Animals↗
PubMed2006

Anthramycin-DNA binding explored by molecular simulations.

The anticancer drug anthramycin inhibits replication and transcription processes by covalently binding to DNA. Here, we use molecular simulations to investigate the interaction between this ligand and the dodecanucleotide d[GCCAACGTTGGC](2). We start from the X-ray structure of the adduct anthramycin-d[CCAACGTTG*G](2), in which the drug binds covalently to guanine.1 We focus on the noncovalent complexes between the oligonucleotide and the anhydro and hydroxy forms of the drug. Molecular dynamics (MD) simulations show that only the hydroxy form lies in front of the reactive center for the whole simulation ( approximately 20 ns), while the anhydro form moves inside the minor groove to the nearest base pair after approximately 10 ns. This sliding process is associated to both energetic and structural relaxations of the complex. The accuracy of our computational setup is established by performing MD simulations of the covalent adduct and of a 14-mer complexed with anhydro-anthramycin. The MD simulations are complemented by hybrid Car-Parrinello quantum mechanics/molecular mechanics (QM/MM) simulations. These show that in the noncovalent complexes the electric field due to DNA polarizes the hydroxy and, even more, the anhydro form of the drug as to favor a nucleophilic attack by the alkylating guanine. This suggests that the binding process may be characterized by a multistep pathway, catalyzed by the electric field of DNA.

Anthramycin↗
PubMed2003

Modeling of adhesion in tablet compression--I. Atomic force microscopy and molecular simulation.

Adhesion problems during tablet manufacturing have been observed to be dependent on many formulation and process factors including the run time on the tablet press. Consequently, problems due to sticking may only become apparent towards the end of the development process when a prolonged run on the tablet press is attempted for the first time. It would be beneficial to predict in a relative sense if a formulation or new chemical entity has the potential for adhesion problems early in the development process. It was hypothesized that favorable intermolecular interaction between the drug molecules and the punch face is the first step or criterion in the adhesion process. Therefore, the rank order of adhesion during tablet compression should follow the rank order of these energies of interaction. The adhesion phenomenon was investigated using molecular simulations and contact mode atomic force microscopy (AFM). Three model compounds were chosen from a family of "profen" compounds. Silicon nitride AFM tips were modified by coating a 20-nm iron layer on the surfaces by sputter coating. Profen flat surfaces were made by melting and recrystallization. The modified AFM probe and each profen surface were immersed in the corresponding profen saturated water during force measurements using AFM. The work of adhesion between iron and ibuprofen, ketoprofen, and flurbiprofen in vacuum were determined to be -184.1, -2469.3, -17.3 mJ. m(-2), respectively. The rank order of the work of adhesion between iron and profen compounds decreased in the order: ketoprofen > ibuprofen > flurbiprofen. The rank order of interaction between the drug molecules and the iron superlattice as predicted by molecular simulation using Cerius(2) is in agreement with the AFM measurements. It has been demonstrated that Atomic Force Microscopy is a powerful tool in studying the adhesion phenomena between organic drug compounds and metal surface. The study has provided insight into the adhesion problems occurring during tablet compression and a direction for continued study.

Adhesiveness↗
PubMed2003

Molecular simulations to determine the chelating mechanisms of various metal ions to the His-tag motif: a preliminary study.

In the present study, molecular simulations were performed to investigate the chelating mechanisms of various metal ions to the His-tag motifs with various His residues. The chelation mostly involved the i and i+2 His residues for Ni(2+), Zn(2+), Cu(2+), and Co(2+), while the cooperation of 3 His residues was necessary when Fe(3+) was involved in chelation with His-tags having more than 4 His residues. Metal ion was best fitted into the pocket formed by the imidazole nitrogens while it was about equally located among these nitrogen atoms. His-tag6 was found to have little effect on the structural integrity while the target protein contains more than 68 amino acid residues. Ni(2+) interacted with the imidazole nitrogen of His3 in the beginning of chelation, and then entered into the pocket formed by His3 and His5 at 4 ns during the 10 ns molecular dynamics simulations. The fast chelating process resulted in successful application of IMAC techniques in efficient protein purification.

Amino Acid Motifs↗
PubMed2006

Understanding the nonfouling mechanism of surfaces through molecular simulations of sugar-based self-assembled monolayers.

This paper presents a molecular simulation study of the interactions of a protein (lysozyme) with self-assembled monolayers (SAMs) of mannitol and sorbitol terminated alkanethiols in the presence of explicit water molecules and ions. The all-atom simulations were performed to calculate the force generated on the protein as a function of its distance above the SAM surfaces. The structural and dynamic properties of water molecules both above the SAM surfaces and around the SAM head groups were analyzed to provide a better understanding of the nonfouling behavior of the sugar-based SAM surfaces. Results from this work suggest that both mannitol and sorbitol SAMs generate a tightly bound, structured water layer around the SAM chains. This hydration layer creates a repulsive force on the protein when it approaches the surface, resulting in a nonfouling surface despite the presence of hydrogen-bond donor groups. This work demonstrates the importance of strong surface-water interactions for surface resistance to nonspecific protein adsorption.

Adsorption↗
PubMed1994

"Simulated molecular evolution" or computer-generated artifacts?

1. The authors define a function with value 1 for the positive examples and 0 for the negative ones. They fit a continuous function but do not deal at all with the error margin of the fit, which is almost as large as the function values they compute. 2. The term "quality" for the value of the fitted function gives the impression that some biological significance is associated with values of the fitted function strictly between 0 and 1, but there is no justification for this kind of interpretation and finding the point where the fit achieves its maximum does not make sense. 3. By neglecting the error margin the authors try to optimize the fitted function using differences in the second, third, fourth, and even fifth decimal place which have no statistical significance. 4. Even if such a fit could profit from more data points, the authors should first prove that the region of interest has some kind of smoothness, that is, that a continuous fit makes any sense at all. 5. "Simulated molecular evolution" is a misnomer. We are dealing here with random search. Since the margin of error is so large, the fitted function does not provide statistically significant information about the points in search space where strings with cleavage sites could be found. This implies that the method is a highly unreliable stochastic search in the space of strings, even if the neural network is capable of learning some simple correlations. 6. Classical statistical methods are for these kind of problems with so few data points clearly superior to the neural networks used as a "black box" by the authors, which in the way they are structured provide a model with an error margin as large as the numbers being computed.7. And finally, even if someone would provide us with a function which separates strings with cleavage sites from strings without them perfectly, so-called simulated molecular evolution would not be better than random selection.Since a perfect fit would only produce exactly ones or zeros,starting a search in a region of space where all strings in the neighborhood get the value zero would not provide any kind of directional information for new iterations. We would just skip from one point to the other in a typical random walk manner.

Amino Acid Sequence↗
PubMed2004

Monte Carlo molecular simulation of the hydration of K-montmorillonite at 353 K and 625 bar.

Monte Carlo molecular simulations of the hydration of K-saturated Wyoming-type montmorillonite at constant stress in the NPzzT ensemble and at constant chemical potential in the grand canonical muVT ensemble, under basin-like conditions of 353 K and 625 bar, show a strong tendency of the K+ ions to adhere to the siloxane surface, forming predominant inner-sphere complexes with tetrahedral oxygen atoms and adsorbed water molecules. Simulations in the grand canonical ensemble predict that none of the K-montmorillonite hydrates, the one-, two-, and three-layer hydrates, are stable in this environment of high depth, temperature, and pressure. The most nearly stable configuration corresponds to the one-layer hydrate, characterized by a d001 spacing of 12.75 A, the adsorbed water being 60 molecules/layer or 180.83 mg of H2O/g of clay, an internal energy of -22.73 kcal/mol, an interlayer density of 0.365 g/mL, and a pressure tensor, Pzz, of 1999.9 bar. The interlayer structure consists of two close layers of water molecules 0.50 A from the midplane, with broad shoulders on the sides, the protons oriented toward the midplane and the siloxane surfaces, and the K+ ions close to the clay surfaces and on the interlayer midplane.

Adsorption↗
PubMed2006

Molecular simulations of cotranslational protein folding: fragment stabilities, folding cooperativity, and trapping in the ribosome.

Although molecular simulation methods have yielded valuable insights into mechanistic aspects of protein refolding in vitro, they have up to now not been used to model the folding of proteins as they are actually synthesized by the ribosome. To address this issue, we report here simulation studies of three model proteins: chymotrypsin inhibitor 2 (CI2), barnase, and Semliki forest virus protein (SFVP), and directly compare their folding during ribosome-mediated synthesis with their refolding from random, denatured conformations. To calibrate the methodology, simulations are first compared with in vitro data on the folding stabilities of N-terminal fragments of CI2 and barnase; the simulations reproduce the fact that both the stability and thermal folding cooperativity increase as fragments increase in length. Coupled simulations of synthesis and folding for the same two proteins are then described, showing that both fold essentially post-translationally, with mechanisms effectively identical to those for refolding. In both cases, confinement of the nascent polypeptide chain within the ribosome tunnel does not appear to promote significant formation of native structure during synthesis; there are however clear indications that the formation of structure within the nascent chain is sensitive to location within the ribosome tunnel, being subject to both gain and loss as the chain lengthens. Interestingly, simulations in which CI2 is artificially stabilized show a pronounced tendency to become trapped within the tunnel in partially folded conformations: non-cooperative folding, therefore, appears in the simulations to exert a detrimental effect on the rate at which fully folded conformations are formed. Finally, simulations of the two-domain protease module of SFVP, which experimentally folds cotranslationally, indicate that for multi-domain proteins, ribosome-mediated folding may follow different pathways from those taken during refolding. Taken together, these studies provide a first step toward developing more realistic methods for simulating protein folding as it occurs in vivo.

Binding Sites↗
PubMed2006

Molecular simulations of benzene and PAH interactions with soot.

Molecular mechanics simulations and ab initio calculations were performed to investigate the mechanism of PAH-soot adsorption. Partitioning of benzene, naphthalene, fluorene, phenanthrene, anthracene, pyrene, fluoranthene, benzo[a]anthracene, benzo[k]fluoranthene, benzo[a]pyrene, and benzo[g,h,i]perylene between water and soot was modeled with classical and quantum mechanical calculations in order to determine a method for predicting log(K(d)) values. In both cases, the predicted mechanism of adsorption is interaction of the pi-electrons in the PAH and soot (i.e., pi-pi van der Waals forces). Solvation energies, the energy difference between the solute in the gas phase and in the model aqueous phase, calculated with molecular mechanics did not follow the observed solubilities of the PAHs. Molecular dynamics simulations overestimate the favorability of PAHs in the aqueous phase. Hence, the partitioning between the aqueous phase and soot does not accurately correlate with observed log(K(d)) values. Models of PAH adsorption using structures from molecular mechanics and energies from ab initio calculations do produce water-soot partitioning energies that correlate well with observed log(K(d)) values. The log(K(d)) values for benzene, anthracene, fluorene, and benzo[a]pyrene were predicted based on the correlation between calculated partitioning energies and observed log(K(d)) values. Results presented here suggest that partitioning of PAHs onto soot should depend on the size of the PAH, the planarity of PAH molecule, and the aromaticity of the compound. The methodology developed bythis research can be used to predict PAH K(d) values that have not yet been measured.

Benzene↗
PubMed2006

Vibrational spectra of methane clathrate hydrates from molecular dynamics simulation.

Molecular dynamics simulations were performed on methane clathrate hydrates at ambient conditions. Thermal expansion results over the temperature range 60-300 K show that the unit cell volume increases with temperature in agreement with experiment. Power spectra were obtained at 273 K from velocity autocorrelation functions for selected atoms, and normal modes were assigned. The spectra were further classified according to individual atom types, allowing the assignment of contributions from methane molecules located in small and large cages within the structure I unit cell. The symmetric C-H stretch of methane in the small cages occurs at a higher frequency than for methane located in the large cages, with a peak separation of 14 cm(-1). Additionally, we determined that the symmetric C-H stretch in methane gas occurs at the same frequency as methane in the large cages. Results of molecular dynamics simulations indicate the use of power spectra obtained from the velocity autocorrelation function is a reliable method to investigate the vibrational behavior of guest molecules in clathrate hydrates.

Letter↗
PubMed2005

Molecular simulation of the homogeneous crystal nucleation of carbon dioxide.

We report on a molecular simulation study of the homogeneous nucleation of CO2 in the supercooled liquid at low pressure (P = 5 MPa) and for degrees of supercooling ranging from 32% to 60%. In all cases, regardless of the degree of supercooling, the structure of the crystal nuclei is that of the Pa3 phase, the thermodynamically stable phase. For the more moderate degree of supercooling of 32%, the nucleation is an activated process and requires a method to sample states of high free energy. In this work, we apply a series of bias potentials, which promote the ordering of the centers of mass of the molecules and allow us to gradually grow crystal nuclei. The reliability of the results so obtained is assessed by studying the evolution of the nuclei in the absence of any bias potential, and by determining their probability of growth. We estimate that the size of the critical nucleus, for which the probability of growth is 0.5, is approximately 240 molecules. Throughout the nucleation process, the crystal nuclei clearly exhibit a Pa3 structure, in apparent contradiction with Ostwald's rule of stages. The other polymorphs have a much larger free energy. This makes their formation highly unlikely and accounts for the fact that the nucleation of CO2 proceeds directly in the stable Pa3 structure.

Journal Article↗
PubMed2005

Molecular simulation study of the influence of small molecules on the dynamic and structural properties of phospholipid bilayers.

Molecular simulations were performed for phospholipid bilayers as model cell membranes to understand their dynamic and structural properties, in particular how those are affected in the presence of the cell-preservation agents trehalose and DMSO. Simulations covered a large range of temperatures and concentrations of trehalose and DMSO. A number of properties were examined, including the area per head group, lipid-tail order parameter, density profile, conformation, and binding of trehalose. Our results indicate that trehalose specifically interacts with the phospholipid bilayer by binding to the polar head groups (H-bonds). This binding serves as a bridge between the lipids and prevents their aggregation. The hydrophilic nature of trehalose also makes it a water substitute when binding to the phospholipids. At high trehalose concentrations, trehalose forms a glassy matrix and preserves the bilayer structure by 'freezing' it in its natured state, as the core of the bilayer is found to remain fluid. Simulations with DMSO revealed that DMSO can penetrate and freely diffuse through the bilayer, which is attributed to its unique characteristic of interacting favorably with both hydrophilic and hydrophobic domains.

Computer Simulation↗
PubMed2003

A computationally inexpensive modification of the point dipole electrostatic polarization model for molecular simulations.

We present an approximation, which allows reduction of computational resources needed to explicitly incorporate electrostatic polarization into molecular simulations utilizing empirical force fields. The proposed method is employed to compute three-body energies of molecular complexes with dipolar electrostatic probes, gas-phase dimerization energies, and pure liquid properties for five systems that are important in biophysical and organic simulations-water, methanol, methylamine, methanethiol, and acetamide. In all the cases, the three-body energies agreed with high level ab initio data within 0.07 kcal/mol, dimerization energies-within 0.43 kcal/mol (except for the special case of the CH(3)SH), and computed heats of vaporization and densities differed from the experimental results by less than 2%. Moreover, because the presented method allows a significant reduction in computational cost, we were able to carry out the liquid-state calculations with Monte Carlo technique. Comparison with the full-scale point dipole method showed that the computational time was reduced by 3.5 to more than 20 times, depending on the system in hand and on the desired level of the full-scale model accuracy, while the difference in energetic results between the full-scale and the presented approximate model was not great in the most cases. Comparison with the nonpolarizable OPLS-AA force field for all the substances involved and with the polarizable POL3 and q90 models for water and methanol, respectively, demonstrates that the presented technique allows reduction of computational cost with no sacrifice of accuracy. We hope that the proposed method will be of benefit to research employing molecular modeling technique in the biophysical and physical organic chemistry areas.

Algorithms↗
PubMed2007

Solubility of KF and NaCl in water by molecular simulation.

The solubility of two ionic salts, namely, KF and NaCl, in water has been calculated by Monte Carlo molecular simulation. Water has been modeled with the extended simple point charge model (SPC/E), ions with the Tosi-Fumi model and the interaction between water and ions with the Smith-Dang model. The chemical potential of the solute in the solution has been computed as the derivative of the total free energy with respect to the number of solute particles. The chemical potential of the solute in the solid phase has been calculated by thermodynamic integration to an Einstein crystal. The solubility of the salt has been calculated as the concentration at which the chemical potential of the salt in the solution becomes identical to that of the pure solid. The methodology used in this work has been tested by reproducing the results for the solubility of KF determined previously by Ferrario et al. [J. Chem. Phys. 117, 4947 (2002)]. For KF, it was found that the solubility of the model is only in qualitative agreement with experiment. The variation of the solubility with temperature for KF has also been studied. For NaCl, the potential model used predicts a solubility in good agreement with the experimental value. The same is true for the hydration chemical potential at infinite dilution. Given the practical importance of solutions of NaCl in water the model used in this work, whereas simple, can be of interest for future studies.

Computer Simulation↗
PubMed2001

Molecular simulation of the interaction of kappa-conotoxin-PVIIA with the Shaker potassium channel pore.

Molecular simulation techniques were appplied to predict the interaction of the voltage-dependent Shaker potassium channel with the channel-blocking toxin kappa-conotoxin-PVIIA (PVIIA). A structural thee-dimensional model of the extracellular vestibule of the potassium channel was constructed based on structural homologies with the bacterial potassium channel Kcsa, whose structure has been solved by X-ray crystallography. The docking of the PVIIA molecule was obtained by a geometric recognition algorithm, yielding 100 possible conformations. A series of residue-residue distance restraints, predicted from mutation-cycle experiments, were used to select a small set of a plausible channel-toxin complex models among the resulting possible conformations. The four final conformations, with similar characteristics, can explain most of the single-point mutation experiments done with this system. The models of the Shaker-PVIIA interaction predict two clusters of amino acids, critical for the binding of the toxin to the channel. The first cluster is the amino acids R2, I3, Q6 and K7 that form the plug of the toxin that interacts with the entrance to the selectivity filter of the channel. The second cluster of residues, R22, F23, N24 and K25, interacts with a channel region near to the external entrance of the pore vestibule. The consistency of the obtained models and the experimental data indicate that the Shaker-PVIIA complex model is reasonable and can be used in further biological studies such as the rational design of blocking agents of potassium channels and the mutagenesis of both toxins and potassium channels.

Algorithms↗