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 73 records · Page 4Linked to original sources

PubMed2003

Targeting aurora2 kinase in oncogenesis: a structural bioinformatics approach to target validation and rational drug design.

The aurora kinases are a novel oncogenic family of mitotic serine/threonine kinases (S/T kinases) that are overexpressed in a number of solid tumors, including pancreas and colorectal cancer. A PSI-BLAST search [National Center for Biotechnology Information (NCBI)] with the sequence of the S/T kinase domain of human aurora1 kinase [also known as AUR1, ARK2, AIk2, AIM-1, and STK12] and human aurora2 kinase (also known as AUR2, ARK1, AIK, BTAK, and STK15) showed a high sequence similarity to the three-dimensional structures of bovine cAMP-dependent kinase [Brookhaven Protein Data Bank code 1CDK], murine cAMP-dependent kinase (1APM), and Caenorhabditis elegans twitchin kinase (1KOA). When the aurora1 or aurora2 sequence was input into the tertiary structure prediction programs THREADER and 3D-PSSM (three-dimensional position-sensitive scoring matrix), the top structural matches were 1CDK, 1APM, and 1KOA, confirming that these domains are structurally conserved. The structural models of aurora1 and aurora2 were built using 1CDK as the template structure. Molecular dynamics and docking simulations, targeting the ATP binding site of aurora2 with adenylyl imidodiphosphate (AMP-PNP), staurosporine, and six small molecular S/T kinase inhibitors, identified active-site residues that interact with these inhibitors differentially. The docked structures of the aurora2-AMP-PNP and aurora2-staurosporine complexes indicated that the adenine ring of AMP-PNP and the indolocarbazole moiety of staurosporine have similar positions and orientations and provided the basis for the docking of the other S/T kinase inhibitors. Inhibitors with isoquinoline and quinazoline moieties were recognized by aurora2 in which H-89 and 6,7-dimethoxyquinazoline compounds exhibited high binding energies compared with that of staurosporine. The calculated binding energies for the docked small-molecule inhibitors were qualitatively consistent with the IC(50) values generated using an in vitro kinase assay. The aurora2 structural model provides a rational basis for site-directed mutagenesis of the active site; design of novel H-89, staurosporine, and quinazoline analogues; and the screening of the available chemical database for the identification of other novel, small-molecular entities.

Amino Acid Sequence↗

PubMed2003

Studying the roles of W86, E202, and Y337 in binding of acetylcholine to acetylcholinesterase using a combined molecular dynamics and multiple docking approach.

A combined molecular dynamics simulation and multiple ligand docking approach is applied to study the roles of the anionic subsite residues (W86, E202, Y337) in the binding of acetylcholine (ACh) to acetylcholinesterase (AChE). We find that E202 stabilizes docking of ACh via electrostatic interactions. However, we find no significant electrostatic contribution from the aromatic residues. Docking energies of ACh to mutant AChE show a more pronounced effect because of size/shape complementarity. Mutating to smaller residues results in poorer binding, both in terms of docking energy and statistical docking probability. Besides separating out electrostatics by turning off the partial charges from each residue and comparing it with the native, the mutations in this study are W86F, W86A, E202D, E202Q, E202A, Y337F, and Y337A. We also find that all perturbations result in a significant reduction in binding of extended ACh in the catalytically productive orientation. This effect is primarily caused by a small shift in preferred position of the quaternary tail.

Acetylcholine↗

PubMed2006

Possible dynamic anchor points in a benzoxazinone derivative-human oxytocin receptor system--a molecular docking and dynamics calculation.

In this study, we performed a molecular docking and dynamics simulation for a benzoxazinone-human oxytocin receptor system to determine the possible hydrophobic and electrostatic interaction points in the dynamic complex. After the homology modeling, the ligand was docked into the putative active using AutoDock 3.05. After the application of energetic and structural filters, the complexes obtained were further refined with a simulated annealing protocol (AMBER8) to remove steric clashes. Three complexes were selected for subjection to the molecular dynamics simulation (5 ns), and the results on the occurrence of average anchor points showed a stable complex between the benzoxazinone derivative and the receptor. The complex could be used as a good starting point for further analysis with site-directed mutagenesis, or further computational research.

Animals↗

PubMed2002

Studying enzyme binding specificity in acetylcholinesterase using a combined molecular dynamics and multiple docking approach.

