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

The conformation of cyclo(-D-Pro-Ala4-) as a model for cyclic pentapeptides of the DL4 type.

The conformation of the cyclic pentapeptide cyclo(-D-Pro-Ala(4)-) in solution and in the solid state was reinvestigated using modern NMR techniques. To allow unequivocal characterization of hydrogen bonds, relaxation behavior, and intramolecular distances, differently labeled isotopomers were synthesized. The NMR results, supported by extensive MD simulations, demonstrate unambiguously that the preferred conformation previously described by us, but recently questioned, is indeed correct. The validation of the conformational preferences of this cyclic peptide is important given that this system is a template for several bioactive compounds and for controlled "spatial screening" for the search of bioactive conformations.

Alanine↗
PubMed2006

Microsolvation effects on the excited-state dynamics of protonated tryptophan.

To better understand the complex photophysics of the amino acid tryptophan, which is widely used as a probe of protein structure and dynamics, we have measured electronic spectra of protonated, gas-phase tryptophan solvated with a controlled number of water molecules and cooled to approximately 10 K. We observe that, even at this temperature, the bare molecule exhibits a broad electronic spectrum, implying ultrafast, nonradiative decay of the excited state. Surprisingly, the addition of two water molecules sufficiently lengthens the excited-state lifetime that we obtain a fully vibrationally resolved electronic spectrum. Quantum chemical calculations at the RI-CC2/aug-cc-pVDZ level, together with TDDFT/pw based first-principles MD simulations of the excited-state dynamics, clearly demonstrate how interactions with water destabilize the photodissociative states and increase the excited-state lifetime.

Models, Molecular↗
PubMed2007

Proton shuttles and phosphatase activity in soluble epoxide hydrolase.

Recently, a novel metal Mg2+-dependent phosphatase activity has been discovered in the N-terminal domain of the soluble epoxide hydrolase (sEH), opening a new branch of fatty acid metabolism and providing an additional site for drug targeting. Importantly, the sEH N-terminal fold belongs to the haloacid dehalogenase (HAD) superfamily, which comprises a vast majority of phosphotransferases. Herein, we present the results of a computational study of the sEH phosphatase activity, which includes classical molecular dynamics (MD) simulations and mixed quantum mechanical/molecular mechanics (QM/MM) calculations. On the basis of experimental results, a two-step mechanism has been proposed and herein investigated: (1) phosphoenzyme intermediate formation and (2) phosphoenzyme intermediate hydrolysis. Building on our earlier work, we now provide a detailed description of the reaction mechanism for the whole catalytic cycle along with its free energy profile. The present computations suggest metaphosphate-like transition states for these phosphoryl transfers. They also reveal that the enzyme promotes water deprotonation and facilitates shuttling of protons via a metal-ligand connecting water bridge (WB). These WB-mediated proton shuttles are crucial for the activation of the solvent nucleophile and for the stabilization of the leaving group. Moreover, due to the conservation of structural features in the N-terminal catalytic site of sEH and other members of the HAD superfamily, we suggest a generalization of our findings to these other metal-dependent phosphatases.

Catalysis↗
PubMed2001

Molecular dynamics study of the IIA binding site in human serum albumin: influence of the protonation state of Lys195 and Lys199.

The IIA binding site of human serum albumin (HSA) preferentially binds hydrophobic organic anions of medium size (e.g., aspirin, benzylpenicillin, warfarin, etc.) and bilirubin. This binding ability is particularly important for the distribution, metabolism, and efficacy of drugs. In addition, HSA can also covalently link to different IIA substrates owing to the presence of a highly reactive residue, Lys199, which is strategically located in the IIA site. Herein, we present results of three restrained molecular dynamics (MD) simulations of the IIA binding site on the HSA protein. From these simulations, we have determined the influence that the ionization state of the key residue, Lys199, and the nearby Lys195 has on the structure and dynamics of the IIA binding site. When Lys199 is neutral the computed average distances for the most significant interresidue contacts are in good agreement with those estimated from the X-ray coordinates. The analysis of the solvent structure and dynamics indicates that the basic form of Lys199 is likely connected to the acid form of Lys195 through a network of H-bonding water molecules with a donor --> acceptor character. The presence of these water bridges can be important for stabilizing the configuration of the IIA binding site and/or promoting a potential Lys195 --> Lys199 proton-transfer process. These results suggest that both lysine residues located in the IIA binding site of HSA, Lys195 and Lys199, could play a combined and comparable chemical role. Our simulations also give insight into the binding of bilirubin to HSA.

