Acetylcholinesterase inhibition by fasciculin: crystal structure of the complex

Molecular Basis for Mambalgin-2 Interaction with Heterotrimeric α-ENaC/ASIC1a/γ-ENaC Channels in Cancer Cells

Cancer progression is characterized by microenvironmental acidification. Tumor cells adapt to low environmental pH by activating acid-sensing trimeric ion channels of the DEG/ENaC family. The α-ENaC/ASIC1a/γ-ENaC heterotrimeric channel is a tumor-specific acid-sensing channel, and its targeting can be considered a new strategy for cancer therapy. Mambalgin-2 from the Dendroaspis polylepis venom inhibits the α-ENaC/ASIC1a/γ-ENaC heterotrimer more effectively than the homotrimeric ASIC1a channel, initially proposed as the target of mambalgin-2. Although the molecular basis of such mambalgin selectivity remained unclear. Here, we built the models of the complexes of mambalgin-2 with the α-ENaC/ASIC1a/γ-ENaC and ASIC1a channels, performed MD and predicted the difference in the binding modes. The importance of the 'head' loop region of mambalgin-2 for the interaction with the hetero-, but not with the homotrimeric channel was confirmed by site-directed mutagenesis and electrophysiology. A new mode of allosteric regulation of the ENaC channels by linking the thumb domain of the ASIC1a subunit with the palm domain of the γ-ENaC subunit was proposed. The data obtained provide new insights into the regulation of various types of acid-sensing ion channels and the development of new strategies for cancer treatment.

Toxic Effects of Koumine on the Early-Life Development Stage of Zebrafish

Koumine is one of the most abundant alkaloids found in Gelsemium elegans, and it has a wide range of pharmacological effects including antitumor, anti-inflammatory, analgesic treatment effects, and antianxiety. However, its high toxicity and unclear mechanism of action have greatly limited the medicinal development and use of koumine. We investigated the toxic effects of koumine on the developmental toxicity and behavioral neurotoxicity of zebrafish embryos and larvae. Embryos at 6 h postfertilization (hpf) were exposed to 12.5, 25, 50, 75, and 100 mg/L of koumine until 120 hpf. Koumine affected the hatching and heartbeats of the embryos. The morphological analysis also revealed many abnormalities, such as shortened bodies, yolk sac edemas, tail malformations, and pericardial edemas. To identify the neurotoxicity of koumine, the behavior of the larvae was measured. Koumine at 50 and 100 mg/L affect the escape response. The embryos exhibited uncoordinated muscle contractions along the body axis in response to touch at 36 hpf. More importantly, we found that the neurotoxicity of koumine is mainly caused by influencing the ACh content and the activity of AChE without impairing motor neuron development. A comprehensive analysis shows that a high concentration of koumine has obvious toxic effects on zebrafish, and the safe concentration of koumine for zebrafish should be less than 25 mg/L. These results will be valuable for better understanding the toxicity of koumine and provide new insights into the application of koumine.

MT9, a natural peptide from black mamba venom antagonizes the muscarinic type 2 receptor and reverses the M2R-agonist-induced relaxation in rat and human arteries

All five muscarinic receptors have important physiological roles. The endothelial M2 and M3 subtypes regulate arterial tone through direct coupling to Gq or Gi/o proteins. Yet, we lack selective pharmacological drugs to assess the respective contribution of muscarinic receptors to a given function. We used mamba snake venoms to identify a selective M2R ligand to investigate its contribution to arterial contractions. Using a bio-guided screening binding assay, we isolated MT9 from the black mamba venom, a three-finger toxin active on the M2R subtype. After sequencing and chemical synthesis of MT9, we characterized its structure by X-ray diffraction and determined its pharmacological characteristics by binding assays, functional tests, and ex vivo experiments on rat and human arteries. Although MT9 belongs to the three-finger fold toxins family, it is phylogenetically apart from the previously discovered muscarinic toxins, suggesting that two groups of peptides evolved independently and in a convergent way to target muscarinic receptors. The affinity of MT9 for the M2R is 100 times stronger than that for the four other muscarinic receptors. It also antagonizes the M2R/G_i pathways in cell-based assays. MT9 acts as a non-competitive antagonist against acetylcholine or arecaine, with low nM potency, for the activation of isolated rat mesenteric arteries. These results were confirmed on human internal mammary arteries. In conclusion, MT9 is the first fully characterized M2R-specific natural toxin. It should provide a tool for further understanding of the effect of M2R in various arteries and may position itself as a new drug candidate in cardio-vascular diseases.

Structural and Functional Diversity of Animal Toxins Interacting With GPCRs

Peptide toxins from venoms have undergone a long evolutionary process allowing host defense or prey capture and making them highly selective and potent for their target. This has resulted in the emergence of a large panel of toxins from a wide diversity of species, with varied structures and multiple associated biological functions. In this way, animal toxins constitute an inexhaustible reservoir of druggable molecules due to their interesting pharmacological properties. One of the most interesting classes of therapeutic targets is the G-protein coupled receptors (GPCRs). GPCRs represent the largest family of membrane receptors in mammals with approximately 800 different members. They are involved in almost all biological functions and are the target of almost 30% of drugs currently on the market. Given the interest of GPCRs in the therapeutic field, the study of toxins that can interact with and modulate their activity with the purpose of drug development is of particular importance. The present review focuses on toxins targeting GPCRs, including peptide-interacting receptors or aminergic receptors, with a particular focus on structural aspects and, when relevant, on potential medical applications. The toxins described here exhibit a great diversity in size, from 10 to 80 amino acids long, in disulfide bridges, from none to five, and belong to a large panel of structural scaffolds. Particular toxin structures developed here include inhibitory cystine knot (ICK), three-finger fold, and Kunitz-type toxins. We summarize current knowledge on the structural and functional diversity of toxins interacting with GPCRs, concerning first the agonist-mimicking toxins that act as endogenous agonists targeting the corresponding receptor, and second the toxins that differ structurally from natural agonists and which display agonist, antagonist, or allosteric properties.

A new Kunitz-type snake toxin family associated with an original mode of interaction with the vasopressin 2 receptor

BACKGROUND AND PURPOSE

Venomous animals express numerous Kunitz-type peptides. The mambaquaretin-1 (MQ1) identified from the Dendroaspis angusticeps venom is the most selective antagonist of the arginine-vasopressin V2 receptor (V2R) and the unique Kunitz-type peptide active on a GPCR. We aimed to exploit other mamba venoms to enlarge the V2R-Kunitz peptide family and gain insight into the MQ1 molecular mode of action.

EXPERIMENTAL APPROACH

We used a bio-guided screening assay to identify novel MQs and placed them phylogenetically. MQs were produced by solid phase peptide synthesis and characterized in vitro by binding and functional tests and in vivo by diuresis measurement in rats.

KEY RESULTS

Eight additional MQs were identified with nanomolar affinities for the V2R, all antagonists. MQs form a new subgroup in the Kunitz family, close to the V2R non-active dendrotoxins and to 2 V2R active cobra toxins. Sequence comparison between active and non-active V2R Kunitz peptides highlighted 5 positions, belong which, four are involved in V2R interaction and which belong to the 2 large MQ1 loops. We finally determined that 8 positions, part of these 2 loops, interact with the V2R. The variant MQ1-K39A showed higher affinity for the hV2R but not for the rat V2R.

