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 19 recordsLinked to original sources
PubMed1992

MD simulation of subtilisin BPN' in a crystal environment.

In this paper we present a molecular dynamics (MD) simulation of subtilisin BPN' in a crystalline environment containing four protein molecules and solvent. Conformational and dynamic properties of the molecules are compared with each other and with respect to the X-ray structure to test the validity of the force field. The agreement between simulated and experimental structure using the GROMOS force field is better than that obtained in the literature using other force fields for protein crystals. The overall shape of the molecule is well preserved, as is the conformation of alpha-helices and beta-strands. Structural differences are mainly found in loop regions. Solvent networks found in the X-ray structure were reproduced by the simulation, which was unbiased with respect to the crystalline hydration structure. These networks seem to play an important role in the stability of the protein; evidence of this is found in the structure of the active site. The weak ion binding site in the X-ray structure of subtilisin BPN' is occupied by a monovalent ion. When a calcium ion is placed in the initial structure, three peptide ligands are replaced by 5 water ligands, whereas a potassium ion retains (in part) its original ligands. Existing force fields yield a reliable method to probe local structure and short-time dynamics of proteins, providing an accuracy of about 0.1 nm.

Computer Simulation↗
PubMed2025

Computational identification of potential antifungal targets against Claviceps purpurea via MD simulation and MM/GBSA.

Ergot alkaloids produced by the fungus Claviceps purpurea pose significant risks to agriculture and human health. This study systematically investigates the pathogenicity of C. purpurea, analyzing five strains for their proteomic profiles, which revealed genetic variability in size and GC content. We identified proteins localized in various cellular compartments, contributing to our understanding of essential cellular processes. A focus on potential drug targets led to the identification of Alpha-N-acetylglucosaminidase, a hydrolase with significant mass and functional relevance, despite not matching a UniProt entry. The 3D structure prediction confirmed its integrity, making it a suitable target for further analysis. Molecular docking identified ligands CID:51,535,944 and CID:145,242,255 with strong binding affinities to Alpha-N-acetylglucosaminidase, highlighting interactions with key residues like TRP138, ARG651. Molecular docking interactions were validated and showed consistency through MD simulation analyses with greater RMSD, RMSF and PL contacts. This research enhances our understanding of C. purpurea, offering insights into its genetic diversity and cellular mechanisms while identifying promising therapeutic targets. The findings contribute to strategies for mitigating the economic and health impacts of C. purpurea infections, paving the way for innovative interventions in sustainable agriculture.

Molecular Docking Simulation↗
PubMed1988

High propeller twist and unusual hydrogen bonding patterns from the MD simulation of (dG)6.(dC)6.

A molecular dynamics (MD) study of (dG)6.(dC)6 including counter ions and 292 water molecules was made. The hydrogen bonding pattern and propeller twist angles for the mini-helix are reported as averages for times spanning 21-30, 31-40, 41-50, and 51-60 ps. The propeller twist angles range from 18 degrees to 38 degrees. Bifurcated and interstrand neighboring base (twisted) hydrogen bonding patterns were found.

Hydrogen Bonding↗
PubMed1992

Molecular modeling of an antigenic complex between a viral peptide and a class I major histocompatibility glycoprotein.

Computer simulation of the conformations of short antigenic peptides (5-10 residues) either free or bound to their receptor, the major histocompatibility complex (MHC)-encoded glycoprotein H-2 Ld, was employed to explain experimentally determined differences in the antigenic activities within a set of related peptides. Starting for each sequence from the most probable conformations disclosed by a pattern-recognition technique, several energy-minimized structures were subjected to molecular dynamics simulations (MD) either in vacuo or solvated by water molecules. Notably, antigenic potencies were found to correlate to the peptides propensity to form and maintain an overall alpha-helical conformation through regular i,i + 4 hydrogen bonds. Accordingly, less active or inactive peptides showed a strong tendency to form i,i + 3 hydrogen bonds at their N-terminal end. Experimental data documented that the C-terminal residue is critical for interaction of the peptide with H-2 Ld. This finding could be satisfactorily explained by a 3-D Q.S.A.R. analysis postulating interactions between ligand and receptor by hydrophobic forces. A 3-D model is proposed for the complex between a high-affinity nonapeptide and the H-2 Ld receptor. First, the H-2 Ld molecule was built from X-ray coordinates of two homologous proteins: HLA-A2 and HLA-Aw68, energy-minimized and studied by MD simulations. With HLA-A2 as template, the only realistic simulation was achieved for a solvated model with minor deviations of the MD mean structure from the X-ray conformation. Water simulation of the H-2 Ld protein in complex with the antigenic nonapeptide was then achieved with the template-derived optimal parameters. The bound peptide retains mainly its alpha-helical conformation and binds to hydrophobic residues of H-2 Ld that correspond to highly polymorphic positions of MHC proteins. The orientation of the nonapeptide in the binding cleft is in accordance with the experimentally determined distribution of its MHC receptor-binding residues (agretope residues). Thus, computer simulation was successfully employed to explain functional data and predicts alpha-helical conformation for the bound peptide.

