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 1,171 records · Page 65Linked to original sources
PubMed2006

Modeling evolution of hydrogen bonding and stabilization of transition states in the process of cocaine hydrolysis catalyzed by human butyrylcholinesterase.

Molecular dynamics (MD) simulations and quantum mechanical/molecular mechanical (QM/MM) calculations were performed on the prereactive enzyme-substrate complex, transition states, intermediates, and product involved in the process of human butyrylcholinesterase (BChE)-catalyzed hydrolysis of (-)-cocaine. The computational results consistently reveal a unique role of the oxyanion hole (consisting of G116, G117, and A199) in BChE-catalyzed hydrolysis of cocaine, compared to acetylcholinesterase (AChE)-catalyzed hydrolysis of acetylcholine. During BChE-catalyzed hydrolysis of cocaine, only G117 has a hydrogen bond with the carbonyl oxygen (O31) of the cocaine benzoyl ester in the prereactive BChE-cocaine complex, and the NH groups of G117 and A199 are hydrogen-bonded with O31 of cocaine in all of the transition states and intermediates. Surprisingly, the NH hydrogen of G116 forms an unexpected hydrogen bond with the carboxyl group of E197 side chain and, therefore, is not available to form a hydrogen bond with O31 of cocaine in the acylation. The NH hydrogen of G116 is only partially available to form a weak hydrogen bond with O31 of cocaine in some structures involved in the deacylation. The change of the estimated hydrogen-bonding energy between the oxyanion hole and O31 of cocaine during the reaction process demonstrates how the protein environment can affect the energy barrier for each step of the BChE-catalyzed hydrolysis of cocaine. These insights concerning the effects of the oxyanion hole on the energy barriers provide valuable clues on how to rationally design BChE mutants with a higher catalytic activity for the hydrolysis of (-)-cocaine.

Butyrylcholinesterase↗
PubMed2006

Molecular dynamics study of the stabilities of consensus designed ankyrin repeat proteins.

Two designed ankyrin repeat (AR) proteins (E3_5 and E3_19) are high homologous (with about 87% sequence identity) and their crystal structures have a Calpha atom-positional root-mean-square difference of about 0.14 nm. However, it was found that E3_5 is considerably more stable than E3_19 in guanidinium hydrochloride and thermal denaturation experiments. With the goal of providing insights into the various factors contributing to the stabilities of the designed AR proteins and suggesting possible mutations to enhance their stabilities, homology modeling and molecular dynamics (MD) simulations with explicit solvent have been performed. Because the crystal structure of E3_19 was solved later than that of E3_5, a homology model of E3_19 based on the crystal structure of E3_5 was also used in the simulations. E3_5 shows a very stable trajectory in both crystal and solution simulations. In contrast, the C-terminal repeat of E3_19 unfolds in the simulations starting from either the modeled structure or the crystal structure, although it has a sequence identical to that of E3_5. A continuum electrostatic model was used to estimate the effect of single mutations on protein stability and to study the interaction between the internal ARs and the C-terminal capping AR. Mutations involving charged residues were found to have large effects on stability. Due to the difference in charge distribution in the internal ARs of E3_19 and E3_5, their interaction with the C-terminal capping AR is less favorable in E3_19. The simulation trajectories suggest that the stability of the designed AR proteins can be increased by optimizing the electrostatic interactions within and between the different repeats.

Ankyrin Repeat↗
PubMed2006

Homology modeling and molecular dynamics study of West Nile virus NS3 protease: a molecular basis for the catalytic activity increased by the NS2B cofactor.

