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,297 records · Page 72Linked to original sources
PubMed2005

Anesthetic interaction with ketosteroid isomerase: insights from molecular dynamics simulations.

The nature and the sites of interactions between anesthetic halothane and homodimeric Delta5-3-ketosteroid isomerase (KSI) are characterized by flexible ligand docking and confirmed by 1H-15N NMR. The dynamics consequence of halothane interaction and the implication of the dynamic changes to KSI function are studied by multiple 5-ns molecular dynamics simulations in the presence and absence of halothane. Both docking and MD simulations show that halothane prefer the amphiphilic dimeric interface to the hydrophobic active site of KSI. Halothane occupancy at the dimer interface disrupted the intersubunit hydrogen bonding formed either directly through side chains of polar residues or indirectly through the mediation of the interfacial water molecules. Moreover, in the presence of halothane, the exchange rate of the bound waters with bulk water was increased. Halothane perturbation to the dimer interface affected the overall flexibility of the active site. This action is likely to contribute to the halothane-induced reduction of the KSI activity. The allosteric halothane modulation of the dynamics-function relationship of KSI without direct competition at the enzymatic active sites may be generalized to offer a unifying explanation of anesthetic action on a diverse range of multidomain neuronal proteins that are potentially relevant to clinical general anesthesia.

Anesthetics, Inhalation↗
PubMed2004

Homology modelling and molecular dynamics simulations: comparative studies of human aquaporin-1.

The structures of the mammalian water transport protein Aqp1 and of its bacterial homologue GlpF enables us to test whether homology models can be used to explore relationships between structure, dynamics and function in mammalian transport proteins. Molecular dynamics simulations (totalling almost 40 ns) were performed starting from: the X-ray structure of Aqp1; a homology model of Aqp1 based on the GlpF structure; and intermediate resolution structures of Aqp1 derived from electron microscopy. Comparisons of protein RMSDs vs. time suggest that the homology models are of comparable conformational stability to the X-ray structure, whereas the intermediate resolution structures exhibit significant conformation drift. For simulations based on the X-ray structure and on homology models, the flexibility profile vs. residue number correlates well with the crystallographic B-values for each residue. In the simulations based on intermediate resolution structures, mobility of the highly conserved NPA loops is substantially higher than in the simulations based on the X-ray structure or the homology models. Pore radius profiles remained relatively constant in the X-ray and homology model simulations but showed substantial fluctuations (reflecting the higher NPA loop mobility) in the intermediate resolution simulations. The orientation of the dipoles of water molecules within the pore is of key importance in maintaining low proton permeability through Aqp1. This property seems to be quite robust to the starting model used in the simulation. These simulations suggest that homology models based on bacterial homologues may be used to derive functionally relevant information on the structural dynamics of mammalian transport proteins.

Aquaporin 1↗
PubMed2003

Excited state charge-transfer dynamics study of poplar plastocyanin by ultrafast pump-probe spectroscopy and molecular dynamics simulation.

We have applied ultrafast pump-probe spectroscopy to investigate the excited state dynamics of the blue copper protein poplar plastocyanin, by exciting in the blue side of its 600-nm absorption band. The decay of the charge-transfer excited state occurs exponentially with a time constant of approximately 280 fs and is modulated by well visible oscillations. The Fourier transform of the oscillatory component, besides providing most of the vibrational modes found by conventional resonance Raman, presents additional bands in the low frequency region modes, which are reminiscent of collective motions of biological relevance. Notably, a high frequency mode at approximately 508 cm(-1), whose dynamics are consistent with that of the excited state and already observed for other blue copper proteins, is shown to be present also in poplar plastocyanin. This vibrational mode is reproduced by a molecular dynamics simulation involving the excited state of the copper site.

Computer Simulation↗
PubMed2006

X-ray diffraction "fingerprinting" of DNA structure in solution for quantitative evaluation of molecular dynamics simulation.

Solution state x-ray diffraction fingerprinting is demonstrated as a method for experimentally assessing the accuracy of molecular dynamics (MD) simulations. Fourier transforms of coordinate data from MD simulations are used to produce reciprocal space "fingerprints" of atomic pair distance correlations that are characteristic of the ensemble and are the direct numerical analogues of experimental solution x-ray diffraction (SXD). SXD experiments and MD simulations were carried out to test the ability of experiment and simulation to resolve sequence-dependent modifications in helix conformation for B-form DNA. SXD experiments demonstrated that solution-state poly(AT) and poly(A)-poly(T) duplex DNA sequences exist in ensembles close to canonical B-form and B'-form structures, respectively. In contrast, MD simulations analyzed in terms of SXD fingerprints are shown to deviate from experiment, most significantly for poly(A)-poly(T) duplex DNA. Compared with experiment, MD simulation shortcomings were found to include both mismatches in simulated conformer structures and number population within the ensembles. This work demonstrates an experimental approach for quantitatively evaluating MD simulations and other coordinate models to simulate biopolymer structure in solution and suggests opportunities to use solution diffraction data as experimental benchmarks for developing supramolecular force fields optimized for a range of in situ applications.

