Mechanism of acetylcholinesterase inhibition by fasciculin: a 5-ns molecular dynamics simulation

AChE activity self-breathing control mechanisms regulated by H2Sn and GSH: Persulfidation and glutathionylation on sulfhydryl after disulfide bonds cleavage

Hydrogen sulfide (H2S), or dihydrogen sulfane (H2Sn), acts as a signal molecule through the beneficial mechanism of persulfidation, known as the post-translational transformation of cysteine residues to persulfides. We previously reported that Glutathione (GSH) could regulate enzyme activity through S-desulfurization or glutathionylation of residues to generate protein-SG or protein-SSG, releasing H2S. However, little is known about the mechanisms by which H2Sn and GSH affect the disulfide bonds. In this study, we provide direct evidences that H2Sn and GSH modify the sulfhydryl group on Cys272, which forms disulfide bonds in acetylcholinesterase (AChE), to generate Cys-SSH and Cys-SSG, respectively. Glutathionylation of disulfide is a two-step reaction based on nucleophilic substitution, in which the first CS bond is broken, then the SS bond is broken to release H2S. H2Sn and GSH controlled self-breathing motion in enzyme catalysis by disconnecting specific disulfide bonds and modifying cysteine residues, thereby regulating AChE activity. Here, we elucidated H2Sn and GSH mechanisms on disulfide in the AChE system and proposed a self-breathing control theory induced by H2Sn and GSH. These theoretical findings shed light on the biological functions of H2Sn and GSH on sulfhydryl and disulfide bonds and enrich the theory of enzyme activity regulation.

Disclosing gate-opening/closing events inside a flexible metal-organic framework loaded with CO2 by reactive and essential dynamics

We have combined reactive molecular dynamics simulations with principal component analysis to provide a clearer view of the interactions and motion of the CO2 molecules inside a metal-organic framework and the movements of the MOF components that regulate storage, adsorption, and diffusion of the guest species. The tens-of-nanometer size of the simulated model, the capability of the reactive force field tuned to reproduce the inorganic-organic material confidently, and the unconventional use of essential dynamics have effectively disclosed the gate-opening/closing phenomenon, possible coordinations of CO2 at the metal centers, all the diffusion steps inside the MOF channels, the primary motions of the linkers, and the effects of their concerted rearrangements on local CO2 relocations.

Computational Studies on Acetylcholinesterase Inhibitors: From Biochemistry to Chemistry

Acetylcholinesterase inhibitors are the most promising therapeutics for Alzheimer's disease treatment as these prevent the loss of acetylcholine and slows the progression of the disease. The drugs approved for the management of Alzheimer's disease by the FDA are acetylcholinesterase inhibitors but are associated with side effects. Consistent and stringent efforts by the researchers with the help of computational methods opened new ways of developing novel molecules with good acetylcholinesterase inhibitory activity. In this manuscript, we reviewed the studies that identified the essential structural features of acetylcholinesterase inhibitors at the molecular level as well as the techniques like molecular docking, molecular dynamics, quantitative structure-activity relationship, virtual screening, and pharmacophore modelling that were used in designing these inhibitors.

A Comprehensive Review of Cholinesterase Modeling and Simulation

The present article reviews published efforts to study acetylcholinesterase and butyrylcholinesterase structure and function using computer-based modeling and simulation techniques. Structures and models of both enzymes from various organisms, including rays, mice, and humans, are discussed to highlight key structural similarities in the active site gorges of the two enzymes, such as flexibility, binding site location, and function, as well as differences, such as gorge volume and binding site residue composition. Catalytic studies are also described, with an emphasis on the mechanism of acetylcholine hydrolysis by each enzyme and novel mutants that increase catalytic efficiency. The inhibitory activities of myriad compounds have been computationally assessed, primarily through Monte Carlo-based docking calculations and molecular dynamics simulations. Pharmaceutical compounds examined herein include FDA-approved therapeutics and their derivatives, as well as several other prescription drug derivatives. Cholinesterase interactions with both narcotics and organophosphate compounds are discussed, with the latter focusing primarily on molecular recognition studies of potential therapeutic value and on improving our understanding of the reactivation of cholinesterases that are bound to toxins. This review also explores the inhibitory properties of several other organic and biological moieties, as well as advancements in virtual screening methodologies with respect to these enzymes.

1-(3-Tert-Butylphenyl)-2,2,2-Trifluoroethanone as a Potent Transition-State Analogue Slow-Binding Inhibitor of Human Acetylcholinesterase: Kinetic, MD and QM/MM Studies

Kinetic studies and molecular modeling of human acetylcholinesterase (AChE) inhibition by a fluorinated acetophenone derivative, 1-(3-tert-butylphenyl)-2,2,2-trifluoroethanone (TFK), were performed. Fast reversible inhibition of AChE by TFK is of competitive type with K_i = 5.15 nM. However, steady state of inhibition is reached slowly. Kinetic analysis showed that TFK is a slow-binding inhibitor (SBI) of type B with K_i* = 0.53 nM. Reversible binding of TFK provides a long residence time, τ = 20 min, on AChE. After binding, TFK acylates the active serine, forming an hemiketal. Then, disruption of hemiketal (deacylation) is slow. AChE recovers full activity in approximately 40 min. Molecular docking and MD simulations depicted the different steps. It was shown that TFK binds first to the peripheral anionic site. Then, subsequent slow induced-fit step enlarged the gorge, allowing tight adjustment into the catalytic active site. Modeling of interactions between TFK and AChE active site by QM/MM showed that the “isomerization” step of enzyme-inhibitor complex leads to a complex similar to substrate tetrahedral intermediate, a so-called “transition state analog”, followed by a labile covalent intermediate. SBIs of AChE show prolonged pharmacological efficacy. Thus, this fluoroalkylketone intended for neuroimaging, could be of interest in palliative therapy of Alzheimer’s disease and protection of central AChE against organophosphorus compounds.

Modeling Epac1 interactions with the allosteric inhibitor AM-001 by co-solvent molecular dynamics

