SEARCH · PubMed Health

Results for “Molecular Dynamic simulation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2Linked to original sources

PubMed1990

Studies on somatostatin with time-resolved spectroscopy and molecular dynamics simulations.

Somatostatin was studied by time-resolved fluorescence spectroscopy and by molecular dynamics simulations. The results obtained indicate the existence of mainly one conformation in DMSO, and the existence of several conformations in other solvents. Molecular dynamics simulations showed that the most important region for activity (residue 7-9) is also the most flexible region in the peptide.

Computer Simulation↗

PubMed1991

Molecular dynamics simulation of Lewis blood groups and related oligosaccharides.

Molecular dynamics simulations without explicit inclusion of solvent molecules have been performed to study the motions of Lewisa and Lewisb blood group oligosaccharides, and two blood group A tetrasaccharides having type I and type II core chains. The blood group H trisaccharide has also been studied and compared with the blood group A type II core chain. The potential energy surface developed by Rasmussen and co-workers was used with the molecular mechanics code CHARMM. The lowest energy minima of the component disaccharide fragments were obtained from conformational energy mapping. The lowest energy minima of these disaccharide fragments were used to build the tri- and tetrasaccharides that were further minimized before the actual heating/equilibration and dynamics simulations. The trajectories of the disaccharide fragments, e.g., Fuc alpha- (1----4)GlcNAc, Gal beta-(1----4)GlcNAc, etc., show transitions among various minima. However, the oligosaccharides were found to be dynamically stable and no transitions to other minimum energy conformations were observed in the time series of the glycosidic dihedral angles even during trajectories as long as 300 ps. The stable conformations of the glycosidic linkages in the oligosaccharides are not necessarily the same as the minimum energy conformation of the corresponding isolated disaccharides. The average fluctuations of the glycosidic angles in the oligosaccharides were well within the range of +/- 15 degrees. The results of these trajectory calculations were consistent with the relatively rigid single-conformation models derived for these oligosaccharides from 1H-nmr data.

Carbohydrate Conformation↗

PubMed1990

Molecular dynamics simulations of active site mutants of triosephosphate isomerase.

Molecular dynamics simulations of triosephosphate isomerase (TIM) and of some active site TIM mutants were performed in an attempt to elucidate possible interactions important for catalytic activity and binding. A variety of active site residues in TIM have been altered, resulting in all cases in decreases in catalytic activity. Second-site suppressor mutants were characterized for two of these active site mutants. The pseudorevertants have increased activity compared to the single mutant from which they were derived and, surprisingly, in both cases the increase in activity is a result of the replacement of an active site serine for proline. We performed simulations of wild-type TIM and the active site mutants with the substrate dihydroxyacetone phosphate bound both noncovalently and covalently. The noncovalent complexes were used to examine interactions important to binding while the covalent complexes are models of the transition state structure for enolization, which is the rate-determining step for the mutants. The difference between these two states, then, is related to the catalytic activity. We found various protein-substrate interactions that improved in the noncovalent mutant complexes, which correlates with the experimentally observed increase in binding affinity upon mutation. In the covalent complexes we observed improved electrostatic stabilization of the transition state upon introduction of Pro, which is also consistent with the experimental data. Our simulations reproduce the highly co-operative nature of the interactions in the active site and suggest that this approach may be useful for identifying particularly promising sites for mutation.

Binding Sites↗

PubMed1987

Lipid bilayer polypeptide interactions studied by molecular dynamics simulation.

A model membrane with a polypeptide alpha-helix inserted has been simulated by molecular dynamics at a temperature well above the gel/liquid crystalline phase transition temperature. Order parameters of the lipids and other equilibrium and dynamic quantities have been calculated. Three systems, polyglycine constrained into an alpha-helical configuration, glycophorin with similarly conformationally constrained backbone and finally glycophorin free to change its backbone conformation, have been studied. In all cases there was an ordering of the chains close to the helix. This effect was, however, much smaller for glycophorin with its rather bulky side chains than for polyglycine. The dynamics of the lipids were affected by the neighbouring helix, not drastically however. Lateral diffusion and reorientational time correlations of lipids close to the helix were slower than for the bulk ones, but not more than two or three times. Thus, we did not find any evidence of bound or frozen boundary lipids.

