SEARCH · PubMed Health

Results for “Molecular docking and molecular dynamics 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 127 records · Page 7Linked to original sources

PubMed2005

Construction of a 3D model of nattokinase, a novel fibrinolytic enzyme from Bacillus natto. A novel nucleophilic catalytic mechanism for nattokinase.

A three-dimensional structural model of nattokinase (NK) from Bacillus natto was constructed by homology modeling. High-resolution X-ray structures of Subtilisin BPN' (SB), Subtilisin Carlsberg (SC), Subtilisin E (SE) and Subtilisin Savinase (SS), four proteins with sequential, structural and functional homology were used as templates. Initial models of NK were built by MODELLER and analyzed by the PROCHECK programs. The best quality model was chosen for further refinement by constrained molecular dynamics simulations. The overall quality of the refined model was evaluated. The refined model NKC1 was analyzed by different protein analysis programs including PROCHECK for the evaluation of Ramachandran plot quality, PROSA for testing interaction energies and WHATIF for the calculation of packing quality. This structure was found to be satisfactory and also stable at room temperature as demonstrated by a 300ps long unconstrained molecular dynamics (MD) simulation. Further docking analysis promoted the coming of a new nucleophilic catalytic mechanism for NK, which is induced by attacking of hydroxyl rich in catalytic environment and locating of S221.

Amino Acid Sequence↗

PubMed2003

Improved prediction of HIV-1 protease-inhibitor binding energies by molecular dynamics simulations.

BACKGROUND: The accurate prediction of enzyme-substrate interaction energies is one of the major challenges in computational biology. This study describes the improvement of protein-ligand binding energy prediction by incorporating protein flexibility through the use of molecular dynamics (MD) simulations. RESULTS: Docking experiments were undertaken using the program AutoDock for twenty-five HIV-1 protease-inhibitor complexes determined by x-ray crystallography. Protein-rigid docking without any dynamics produced a low correlation of 0.38 between the experimental and calculated binding energies. Correlations improved significantly for all time scales of MD simulations of the receptor-ligand complex. The highest correlation coefficient of 0.87 between the experimental and calculated energies was obtained after 0.1 picoseconds of dynamics simulation. CONCLUSION: Our results indicate that relaxation of protein complexes by MD simulation is useful and necessary to obtain binding energies that are representative of the experimentally determined values.

Computer Simulation↗

PubMed2026

Sanguinarine as a multi-target therapeutic candidate for laryngeal cancer: insights from network pharmacology, molecular dynamics and in vitro validation.

OBJECTIVE: To identify the core targets and elucidate the potential molecular mechanisms of sanguinarine (SA) against laryngeal squamous cell carcinoma (LSCC), and to validate its antitumor effects in vitro. METHODS: Potential targets of SA were predicted using SwissTargetPrediction, TargetNet, and SuperPred and intersected with LSCC-related targets obtained from the GeneCards, OMIM, and DISEASES databases. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed. A protein-protein interaction (PPI) network was constructed using the STRING database (combined score > 0.900), and topological parameters including degree centrality (DC), betweenness centrality (BC), closeness centrality (CC), eigenvector centrality (EC), and local average connectivity (LAC) were calculated in Cytoscape to identify core genes based on median thresholds. Molecular docking and 100-ns molecular dynamics (MD) simulations were conducted for epidermal growth factor receptor (EGFR), Phosphatidylinositide-3-kinase catalytic subunit alpha (PIK3CA), phosphatidylinositol-4,5-biphosphate 3-kinase catalytic subunit β (PIK3CB), phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit delta (PIK3CD), and Non-Receptor Tyrosine Kinase (SRC). The effects of SA on LSCC were evaluated using CCK-8, colony formation, Transwell migration, and wound-healing assays in TU177 cells and TU212. RESULTS: A total of 213 common targets were identified, which were significantly enriched in PI3K-Akt signaling, EGFR tyrosine kinase inhibitor resistance, and adhesion- and migration-related pathways. The PPI network comprised 259 nodes and 259 edges, from which five core genes-PIK3CA, PIK3CB, PIK3CD, EGFR, and SRC-were identified. Molecular docking revealed strong binding affinities between SA and the PI3K family proteins (- 9.79 to - 10.96 kcal/mol), as well as EGFR (- 8.58 kcal/mol) and SRC (- 6.77 kcal/mol). MD simulations indicated greater stability of SA complexes with EGFR and PI3K family members compared with SRC. In vitro assays demonstrated that SA significantly inhibited TU177 cell and TU212 cell proliferation, colony formation, and migration. CONCLUSION: SA may exert anti-laryngeal cancer effects through synergistic multi-target inhibition centered on the EGFR/SRC/PI3K signaling axis, highlighting its potential as a promising therapeutic candidate for LSCC.