The exchange proteins activated by cAMP (EPAC) are implicated in a large variety of physiological processes and they are considered as promising targets for a wide range of therapeutic applications. Several recent reports provided evidence for the therapeutic effectiveness of the inhibiting EPAC1 activity cardiac diseases. In that context, we recently characterized a selective EPAC1 antagonist named AM-001. This compound was featured by a non-competitive mechanism of action but the localization of its allosteric site to EPAC1 structure has yet to be investigated. Therefore, we performed cosolvent molecular dynamics with the aim to identify a suitable allosteric binding site. Then, the docking and molecular dynamics were used to determine the binding of the AM-001 to the regions highlighted by cosolvent molecular dynamics for EPAC1. These analyses led us to the identification of a suitable allosteric AM-001 binding pocket at EPAC1. As a model validation, we also evaluated the binding poses of the available AM-001 analogues, with a different biological potency. Finally, the complex EPAC1 with AM-001 bound at the putative allosteric site was further refined by molecular dynamics. The principal component analysis led us to identify the protein motion that resulted in an inactive like conformation upon the allosteric inhibitor binding.

Molecular modeling, simulation and principal component analysis of binding of resveratrol and its analogues with DNA

Structure-based drug designing has become a significant subject of research, and several clinically promising DNA binding compounds were evolved using this technique. The interaction of an octamer DNA sequence d(CCAATTGG)₂ with a natural stilbene, resveratrol and its analogues have been studied using molecular docking method. Out of the ten compounds studied, seven compounds were found to bind to the minor groove of AATT segment of the sequence. Pterostilbene, a natural analogue of resveratrol, showed the lowest binding energy. Rhaponticin, a natural analogue of resveratrol and digalloylresveratrol, a synthetic ester of resveratrol bind to the major groove of the AATT segment while dihydroresveratrol binds to the minor groove of GC terminal base pair. ADMET (Absorption, distribution, metabolism, excretion and toxicity) study showed that all compounds obey Lipinski rule and are accepted as orally active drugs based on different physicochemical descriptors. Molecular dynamics simulations were performed for the complex with lowest binding energy and trajectory analysis were performed. Principal component analysis has been performed to underline the prominent motions in alone DNA and when it is bound to pterostilbene. AbbreviationsADMETAbsorption, distribution, metabolism, excretion and toxicityDIGDigalloyl resveratrolDNADeoxyribonucleic acidELElectrostatic energyENPOLARNonpolar solvation energyESURFSurface areaGBGeneralized BornHBAHydrogen bond acceptorsHBDHydrogen bond donorsLGALamarckian genetic algorithmMDMolecular dynamicsPBPoisson-BoltzmannPCAPrincipal component analysisPTPterostilbeneRMSDRoot mean square deviationSASimulated annealingTLX3T-cell leukemia homeobox 3VDWvan der WaalsCommunicated by Ramaswamy H. Sarma.

Dynamic Characterization of the Human Heme Nitric Oxide/Oxygen (HNOX) Domain under the Influence of Diatomic Gaseous Ligands

Soluble guanylate cyclase (sGC) regulates numerous physiological processes. The β subunit Heme Nitric Oxide/Oxygen (HNOX) domain makes this protein sensitive to small gaseous ligands. The structural basis of the activation mechanism of sGC under the influence of ligands (NO, O₂, CO) is poorly understood. We examine the effect of different ligands on the human sGC HNOX domain. HNOX systems with gaseous ligands were generated and explored using Molecular Dynamics (MD). The distance between heme Fe2+ and histidine in the NO-ligated HNOX (NO-HNOX) system is larger compared to the O₂, CO systems. NO-HNOX rapidly adopts the conformation of the five-group metal coordination system. Loops α, β, γ and helix-f exhibit increased mobility and different hydrogen bond networks in NO-HNOX compared to the other systems. The removal of His from the Fe coordination sphere in NO-HNOX is assisted by interaction of the imidazole ring with the surrounding residues which in turn leads to the release of signaling helix-f and activation of the sGC enzyme. Insights into the conformational dynamics of a human sGC HNOX domain, especially for regions which are functionally critical for signal transduction, are valuable in the understanding of cardiovascular diseases.

Identification of new allosteric sites and modulators of AChE through computational and experimental tools

Allosteric sites on proteins are targeted for designing more selective inhibitors of enzyme activity and to discover new functions. Acetylcholinesterase (AChE), which is most widely known for the hydrolysis of the neurotransmitter acetylcholine, has a peripheral allosteric subsite responsible for amyloidosis in Alzheimer's disease through interaction with amyloid β-peptide. However, AChE plays other non-hydrolytic functions. Here, we identify and characterise using computational tools two new allosteric sites in AChE, which have allowed us to identify allosteric inhibitors by virtual screening guided by structure-based and fragment hotspot strategies. The identified compounds were also screened for in vitro inhibition of AChE and three were observed to be active. Further experimental (kinetic) and computational (molecular dynamics) studies have been performed to verify the allosteric activity. These new compounds may be valuable pharmacological tools in the study of non-cholinergic functions of AChE.

Computational Studies on Acetylcholinesterases

Functions of biomolecules, in particular enzymes, are usually modulated by structural fluctuations. This is especially the case in a gated diffusion-controlled reaction catalyzed by an enzyme such as acetylcholinesterase. The catalytic triad of acetylcholinesterase is located at the bottom of a long and narrow gorge, but it catalyzes the extremely rapid hydrolysis of the neurotransmitter, acetylcholine, with a reaction rate close to the diffusion-controlled limit. Computational modeling and simulation have produced considerable advances in exploring the dynamical and conformational properties of biomolecules, not only aiding in interpreting the experimental data, but also providing insights into the internal motions of the biomolecule at the atomic level. Given the remarkably high catalytic efficiency and the importance of acetylcholinesterase in drug development, great efforts have been made to understand the dynamics associated with its functions by use of various computational methods. Here, we present a comprehensive overview of recent computational studies on acetylcholinesterase, expanding our views of the enzyme from a microstate of a single structure to conformational ensembles, strengthening our understanding of the integration of structure, dynamics and function associated with the enzyme, and promoting the structure-based and/or mechanism-based design of new inhibitors for it.

The plastoquinol-plastoquinone exchange mechanism in photosystem II: insight from molecular dynamics simulations