CONCLUSIONS AND IMPLICATIONS

A new function and mode of action is now associated with the Kunitz-peptides. The number of MQ1 residues involved in V2R binding is large and may explain its absolute selectivity. MQ1-K39A represents the first step in the improvement of the MQ1 design for medicinal perspective.

Design and synthesis of new carbamates as inhibitors for fatty acid amide hydrolase and cholinesterases: Molecular dynamic, in vitro and in vivo studies

As anandamide (N-arachidonoylethanolamine, AEA) shows neuroprotective effects, the inhibition of its degradative enzyme, fatty acid amide hydrolase (FAAH) has been considered as a hopeful avenue for the treatment of neurodegenerative diseases, like Alzheimer's disease (AD). Memory loss, cognitive impairment and diminution of the cholinergic tone, due to the dying cholinergic neurons in the basal forebrain, are common hallmarks in patients with AD. By taking advantage of cholinesterase inhibitors (ChEIs), the degradation of acetylcholine (ACh) is decreased leading to enhanced cholinergic neurotransmission in the aforementioned region and ultimately improves the clinical condition of AD patients. In this work, new carbamates were designed as inhibitors of FAAH and cholinestrases (ChEs) (acetylcholinestrase (AChE), butyrylcholinestrase (BuChE)) inspired by the structure of the native substrates, structure of active sites and the SARs of the well-known inhibitors of these enzymes. All the designed compounds were synthesized using different reactions. All the target compounds were tested for their inhibitory activity against FAAH and ChEs by employing the Cayman assay kit and Elman method respectively. Generally, compounds possessing aminomethyl phenyl linker was more potent compared to their corresponding compounds possessing piperazinyl ethyl linker. The inhibitory potential of the compounds 3a-q extended from 0.83 ± 0.03 μM (3i) to ˃100 μM (3a) for FAAH, 0.39 ± 0.02 μM (3i) to 24% inhibition in 113 ± 4.8 μM (3b) for AChE, and 1.8 ± 3.2 μM (3i) to 23.2 ± 0.2 μM (3b) for BuChE. Compound 3i a heptyl carbamate analog possessing 2-oxo-1,2-dihydroquinolin ring and aminomethyl phenyl linker showed the most inhibitory activity against three enzymes. Also, compound 3i was investigated for memory improvement using the Morris water maze test in which the compound showed better memory improvement at 10 mg/kg compared to reference drug rivastigmine at 2.5 mg/kg. Molecular docking and molecular dynamic studies of compound 3i into the enzymes displayed the possible interactions of key residues of the active sites with compound 3i. Finally, kinetic study indicated that 3i inhibits AChE through the mixed- mode mechanism and non-competitive inhibition mechanism was revealed for BuChE.

Analogues of cannabinoids as multitarget drugs in the treatment of Alzheimer's disease

Given that neuronal degeneration in Alzheimer's disease (AD) is caused by the combination of multiple neurotoxic insults, current directions in the research of novel therapies to treat this disease attempts to design multitarget strategies that could be more effective than the simply use of acetylcholinesterase inhibitors; currently, the most used therapy for AD. One option, explored recently, is the synthesis of new analogues of cannabinoids that could competitively inhibit the acetylcholinesterase (AChE) enzyme and showing the classic neuroprotective profile of cannabinoid compounds. In this work, molecular docking has been used to design some cannabinoid analogues with such multitarget properties, based on the similarities of donepezil and Δ⁹-tetrahydrocannabinol. The analogues synthesized, compounds 1 and 2, demonstrated to have two interesting characteristics in different in vitro assays: competitive inhibition of AChE and competitive antagonism at the CB₁/CB₂ receptors. They are highly lipophilic, highlighting that they could easily reach the CNS, and apparently presented a low toxicity. These results open the door to the synthesis of new compounds for a more effective treatment of AD.

The multiple biological roles of the cholinesterases

It is tacitly assumed that the biological role of acetylcholinesterase is termination of synaptic transmission at cholinergic synapses. However, together with its structural homolog, butyrylcholinesterase, it is widely distributed both within and outside the nervous system, and, in many cases, the role of both enzymes remains obscure. The transient appearance of the cholinesterases in embryonic tissues is especially enigmatic. The two enzymes' extra-synaptic roles, which are known as 'non-classical' roles, are the topic of this review. Strong evidence has been presented that AChE and BChE play morphogenetic roles in a variety of eukaryotic systems, and they do so either by acting as adhesion proteins, or as trophic factors. As trophic factors, one mode of action is to directly regulate morphogenesis, such as neurite outgrowth, by poorly understood mechanisms. The other mode is by regulating levels of acetylcholine, which acts as the direct trophic factor. Alternate substrates have been sought for the cholinesterases. Quite recently, it was shown that levels of the aggression hormone, ghrelin, which also controls appetite, are regulated by butyrylcholinesterase. The rapid hydrolysis of acetylcholine by acetylcholinesterase generates high local proton concentrations. The possible biophysical and biological consequences of this effect are discussed. The biological significance of the acetylcholinesterases secreted by parasitic nematodes is reviewed, and, finally, the involvement of acetylcholinesterase in apoptosis is considered.

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.

Comparative mapping of selected structural determinants on the extracellular domains of cholinesterase-like cell-adhesion molecules

Cell adhesion generally involves formation of homophilic or heterophilic protein complexes between two cells to form transcellular junctions. Neural cell-adhesion members of the α/β-hydrolase fold superfamily of proteins use their extracellular or soluble cholinesterase-like domain to bind cognate partners across cell membranes, as illustrated by the neuroligins. These cell-adhesion molecules currently comprise the synaptic organizers neuroligins found in all animal phyla, along with three proteins found only in invertebrates: the guidance molecule neurotactin, the glia-specific gliotactin, and the basement membrane protein glutactin. Although these proteins share a cholinesterase-like fold, they lack one or more residues composing the catalytic triad responsible for the enzymatic activity of the cholinesterases. Conversely, they are found in various subcellular localisations and display specific disulfide bonding and N-glycosylation patterns, along with individual surface determinants possibly associated with recognition and binding of protein partners. Formation of non-covalent dimers typical of the cholinesterases is documented for mammalian neuroligins, yet whether invertebrate neuroligins and their neurotactin, gliotactin and glutactin relatives also form dimers in physiological conditions is unknown. Here we provide a brief overview of the localization, function, evolution, and conserved versus individual structural determinants of these cholinesterase-like cell-adhesion proteins. This article is part of the special issue entitled 'Acetylcholinesterase Inhibitors: From Bench to Bedside to Battlefield'.

Snake three-finger α-neurotoxins and nicotinic acetylcholine receptors: molecules, mechanisms and medicine