The West Nile virus (WNV) NS3 serine protease, which plays an important role in assembly of infective virion, is an attractive target for anti-WNV drug development. Cofactors NS2B and NS4A increase the catalytic activity of NS3 in dengue virus and Hepatitis C virus, respectively. Recent studies on the WNV-NS3 characterize the catalytically active form of NS3 by tethering the 40-residue cofactor NS2B. It is suggested that NS2B is essential for the NS3 activity in WNV, while there is no information of the WNV-NS3-related crystal structure. To understand the role of NS2B/substrate in the NS3 catalytic activity, we built a series of models: WNV-NS3 and WNV-NS3-NS2B and WNV-NS3-NS2B-substrate using homology modeling and molecular modeling techniques. Molecular dynamics (MD) simulations were performed for 2.75 ns on each model, to investigate the structural stabilization and catalytic triad motion of the WNV NS3 protease with and without NS2B/substrate. The simulations show that the NS3 rearrangement occurs upon the NS2B binding, resulting in the stable D75-OD1...H51-NH hydrogen bonding. After the substrate binds to the NS3-NS2B active site, the NS3 protease becomes more stable, and the catalytic triad is formed. These results provide a structural basis for the activation and stabilization of the enzyme by its cofactor and substrate.

Amino Acid Sequence↗
PubMed1993

Essential dynamics of proteins.

Analysis of extended molecular dynamics (MD) simulations of lysozyme in vacuo and in aqueous solution reveals that it is possible to separate the configurational space into two subspaces: (1) an "essential" subspace containing only a few degrees of freedom in which anharmonic motion occurs that comprises most of the positional fluctuations; and (2) the remaining space in which the motion has a narrow Gaussian distribution and which can be considered as "physically constrained." If overall translation and rotation are eliminated, the two spaces can be constructed by a simple linear transformation in Cartesian coordinate space, which remains valid over several hundred picoseconds. The transformation follows from the covariance matrix of the positional deviations. The essential degrees of freedom seem to describe motions which are relevant for the function of the protein, while the physically constrained subspace merely describes irrelevant local fluctuations. The near-constraint behavior of the latter subspace allows the separation of equations of motion and promises the possibility of investigating independently the essential space and performing dynamic simulations only in this reduced space.

Animals↗
PubMed1994

The enzymatic mechanism of carboxypeptidase: a molecular dynamics study.

An MD simulation of the system carboxypeptidase A (CPA) with the tetrapeptide Val-Leu-Phe-Phe has been performed in order to learn about the substrate disposition just prior to nucleophilic attack. We have explored the model in which the substrate does not substitute the zinc-coordinated water (the "water" mechanism). The simulations do suggest as feasible that the Zn-OH2 group performs a nucleophilic attack on the Phe-Phe peptidic bond. We have also investigated the model in which the carbonyl oxygen displaces the zinc-coordinated water. In this case the substrate and Glu-270 orient themselves to allow an anhydride intermediate during the peptidic bond cleavage (the "anhydride" mechanism). Based on the results of the simulations, both "water" and "anhydride" mechanisms are structurally feasible, although the former model seems more probable on chemical grounds.

Amino Acid Sequence↗
PubMed2005

Conformational study on glycosylated asparagine-oligopeptides by NMR spectroscopy and molecular dynamics calculations.

The conformational properties of the homo oligomers of increasing chain length Boc-(Asn)(n)-NHMe (n = 2, 4, 5), (GlcNAc-beta-Asn)(n)-NHMe (n = 2, 4, 5, 8) and Boc-[GlcNAc(Ac)(3)-beta-Asn](n)-NHMe (n = 2, 4, 5) were studied by using NOE experiments and molecular dynamic calculations (MD). Sequential NOEs and medium range NOEs, including (i,i+2) interactions, were detected by ROESY experiments and quantified. The calculated inter-proton distances are longer than those characteristic of beta-turn secondary structures. Owing to the large conformational motions expected for linear peptides, MD simulations were performed without NMR constraints, with explicit water and by applying different treatments of the electrostatic interactions. In agreement with the NOE results, the simulations showed, for all peptides, the presence of both folded and unfolded structures. The existence of significant populations of beta-turn structures can be excluded for all the examined compounds, but two families of structures were more often recognized. The first one with sinusoidal or S-shaped forms, and another family of large turns together with some more extended conformations. Only the glycosylated pentapeptide shows in vacuo a large amount of structures with helical shaped form. The results achieved in water and in DMSO are compared and discussed, together with the effect of the glycosylation.

Asparagine↗
PubMed2006