In the photosystem II (PSII) of oxygenic photosynthetic organisms, the reaction center (RC) core mediates the light-induced electron transfer leading to water splitting and production of reduced plastoquinone molecules. The reduction of plastoquinone to plastoquinol lowers PSII affinity for the latter and leads to its release. However, little is known about the role of protein dynamics in this process. Here, molecular dynamics simulations of the complete PSII complex embedded in a lipid bilayer have been used to investigate the plastoquinol release mechanism. A distinct dynamic behavior of PSII in the presence of plastoquinol is observed which, coupled to changes in charge distribution and electrostatic interactions, causes disruption of the interactions seen in the PSII-plastoquinone complex and leads to the "squeezing out" of plastoquinol from the binding pocket. Displacement of plastoquinol closes the second water channel, recently described in a 2.9 Å resolution PSII structure (Guskov et al. in Nat Struct Mol Biol 16:334-342, 2009), allowing to rule out the proposed "alternating" mechanism of plastoquinol-plastoquinone exchange, while giving support to the "single-channel" one. The performed simulations indicated a pivotal role of D₁-Ser264 in modulating the dynamics of the plastoquinone binding pocket and plastoquinol-plastoquinone exchange via its interaction with D₁-His252 residue. The effects of the disruption of this hydrogen bond network on the PSII redox reactions were experimentally assessed in the D₁ site-directed mutant Ser264Lys.

Molecular investigations of protriptyline as a multi-target directed ligand in Alzheimer's disease

Alzheimer's disease (AD) is a complex neurodegenerative disorder involving multiple cellular and molecular processes. The discovery of drug molecules capable of targeting multiple factors involved in AD pathogenesis would greatly facilitate in improving therapeutic strategies. The repositioning of existing non-toxic drugs could dramatically reduce the time and costs involved in developmental and clinical trial stages. In this study, preliminary screening of 140 FDA approved nervous system drugs by docking suggested the viability of the tricyclic group of antidepressants against three major AD targets, viz. Acetylcholinesterase (AChE), β-secretase (BACE-1), and amyloid β (Aβ) aggregation, with one member, protriptyline, showing highest inhibitory activity. Detailed biophysical assays, together with isothermal calorimetry, fluorescence quenching experiments, kinetic studies and atomic force microscopy established the strong inhibitory activity of protriptyline against all three major targets. The molecular basis of inhibition was supported with comprehensive molecular dynamics simulations. Further, the drug inhibited glycation induced amyloid aggregation, another important causal factor in AD progression. This study has led to the discovery of protriptyline as a potent multi target directed ligand and established its viability as a promising candidate for AD treatment.

An investigation of G-quadruplex structural polymorphism in the human telomere using a combined approach of hydrodynamic bead modeling and molecular dynamics simulation

Guanine-rich oligonucleotides can adopt noncanonical tertiary structures known as G-quadruplexes, which can exist in different forms depending on experimental conditions. High-resolution structural methods, such as X-ray crystallography and NMR spectroscopy, have been of limited usefulness in resolving the inherent structural polymorphism associated with G-quadruplex formation. The lack of, or the ambiguous nature of, currently available high-resolution structural data, in turn, has severely hindered investigations into the nature of these structures and their interactions with small-molecule inhibitors. We have used molecular dynamics in conjunction with hydrodynamic bead modeling to study the structures of the human telomeric G-quadruplex-forming sequences at the atomic level. We demonstrated that molecular dynamics can reproduce experimental hydrodynamic measurements and thus can be a powerful tool in the structural study of existing G-quadruplex sequences or in the prediction of new G-quadruplex structures.

Nature: a substantial source of auspicious substances with acetylcholinesterase inhibitory action

Acetylcholinesterase (AChE) (EC 3.1.1.7) is an important enzyme that breaks down of acetylcholine in synaptic cleft in neuronal junctions. Inhibition of AChE is associated with treatment of several diseases such as Alzheimer's disease (AD), myasthenia gravis, and glaucoma as well as the mechanisms of insecticide and anthelmintic drugs. Several AChE inhibitors are available in clinical use currently for the treatment of AD; however, none of them has ability, yet, to seize progress of the disease. Consequently, an extensive research has been going on finding new AChE inhibitors. In this sense, natural inhibitors have gained great attention due to their encouraging effects toward AChE. In this review, promising candidate molecules with marked AChE inhibition from both plant and animal sources will be underlined.

Molecular dynamics of HIV1-integrase in complex with 93del - a structural perspective on the mechanism of inhibition

HIV1 integrase is an important target for the antiviral therapy. Guanine-rich quadruplex, such as 93del, have been shown to be potent inhibitors of this enzyme and thus representing a new class of antiviral agents. Although X-ray and NMR structures of HIV1 integrase and 93del have been reported, there is no structural information of the complex and the mechanism of inhibition still remains unexplored. A number of computational methods including automated protein-DNA docking and molecular dynamics simulation in explicit solvent were used to model the binding of 93del to HIV1 integrase. Analysis of the dynamic behaviour of the complex using principal components analysis and elastic network modelling techniques allow us to understand how the binding of 93del aptamer and its interactions with key residues affect the intrinsic motions of the catalytic loops by stabilising them in catalytically inactive conformations. Such insights into the structural mechanism of inhibition can aid in improving the design of anti-HIV aptamers.

Insight into TPMT(∗)23 mutation mis-folding using molecular dynamics simulation and protein structure analysis

Thiopurine S-methyltransferase (TPMT) is an important enzyme that metabolizes thiopurine drugs. This enzyme exhibits a large number of interindividual polymorphism. TPMT(∗)23 polymorphism has been reported in a few cases in the world in co-dominance with TPMT(∗)3A. The phenotype has been reported to affect enzyme activity in vivo and in vitro. Its underlying structural basis is not clarified yet. In our study, the wild type (WT) protein structure was analyzed and the amino acids bordering water channels in thiopurine sites were identified. Molecular dynamics of both the WT and TPMT(∗)23 mutation was carried out. In addition, the effects of this mutation, especially on the thiopurine site which is closed with a pincer like mechanism, were investigated. We focused on explaining how a locally occurred A167G substitution propagated through hydrogen bonds alteration to induce structural modification which affects both thiopurine and S-adenosylmethionine receptors. Finally, a genetic prediction of mutation functional consequences has been conducted confirming altered activity. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:20.

