SEARCH · PubMed Health

Results for “molecular docking studies”

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 235 records · Page 13Linked to original sources

PubMed2026

Proteomic regulation of anti-proliferative and anti-migratory activity by potent phytochemicals from Pistacia integerrima J.L. Steward Ex Brandis via PI3K, AKT1, and KRAS for Lung Cancer.

BACKGROUND: Non-small cell lung cancer (NSCLC) is the leading cause of mortality worldwide and remains a major therapeutic challenge due to high metastasis, drug resistance and limited treatments. Pistacia integerrima J.L. Steward Ex Brandis (PI) consists of flavonoids, steroids, terpenoids and phenolic compounds reported for pharmacological activities. The efficacy of potent bioactives from P. integerrima may be ascertained employing cytotoxic, antiproliferative, anti-migratory, and anti-metastatic evaluations in A549 NSCLC cells with proteomic profiling, molecular docking, and dynamics simulation study. METHODS AND RESULTS: PI EtAc produced significant dose-dependent cytotoxicity in A549 cells (100 µg/mL, p < 0.0001 in the MTT assay. There was a pronounced decrease in colony formation after treatment with EtAc, with 18.41% (p < 0.002), and markedly. Furthermore. PI EtAC markedly inhibited cell migration emphasized by wound healing and Transwell migration (p < 0.01) assays, indicating reduced metastatic migratory potential. Proteomic analysis demonstrated significant downregulation of Endoglin (CD105), KLK5 and MMP-2, indicating suppression of angiogenic and metastatic signalling pathways in the Human XL Oncology protein array. The interaction of major PI phytochemicals with key NSCLC-associated targets was recorded in Molecular docking, revealing favourable binding affinities of kaempferol, β- sitosterol, luteolin, and quercetin towards several oncogenic targets, including AKT1(- 7.6 kcal/mol), PI3K(- 9.4 kcal/mol), KRAS (- 8.5 kcal/mol) and MMP9 (- 8.1 kcal/mol). Molecular dynamics simulation confirmed the structural stability of the kaempferol -AKT1 complex throughout the 100 ns simulation. CONCLUSION: Pistacia integerrima bioactives exhibited significant anti-proliferative, anti-migratory, and anti-metastatic activities in vitro, which may provide scientific rationale identifying newer promising candidates for NSCLC.

Humans↗

PubMed2001

DNA repair mechanism by photolyase: electron transfer path from the photolyase catalytic cofactor FADH(-) to DNA thymine dimer.

Photolyase is an enzyme that catalyses photorepair of thymine dimers in UV damaged DNA by electron transfer reaction. The structure of the photolyase/DNA complex is unknown at present. Using crystal structure coordinates of the substrate-free enzyme from E. coli, we have recently built a computer molecular model of a thymine dimer docked to photolyase catalytic site and studied molecular dynamics of the system. In this paper, we present analysis of the electronic coupling and electron transfer pathway between the catalytic cofactor FADH(-) and the pyrimidine dimer by the method of interatomic tunneling currents. Electronic structure is treated in the extended Hückel approximation. The root mean square transfer matrix element is about 6 cm(-1), which is consistent with the experimentally determined rate of transfer. We find that electron transfer mechanism responsible for the repair utilizes an unusual folded conformation of FADH(-) in photolyases, in which the isoalloxazine ring of the flavin and the adenine are in close proximity, and the peculiar features of the docked orientation of the dimer. The tunneling currents show explicitly that despite of the close proximity between the donor and acceptor complexes, the electron transfer mechanism between the flavin and the thymine bases is not direct, but indirect, with the adenine acting as an intermediate. These calculations confirm the previously made conclusion based on an indirect evidence for such mechanism.

Animals↗

PubMed1996

Gbeta subunit interacts with a peptide encoding region 956-982 of adenylyl cyclase 2. Cross-linking of the peptide to free Gbetagamma but not the heterotrimer.

The region encoded by amino acids 956-982 of adenylyl cyclase 2 is important for Gbetagamma stimulation. Interactions of a peptide encoding the 956-982 region of adenylyl cyclase 2 (QEHAQEPERQYMHIGTMVEFAYALVGK (QEHA peptide)) with Gbetagamma subunits were studied. QEHA peptide was covalently attached to beta subunit of free Gbetagamma by the cross-linker N-succinimidyl(4-iodoacetyl)aminobenzoate. Cross-linking was proportional to the amount of QEHA peptide added; other control peptides cross-linked minimally. When Go was used, very little cross-linking was observed with GDP and EDTA, but upon activation by guanosine 5'-3-O-(thio)triphosphate and Mg2+, specific cross-linking of the QEHA peptide to Gbeta was observed. We conclude that beta subunits of G proteins contain effector interaction domains that are occluded by Galpha subunits in the heterotrimer. Molecular modeling studies used to dock the QEHA peptide on to Gbeta indicate that amino acids 75-165 of Gbeta may be involved in effector interactions.

Adenylyl Cyclases↗

PubMed2025

Integrated data mining and network pharmacology to explore the prescription patterns from a senior TCM oncologist's clinical practice in treating chemotherapy-induced hand-foot syndrome.

Hand-foot syndrome (HFS) is a common and refractory adverse effect of chemotherapy lacking specific therapeutic strategies currently. Traditional Chinese medicine (TCM) has shown empirical efficacy in clinical HFS management. This study integrated data mining and network pharmacology to systematically elucidate the medication principles and molecular mechanisms underlying Professor Gang Xie's prescriptions for HFS. All medical records from Professor Xie's specialist clinic (January 2020 to March 2025) were retrospectively collected and standardized in Excel. Prescriptions were analyzed through frequency statistics, association and clustering. Active ingredients of core herb pairs and their disease-related targets were identified using TCMSP, HERB, GeneCards, PharmGKB and GEO databases. Protein-protein interaction (PPI) networks, gene ontology (GO), and Kyoto encyclopedia of genes and genomes (KEGG) pathway analyses were performed. Molecular docking validated interactions between key bioactive compounds and targets. This study involved 217 prescriptions containing 150 herbs. Core herb combinations comprised Radix Astragali (Huangqi), Poria (Fuling), and Radix Pseudostellariae (Taizishen), predominantly classified as spleen-tonifying agents with warm properties, targeting lung, spleen, and stomach meridians. Network analysis identified 67 bioactive compounds and 899 disease targets. Quercetin, kaempferol, acacetin and luteolin were identified the key ingredients. The core targets (TP53, STAT3, PIK3CA, HSP90AA1, AKT1, CTNNB1, PI3KR1, MAPK1) were enriched in MAPK and PI3K-Akt signaling pathways. Molecular docking confirmed strong binding affinity between key compounds and targets. Professor Xie's therapeutic strategy for HFS emphasizes "spleen fortification, phlegm elimination, and stasis resolution." The core herb combination likely exerts anti-HFS effects via modulation of MAPK and PI3K-Akt pathways, providing a pharmacological basis for TCM-driven HFS management.

Network Pharmacology↗

PubMed2005

Structure-activity relationships for 1-alkyl-3-(1-naphthoyl)indoles at the cannabinoid CB(1) and CB(2) receptors: steric and electronic effects of naphthoyl substituents. New highly selective CB(2) receptor agonists.

In an effort to improve indole-based CB(2) cannabinoid receptor ligands and also to develop SAR for both the CB(1) and CB(2) receptors, 47 indole derivatives were prepared and their CB(1) and CB(2) receptor affinities were determined. The indole derivatives include 1-propyl- and 1-pentyl-3-(1-naphthoyl)indoles both with and without a 2-methyl substituent. Naphthoyl substituents include 4- and 7-alkyl groups as well as 2-, 4-, 6-, 7-methoxy and 4-ethoxy groups. The effects of these substituents on receptor affinities are discussed and structure-activity relationships are presented. In the course of this work three new highly selective CB(2) receptor agonists were identified, 1-propyl-3-(4-methyl-1-naphthoylindole (JWH-120), 1-propyl-2-methyl-3-(6-methoxy-1-naphthoylindole (JWH-151), and 1-pentyl-3-(2-methoxy-1-naphthoylindole (JWH-267). GTPgammaS assays indicated that JWH-151 is a full agonist at CB(2), while JWH-120 and JWH-267 are partial agonists. Molecular modeling and receptor docking studies were carried out on a set of 3-(4-propyl-1-naphthoyl)indoles, a set of 3-(6-methoxy-1-naphthoyl)indoles and the pair of N-pentyl-3-(2-methoxy-1-naphthoyl)indoles. Docking studies indicated that the CB(1) receptor affinities of these compounds were consistent with their aromatic stacking interactions in the aromatic microdomain of the CB(1) receptor.

Alkylation↗

PubMed2003

3-Indolyl-1-naphthylmethanes: new cannabimimetic indoles provide evidence for aromatic stacking interactions with the CB(1) cannabinoid receptor.

A series of 1-pentyl-1H-indol-3-yl-(1-naphthyl)methanes (9-11) and 2-methyl-1-pentyl-1H-indol-3-yl-(1-naphthyl)methanes (12-14) have been synthesized to investigate the hypothesis that cannabimimetic 3-(1-naphthoyl)indoles interact with the CB(1) receptor by hydrogen bonding to the carbonyl group. Indoles 9-11 have significant (K(i)=17-23nM) receptor affinity, somewhat less than that of the corresponding naphthoylindoles (5, 15, 16). 2-Methyl-1-indoles 12-14 have little affinity for the CB(1) receptor, in contrast to 2-methyl-3-(1-naphthoyl)indoles 17-19, which have affinities comparable to those of 5, 15, 16. A cannabimimetic indene hydrocarbon (26) was synthesized and found to have K(i)=26+/-4nM. Molecular modeling and receptor docking studies of naphthoylindole 16, its 2-methyl congener (19) and indolyl-1-naphthylmethanes 11 and 14, combined with the receptor affinities of these cannabimimetic indoles, strongly suggest that these cannabinoid receptor ligands bind primarily by aromatic stacking interactions in the transmembrane helix 3-4-5-6 region of the CB(1) receptor.

Cannabinoids↗

PubMed1995

Characterization of dynamic state inhibitors of HIV-1 protease.

The inhibitory behavior of simple N-protected amino acids containing a dysfunctional enol group was analyzed. Pseudo-Hill coefficients below 0.5 were obtained for the enol compounds. The active moiety of the enols described herein is involved in an equilibrium between keto and enol forms. The enol in its E form interacts strongly with the enzyme and we suggest that the Z form is expelled from the active site, as suggested by molecular modeling studies using the docking technique. Also, tetrahedral intermediate mimicry is achieved by the enol form, but not by the corresponding keto form. A pseudo-Hill coefficient of 0.25 has been determined for the Boc-phenylalanine enol, and by multiplying the concentration of inhibitor by its pseudo-Hill coefficient, a Ki of 0.82 microM has been determined for the inhibitory E-enol form. Boc-tyrosine enol generated a Ki of 0.39 microM. The requirements at the P1' position to achieve good inhibition are also described.

Amino Acid Sequence↗

PubMed2004

Study on affinity profile toward native human and bovine adenosine receptors of a series of 1,8-naphthyridine derivatives.

A new series of 1,8-naphthyridine derivatives (29-44 and 46-52) bearing various substituents in different positions on the heterocyclic nucleus were synthesized in order to analyze the effects produced on the affinity toward the bovine adenosine receptors. These derivatives represent an extension of our previous work on this class of compounds with high affinity toward A(1) adenosine receptors.(19) The results of radioligand binding assays indicate that a large number of the 1,8-naphthyridine derivatives proved to be A(1) selective, with a high affinity toward bovine adenosine receptors in the low nanomolar range, and one (29) in the subnanomolar range. Furthermore, the new series of 1,8-naphthyridine derivatives (29-44 and 46-52), together with the analogous derivatives 1-28 previously studied,(19) were tested to evaluate their affinity toward human cortical A(1) receptors and human striatal A(2A) receptors. The results indicate that all the 1,8-naphthyridine compounds generally possess a higher affinity toward the bovine A(1) receptor compared with the human A(1) receptor. As regards the affinity toward the A(2A) bovine receptor, only a few compounds possess a moderate affinity, which for some compounds remained approximately the same toward the A(2A) human receptor. A molecular modeling study of the docking of the 1,8-naphthyridine compounds with both the bovine and the human A(1) adenosine receptors was carried out with the aim of explaining the marked decrease in the affinity toward human A(1) adenosine receptors in comparison with bovine A(1) adenosine receptors. This study indicated that the structural differences, albeit small, of the active sites of the two receptors make differences in the dimensions of the site and this influenced the ability of the title compounds to interact with the two A(1) receptors.

Animals↗

PubMed1999

Conformationally constrained butyrophenones with mixed dopaminergic (D(2)) and serotoninergic (5-HT(2A), 5-HT(2C)) affinities: synthesis, pharmacology, 3D-QSAR, and molecular modeling of (aminoalkyl)benzo- and -thienocycloalkanones as putative atypical antipsychotics.

A series of novel conformationally restricted butyrophenones (2-(aminoethyl)- and 3-(aminomethyl)thieno- or benzocycloalkanones bearing (6-fluorobenzisoxazolyl)piperidine, (p-fluorobenzoyl)piperidine, (o-methoxyphenyl)piperazine, or linear butyrophenone fragments) were prepared and evaluated as atypical antipsychotic agents by in vitro assays of affinity for dopamine receptors (D(1), D(2)) and serotonin receptors (5-HT(2A), 5-HT(2C)) and by in vivo assays of antipsychotic potential and the risk of inducing extrapyramidal side effects. Potency and selectivity depended mainly on the amine fragment connected to the cycloalkanone structure. As a group, compounds with a benzisoxazolyl fragment had the highest 5-HT(2A) activities, followed by the benzoylpiperidine derivatives; in general, alpha-substituted cycloalkanone derivatives were more active than the corresponding beta-substituted congeners. CoMFA (comparative molecular field analysis) and docking studies showed electrostatic, steric, and lipophilic determinants of 5-HT(2A) and D(2) affinities and 5-HT(2A)/D(2) selectivity. The in vitro and in vivo pharmacological profiles of N-[(4-oxo-4H-5, 6-dihydrocyclopenta[b]thiophene-5-yl)ethyl]-4-(6-fluorobenzisox azol-3 -yl)piperidine (23b, QF 0510B), N-[(4-oxo-4,5,6, 7-tetrahydrobenzo[b]thiophene-5-yl)ethyl]-4-(6-fluorobenzisoxazol- 3-y l)piperidine (24b, QF 0610B), and N-[(7-oxo-4,5,6, 7-tetrahydrobenzo[b]thiophene-6-yl)ethyl]-4-(6-fluorobenzisoxazol- 3-y l)piperidine (29b, QF 0902B) suggest that they may be effective antipsychotic drugs with low propensity to induce extrapyramidal side effects.

Animals↗

PubMed1999

Deguanylation of human immunodeficiency virus (HIV-1) RNA by recombinant pokeweed antiviral protein.

Modeling studies, combined with the molecular docking of the trinucleotide GGG into the active site of the deadenylating RNA N-glycosidase pokeweed antiviral protein (PAP), indicated that a guanine base can fit into the active site pocket of PAP without disturbing its unique geometry and is sandwiched between residues Tyr(72) and Tyr(123) very much like an adenine base. The guanine base can form two specific hydrogen bonds with the active site residues Ser(121) and Val(73) and the attached negatively charged phosphate groups can entertain stabilizing electrostatic interactions with two clusters of positively charged patches on the PAP surface formed by Lys(210) and Arg(179) from one side and Arg(122) and Arg(135) from the other side of the active site. These observations prompted the hypothesis that the RNA depurinating activity of PAP may not be restricted to adenine residues and PAP should be capable of deguanylating ribosomal and viral RNA as well. This hypothesis was experimentally confirmed by direct demonstration that guanine base is released from both ribosomal and HIV-1 RNA after treatment with purified recombinant PAP using quantitative high performance liquid chromatography. Recombinant PAP released adenine and guanine residues at a 1:1 ratio from HIV-1 RNA and at an approximately 3:1 (adenine:guanine) ratio from Escherichia coli ribosomal RNA. At a concentration of 5 microM, recombinant PAP released 263 +/- 10 pmol of adenine and 100 +/- 11 pmol of guanine from 1 microgram of E. coli ribosomal RNA (16S + 23S) within 4 h of treatment. By comparison, 138 +/- 12 pmol of adenine and 143 +/- 10 pmol of guanine were released from 1 microgram of HIV-1 RNA under identical treatment conditions (5 microM recombinant PAP, 4 h treatment). The deguanylation of the ribosomal and viral RNA targets by recombinant PAP was concentration-dependent and is abolished by alanine substitutions of the catalytic active site residues Tyr(72) and Tyr(123). To our knowledge, these findings provide the first evidence that PAP can deguanylate both ribosomal and viral RNA.

Adenine Nucleotides↗

PubMed2003

Computational modelling of the open-state Kv 1.5 ion channel block by bupivacaine.

Binding of R(+)-bupivacaine to open-state homology models of the mammalian K(v)1.5 membrane ion channel is studied using automated docking and molecular dynamics (MD) methods. Homology models of K(v)1.5 are built using the 3D structures of the KcsA and MthK channels as a template. The packing of transmembrane (TM) alpha-helices in the KcsA structure corresponds to a closed channel state. Opening of the channel may be reached by a conformational transition yielding a bent structure of the internal S6 helices. Our first model of the K(v) open state involves a PVP-type of bending hinge in the internal helices, while the second model corresponds to a Gly-type of bending hinge as found in the MthK channel. Ligand binding to these models is probed using the common local anaesthetic bupivacaine, where blocker binding from the intracellular side of the channel is considered. Conformational properties and partial atomic charges of bupivacaine are determined from quantum mechanical HF/6-31G* calculations with inclusion of solvent effects. The automated docking and MD calculations for the PVP-bend model predict that bupivacaine could bind either in the central cavity or in the PVP region of the channel pore. Linear interaction energy (LIE) estimates of the binding free energies for bupivacaine predict strongest binding to the PVP region. Surprisingly, no binding is predicted for the Gly-bend model. These results are discussed in light of electrophysiological data which show that the K(v)1.5 channel is unable to close when bupivacaine is bound.

Amino Acid Sequence↗

PubMed2003

Design, synthesis and biological evaluation of ambenonium derivatives as AChE inhibitors.

Ambenonium (1), an old AChE inhibitor, is endowed with an outstanding affinity and a peculiar mechanism of action that, taken together, make it a very promising pharmacological tool for the treatment of Alzheimer's disease (AD). Unfortunately, the bisquaternary structure of 1 prevents its passage through the blood brain barrier. In a search of centrally active ambenonium derivatives, we planned to synthesize tertiary amines of 1, such as 2 and 3. In addition, to add new insights into the binding mechanism of the inhibitor, we designed constrained analogues of ambenonium by incorporating the diamine functions into cyclic moieties (4-12). The biological evaluation of the new compounds has been assessed in vitro against human AChE and BChE. All tertiary amine derivatives resulted more than 1000-fold less potent than 1 and, unlike prototype, did not show any selectivity between the two enzymes. This result, because of recent findings concerning the role of BChE in AD, makes our compounds, endowed with a well-balanced profile of AChE/BChE inhibition, valuable candidates for further development. To better clarify the interactions that account for the high affinity of 1, docking simulations and molecular dynamics studies on the AChE-1 complex were also carried out.

Acetylcholinesterase↗

PubMed2006

Metabolism of bergamottin by cytochromes P450 2B6 and 3A5.

Cytochromes P450 (P450) 2B6 and 3A5 are inactivated by bergamottin (BG). P450 2B6 metabolized BG primarily to M3 and M4 and one minor metabolite (M1). The metabolites were analyzed, and the data indicated that M1 was bergaptol, M3 was 5'-OH-BG, and M4 was a mixture of 6'- and 7'-OH-BG. Because 6'- and 7'-OH-BG were the primary metabolites, it suggested that P450 2B6 preferentially oxidized the geranyloxy chain of BG. Metabolism of BG by P450 3A5 resulted in three major metabolites: [bergaptol, M3 (5'-OH-BG), and M5 (2'-OH-BG)], and two minor metabolites [M2 (6',7'-dihydroxy-BG) and M4 (6'- and 7'-OH-BG)]. Because bergaptol was the most abundant metabolite formed, it suggested that P450 3A5 metabolized BG mainly by cleaving the geranyl-oxy chain. Molecular modeling studies confirmed that docking of BG in the P450 2B6 active site favors oxidation in the terminal region of the geranyl-oxy chain, whereas positioning the 2'-carbon of BG nearest the heme iron is preferred by P450 3A5. Glutathione (GSH)-BG conjugates were formed by both P450. Each enzyme predominantly formed conjugates with m/z values of 662. Tandem mass spectrometry analysis of the GSH conjugates indicated that the oxidation forming a reactive intermediate occurred on the furan moiety of BG, presumably through the initial formation of an epoxide at the furan double bond. The data indicate that oxidation of the geranyl-oxy chain resulted in the formation of stable metabolites of BG, whereas oxidation of the furan ring produced reactive intermediates that may be responsible for binding to and inactivating P450 2B6 and 3A4.

Aryl Hydrocarbon Hydroxylases↗

PubMed2006

Discovery of potential sorbitol dehydrogenase inhibitors from virtual screening.

Sorbitol dehydrogenase (SDH) is the second enzyme in the polyol pathway of glucose metabolism and is a possible target for the treatment of the complications of diabetes. In this study the molecular modelling program DOCK was used to analyse 249,071 compounds from the National Cancer Institute Database and predict those with high affinity for SDH. From a total of 21 tested the 7 compounds including flavin adenine dinucleotide disodium hydrate, (+)-Amethopterin, 3-hydroxy-2-napthoic(2-hydroxybenzylidene) hydrazide, folic acid, N-2,4-dinitrophenyl-L-cysteic acid, Vanillin azine and 1H-indole-2,3-dione,5-bromo-6-nitro-1-(2,3,4-tri-O-acetyl-alpha-L-arabinopyranosyl)-(9Cl), were shown to inhibit SDH and displayed IC50 values of 0.192 microM, 1.1 microM, 1.2 microM, 4.5 microM, 5.3 microM, 7 microM and 28 microM, respectively. These compounds may aid the design of pharmaceutical agents for the treatment of diabetes complications.

Database Management Systems↗

PubMedDate not supplied

Computer aids in drug design--highlights.

An outline of some procedures used for computer aided drug design has been given. The emphasis is on lead generation, both primary and secondary, and on current lead optimization processes which rely on an appreciation of ligand-receptor interactions. The areas covered include techniques for structure based drug design and lead optimization, ligand-receptor interaction studies using molecular surfaces and docking, extraction of pharmacophore pattern to obtain new leads, expert system application etc.

Computer-Aided Design↗

PubMed2005

Structure-based discovery of human L-xylulose reductase inhibitors from database screening and molecular docking.

Human L-xylulose reductase (XR) is an enzyme of the glucuronic acid/uronate cycle of glucose metabolism and is a possible target for treatment of the long-term complications of diabetes. In this study we utilised the molecular modelling program DOCK to analyse the 249,071 compounds of the National Cancer Institute Database and retrieved those compounds with high predicted affinity for XR. Several carboxylic acid-based compounds were tested and shown to inhibit XR. These included nicotinic acid (IC50=100 microM), benzoic acid (IC50=29 microM) and their derivatives. These results extend and improve upon the activities of known, commercially available inhibitors of XR such as the aliphatic fatty acid n-butyric acid (IC50=64 microM). To optimise the interaction between the inhibitor and the holoenzyme, the program GRID was used to design de novo compounds based on the inhibitor benzoic acid. The inclusion of a hydroxy-phenyl group and a phosphate to the benzoic acid molecule increased the net binding energy by 1.3- and 2.4-fold, respectively. The resultant compounds may produce inhibitors with improved specificity for XR.

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↗

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↗