Snake venom three-finger α-neurotoxins (α-3FNTx) act on postsynaptic nicotinic acetylcholine receptors (nAChRs) at the neuromuscular junction (NMJ) to produce skeletal muscle paralysis. The discovery of the archetypal α-bungarotoxin (α-BgTx), almost six decades ago, exponentially expanded our knowledge of membrane receptors and ion channels. This included the localisation, isolation and characterization of the first receptor (nAChR); and by extension, the pathophysiology and pharmacology of neuromuscular transmission and associated pathologies such as myasthenia gravis, as well as our understanding of the role of α-3FNTxs in snakebite envenomation leading to novel concepts of targeted treatment. Subsequent studies on a variety of animal venoms have yielded a plethora of novel toxins that have revolutionized molecular biomedicine and advanced drug discovery from bench to bedside. This review provides an overview of nAChRs and their subtypes, classification of α-3FNTxs and the challenges of typifying an increasing arsenal of structurally and functionally unique toxins, and the three-finger protein (3FP) fold in the context of the uPAR/Ly6/CD59/snake toxin superfamily. The pharmacology of snake α-3FNTxs including their mechanisms of neuromuscular blockade, variations in reversibility of nAChR interactions, specificity for nAChR subtypes or for distinct ligand-binding interfaces within a subtype and the role of α-3FNTxs in neurotoxic envenomation are also detailed. Lastly, a reconciliation of structure-function relationships between α-3FNTx and nAChRs, derived from historical mutational and biochemical studies and emerging atomic level structures of nAChR models in complex with α-3FNTxs is discussed.

Assessing the performance of the MM/PBSA and MM/GBSA methods. 10. Impacts of enhanced sampling and variable dielectric model on protein-protein Interactions

Enhanced sampling has been extensively used to capture the conformational transitions in protein folding, but it attracts much less attention in the studies of protein-protein recognition. In this study, we evaluated the impact of enhanced sampling methods and solute dielectric constants on the overall accuracy of the molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) and molecular mechanics/generalized Born surface area (MM/GBSA) approaches for the protein-protein binding free energy calculations. Here, two widely used enhanced sampling methods, including aMD and GaMD, and conventional molecular dynamics (cMD) simulations with two AMBER force fields (ff03 and ff14SB) were used to sample the conformations for 21 protein-protein complexes. The MM/PBSA and MM/GBSA calculation results illustrate that the standard MM/GBSA based on the cMD simulations yields the best Pearson correlation (r_p = -0.523) between the predicted binding affinities and the experimental data, which is much higher than that given by MM/PBSA (r_p = -0.212). Two enhanced sampling methods (aMD and GaMD) are indeed more efficient for conformational sampling, but they did not improve the binding affinity predictions for protein-protein systems, suggesting that the aMD or GaMD sampling (at least in short timescale simulations) may not be a good choice for the MM/PBSA and MM/GBSA predictions of protein-protein complexes. The solute dielectric constant of 1.0 is recommended to MM/GBSA, but a higher solute dielectric constant is recommended to MM/PBSA, especially for the systems with higher polarity on the protein-protein binding interfaces. Then, a preliminary assessment of the MM/GBSA calculations based on a variable dielectric generalized Born (VDGB) model was conducted. The results highlight the potential power of VDGB in the free energy predictions for protein-protein systems, but more thorough studies should be done in the future.

Metal based donepezil analogues designed to inhibit human acetylcholinesterase for Alzheimer's disease

Among neurodegenerative disorders, Alzheimer's disease (AD) is one of the most common disorders showing slow progressive cognitive decline. Targeting acetylcholinesterase (AChE) is one of the major strategies for AD therapeutics, as cholinergic pathways in the cerebral cortex and basal forebrain are compromised. Herein, we report the design of some copper and other metal based donepezil derivatives, employing density functional theory (DFT). All designed compounds are optimized at the B3LYP/SDD level of theory. Dipole moments, electronic energie, enthalpies, Gibbs free energies, and HOMO-LUMO gaps of these modified compounds are also investigated in the subsequent analysis. The molecules were then subjected to molecular docking analysis with AChE to study the molecular interactions broadly. Ensemble based docking and molecular dynamics (MD) simulations of the best candidates were also performed. Docking and MD simulation reveal that modified drugs are more potent than unmodified donepezil, where Trp86, Tyr337, Phe330 residues play some important roles in drug-receptor interactions. According to ensemble based docking, D9 shows greater binding affinity compared to the parent in most conformations obtained from protein data bank and MD simulation. In addition, it is observed that the π- π stacking with the residues of Trp86, Tyr337, Tyr341, Tyr124 and Trp286 may be required for strong ligand binding. Moreover, ADME/T analysis suggests that modified derivatives are less toxic and have improved pharmacokinetic properties than those of the parent drug. These results further confirm the ability of metal-directed drugs to bind simultaneously to the active sites of AChE and support them as potential candidates for the future treatment of Alzheimer's disease.

Characterization of the interactions between coumarin-derivatives and acetylcholinesterase: Examination by NMR and docking simulations

Alzheimer's disease (AD) is one of the most common forms of dementia and a significant threat to the elderly populations, especially in the Western world. The rapid hydrolysis of the principal neurotransmitter into choline and acetate by acetylcholinesterase (AChE) at synapses causes the loss of cognitive response that becomes the real cause of AD. Therefore, inhibition of AChE is the most fundamental therapy among currently available treatments for AD. In this context, we designed and performed molecular recognitions studies of coumarin-based inhibitors towards AChE. STD NMR and Tr-NOESY applications were utilized to evaluate the binding epitope, the dissociation constant (K_D) and bound conformations of these inhibitors within this inhibitor-AChE complex. Compound 1, which has a similar inhibition activity to tacrine (a current drug) led in this study as a stronger binder with K_D = 30 μM ,even greater than tacrine (K_D = 140 μM). Moreover, docking simulations mimic NMR results and provided evidence of synchronizing binding of compound 1 with three sites; the peripheral anionic site, the bottom of the gorge, and the catalytic site. Therefore, we envisioned from our experimental and theoretical results that coumarin-based inhibitors containing a piperidinyl scaffold might be a potential drug candidates for AD in the future.

Wireless control of cellular function by activation of a novel protein responsive to electromagnetic fields

The Kryptopterus bicirrhis (glass catfish) is known to respond to electromagnetic fields (EMF). Here we tested its avoidance behavior in response to static and alternating magnetic fields stimulation. Using expression cloning we identified an electromagnetic perceptive gene (EPG) from the K. bicirrhis encoding a protein that responds to EMF. This EPG gene was cloned and expressed in mammalian cells, neuronal cultures and in rat’s brain. Immunohistochemistry showed that the expression of EPG is confined to the mammalian cell membrane. Calcium imaging in mammalian cells and cultured neurons expressing EPG demonstrated that remote activation by EMF significantly increases intracellular calcium concentrations, indicative of cellular excitability. Moreover, wireless magnetic activation of EPG in rat motor cortex induced motor evoked responses of the contralateral forelimb in vivo. Here we report on the development of a new technology for remote, non-invasive modulation of cell function.

Three-Finger Proteins from the Ly6/uPAR Family: Functional Diversity within One Structural Motif