Biochemistry of envenomation

Venoms and toxins are of significant interest due to their ability to cause a wide range of pathophysiological conditions that can potentially result in death. Despite their wide distribution among plants and animals, the biochemical pathways associated with these pathogenic agents remain largely unexplored. Impoverished and underdeveloped regions appear especially susceptible to increased incidence and severity due to poor socioeconomic conditions and lack of appropriate medical treatment infrastructure. To facilitate better management and treatment of envenomation victims, it is essential that the biochemical mechanisms of their action be elucidated. This review aims to characterize downstream envenomation mechanisms by addressing the major neuro-, cardio-, and hemotoxins as well as ion-channel toxins. Because of their use in folk and traditional medicine, the biochemistry behind venom therapy and possible implications on conventional medicine will also be addressed.

Molecular dynamics simulations of the coenzyme induced conformational changes of Mycobacterium tuberculosis L-alanine dehydrogenase

Mycobacterium tuberculosis L-alanine dehydrogenase (L-MtAlaDH) catalyzes the NADH-dependent reversible oxidative deamination of L-alanine to pyruvate and ammonia. L-MtAlaDH has been proposed to be a potential target in the treatment of tuberculosis. Based on the crystal structures of this enzyme, molecular dynamics simulations were performed to investigate the conformational changes of L-MtAlaDH induced by coenzyme NADH. The results show that the presence of NADH in the binding domain restricts the motions and conformational distributions of L-MtAlaDH. There are two loops (residues 94-99 and 238-251) playing important roles for the binding of NADH, while another loop (residues 267-293) is responsible for the binding of substrate. The opening/closing and twisting motions of two domains are closely related to the conformational changes of L-MtAlaDH induced by NADH.

Active site gating and substrate specificity of butyrylcholinesterase and acetylcholinesterase: insights from molecular dynamics simulations

Butyrylcholinesterase (BChE) and acetylcholinesterase (AChE) are highly homologous proteins with distinct substrate preferences. In this study we compared the active sites of monomers and tetramers of human BChE and human AChE after performing molecular dynamics (MD) simulations in water-solvated systems. By comparing the conformational dynamics of gating residues of AChE and BChE, we found that the gating mechanisms of the main door of AChE and BChE are responsible for their different substrate specificities. Our simulation of the tetramers of AChE and BChE indicates that both enzymes could have two dysfunctional active sites due to their restricted accessibility to substrates. The further study on catalytic mechanisms of multiple forms of AChE and BChE would benefit from our comparison of the active sites of the monomers and tetramers of both enzymes.

Long route or shortcut? A molecular dynamics study of traffic of thiocholine within the active-site gorge of acetylcholinesterase

The principal role of acetylcholinesterase is termination of nerve impulse transmission at cholinergic synapses, by rapid hydrolysis of the neurotransmitter acetylcholine to acetate and choline. Its active site is buried at the bottom of a deep and narrow gorge, at the rim of which is found a second anionic site, the peripheral anionic site. The fact that the active site is so deeply buried has raised cogent questions as to how rapid traffic of substrate and products occurs in such a confined environment. Various theoretical and experimental approaches have been used to solve this problem. Here, multiple conventional molecular dynamics simulations have been performed to investigate the clearance of the product, thiocholine, from the active-site gorge of acetylcholinesterase. Our results indicate that thiocholine is released from the peripheral anionic site via random pathways, while three exit routes appear to be favored for its release from the active site, namely, along the axis of the active-site gorge, and through putative back- and side-doors. The back-door pathway is that via which thiocholine exits most frequently. Our results are in good agreement with kinetic and kinetic-crystallography studies. We propose the use of multiple molecular dynamics simulations as a fast yet accurate complementary tool in structural studies of enzymatic trafficking.

Dynamic conformational changes in munc18 prevent syntaxin binding

The Sec1/munc18 protein family is essential for vesicle fusion in eukaryotic cells via binding to SNARE proteins. Protein kinase C modulates these interactions by phosphorylating munc18a thereby reducing its affinity to one of the central SNARE members, syntaxin-1a. The established hypothesis is that the reduced affinity of the phosphorylated munc18a to syntaxin-1a is a result of local electrostatic repulsion between the two proteins, which interferes with their compatibility. The current study challenges this paradigm and offers a novel mechanistic explanation by revealing a syntaxin-non-binding conformation of munc18a that is induced by the phosphomimetic mutations. In the present study, using molecular dynamics simulations, we explored the dynamics of the wild-type munc18a versus phosphomimetic mutant munc18a. We focused on the structural changes that occur in the cavity between domains 3a and 1, which serves as the main syntaxin-binding site. The results of the simulations suggest that the free wild-type munc18a exhibits a dynamic equilibrium between several conformations differing in the size of its cavity (the main syntaxin-binding site). The flexibility of the cavity's size might facilitate the binding or unbinding of syntaxin. In silico insertion of phosphomimetic mutations into the munc18a structure induces the formation of a conformation where the syntaxin-binding area is rigid and blocked as a result of interactions between residues located on both sides of the cavity. Therefore, we suggest that the reduced affinity of the phosphomimetic mutant/phosphorylated munc18a is a result of the closed-cavity conformation, which makes syntaxin binding energetically and sterically unfavorable. The current study demonstrates the potential of phosphoryalation, an essential biological process, to serve as a driving force for dramatic conformational changes of proteins modulating their affinity to target proteins.

Protein phosphorylation plays a significant regulatory role in multi-component systems engaged in signal transduction or coordination of cellular processes, by activating or deactivating proteins. The potential of phosphorylation to induce substantial conformational changes in proteins, thereby changing their affinity to target proteins, has already been shown but the dynamics of the process is not fully elucidated. In the present study, we investigated, by molecular dynamics simulations, the dynamic conformational changes in munc18a, a protein that is crucial for neurotransmitter release and interacts tightly with the SNARE syntaxin-1. We further investigated the conformational changes that occur in munc18a when it is phosphorylated, reducing its affinity to syntaxin-1a. The results of the simulations suggest that there is a conformational flexibility of the syntaxin-unbounded munc18a that allows changes in the shape of the syntaxin-1a binding cavity. In silico insertion of phosphomimetic mutations into munc18a led to a reduction in the flexibility and closure of the syntaxin-binding site. We suggest that the reduced affinity of phosphorylated munc18a to syntaxin-1a stems from the difficulty of syntaxin-1a to bind to the munc18a closed-cavity conformation, induced by the PKC phosphorylation of munc18a.