Computer Simulation↗
PubMed1990

Insights into the function of the zinc hydroxide-Thr199-Glu106 hydrogen bonding network in carbonic anhydrases.

The exact functional role of the zinc hydroxide (water)-Thr199-Glu106 hydrogen bond network in the carbonic anhydrases is unknown. However, from the results of molecular dynamics simulations (MD) we are able to better define its function. From computer graphics analysis and MD simulations on the zinc hydroxide form of human carbonic anhydrase II we find that this interaction forces the hydroxide hydrogen atom to be in a "down" position relative to the deep water-binding pocket. From previous work we have found that this pocket is a high-affinity binding site for CO2. We also note that during the timescale of our simulation (126 ps) the hydrogen bonds between the hydroxide hydrogen atom and Thr199 and the one between Thr199 and Glu106 are not fluxional. We propose that the role of the zinc hydroxide (water)-Thr199-Glu106 hydrogen bond network is to lock the hydrogen atom in the down position in order to expose the CO2 molecule bound in the deep water pocket to a lone pair of the hydroxide oxygen atom. This would allow for the rapid reaction of the CO2 molecule around the zinc ion. Furthermore, if the hydroxide hydrogen atom were not locked in the down position the binding of CO2 to the deep water pocket could be interfered with by the unrestrained hydroxide hydrogen atom (e.g. the N-Zn-O-H torsion could undergo rotational transitions that would partially block the deep water pocket). In summary, the roles we ascribe to this hydrogen bonding network are (1) to allow for facile access of CO2 to the deep water pocket and (2) to allow for maximal exposure of a hydroxide oxygen lone pair to the CO2 carbon atom.

Binding Sites↗
PubMed1991

Mass-weighted molecular dynamics simulation of cyclic polypeptides.

A modified molecular dynamics (MD) method in which atomic masses are weighted was developed previously for studying the conformational flexibility of neuroregulating tetrapeptide Phe-Met-Arg-Phe-amide (FMRF-amide). The method has now been applied to longer and constrained molecules, namely a disulfide-linked cyclic hexapeptide, c[CYFQNC], and its linear and "pseudo-cyclic" analogues. The sampling of dehedral conformational space of teh linear hexapeptide in mass-weighted MD simulations was found to be improved significantly over conventional MD simulations, as in the case of the shorter FMRF-amide molecule studied previously. In the cyclic hexapeptide, the internal constraint of the molecule due to the intramolecular disulfide bond (hence the absence of free terminals in the molecule) does not adversely affect the significant improvement of conformational sampling in mass-weighted MD simulations over normal MD simulations. The pseudo-cyclic polypeptide is identical to the linear CYFQNC molecule in amino acid sequence (i.e., side chains of the cysteine residues are reduced), but the positions of its two terminal heavy atoms were held fixed in space such that the molecule has a nearly cyclic conformation. For this molecule, the mass-weighted MD simulation generated a wide range of polypeptide backbone conformations covering the internal dihedral degrees of freedom; moreover, the physical space of the pseudo-cyclic structure was also sampled in a complete revolution of the entire molecular fragment about the two fixed termini during the simulation. These characteristics suggest that mass-weighted MD can also be an extremely useful method for conformational analyses of constrained molecules and, in particular, for modeling loops on protein surfaces.

Amino Acid Sequence↗
PubMedDate not supplied

On the dependence of molecular conformation on the type of solvent environment: a molecular dynamics study of cyclosporin A.

