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 973 records · Page 54Linked to original sources
PubMed2005

Protonation of the acidic residues in the transmembrane cation-binding sites of the ca(2+) pump.

The ionization states of the acidic residues around the Ca2+-binding sites of sarcoplasmic reticulum Ca2+ ATPase are studied by continuum electrostatic calculations and all-atom molecular dynamics simulations with explicit solvent and phospholipids. The two methods consistently indicate that Glu58 and Glu908 are protonated at neutral pH. The Ca2+ coordination and the H-bonds formed by the protonation of Glu58 and Glu908 are stable in an MD simulation, whereas the H-bonds are disrupted and the Ca2+ coordination geometry is severely altered in another simulation treating these residues unprotonated. The results clearly indicate that the H-bonds formed by protonation of Glu58 and Glu908 provide extra stability for the Ca2+-binding sites of Ca2+ ATPase.

Binding Sites↗
PubMedDate not supplied

Molecular dynamics of a vasopressin V2 receptor in a phospholipid bilayer membrane.

Molecular dynamics simulations were carried out for a V2 receptor (V2R) model embedded in a dimyristoylphosphatidylcholine (DMPC) bilayer. Both free and ligand-bound states of V2R were modeled. Our initial V2R model was obtained using a rule-based automated method for GPCR modeling and refined using constrained simulated annealing in vacuo. The docking site of the native vasopressin ligand was selected and justified upon consideration of ligand-receptor interactions and structure-activity data. The primary purpose of this work was to investigate the usefulness of MD simulation of an integral membrane protein like a GPCR receptor, upon inclusion of a carefully parameterized surrounding lipid membrane and water. Physical properties of the system were evaluated and compared with the fully hydrated pure DMPC bilayer membrane. The solvation interactions, individual lipid-protein interaction and fluctuations of the protein, the lipid, and water were analyzed in detail. As expected, the membrane-spanning helices of the protein fluctuate less than the peripheral loops do. The protein appears to disturb the local lipid structure. Simulations were carried out using AMBER 4.1 package upon constant number-pressure-temperature (NPT) conditions on massively parallel computers Cray T3E and IBM SP2.

Animals↗
PubMed2004

Vibrational spectra from atomic fluctuations in dynamics simulations. II. Solvent-induced frequency fluctuations at femtosecond time resolution.

The midinfrared (MIR) spectra of molecules in polar solvents exhibit inhomogeneously broadened bands whose spectral positions are shifted as compared to the gas phase. The shifts are caused by interactions with structured solvation shells and the broadenings by fluctuations of these interactions. The MIR spectra can be calculated from hybrid molecular dynamics (MD) simulations, which treat the solute molecule by density functional theory and the solvent by molecular mechanics by the so-called instantaneous normal mode analysis (INMA) or by Fourier transforming the time correlation function (FTTCF) of the molecular dipole moment. In Paper I of this work [M. Schmitz and P. Tavan, J. Chem. Phys. 121, 12233 (2004)] we explored an alternative method based on generalized virial (GV) frequencies noting, however, that GV systematically underestimates frequencies. As shown by us these artifacts are caused by solvent-induced fluctuations of the (i) equilibrium geometry, (ii) force constants, and (iii) normal mode directions as well as by (iv) diagonal and (v) off-diagonal anharmonicities. Here we now show, by analyzing the time scales of fluctuations and sample MD trajectories of formaldehyde in the gas phase and in water, that all these sources of computational artifacts can be made visible by a Fourier analysis of the normal coordinates. Correspondingly, the error sources (i) and (iii)-(v) can be removed by bandpass filtering, as long as the spectral signatures of the respective effects are well separated from the fundamental band. Furthermore, the artifacts arising from effect (ii) can be strongly diminished by a time-resolved version of the GV approach (TF-GV). The TF-GV method then yields for each mode j a trajectory of the vibrational frequency omega(j)(tmid R:tau) at a time resolution tau>tau(j), which is only limited by the corresponding oscillation time tau(j)=2pi/omega(j) and, thus, is in the femtosecond range. A correlation analysis of these trajectories clearly separates the librational motions from the conformational dynamics of the solvation shells and yields the inhomogeneously broadened MIR spectra, if the theory of motional narrowing is properly included. The MIR spectrum of formaldehyde in solution obtained by TF-GV agrees very well with the FTTCF result, if one applies the so-called "harmonic approximation" quantum correction factor and a temperature scaling to the FTTCF intensities. Also for INMA an excellent agreement is achieved if one disregards a slight INMA overestimate of linewidths.