Glycophorins↗

PubMed1992

Molecular dynamics simulation in vacuo and in solution of [Aib5,6-D-Ala8] cyclolinopeptide A: a conformational and comparative study.

The conformation of a Cyclolinopeptide A analogue, c-(Pro-Pro-Phe-Phe-Aib-Aib-Ile-D-Ala-Val), has been investigated by means of molecular dynamics simulations, in various molecular environments. The molecular dynamics results are compared with that obtained for Cyclolinopeptide A and a detailed analysis of the different behaviour for the two compounds is reported. A complete analysis of hydrogen bonds is presented.

Amino Acid Sequence↗

PubMed1991

Molecular dynamics simulations of dinucleoside and dinucleoside-drug crystal hydrates.

Molecular dynamics simulations have been performed on the dinucleoside monophosphates rGpC and dCpG, the latter in its intercalation complex with the acridine drug proflavine. The simulations were performed on the crystal structures, with crystallographically-located solvent molecules. It was found that satisfactory results were best obtained with restraints placed on the movements of the water molecules. Motions of individual atoms have been examined in terms of rms fluctuations and anisotropy and correlation functions. Relative motions of groups (phosphates, sugars, bases and proflavine molecules) have been analysed.

Computer Simulation↗

PubMed1991

Molecular dynamics simulation of the aqueous solvation of sugars.

Recently, several molecular dynamics simulations of the aqueous solvation of carbohydrates have been reported. These studies represent the first theoretical picture available of the microscopic character of sugar solutions, and may provide explanations of the unusual and complex behavior of this class of molecules in solution. This paper will discuss two MD simulations of D-glucopyranose, including a free energy perturbation calculation of the anomeric free energy difference. Solvation was found to have little effect upon the mean conformational structure of the pyranoid rings, but the presence of solvent significantly affected the motions and orientations of the exocyclic groups. Adjacent functional groups of the sugar rings were found to mutually perturb one another's hydration, depending upon the local stereochemistry, which may prove to play a part in the observed anomeric preferences of the sugars. From a component analysis of the free energy of solvation of the two anomers of D-glucopyranose, it was found that a large solvation term favors the beta anomer, which is the form found to be preferred in aqueous solution.

Carbohydrates↗

PubMedDate not supplied

Molecular dynamics simulations and the conformational mobility of blood group oligosaccharides.

Molecular dynamics simulations were carried out without explicit consideration of solvent to explore the conformational mobility of blood group A and H oligosaccharides. The potential energy force field of Rasmussen and co-workers was used with the CHARMM program on a number of disaccharide and trisaccharide models composed of fucose, galactose, glucose, N-acetyl glucosamine, and N-acetyl galactosamine chosen to represent various fragments of blood group oligosaccharides. In agreement with results of earlier studies, stable chair conformations were found for each pyranoside from which no transitions were detected in simulations as long as 800 ps. Exocyclic dihedral angles, including that of C5-C6, generally show numerous transitions on a time scale of approximately 5-30 ps. The dihedral angles of some but not all glycosidic linkages of blood group oligosaccharides show transitions on the time scale of 30-50 ps, implying that the extent of internal motion in blood group oligosaccharides depends strongly on linkage stereochemistry. For certain blood group A and H oligosaccharides that show limited internal motion in these simulations, we argue that the calculations are consistent with our previous analysis of 1H nuclear Overhauser enhancement (NOE) data that imply single conformations over a wide range of temperature and solvent conditions. While the trajectories are consistent with 13C T1 values that have been interpreted as indicating rigid conformations, measurements of 13C-NOE and T1 as a function of magnetic field strength are proposed as an improved method for experimental detection of the internal motion that is suggested for certain oligosaccharides in these simulations. The results of these simulations differ substantially from those of peptides of a similar molecular weight in that the oligosaccharides show much less internal motion.

Blood Group Antigens↗

