Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “MD 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 541 records · Page 30Linked to original sources
PubMed2004

Free energy and configurational entropy of liquid silica: fragile-to-strong crossover and polyamorphism.

Recent molecular dynamics (MD) simulations of liquid silica, using the "BKS" model [Phys. Rev. Lett. 64, 1955 (1990)]], have demonstrated that the liquid undergoes a dynamical crossover from super-Arrhenius, or "fragile" behavior, to Arrhenius, or "strong" behavior, as temperature T is decreased. From extensive MD simulations, we show that this fragile-to-strong crossover (FSC) can be connected to changes in the properties of the potential energy landscape, or surface (PES), of the liquid. To achieve this, we use thermodynamic integration to evaluate the absolute free energy of the liquid over a wide range of density and T. We use this free energy data, along with the concept of "inherent structures" of the PES, to evaluate the absolute configurational entropy S(c) of the liquid. We find that the temperature dependence of the diffusion coefficient and of S(c) are consistent with the prediction of Adam and Gibbs, including in the region where we observe the FSC to occur. We find that the FSC is related to a change in the properties of the PES explored by the liquid, specifically an inflection in the T dependence of the average inherent structure energy. In addition, we find that the high T behavior of S(c) suggests that the liquid entropy might approach zero at finite T, behavior associated with the so-called Kauzmann paradox. However, we find that the change in the PES that underlies the FSC is associated with a change in the T dependence of S(c) that elucidates how the Kauzmann paradox is avoided in this system. Finally, we also explore the relation of the observed PES changes to the recently discussed possibility that BKS silica exhibits a liquid-liquid phase transition, a behavior that has been proposed to underlie the observed polyamorphism of amorphous solid silica.

Journal Article↗
PubMed2000

X-Ray crystal structure and molecular dynamics simulations of silver hake parvalbumin (Isoform B).

Parvalbumins constitute a class of calcium-binding proteins characterized by the presence of several helix-loop-helix (EF-hand) motifs. In a previous study (Revett SP, King G, Shabanowitz J, Hunt DF, Hartman KL, Laue TM, Nelson DJ, 1997, Protein Sci 7:2397-2408), we presented the sequence of the major parvalbumin isoform from the silver hake (Merluccius bilinearis) and presented spectroscopic and structural information on the excised "EF-hand" portion of the protein. In this study, the X-ray crystal structure of the silver hake major parvalbumin has been determined to high resolution, in the frozen state, using the molecular replacement method with the carp parvalbumin structure as a starting model. The crystals are orthorhombic, space group C2221, with a = 75.7 A, b = 80.7 A, and c = 42.1 A. Data were collected from a single crystal grown in 15% glycerol, which served as a cryoprotectant for flash freezing at -188 degrees C. The structure refined to a conventional R-value of 21% (free R 25%) for observed reflections in the range 8 to 1.65 A [1 > 2sigma(I)]. The refined model includes an acetylated amino terminus, 108 residues (characteristic of a beta parvalbumin lineage), 2 calcium ions, and 114 water molecules per protein molecule. The resulting structure was used in molecular dynamics (MD) simulations focused primarily on the dynamics of the ligands coordinating the Ca2+ ions in the CD and EF sites. MD simulations were performed on both the fully Ca2+ loaded protein and on a Ca2+ deficient variant, with Ca2+ only in the CD site. There was substantial agreement between the MD and X-ray results in addressing the issue of mobility of key residues in the calcium-binding sites, especially with regard to the side chain of Ser55 in the CD site and Asp92 in the EF site.

Animals↗
PubMed1995

A model for transmembrane helix with a cis-proline in the middle.

The presence of a higher percentage of Proline in the transmembrane helices of transport proteins indicates that they are involved in the function of these integral membrane proteins (IMPs). In many cases, the possible involvement of cis-trans isomerization in function/folding of IMPs has been suggested. The introduction of cis-Pro in an ideal alpha-helix results in a helix-turn-helix motif. A molecular dynamics (MD) simulation is carried out on the sequence ACE-(ALA)10-cis-Pro-(ALA)10-NME with ideal alpha-helical structure to investigate if and how a straight helix can accommodate a cis-Pro. The analysis of the conformations accessed during MD simulation showed that the residues near cis-Pro can adopt alternate conformations other than the right-handed helical conformation such that an almost straight helix is obtained. This may have implications in the involvement of cis-trans isomerization in folding and/or function of IMPs.