The discovery in higher animals of proteins from the Ly6/uPAR family, which have structural homology with snake "three-finger" neurotoxins, has generated great interest in these molecules and their role in the functioning of the organism. These proteins have been found in the nervous, immune, endocrine, and reproductive systems of mammals. There are two types of the Ly6/uPAR proteins: those associated with the cell membrane by GPI-anchor and secreted ones. For some of them (Lynx1, SLURP-1, SLURP-2, Lypd6), as well as for snake α-neurotoxins, the target of action is nicotinic acetylcholine receptors, which are widely represented in the central and peripheral nervous systems, and in many other tissues, including epithelial cells and the immune system. However, the targets of most proteins from the Ly6/uPAR family and the mechanism of their action remain unknown. This review presents data on the structural and functional properties of the Ly6/uPAR proteins, which reveal a variety of functions within a single structural motif.

Assessing the performance of the MM/PBSA and MM/GBSA methods. 6. Capability to predict protein-protein binding free energies and re-rank binding poses generated by protein-protein docking

Understanding protein-protein interactions (PPIs) is quite important to elucidate crucial biological processes and even design compounds that interfere with PPIs with pharmaceutical significance. Protein-protein docking can afford the atomic structural details of protein-protein complexes, but the accurate prediction of the three-dimensional structures for protein-protein systems is still notoriously difficult due in part to the lack of an ideal scoring function for protein-protein docking. Compared with most scoring functions used in protein-protein docking, the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) and Molecular Mechanics/Poisson Boltzmann Surface Area (MM/PBSA) methodologies are more theoretically rigorous, but their overall performance for the predictions of binding affinities and binding poses for protein-protein systems has not been systematically evaluated. In this study, we first evaluated the performance of MM/PBSA and MM/GBSA to predict the binding affinities for 46 protein-protein complexes. On the whole, different force fields, solvation models, and interior dielectric constants have obvious impacts on the prediction accuracy of MM/GBSA and MM/PBSA. The MM/GBSA calculations based on the ff02 force field, the GB model developed by Onufriev et al. and a low interior dielectric constant (εin = 1) yield the best correlation between the predicted binding affinities and the experimental data (rp = -0.647), which is better than MM/PBSA (rp = -0.523) and a number of empirical scoring functions used in protein-protein docking (rp = -0.141 to -0.529). Then, we examined the capability of MM/GBSA to identify the possible near-native binding structures from the decoys generated by ZDOCK for 43 protein-protein systems. The results illustrate that the MM/GBSA rescoring has better capability to distinguish the correct binding structures from the decoys than the ZDOCK scoring. Besides, the optimal interior dielectric constant of MM/GBSA for re-ranking docking poses may be determined by analyzing the characteristics of protein-protein binding interfaces. Considering the relatively high prediction accuracy and low computational cost, MM/GBSA may be a good choice for predicting the binding affinities and identifying correct binding structures for protein-protein systems.

Hot Spots for Protein Partnerships at the Surface of Cholinesterases and Related α/β Hydrolase Fold Proteins or Domains-A Structural Perspective

The hydrolytic enzymes acetyl- and butyryl-cholinesterase, the cell adhesion molecules neuroligins, and the hormonogenic macromolecule thyroglobulin are a few of the many members of the α/β hydrolase fold superfamily of proteins. Despite their distinctive functions, their canonical subunits, with a molecular surface area of ~20,000 Ų, they share binding patches and determinants for forming homodimers and for accommodating structural subunits or protein partners. Several of these surface regions of high functional relevance have been mapped through structural or mutational studies, while others have been proposed based on biochemical data or molecular docking studies. Here, we review these binding interfaces and emphasize their specificity versus potentially multifunctional character.

Structural Mechanism for Modulation of Synaptic Neuroligin-Neurexin Signaling by MDGA Proteins

Neuroligin-neurexin (NL-NRX) complexes are fundamental synaptic organizers in the central nervous system. An accurate spatial and temporal control of NL-NRX signaling is crucial to balance excitatory and inhibitory neurotransmission, and perturbations are linked with neurodevelopmental and psychiatric disorders. MDGA proteins bind NLs and control their function and interaction with NRXs via unknown mechanisms. Here, we report crystal structures of MDGA1, the NL1-MDGA1 complex, and a spliced NL1 isoform. Two large, multi-domain MDGA molecules fold into rigid triangular structures, cradling a dimeric NL to prevent NRX binding. Structural analyses guided the discovery of a broad, splicing-modulated interaction network between MDGA and NL family members and helped rationalize the impact of autism-linked mutations. We demonstrate that expression levels largely determine whether MDGAs act selectively or suppress the synapse organizing function of multiple NLs. These results illustrate a potentially brain-wide regulatory mechanism for NL-NRX signaling modulation.

•

The MDGA1 extracellular region has an unusual triangular multi-domain arrangement

•

The NL1-MDGA1 complex structure reveals how MDGA proteins block neurexin binding

•

MDGA1 and MDGA2 bind all NL isoforms, a process fine-tuned by alternative splicing

•

MDGA1 and MDGA2 suppress NL synaptogenic activity in a concentration-dependent manner

Schwann Cells in Neuromuscular Junction Formation and Maintenance

The neuromuscular junction (NMJ) is a tripartite synapse that is formed by motor nerve terminals, postjunctional muscle membranes, and terminal Schwann cells (TSCs) that cover the nerve-muscle contact. NMJ formation requires intimate communications among the three different components. Unlike nerve-muscle interaction, which has been well characterized, less is known about the role of SCs in NMJ formation and maintenance. We show that SCs in mice lead nerve terminals to prepatterned AChRs. Ablating SCs at E8.5 (i.e., prior nerve arrival at the clusters) had little effect on aneural AChR clusters at E13.5, suggesting that SCs may not be necessary for aneural clusters. SC ablation at E12.5, a time when phrenic nerves approach muscle fibers, resulted in smaller and fewer nerve-induced AChR clusters; however, SC ablation at E15.5 reduced AChR cluster size but had no effect on cluster density, suggesting that SCs are involved in AChR cluster maturation. Miniature endplate potential amplitude, but not frequency, was reduced when SCs were ablated at E15.5, suggesting that postsynaptic alterations may occur ahead of presynaptic deficits. Finally, ablation of SCs at P30, after NMJ maturation, led to NMJ fragmentation and neuromuscular transmission deficits. Miniature endplate potential amplitude was reduced 3 d after SC ablation, but both amplitude and frequency were reduced 6 d after. Together, these results indicate that SCs are not only required for NMJ formation, but also necessary for its maintenance; and postsynaptic function and structure appeared to be more sensitive to SC ablation.

SIGNIFICANCE STATEMENT

Neuromuscular junctions (NMJs) are critical for survival and daily functioning. Defects in NMJ formation during development or maintenance in adulthood result in debilitating neuromuscular disorders. The role of Schwann cells (SCs) in NMJ formation and maintenance was not well understood. We genetically ablated SCs during development and after NMJ formation to investigate the consequences of the ablation. This study reveals a critical role of SCs in NMJ formation as well as maintenance.

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 three-finger toxin fold: a multifunctional structural scaffold able to modulate cholinergic functions