PubMed2026

Elucidating the Mechanism of Xiaoqinglong Decoction in Chronic Urticaria Treatment: An Integrated Approach of Network Pharmacology, Bioinformatics Analysis, Molecular Docking, and Molecular Dynamics Simulations.

INTRODUCTION: Xiaoqinglong Decoction (XQLD) is a traditional Chinese medicinal formula commonly used to treat chronic urticaria (CU). However, its underlying therapeutic mechanisms remain incompletely characterized. This study employed an integrated approach combining network pharmacology, bioinformatics, molecular docking, and molecular dynamics simulations to identify the active components, potential targets, and related signaling pathways involved in XQLD's therapeutic action against CU, thereby providing a mechanistic foundation for its clinical application. METHODS: The active components of XQLD and their corresponding targets were identified using the Traditional Chinese Medicine Systems Pharmacology (TCMSP) database. CU-related targets were retrieved from the OMIM and GeneCards databases. Subsequently, core components and targets were determined via protein-protein interaction (PPI) network analysis and component-target-pathway network construction. Topological analyses were performed using Cytoscape software to prioritize core nodes within these networks. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted via the DAVID database to identify enriched biological processes and signaling pathways. Molecular docking was performed to evaluate binding interactions between key components and core targets, while molecular dynamics (MD) simulations were employed to assess the stability of the component-target complexes with the lowest binding energy. Finally, CU-related targets of XQLD were validated using datasets from the Gene Expression Omnibus (GEO) database. RESULTS: A total of 135 active components and 249 potential targets of XQLD were identified, alongside 1,711 CU-related targets. Core components, such as quercetin, kaempferol, beta-sitosterol, naringenin, stigmasterol, and luteolin, exhibited high degree values in the constructed networks. The core targets identified included AKT1, TNF, IL6, TP53, PTGS2, CASP3, BCL2, ESR1, PPARG, and MAPK3. GO and KEGG pathway enrichment analyses revealed the PI3K-Akt signaling pathway as a central regulatory mechanism. Molecular docking studies demonstrated strong binding affinities between active components and core targets, with the stigmasterol-AKT1 complex exhibiting the lowest binding energy (-11.4 kcal/mol) and high stability in MD simulations. Validation using GEO datasets identified 12 core genes shared between CU-related targets and XQLD-associated targets, including PTGS2 and IL6, which were also prioritized as core targets in the network pharmacology analyses. DISCUSSION: This study comprehensively integrates multidisciplinary approaches to clarify the potential molecular mechanisms of XQLD in treating CU, highlighting its multitarget and multipathway synergistic effects. Molecular docking and dynamics simulations confirm the stable interaction between stigmasterol and the core target AKT1. Additionally, GEO dataset analysis verifies the pathogenic relevance of targets such as PTGS2 and IL6, significantly enhancing the credibility of our findings. These results provide a modern scientific basis for the traditional therapeutic effects of XQLD on CU and have important implications for developing multitarget treatments for this condition. However, this study mainly relies on database mining and computational simulations. Further in vitro and in vivo experimental validations are needed to confirm the predicted component-target-pathway interactions. CONCLUSION: This study identifies the active components, potential targets, and pathways through which XQLD exerts therapeutic effects on CU. These findings provide a theoretical foundation for further mechanistic studies and support their clinical application in the treatment of CU.

Molecular Docking Simulation↗

PubMedDate not supplied

Molecular dynamics simulations of opioid peptide analogs containing multiple conformational restrictions.