Membrane Proteins↗
PubMed1995

Quantification and visualization of molecular surface flexibility.

Two new methods for the quantification and visualization of the flexibility of molecular surfaces are presented. Both methods rely on results of molecular dynamics (MD) simulations. Whereas method I is based on a simple but fast grid-counting algorithm, method II uses a mapping function that allows for a sharp and clear visualization of atomic RMS fluctuations on a molecular surface. To demonstrate the scope of the methods, MD simulations of two proteins, PTI and ubiquitin, were performed. The flexibility data are mapped onto the molecular surfaces of the proteins and visualized using texture mapping technology available on modern workstations.

Algorithms↗
PubMed1992

Solvent dielectric effects on protein dynamics.

Electron paramagnetic resonance (EPR) spectroscopy and molecular dynamics (MD) simulations were used to investigate the dynamics of alpha-chymotrypsin in solvents ranging in dielectric constant from 72 to 1.9. EPR measurements showed that motions in the vicinity of two spin-labeled amino acids (Met-192 and Ser-195) decreased dramatically with decreasing solvent dielectric constant, a trend consistent with changes in the electrostatic force between charged residues of the protein. EPR results and MD simulations revealed a very similar functional dependence between rates of motion in the protein and the dielectric constant of the bulk solvent; however, predicted motions of protein atoms were markedly faster than measured motions of the spin labels. MD calculations for dielectric constants of 5 and 72 showed the greatest differences near the outer surface of the protein. In general, at the lower dielectric constant many atoms of the protein move more slowly, and many of the slowest residues are near the exterior. These results suggest that altered dynamics may contribute to the unusual properties--e.g., modified stereoselectivities--of enzymes in nearly dry organic solvents.

Chymotrypsin↗
PubMed2006

Comparison of homology models and X-ray structures of the nuclear receptor CAR: assessing the structural basis of constitutive activity.

The constitutive androstane receptor (CAR) possesses an intrinsic basal activity whose structural basis has been analysed during the last decade. Recently, we published a homology model of the CAR ligand binding domain (LBD) based on the X-ray structures of the closely related pregnane X (PXR) and vitamin D (VDR) receptor. A detailed analysis of the homology model and molecular dynamics (MD) simulations afforded us to propose a potential mechanism underlying the constitutive activity of CAR. Almost simultaneously, X-ray structures of human and mouse CAR LBD were released. In the present study, a detailed analysis and comparison of homology model and X-ray structures is carried out in order to evaluate the quality and reliability of our homology modelling procedure. The hypothesis of the constitutive activity which we proposed on the basis of our modelling results was tested for consistency with the crystal structures. In addition, the features stated to be essential for the basal activity based on the X-ray data were investigated by means of molecular dynamics simulations. Our results show that the homology modelling procedure was able to predict the CAR LBD structure with high accuracy. Structural features that have been revealed as critical for constitutive activity in the model are also observed in the X-ray structures. Furthermore, the MD simulations of the CAR X-ray structures and a detailed analysis of other NRs clarify the role of distinct structural features that have been assigned an important role for the constitutive activity.

Amino Acid Sequence↗
PubMed2005

Molecular dynamics simulation of clustered DNA damage sites containing 8-oxoguanine and abasic site.