Influence of the acetylcholinesterase active site protonation on omega loop and active site dynamics

Existence of alternative entrances in acetylcholinesterase (AChE) could explain the contrast between the very high AChE catalytic efficiency and the narrow and long access path to the active site revealed by X-ray crystallography. Alternative entrances could facilitate diffusion of the reaction products or at least water and ions from the active site. Previous molecular dynamics simulations identified side door and back door as the most probable alternative entrances. The simulations of non-inhibited AChE suggested that the back door opening events occur only rarely (0.8% of the time in the 10ns trajectory). Here we present a molecular dynamics simulation of non-inhibited AChE, where the back door opening appears much more often (14% of the time in the 12ns trajectory) and where the side door opening was observed quite frequently (78% of trajectory time). We also present molecular dynamics, where the back door does not open at all, or where large conformational changes of the AChE omega loop occur together with alternative passage opening events. All these differences in AChE dynamical behavior are caused by different protonation states of two glutamate residues located on bottom of the active site gorge (Glu202 and G450 in Mus musculus AChE). Our results confirm the results of previous molecular dynamics simulations, expand the view and suggest the probable reasons for the overall conformational behavior of AChE omega loop.

Understanding the specificity of a docking interaction between JNK1 and the scaffolding protein JIP1

The up-regulation of JNK activity is associated with a number of disease states. The JNK-JIP1 interaction represents an attractive target for the inhibition of JNK-mediated signaling. In this study, molecular dynamics simulations have been performed on the apo-JNK1 and the JNK1•L-pepJIP1 and JNK1•D-pepJIP1 complexes to investigate the interaction between the JIP1 peptides and JNK1. Dynamic domain studies based on essential dynamics (ED) analysis of apo-JNK1 and the JNK1•L-pepJIP1 complex have been performed to analyze and compare details of conformational changes, hinge axes, and hinge bending regions in both structures. The activation loop, the αC helix, and the G loop are found to be highly flexible and to exhibit significant changes in dynamics upon L-pepJIP1 binding. The conformation of the activation loop for the apo state is similar to that of inactive apo-ERK2, while the activation loop in JNK1•L-pepJIP1 complex resembles that of the inactive ERK2 bound with pepHePTP. ED analysis shows that, after the binding of l-pepJIP1, the N- and C-terminal domains of JNK1 display both a closure and a twisting motion centered around the activation loop, which functions as a hinge. In contrast, no domain motion is detected for the apo state for which an open conformation is favored. The present study suggests that L-pepJIP1 regulates the interdomain motions of JNK1 and potentially the active site via an allosteric mechanism. The binding free energies of L-pepJIP1 and D-pepJIP1 to JNK1 are estimated using the molecular mechanics Poisson-Boltzmann and generalized-Born surface area (MM-PB/GBSA) methods. The contribution of each residue at the interaction interface to the binding affinity of L-pepJIP1 with JNK1 has been analyzed by means of computational alanine-scanning mutagenesis and free energy decomposition. Several critical interactions for binding (e.g., Arg156/L-pepJIP1 and Glu329/JNK1) have been identified. The binding free energy calculation indicates that the electrostatic interaction contributes critically to specificity, rather than to binding affinity between the peptide and JNK1. Notably, the binding free energy calculations predict that D-pepJIP1 binding to JNK1 is significantly weaker than the L form, contradicting the previous suggestion that D-pepJIP1 acts as an inhibitor toward JNK1. We have performed experiments using purified JNK1 to confirm that, indeed, D-pepJIP1 does not inhibit the ability of JNK1 to phosphorylate c-Jun in vitro.

Molecular dynamics simulations of the rotary motor F(0) under external electric fields across the membrane

The membrane-bound component F(0), which is a major component of the F(0)F(1)-ATP synthase, works as a rotary motor and plays a central role in driving the F(1) component to transform chemiosmotic energy into ATP synthesis. We conducted molecular dynamics simulations of b(2)-free F(0) in a 1-palmitoyl-2-oleoyl-phosphatidylcholine lipid bilayer for tens of nanoseconds with two different protonation states of the cAsp-61 residue at the interface of the a-c complex in the absence of electric fields and under electric fields of +/-0.03 V/nm across the membrane. To our surprise, we observed that the upper half of the N-terminal helix of the c(1) subunit rotated about its axis clockwise by 30 degrees . An energetic analysis revealed that the electrostatic repulsion between this N-terminal helix and subunit c(12) was a major contributor to the observed rotation. A correlation map analysis indicated that the correlated motions of residues in the interface of the a-c complex were significantly reduced by external electric fields. The deuterium order parameter (S(CD)) profile calculated by averaging all the lipids in the F(0)-bound bilayer was not very different from that of the pure bilayer system, in agreement with recent (2)H solid-state NMR experiments. However, by delineating the lipid properties according to their vicinity to F(0), we found that the S(CD) profiles of different lipid shells were prominently different. Lipids close to F(0) formed a more ordered structure. Similarly, the lateral diffusion of lipids on the membrane surface also followed a shell-dependent behavior. The lipids in the proximity of F(0) exhibited very significantly reduced diffusional motion. The numerical value of S(CD) was anticorrelated with that of the diffusion coefficient, i.e., the more ordered lipid structures led to slower lipid diffusion. Our findings will help elucidate the dynamics of F(0) depending on the protonation state and electric field, and may also shed some light on the interactions between the motor F(0) and its surrounding lipids under physiological conditions, which could help to rationalize its extraordinary energy conversion efficiency.

Differential binding of bispyridinium oxime drugs with acetylcholinesterase

AIM