Molecular dynamics simulations were performed on the potent and slightly mu-receptor selective cyclic dermorphin analog H-Tyr-D-Orn-Phe-Glu-NH2 as well as on analogs containing a conformationally restricted phenylalanine derivative in place of Phe in the 3 position of the peptide sequence. Peptides studied included the potent and highly mu-selective analogs H-Tyr-D-Orn-Aic-Glu-NH2 (Aic = 2-aminoindan-2-carboxylic acid), H-Tyr-D-Orn-Atc-Glu-NH2 (Atc = 2-aminotetralin-2-carboxylic acid) and H-Tyr-D-Orn-D-Atc-Glu-NH2, and the weakly active analog H-Tyr-D-Orn-Tic-Glu-NH2 (Tic = tetrahydroisoquinoline-3-carboxylic acid). Four different starting conformations were chosen for each peptide, and after equilibration each simulation was allowed to proceed for 100 picoseconds at 600 degrees K. The 14-membered ring structures in the Phe-, Aic-, L- and D-Atc-containing analogs showed moderate structural flexibility, while the peptide ring in the Tic-containing analog was more rigid. As theoretically predicted, the phi 3 and psi 3 angles of the Aic-, L- and D-Atc-containing analogs were limited to values of either about +50 degrees or -50 degrees during almost the entire period of the simulations. In the Tic-containing analog the phi 3 and psi 3 angles were 0 degrees and 90 degrees, respectively, and did not change for the entire duration of the simulation. The side chains of the constrained amino acids showed limited movement, but transitions between the allowed conformations did occur on the time scale of the simulations.(ABSTRACT TRUNCATED AT 250 WORDS)

Computer Simulation↗

PubMed1985

Molecular dynamics simulation of photodissociation of carbon monoxide from hemoglobin.

A molecular dynamics simulation of the photodissociation of carbon monoxide from the alpha subunit of hemoglobin is described. To initiate photodissociation, trajectories of the liganded molecule were interrupted, the iron-carbon monoxide bond was broken, and the parameters of the iron-nitrogen bonds were simultaneously altered to produce a deoxyheme conformation. Heme potential functions were used that reproduce the energies and forces for the iron out-of-plane motion obtained from quantum mechanical calculations. The effect of the protein on the rate and extent of the displacement of the iron from the porphyrin plane was assessed by comparing the results with those obtained for an isolated complex of heme with imidazole and carbon monoxide. The half-time for the displacement of the iron from the porphyrin plane was found to be 50-150 fs for both the protein and the isolated complex. These results support the interpretation of optical absorption studies using 250-fs laser pulses that the iron is displaced from the porphyrin plane within 350 fs in both hemoglobin and a free heme complex in solution.

Animals↗

PubMed1992

Molecular dynamics simulation of bovine prothrombin fragment 1 in the presence of calcium ions.

Early solvation-induced structural reorganization of calcium prothrombin fragment 1 is simulated with molecular dynamics. Initial coordinates are those of the 2.2-A resolution crystal structure [Soriano-Garcia, M., Padmanabhan, K., de Vos, A. M., & Tulinsky, A. (1992) Biochemistry 31, 2554-2556]. The molecular dynamics code AMBER, appropriately modified to include long-range (less than or equal to 22.0 A) ionic forces, was employed. The solution structure appears to equilibrate within 100 ps. Although minor changes are seen in various structural domains, the early solution structure basically maintains an intricate network of nine gamma-carboxyglutamic acid (Gla) residues encapsulating seven calcium ions. However, the Gla domain moves with respect to the kringle domain. This motion is mainly due to the movement of Ser34-Leu35 that appears to be a flexible hinge between the domains. The N-terminus of Ala 1 is in a tightly bound complex with three Gla residues that remains stable in the solution structure when the long-range electrostatic cutoff is employed and the near planar alignment of the seven calcium ions is only slightly distorted. The simulation structure is discussed in terms of experiments that studied calcium ion-induced quenching of the intrinsic fluorescence, protection of the N-terminal amino group from acetylation by calcium ions, chemical modification of the N-terminus to a trinitrophenyl derivative, and the possibility of a calcium-binding site(s) in the kringle domain.

Amino Acid Sequence↗

PubMed1991

Molecular dynamics simulations of the cytochrome c3-rubredoxin complex from Desulfovibrio vulgaris.