Clustered DNA damage sites induced by ionizing radiation have been suggested to have serious consequences to organisms, such as cancer, due to their reduced probability to be repaired by the enzymatic repair machinery of the cell. Although experimental results have revealed that clustered DNA damage sites effectively retard the efficient function of repair enzymes, it remains unclear as to what particular factors influence this retardation. In this study, approaches based on molecular dynamics (MD) simulation have been applied to examine conformational changes and energetic properties of DNA molecules containing clustered damage sites consisting of two lesioned sites, namely 7,8-dihydro-8-oxoguanine (8-oxoG) and apurinic/apyrimidinic (AP) site, located within a few base pairs of each other. After 1 ns of MD simulation, one of the six DNA molecules containing a clustered damage site develops specific characteristic features: sharp bending at the lesioned site and weakening or complete loss of electrostatic interaction energy between 8-oxoG and bases located on the complementary strand. From these results it is suggested that these changes would make it difficult for the repair enzyme to bind to the lesions within the clustered damage site and thereby result in a reduction of its repair capacity.

Algorithms↗
PubMed2006

Li+ transport in lithium sulfonylimide-oligo(ethylene oxide) ionic liquids and oligo(ethylene oxide) doped with LiTFSI.

The Li+ environment and transport in an ionic liquid (IL) comprised of Li+ and an anion of bis(trifluoromethanesulfonyl)imide anion (TFSI-) tethered to oligoethylene oxide (EO) (EO(12)TFSI-/Li+) were determined and compared to those in a binary solution of the oligoethylene oxide with LiTFSI salt (EO(12)/LiTFSI) by using molecular dynamics (MD) simulations and AC conductivity measurements. The latter revealed that the AC conductivity is 1 to 2 orders of magnitude less in the IL compared to the oligoether/salt binary electrolyte with greater differences being observed at lower temperatures. The conductivity of these electrolytes was accurately predicted by MD simulations, which were used in conjunction with a microscopic model to determine mechanisms of Li+ transport. It was discerned that structure-diffusion of the Li+ cation in the binary electrolyte (EO(12)/LiTFSI-) was similar to that in EO(12)TFSI-/Li+ IL at high temperature (>363 K), thus, one can estimate conductivity of IL at this temperature range if one knows the structure-diffusion of Li+ in the binary electrolyte. However, the rate of structure-diffusion of Li+ in IL was found to slow more dramatically with decreasing temperature than in the binary electrolyte. Lithium motion together with EO(12) solvent accounted for 90% of Li+ transport in EO(12)/LiTFSI-, while the Li+ motion together with the EO(12)TFSI- anion contributed approximately half to the total Li+ transport but did not contribute to the charge transport in IL.

Journal Article↗
PubMed1997

Molecular modeling of the neurophysin I/oxytocin complex.

Neurophysins I and II (NPI and NPII) act in the neurosecretory granules as carrier proteins for the neurophyseal hormones oxytocin (OT) and vasopressin (VP), respectively. The NPI/OT functional unit, believed to be an (NPI/OT)2 heterotetramer, was modeled using low-resolution structure information, viz. the C alpha carbon atom coordinates of the homologous NPII/dipeptide complex (file 1BN2 in the Brookhaven Protein Databank) as a template. Its all-atom representation was obtained using standard modeling tools available within the INSIGHT/Biopolymer modules supplied by Biosym Technologies Inc. A conformation of the NPI-bound OT, similar to that recently proposed in a transfer NOE experiment, was docked into the ligand-binding site by a superposition of its Cys1-Tyr2 fragment onto the equivalent portion of the dipeptide in the template. The starting complex for the initial refinements was prepared by two alternative strategies, termed Model I and Model II, each ending with a approximately 100 ps molecular dynamics (MD) simulation in water using the AMBER 4.1 force field. The free homodimer NPI2 was obtained by removal of the two OT subunits from their sites, followed by a similar structure refinement. The use of Model I, consisting of a constrained simulated annealing, resulted in a structure remarkably similar to both the NPII/dipeptide complex and a recently published solid-state structure of the NPII/OT complex. Thus, Model I is recommended as the method of choice for the preparation of the starting all-atom data for MD. The MD simulations indicate that, both in the homodimer and in the heterotetramer, the 3(10)-helices demonstrate an increased mobility relative to the remaining body of the protein. Also, the C-terminal domains in the NPI2 homodimer are more mobile than the N-terminal ones. Finally, a distinct intermonomer interaction is identified, concentrated around its most prominent, although not unique, contribution provided by an H-bond from Ser25 O gamma in one NPI unit to Glu81 O epsilon in the other unit. This interaction is present in the heterotetramer (NPI/OT)2 and absent or weak in the NPI2 homodimer. We speculate that this interaction, along with the increased mobility of the 3(10)-helices and the carboxy domains, may contribute to the allosteric communication between ligand binding and NPI dimerization.