Three-finger fold toxins are miniproteins frequently found in Elapidae snake venoms. This fold is characterized by three distinct loops rich in β-strands and emerging from a dense, globular core reticulated by four highly conserved disulfide bridges. The number and diversity of receptors, channels, and enzymes identified as targets of three-finger fold toxins is increasing continuously. Such manifold diversity highlights the specific adaptability of this fold for generating pleiotropic functions. Although this toxin superfamily disturbs many biological functions by interacting with a large diversity of molecular targets, the most significant target is the cholinergic system. By blocking the activity of the nicotinic and muscarinic acetylcholine receptors or by inhibiting the enzyme acetylcholinesterase, three-finger fold toxins interfere most drastically with neuromuscular junction functioning. Several of these toxins have become powerful pharmacological tools for studying the function and structure of their molecular targets. Most importantly, since dysfunction of these receptors/enzyme is involved in many diseases, exploiting the three-finger scaffold to create novel, highly specific therapeutic agents may represent a major future endeavor. This is an article for the special issue XVth International Symposium on Cholinergic Mechanisms.

Green synthesis of novel spiro-indenoquinoxaline derivatives and their cholinesterases inhibition activity

A convenient synthesis of substituted spiroindenoquinoxalines at mild and green conditions was developed. Multicomponent reaction of substituted phenylene diamines, ninhydrin, malononitrile and N,N'-substituted-2-nitroethene-1,1-diamines produced the target compounds. Twelve new spiroindenoquinoxalines were obtained, and their ability in inhibition of acetyl and butyrylcholinesterases were investigated both in vitro and in silico. All compounds showed moderate level activity against both acetyl and butyrylcholinesterases.

JET2 Viewer: a database of predicted multiple, possibly overlapping, protein-protein interaction sites for PDB structures

The database JET2 Viewer, openly accessible at http://www.jet2viewer.upmc.fr/, reports putative protein binding sites for all three-dimensional (3D) structures available in the Protein Data Bank (PDB). This knowledge base was generated by applying the computational method JET² at large-scale on more than 20 000 chains. JET² strategy yields very precise predictions of interacting surfaces and unravels their evolutionary process and complexity. JET2 Viewer provides an online intelligent display, including interactive 3D visualization of the binding sites mapped onto PDB structures and suitable files recording JET² analyses. Predictions were evaluated on more than 15 000 experimentally characterized protein interfaces. This is, to our knowledge, the largest evaluation of a protein binding site prediction method. The overall performance of JET² on all interfaces are: Sen = 52.52, PPV = 51.24, Spe = 80.05, Acc = 75.89. The data can be used to foster new strategies for protein-protein interactions modulation and interaction surface redesign.

Unbinding of fluorinated oxime drug from the AChE gorge in polarizable water: a well-tempered metadynamics study

A well-tempered metadynamics study reveals that fluorinated obidoxime is held more firmly in the AChE gorge in comparison to obidoxime.

Probing the origins of human acetylcholinesterase inhibition via QSAR modeling and molecular docking

Alzheimer's disease (AD) is a chronic neurodegenerative disease which leads to the gradual loss of neuronal cells. Several hypotheses for AD exists (e.g., cholinergic, amyloid, tau hypotheses, etc.). As per the cholinergic hypothesis, the deficiency of choline is responsible for AD; therefore, the inhibition of AChE is a lucrative therapeutic strategy for the treatment of AD. Acetylcholinesterase (AChE) is an enzyme that catalyzes the breakdown of the neurotransmitter acetylcholine that is essential for cognition and memory. A large non-redundant data set of 2,570 compounds with reported IC50 values against AChE was obtained from ChEMBL and employed in quantitative structure-activity relationship (QSAR) study so as to gain insights on their origin of bioactivity. AChE inhibitors were described by a set of 12 fingerprint descriptors and predictive models were constructed from 100 different data splits using random forest. Generated models afforded R (2), [Formula: see text] and [Formula: see text] values in ranges of 0.66-0.93, 0.55-0.79 and 0.56-0.81 for the training set, 10-fold cross-validated set and external set, respectively. The best model built using the substructure count was selected according to the OECD guidelines and it afforded R (2), [Formula: see text] and [Formula: see text] values of 0.92 ± 0.01, 0.78 ± 0.06 and 0.78 ± 0.05, respectively. Furthermore, Y-scrambling was applied to evaluate the possibility of chance correlation of the predictive model. Subsequently, a thorough analysis of the substructure fingerprint count was conducted to provide informative insights on the inhibitory activity of AChE inhibitors. Moreover, Kennard-Stone sampling of the actives were applied to select 30 diverse compounds for further molecular docking studies in order to gain structural insights on the origin of AChE inhibition. Site-moiety mapping of compounds from the diversity set revealed three binding anchors encompassing both hydrogen bonding and van der Waals interaction. Molecular docking revealed that compounds 13, 5 and 28 exhibited the lowest binding energies of -12.2, -12.0 and -12.0 kcal/mol, respectively, against human AChE, which is modulated by hydrogen bonding, π-π stacking and hydrophobic interaction inside the binding pocket. These information may be used as guidelines for the design of novel and robust AChE inhibitors.

Perspectives for the structure-based design of acetylcholinesterase reactivators

Rational design of active molecules through structure-based methods has been gaining adepts during the last decades due to the wider availability of protein structures, most of them conjugated with relevant ligands. Acetylcholinesterase (AChE) is a molecular target with a considerable amount of data related to its sequence and 3-dimensional structure. In addition, there are structural insights about the mechanism of action of the natural substrate and drugs used in Alzheimer's disease, organophosphorus compounds, among others. We looked for AChE structural data useful for in silico design of potential interacting molecules. In particular, we focused on information regarding the design of ligands aimed to reactivate AChE catalytic activity. The structures of 178 AChE were annotated and categorized on different subsets according to the nature of the ligand, source organisms and experimental details. We compared sequence homology among the active site from Torpedo californica, Mus musculus and Homo sapiens with the latter two species having the closest relationship (88.9% identity). In addition, the mechanism of organophosphorus binding and the design of effective reactivators are reviewed. A curated data collection obtained with information from several sources was included for researchers working on the field. Finally, a molecular dynamics simulation with human AChE indicated that the catalytic pocket volume stabilizes around 600 Å(3), providing additional clues for drug design.

Efforts toward treatments against aging of organophosphorus-inhibited acetylcholinesterase

Aging is a dealkylation reaction of organophosphorus (OP)-inhibited acetylcholinesterase (AChE). Despite many studies to date, aged AChE cannot be reactivated directly by traditional pyridinium oximes. This review summarizes strategies that are potentially valuable in the treatment against aging in OP poisoning. Among them, retardation of aging seeks to lower the rate of aging through the use of AChE effectors. These drugs should be administered before AChE is completely aged. For postaging treatment, realkylation of aged AChE by appropriate alkylators may pave the way for oxime treatment by neutralizing the oxyanion at the active site of aged AChE. The other two strategies, upregulation of AChE expression and introduction of exogenous AChE, cannot resurrect aged AChE but may compensate for lowered active AChE levels by in situ production or external introduction of active AChE. Upregulation of AChE expression can be triggered by some peptides. Sources of exogenous AChE can be whole blood or purified AChE, either from human or nonhuman species.

Top-down venomics of the East African green mamba, Dendroaspis angusticeps, and the black mamba, Dendroaspis polylepis, highlight the complexity of their toxin arsenals