To performe a time-dependent topographical delineation of protein-drug interactions to gain molecular insight into the supremacy of Ortho-7 over HI-6 in reactivating tabun-conjugated mouse acetylcholinesterase (mAChE).

METHODS

We conducted all-atom steered molecular dynamics simulations of the two protein-drug complexes. Through a host of protein-drug interaction parameters (rupture force profiles, hydrogen bonds, water bridges, hydrophobic interactions), geometrical, and orientation ordering of the drugs, we monitored the enzyme's response during the release of the drugs from its active-site.

RESULTS

The results show the preferential binding of the drugs with the enzyme. The pyridinium ring of HI-6 shows excellent complementary binding with the peripheral anionic site, whereas one of two identical pyridinium rings of Ortho-7 has excellent binding compatibility in the enzyme active-site where it can orchestrate the reactivation process. We found that the active pyridinium ring of HI-6 undergoes a complete turn along the active site axis, directed away from the active-site region during the course of the simulation.

CONCLUSION

Due to excellent cooperative binding of Ortho-7, as rendered by several cation-pi interactions with the active-site gorge of the enzyme, Ortho-7 may be a more efficient reactivator than HI-6. Our work supports the growing body of evidence that the efficacy of the drugs is due to the differential bindings of the oximes with AChE and can aid to the rational design of oxime drugs.

Molecular dynamics and principal components analysis of human telomeric quadruplex multimers

Guanine-rich DNA repeat sequences located at the terminal ends of chromosomal DNA can fold in a sequence-dependent manner into G-quadruplex structures, notably the terminal 150-200 nucleotides at the 3' end, which occur as a single-stranded DNA overhang. The crystal structures of quadruplexes with two and four human telomeric repeats show an all-parallel-stranded topology that is readily capable of forming extended stacks of such quadruplex structures, with external TTA loops positioned to potentially interact with other macromolecules. This study reports on possible arrangements for these quadruplex dimers and tetramers, which can be formed from 8 or 16 telomeric DNA repeats, and on a methodology for modeling their interactions with small molecules. A series of computational methods including molecular dynamics, free energy calculations, and principal components analysis have been used to characterize the properties of these higher-order G-quadruplex dimers and tetramers with parallel-stranded topology. The results confirm the stability of the central G-tetrads, the individual quadruplexes, and the resulting multimers. Principal components analysis has been carried out to highlight the dominant motions in these G-quadruplex dimer and multimer structures. The TTA loop is the most flexible part of the model and the overall multimer quadruplex becoming more stable with the addition of further G-tetrads. The addition of a ligand to the model confirms the hypothesis that flat planar chromophores stabilize G-quadruplex structures by making them less flexible.

Nucleotide-mediated conformational changes of monomeric actin and Arp3 studied by molecular dynamics simulations

Members of the actin family of proteins exhibit different biochemical properties when ATP, ADP-P(i), ADP, or no nucleotide is bound. We used molecular dynamics simulations to study the effect of nucleotides on the behavior of actin and actin-related protein 3 (Arp3). In all of the actin simulations, the nucleotide cleft stayed closed, as in most crystal structures. ADP was much more mobile within the cleft than ATP, despite the fact that both nucleotides adopt identical conformations in actin crystal structures. The nucleotide cleft of Arp3 opened in most simulations with ATP, ADP, and no bound nucleotide. Deletion of a C-terminal region of Arp3 that extends beyond the conserved actin sequence reduced the tendency of the Arp3 cleft to open. When the Arp3 cleft opened, we observed multiple instances of partial release of the nucleotide. Cleft opening in Arp3 also allowed us to observe correlated movements of the phosphate clamp, cleft mouth, and barbed-end groove, providing a way for changes in the nucleotide state to be relayed to other parts of Arp3. The DNase binding loop of actin was highly flexible regardless of the nucleotide state. The conformation of Ser14/Thr14 in the P1 loop was sensitive to the presence of the gamma-phosphate, but other changes observed in crystal structures were not correlated with the nucleotide state on nanosecond timescales. The divalent cation occupied three positions in the nucleotide cleft, one of which was not previously observed in actin or Arp2/3 complex structures. In sum, these simulations show that subtle differences in structures of actin family proteins have profound effects on their nucleotide-driven behavior.

Exploration structurale du paysage conformationnel de l’acétylcholinestérase par cristallographie cinétique

Acetylcholinesterase is a very rapid enzyme, essential in the process of nerve impulse transmission at cholinergic synapses. It is the target of all currently approved anti-Alzheimer drugs and further progress in the modulation of its activity requires structural as well as dynamical information. Exploration of the conformational energy landscape of a protein by means of X-ray crystallography requires the use of experimental tricks, to overcome the inherently static nature of crystallographic structures. Here we report three experimental approaches that allowed to gain structural insight into the dynamics of acetylcholinesterase, which is relevant for structure-based drug design.

A model study of interactions between TcAChE peripheral site segment Tyr70Val71 and loop 1 of Fasciculin 2

The continuous chain of residues (Thr7 to Ala12) of Loop1 of Fas2 (F1) and its interaction with the peripheral binding sites (Tyr70-Val71) of AChE (P1) has been studied. Our results suggest that the flexibility of Loop1 might be caused by either the partially protonated guanidine group of Arg11 under experimental conditions or by the interaction with the negatively charged center of substrates. The binding energy of F1-P1 is predicted to be -16.6 kcal/mol at the B3LYP/6-311G(d,p) level, which is assumed to originate from one isolated O7...HN10 H-bond, one possible O10...HC71 unconventional O...HC type H-bonding, and the improved pi-bonding cooperativity around the peptide group of the AChE segment Tyr70-Val71. The classical Kitaura-Morokuma energy decomposition analysis, the NPA charge analysis, and the AIM analysis consistently reveal that the peptide group in segment P1 is more polarizable, which might play the key role in the interactions between F1 and P1. The PCM solvent effect corrected results reveal decrease of the interaction energy of the considered model. The importance of Thr8 of Fas2 in the P-site binding of AChE is also concluded. Site-directed mutations on either the Fas2 residue of Thr8 or the AChE residue of Tyr70 are expected to alter the binding behavior of the Loop1 of Fas2 with AChE.