Binding Sites↗
PubMed2001

Improved scoring of ligand-protein interactions using OWFEG free energy grids.

A new approach to rapidly score protein-ligand interactions is tested on several protein-ligand systems. Results using this approach - the OWFEG free energy grid - are quite promising and are generally in better agreement with experiment (in some cases much better) than those obtained employing scoring techniques currently in wide use. The OWFEG free energy grid is generated from a one-window free energy perturbation MD simulation (Pearlman, D. A. J. Med. Chem. 1999, 42, 4313-4324). The OWFEG approach is applied to three protein systems: IMPDH, MAP kinase p38, and HIV-1 aspartyl protease. OWFEG scores are compared to experimental K(i) and IC50 data in each case. Empirical scoring functions applied to the same systems for comparison include ChemScore, Piecewise Linear Potential (PLP), and Dock energy score.

Databases, Factual↗
PubMed1994

Computer modeling analysis for enediyne chromophore-apoprotein complex of macromolecular antitumor antibiotic C-1027.

On the basis of the van der Waals and electrostatic potential energy calculations and energy minimization (EM), a model for the apoprotein-chromophore complex of C-1027 has been constructed by docking the apoprotein (Apo-C1027) and the chromophore (C1027-Chr). The model of the apoprotein-chromophore complex also suggests the most probable stereochemistry of the chromophore (8R,9S,13R,17R), specific interactions between Apo-C1027 and C1027-Chr, and the origin of the stabilization of the chromophore by the apoprotein. The acetylenic bond of C1027-Chr appears to be stabilized by (1) hydrophobic interactions with the bottom of the pocket of Apo-C1027, (2) orbital interaction with Cys36-Cys45 disulfide bond, (3) van der Waals contact with Pro76, and (4) pi-pi stacking with the benzene moiety of the chromophore. In this model, molecular dynamics (MD) simulation of the chromophore indicates that the acetylene bond distance (C2-C7) of the apoprotein-bound chromophore is longer than that of the unbound chromophore.

Amino Acid Sequence↗
PubMed2002

Molecular dynamics and free energy analyses of cathepsin D-inhibitor interactions: insight into structure-based ligand design.

In this study, we compare the calculated and experimental binding free energies for a combinatorial library of inhibitors of cathepsin D (CatD), an aspartyl protease. Using a molecular dynamics (MD)-based, continuum solvent method (MM-PBSA), we are able to reproduce the experimental binding affinity for a set of seven inhibitors with an average error of ca. 1 kcal/mol and a correlation coefficient of 0.98. By comparing the dynamical conformations of the inhibitors complexed with CatD with the initial conformations generated by CombiBuild (University of California, San Francisco, CA, 1995), we have found that the docking conformation observed in an X-ray structure of one peptide inhibitor bound to CatD (Proc. Natl. Acad. Sci. U.S.A. 1993, 90, 6796-6800) is in good agreement with our MD simulation. However, the DOCK scoring function, based on intermolecular van der Waals and electrostatics, using a distance-dependent dielectric constant (J. Comput. Chem. 1992, 13, 505-524), poorly reproduces the trend of experimental binding affinity for these inhibitors. Finally, the use of PROFEC (J. Comput.-Aided Mol. Des. 1998, 12, 215-227) analysis allowed us to identify two possible substitutions to improve the binding of one of the better inhibitors to CatD. This study offers hope that current methods of estimating the free energy of binding are accurate enough to be used in a multistep virtual screening protocol.

Algorithms↗
PubMed2006

In-situ synthesis of a tacrine-triazole-based inhibitor of acetylcholinesterase: configurational selection imposed by steric interactions.