Journal Article↗
PubMed2005

Receptor rigidity and ligand mobility in trypsin-ligand complexes.

The trypsin-like serine proteases comprise a structurally similar family of proteins with a wide diversity of biological functions. Members of this family play roles in digestion, hemostasis, immune responses, and cancer metastasis. Bovine trypsin is an archetypical member of this family that has been extensively characterized both functionally and structurally, and that preferentially hydrolyzes Arg/Lys-Xaa peptide bonds. We have used molecular dynamics (MD) simulations to study bovine trypsin complexed with the two noncovalent small-molecule ligands, benzamidine and tranylcypromine, that have the same hydrogen-bond donating moieties as Arg and Lys side-chains, respectively. Multiple (10) simulations ranging from 1 ns to 2.2 ns, with explicit water molecules and periodic boundary conditions, were performed. The simulations reveal that the trypsin binding pocket residues are relatively rigid regardless of whether there is no ligand, a high-affinity ligand (benzamidine), or a low-affinity ligand (tranylcypromine). The thermal average of the conformations sampled by benzamidine bound to trypsin is planar and consistent with the planar internal geometry of the benzamidine crystallographic model coordinates. However, the most probable bound benzamidine conformations are +/-25 degrees out of plane, implying that the observed X-ray electron density represents an average of densities from two mirror symmetric, nonplanar conformations. Solvated benzamidine has free energy minima at +/-45 degrees , and the induction of a more planar geometry upon binding is associated with approximately 1 kcal/mol of intramolecular strain. Tranylcypromine's hydrogen-bonding pattern in the MD differs substantially from that inferred from the X-ray electron density. Early in simulations of this system, tranylcypromine adopts an alternative binding conformation, changing from the crystallographic conformation, with a direct hydrogen bond between its amino moiety and the backbone oxygen of Gly219, to one having a bridging water molecule. This result is consistently seen with the CHARMM22, Amber, or OPLS-AA force fields. The trypsin-tranylcypromine hydrogen-bonding pattern observed in the simulations also occurs as the crystallographic binding mode of the Lys15 side-chain of bovine pancreatic trypsin inhibitor bound to trypsin. In this latter cocrystal, a bridging crystallographic water does reside between the side-chain's amino group and the trypsin Gly219 backbone oxygen. Furthermore, the trypsin-tranylcypromine simulations sample two different stable noncrystallographic binding poses. These data suggest that some of the electron density ascribed to tranylcypromine in the X-ray model is rather due to a bound water molecule, and that multiple tranylcypromine binding conformations (crystallographic disorder) may be the cause of ambiguous electron density. The combined trypsin-benzamidine and trypsin- tranylcypromine results highlight the ability of simulations to augment protein-ligand complex structural data by deconvoluting the effects of thermal and structural averaging, and by finding energetically optimal ligand and bound water positions for weakly bound ligands.

Animals↗
PubMed1999

Cationic 5,10,15,20-tetrakis(N-methylpyridinium-4-yl)porphyrin fully intercalates at 5'-CG-3' steps of duplex DNA in solution.

The interaction of 5,10,15, 20-tetrakis(N-methylpyridinium-4-yl)porphyrin (T4MPyP(4+)) with the oligonucleotide DNA duplex [d(GCACGTGC)](2) was studied by two-dimensional (1)H NMR spectroscopy, optical absorbance, circular dichroism, and molecular dynamics simulation employing particle mesh Ewald methods. T4MPyP(4+) is one of the largest aromatic molecules for which intercalative binding to DNA has been proposed, although this has been called into question by recent X-ray crystallographic work [Lipscomb et al. (1996) Biochemistry 35, 2818-2823]. T4MPyP(4+) binding to [d(GCACGTGC)](2) produced a single set of (mostly) upfield-shifted DNA resonances in slow exchange with the resonances of the free DNA. Intra- and intermolecular NOEs observed in the complex showed that the porphyrin intercalates at the central 5'-CG-3' step of the DNA duplex without disrupting the flanking base pairs. Absorption and circular dichroism spectra of the complex also support intercalative binding. Molecular dynamics simulations (using explicit solvent and PME methods), carried out for fully and partially intercalated complexes, yielded stable trajectories and plausible structures, but only the symmetrical, fully intercalated model agreed with NOESY data. Stable hydrogen bonding was observed during 600 ps of MD simulation for the base pairs flanking the binding site.