Humans↗

PubMed2005

Calpha-H...O = C hydrogen bonds contribute to the specificity of RGD cell-adhesion interactions.

BACKGROUND: The Arg-Gly-Asp (RGD) cell adhesion sequence occurs in several extracellular matrix molecules known to interact with integrin cell-surface receptors. Recently published crystal structures of the extracellular regions of two integrins in complex with peptides containing or mimicking the RGD sequence have identified the Arg and Asp residues as key specificity determinants for integrin recognition, through hydrogen bonding and metal coordination interactions. The central Gly residue also appears to be in close contact with the integrin surface in these structures. RESULTS: When hydrogen atoms are modelled on the central Gly residue with standard stereochemistry, the interaction between this residue and a carbonyl group in the integrin surface shows all the hallmarks of Calpha-H...O = C hydrogen bonding, as seen in the collagen triple helix and in many crystal structures of small organic molecules. Moreover, molecular dynamic simulations of the docking of RGD-containing fragments on integrin surfaces support the occurrence of these interactions. There appears to be an array of four weak and conventional hydrogen bonds lining up the RGD residues with main chain carbonyl groups in the integrin surface. CONCLUSIONS: The occurrence of weak Calpha-H...O = C hydrogen bonds in the RGD-integrin interaction highlights the importance of the conserved Gly residue in the RGD motif and its contribution to integrin-ligand binding specificity. Our analysis shows how weak hydrogen bonds may also play important biological roles by contributing to the specificity of macromolecular recognition.

Arginine↗

PubMedDate not supplied

Development and validation of empirical force field parameters for netropsin.

The netropsin molecule preferentially binds to the four consecutive A.T base pairs of the DNA minor groove and could therefore inhibit the expression of specific genes. The understanding of its binding on a molecular level is indispensable for computer-aided design of new antitumor agents. This knowledge could be obtained via molecular dynamics (MD) and docking simulations, but in this case appropriate force field parameters for the netropsin molecule should be explicitly defined. Our parametrization was based on the results of quantum chemical calculations. The resulting set of parameters was able to reproduce bond lengths, bond angles, torsional angles of the ab initio minimized geometry within 0.03 A, 3 deg and 5 deg, respectively, and its vibrational frequencies with a relative error of 4.3% for low and 2.8% for high energy modes. To show the accuracy of the developed parameters we calculated an IR spectrum of the netropsin molecule using MD simulation and found it to be in good agreement with the experimental one. Finally, we performed a 10 ns long MD simulation of the netropsin-DNA complex immersed in explicit water. The overall complex conformation remained stable at all times, and its secondary structure was well retained.

Algorithms↗

PubMed2005

Antagonist and agonist binding models of the human gonadotropin-releasing hormone receptor.