DNA↗
PubMed2006

Molecular dynamics simulations of the ligand-binding domain of an N-methyl-D-aspartate receptor.

The mechanism of partial agonism at N-methyl-D-aspartate receptors is an unresolved issue, especially with respect to the role of protein dynamics. We have performed multiple molecular dynamics simulations (7 x 20 ns) to examine the behavior of the ligand-binding core of the NR1 subunit with a series of ligands. Our results show that water plays an important role in stabilizing different conformations of the core and how a closed cleft conformation of the protein might be stabilized in the absence of ligands. In the case of ligand-bound simulations with both full and partial agonists, we observed that ligands within the binding cleft may undergo distinct conformational changes, without grossly influencing the degree of cleft closure within the ligand-binding domain. In agreement with recently published crystallographic data, we also observe similar changes in backbone torsions corresponding to the hinge region between the two lobes for the partial agonist, D-cycloserine. This observation rationalizes the classification of D-cycloserine as a partial agonist and should provide a basis with which to predict partial agonism in this class of receptor by analyzing the behavior of these torsions with other potential ligands.

Animals↗
PubMed2004

Characterization of molecular structures and properties of polyurethanes using molecular dynamics simulations.

Thermotropic polyurethanes with mesogenic groups in side chains were prepared from two diisocyanates and four diols with stoichiometric ratios of reactive isocyanate (NCO) and hydroxy (OH) groups. Their thermal behavior was determined by differential scanning calorimetry. The effect of structure modifications of the diisocyanates and diols, in particular changes in the mesogen, were investigated. Introduction of mesogenic segments into the polymers suppresses the ordering. Stiff end substituents (phenyl and alkoxy groups) of the mesogens stabilize the mesophases to such an extent that the negative influence of long polymer chains is compensated and the liquid-crystalline properties are recovered. All-atom molecular dynamics simulations in the Cerius2 modeling environment were carried out to characterize the structures of the polymers. Analysis of the dynamic trajectories at 20, 100, 120 and 170 degrees C revealed changes in conformation of macromolecules, which correlate with DSC measurements.

Calorimetry↗
PubMed2006

Structural instability of the prion protein upon M205S/R mutations revealed by molecular dynamics simulations.

The point mutations M205S and M205R have been demonstrated to severely disturb the folding and maturation process of the cellular prion protein (PrP(C)). These disturbances have been interpreted as consequences of mutation-induced structural changes in PrP, which are suggested to involve helix 1 and its attachment to helix 3, because the mutated residue M205 of helix 3 is located at the interface of these two helices. Furthermore, current models of the prion protein scrapie (PrP(Sc)), which is the pathogenic isoform of PrP(C) in prion diseases, imply that helix 1 disappears during refolding of PrP(C) into PrP(Sc). Based on molecular-dynamics simulations of wild-type and mutant PrP(C) in aqueous solution, we show here that the native PrP(C) structure becomes strongly distorted within a few nanoseconds, once the point mutations M205S and M205R have been applied. In the case of M205R, this distortion is characterized by a motion of helix 1 away from the hydrophobic core into the aqueous environment and a subsequent structural decay. Together with experimental evidence on model peptides, this decay suggests that the hydrophobic attachment of helix 1 to helix 3 at M205 is required for its correct folding into its stable native structure.

Computer Simulation↗
PubMed2003

Human factor Xa bound amidine inhibitor conformation by double rotational-echo double resonance nuclear magnetic resonance and molecular dynamics simulations.

Double rotational-echo double resonance (double REDOR) NMR was used to investigate the conformation of a (13)C-, (15)N-, and (19)F-labeled inhibitor (Berlex Biosciences compound no. ZK-806299) bound to human factor Xa. Conformationally dependent carbon-fluorine dipolar couplings were measured by (13)C[(19)F] REDOR. Natural abundance carbon signals in the full-echo spectra were removed by (13)C[(15)N] REDOR. Major and minor binding modes were suggested by the NMR data, but only the former had adequate signal to noise for distance determinations. Molecular dynamics simulations restrained by double-REDOR-determined intramolecular (13)C-(19)F distances revealed two models for the dominant binding mode that are consistent with the NMR data. We conclude that ZK-806299 binds similarly to both FXa. Moreover, it appears to bind to FXa in a fashion previously demonstrated for ZK-807834, a more selective FXa inhibitor.