Binding Sites↗
PubMed2006

Influence of point mutations on the flexibility of cytochrome b5: molecular dynamics simulations of holoproteins.

Two membrane-bound isoforms of cytochrome b5 have been identified in mammals, one associated with the outer mitochondrial membrane (OM b5) and the other with the endoplasmic reticulum (microsomal, or Mc b5). The soluble heme binding domains of OM and Mc b5 have highly similar three-dimensional structures but differ significantly in physical properties, with OM b5 exhibiting higher stability due to stronger heme association. In this study, we present results of 8.5-ns length molecular dynamics simulations for rat Mc b5, bovine Mc b5, and rat OM b5, as well as for two rat OM b5 mutants that were anticipated to exhibit properties intermediate between those of rat OM b5 and the two Mc proteins: the A18S/I32L/L47R triple mutant (OM3M) and the A18S/I25L/I32L/L47R/L71S quintuple mutant (OM5M). Analysis of the structure, fluctuations, and interactions showed that the five b5 variants used in this study differed in organization of their molecular surfaces and heme binding cores in a way that could be used to explain certain experimentally observed physical differences. Overall, our simulations provided qualitative microscopic explanations of many of the differences in physical properties between OM and Mc b5 and two mutants in terms of localized changes in structure and flexibility. They also reveal that opening of a surface cleft between hydrophobic cores 1 and 2 in bovine Mc b5, observed in two previously reported simulations (E. M. Storch and V. Daggett, Biochemistry, 1995, Vol. 34, pp. 9682-9693; A. Altuve, Biochemistry, 2001, Vol. 40, pp. 9469-9483), probably resulted from removal of crystal contacts and likely does not occur on the nanosecond time scale. Finally, the MD simulations of OM5M b5 verify that stability and dynamic properties of cytochrome b5 are remarkably resistant to mutations that dramatically alter the stability and structure of the apoprotein.

Animals↗
PubMed2005

Structural dynamics of precursor and product of the RNA enzyme from the hepatitis delta virus as revealed by molecular dynamics simulations.

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme involved in the replication of a human pathogen, the hepatitis delta virus. Recent crystal structures of the precursor and product of self-cleavage, together with detailed kinetic analyses, have led to hypotheses on the catalytic strategies employed by the HDV ribozyme. We report molecular dynamics (MD) simulations (approximately 120 ns total simulation time) to test the plausibility that specific conformational rearrangements are involved in catalysis. Site-specific self-cleavage requires cytidine in position 75 (C75). A precursor simulation with unprotonated C75 reveals a rather weak dynamic binding of C75 in the catalytic pocket with spontaneous, transient formation of a H-bond between U-1(O2') and C75(N3). This H-bond would be required for C75 to act as the general base. Upon protonation in the precursor, C75H+ has a tendency to move towards its product location and establish a firm H-bonding network within the catalytic pocket. However, a C75H+(N3)-G1(O5') H-bond, which would be expected if C75 acted as a general acid catalyst, is not observed on the present simulation timescale. The adjacent loop L3 is relatively dynamic and may serve as a flexible structural element, possibly gated by the closing U20.G25 base-pair, to facilitate a conformational switch induced by a protonated C75H+. L3 also controls the electrostatic environment of the catalytic core, which in turn may modulate C75 base strength and metal ion binding. We find that a distant RNA tertiary interaction involving a protonated cytidine (C41) becomes unstable when left unprotonated, leading to disruptive conformational rearrangements adjacent to the catalytic core. A Na ion temporarily compensates for the loss of the protonated hydrogen bond, which is strikingly consistent with the experimentally observed synergy between low pH and high Na+ concentrations in mediating residual self-cleavage of the HDV ribozyme in the absence of divalents.

Base Sequence↗
PubMed1999

Conformational dynamics of cytochrome c: correlation to hydrogen exchange.