G-protein-coupled receptors (GPCRs) constitute one of the most important classes of drug targets. Since the first high-resolution structure of a GPCR was determined by Palczewski and co-workers [K. Palczewski, T. Kumasaka, T. Hori, C.A. Behnke, H. Motoshima, B.A. Fox, I. Le Trong, D.C. Teller, T. Okada, R.E. Stenkamp, M. Yamamoto, M. Miyano, Crystal structure of rhodopsin: a G-protein-coupled receptor, Science 289 (2000) 739-745], development of in silico models of rhodopsin-like GPCRs could be rationally founded. In this work, we present a model of the human gonadotropin-releasing hormone receptor based on the rhodopsin structure. The transmembrane helices are modeled by homology, while the extra- and intra-cellular loops are modeled in such a way that experimentally determined interactions and microdomains (e.g., hydrophobic cores) are retained. We conclude that specifically tailored models, compared to more automatic approaches, have the benefit that known interactions are easily introduced early in the homology modeling. Furthermore, tailored models, although more tedious to construct, are better suited for drug lead finding and for compound optimization. To test the stability of the receptor, we performed a 1 ns molecular dynamics simulation. Moreover, we docked two agonists (native GnRH and Triptorelin, [dTrp(6)]-GnRH) and two antagonists (Cetrorelix, dNal(1)-dCpa(2)-dPal(3)-Ser(4)-Tyr(5)-dCit(6)-Leu(7)-Arg(8)-Pro(9)-dAla(10)), and the covalently constrained dicyclic decapeptide dicyclo(1,1'-5/4-10)[Ac-Glu(1)(Gly(1)')-dCpa(2)-dTrp(3)-Asp(4)-dbu(5)-dNal(6)-Leu(7)-Arg(8)-Pro(9)-dpr(10)-NH(2)] into the putative receptor binding site. The docked ligand conformations result in ligand-receptor interactions that are generally in good agreement with site-directed mutagenesis and ligand-binding studies presented in the literature. Our results indicate that the binding conformation of the antagonists differs from that of the agonists. This difference can be linked to the activation or inhibition of the receptor.

Binding Sites↗

PubMed2026

Discovery of NAT-6-321056 as a novel modulator of VEGFR2 signaling to suppress tumor angiogenesis.

Vascular endothelial growth factor receptor 2 (VEGFR2) is a master regulator of angiogenesis and cancer progression. However, current VEGFR2 modulators face significant challenges, including off-target toxicity and acquired resistance, underscoring the urgent need for novel therapeutic agents with improved efficacy and safety profiles. Here, we reported that virtual screening of 39,442 natural products from the ZINC natural products-derived library, coupled with molecular docking and molecular dynamics (MD) simulations to evaluate the binding stability of candidate compounds, identified NAT-6-321056 as a highly promising modulator of VEGFR2 signaling. Biological evaluations demonstrated that NAT-6-321056 exerted potent inhibition on the growth of a broad spectrum of cancer cells, including both solid tumors and hematological malignancies. In EA.hy 926 endothelial cells and SK-N-DZ neuroblast cells, the compound significantly suppressed proliferation, migration, and invasion. Microscale thermophoresis (MST) confirmed direct binding of NAT-6-321056 to VEGFR2 with favorable affinity. Kinase profiling against a panel of 33 kinases indicated that NAT-6-321056 exhibited a multi-kinase modulation profile. Mechanistic studies revealed that NAT-6-321056 suppressed the expression of hypoxia-inducible factor 1-alpha (HIF-1α) and was associated with reduced VEGFR2 phosphorylation and attenuation of the downstream ERK/JNK/AKT signaling pathways. Moreover, NAT-6-321056 exhibited robust in vivo anti-angiogenic effects in both the chick chorioallantoic membrane (CAM) assay and transgenic zebrafish vascular fluorescence imaging models. Computational absorption, distribution, metabolism, excretion, and toxicity (ADMET) prediction suggested acceptable drug-like properties. Collectively, these findings demonstrated that NAT-6-321056 is a promising modulator of VEGFR2 signaling with potent anti-angiogenic activity and represents a viable candidate for cancer therapy.

Vascular Endothelial Growth Factor Receptor-2↗

PubMed2005

Exploring the active-site of a rationally redesigned lipase for catalysis of Michael-type additions.