A combined molecular dynamics simulation and multiple ligand docking approach is applied to study the binding specificity of acetylcholinesterase (AChE) with its natural substrate acetylcholine (ACh), a family of substrate analogues, and choline. Calculated docking energies are well correlated to experimental k(cat)/K(M) values, as well as to experimental binding affinities of a related series of TMTFA inhibitors. The "esteratic" and "anionic" subsites are found to act together to achieve substrate binding specificity. We find that the presence of ACh in the active site of AChE not only stabilizes the setup of the catalytic triad but also tightens both subsites to achieve better binding. The docking energy gained from this induced fit is 0.7 kcal/mol for ACh. For the binding of the substrate tailgroup to the anionic subsite, both the size and the positive charge of the tailgroup are important. The removal of the positive charge leads to a weaker binding of 1 kcal/mol loss in docking energy. Substituting each tail methyl group with hydrogen results in both an incremental loss in docking energy and also a decrease in the percentage of structures docked in the active site correctly set up for catalysis.

Acetophenones↗

PubMed2025

Exploring effector protein dynamics and natural fungicidal potential in rice blast pathogen Magnaporthe oryzae.

Rice blast, caused by Magnaporthe oryzae, is one of the most destructive fungal diseases in rice, resulting in major economic losses worldwide. Genetic and genomic studies have identified key genes and proteins, such as AvrPik variants and MAX proteins, that are crucial for the pathogen's virulence. These effector proteins interact with specific alleles of the Pik gene family on rice chromosome 11, modulating the host's immune response. In this study, we investigated 35 plant-derived metabolites known for their antifungal properties as potential fungicides against M. oryzae. Using molecular docking, we identified Hecogenin and Cucurbitacin E as strong binders to MAX40 and APIKL2A proteins, which are essential for the fungus's immune evasion and pathogenicity. Molecular dynamics simulations further confirmed that these compounds form stable, strong interactions with the target proteins, validating their potential as therapeutic agents. Additionally, the compounds were evaluated based on Lipinski's rule of five and toxicity predictions, indicating their suitability for agricultural use. These results suggest that Hecogenin and Cucurbitacin E could serve as promising lead candidates in the development of novel fungicides for rice blast, offering new strategies for crop protection and sustainable agricultural practices.

Oryza↗

PubMed2026

Integrative computational analysis combining network pharmacology, regulatory network modeling, and molecular dynamics reveals the mechanisms of Quanshen compound in ITP.

UNLABELLED: Immune thrombocytopenia (ITP) is a hemorrhagic disorder caused by immune dysfunction. Quanshen Compound (QSC) is an in-house preparation developed by the Uyghur Hospital in Hotan Prefecture. This study primarily investigates and validates the potential pharmacological basis and mechanism of action of QSC in modulating immune thrombopoiesis. Based on the multi-database screening of the QSC and the related targets of ITP, the intersection was obtained to construct a protein-protein interaction (PPI) network and screen the core targets; the intersection targets were analyzed for gene ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis using R packages; a component-target-pathway network was constructed to screen the key active components and their mechanisms of action. At the same time, the TF-mRNA-miRNA regulatory network of the core targets was constructed, and chromosome localization and subcellular localization analysis were performed; further, the binding stability of key components and core targets was verified through molecular docking and molecular dynamics simulation. A total of 227 potential target sites were screened out, among which TNF, IL6, AKT1, TP53 and IL1B were the core targets. The enrichment results indicated that these intersecting target sites mainly participated in inflammatory responses, immune regulation and hemostasis-related biological processes, and were significantly enriched in the PI3K-Akt signaling pathway, Toll-like receptor signaling pathway, Th17 cell differentiation and PD-1/PD-L1 signaling pathway. The core target TF-mRNA-miRNA regulatory network contained 184 nodes and 200 edges, suggesting that the core targets were subject to multi-level regulation. Molecular docking results showed that the main active components had good binding activity with the core targets, and molecular dynamics simulation further verified the stability of the complex. QSC may improve ITP through a multi-component, multi-target, and multi-pathway synergistic mechanism involving key targets such as TNF, IL6, AKT1, TP53, and IL1B, as well as the PI3K-Akt signaling pathway. These findings provide new insights into the potential therapeutic mechanisms of QSC against ITP and warrant further experimental validation. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at https://doi.org/10.1007/s40203-026-00718-0.

Immune thrombocytopenia↗

PubMed2004

Development of the force field parameters for phosphoimidazole and phosphohistidine.