We study the dynamical fluctuations of horse heart cytochrome c by molecular dynamics (MD) simulations in aqueous solution, at four temperatures: 300 K, 360 K, 430 K, and 550 K. Each simulation covers a production time of at least 1.5 nanoseconds (ns). The conformational dynamics of the system is analyzed in terms of collective motions that involve the whole protein, and local motions that involve the formation and breaking of intramolecular hydrogen bonds. The character of the MD trajectories can be described within the framework of rugged energy landscape dynamics. The MD trajectories sample multiple conformational minima, with basins in protein conformational space being sampled for a few hundred picoseconds. The trajectories of the system in configurational space can be described in terms of diffusion of a particle in real space with a waiting time distribution due to partial trapping in shallow minima. As a consequence of the hierarchical nature of the dynamics, the mean square displacement autocorrelation function, <|x(t) - x(0)|2>, exhibits a power law dependence on time, with an exponent of around 0.5 for times shorter than 100 ps, and an exponent of 1.75 for longer times. This power law behavior indicates that the system exhibits suppressed diffusion (sub-diffusion) in sampling of configurational space at time scales shorter than 100 ps, and enhanced (super-diffusion) at longer time scales. The multi-basin feature of the trajectories is present at all temperatures simulated. Structural changes associated with inter-basin displacements correspond to collective motions of the Omega loops and coiled regions and relative motions of the alpha-helices as rigid bodies. Similar motions may be involved in experimentally observed amide hydrogen exchange. However, some groups showing large correlated motions do not expose the amino hydrogens to the solvent. We show that large fluctuations are not necessarily correlated to hydrogen exchange. For example, regions of the proteins forming alpha helices and turns show significant fluctuations, but as rigid bodies, and the hydrogen bonds involved in the formation of these structures do not break in proportion to these fluctuations. Proteins 1999;36:175-191. Published 1999 Wiley-Liss, Inc.

Amides↗
PubMed2005

Polyoxometalates in solution: Molecular dynamics simulations on the alpha-PW12O40(3-) Keggin anion in aqueous media.

The Keggin anion, PW12O40(3-), is one of the most representative polyoxometalates (POMs). In recent years increasing theoretical work focused on this family of compounds has explained or even predicted some of their properties using quantum mechanics methods. In this report we applied for the first time molecular dynamics (MD) to the title compound to analyze its interactions with water. We used three atomic charge definitions (Mulliken, ChelpG, and formal charges) to carry out MD simulations. The results show that the terminal oxygens of the cluster are invariably most effectively solvated by water because of their prominent position within the framework. On the other hand, bridging oxygens, which are confined in more internal positions, concentrate a smaller portion of the whole solvation. We investigated the hydrogen bonds existing between water and the cluster, confirming that the terminal positions form more contacts with H2O than any other site of the cluster. In the end, the lifetime of such contacts is longer with bridging oxygens, presumably due to their higher atomic charge.

Journal Article↗
PubMed2005

Multiple loop conformations of peptides predicted by molecular dynamics simulations are compatible with nuclear magnetic resonance.

The affinity and selectivity of protein-protein interactions can be fine-tuned by varying the size, flexibility, and amino acid composition of involved surface loops. As a model for such surface loops, we study the conformational landscape of an octapeptide, whose flexibility is chemically steered by a covalent ring closure integrating an azobenzene dye into and by a disulfide bridge additionally constraining the peptide backbone. Because the covalently integrated azobenzene dyes can be switched by light between a bent cis state and an elongated trans state, six cyclic peptide models of strongly different flexibilities are obtained. The conformational states of these peptide models are sampled by NMR and by unconstrained molecular dynamics (MD) simulations. Prototypical conformations and the free-energy landscapes in the high-dimensional space spanned by the phi/psi angles at the peptide backbone are obtained by clustering techniques from the MD trajectories. Multiple open-loop conformations are shown to be predicted by MD particularly in the very flexible cases and are shown to comply with the NMR data despite the fact that such open-loop conformations are missing in the refined NMR structures.

Aminobenzoates↗
PubMed2000

Sequence-dependent variation in DNA minor groove width dictates orientational preference of Hoechst 33258 in A-tract recognition: solution NMR structure of the 2:1 complex with d(CTTTTGCAAAAG)(2).