Allosteric Regulation↗
PubMed2005

Integral equation theory for atactic polystyrene melt with a coarse-grained model.

In this work, an integral equation approach to investigate the atactic polystyrene (aPS) melt based on polymer reference interaction site model (PRISM) theory is proposed. The intramolecular structure factors, required as input to PRISM theory, are obtained from the semiflexible chain model. With a novel coarse-graining procedure and the explicit-atom molecular-dynamics (MD) simulations for aPS, the parameters needed for the coarse-grained model are obtained by using an automatic simplex optimization. These parameters can be used to describe the structure and thermodynamic properties of the complex aPS melt and good agreement is obtained between the theory and MD simulations. The proposed integral equation approach provides a basis for describing the structure and properties of PS nanocomposites where the application of molecular simulation is difficult.

Algorithms↗
PubMed2004

Reorganization in apo- and holo-beta-lactoglobulin upon protonation of Glu89: molecular dynamics and pKa calculations.

Molecular dynamics (MD) simulations starting from crystallographic data allowed us to directly account for the effects of the protonation state of Glu89 on the conformational stability of apo- and holo-beta-lactoglobulin (BLG). In apo-BLG simulations starting from the protonated crystal structure, we observe a long-lived H-bond interaction between the protonated Glu89 and Ser116. This interaction, sequestering the proton from the aqueous medium, explains a pK(half) value evaluated at pH 7.3 by continuum electrostatics/Monte Carlo computation on MD data, using linear response approximation. A very large root-mean-square deviation (RMSD; 5.11 A) is observed for the EF loop between protonated and unprotonated apo-BLG. This results from a quite different orientation of the EF loop that acts either as a closed or as an open lid above the protein calyx. Proton exchange by Glu89 in apo- but not in holo-BLG is associated with a reorganization energy of 4.7 kcal/mol. A 3-ns MD simulation starting from the crystal structure of protonated apo-BLG, but considering the Glu89 as unprotonated, shows the progressive opening of the lid giving rise to the Tanford transition. In both holo-BLG forms, the lid is most probably held in place by hydrophobic interactions of amino acid side-chains of the EF loop with the palmitate hydrocarbon tail.

Animals↗
PubMed1998

Solution conformation of methylated macrolide antibiotics roxithromycin and erythromycin using NMR and molecular modelling. Ribosome-bound conformation determined by TRNOE and formation of cytochrome P450-metabolite complex.

Conformational study of methylated derivatives of macrolide antibiotics roxithromycin (6-OMe-roxithromycin and 6,11-OMe-roxithromycin) has been achieved by NMR in solution and molecular dynamics (MD) simulations and compared to 6-OMe-erythromycin (clarithromycin). A complete conformational study by NMR has been led by determination of homonuclear coupling constants and NOEs. Heteronuclear 1H-13C coupling constants were also measured to investigate the orientation of the sugar moieties with respect to the erythronolide. MD simulations were performed using the crystallographic coordinates as the starting conformation. For each compound, experimental results were compared to calculated conformations in order to identify eventual conformational equilibrium in solution. It is shown that the effect of the methylation is opposite for roxithromycin compared to erythromycin especially on motional properties as the roxithromycin derivatives gain in mobility while the erythromycin derivatives behaves as a more restrained molecule. The study of macrolide-ribosome interactions has been investigated using transferred NOESY 1H NMR experiments and the conformations weakly bound to bacterial ribosomes were determined. Biological interactions of these compounds with membranar liver protein cytochrome P450 was also discussed with regard to their structural properties.

Anti-Bacterial Agents↗
PubMed2002

Intermolecular interactions between cholecystokinin-8 and the third extracellular loop of the cholecystokinin-2 receptor.