Phosphorylation of histidine-containing proteins is a key step in the mechanism of many phosphate transfer enzymes (kinases, phosphatases) and is the first stage in a wide variety of signal transduction cascades in bacteria, yeast, higher plants, and mammals. Studies of structural and dynamical aspects of such enzymes in the phosphorylated intermediate states are important for understanding the intimate molecular mechanisms of their functioning. Such information may be obtained via molecular dynamics and/or docking simulations, but in this case appropriate force field parameters for phosphohistidine should be explicitly defined. In the present article we describe development of the GROMOS96 force field parameters for phosphoimidazole molecule--a realistic model of the phosphohistidine side chain. The parameterization is based on the results of ab initio quantum chemical calculations with subsequent refinement and testing using molecular mechanics and molecular dynamics simulations. The set of force constants and equilibrium geometry is employed to derive force field for the phosphohistidine moiety. Resulting parameters and topology are incorporated into the molecular modeling package GROMACS and used in molecular dynamics simulations of a phosphohistidine-containing protein in explicit solvent.

Algorithms↗

PubMed2004

Homology modeling, force field design, and free energy simulation studies to optimize the activities of histone deacetylase inhibitors.

As an effort to develop therapeutics for cancer treatments, a number of effective histone deacetylase inhibitors with structural diversity have been discovered. To gain insight into optimizing the activity of an identified lead compound, a computational protocol sequentially involving homology modeling, docking experiments, molecular dynamics simulation, and free energy perturbation calculations was applied for rationalizing the relative activities of known histone deacetylase inhibitors. With the newly developed force field parameters for the coordination environment of the catalytic zinc ion in hand, the computational strategy proved to be successful in predicting the rank orders for 12 derivatives of three hydroxamate-based inhibitor scaffolds with indole amide, pyrrole, and sulfonamide moieties. The results showed that the free energy of an inhibitor in aqueous solution should be an important factor in determining the binding free energy. Hence, in order to enhance the inhibitory activity by adding or substituting a chemical group, the increased stabilization in solution due to the structural changes must be overcome by a stronger enzyme-inhibitor interaction. It was also found that to optimize inhibitor potency, the hydrophobic head of an inhibitor should be elongated or enlarged so that it can interact with Pro29 and His28 that are components of the flexible loop at the top of the active site.

Amino Acid Sequence↗

PubMed2005

Modeling Polarization in Proteins and Protein-ligand Complexes: Methods and Preliminary Results.

This chapter discusses methods for modeling electronic polarization in proteins and protein-ligand complexes. Two different approaches are considered: explicit incorporation of polarization into a molecular mechanics force field and the use of mixed quantum mechanics/molecular mechanics methods to model polarization in a restricted region of the protein or protein-ligand complex. A brief description is provided of the computational methodology and parameterization protocols and then results from two preliminary studies are presented. The first study employs quantum mechanics/molecular mechanics (QM/MM) methods to improve the accuracy of protein-ligand docking; here, incorporation of polarization is shown to dramatically improve the robustness of the accuracy of structural prediction of the protein-ligand docking by enabling qualitative improvement in the selection of the correct hydrogen bonding patterns of the docked ligand. The second study discusses a 2-ns simulation of bovine pancreatic trypsin inhibitor (BPTI) in water using a variety of fixed charge and polarizable models for both the protein and the solvent, analyzing observed root mean square deviations (RMSD), intraprotein hydrogen bonding, and water structure and dynamics. All of these efforts are in a relatively early stage of development, the results are encouraging in that stable methods have been developed, and significant effects of polarization are seen and (in the case of the QM/MM-based docking) improvements have been validated as compared to experiment. With regard to accuracy and robustness of full simulations, a great deal more work needs to be done to quantitate and improve the present models.

Animals↗

PubMed2004

Modelling and simulation: a computational perspective in anticancer drug discovery.

The availability of high-quality molecular graphics tools in the public domain is changing the way macromolecular structure is perceived by researchers, educators and students alike. Computational methods have become increasingly important in a number of areas such as comparative or homology modelling, functional site location, characterization of ligand-binding sites in proteins, docking of small molecules into protein binding sites, protein-protein docking, and molecular dynamics simulations. The results obtained yield information that sometimes is beyond current experimental possibilities and can be used to guide and improve a vast array of experiments. On the basis of our improved level of understanding of molecular recognition and the widespread availability of target structures, it is reasonable to assume that computational methods will continue aiding not only in the design and interpretation of hypothesis-driven experiments in the field of cancer research but also in the rapid generation of new hypotheses.