The solution structure of the dodecamer duplex d(CTTTTGCAAAAG)(2)and its 2:1 complex with the bis -benzimidazole Hoechst 33258 has been investigated by NMR and NOE-restrained molecular dynamics (rMD) simulations. Drug molecules are bound in each of the two A-tracts with the bulky N-methylpiperazine ring of each drug located close to the central TG (CA) step, binding essentially to the narrow minor groove of each A-tract. MD simulations over 1 ns, using an explicit solvation model, reveal time-averaged sequence-dependent narrowing of the minor groove from the 3'-end towards the 5'-end of each TTTT sequence. Distinct junctions at the TpG (CpA) steps, characterised by large positive roll, low helical and propeller twists and rapid AT base pair opening rates, add to the widening of the groove at these sites and appear to account for the bound orientation of the two drug molecules with the N-methylpiperazine ring binding in the wider part of the groove close to the junctions. Comparisons between the free DNA structure and the 2:1 complex (heavy atom RMSD 1.55 A) reveal that these sequence-dependent features persist in both structures. NMR studies of the sequence d(GAAAAGCTTTTC)(2), in which the A-tracts have been inverted with the elimination of the TpG junctions, results in loss of orientational specificity of Hoechst 33258 and formation of multiple bound species in solution, consistent with the drug binding in a number of different orientations.

Allosteric Regulation↗
PubMed1999

A kinetic model for the internal motions of proteins: diffusion between multiple harmonic wells.

The dynamics of collective protein motions derived from Molecular Dynamics simulations have been studied for two small model proteins: initiation factor I and the B1 domain of Protein G. First, we compared the structural fluctuations, obtained by local harmonic approximations in different energy minima, with the ones revealed by large scale molecular dynamics (MD) simulations. It was found that a limited set of harmonic wells can be used to approximate the configurational fluctuations of these proteins, although any single harmonic approximation cannot properly describe their dynamics. Subsequently, the kinetics of the main (essential) collective protein motions were characterized. A dual-diffusion behavior was observed in which a fast type of diffusion switches to a much slower type in a typical time of about 1-3 ps. From these results, the large backbone conformational fluctuations of a protein may be considered as "hopping" between multiple harmonic wells on a basically flat free energy surface.

Kinetics↗
PubMed2003

Molecular mechanics and dynamics simulations of various dispersant models on the water surface (001).

Five dispersant-molecule models of succinimide, acrylate, imide, phenylsulfonic and salicyl were used to study their interactions with the water surface (001). The interaction energy, molecular configuration, charge distribution and radial distribution function (RDF) curve for each of the dispersant molecules were analyzed from the molecular mechanics (MM) and molecular dynamics (MD) simulation results. It can be seen that the system energies, mostly electrostatic and hydrogen bond energies, were reduced significantly when the dispersant molecules interacted with the water surface. The hydrophilic group of a dispersant molecule can attach itself to the water surface firmly and reach for a stable energy-minimized configuration, which is helpful to the dispersants' dispersancy. The influence exerted by the hydrophobic group of the dispersant molecule, which was the substituted hydrocarbon chain of n-octadecanyl in this paper, is discussed in comparison with the naked polar headgroup. Steric configuration, charge distribution and substitute hydrocarbon chain of the dispersant molecule influenced the interaction between dispersants and polar water surface.

Acrylates↗
PubMed2001

A nanosecond molecular dynamics study of antiparallel d(G)7 quadruplex structures: effect of the coordinated cations.

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G7) quadruplex structures with different coordinated ions, namely Na+ and K+ ion, water and Na+ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of approximately 1.5 A from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K+ ions are quite immobile. MD studies at 400K indicate that K+ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K+ ion coordinated structure is more favourable than that for the Na+ ion coordinated antiparallel quadruplex structure.

Cations, Monovalent↗
PubMed1995

Molecular dynamics and structure-based drug design for predicting non-natural nonapeptide binding to a class I MHC protein.

Starting from the known three-dimensional structure of the class I major histocompatibility complex-encoded HLA-B*2705 protein, three non-natural nonapeptides were designed to fit optimally the HLA-B*2705-binding groove. The optimization was performed using structure-based drug design methods and the fact that all the possible interactions of the secondary anchor residue (position 3) with its human leukocyte antigen-binding pocket (pocket D) in nature are not entirely utilized. 150 ps molecular-dynamics (MD) simulation in water was employed to study the stability of the bimolecular complexes with three non-natural peptides (P3 = homophenylalanine, beta-naphthylalanine, alpha-naphthylalanine) as well as with the two natural homologues (P3 = Gly, Leu). Various structural and dynamical properties (atomic fluctuations, solvent-accessible surface areas, peptide Calpha-atom positions) of the simulated bimolecular complexes were used to compare the three non-natural with the two natural ligands. Since the various molecular properties have been shown previously to be related to the binding affinity of nonapeptide ligands to the major histocompatibility complex (MHC) HLA-B*2705 protein, the MD data predict a rather higher stability of MHC-ligand complexes with the three non-natural peptides, suggesting that the unnatural peptides studied show an enhanced binding affinity to the HLA-B*2705 protein. These results are in agreement with the experimental values of a semi-quantitative in vitro assembly assay, performed on the five nonapeptides (P3 = Gly, Leu, homophenylalanine, beta-naphthylalanine, alpha-naphthylalanine), which shows their ability to stabilize the native conformation of the HLA-B*2705 heavy chain and also shows that the three non-natural ligands bind with higher affinity (0.5 micro M) to the MHC protein than the two natural homologues (40 micro M). Thus, this study demonstrates that structural information combined with rational design and molecular-dynamics simulations can illustrate and predict MHC binding and potential T-cell epitope properties as well as contribute to the design of new non-peptidic MHC inhibitors that may be useful for the selective immunotherapy of autoimmune diseases to which HLA alleles are directly associated.