The structure of the third extracellular loop of the human cholecystokinin-2 receptor, CCK2-R(352-379), and its interactions with the C-terminal octapeptide of cholecystokinin (CCK-8) have been determined by high-resolution NMR and computer simulations. In the presence of dodecylphosphocholine micelles, the structure of the receptor fragment consisted of three helices, with the first and third corresponding to residues of the extracellular ends of transmembrane helices (TM) 6 and 7, respectively. The central, extracellular helix, consisting of residues 363-368, was found to be closely associated with the membrane mimetic used during the spectroscopic studies and molecular dynamics (MD) simulations. Upon titration of CCK-8 to the receptor domain, chemical shift perturbation and intermolecular NOEs (Trp30, Met31 of CCK-8 and P371, F374 of CCK2-R) indicated the formation of a stable complex and specific ligand/receptor interactions. Using the NOE-generated intermolecular contact points, extensive MD simulations of CCK-8 bound to the CCK2 receptor were carried out. The results, with CCK-8 in close proximity to TM7, differ from previous structural studies of CCK-8 association with CCK1-R, in which the ligand formed a number of interactions with TM6. These differences may play a role in the ligand specificity displayed by the CCK1 and CCK2 receptor subtypes.

Amino Acid Sequence↗
PubMed2004

Structural and dynamic effects of alpha-helix deletion in Sso7d: implications for protein thermal stability.

Sso7d is a 62-residue protein from the hyperthemophilic archaeon Sulfolobus solfataricus with a denaturation temperature close to 100 degrees C around neutral pH. An engineered form of Sso7d truncated at leucine 54 (L54Delta) is significantly less stable, with a denaturation temperature of 53 degrees C. Molecular dynamics (MD) studies of Sso7d and its truncated form at two different temperatures have been performed. The results of the MD simulations at 300 K indicate that: (1) the flexibility of Sso7d chain at 300 K agrees with that detected from X-ray and NMR structural studies; (2) L54Delta remains stable in the native folded conformation and possesses an overall dynamic behavior similar to that of the parent protein. MD simulations performed at 500 K, 10 ns long, indicate that, while Sso7d is in-silico resistant to high temperature, the truncated variant partially unfolds, revealing the early phases of the thermal unfolding pathway of the protein. Analysis of the trajectories of L54Delta suggests that the unzipping of the N-terminal and C-terminal beta-strands should be the first event of the unfolding pathway, and points out the regions more resistant to thermal unfolding. These findings allow one to understand the role played by specific interactions connecting the two ends of the chain for the high thermal stability of Sso7d, and support recent hypotheses on its folding mechanism emerged from site-directed mutagenesis studies.

Amino Acids↗
PubMed2005

Molecular dynamics investigations of the polysaccharide scleroglucan: first study on the triple helix structure.

Explicit solvent molecular dynamics (MD) simulations on the triple helix of the polysaccharide Scleroglucan (Sclg) at two temperatures (273 and 300 K) were carried out. Owing to the complexity of the system, a united-atom force field, based on the properly modified GROMACS parameters, was adopted. To test these parameters for our system, MD simulations of the two disaccharidic units, representing the main chain and the side-chain linkages of the Sclg repeating unit, were performed and the results were compared with the literature data. The simulated triple helix of Sclg retained the main experimentally determined features of the polymer. The residence times of the solvent molecules at 273 and 300 K were analyzed. The results show that the more internal water molecules, interacting with the core of the Sclg triplex are not influenced substantially by changing the temperature, on the contrary the water molecules, interacting with the side-chain glucose residues show more significant differences. These data suggest that the more external water molecules, interacting with the side chain, play a major role in the conformational transition experimentally observed at low temperature.

Carbohydrate Conformation↗
PubMed2004

Searching for nanostructures of the cubic mesophase of liquid crystal molecules, BABH8.