We report the characterization, by combination of high-resolution on-line molecular mass and disulfide bond profiling and top-down MS/MS analysis, of the venom proteomes of two congeneric African snake species of medical importance, Dendroaspis angusticeps (green mamba) and D. polylepis (black mamba). Each of these mamba venoms comprised more than two-hundred polypeptides belonging to just a few toxin families. Both venom proteomes are overwhelmingly composed of post-synaptically-acting short- and long-chain neurotoxins that potently inhibit muscle- and neuronal-type nicotinic acetylcholine receptors; muscarinic cardiotoxins; and dendrotoxins, that block some of the Kv1, n-class of K+ channels. However, the identity of the major proteins and their relative abundances exhibit marked interspecific variation. In addition, the greater resolution of the top-down venomic analytical approach revealed previously undetected protein species, isoforms and proteoforms, including the identification and precise location of modified lysine residues in a number of proteins in both venoms, but particularly in green mamba toxins. This comparative top-down venomic analysis unveiled the untapped complexity of Dendroaspis venoms and lays the foundations for rationalizing the notably different potency of green and black mamba lethal arsenals at locus resolution.

SIGNIFICANCE PARAGRAPH

We report the characterization, by combination of high-resolution on-line molecular mass and disulfide bond profiling and top-down MS/MS analysis, of the venom proteomes of two congeneric African snake species of medical importance, Dendroaspis angusticeps (green mamba) and D. polylepis (black mamba). Each of these mamba venoms comprised more than two-hundred polypeptides belonging to just a few toxin families. Both venom proteomes are overwhelmingly composed of post-synaptically-acting short- and long-chain neurotoxins that potently inhibit muscle- and neuronal-type nicotinic acetylcholine receptors; muscarinic cardiotoxins; and dendrotoxins, that block some of the Kv1, n-class of K+ channels. However, the identity of the major proteins and their relative abundances exhibit marked interspecific variation. In addition, the greater resolution of the top-down venomic analytical approach revealed previously undetected protein species, isoforms and proteoforms, including the identification and precise location of modified lysine residues in a number of proteins in both venoms, but particularly in green mamba toxins. This comparative top-down venomic analysis unveiled the untapped complexity of Dendroaspis venoms and lays the foundations for rationalizing the notably different potency of green and black mamba lethal arsenals at locus resolution.

Hopeahainol A binds reversibly at the acetylcholinesterase (AChE) peripheral site and inhibits enzyme activity with a novel higher order concentration dependence

Natural product inhibitors of AChE are of interest both because they offer promise as inexpensive drugs for symptomatic relief in Alzheimer's disease and because they may provide insights into the structural features of the AChE catalytic site. Hopeahainol A is an uncharged polyphenol AChE inhibitor from the stem bark of Hopea hainanensis with a constrained, partially dearomatized bicyclic core. Molecular modeling indicates that hopeahainol A binds at the entrance of the long but narrow AChE active site gorge because it is too bulky to be accommodated within the gorge without severe distortion of the gorge as depicted in AChE crystal structures. We conducted inhibitor competition experiments in which AChE inhibition was measured with hopeahainol A together with either edrophonium (which binds at the base of the gorge) or thioflavin T (which binds to the peripheral or P-site near the gorge mouth). The results agreed with the molecular modeling and indicated that hopeahainol A at lower concentrations (<200 μM) bound only to the P-site, as hopeahainol A and thioflavin T were unable to form a ternary complex with AChE while hopeahainol A and edrophonium did form a ternary complex with essentially no competition between them. Inhibition increased to a striking extent at higher concentrations of hopeahainol A, with plots analogous to classic Dixon plots showing a dependence on hopeahainol A concentrations to the third- or fourth order. The inhibition at higher hopeahainol A concentrations was completely reversed on dilution and blocked by bound edrophonium. We hypothesize that bound hopeahainol A induces conformational changes in the AChE active site that allow binding of additional hopeahainol A molecules, a phenomenon that would be unprecedented for a reversible inhibitor that apparently forms no covalent bonds with AChE.

Design and synthesis of N-substituted-2-hydroxyiminoacetamides and interactions with cholinesterases

Within this study, we designed and synthesized four new oxime compounds of the N-substituted 2-hydroxyiminoacetamide structure and evaluated their interactions with acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). Our aim was to explore the possibility of extending the dual-binding mode of interaction between the enzyme and the inhibitor to a so-called triple-binding mode of interaction through the introduction of an additional binding moiety. N-substituted 2-hydroxyiminoacetamide 1 was prepared via BOP catalyzed amidation of hydroxyiminoacetic acid with 3-azido-1-phenylpropylamine. An azide group enabled us to prepare more elaborate structures 2-4 by the copper-catalyzed azide-alkyne cycloaddition. The new compounds 1-4 differed in their presumed AChE peripheral site binding moiety, which ranged from an azide group to functionalized heterocycles. Molecular docking studies revealed that all three binding moieties are involved in the non-covalent interactions with ChEs for all of the four compounds, albeit not always in the complete accordance with the proposed hypothesis. All of the four compounds reversibly inhibited the ChEs with their inhibition potency increasing in the same order for both enzymes (1 < 2 < 4 < 3). A higher preference for binding to BChE (K_I from 0.30 μmol/L to 130 μmol/L) over AChE (K_I from 50 μmol/L to 1200 μmol/L) was observed for all of the compounds. Compounds were screened for reactivation of cyclosarin-, sarin- and VX-inhibited AChE and BChE.

The impact of crystallization conditions on structure-based drug design: A case study on the methylene blue/acetylcholinesterase complex

Structure-based drug design utilizes apoprotein or complex structures retrieved from the PDB. >57% of crystallographic PDB entries were obtained with polyethylene glycols (PEGs) as precipitant and/or as cryoprotectant, but <6% of these report presence of individual ethyleneglycol oligomers. We report a case in which ethyleneglycol oligomers' presence in a crystal structure markedly affected the bound ligand's position. Specifically, we compared the positions of methylene blue and decamethonium in acetylcholinesterase complexes obtained using isomorphous crystals precipitated with PEG200 or ammonium sulfate. The ligands' positions within the active-site gorge in complexes obtained using PEG200 are influenced by presence of ethyleneglycol oligomers in both cases bound to W84 at the gorge's bottom, preventing interaction of the ligand's proximal quaternary group with its indole. Consequently, both ligands are ∼3.0Å further up the gorge than in complexes obtained using crystals precipitated with ammonium sulfate, in which the quaternary groups make direct π-cation interactions with the indole. These findings have implications for structure-based drug design, since data for ligand-protein complexes with polyethylene glycol as precipitant may not reflect the ligand's position in its absence, and could result in selecting incorrect drug discovery leads. Docking methylene blue into the structure obtained with PEG200, but omitting the ethyleneglycols, yields results agreeing poorly with the crystal structure; excellent agreement is obtained if they are included. Many proteins display features in which precipitants might lodge. It will be important to investigate presence of precipitants in published crystal structures, and whether it has resulted in misinterpreting electron density maps, adversely affecting drug design.

The reactivation of tabun-inhibited mutant AChE with Ortho-7: steered molecular dynamics and quantum chemical studies