Recently, researchers have used acetylcholinesterase (AChE) as a reaction vessel to synthesize its own inhibitors. Thus, 1 (syn-TZ2PA6), a femtomolar AChE inhibitor, which is formed in a 1:1 mixture with its anti-isomer by solution phase reaction from 3 (TZ2) and 4 (PA6), can be synthesized exclusively inside the AChE gorge. Our computational approach based on quantum mechanical/molecular mechanical (QM/MM) calculations, molecular dynamics (MD), and targeted molecular dynamics (TMD) studies answers why 1 is the sole product in the AChE environment. Ab initio QM/MM results show that the reaction in the AChE gorge occurs when 3/azide and 4/acetylene are extended in a parallel orientation. An MD simulation started from the final structure of QM/MM calculations keeps the azide's and acetylene's parallel orientations intact for 10 ns of simulation time. A TMD simulation applied on an antiparallel azide-acetylene conformation flips the acetylene easily to bring it to a position that is parallel to azide. A second set of QM/MM calculations performed on this flipped structure generates a similar minimum-energy path as obtained previously. Even a TMD simulation carried out on a parallel azide-acetylene conformation could not deform their parallel arrangement. All of these results, thus, imply that inside the AChE gorge, the azide group of 3 and the acetylene group of 4 always remain parallel, with the consequence that 1 is the only product. The architecture of the gorge plays an important role in this selective formation of 1.

Acetylcholinesterase↗
PubMed2006

Proposal of a new binding orientation for non-peptide AT1 antagonists: homology modeling, docking and three-dimensional quantitative structure-activity relationship analysis.

A three-dimensional model of the AT1 receptor was constructed by means of a homology modeling procedure, using the X-ray structure of bovine rhodopsin as the initial template and taking into account the available site-directed mutagenesis data. The docking of losartan and its active metabolite EXP3174, followed by 1 ns of molecular dynamics (MD) simulation inserted into the phospholipid bilayer, suggested a different binding orientation for these antagonists from those previously proposed. Furthermore, the docking of several non-peptide antagonists was used as an alignment tool for the development of a three-dimensional quantitative structure-activity relationship (3D-QSAR) model, and the good results confirmed our binding hypothesis and the reliability of the model.

Amino Acid Sequence↗
PubMed2006

Insights on resistance to reverse transcriptase: the different patterns of interaction of the nucleoside reverse transcriptase inhibitors in the deoxyribonucleotide triphosphate binding site relative to the normal substrate.

It is presently known that the long-term failure in the treatment of AIDS with the currently available nucleotide reverse transcriptase inhibitors (NRTIs) is related to the development of resistance by reverse transcriptase (RT) at the binding or incorporation level or both, or subsequent to the nucleotide incorporation (excision). To achieve greater insight on the differential interactions of two NRTIs that are mainly discriminated by different mechanisms, 2',3'-didehydro-2',3'-dideoxythymidine-5'-triphosphate (d4TTP, that is, phosphorylated stavudine) and 2',3'-dideoxycytidine-5'-triphosphate (ddCTP, that is, phosphorylated zalcitabine), with the primer/template (p/t) and with the N binding site of reverse transcriptase (RT) in relation to the normal substrate (dNTP), we have conducted a series of molecular dynamics (MD) simulations. We propose that the different resistance profiles arise from the different conformations adopted by the inhibitors at the N site. d4TTP adopts an ideal conformation for catalysis because it forms an ion-dipole intramolecular interaction with the beta-phosphate oxygen of the triphosphate, as does the normal substrate. In ddCTP, the lack of this essential interaction results in a different, noncatalytic conformation.

Binding Sites↗
PubMed2006

Synthesis of 2-methyl-3-indolylacetic derivatives as anti-inflammatory agents that inhibit preferentially cyclooxygenase 1 without gastric damage.

Novel substituted 2-methyl-3-indolylacetic derivatives were synthesized and evaluated for their activity in vitro and in vivo on COX-1 and COX-2. Active compounds were screened to determine their gastrointestinal tolerability in vivo in the rat. Results showed that 3 and 4 preferentially inhibited COX-1 in vitro and in vivo. MD simulations indicated an induced fit for COX-1 but not for COX-2, probably because of a lower plasticity of the latter.

Acute Disease↗
PubMed2000

Synthesis and conformational studies of peptidomimetics containing furanoid sugar amino acids and a sugar diacid.