Amino Acid Sequence↗

PubMed1999

Catechol(amine)s as probes of lactoperoxidase catalytic site structure: spectroscopic and modeling studies.

Binding affinities to lactoperoxidase (LPO) of a homologous series of substituted catechol(amine)s [such as catechol, 4-methylcatechol, 3,4-dihydroxybenzoic acid, 3,4-dihydroxyphenylacetic acid, 3-(3,4-dihydroxyphenyl)propionic acid; dopamine, noradrenaline, adrenaline; L-3,4-dihydroxyphenylalanine] were studied by UV-visible spectroscopy and docking simulations. Dissociation constant (Kd) values were calculated by direct fitting of the experimental data and fall in a range of 3-95 mM. Thermodynamic parameters are comparable with those reported for the interaction of LPO with p-substituted phenols, suggesting a similar general mode of binding. Furthermore, the relative contributions to binding energy, described by the unimolecular constant Ku, show that interaction between protein and ligands originates from a relatively large number of groups. Docking and molecular dynamics simulations, in agreement with experimental evidence, predict that the substrate is localized into the access channel in the vicinity of heme distal pocket. This channel is characterized by a hydrophobic patch (six Phe residues) and by a charged contribution (two Glu and one His residues). All of the substrates, except caffeic acid, may approach the protein active site. Positively charged Arg372 acts as a gate above the heme distal pocket and seems to address substrate orientation in relation to the side-chain terminal group.

3,4-Dihydroxyphenylacetic Acid↗

PubMed2026

Resveratrol Attenuates Gemcitabine Resistance in Hepatocellular Carcinoma Cells by Inhibiting Thymidylate Synthase.

BACKGROUND: Hepatocellular carcinoma (HCC) is a leading cause of cancer death worldwide. Gemcitabine (Gem) is a commonly used drug against HCC, but its efficacy is limited by the development of resistance. Resveratrol (Res), a natural polyphenol with antitumor activity, may reverse Gem resistance in HCC, although the mechanism remains unclear. METHODS: The effects of Res on the proliferation, apoptosis, cell cycle, and invasion of Hep3B and HuH-7 cells were assessed via cell counting kit-8 (CCK-8), clonogenic, flow cytometry, and Transwell assays, respectively. Potential Res targets were predicted by network pharmacology, and markers of HCC prognosis were identified from the cancer genome atlas (TCGA) data. The interaction between Res and thymidylate synthase (TYMS) was validated by molecular docking and dynamics simulation. A Gem-resistant HuH-7 cell line (HuH-7/GR) was established, and when these cells were treated with Res combined with Gem, the effect on Gem sensitivity was detected by CCK-8 assay, clonogenic assay, and flow cytometry. Finally, a subcutaneous nude mouse model of HCC was used to evaluate the in vivo effects of Res combined with Gem. RESULTS: Res inhibited HCC cell proliferation, induced apoptosis and G2/M arrest, and suppressed invasion in a concentration-dependent manner. Network pharmacology and TCGA analysis identified TYMS as an important target gene for Res. TYMS was highly expressed in HCC tissues and correlated with poor prognosis. Res treatment reduced TYMS expression, while molecular docking and simulation showed stable binding of Res to TYMS. TYMS levels were elevated in HuH-7/GR resistant cells. Res combined with Gem was found to reverse drug resistance, inhibit proliferation and colony formation, and induce apoptosis. The Res + Gem combination group showed the smallest tumor volume in the in vivo model. CONCLUSION: By attenuating Gem resistance through TYMS inhibition, Res holds promise as a clinically viable adjunct to Gem-based chemotherapy, offering a potential strategy to improve outcomes in HCC patients.

Resveratrol↗

PubMed2003

Adenosine A2A-dopamine D2 receptor-receptor heteromerization: qualitative and quantitative assessment by fluorescence and bioluminescence energy transfer.