Protein complex formation by acetylcholinesterase and the neurotoxin fasciculin-2 appears to involve an induced-fit mechanism

Specific, rapid association of protein complexes is essential for all forms of cellular existence. The initial association of two molecules in diffusion-controlled reactions is often influenced by the electrostatic potential. Yet, the detailed binding mechanisms of proteins highly depend on the particular system. A complete protein complex formation pathway has been delineated by using structural information sampled over the course of the transformation reaction. The pathway begins at an encounter complex that is formed by one of the apo forms of neurotoxin fasciculin-2 (FAS2) and its high-affinity binding protein, acetylcholinesterase (AChE), followed by rapid conformational rearrangements into an intermediate complex that subsequently converts to the final complex as observed in crystal structures. Formation of the intermediate complex has also been independently captured in a separate 20-ns molecular dynamics simulation of the encounter complex. Conformational transitions between the apo and liganded states of FAS2 in the presence and absence of AChE are described in terms of their relative free energy profiles that link these two states. The transitions of FAS2 after binding to AChE are significantly faster than in the absence of AChE; the energy barrier between the two conformational states is reduced by half. Conformational rearrangements of FAS2 to the final liganded form not only bring the FAS2/AChE complex to lower energy states, but by controlling transient motions that lead to opening or closing one of the alternative passages to the active site of the enzyme also maximize the ligand's inhibition of the enzyme.

Three hydrolases and a transferase: comparative analysis of active-site dynamics via the BioSimGrid database