Conformational analysis of a potent SSTR3-selective somatostatin analogue by NMR in water solution.

The three-dimensional structure of a potent SSTR3-selective analogue of somatostatin, cyclo(3-14)H-Cys(3)-Phe(6)-Tyr(7)-D-Agl(8)(N(beta) Me, 2-naphthoyl)-Lys(9)-Thr(10)-Phe(11)-Cys(14)-OH (des-AA(1, 2, 4, 5, 12, 13)[Tyr(7), D-Agl(8)(N(beta) Me, 2-naphthoyl)]-SRIF) (peptide 1) has been determined by (1)H NMR in water and molecular dynamics (MD) simulations. The peptide exists in two conformational isomers differing mainly by the cis/trans isomerization of the side chain in residue 8. The structure of 1 is compared with the consensus structural motifs of other somatostatin analogues that bind predominantly to SSTR1, SSTR2/SSTR5 and SSTR4 receptors, and to the 3D structure of a non-selective SRIF analogue, cyclo(3-14)H-Cys(3)-Phe(6)-Tyr(7)-D-2Nal(8)-Lys(9)-Thr(10)-Phe(11)-Cys(14)-OH (des-AA(1, 2, 4, 5, 12, 13)[Tyr(7), D-2Nal(8)]-SRIF) (peptide 2). The structural determinant factors that could explain selectivity of peptide 1 for SSTR3 receptors are discussed.

Amino Acid Sequence↗
PubMed2007

Bombesin analogs containing alpha-amino-isobutyric acid with potent anticancer activity.

Six octapeptide bombesin (BN) analogs were synthesized by substituting alpha-aminoisobutyric acid (Aib), in place of Ala9 or Gly11, or both, in the [D-Phe6, desMet14]-BN (6-14) sequence: D-Phe6-Gln7-Trp8-Ala9-Val10-Gly11-His12-Leu13-NH2 (P0). Additionally, Leu13 was replaced with isoleucine in two analogs and one of the analogs was butanoylated at the N-terminus. The antiproliferative activity of the analogs was tested in vitro on human pancreatic (MiaPaCa-2) and colon cancer (SW620, HT29 and PTC) cell lines using the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. The analogs demonstrated anticancer activity in the above cell lines at concentrations ranging from 0.01 nM to 1 microM. One of the analogs, P6, was evaluated for in vivo tumor regression in a xenograft model of human primary colon cancer in athymic nude mice and was found to cause significant reduction in tumor volume. NMR and molecular dynamics (MD) simulation studies for this analog revealed the presence of a mixed 3(10)/alpha-helical structure. This study demonstrates that the designed BN analogs retain their anticancer activity after the incorporation of the constrained amino acid, Aib, and are potential molecules for future use in cancer therapy and drug targeting.

Amino Acid Sequence↗
PubMed2006

Molecular dynamics study of amyloid formation of two Abl-SH3 domain peptides.

Molecular dynamics (MD) simulations were carried out for two-strand and ten-strand beta-sheets constructed from two peptides corresponding to the diverging turn of two homologous Abl-SH3 domains, DLSFMKGE (MK; from Drosophila) and DLSFKKGE (KK; from man), in explicit water at the temperatures of 30, 170/190 and 300 K. It was found that the 2 x MK beta-sheet is more stable than the 2 x KK beta-sheet, and that the 10 x MK beta-sheet is more stable than the 10 x KK beta-sheet; this suggests that the MK systems are fibril-creating and the KK systems are not. These results might explain why most SH3 domains possess two conserved basic residues at the diverging turn, which may act as gatekeepers in order to avoid aggregation.

Amino Acid Sequence↗
PubMed2000

Structure and thermodynamics of RNA-protein binding: using molecular dynamics and free energy analyses to calculate the free energies of binding and conformational change.