Journal Article↗
PubMed2004

Structure and dynamics of a predicted ferredoxin-like selenoprotein in Japanese encephalitis virus.

Homologues of the selenoprotein glutathione peroxidase (GPx) have been previously identified in poxviruses and in RNA viruses including HIV-1 and hepatitis C virus (HCV). Sequence analysis of the NS4 region of Japanese encephalitis virus (JEV) suggests it may encode a structurally related but functionally distinct selenoprotein gene, more closely related to the iron-binding protein ferredoxin than to GPx, with three highly conserved UGA codons that align with essential Cys residues of ferredoxin. Comparison of the probe JEV sequence to an aligned family of ferredoxin sequences gave an overall 30.3% identity and 45.8% similarity, and was statistically significant at 4.9 S.D. (P < 10(-6)) above the average score computed for randomly shuffled sequences. A 3-dimensional model of the hypothetical JEV protein (JEV model) was constructed by homology modeling using SYBYL, based upon a high resolution X-ray structure of ferredoxin (PDB code: 1awd). The JEV model and the model from 1awd were subsequently subjected to molecular dynamics simulations in aqueous medium using AMBER 6. The solution structure of the JEV model indicates that it could fold into a tertiary structure globally similar to ferredoxin 1awd, with RMSD between the averaged structures of 1.8 A for the aligned regions. The modeling and MD simulations data also indicate that this structure for the JEV protein is energetically favorable, and that it could be quite stable at room temperature. This protein might play a role in JEV infection and replication via TNF and other cellular stimuli mediated via redox mechanisms.

Amino Acid Sequence↗
PubMed2006

[Screening of new HIV inhibitors based on the database of traditional Chinese medicine].

AIM: To report the preliminary result of the HIV inhibitor screening based on cheminformatics tools and the traditional Chinese medicine database. METHODS: Database search was carried out with saquinavir molecule as a template, further screening was made with docking. Detailed studies using molecular dynamics simulation of 50 ps and 200 ps were made with respect to a potential leading compound, leucovorin. RESULTS: The leucovorin molecule distinguished from other molecules as a potential drug candidate and is subject to extensive studies. The bonding profile and energy were calculated with MD simulations. CONCLUSION: Our results could be very helpful when we modify leucovorin or design new inhibitors against HIV.

Anti-HIV Agents↗
PubMed2006

A study of collective atomic fluctuations and cooperativity in the U1A-RNA complex based on molecular dynamics simulations.

Cooperative interactions play an important role in recognition and binding in macromolecular systems. In this study, we find that cross-correlated atomic fluctuations can be used to identify cooperative networks in a protein-RNA system. The dynamics of the RRM-containing protein U1A-stem loop 2 RNA complex have been calculated theoretically from a 10 ns molecular dynamics (MD) simulation. The simulation was analyzed by calculating the covariance matrix of all atomic fluctuations. These matrix elements are then presented in the form of a two-dimensional grid, which displays fluctuations on a per residue basis. The results indicate the presence of strong, selective cross-correlated fluctuations throughout the RRM in U1A-RNA. The atomic fluctuations correspond well with previous biophysical studies in which a multiplicity of cooperative networks have been reported and indicate that the various networks identified in separate individual experiments are fluctuationally correlated into a hyper-network encompassing most of the RRM. The calculated results also correspond well with independent results from a statistical covariance analysis of 330 aligned RRM sequences. This method has significant implications as a predictive tool regarding cooperativity in the protein-nucleic acid recognition process.

Amino Acid Sequence↗