Michael-type additions of various thiols and alpha,beta-unsaturated carbonyl compounds were performed in organic solvent catalyzed by wild-type and a rationally redesigned mutant of Candida antarctica lipase B. The mutant lacks the nucleophilic serine 105 in the active-site; this results in a changed catalytic mechanism of the enzyme. The possibility of utilizing this mutant for Michael-type additions was initially explored by quantum-chemical calculations on the reaction between acrolein and methanethiol in a model system. The model system was constructed on the basis of docking and molecular-dynamics simulations and was designed to simulate the catalytic properties of the active site. The catalytic system was explored experimentally with a range of different substrates. The kca values were found to be in the range of 10(-3) to 4 min(-1), similar to the values obtained with aldolase antibodies. The enzyme proficiency was 10(7). Furthermore, the Michael-type reactions followed saturation kinetics and were confirmed to take place in the enzyme active site.

Binding Sites↗

PubMed2007

Does tetracycline bind helix 2 of prion? An integrated spectroscopical and computational study of the interaction between the antibiotic and alpha helix 2 human prion protein fragments.

We demonstrate here that tetracycline (TC) can strongly interact (KD' = 189 +/- 7 nM) with model peptides derived from the C-terminal globular domain of the prion protein, hPrP [173-195], and that interaction concerns residues within the C-terminal half of the helix 2, a short region previously indicated as endowed with ambivalent conformational behavior and implicated in PrP conversion to the beta-sheet-rich, infective scrapie variant. Data have been confirmed by binding studies with the N-terminal truncated 180-195 variant that displays a dissociation constant of 483 +/- 30 nM. Remarkably, TC does not influence the structure of the N-terminally fluoresceinated peptides that both show alpha-helical conformations. Docking calculations and molecular dynamics simulations suggest a direct, strong interaction of the antibiotic with exposed side chain functional groups of threonines 190-193 on the solvent-exposed surface of helix 2.

Amino Acid Sequence↗

PubMed2003

Rational design of a lipase to accommodate catalysis of Baeyer-Villiger oxidation with hydrogen peroxide.

The mechanism and potential energy surface for the Baeyer-Villiger oxidation of acetone with hydrogen peroxide catalyzed by a Ser105-Ala mutant of Candida antarctica Lipase B has been determined using ab initio and density functional theories. Initial substrate binding has been studied using an automated docking procedure and molecular dynamics simulations. Substrates were found to bind to the active site of the mutant. The activation energy for the first step of the reaction, the nucleophilic attack of hydrogen peroxide on the carbonyl carbon of hydrogen peroxide, was calculated to be 4.4 kcal x mol(-1) at the B3LYP/6-31+G* level. The second step, involving the migration of the alkyl group, was found to be the rate-determining step with a computed activation energy of 19.9 kcal x mol(-1) relative the reactant complex. Both steps were found to be lowered considerably in the reaction catalyzed by the mutated lipase, compared to the uncatalyzed reaction. The first step was lowered by 36.0 kcal x mol(-1) and the second step by 19.5 kcal x mol(-1). The second step of the reaction, the rearrangement step, has a high barrier of 27.7 kcal x mol(-1) relative to the Criegee intermediate. This could lead to an accumulation of the intermediate. It is not clear whether this result is an artifact of the computational procedure, or an indication that further mutations of the active site are required. Figure Second TS (18TS) in the Baeyer-Villiger oxidation in a mutant of CALB. Distances in A

Binding Sites↗

PubMed2003

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

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

Computer Simulation↗

PubMed2006

Computational studies of H5N1 hemagglutinin binding with SA-alpha-2, 3-Gal and SA-alpha-2, 6-Gal.

For influenza H5N1 hemagglutinin, a switch from SA-alpha-2, 3-Gal to SA-alpha-2, 6-Gal receptor specificity is a critical step leading to the conversion from avian-to-human to human-to-human infection. Therefore, the understanding of the binding modes of SA-alpha-2, 3-Gal and SA-alpha-2, 6-Gal to H5N1 hemagglutinin will be very important for the examination of possible mutations needed for going from an avian to a human flu virus. Based on the available H5N1 hemagglutinin crystal structure, the binding profiles between H5N1 hemagglutinin and two saccharide ligands, SA-alpha-2, 3-Gal and SA-alpha-2, 6-Gal, were investigated by ab initio quantum mechanics, molecular docking, molecular mechanics, and molecular dynamics simulations. It was found that SA-alpha-2, 3-Gal has strong multiple hydrophobic and hydrogen bond interactions in its trans conformation with H5N1 hemagglutinin, whereas the SA-alpha-2, 6-Gal only shows weak interactions in a different conformation (cis type).

Computational Biology↗

PubMed2003

Identifying inhibitors of the SARS coronavirus proteinase.

The Severe Acute Respiratory Syndrome (SARS) is a serious respiratory illness that has recently been reported in parts of Asia and Canada. In this study, we use molecular dynamics (MD) simulations and docking techniques to screen 29 approved and experimental drugs against the theoretical model of the SARS CoV proteinase as well as the experimental structure of the transmissible gastroenteritis virus (TGEV) proteinase. Our predictions indicate that existing HIV-1 protease inhibitors, L-700,417 for instance, have high binding affinities and may provide good starting points for designing SARS CoV proteinase inhibitors.

Antiviral Agents↗

PubMed2005

Exploring blocker binding to a homology model of the open hERG K+ channel using docking and molecular dynamics methods.

Binding of blockers to the human voltage-gated hERG potassium channel is studied using a combination of homology modelling, automated docking calculations and molecular dynamics simulations, where binding affinities are evaluated using the linear interaction energy method. A homology model was constructed based on the available crystal structure of the bacterial KvAP channel and the affinities of a series of sertindole analogues predicted using this model. The calculations reproduce the relative binding affinities of these compounds very well and indicate that both polar interactions near the intracellular opening of the selectivity filter as well as hydrophobic complementarity in the region around F656 are important for blocker binding. These results are consistent with recent alanine scanning mutation experiments on the blocking of the hERG channel by other compounds.

Amino Acid Sequence↗

PubMed2001

Mechanisms of tetraethylammonium ion block in the KcsA potassium channel.

We report results from automated docking and microscopic molecular dynamics simulations of the tetraethylammonium (TEA) complexes with KcsA. Binding modes and energies for TEA binding at the external and internal sides of the channel pore are examined utilising the linear interaction energy method. Effects of the channel ion occupancy (based on our previous results for the ion permeation mechanisms) on the binding energies are considered. Calculations show that TEA forms stable complexes at both the external and internal entrances of the selectivity filter. Furthermore, the effects of the Y82V mutation are evaluated and the results show, in agreement with experimental data, that the mutant has a significantly reduced binding affinity for TEA at the external binding site, which is attributed to stabilising hydrophobic interactions between the ligand and the tyrosines.

Bacterial Proteins↗

PubMed2004

Cyclohexylcarbamic acid 3'- or 4'-substituted biphenyl-3-yl esters as fatty acid amide hydrolase inhibitors: synthesis, quantitative structure-activity relationships, and molecular modeling studies.

Fatty acid amide hydrolase (FAAH) is a promising target for modulating endocannabinoid and fatty acid ethanolamide signaling, which may have important therapeutic potential. We recently described a new class of O-arylcarbamate inhibitors of FAAH, including the cyclohexylcarbamic acid biphenyl-3-yl ester URB524 (half-maximal inhibitory concentration, IC(50) = 63 nM), which have significant anxiolytic-like properties in rats. In the present study, by introducing a selected group of substituents at the meta and para positions of the distal phenyl ring of URB524, we have characterized structure-activity profiles for this series of compounds and shown that introduction of small polar groups in the meta position greatly improves inhibitory potency. Most potent in the series was the m-carbamoyl derivative URB597 (4i, IC(50) = 4.6 nM). Furthermore, quantitative structure-activity relationship (QSAR) analysis of an extended set of meta-substituted derivatives revealed a negative correlation between potency and lipophilicity and suggested that small-sized substituents may undertake polar interactions with the binding pocket of the enzyme. Docking studies and molecular dynamics simulations, using the crystal structure of FAAH, indicated that the O-biphenyl scaffold of the carbamate inhibitors can be accommodated within a lipophilic region of the substrate-binding site, where their folded shape mimics the initial 10-12 carbon atoms of the arachidonyl moiety of anandamide (a natural FAAH substrate) and methyl arachidonyl fluorophosphonate (a nonselective FAAH inhibitor). Moreover, substituents at the meta position of the distal phenyl ring can form hydrogen bonds with atoms located on the polar section of a narrow channel pointing toward the membrane-associated side of the enzyme. The structure-activity characterization reported here should help optimize the pharmacodynamic and pharmacokinetic properties of this class of compounds.

Amidohydrolases↗

PubMed2026

In silico identification of DNMT1 inhibitors from the PlantCyc database through computational approach to assess the anti-cancer potential of nutraceutical compounds in breast cancer.

Breast cancer accounts for a disproportionate share of global cancer-related deaths, with 670,000 fatalities and 2.3 million new diagnoses recorded in women during 2022 alone. Existing treatment modalities carry considerable toxicity burdens, and resistance to available agents remains an unresolved clinical problem. DNA methyltransferase 1 (DNMT1), the enzyme chiefly responsible for maintaining genome-wide methylation patterns during DNA replication, has been mapped out as a high-value target in breast cancer because its dysregulation silences tumour suppressor genes through promoter hypermethylation. The present work involves hierarchical in silico workflow to screen 4549 plant-derived compounds from the PlantCyc database (v16.0.3) against the human DNMT1 catalytic domain (PDB ID: 4WXX). Ten top-scoring compounds were taken forward for molecular docking via AutoDock Vina; Quercetin and Kaempferol both recorded the highest binding affinities at -9.5 kcal/mol, Wogonin (-9.3 kcal/mol) and Xanthohumol (-8.1 kcal/mol) also emerged as strong binders. Pharmacokinetic evaluation using ADMET-AI confirmed that all 10 compounds met Lipinski's rule of five, with human intestinal absorption values at or above 0.98. Wogonin and Xanthohumol were selected for a 100 ns all-atom molecular dynamics (MD) simulation in GROMACS due to their well-rounded ADMET profiles and limited existing data on their specific interactions with DNMT1 in breast cancer. Across all measured trajectory metrics, backbone RMSD, residue fluctuation, radius of gyration, solvent-accessible surface area, and intermolecular hydrogen bond count, Wogonin formed a more stable, compact complex. These findings suggest that Wogonin and Xanthohumol are non-toxic nutraceutical candidates suitable for DNMT1 targeted epigenetic therapy, with computational foundation strong enough to facilitate future in vitro and in vivo validation work.

Humans↗

PubMed1990

Molecular modelling study of antigen binding to oxazolone-specific antibodies: the Ox1 idiotypic IgG and its mature variant with increased affinity to 2-phenyloxazolone.

Structural models of the variable domains of the murine anti-2-phenyloxazolone IgG (Ox1 idiotype) and its somatic variant, which has higher affinity to the hapten 2-phenyloxazolone, were constructed by computer-aided model building using known structures of highly homologous immunoglobulins as templates. Molecular dynamics simulations were used to dock the hapten between the VL and VH domains. The hapten is predicted to bind to slightly different sites in the two models. Hypotheses concerning the role of a number of preferred mutations in anti-oxazolone variants are presented. These can be tested by mutagenesis and crystallography. In particular, the higher binding affinities of the different antibody variants are shown to correlate with better complementarity of electrostatics. The molecular dynamic simulations also suggest that two mobile tryptophans at the mouth of the pocket may play an important role in the binding of hapten.

Amino Acid Sequence↗