Molecular dynamics simulations have been carried out on the complex formed between the tetraheme cytochrome c3 and the iron protein rubredoxin from the sulfate-reducing bacterium Desulfovibrio vulgaris. These simulations were performed both with explicit solvent water molecules included, and without solvent molecules using a distance-dependent dielectric constant to approximate the screening effects of solvent. The results of both simulations are strikingly different, indicating that the representation of environmental effects is important in such simulations. For example, a striking adaptation of the two proteins seen in the nonsolvated simulation is not seen when explicit solvent water is included; in fact, the complex appears to become weaker in the solvated simulation. Nonetheless, the iron-iron distance decreases more significantly in the solvated simulation than in the nonsolvated simulation. It was found that in both cases molecular dynamics optimized the structures further than energy minimization alone.

Computer Simulation↗

PubMed1988

Importance of folded monomer and extended antiparallel dimer structures as enkephalin active conformation. Molecular dynamics simulations of [Met5]enkephalin in water.

Simulations of the molecular dynamics of the [Met5]enkephalin monomer and dimer structures in water have been carried out. The dynamic trajectories have been analyzed in terms of the distances between intra- or intermolecular polar atoms. The time-correlated conformational transitions of an extended monomer structure have been converged into a stationary state among the beta-bend folded forms. However, the dynamics simulation of an extended antiparallel dimer structure has shown no noticeable conformation change. These results imply that both the beta-bend monomer and the extended dimer structures exist together as the fundamental conformation of enkephalins.

Calorimetry↗

PubMed2026

Exploring the mechanism of Shengmai San in treating lung adenocarcinoma based on bioinformatics and molecular dynamics simulation.

To investigate the mechanism of Shengmai San (SMS) in the treatment of lung adenocarcinoma (LUAD) based on an integrated strategy combining "network pharmacology, bioinformatics, molecular docking, and molecular dynamics simulation," aiming to provide a precise combination therapy strategy and identify potential bioactive compounds. Differentially expressed genes in LUAD were identified from the Gene Expression Omnibus database using R (originally developed at Bell Laboratories and currently managed by Lucent Technologies). SMS components (ginseng, Ophiopogon japonicus, and Schisandra chinensis) were retrieved from encyclopaedia of traditional Chinese medicine, with Lipinski-compliant compounds selected. Compound targets were predicted via SwissTargetPrediction and Similarity Ensemble Approach. Intersecting targets between differentially expressed genes and compound targets were identified for "herbs-compounds-targets-disease" network construction. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed. Hub targets were identified by analyzing the protein-protein interaction network. High-prognostic relevance targets were screened from The Cancer Genome Atlas. Compounds targeting these were identified through the herbs-compounds-targets-disease network, and absorption, distribution, metabolism, excretion, and toxicity-compliant compounds were selected using SwissADME (a web-based tool provided by the Molecular Modeling Group of the Swiss Institute of Bioinformatics). Core regulatory targets were identified through molecular docking, with complex stability assessed by molecular dynamics simulations. The key bioactive compounds of SMS for treating LUAD were identified as 7-hydroxy-2,5-dimethyl-4H-1-benzopyran-4-one, N-trans-feruloyltyramine, paprazine, and (E)-N-[(2S)-2-hydroxy-2-(4-hydroxyphenyl)ethyl]-3-(4-hydroxyphenyl)prop-2-enamide. Hub targets included AURKA, CCNA2, CCNB1, CDK1, CHEK1, KIF11, NEK2, PLK1, TTK, and TYMS. Among these, CDK1, CHEK1, and PLK1 demonstrated both high-prognostic relevance and strong binding affinity with SMS, emerging as core regulatory targets for SMS in LUAD treatment. Mechanistically, SMS exerts its anticancer effects primarily by modulating the tumor necrosis factor, interleukin-17, cell cycle, and Lipid and atherosclerosis signaling pathways. The active components of SMS, such as paprazine, may exert antitumor effects partly through downregulating CDK1, CHEK1, and PLK1 expression. Although the present study did not examine drug-resistance models or combination regimens, our findings raise the possibility that, in patients with high expression of these genes, combining SMS with standard chemotherapy or targeted therapy could potentially enhance chemosensitivity and mitigate the development of resistance. This hypothesis, however, requires formal testing in appropriate preclinical models and functional validation studies.

Molecular Dynamics Simulation↗