An adaptive binding mechanism, requiring large conformational rearrangements, occurs commonly with many RNA-protein associations. To explore this process of reorganization, we have investigated the conformational change upon spliceosomal U1A-RNA binding with molecular dynamics (MD) simulations and free energy analyses. We computed the energetic cost of conformational change in U1A-hairpin and U1A-internal loop binding using a hybrid of molecular mechanics and continuum solvent methods. Encouragingly, in all four free energy comparisons (two slightly different proteins, two different RNAs), the free macromolecule was more stable than the bound form by the physically reasonable value of approximately 10 kcal/mol. We calculated the absolute binding free energies for both complexes to be in the same range as that found experimentally.

Amino Acid Sequence↗
PubMed1996

Hydroxyl and water molecule orientations in trypsin: comparison to molecular dynamic structures.

A comparison is presented of experimentally observed hydroxyl and water hydrogens in trypsin determined from neutron density maps with the results of a 140ps molecular dynamics (MD) simulation. Experimental determination of hydrogen and deuterium atom positions in molecules as large as proteins is a unique capability of neutron diffraction. The comparison addresses the degree to which a standard force-field approach can adequately describe the local electrostatic and van der Waals forces that determine the orientations of these hydrogens. The molecular dynamics simulation, based on the all-atom AMBER force-field, allowed free rotation of all hydroxyl groups and movement of water molecules making up a bath surrounding the protein. The neutron densities, derived from 2.1A D2O-H2O difference Fourier maps, provide a database of 27 well-ordered hydroxyl hydrogens. Virtually all of the simulated hydroxyl orientations are within a standard deviation of the experimentally-observed positions, including several examples in which both the simulation and the neutron density indicate that a hydroxyl group is shifted from a 'standard' rotamer. For the most highly ordered water molecules, the hydrogen distributions calculated from the trajectory were in good agreement with neutron density; simulated water molecules that displayed multiple hydrogen bonding networks had correspondingly broadened neutron density profiles. This comparison was facilitated by development of a method to construct a pseudo 2A density map based on the hydrogen atom distributions from the simulation. The degree of internal water molecules is shown to result primarily from the electrostatic environment surrounding that water molecule as opposed to the cavity size available to the molecule. A method is presented for comparing the discrete observations sampled in a dynamics trajectory with the time-averaged data obtained from X-ray or neutron diffraction studies. This method is particularly useful for statically-disordered water molecules, in which the average location assigned from a trajectory may represent a site of relatively low occupancy.

Crystallography↗
PubMed1994

Pharmacophore refinement of gpIIb/IIIa antagonists based on comparative studies of antiadhesive cyclic and acyclic RGD peptides.

Structurally guided design approaches to low-molecular-weight platelet aggregation antagonists addressing the platelet-associated heterodimeric cell surface receptor gpIIb/IIIa rely on comparative studies of an ensemble of conformationally and biologically characterized compounds, since no high-resolution structure of the receptor system is available. We report a classical indirect and comparative pharmacophore refinement approach based on a series of small cyclic Arg-Gly-Asp (RGD) peptides as gpIIb/IIIa-fibrinogen interaction antagonists. These peptides have previously been investigated as potent and selective tumor cell adhesion inhibitors. The definition of geometrical descriptors classifying the RGD peptide conformations and their subsequent analysis over selected RGD- and RXD-containing protein structures allows for a correlation of distinct structural features for platelet aggregation inhibition. An almost parallel alignment of the Arg and Asp side chains was identified by a vector analysis as being present in all active cyclic hexa- and pentapeptides. This orientation is induced mainly by the constraint of backbone cyclization and is not of any covalent tripeptide-inherent origin, which was rationalized by a 500 ps high-energy MD simulation of a sequentially related linear model peptide. The incorporation of the recognition tripeptide Arg-Gly-Asp into the cyclic peptide templates acted as a filter mechanism, restricting the otherwise free torsional relation of both side chains to a parallel orientation. Based on the derived results, several detailed features of the receptor binding site could be deduced in terms of receptor complementarity. These findings should govern the design of next-generation compounds with enhanced activities. Furthermore, the complementary stereochemical characteristics of the substrate can be used as boundary conditions for pseudoreceptor modelling studies that are capable of reconstructing a hypothetical binding pocket, qualitatively resembling the steric and electronic demands of gpIIb/IIIa. It is interesting to note that these features provide clear differentiation to requirements for inhibition of alpha v beta 3 substrate binding. This can account for the extremely high selectivity and activity of some of our constrained peptides for either the alpha IIb beta 3 or the alpha v beta 3 receptor.