Amidines↗
PubMed2006

Molecular dynamics simulations of arachidonic acid-derived pentadienyl radical intermediate complexes with COX-1 and COX-2: insights into oxygenation regio- and stereoselectivity.

The two cyclooxygenase enzymes, COX-1 and COX-2, are responsible for the committed step in prostaglandin biosynthesis and are the targets of the nonsteroidal antiinflammatory drugs aspirin and ibuprofen and the COX-2 selective inhibitors, Celebrex, Vioxx, and Bextra. The enzymes are remarkable in that they catalyze two dioxygenations and two cyclizations of the native substrate, arachidonic acid, with near absolute regio- and stereoselectivity. Several theories have been advanced to explain the nature of enzymatic control over this series of reactions, including suggestions of steric shielding and oxygen channeling. As proposed here, selective radical trapping and spin localization in the substrate-derived pentadienyl radical intermediate can also be envisioned. Herein we describe the results of explicit, 10 ns molecular dynamics simulations of both COX-1 and COX-2 with the substrate-derived pentadienyl radical intermediate bound in the active site. The enzymes' influence on the conformation of the pentadienyl radical was investigated, along with the accessible space above and below the radical plane and the width of several channels to the active site that could function as access routes for molecular oxygen. Additional simulations demonstrated the extent of molecular oxygen mobility within the active site. The results suggest that spin localization is unlikely to play a role in enzymatic control of this reaction. Instead, a combination of oxygen channeling, steric shielding, and selective radical trapping appears to be responsible. This work adds a dynamic perspective to the strong foundation of static structural data available for these enzymes.

Arachidonic Acid↗
PubMed2003

Molecular dynamics simulations on the first two helices of Vpu from HIV-1.

Vpu is an 81 amino acid protein of HIV-1 with two phosphorylation sites. It consists of a short N-terminal end traversing the bilayer and a longer cytoplasmic part. The dual functional role of Vpu is attributed to these topological distinct regions of the protein. The first 52 amino acids of Vpu (HV1H2) have been simulated, which are thought to be embedded in a fully hydrated lipid bilayer and to consist of a transmembrane helix (helix-1) connected via a flexible linker region, including a Glu-Tyr-Arg (EYR) motif, with a second helix (helix-2) residing with its helix long axis on the bilayer surface. Repeated molecular dynamics simulations show that Glu-28 is involved in salt bridge formation with Lys-31 and Arg-34 establishing a kink between the two helices. Helix-2 remains in a helical conformation indicating its stability and function as a "peptide float," separating helix-1 from the rest of the protein. This leads to the conclusion that Vpu consists of three functional modules: helix-1, helix-2, and the remaining residues toward the C-terminal end.

Amino Acid Sequence↗
PubMed1996

260 ps molecular dynamics simulation of substance P with hydrated dimyristoyl phosphatidyl choline bilayer.

We present here results on 260 pico seconds (ps) molecular dynamics (MD) simulation of substance P (SP) in hydrated bilayer of dimyristoyl phosphatidyl choline (DMPC) (39 molecules of DMPC with 776 water molecules). 260 ps MD simulation has been carried out in 0.001 ps time interval with united atom force field, using AMBER 4.0 package. Non bonded pair list was updated every 20 cycles using 12.5 Angstrom cut off distance. Analysis of MD data is done using our package ANALMD. The obtained models are presented using graphics package RASMOL. All simulations, analysis of MD data and graphics is done on INDIGO-2, R-4400 extreme graphics work station. Our results show no systematic change in order parameter, but reduction in transfraction of the chain torsional angles, compared to our earlier results on MD simulation on hydrated DMPC bilayer without SP. C-terminal and central peptide residues adopt partial helical conformation. Helix type as classified on the basis of H-bonds is between alpha and 3(10). The peptide backbone shows flexibility during heating runs. Later, it is stabilized and there was not much change in the spatial position of the backbone. Lipid matrix serves the role of immobilization of the peptide backbone in a preferred conformation.

Amino Acid Sequence↗
PubMed1994

Rational design of a three-heptad coiled-coil protein and comparison by molecular dynamics simulation with the GCN4 coiled coil: presence of interior three-center hydrogen bonds.