Furanoid sugar amino acids (1) were synthesized and used as dipeptide isosteres to induce interesting turn structures in small linear peptides. They belong to a new variety of designed hybrid structures that carry both amino and carboxyl groups on rigid furanose sugar rings. Four such molecules, 6-amino-2,5-anhydro-6-deoxy-D-gluconic acid (3, Gaa) and its mannonic (4, Maa), idonic (5, Iaa), and a 3,4-dideoxyidonic (6, ddIaa) congeners were synthesized. The synthesis followed a novel reaction path in which an intramolecular 5-exo S(N)2 opening of the hexose-derived terminal aziridine ring in 2 by the gamma-benzyloxy oxygen with concomitant debenzylation occurred during pyridinium dichromate oxidation of the primary delta-hydroxyl group to carboxyl function, leading to the formation of furanoid sugar amino acid frameworks in a single step. Incorporation of these furanoid sugar amino acids into Leu-enkephalin replacing its Gly-Gly portion gave analogues 8-11. Detailed structural analysis of these molecules by circular dichroism (CD) and various NMR techniques in combination with constrained molecular dynamics (MD) simulations revealed that two of these analogues, 8a and 10a, have folded conformations composed of an unusual nine-membered pseudo beta-turn-like structure with a strong intramolecular H-bond between LeuNH --> sugarC3-OH. This, in turn, brings the two aromatic rings of Tyr and Phe in close proximity, a prerequisite for biological activities of opioid peptides. The analgesic activities of 8a,b determined by mouse hot-plate and tail-clip methods were similar to that of Leu-enkephalin methyl ester. The syn disposition of the beta-hydroxycarboxyl motif on the sugar rings appears to be the driving force to nucleate the observed turn structures in some of these molecules (8 and 10). Repetition of the motif on both sides of a furanose ring resulted in a novel molecular design of sugar diacid, 2,5-anhydro-D-idaric acid (7, Idac). Bidirectional elongation of the diacid moieties of 7 with identical peptide strands led to the formation of a C2-symmetric reverse-turn mimetic 12 which displayed a very ordered structure consisting of identical intramolecular H-bonds at two ends between LeuNH --> sugar-OH, the same as in 8 and 10.

Amino Acids↗
PubMed2003

Reverse turn induced pi-facial selectivity during polyaniline-supported cobalt(II) salen catalyzed aerobic epoxidation of N-cinnamoyl L-proline derived peptides.

A novel chemo- and diastereoselective aerobic epoxidation of the N-cinnamoyl peptides catalyzed by polyaniline-supported cobalt(II) salen (PASCOS) is described. The N-cinnamoyl proline derived peptides 1 show a high pi-facial selectivity during these epoxidations. The origin of this diastereoselectivity in 1 has been attributed to (i) the propensity of the N-cinnamoyl proline amide to exist predominantly as trans rotamer in CDCl3, DMSO-d6, and CH3CN medium and (ii) existence of these peptides as organized structures (gamma- and beta-turns) due to the presence of intramolecular hydrogen bonds. An extensive solution NMR and MD simulation study on 1d and 1f indicates that the origin of the high pi-facial selectivity is due to the well-defined gamma- and beta-turns which result in the hindrance of one face of the cinnamoyl double bond in the transition state of the epoxidation reaction.

Amino Acids↗
PubMed2006

Synthesis and structure of novel conformationally constrained 1',2'-azetidine-fused bicyclic pyrimidine nucleosides: their incorporation into oligo-DNAs and thermal stability of the heteroduplexes.

[structures: see text] The synthesis of novel 1',2'-aminomethylene bridged (6-aza-2-oxabicyclo[3.2.0]heptane) "azetidine" pyrimidine nucleosides and their transformations to the corresponding phosphoramidite building blocks (20, 39, and 42) for automated solid-phase oligonucleotide synthesis is reported. The novel bicyclonucleoside "azetidine" monomers were synthesized by two different strategies starting from the known sugar intermediate 6-O-benzyl-1,2:3,4-bis-O-isopropylidene-D-psicofuranose. Conformational analysis performed by molecular modeling (ab initio and MD simulations) and NMR showed that the azetidine-fused furanose sugar is locked in a North-East conformation with pseudorotational phase angle (P) in the range of 44.5-53.8 degrees and sugar puckering amplitude (phi(m)) of 29.3-32.6 degrees for the azetidine-modified T, U, C, and 5-Me-C nucleosides. Thermal denaturation studies of azetidine-modified oligo-DNA/RNA heteroduplexes show that the azetidine-fused nucleosides display improved binding affinities when compared to that of previously synthesized North-East sugar constrained oxetane fused analogues.

Azetidines↗
PubMed2006

Conformational analysis of 2-substituted cyclobutane-alpha-amino acid derivatives. A synergistic experimental and computational study.