PubMed1988

Molecular dynamics simulations of ribonuclease T1. Effect of solvent on the interaction with 2'GMP.

Molecular dynamics simulations in vacuum and with a water sphere around the active site were performed on the 2'GMP-RNase T1 complex. The presence of water led to the maintenance of the 2'-GMP-RNase T1 interactions as compared to the X-ray structure, including the hydrogen bonds implicated in the enzyme-inhibitor recognition process. The sidechain of His92 in the molecular dynamics water simulation, however, hydrogen bonds directly to the phosphate of 2'GMP in contrast to the X-ray structure but in support of the role of that residue in the enzyme's catalytic mechanism. Fluctuations of active-site residues are not strongly influenced by water, possibly owing to the exclusion of water by the bound 2'GMP, which did show an increase in mobility. Analysis of the 2'GMP-RNase T1 interactions versus time reveal an equilibrium fluctuation in the presence of water, leading to a less favorable 2'GMP-RNase T1 interaction energy, suggesting a possible relationship between picosecond fluctuations and inhibitor dissociation occurring in the millisecond time domain.

Endoribonucleases↗

PubMed1992

A molecular dynamics simulation of bacteriophage T4 lysozyme.

An analysis of a 400 ps molecular dynamics simulation of the 164 amino acid enzyme T4 lysozyme is presented. The simulation was carried out with all hydrogen atoms modeled explicitly, the inclusion of all 152 crystallographic waters and at a temperature of 300 K. Temporal analysis of the trajectory versus energy, hydrogen bond stability, r.m.s. deviation from the starting crystal structure and radius of gyration, demonstrates that the simulation was both stable and representative of the average experimental structure. Average structural properties were calculated from the enzyme trajectory and compared with the crystal structure. The mean value of the C alpha displacements of the average simulated structure from the X-ray structure was 1.1 +/- 0.1 A; differences of the backbone phi and psi angles between the average simulated structure and the crystal structure were also examined. Thermal-B factors were calculated from the simulation for heavy and backbone atoms and both were in good agreement with experimental values. Relationships between protein secondary structure elements and internal motions were studied by examining the positional fluctuations of individual helix, sheet and turn structures. The structural integrity in the secondary structure units was preserved throughout the simulation; however, the A helix did show some unusually high atomic fluctuations. The largest backbone atom r.m.s. fluctuations were found in non-secondary structure regions; similar results were observed for r.m.s. fluctuations of non-secondary structure phi and psi angles. In general, the calculated values of r.m.s. fluctuations were quite small for the secondary structure elements. In contrast, surface loops and turns exhibited much larger values, being able to sample larger regions of conformational space. The C alpha difference distance matrix and super-positioning analyses comparing the X-ray structure with the average dynamics structure suggest that a 'hinge-bending' motion occurs between the N- and C-terminal domains.

Bacteriophage T4↗

PubMed1990

Molecular dynamics simulation of papain-E-64 (N-[N-(L-3-trans-carboxyoxirane-2-carbonyl)-L-leucyl]agmatine) complex.

To investigate the possible binding mode of E-64 (N-[N-(L-3-trans-carboxyoxirane-2-carbonyl)-L-leucyl]agmatine), a potent cysteine protease inhibitor, to papain active site, molecular dynamics simulations were applied to two complex forms: R- and S- configurational forms of E-64 C2 atom for the covalent bond formation with the papain Cys-25 SH group. The tertiary structures of the papain-E-64 complexes were built by visual interactive modelling and the energy minimization technique, and were subjected to the dynamics simulations of 10 ps. Although no significant difference was observed between the potential energies of energy-minimized R- and S-complex forms, the molecular dynamics simulations suggested that the hydrogen bonding mode of the former form is more advantageous than that of the latter one. Comparing with the hydrogen bonds observed in the papain-E-64 complex crystal, it could be concluded that the present molecular dynamics simulation reflects well the three-dimensional structure concerning the interaction of E-64 with the papain active site. The conformational characteristics of E-64 and its possible interaction mode with papain were also discussed.

Computer Simulation↗