alpha-Helical coiled coils have a 7-residue repeating pattern (abcdefg) where a and d are usually hydrophobic. We have designed a 2-stranded 44-residue coiled-coil protein (P44) consisting of 2 22-residue alpha-helices linked by 2 terminal disulfide groups to test whether the disulfide bridges could stabilize a 3-heptad coiled coil. P44 should be stabilized by intrahelical hydrogen bonds, interhelical disulfide and salt bridges, and interior hydrophobic interactions. A computer model of P44 was built and its stability was studied by molecular dynamics simulation with explicit water. This doubly crosslinked 3-heptad coiled coil did not unfold during a 300-ps simulation with explicit water. This doubly crosslinked 3-heptad coiled coil did not unfold during a 300-ps simulation. But reduced P44 with 4 thiol groups did unfold. For comparison, the 62-residue crystal structure of the 4-heptad coiled coil of transcription activator GCN4 did not unfold during a 300-ps simulation. Thus P44 may be a stable folded protein in aqueous solution. These simulations revealed the presence of 2 local hydrogen bond networks involving intra-helical 3-center hydrogen bonds in the hydrophobic interior of the coiled coils of GCN4 and P44. The NH hydrogen at d makes a 3-center hydrogen bond whose major component is to the i - 4 C = O oxygen at g and minor component is to the solvent-inaccessible i - 3 C = O oxygen at a. Likewise, the NH hydrogen at g makes a 3-center hydrogen bond with the i - 4 C = O oxygen at c and the buried i - 3 C = O oxygen at d.

Amino Acid Sequence↗
PubMedDate not supplied

Homology model and molecular dynamics simulation of carp ovum cystatin.

In this study, a homology model of carp ovum cystatin was constructed based on the crystal structure of chicken egg white cystatin. The results of amino acid sequence alignment indicate that these two proteins exhibit 36.11% of sequence identity. The resultant homology model reveals that carp ovum cystatin shares similar folds as chicken egg white cystatin, particularly in the conserved regions of Q48-V49-G52 and P98-W99 and the locations of two disulfide bonds, C67-C76 and C90-C110. However, the results of 1 ns molecular dynamics simulations show that carp ovum cystatin exhibits less structural integrity than chicken egg white cystatin in explicit water at 300 K. The relatively hydrophilic Met62 of carp ovum cystatin, corresponding to the hydrophobic Leu68 of human cystatin C and Ile66 of chicken egg white cystatin, may destabilize the hydrophobic core and form a dimeric structure more easily through domain swapping. A total of 16 positively charged residues are equally distributed on the surface of carp ovum cystatin, resulting in agglutination with the negatively charged spermatozoa via electrostatic interaction. Thus, carp ovum cystatin is considered to be important in preventing carp eggs from polyspermy.

Agglutination↗
PubMed2004

Rapid protein-ligand docking using soft modes from molecular dynamics simulations to account for protein deformability: binding of FK506 to FKBP.

Most current docking methods to identify possible ligands and putative binding sites on a receptor molecule assume a rigid receptor structure to allow virtual screening of large ligand databases. However, binding of a ligand can lead to changes in the receptor protein conformation that are sterically necessary to accommodate a bound ligand. An approach is presented that allows relaxation of the protein conformation in precalculated soft flexible degrees of freedom during ligand-receptor docking. For the immunosuppressant FK506-binding protein FKBP, the soft flexible modes are extracted as principal components of motion from a molecular dynamics simulation. A simple penalty function for deformations in the soft flexible mode is used to limit receptor protein deformations during docking that avoids a costly recalculation of the receptor energy by summing over all receptor atom pairs at each step. Rigid docking of the FK506 ligand binding to an unbound FKBP conformation failed to identify a geometry close to experiment as favorable binding site. In contrast, inclusion of the flexible soft modes during systematic docking runs selected a binding geometry close to experiment as lowest energy conformation. This has been achieved at a modest increase of computational cost compared to rigid docking. The approach could provide a computationally efficient way to approximately account for receptor flexibility during docking of large numbers of putative ligands and putative docking geometries.

Binding Sites↗
PubMed2006

Assessing the acid-base and conformational properties of histidine residues in human prion protein (125-228) by means of pK(a) calculations and molecular dynamics simulations.