The dependence of the conformation of cyclosporin A (CPA), a cyclic undecapeptide with potent immunosuppressive activity, on the type of solvent environment is examined using the computer simulation method of molecular dynamics (MD). Conformational and dynamic properties of CPA in aqueous solution are obtained from MD simulations of a CPA molecule dissolved in a box with water molecules. Corresponding properties of CPA in apolar solution are obtained from MD simulations of CPA in a box with carbontetrachloride. The results of these simulations in H2O and in CCl4 are compared to each other and to those of previous simulations of crystalline CPA and of an isolated CPA molecule. The conformation of the backbone of the cyclic polypeptide is basically independent of the type of solvent. In aqueous solution the beta-pleated sheet is slightly weaker and the gamma-turn is a bit less pronounced than in apolar solution. Side chains may adopt different conformations in different solvents. In apolar solution the hydrophobic side chain of the MeBmt residue is in an extended conformation with its hydroxyl group hydrogen bonded to the backbone carbonyl group. In aqueous solution this hydrophobic side chain folds over the core of the molecule and the mentioned hydrogen bond is broken in favor of hydrogen bonding to water molecules. The conformation obtained from the MD simulation in CCl4 nicely agrees with experimental atom-atom distance data as obtained from nmr experiments in chloroform. In aqueous solution the relaxation of atomic motion tends to be slower than in apolar solution.

Chemical Phenomena↗
PubMed1992

Conformational states of a TT mismatch from molecular dynamics simulation of duplex d (CGCGATTCGCG).

The TT mismatch region in duplex d (CGCGATTCGCG) was studied using a 500-ps molecular dynamics (MD) simulation in water, and a series of 1-ps MD simulations and energy minimizations in vacuum. The DNA maintained its duplex structure, although the mismatch region showed significantly higher flexibility than the GC regions. The predominant conformation in the 500-ps MD simulation involved an average -42 degrees propeller twist between T6 and T'6, and a -22 degree buckle between A5 and T'7. One hydrogen bond was formed between T6 and T'6, and another between T6 and the O2 of T'7, with both Watson-Crick hydrogen bonds between A5 and T'7 remaining intact. The minimizations resulted in conformations with the equivalent hydrogen-bonding pattern, as well as ones with "wobble pair" hydrogen bonds between T6 and T'6. However, the wobble pair conformation was found to be unstable in the water simulation.

Base Sequence↗
PubMed1989

Superoxide dismutase: fluctuations in the structure and solvation of the active site channel studied by molecular dynamics simulation.

The molecular dynamics (MD) simulation of superoxide dismutase (SOD) in water is carried out for a total of 23 ps. The simulation system is a 26 A sphere centered at the active site of SOD, including 1602 atoms from SOD and 1761 water molecules. There is no gross deviation from the x-ray structure for the average MD structure. The structure and potential fluctuations around the active site are examined. The results provide new insight to the interactions between SOD and its substrate superoxide.

Binding Sites↗
PubMed1991

Mass-weighted molecular dynamics simulation and conformational analysis of polypeptide.

Atomic motions in protein molecules have been studied by molecular dynamics (MD) simulations; dynamics simulation methods have also been employed in conformational studies of polypeptide molecules. It was found that when atomic masses are weighted, the molecular dynamics method can significantly increase the sampling of dihedral conformation space in such studies, compared to a conventional MD simulation of the same total simulation time length. Herein the theoretical study of molecular conformation sampling by the molecular dynamics-based simulation method in which atomic masses are weighted is reported in detail; moreover, a numerical scheme for analyzing the extensive conformational sampling in the simulation of a tetrapeptide amide molecule is presented. From numerical analyses of the mass-weighted molecular dynamics trajectories of backbone dihedral angles, low-resolution structures covering the entire backbone dihedral conformation space of the molecule were determined, and the distribution of rotationally stable conformations in this space were analyzed quantitatively. The theoretical analyses based on the computer simulation and numerical analytical methods suggest that distinctive regimes in the conformational space of the peptide molecule can be identified.

Amino Acid Sequence↗
PubMed1987

Molecular dynamics simulations of cyclosporin A: the crystal structure and dynamic modelling of a structure in apolar solution based on NMR data.