A highly toxic nerve agent, tabun, can inhibit acetylcholinesterase (AChE) at cholinergic sites, which leads to serious cardiovascular complications, respiratory compromise and death. We have examined the structural features of the tabun-conjugated AChE complex with an oxime reactivator, Ortho-7, to provide a strategy for designing new and efficient reactivators. Mutation of mAChE within the choline binding site by Y337A and F338A and its interaction with Ortho-7 has been investigated using steered molecular dynamics (SMD) and quantum chemical methods. The overall study shows that after mutagenesis (Y337A), the reactivator can approach more freely towards the phosphorylated active site of serine without any significant steric hindrance in the presence of tabun compared to the wild type and double mutant. Furthermore, the poor binding of Ortho-7 with the peripheral residues of mAChE in the case of the single mutant compared to that of the wild-type and double mutant (Y337A/F338A) can contribute to better efficacy in the former case. Ortho-7 has formed a greater number of hydrogen bonds with the active site surrounding residues His447 and Phe295 in the case of the single mutant (Y337A), and that stabilizes the drug molecule for an effective reactivation process. The DFT M05-2X/6-31+G(d) level of theory shows that the binding energy of Ortho-7 with the single mutant (Y337A) is energetically more preferred (-19.8 kcal mol(-1)) than the wild-type (-8.1 kcal mol(-1)) and double mutant (Y337A/F338A) (-16.0 kcal mol(-1)). The study reveals that both the orientation of the oxime reactivator for nucleophilic attack and the stabilization of the reactivator at the active site would be crucial for the design of an efficient reactivator.

Thyroglobulin From Molecular and Cellular Biology to Clinical Endocrinology

Thyroglobulin (Tg) is a vertebrate secretory protein synthesized in the thyrocyte endoplasmic reticulum (ER), where it acquires N-linked glycosylation and conformational maturation (including formation of many disulfide bonds), leading to homodimerization. Its primary functions include iodide storage and thyroid hormonogenesis. Tg consists largely of repeating domains, and many tyrosyl residues in these domains become iodinated to form monoiodo- and diiodotyrosine, whereas only a small portion of Tg structure is dedicated to hormone formation. Interestingly, evolutionary ancestors, dependent upon thyroid hormone for development, synthesize thyroid hormones without the complete Tg protein architecture. Nevertheless, in all vertebrates, Tg follows a strict pattern of region I, II-III, and the cholinesterase-like (ChEL) domain. In vertebrates, Tg first undergoes intracellular transport through the secretory pathway, which requires the assistance of thyrocyte ER chaperones and oxidoreductases, as well as coordination of distinct regions of Tg, to achieve a native conformation. Curiously, regions II-III and ChEL behave as fully independent folding units that could function as successful secretory proteins by themselves. However, the large Tg region I (bearing the primary T4-forming site) is incompetent by itself for intracellular transport, requiring the downstream regions II-III and ChEL to complete its folding. A combination of nonsense mutations, frameshift mutations, splice site mutations, and missense mutations in Tg occurs spontaneously to cause congenital hypothyroidism and thyroidal ER stress. These Tg mutants are unable to achieve a native conformation within the ER, interfering with the efficiency of Tg maturation and export to the thyroid follicle lumen for iodide storage and hormonogenesis.

Rapid Estimation of Solvation Energy for Simulations of Protein-Protein Association

We have formulated the Energy by Linear Superposition of Corrections Approximation (ELSCA) for estimating the electrostatic and apolar solvation energy of bringing two proteins into close proximity or into contact as defined by the linearized Poisson-Boltzmann model and a linear function of the solvent-accessible surface area. ELSCA utilizes potentials of mean force between atom types found in the AMBER ff99 force field, a uniform distance-dependent dielectric, and a potential that mimics the change in solvent accessible surface area for bringing two solvated spheres into contact. ELSCA was trained by a linear least-squares fit on more than 39 000 putative complexes, each formed from pairs of nonhomologous proteins with a range of shapes, sizes, and charges. The training set was also designed to capture various stages of complex formation. ELSCA was tested against over 8000 non-native complexes of 45 enzyme/inhibitor, antibody/antigen, and other systems that are known to form complexes and gives an overall correlation of 0.962 with PBSA-derived energies for these complexes. The predictions have a slope of 0.89 on the actual values with a bias of 11.1 kcal/mol. When applied to native complexes of these 45 protein systems, ELSCA reproduces PBSA results with a correlation of 0.787, a slope of 1.13, and a bias of 13.0 kcal/mol. We report parameters for ELSCA in the context of the AMBER ff99 parameter set. Our model is most useful in macromolecular docking and protein association simulations, where large portions of each molecule may be considered rigid.

Ortho-7 bound to the active-site gorge of free and OP-conjugated acetylcholinesterase: cation-π interactions

Despite the immense importance of cation-π interactions prevailing in bispyridinium drug acetylcholinesterase (AChE) complexes, a precise description of cation-π interactions at molecular level has remained elusive. Here, we consider a bispyridinium drug, namely, ortho-7 in three different structures of AChE, with and without complexation with organophosphorus (OP) compounds for detailed investigation using all atom molecular dynamics simulation. By quantum mechanical calculations, Y72, W86, Y124, W286, Y337, and Y341 aromatic residues of the enzyme are investigated for possible cation-π interactions with ortho-7. The cation-π interactions in each of the protein-drug complexes are studied using distance, angle, a suitable functional form of them, and electrostatic criteria. The variation of cation-π functional is remarkably consistent with that of the Columbic variation. It is clearly observed that cation-π interactions for some of the residues in the catalytic active site (CAS) and peripheral anionic site (PAS) of the enzyme are either enhanced or reduced based on the nature of OP conjugation (i.e., nerve gas, tabun or pesticide, fenamiphos) when compared with the OP-free enzyme. The strength of cation-π interaction is strongly dependent on the type OP conjugation. The effect of conjugation at CAS is also seen to influence the cation-π interaction at the PAS region. The variation of cation-π interactions on the type of conjugating OP compounds might be suggestive of a reason as to why wide spectrum drug against any OP poisoning is yet to arrive in the market.

Crystal structure of snake venom acetylcholinesterase in complex with inhibitory antibody fragment Fab410 bound at the peripheral site: evidence for open and closed states of a back door channel

The acetylcholinesterase found in the venom of Bungarus fasciatus (BfAChE) is produced as a soluble, non-amphiphilic monomer with a canonical catalytic domain but a distinct C terminus compared with the other vertebrate enzymes. Moreover, the peripheral anionic site of BfAChE, a surface site located at the active site gorge entrance, bears two substitutions altering sensitivity to cationic inhibitors. Antibody Elec410, generated against Electrophorus electricus acetylcholinesterase (EeAChE), inhibits EeAChE and BfAChE by binding to their peripheral sites. However, both complexes retain significant residual catalytic activity, suggesting incomplete gorge occlusion by bound antibody and/or high frequency back door opening. To explore a novel acetylcholinesterase species, ascertain the molecular bases of inhibition by Elec410, and document the determinants and mechanisms for back door opening, we solved a 2.7-Å resolution crystal structure of natural BfAChE in complex with antibody fragment Fab410. Crystalline BfAChE forms the canonical dimer found in all acetylcholinesterase structures. Equally represented open and closed states of a back door channel, associated with alternate positions of a tyrosine phenol ring at the active site base, coexist in each subunit. At the BfAChE molecular surface, Fab410 is seated on the long Ω-loop between two N-glycan chains and partially occludes the gorge entrance, a position that fully reflects the available mutagenesis and biochemical data. Experimentally based flexible molecular docking supports a similar Fab410 binding mode onto the EeAChE antigen. These data document the molecular and dynamic peculiarities of BfAChE with high frequency back door opening, and the mode of action of Elec410 as one of the largest peptidic inhibitors targeting the acetylcholinesterase peripheral site.