An extensive conformational study of different 2-substituted cyclobutane-alpha-amino acid derivatives in the solid state, in the gas phase, and in solution has been carried out. The study combines experimental techniques, such as X-ray diffraction and NMR spectroscopy, and computational methods, such as DFT calculations and molecular dynamics (MD) simulations, in a set solvent. The study reveals that the substituent at C2 in the cyclobutane ring, when fixed in an equatorial position, modulates the conformational preference of the ring-puckering.

Journal Article↗
PubMed2000

Synthesis and conformational properties of oligonucleotides incorporating 2'-O-phosphorylated ribonucleotides as structural motifs of pre-tRNA splicing intermediates.

To synthesize oligonucleotides containing 2'-O-phosphate groups, four kinds of ribonucleoside 3'-phosphoramidite building blocks 6a-d having the bis(2-cyano-1,1-dimethylethoxy)thiophosphoryl (BCMETP) group were prepared according to our previous phosphorylation procedure. These phosphoramidite units 6a-d were not contaminated with 3'-regioisomers and were successfully applied to solid-phase synthesis to give oligodeoxyuridylates 15, 16 and oligouridylates 21, 22. Self-complementary Drew-Dickerson DNA 12mers 24-28 replaced by a 2'-O-phosphorylated ribonucleotide at various positions were similarly synthesized. In these syntheses, it turned out that KI(3) was the most effective reagent for oxidative desulfurization of the initially generated thiophosphate group to the phosphate group on polymer supports. Without using this conversion step, a tridecadeoxyuridylate 17 incorporating a 2'-O-thiophosphorylated uridine derivative was also synthesized. To investigate the effect of the 2'-phosphate group on the thermal stability and 3D-structure of DNA(RNA) duplexes, T(m) measurement of the self-complementary oligonucleotides obtained and MD simulation of heptamer duplexes 33-36 were carried out. According to these analyses, it was suggested that the nucleoside ribose moiety phosphorylated at the 2'-hydroxyl function predominantly preferred C2'-endo to C3'-endo conformation in DNA duplexes so that it did not significantly affect the stability of the DNA duplex. On the other hand, the 2'-modified ribose moiety was expelled to give a C3'-endo conformation in RNA duplexes so that the RNA duplexes were extremely destabilized.

Base Sequence↗
PubMed2005

Predictive assessment of surface resistances in zeolite membranes using atomically detailed models.

The diffusive transport of molecules through nanoporous membranes is determined by both intracrystalline diffusion and mass transport resistances associated with entering and leaving the membrane material. We compare two methods for assessing the relative importance of these resistances based on atomically detailed descriptions of the membrane material. For extremely thin membranes, net transport can be assessed using dual control volume grand canonical molecular dynamics (DCV-GCMD). We show that previous implementations of this technique may have been influenced by nonisothermal effects in interfacial regions and suggest a simple remedy to this situation. We also introduce an approximate method that uses information only from equilibrium MD simulations, which avoids the significant computational expense associated with DCV-GCMD. This approximate method can be used to rapidly assess the importance of interface-related resistances to mass transport over broad ranges of membrane operating conditions. This method will be useful in allowing a rapid determination of whether these interface resistances are significant in practical experimental situations. These two methods are compared by considering the transport of CH(4) and CF(4) through defect-free silicalite membranes.

Journal Article↗
PubMed2005

Insertion of dipolar molecules in channels of a centrosymmetric organic zeolite: molecular modeling and experimental investigations on diffusion and polarity formation.

The mechanism of insertion of p-nitroaniline (PNA) and its diffusion behavior in channels of the hexagonal host structure of tris(o-phenylenedioxy)cyclotriphosphazene (TPP) was investigated by means of molecular modeling tools. Strong preferential sites in the bulk were found to be due to pi-pi and NH-pi interactions between PNA and channel walls of TPP. MD simulations showed that diffusion can be characterized by jumps from one site to the next, occurring mainly because of the dynamic flexibility of the host structure. Calculations of host-guest interactions between the TPP surface and PNA approaching the entrance of channels with its terminal H2N-first or O2N-first revealed that the H2N-first insertion is clearly preferred. Preferential insertion of PNA is found to be the reason for polar effects, observed experimentally. Because of a distinct guest-host recognition at the surface, guest-guest interactions were found to have a minor influence on polarity.

Journal Article↗