The conformation of the immunosuppressive drug cyclosporin A (CPA), both in apolar solution and in crystalline state, has been studied by computer simulation techniques. Three molecular dynamics (MD) simulations have been performed: one modelling the crystal structure and two modelling the structure in apolar solution, using a restrained MD approach in which data from nuclear magnetic resonance (NMR) and infrared (IR) spectroscopy are taken into account. The simulation of the crystalline state (MDC) concerns a system of 4 unit cells containing 16 cyclosporin A molecules and 22 water molecules, which is simulated using crystalline periodic boundary conditions. The simulations modelling the apolar solvent conformation (MDS) concern one isolated cyclosporin A molecule. In these simulations an extra term in the interatomic potential function is used, which forces the molecule to satisfy a set of 57 atom-atom distance constraints originating from nuclear Overhauser effects (NOEs) obtained from NMR spectroscopy and one distance constraint deduced from IR spectroscopy. From a comparison of the results of the crystal simulation to those of the X-ray experiment in terms of structure, atomic fluctuations, hydrogen bond pattern, etc., it is concluded that the force field that is used yields an adequate representation of crystalline cyclosporin A. Secondly, it is shown that the dynamic modelling technique that is used to obtain a structure in a polar solution from NMR distance information works well. Starting from initial conformations which have a root mean square difference of 0.14 nm both distance restrained MD simulations converge to the same final solution structure. A comparison of the crystal structure of cyclosporin A and the one in apolar solution shows that there are significant differences. The overall difference in atomic positions is 0.09 nm for the C alpha atoms and 0.17 nm for all atoms. In apolar solution, the molecule is slightly more bent and the side chains of 1 MeBmt and 10 MeLeu adopt a different conformation.

Computer Simulation↗
PubMed2026

Genomic and structural insights into the atpB L173I substitution: modulation of the F₀ rotor architecture in Mycobacterium tuberculosis ATP synthase and altered Bedaquiline binding dynamics.

The F₀F₁ ATP synthase of Mycobacterium tuberculosis (M. tuberculosis) is an essential membrane-embedded rotary motor responsible for ATP synthesis and maintenance of the proton motive force in bacteria. The transmembrane F₀ domain comprises the c-subunit (atpE) and the a-subunit (atpB). Their coordinated interactions are needed for proton translocation and torque generation. Bedaquiline (BDQ), FDA-approved diarylquinoline for the treatment of multidrug-resistant tuberculosis (MDR-TB), targets the F₀ motor by binding at the a-c interface and inhibiting rotary catalysis. To the best of our knowledge, this study represents the first attempt to analyze the effects of mutations in the atpB protein on its structural stability in the F₀ domain, thereby highlighting the novelty of this work. In this study, we integrated Indian whole-genome sequencing (WGS) datasets (PRJNA37907) with long-timescale (1000 ns) membrane-embedded molecular dynamics (MD) simulations. Among 57 atpB mutations identified from WGS analysis, L173I was selected for structural and MD analysis. L173I is located at the atpB-atpE interface near the BDQ-binding region, despite V177L and S184A showing higher prevalence. Comparative MD simulations encompassed four systems: wild-type apo, wild-type with BDQ, L173I apo, and L173I with BDQ. Structural interrogation revealed that the L173I substitution induces subtle destabilization of the global fold of the atpB-atpE complex relative to the apo state, while more critically attenuating inter-subunit contacts between the a-subunit and the c-ring. These perturbations provide a mechanistic rationale for reduced BDQ susceptibility, arising from altered interfacial dynamics rather than complete abrogation of drug binding. This integrative genomic-structural framework advances our understanding of ATP synthase-mediated resistance in M. tuberculosis.

Diarylquinolines↗
PubMed1987

Conformational analysis of a cyclic thymopoietin-analogue by 1H n.m.r. spectroscopy and restrained molecular dynamics simulations.

The internuclear distances of the cyclic thymopoietin derivative c[D-Val-Tyr-Arg-Lys-Glu] have been determined using two-dimensional nuclear Overhauser n.m.r. spectroscopy. These distances are used as constraints for a restrained Molecular Dynamics (MD) simulation. The two starting structures used for the calculations consist of a beta and gamma turn for model 1 and two gamma turns for model 2. The rms difference in atomic positions of the two conformations is 0.242 nm. They converge during the restrained MD simulation to the same final structure. The positional rms difference of the time averaged (5-14 ps) conformations is 0.011 nm. The hydrogen bond pattern is similar to that of model 1, but in addition we find three more gamma turns. The vicinal NH-C alpha H couplings agree well with those calculated from the time averaged structures.