There is evidence for strong functional antagonistic interactions between adenosine A2A receptors (A2ARs) and dopamine D2 receptors (D2Rs). Although a close physical interaction between both receptors has recently been shown using co-immunoprecipitation and co-localization assays, the existence of a A2AR-D2R protein-protein interaction still had to be demonstrated in intact living cells. In the present work, fluorescence resonance energy transfer (FRET) and bioluminescence resonance energy transfer (BRET) techniques were used to confirm the occurrence of A2AR-D2R interactions in co-transfected cells. The degree of A2AR-D2R heteromerization, measured by BRET, did not vary after receptor activation with selective agonists, alone or in combination. BRET competition experiments were performed using a chimeric D2R-D1R in which helices 5 and 6, the third intracellular loop (I3), and the third extracellular loop (E3) of the D2R were replaced by those of the dopamine D1 receptor (D1R). Although the wild type D2R was able to decrease the BRET signal, the chimera failed to achieve any effect. This suggests that the helix 5-I3-helix 6-E3 portion of D2R holds the site(s) for interaction with A2AR. Modeling of A2AR and D2R using a modified rhodopsin template followed by molecular dynamics and docking simulations gave essentially two different possible modes of interaction between D2R and A2AR. In the most probable one, helix 5 and/or helix 6 and the N-terminal portion of I3 from D2R approached helix 4 and the C-terminal portion of the C-tail from the A2AR, respectively.

Binding Sites↗

PubMed2004

Severe acute respiratory syndrome coronavirus 3C-like proteinase N terminus is indispensable for proteolytic activity but not for enzyme dimerization. Biochemical and thermodynamic investigation in conjunction with molecular dynamics simulations.

Severe acute respiratory syndrome (SARS) coronavirus is a novel human coronavirus and is responsible for SARS infection. SARS coronavirus 3C-like proteinase (SARS 3CL(pro)) plays key roles in viral replication and transcription and is an attractive target for anti-SARS drug discovery. In this report, we quantitatively characterized the dimerization features of the full-length and N-terminal residues 1-7 deleted SARS 3CL(pro)s by using glutaraldehyde cross-linking SDS-PAGE, size-exclusion chromatography, and isothermal titration calorimeter techniques. Glutaraldehyde cross-linking SDS-PAGE and size-exclusion chromatography results show that, similar to the full-length SARS 3CL(pro), the N-terminal deleted SARS 3CL(pro) still remains a dimer/monomer mixture within a wide range of protein concentrations. Isothermal titration calorimeter determinations indicate that the equilibrium dissociation constant (K(d)) of the N-terminal deleted proteinase dimer (262 microm) is very similar to that of the full-length proteinase dimer (227 microm). Enzymatic activity assay using the fluorescence resonance energy transfer method reveals that N-terminal deletion results in almost complete loss of enzymatic activity for SARS 3CL(pro). Molecular dynamics and docking simulations demonstrate the N-terminal deleted proteinase dimer adopts a state different from that of the full-length proteinase dimer, which increases the angle between the two protomers and reduces the binding pocket that is not beneficial to the substrate binding. This conclusion is verified by the surface plasmon resonance biosensor determination, indicating that the model substrate cannot bind to the N-terminal deleted proteinase. These results suggest the N terminus is not indispensable for the proteinase dimerization but may fix the dimer at the active state and is therefore vital to enzymatic activity.

Binding Sites↗

PubMed2006

Subtype selectivity and flexibility of ionotropic glutamate receptors upon antagonist ligand binding.

The binding modes of a set of known ionotropic glutamate receptor antagonist-ligands have been studied using homology modeling, molecular docking, molecular dynamics (MD) simulations and ab initio quantum mechanical calculations. The core structure of the studied ligands is the decahydroisoquinoline ring, which has a carboxylic acid group at position three and different negatively-charged substituents (R) at position six. The binding affinities of these molecules have been reported earlier. From the current study, the carboxylate group of the decahydroisoquinoline ring hydrogen bonds with Arg485, the amino group with Pro478 and Thr480, and the negatively charged substituent R interacts with the positively charged N-terminus of helix-F. The subtype selectivity of these ligands seems to be strongly dependent on the amino acid at position 650 (GluR2: leucine, GluR5: valine), which affects the conformation of the ligand and ligand-receptor interactions, but depends considerably on the size of the R-group of the ligand. In addition, the MD simulations also revealed that the relative positions of the S1 and S2 domains can alter significantly showing different "closure" and "rotational movements" depending on the antagonist-ligand that is bound. Accordingly, molecular docking of antagonist ligands into static crystal structures cannot sufficiently explain ligand binding and subtype selectivity.

Amino Acid Sequence↗

PubMed2026

Toxicological effects of propyl 4-hydroxybenzoate on gallstone pathogenesis: An integrated mendelian randomization, network toxicology, and experimental study.