Nanostructures of a thermotropic cubic phase forming liquid crystal compound, 1,2-bis-[4-n-octyloxy-benzoyl]-hydrazine was studied by molecular dynamics (MD) simulations. A model for its cubic phase structure was proposed, which was constructed from the building unit of a locally orientational ordered "bundle" consistent with the cell parameter and the space group (Ia3d) from the recent x-ray results. Stability of the model structure was studied by multinanosecond MD simulations. A periodic nanostructure with 2.6 nm periodicity [coincides with the Ia3d (211) reflection] was kept up to 60 ns in the reduced pressure simulation which realizes the experimental value of density. However, the calculated fourth-rank order parameter shows that the simulated final state above does not have cubic orientational symmetry but rather isotropic symmetry.

Journal Article↗
PubMed2003

Molecular dynamics and quantum chemical studies on the catalytic mechanism of Delta5-3-ketosteroid isomerase: the catalytic diad versus the cooperative hydrogen bond mechanism.

To further understand Delta(5)-3-ketosteroid isomerase (KSI) catalysis, we carried out molecular dynamics (MD) simulations of the KSI dimer ligated with a substrate and reaction intermediate analogue and high level ab initio calculations on relevant enzymatic reaction models. Simulation of the enzyme-substrate complex dimer systems showed asymmetric dynamics between the two monomers, in which the hydrogen bond pattern between the substrate and active site residues in the first and the second subunits supported the cooperative hydrogen bond (CH) and the catalytic diad (CD) mechanisms, respectively. On the other hand, only the CH mechanism was supported in the MD simulation of the enzyme-intermediate complex dimer. From MP2/6-31+G**//RHF/6-31G** calculations, we found the kinetic barriers for the two reaction mechanisms were similar. The CH route afforded a greater stabilization to the enolate intermediate than did the CD counterpart. Thus, the present computational studies indicate that the CH mechanism would be favored over the CD one in the catalytic action of KSI. However, the latter could not be ruled out conclusively because of the explicit appearance of a CD configuration in the MD trajectories of the enzyme-substrate complex and because of the similar intrinsic activation barrier for the CH and CD mechanisms. The appearance of configurations that favor the CD pathway is rationalized in terms of a model in which the KSI-substrate complex does not have a strong preference for one hydrogen bonding pattern over another, while the KSI-intermediate complex favors a cooperative hydrogen bond pattern in order to stabilize the reaction intermediate. This hypothesis is supported by the ab initio calculations which indicate that the CH intermediate is more stable than the CD one by approximately 6.3 kcal/mol.

Catalysis↗
PubMed2026

Genomic and structural analysis of dacB variants associated with cephalosporin resistance in Pseudomonas aeruginosa.

The rise of resistance to fourth-generation cephalosporin in Pseudomonas aeruginosa (P. aeruginosa) is a global concern. The resistance is largely driven by variants of chromosomally encoded AmpC β-lactamase, known as Pseudomonas-derived cephalosporinase (PDC), which arise from the mutations in the ampC gene. In addition, alteration in dacB, which encode the penicillin-binding protein 4 (PBP4), can lead to the overexpression of ampC, thereby contributing to β-lactam resistance. Present work analyzed 208 clinical isolates of P. aeruginosa using whole-genome sequencing (WGS) and detected multiple nonsynonymous single nucleotide polymorphisms (nsSNPs), such as Y264C, G444D, and a double mutation (A394P-T428P). All nsSNPs were predicted to be deleterious by several prediction program. Molecular dynamics (MD) simulations suggested that these substitutions destabilize PBP4, increase structural flexibility, and contribute to the resistance mechanism, which favored their selection. To determine the effective therapeutics against these mutations, molecular docking was conducted with various antibiotics. Cefoperazone exhibited the highest binding affinity (-7.3 kcal/mol) among multiple PBP4 variants. The Molecular dynamics (MD) simulations and Molecular Mechanics Poisson Boltzmann Surface Area calculations (MMPBSA) further confirmed the favorable interactions between cefoperazone and PBP4 variants. In vitro MIC analyses supported these findings, indicating that cefoperazone displayed significant activity against clinical dacB mutants of P. aeruginosa. The study offers structural insight of dacB variants leading to antibiotic resistance and emphasizes the need to prioritize specific antibiotics to address the challenges arising from protein mutations.

Pseudomonas aeruginosa↗