Comparative molecular dynamics (MD) simulations enable us to explore the conformational dynamics of the active sites of distantly related enzymes. We have used the BioSimGrid (http://www.biosimgrid.org) database to facilitate such a comparison. Simulations of four enzymes were analyzed. These included three hydrolases and a transferase, namely acetylcholinesterase, outer-membrane phospholipase A, outer-membrane protease T, and PagP (an outer-membrane enzyme which transfers a palmitate chain from a phospholipid to lipid A). A set of 17 simulations were analyzed corresponding to a total of approximately 0.1 micros simulation time. A simple metric for active-site integrity was used to demonstrate the existence of clusters of dynamic conformational behaviour of the active sites. Small (i.e. within a cluster) fluctuations appear to be related to the function of an enzymatically active site. Larger fluctuations (i.e. between clusters) correlate with transitions between catalytically active and inactive states. Overall, these results demonstrate the potential of a comparative MD approach to analysis of enzyme function. This approach could be extended to a wider range of enzymes using current high throughput MD simulation and database methods.

Molecular dynamics simulations of human butyrylcholinesterase

Herein, we present results from molecular dynamics (MD) simulations of the human butyrylcholinesterase (BuChE) enzyme in aqueous solution. Two configurations of the unbound form of BuChE differing in the presence or absence of a sodium ion inside the protein gorge were simulated for 10 and 5 ns, respectively. Besides complementing the structural information provided by X-ray data, the MD simulations give insight into the structure of the native BuChE enzyme. For example, it is shown that: the nucleophilic Ser(198) residue and the various binding subsites in the BuChE catalytic cavity are readily accessible from the exterior of the protein; the presence of the sodium ion dynamically explores two different binding sites in the gorge leading to the active site and stabilizes the productive conformation of the Glu(325)/His(438)/Ser(198) catalytic triad; several long-lived water bridges are fully integrated into the architecture of the active site; the positions of the residues at the rim of the gorge region display large deviations with respect to the crystal structure; and two side doors, constituted by residues situated at the tip of the acyl- and Omega-loops, respectively, open wide enough to allow the passage of water molecules. In conclusion, we compare our theoretical results with those from previous work on mouse acetylcholinesterase and discuss their implications for substrate binding and catalysis in BuChE.

Molecular docking study on the “back door” hypothesis for product clearance in acetylcholinesterase

Acetylcholinesterase (AChE) is one of the fastest enzymes known, even though the active site is buried inside the protein at the end of a 20-A deep narrow gorge. Among the great variety of crystal structures of this enzyme, both in the absence and presence of various ligands and proteins, the structure of a complex of AChE with the pseudo-irreversible inhibitor Mf268 is of particular interest, as it assists in the proposal of a back door for product clearance from the active site. Binding of Mf268 to AChE results in the carbamoylation of Ser200 and liberation of an eseroline-fragment as the leaving group. The crystal structure of the AChE-Mf268 complex, however, proves that eseroline has escaped from the enzyme, despite the fact that the Ser-bound inhibitor fragment blocks the gorge entrance. The existence of alternative routes other than through the gorge for product clearance has been postulated but is still controversially discussed in the literature, as an experimental proof for such a back door is still missing. We have used Monte Carlo-based molecular docking methods in order to examine possible alternative pathways that could allow eseroline to be released from the protein after being cleaved from the substrate by Ser200. Based on our results, a short channel at the bottom of the gorge seems to be the most probable back-door site, which begins at amino acid Trp84 and ends at the enzyme surface in a cavity close to amino acid Glu445. [Figure: see text].

Effects of soman inhibition and of structural differences on cholinesterase molecular dynamics: a neutron scattering study

Incoherent elastic neutron scattering experiments on members of the cholinesterase family were carried out to investigate how molecular dynamics is affected by covalent inhibitor binding and by differences in primary and quaternary structure. Tetrameric native and soman-inhibited human butyrylcholinesterase (HuBChE) as well as native dimeric Drosophila melanogaster acetylcholinesterase (DmAChE) hydrated protein powders were examined. Atomic mean-square displacements (MSDs) were found to be identical for native HuBChE and for DmAChE in the whole temperature range examined, leading to the conclusion that differences in activity and substrate specificity are not reflected by a global modification of subnanosecond molecular dynamics. MSDs of native and soman-inhibited HuBChE were identical below the thermal denaturation temperature of the native enzyme, indicating a common mean free-energy surface. Denaturation of the native enzyme is reflected by a relative increase of MSDs consistent with entropic stabilization of the unfolded state. The results suggest that the stabilization of HuBChE phosphorylated by soman is due to an increase in free energy of the unfolded state due to a decrease in entropy.

Conformational dynamics of the ligand-binding domain of inward rectifier K channels as revealed by molecular dynamics simulations: toward an understanding of Kir channel gating

Inward rectifier (Kir) potassium channels are characterized by two transmembrane helices per subunit, plus an intracellular C-terminal domain that controls channel gating in response to changes in concentration of various ligands. Based on the crystal structure of the tetrameric C-terminal domain of Kir3.1, it is possible to build a homology model of the ATP-binding C-terminal domain of Kir6.2. Molecular dynamics simulations have been used to probe the dynamics of Kir C-terminal domains and to explore the relationship between their dynamics and possible mechanisms of channel gating. Multiple simulations, each of 10 ns duration, have been performed for Kir3.1 (crystal structure) and Kir6.2 (homology model), in both their monomeric and tetrameric forms. The Kir6.2 simulations were performed with and without bound ATP. The results of the simulations reveal comparable conformational stability for the crystal structure and the homology model. There is some decrease in conformational flexibility when comparing the monomers with the tetramers, corresponding mainly to the subunit interfaces in the tetramer. The beta-phosphate of ATP interacts with the side chain of K185 in the Kir6.2 model and simulations. The flexibility of the Kir6.2 tetramer is not changed greatly by the presence of bound ATP, other than in two loop regions. Principal components analysis of the simulated dynamics suggests loss of symmetry in both the Kir3.1 and Kir6.2 tetramers, consistent with "dimer-of-dimers" motion of subunits in C-terminal domains of the corresponding Kir channels. This is suggestive of a gating model in which a transition between exact tetrameric symmetry and dimer-of-dimers symmetry is associated with a change in transmembrane helix packing coupled to gating of the channel. Dimer-of-dimers motion of the C-terminal domain tetramer is also supported by coarse-grained (anisotropic network model) calculations. It is of interest that loss of exact rotational symmetry has also been suggested to play a role in gating in the bacterial Kir homolog, KirBac1.1, and in the nicotinic acetylcholine receptor channel.

Molecular dynamics simulation of the M2 helices within the nicotinic acetylcholine receptor transmembrane domain: structure and collective motions

Multiple nanosecond duration molecular dynamics simulations were performed on the transmembrane region of the Torpedo nicotinic acetylcholine receptor embedded within a bilayer mimetic octane slab. The M2 helices and M2-M3 loop regions were free to move, whereas the outer (M1, M3, M4) helix bundle was backbone restrained. The M2 helices largely retain their hydrogen-bonding pattern throughout the simulation, with some distortions in the helical end and loop regions. All of the M2 helices exhibit bending motions, with the hinge point in the vicinity of the central hydrophobic gate region (corresponding to residues alphaL251 and alphaV255). The bending motions of the M2 helices lead to a degree of dynamic narrowing of the pore in the region of the proposed hydrophobic gate. Calculations of Born energy profiles for various structures along the simulation trajectory suggest that the conformations of the M2 bundle sampled correspond to a closed conformation of the channel. Principal components analyses of each of the M2 helices, and of the five-helix M2 bundle, reveal concerted motions that may be relevant to channel function. Normal mode analyses using the anisotropic network model reveal collective motions similar to those identified by principal components analyses.

Molecular dynamics simulations of evolved collective motions of atoms in the myosin motor domain upon perturbation of the ATPase pocket

A crucial point for mechanical force generation in actomyosin systems is how the energy released by ATP hydrolysis in the myosin motor domain gives rise to the movement of the myosin head along the actin filament. We assumed the signal of the ATP hydrolysis to be transmitted as modulated atomic vibrations from the nucleotide-binding site throughout the myosin head, and carried out 1-ns all-atom molecular dynamics simulations for that signal transmission. We distributed the released energy to atoms located around the ATPase pocket as kinetic energies and examined how the effect of disturbance extended throughout the motor domain. The result showed that the disturbance signal extended over the motor domain in 150 ps and induced slowly varying collective motions of atoms at the actin-binding site and the junction with the neck, both of which are relevant to the movement of the myosin head along the actin filament. We also performed a principal component analysis of thermal atomic motions for the motor domain, and the first principal component was consistent with the response to the disturbance given to the ATPase pocket.

The dynamics of ligand barrier crossing inside the acetylcholinesterase gorge

The dynamics of ligand movement through the constricted region of the acetylcholinesterase gorge is important in understanding how the ligand gains access to and is released from the active site of the enzyme. Molecular dynamics simulations of the simple ligand, tetramethylammonium, crossing this bottleneck region are conducted using umbrella potential sampling and activated flux techniques. The low potential of mean force obtained is consistent with the fast reaction rate of acetylcholinesterase observed experimentally. From the results of the activated dynamics simulations, local conformational fluctuations of the gorge residues and larger scale collective motions of the protein are found to correlate highly with the ligand crossing.

Acetylcholinesterase in motion: visualizing conformational changes in crystal structures by a morphing procedure

In order to visualize and appreciate conformational changes between homologous three-dimensional (3D) protein structures or protein/inhibitor complexes, we have developed a user-friendly morphing procedure. It enabled us to detect coordinated conformational changes not easily discernible by analytic methods or by comparison of static images. This procedure was applied to comparison of native Torpedo californica acetylcholinesterase and of complexes with reversible inhibitors and conjugates with covalent inhibitors. It was likewise shown to be valuable for the visualization of conformational differences between acetylcholinesterases from different species. The procedure involves generation, in Cartesian space, of 25 interpolated intermediate structures between the initial and final 3D structures, which then serve as the individual frames in a QuickTime movie.

Structural and dynamic properties of water around acetylcholinesterase

Structural and dynamic properties of water molecules around acetylcholinesterase are examined from a 10-nsec molecular dynamics simulation to help understand how the protein alters water properties. Water structure is broken down into hydration sites constructed from the water density <3.6 A from the protein surface. These sites are characterized according to occupancy, number of water neighbors, hydrogen bonds, dipole moment, and residence time. The site description provides a convenient means to describe the extent and localization of these properties. Determining the network of paths that waters follow from site to site and measuring the rate of flow of waters from the sites to the bulk make it possible to quantitatively study the time scales and paths that water molecules follow as they move around the protein.