A thorough study of the acid-base behavior of the four histidines and the other titratable residues of the structured domain of human prion protein (125-228) is presented. By using multi-tautomer electrostatic calculations, average titration curves have been built for all titratable residues, using the whole bundles of NMR structures determined at pH 4.5 and 7.0. According to our results, (1) only histidine residues are likely to be involved in the first steps of the pH-driven conformational transition of prion protein; (2) the pK(a)'s of His140 and His177 are approximately 7.0, whereas those of His155 and His187 are < 5.5. 10-ns long molecular dynamics simulations have been performed on five different models, corresponding to the most significant combinations of histidine protonation states. A critical comparison between the available NMR structures and our computational results (1) confirms that His155 and His187 are the residues whose protonation is involved in the conformational rearrangement of huPrP in mildly acidic condition, and (2) shows how their protonation leads to the destructuration of the C-terminal part of HB and to the loss of the last turn of HA that represent the crucial microscopic steps of the rearrangement.

Amino Acid Sequence↗
PubMed1995

Solution conformation of mu-selective dermorphin and delta-selective deltorphin-I in phospholipid micelles, studied by NMR spectroscopy and molecular dynamics simulations.

Complete proton resonance assignments of the naturally occurring mu-selective dermorphin (H-Tyr-D-Ala-Phe-Gly-Tyr-Pro-Ser-NH2) and delta-selective deltorphin-I (H-Tyr-D-Ala-Phe-Asp-Val-Val-Gly-NH2) were carried out by two-dimensional 1H-NMR techniques to investigate the conformational features in the membrane-mimetic micelles of perdeuterated dodecylphosphocholine. Fifty possible three-dimensional structures for respective peptides were generated by means of distance geometry calculations, all of which satisfy the proton-proton distances derived from NOE measurements within the allowable range, and 25 of them were subjected to the molecular dynamics simulations for 10 ps, in which the NOE distances were included as the energetic constraints. Although conformers simulated for dermorphin showed relatively large conformational variations because of the limited NOE data, most of them were characterized as an entirely folded structure bent at the Gly4 residue, where each of the N- and C-terminal tetrapeptides took an extended conformation. On the other hand, most conformations of deltorphin-I showed the common feature that the N-terminal Tyr-D-Ala-Phe-Asp and C-terminal Val-Val-Gly-NH2 sequences took respective folded conformations, and these were almost at right angles on the border of the Asp-Val sequence. These conformational characteristics are discussed in terms of the possible relationship with the mu/delta-opioid receptor selectivity.

Amino Acid Sequence↗
PubMed2002

Molecular dynamics simulations on the free and complexed N-terminal SH2 domain of SHP-2.

Signal transduction events are often mediated by small protein domains such as SH2 (Src homology 2) domains that recognize phosphotyrosines (pY) and flanking sequences. In case of the SHP-2 receptor tyrosine phosphatase an N-terminal SH2 domain binds and inactivates the phosphatase (PTP) domain. The pY-peptide-binding site on the N-terminal SH2 domain does not overlap with the PTP binding region. Nevertheless, pY-peptide binding causes domain dissociation and phosphatase activation. Comparative multi-nanosecond molecular dynamics simulations on the N-SH2 domain in ligand-bound and free states have been performed to study the allosteric mechanism that leads to domain dissociation upon pY-peptide binding. Significant ligand-dependent differences in the conformational flexibility of regions that are involved in SH2-PTP domain association have been observed. The results support a mechanism of signal transduction where SH2-peptide binding modulates the domain flexibility and reduces its capacity to fit into the entrance of the PTP catalytic domain of SHP-2.

Allosteric Regulation↗
PubMed1999

The mechanism of channel formation by alamethicin as viewed by molecular dynamics simulations.

Alamethicin is a 20-residue channel-forming peptide that forms a stable amphipathic alpha-helix in membrane and membrane-mimetic environments. This helix contains a kink induced by a central Gly-X-X-Pro sequence motif. Alamethicin channels are activated by a cis positive transbilayer voltage. Channel activation is suggested to correspond to voltage-induced insertion of alamethicin helices in the bilayer. Alamethicin forms multi-conductance channels in lipid bilayers. These channels are formed by parallel bundles of transmembrane helices surrounding a central pore. A change in the number of helices per bundle switches the single channel conductance level. Molecular dynamics simulations of alamethicin in a number of different environments have been used to explore its channel-forming properties. These simulations include: (i) alamethicin in solution in water and in methanol; (ii) a single alamethicin helix at the surface of a phosphatidylcholine bilayer; (iii) single alamethicin helices spanning a phosphatidylcholine bilayer; and (iv) channels formed by bundles of 5, 6, 7 or 8 alamethicin helices spanning a phosphatidylcholine bilayer. The total simulation time is c. 30 ns. Thus, these simulations provide a set of dynamic snapshots of a possible mechanism of channel formation by this peptide.

Alamethicin↗