Amino Acid Sequence↗
PubMed1987

Energy minimization and molecular dynamics studies of Asn-102 elastase.

Four isomeric forms of the Asn-102 PPE (D102N mutant according to the emerging protocol, [Knowles, Science, 236 (1987) 1252-1258]) have been investigated using energy minimization (EM) and molecular dynamics (MD) techniques. MD simulation data for 175 ps are reported for each form (in total 700 ps for about 2500 atoms). The His-57 N epsilon-protonated forms are calculated to be more stable than the N delta-protonated ones. The active site region of the most stable form is very similar to that found in the D102N rat trypsin enzyme [Craik et al., Science, 237 (1987) 909-913]. Conformations of the active sites and their hydrogen bond patterns are presented for each of these forms and are compared with the structure of the native enzyme active site. The pH dependent activity of the D102N derivative is discussed.

Animals↗
PubMed2024

Early Events in β2AR Dimer Dynamics Mediated by Activation-Related Microswitches.

G-Protein-Coupled Receptors (GPCRs) make up around 3-4% of the human genome and are the targets of one-third of FDA-approved drugs. GPCRs typically exist as monomers but also aggregate to form higher-order oligomers, including dimers. β2AR, a pharmacologically relevant GPCR, is known to be targeted for the treatment of asthma and cardiovascular diseases. The activation of β2AR at the dimer level remains under-explored. In the current study, molecular dynamics (MD) simulations have been performed to understand activation-related structural changes in β2AR at the dimer level. The transition from inactive to active and vice versa has been studied by starting the simulations in the apo, agonist-bound, and inverse agonist-bound β2AR dimers for PDB ID: 2RH1 and PDB ID: 3P0G, respectively. A cumulative total of around 21-μs simulations were performed. Residue-based distances, RMSD, and PCA calculations suggested that either of the one monomer attained activation-related features for the apo and agonist-bound β2AR dimers. The TM5 and TM6 helices within the two monomers were observed to be in significant variation in all the simulations. TM5 bulge and proximity of TM2 and TM7 helices may be contributing to one of the early events in activation. The dimeric interface between TM1 and helix 8 were observed to be well maintained in the apo and agonist-bound simulations. The presence of inverse agonists favored inactive features in both the monomers. These key features of activation known for monomers were observed to have an impact on β2AR dimers, thereby providing an insight into the oligomerization mechanism of GPCRs.

Receptors, Adrenergic, beta-2↗
PubMed2002

Molecular dynamics and IR spectroscopy in investigation of phase transitions in molecular crystal 4,4'-bis(11-hydroxy-1-undecyloxy)biphenyl.

Molecular dynamics (MD) simulations combined with temperature-dependent IR spectroscopic measurements were used to study phase transitions in molecular crystals of the mesogenic diol 4,4'-bis(11-hydroxy-1-undecyloxy)biphenyl. DSC measurements revealed four phase transitions in this molecular crystal at approximately 327.1 K, 389.8 K, 419.1 K and 431.9 K. Analysis of the dynamic trajectories at temperatures of 300 K, 360 K, 400 K and 480 K revealed changes in conformation of the mesogenic diol molecules and consequently changes in crystal packing and crystal structure in the temperature range 300-480 K and enabled us to understand the mechanism of the phase transitions.

Biphenyl Compounds↗
PubMed2003

Molecular recognition between 4a S/R-galanthamine diastereoisomers and alpha-cyclodextrin.