Structural differences among subfamilies of EF-hand proteins--a view from the pseudo two-fold symmetry axis

We have analyzed the conformations of EF-lobes, adjacent pairs of EF-hand domains, in a coordinate system based on the approximate two-fold (z) axis that relates the two EF-hands. Two parameters - dE(ø), the azimuthal angle between the y-axis and the projection of the offset vector to helix E onto the yz-plane, and δdF(ø), the difference angle between the two helices (F1 and F2) of odd and even domains--characterize the openness of a single EF-hand domain and of an EF-lobe, respectively. We describe and compare values of dE(ø) and of δdF(ø) for EF-hand proteins of five subfamilies--CTER, CPV, S100, PARV, CALP--in calci- and apo- forms, with and without bound target proteins. Each subfamily has characteristic changes associated with binding calcium and/or target proteins.

Predicting affinity- and specificity-enhancing mutations at protein-protein interfaces

Manipulations of PPIs (protein-protein interactions) are important for many biological applications such as synthetic biology and drug design. Combinatorial methods have been traditionally used for such manipulations, failing, however, to explain the effects achieved. We developed a computational method for prediction of changes in free energy of binding due to mutation that bring about deeper understanding of the molecular forces underlying binding interactions. Our method could be used for computational scanning of binding interfaces and subsequent analysis of the interfacial sequence optimality. The computational method was validated in two biological systems. Computational saturated mutagenesis of a high-affinity complex between an enzyme AChE (acetylcholinesterase) and a snake toxin Fas (fasciculin) revealed the optimal nature of this interface with only a few predicted affinity-enhancing mutations. Binding measurements confirmed high optimality of this interface and identified a few mutations that could further improve interaction fitness. Computational interface scanning of a medium-affinity complex between TIMP-2 (tissue inhibitor of metalloproteinases-2) and MMP (matrix metalloproteinase) 14 revealed a non-optimal nature of the binding interface with multiple mutations predicted to stabilize the complex. Experimental results corroborated our computational predictions, identifying a large number of mutations that improve the binding affinity for this interaction and some mutations that enhance binding specificity. Overall, our computational protocol greatly facilitates the discovery of affinity- and specificity-enhancing mutations and thus could be applied for design of potent and highly specific inhibitors of any PPI.

Mapping of the binding landscape for a picomolar protein-protein complex through computation and experiment

Our understanding of protein evolution would greatly benefit from mapping of binding landscapes, i.e., changes in protein-protein binding affinity due to all single mutations. However, experimental generation of such landscapes is a tedious task due to a large number of possible mutations. Here, we use a simple computational protocol to map the binding landscape for two homologous high-affinity complexes, involving a snake toxin fasciculin and acetylcholinesterase from two different species. To verify our computational predictions, we experimentally measure binding between 25 Fas mutants and the 2 enzymes. Both computational and experimental results demonstrate that the Fas sequence is close to the optimum when interacting with its targets, yet a few mutations could further improve Kd, kon, and koff. Our computational predictions agree well with experimental results and generate distributions similar to those observed in other high-affinity PPIs, demonstrating the potential of simple computational protocols in capturing realistic binding landscapes.

Molecular characterization of monoclonal antibodies that inhibit acetylcholinesterase by targeting the peripheral site and backdoor region

The inhibition properties and target sites of monoclonal antibodies (mAbs) Elec403, Elec408 and Elec410, generated against Electrophorus electricus acetylcholinesterase (AChE), have been defined previously using biochemical and mutagenesis approaches. Elec403 and Elec410, which bind competitively with each other and with the peptidic toxin inhibitor fasciculin, are directed toward distinctive albeit overlapping epitopes located at the AChE peripheral anionic site, which surrounds the entrance of the active site gorge. Elec408, which is not competitive with the other two mAbs nor fasciculin, targets a second epitope located in the backdoor region, distant from the gorge entrance. To characterize the molecular determinants dictating their binding site specificity, we cloned and sequenced the mAbs; generated antigen-binding fragments (Fab) retaining the parental inhibition properties; and explored their structure-function relationships using complementary x-ray crystallography, homology modeling and flexible docking approaches. Hypermutation of one Elec403 complementarity-determining region suggests occurrence of antigen-driven selection towards recognition of the AChE peripheral site. Comparative analysis of the 1.9Å-resolution structure of Fab408 and of theoretical models of its Fab403 and Fab410 congeners evidences distinctive surface topographies and anisotropic repartitions of charges, consistent with their respective target sites and inhibition properties. Finally, a validated, data-driven docking model of the Fab403-AChE complex suggests a mode of binding at the PAS that fully correlates with the functional data. This comprehensive study documents the molecular peculiarities of Fab403 and Fab410, as the largest peptidic inhibitors directed towards the peripheral site, and those of Fab408, as the first inhibitor directed toward the backdoor region of an AChE and a unique template for the design of new, specific modulators of AChE catalysis.

The natural product dihydrotanshinone I provides a prototype for uncharged inhibitors that bind specifically to the acetylcholinesterase peripheral site with nanomolar affinity

Cholinergic synaptic transmission often requires extremely rapid hydrolysis of acetylcholine by acetylcholinesterase (AChE). AChE is inactivated by organophosphates (OPs) in chemical warfare nerve agents. The resulting accumulation of acetylcholine disrupts cholinergic synaptic transmission and can lead to death. A potential long-term strategy for preventing AChE inactivation by OPs is based on evidence that OPs must pass through a peripheral site or P-site near the mouth of the AChE active site gorge before reacting with a catalytic serine in an acylation site or A-site at the base of the gorge. An ultimate goal of this strategy is to design compounds that bind tightly at or near the P-site and exclude OPs from the active site while interfering minimally with the passage of acetylcholine. However, to target the AChE P-site with ligands and potential drugs that selectively restrict access, much more information must be gathered about the structure-activity relationships of ligands that bind specifically to the P-site. We apply here an inhibitor competition assay that can correctly determine whether an AChE inhibitor binds to the P-site, the A-site, or both sites. We have used this assay to examine three uncharged, natural product inhibitors of AChE, including aflatoxin B1, dihydrotanshinone I, and territrem B. The first two of these inhibitors are predicted by the competition assay to bind selectively to the P-site, while territrem B is predicted to span both the P- and A-sites. These predictions have recently been confirmed by X-ray crystallography. Dihydrotanshinone I, with an observed binding constant (KI) of 750 nM, provides a good lead compound for the development of high-affinity, uncharged inhibitors with specificity for the P-site.