Magnetic Resonance Spectroscopy↗
PubMed1988

On the occurrence of three-center hydrogen bonds in cyclodextrins in crystalline form and in aqueous solution: comparison of neutron diffraction and molecular dynamics results.

Three-center (bifurcated) hydrogen bonds may play a role by serving as an intermediate state between different dynamically changing hydrogen bonding patterns. Hydrogen bonding configurations can be studied experimentally by neutron diffraction and theoretically by computer simulation techniques. Here, both methods are used to analyse the occurrence of three-center hydrogen bonds in crystals of cyclodextrins. Almost all experimentally observed three-center hydrogen bonds in the crystal are reproduced in the molecular dynamics (MD) simulations, even as far as the detailed asymmetric geometry is concerned. On the basis of this result a MD simulation of cyclodextrin in aqueous solution is searched for the occurrence of three-center hydrogen bonds. Significant differences are found. In solution more different three-center hydrogen bonds per alpha-cyclodextrin molecule are observed than in the crystal but the population (existence as percent of the simulation period) of each three-center hydrogen bond is lower in solution than in crystal. These may indeed serve as intermediate states in the process of changing one hydrogen bonding pattern into another.

Chemical Phenomena↗
PubMed1992

Molecular dynamics study of the conformational behavior of a representative elastin building block: Boc-Gly-Val-Gly-Gly-Leu-OMe.

The conformational behavior of the synthetic peptide, Boc-Gly-Val-Gly-Gly-Leu-OMe, containing the X-Gly-Gly and Gly-Gly-X (X = Val or Leu) repeating sequences and constituting a fragment of elastin was investigated by molecular mechanics and molecular dynamics (MD) simulation. The results suggest that, irrespective of the approximations used, the molecule shows a manifold of low energy conformations characterized by gamma-turns and type II beta-turns. Furthermore, MD simulations point out a conformational floppiness due to very low barriers between different conformations. Experimental CD measurements in a virtually apolar medium (dioxane--epsilon = 2.209), which better mimics the vacuum conditions of the simulation, support the theoretical results. The general emerging picture, indicating the molecule as characterized by a combination of flexibility with conformational preferences, is in agreement with previous experimental findings and enriches of new aspects the description of the microscopic behavior of this molecule suggesting more detailed interpretation of previous data.

Amino Acid Sequence↗
PubMed1992

Solvent dielectric effects on protein dynamics.

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

Chymotrypsin↗
PubMed2026

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

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

Pseudomonas aeruginosa↗
PubMed2025

Integrated network pharmacology, molecular docking and experimental validation to investigate the mechanism of tannic acid in nasopharyngeal cancer.

Tannic acid (TA) is the primary bioactive component in the gallnut (Galla chinensis) and has exhibited the anticancer effects. However, the mechanism of its anti-cancer activity in nasopharyngeal carcinoma (NPC) remains unclear. This research aims to explore the underlying mechanism of TA in the treatment of nasopharyngeal cancer using network pharmacology, molecular docking and experimental validation. Firstly, the targets of TA and NPC were predicted and collected through databases, and the intersection targets were identified. Subsequently, protein-protein interaction (PPI) network analysis, Gene Ontology (GO) enrichment, Kyoto Encyclopedia of Genes Genomes (KEGG) pathway enrichment analysis, molecular docking and molecular dynamics (MD) simulation were conducted to uncover the potential mechanisms of TA in treatment of NPC. Finally, in vitro experiments were utilized to verify the mechanism of TA with anticancer activity in NPC. The results of network pharmacology revealed 42 intersection targets between NPC-related targets and TA-related targets. The phosphoinositide 3-kinase (PI3K)/protein kinase B (AKT) signaling was identified as the main target pathway of TA against NPC. Additionally, molecular docking and MD simulation confirmed the closely binding affinities of TA with AKT1. Furthermore, the results of in vitro experiments demonstrated that TA exerts anticancer activity against NPC by targeting the PI3K/AKT signaling pathway, leading to the suppression of cell proliferation. TA is a promising therapeutic candidate for NPC through PI3K/AKT signaling pathway. These results provide insights into the clinical application of TA, particularly when considered in combination with other therapeutic modalities.

Molecular Docking Simulation↗