Molecular recognition between 4a S/R-galanthamine diastereoisomers (1: 4a S-galanthamine;2: 4a R-galanthamine) and alpha-cyclodextrin (alpha-CD) were studied by use of docking and molecular dynamics (MD) simulation approaches. The binding energy of constructed 2...alpha-CD complexes is approximately 17 kcal mol(-1) lower than that of 1...alpha-CD, implying a stronger binding ability of 2 with alpha-CD than that of1. The theoretical modeling result is consistent with our previous CZE result, which demonstrated that alpha-CD is an efficient chiral additive for separating 1 and 2. The modeling result also indicates that both hydrophobic interaction and H-bond force may work as major factors for molecular recognition between the galanthamine diastereoisomers and alpha-CD.

Computer Simulation↗
PubMed2005

Dynamic conformational states of DNA containing T.T or BrdU.T mispaired bases: wobble H-bond pairing versus cross-strand inter-atomic contacts.

The dynamic structure of 11-mer DNA duplexes of different sequences with or without homopyrimidine (T.T, or BrdU.T) mismatches was studied by molecular dynamics (MD) simulations on a time scale from 200 ps to 1 ns. The conformational analysis suggests that in mismatched duplexes the formation of classical T.T wobble H-bonding pairing is nearest-neighbor sequence-dependent and, in most cases, three-centered H-bonds and numerous alternative close cross-strand interatomic contacts exist. Thus, in duplex W1, where the central triplet is 5'd(CTA).d(TTG), two wobble conformations W upward arrow (alphabeta) and W downward arrow (betaalpha) are formed and exchange rapidly at 300 K. In contrast, when the central triplet is 5'd(TTT).d(ATA) (W2 duplex) wobble conformations are rarely observed at 300 K, and the T.T mispair most often adopts a "twisted" conformation with one largely persistent normal H-bond, plus a stable cross-strand contact involving a T flanking base. However, at elevated temperature (400 K) the same W2 duplex shows frequent exchange between the two classical wobble conformations (alphabeta<-->betaalpha), as is in the case when the central triplet is 5'd(TBrdUT).d(ATA) (W3 duplex at 300 K). It is suggested that in the W2 sequence, restrictions due to thymine-methyl/pi interactions prevent the formation of wobble pairing and thermal activation energy, and/or the chemical replacement of T by BrdU are required in order for the T(BrdU).T mismatch to adopt and exchange between wobble conformations. The specific short and/or long-lived (double/triple) cross-strand dynamic interactions in W1, W2 and W3 duplexes are throughout characterized. These frequent atomic encounters exemplify possible inter-strand charge transfer pathways in the studied DNA molecules.

Base Composition↗
PubMed2006

Effect of HPr phosphorylation on structure, dynamics, and interactions in the course of transcriptional control.

The serine46-phosphorylated form of the bacterial protein HPr fulfils an essential function in carbon catabolite repression (CCR). Using molecular dynamics (MD) we studied the effect of Ser46 phosphorylation on the molecular properties of HPr and its capability to act as the co-repressor of carbon catabolite protein A (CcpA). The calculated pK (a) values for a representative set of HPr(Ser46P) structures indicate that the phosphate group of HPr(Ser46P) exists predominantly in the unprotonated form under neutral conditions. A hydrogen bond detected in HPr(Ser46P) between one phosphate-group oxygen and a side-chain hydrogen of Asn43-an amino acid conserved in all HPr proteins of Gram-positive bacteria that regulate their carbon consumption by CCR-might fulfil an important role in CcpA-HPr(Ser46P) complex formation. MD simulations show that the Ser46P-Asn43 hydrogen bond present in the unbound structure is replaced by intermolecular interactions upon complex formation. The degree to which amino acids in the CcpA-HPr(Ser46P) interface contribute to cofactor binding was analyzed by in silico alanine scanning. Lys307, Arg303, Asp296, Val300, and Tyr295 of CcpA were identified as important amino acids for the CcpA-HPr(Ser46P) interaction. Three of these residues are directly involved in sensing the correct phosphorylation state at His15(HPr) and Ser46(HPr). A substitution of interface residues Val319, Val314, Ser316, Leu321 and Gln320 by alanine showed that these amino acids, which contact helix alpha2 of HPr(Ser46P), play a less prominent role for complex formation.

Amino Acid Substitution↗