BACKGROUND: Gallstone disease is a prevalent digestive disorder with substantial global socioeconomic burden. Propyl 4-hydroxybenzoate (PP), a widely used paraben preservative, exhibits potential metabolic and hepatic toxicity, yet its role in gallstone pathogenesis remains unclear. This study aimed to explore the causal association between PP exposure and gallstone formation and the underlying mechanism. METHODS: Two-sample Mendelian randomization (MR) was performed using genome-wide association study (GWAS) data. Network toxicology, molecular docking, and molecular dynamics simulation were applied to screen for core targets. In vivo experiments, transcriptome sequencing, Western blot (WB), and ELISA were conducted for mechanistic validation. RESULTS: MR confirmed a causal link between circulating PP levels and an elevated risk of gallstones (P < 0.05), with AKT1 identified as the key target. In mice, PP aggravated gallstone formation by activating the AKT1-NF-κB-CXCL1 pathway, enhancing hepatic inflammation and neutrophil extracellular traps (NETs) formation; these effects were reversed by AKT inhibition. CONCLUSION: PP promotes gallstone formation via the AKT1-NF-κB-CXCL1-NETs axis. Our findings highlight PP as an environmental risk factor for gallstones, providing novel insights into their prevention and targeted therapy.

Animals↗

PubMed2006

Experimental/theoretical electrostatic properties of a styrylquinoline-type HIV-1 integrase inhibitor and its progenitors.

We have established that polyhydroxylated styrylquinolines are potent inhibitors of HIV-1 integrase (IN). Among them, we have identified (E)-8-hydroxy-2-[2-(4,5-dihydroxy-3-methoxyphenyl)-ethenyl]-7-quinolinecarboxylic acid (1) as a promising lead. Previous molecular dynamics simulations and docking procedures have shown that the inhibitory activity involves one or two metal cations (Mg2+), which are present in the vicinity of the active center of the enzyme. However, such methods are generally based on a force-field approach and still remain not as reliable as ab initio calculations with extended basis sets on the whole system. To go further in this area, the aim of the present study was to evaluate the predictive ability of the electron density and electrostatic properties in the structure-activity relationships of this class of HIV-1 antiviral drugs. The electron properties of the two chemical progenitors of 1 were derived from both high-resolution X-ray diffraction experiments and ab initio calculations. The twinning phenomenon and solvent disorder were observed during the crystal structure determination of 1. Molecule 1 exhibits a planar s-trans conformation, and a zwitterionic form in the crystalline state is obtained. This geometry was used for ab initio calculations, which were performed to characterize the electronic properties of 1. The electron densities, electrostatic potentials, and atomic charges of 1 and its progenitors are here compared and analyzed. The experimental and theoretical deformation density bond peaks are very comparable for the two progenitors. However, the experimental electrostatic potential is strongly affected by the crystal field and cannot straightforwardly be used as a predictive index. The weak difference in the theoretical electron densities between 1 and its progenitors reveals that each component of 1 conserves its intrinsic properties, an assumption reinforced by a 13C NMR study. This is also shown through an excellent correlation of the atomic charges for the common fragments. The electrostatic potential minima in zwitterionic and nonzwitterionic forms of 1 are discussed in relation with the localization of possible metal chelation sites.

Carbon Isotopes↗

PubMed2000

Shape complementarity, binding-site dynamics, and transition state stabilization: a theoretical study of Diels-Alder catalysis by antibody 1E9.

Antibody 1E9 is a protein catalyst for the Diels-Alder reaction between tetrachlorothiophene dioxide and N-ethylmaleimide. Quantum mechanical calculations have been employed to study the 1E9-catalyzed Diels-Alder reaction in the gas phase. The transition states and intermediates were all determined at the B3LYP/6-31G*//HF/6-31G* level. The cycloaddition step is predicted to be rate-determining, and the endo reaction pathway is strongly favored. Binding of the reactants and the transition states to antibody 1E9 was investigated by docking and molecular dynamics simulations. The linear interaction energy (LIE) method was adopted to estimate the free energy barrier of the 1E9-catalyzed Diels-Alder reaction. The catalytic efficiency of antibody 1E9 is achieved by enthalpic stabilization of the transition state, near-perfect shape complementarity of the hydrophobic binding site for the transition state, and a strategically placed hydrogen-bonding interaction.

Antibodies, Catalytic↗