Induced-fit or preexisting equilibrium dynamics? Lessons from protein crystallography and MD simulations on acetylcholinesterase and implications for structure-based drug design

Novichok Nerve Agents as Inhibitors of Acetylcholinesterase-In Silico Study of Their Non-Covalent Binding Affinity

In silico studies were performed to assess the binding affinity of selected organophosphorus compounds toward the acetylcholinesterase enzyme (AChE). Quantum mechanical calculations, molecular docking, and molecular dynamics (MD) with molecular mechanics Generalized-Born surface area (MM/GBSA) were applied to assess quantitatively differences between the binding energies of acetylcholine (ACh; the natural agonist of AChE) and neurotoxic, synthetic correlatives (so-called "Novichoks", and selected compounds from the G- and V-series). Several additional quantitative descriptors like root-mean-square fluctuation (RMSF) and the solvent accessible surface area (SASA) were briefly discussed to give-to the best of our knowledge-the first quantitative in silico description of AChE-Novichok non-covalent binding process and thus facilitate the search for an efficient and effective treatment for Novichok intoxication and in a broader sense-intoxication with other warfare nerve agents as well.

Chalcone Scaffolds Exhibiting Acetylcholinesterase Enzyme Inhibition: Mechanistic and Computational Investigations

This study was aimed to perform the mechanistic investigations of chalcone scaffold as inhibitors of acetylcholinesterase (AChE) enzyme using molecular docking and molecular dynamics simulation tools. Basic chalcones (C1–C5) were synthesized and their in vitro AChE inhibition was tested. Binding interactions were studied using AutoDock and Surflex-Dock programs, whereas the molecular dynamics simulation studies were performed to check the stability of the ligand–protein complex. Good AChE inhibition (IC₅₀ = 22 ± 2.8 to 37.6 ± 0.75 μM) in correlation with the in silico results (binding energies = −8.55 to −8.14 Kcal/mol) were obtained. The mechanistic studies showed that all of the functionalities present in the chalcone scaffold were involved in binding with the amino acid residues at the binding site through hydrogen bonding, π–π, π–cation, π–sigma, and hydrophobic interactions. Molecular dynamics simulation studies showed the formation of stable complex between the AChE enzyme and C4 ligand.

Computational evidence of new putative allosteric sites in the acetylcholinesterase receptor

Acetylcholinesterase (AChE), with a rigid structure and buried active site at the end of a deep narrow gorge, is interesting enough to solve the paradox between high catalytic activity and unavailability of the active site in treatment of Alzheimer's disease (AD). In this way, the blind docking process is performed on an ensemble of AChE structures created with molecular dynamics (MD) simulations to survey the whole space of AChE to find multiple access pathways to the active site and ranking them based on their affinity scores. Our results show that there are other allosteric binding sites in the protein structure whose inhibition, can affect protein function by disrupting the release of the Acetylcholine (AC) degradation products. In this study, inhibitory activities of Hybride14 and two natural compounds (Papaverine and Palmatine) were evaluated for all possible allosteric sites via docking method. The results confirmed the non-competitive inhibition mechanism. The best binding mode for these inhibitors and efficacy of hydrogen bonds and hydrophobic interactions on inhibitory activities of ligands were also disclosed. Furthermore, our studies provide significant molecular insight for AChE inhibition that could aid in the development of new drugs for AD's treatment.

In silico Exploration of Bioactive Phytochemicals Against Neurodegenerative Diseases Via Inhibition of Cholinesterases

Neurodegenerative disorders are estimated to become the second leading cause of death worldwide by 2040. Despite the widespread use of diverse allopathic drugs, these brain-associated disorders can only be partially addressed and long term treatment is often linked with dependency and other unwanted side effects. Nature, believed to be an arsenal of remedies for any illness, presents an interesting avenue for the development of novel neuroprotective agents. Interestingly, inhibition of cholinesterases, involved in the breakdown of acetylcholine in the synaptic cleft, has been proposed to be neuroprotective. This review therefore aims to provide additional insight via docking studies of previously studied compounds that have shown potent activity against acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) in vitro. Indeed, the determination of potent plant-based ligands for this purpose through in silico methods enables the elimination of lengthy and costly traditional methods of drug discovery. Herein, a literature search was conducted to identify active phytochemicals which are cholinesterase inhibitors. Following which in silico docking methods were applied to obtain docking scores. Compound structures were extracted from online ZINC database and optimized using AM1 implemented in gaussian09 software. Noteworthy ligands against AChE highlighted in this study include: 19,20-dihydroervahanine A and 19, 20-dihydrotabernamine. Regarding BChE inhibition, the best ligands were found to be 8-Clavandurylkaempferol, Na-methylepipachysamine D; ebeiedinone; and dictyophlebine. Thus, ligand optimization between such phytochemicals and cholinesterases coupled with in vitro, in vivo studies and randomized clinical trials can lead to the development of novel drugs against neurodegenerative disorders.

Impact of Sucrose as Osmolyte on Molecular Dynamics of Mouse Acetylcholinesterase

The enzyme model, mouse acetylcholinesterase, which exhibits its active site at the bottom of a narrow gorge, was investigated in the presence of different concentrations of sucrose to shed light on the protein and water dynamics in cholinesterases. The study was conducted by incoherent neutron scattering, giving access to molecular dynamics within the time scale of sub-nano to nanoseconds, in comparison with molecular dynamics simulations. With increasing sucrose concentration, we found non-linear effects, e.g., first a decrease in the dynamics at 5 wt% followed by a gain at 10 wt% sucrose. Direct comparisons with simulations permitted us to understand the following findings: at 5 wt%, sugar molecules interact with the protein surface through water molecules and damp the motions to reduce the overall protein mobility, although the motions inside the gorge are enhanced due to water depletion. When going to 10 wt% of sucrose, some water molecules at the protein surface are replaced by sugar molecules. By penetrating the protein surface, they disrupt some of the intra-protein contacts, and induce new ones, creating new pathways for correlated motions, and therefore, increasing the dynamics. This exhaustive study allowed for an explanation of the detail interactions leading to the observed non-linear behavior.

Spontaneous Hinge-Bending Motions of Angiotensin I Converting Enzyme: Role in Activation and Inhibition

The inhibition of human angiotensin I converting enzyme (ACE) has been regarded as a promising approach for the treatment of hypertension. Despite research attempts over many years, our understanding the mechanisms of activation and inhibition of ACE is still far from complete. Here, we present results of all atom molecular dynamics simulations of ACE with and without ligands. Two types of inhibitors, competitive and mixed non-competitive, were used to model the ligand bound forms. In the absence of a ligand the simulation showed spontaneous large hinge-bending motions of multiple conversions between the closed and open states of ACE, while the ligand bound forms were stable in the closed state. Our simulation results imply that the equilibrium between pre-existing backbone conformations shifts in the presence of a ligand. The hinge-bending motion of ACE is considered as an essential to the enzyme function. A mechanistic model of activation and the inhibition may provide valuable information for novel inhibitors of ACE.

Gorge Motions of Acetylcholinesterase Revealed by Microsecond Molecular Dynamics Simulations

Acetylcholinesterase, with a deep, narrow active-site gorge, attracts enormous interest due to its particularly high catalytic efficiency and its inhibitors used for treatment of Alzheimer’s disease. To facilitate the massive pass-through of the substrate and inhibitors, “breathing” motions to modulate the size of the gorge are an important prerequisite. However, the molecular mechanism that governs such motions is not well explored. Here, to systematically investigate intrinsic motions of the enzyme, we performed microsecond molecular dynamics simulations on the monomer and dimer of Torpedo californica acetylcholinesterase (TcAChE) as well as the complex of TcAChE bound with the drug E2020. It has been revealed that protein-ligand interactions and dimerization both keep the gorge in bulk, and opening events of the gorge increase dramatically compared to the monomer. Dynamics of three subdomains, S3, S4 and the Ω-loop, are tightly associated with variations of the gorge size while the dynamics can be changed by ligand binding or protein dimerization. Moreover, high correlations among these subdomains provide a basis for remote residues allosterically modulating the gorge motions. These observations are propitious to expand our understanding of protein structure and function as well as providing clues for performing structure-based drug design.

Single Molecule Hydrodynamic Separation Allows Sensitive and Quantitative Analysis of DNA Conformation and Binding Interactions in Free Solution

Limited tools exist that are capable of monitoring nucleic acid conformations, fluctuations, and distributions in free solution environments. Single molecule free solution hydrodynamic separation enables the unique ability to quantitatively analyze nucleic acid biophysics in free solution. Single molecule fluorescent burst data and separation chromatograms can give layered insight into global DNA conformation, binding interactions, and molecular distributions. First, we show that global conformation of individual DNA molecules can be directly visualized by examining single molecule fluorescent burst shapes and that DNA exists in a dynamic equilibrium of fluctuating conformations as it is driven by Poiseuille flow through micron-sized channels. We then show that this dynamic equilibrium of DNA conformations is reflected as shifts in hydrodynamic mobility that can be perturbed using salt and ionic strength to affect packing density. Next, we demonstrate that these shifts in hydrodynamic mobility can be used to investigate hybridization thermodynamics and binding interactions. We distinguish and classify multiple interactions within a single sample, and demonstrate quantification amidst large concentration differences for the detection of rare species. Finally, we demonstrate that these differences can resolve perfect complement, 2 bp mismatched, and 3 bp mismatched sequences. Such a system can be used to garner diverse information about DNA conformation and structure, and potentially be extended to other molecules and mixed-species interactions, such as between nucleic acids and proteins or synthetic polymers.

Anti-choline esterase activity of ceramides from the Red Sea marine sponge Mycale euplectellioides

New ceramides were isolated, and structurally elucidated, from the marine sponge Mycale euplectellioides. The isolated compounds showed moderate anti-choline esterase activity in vitro and tight binding to AChE as confirmed through a docking study.

Sampling of conformational ensemble for virtual screening using molecular dynamics simulations and normal mode analysis

Aim: Molecular dynamics simulations and normal mode analysis are well-established approaches to generate receptor conformational ensembles (RCEs) for ligand docking and virtual screening. Here, we report new fast molecular dynamics-based and normal mode analysis-based protocols combined with conformational pocket classifications to efficiently generate RCEs. Materials & Methods: We assessed our protocols on two well-characterized protein targets showing local active site flexibility, dihydrofolate reductase and large collective movements, CDK2. The performance of the RCEs was validated by distinguishing known ligands of dihydrofolate reductase and CDK2 among a dataset of diverse chemical decoys. Results & discussion: Our results show that different simulation protocols can be efficient for generation of RCEs depending on different kind of protein flexibility.[Formula: see text]

A wrench in the works of human acetylcholinesterase: soman induced conformational changes revealed by molecular dynamics simulations

Irreversible inactivation of human acetylcholinesterase (hAChE) by organophosphorous pesticides (OPs) and chemical weapon agents (CWA) has severe morbidity and mortality consequences. We present data from quantum mechanics/molecular mechanics (QM/MM) and 80 classical molecular dynamics (MD) simulations of the apo and soman-adducted forms of hAChE to investigate the effects on the dynamics and protein structure when the catalytic Serine 203 is phosphonylated. We find that the soman phosphonylation of the active site Ser203 follows a water assisted addition-elimination mechanism with the elimination of the fluoride ion being the highest energy barrier at 6.5 kcal/mole. We observe soman-dependent changes in backbone and sidechain motions compared to the apo form of the protein. These alterations restrict the soman-adducted hAChE to a structural state that is primed for the soman adduct to be cleaved and removed from the active site. The altered motions and resulting structures provide alternative pathways into and out of the hAChE active site. In the soman-adducted protein both side and back door pathways are viable for soman adduct access. Correlation analysis of the apo and soman adducted MD trajectories shows that the correlation of gorge entrance and back door motion is disrupted when hAChE is adducted. This supports the hypothesis that substrate and product can use two different pathways as entry and exit sites in the apo form of the protein. These alternative pathways have important implications for the rational design of medical countermeasures.

Correlation of the dynamics of native human acetylcholinesterase and its inhibited huperzine A counterpart from sub-picoseconds to nanoseconds

It is a long debated question whether catalytic activities of enzymes, which lie on the millisecond timescale, are possibly already reflected in variations in atomic thermal fluctuations on the pico- to nanosecond timescale. To shed light on this puzzle, the enzyme human acetylcholinesterase in its wild-type form and complexed with the inhibitor huperzine A were investigated by various neutron scattering techniques and molecular dynamics simulations. Previous results on elastic neutron scattering at various timescales and simulations suggest that dynamical processes are not affected on average by the presence of the ligand within the considered time ranges between 10 ps and 1 ns. In the work presented here, the focus was laid on quasi-elastic (QENS) and inelastic neutron scattering (INS). These techniques give access to different kinds of individual diffusive motions and to the density of states of collective motions at the sub-picoseconds timescale. Hence, they permit going beyond the first approach of looking at mean square displacements. For both samples, the autocorrelation function was well described by a stretched-exponential function indicating a linkage between the timescales of fast and slow functional relaxation dynamics. The findings of the QENS and INS investigation are discussed in relation to the results of our earlier elastic incoherent neutron scattering and molecular dynamics simulations.

Evaluation of binding competition and group epitopes of acetylcholinesterase inhibitors by STD NMR, Tr-NOESY, DOSY and molecular docking: an old approach but new findings

We suggest that the small inhibitors of AChE like coumarin and gallic acid's structure-based scaffolds may help in controlling Alzheimer's disease.

Identification of potential herbal inhibitor of acetylcholinesterase associated Alzheimer's disorders using molecular docking and molecular dynamics simulation

Cholinesterase inhibitors (ChE-Is) are the standard for the therapy of AD associated disorders and are the only class of approved drugs by the Food and Drug Administration (FDA). Additionally, acetylcholinesterase (AChE) is the target for many Alzheimer's dementia drugs which block the function of AChE but have some side effects. Therefore, in this paper, an attempt was made to elucidate cholinesterase inhibition potential of secondary metabolite from Cannabis plant which has negligible or no side effect. Molecular docking of 500 herbal compounds, against AChE, was performed using Autodock 4.2 as per the standard protocols. Molecular dynamics simulations have also been carried out to check stability of binding complex in water for 1000 ps. Our molecular docking and simulation have predicted high binding affinity of secondary metabolite (C28H34N2O6) to AChE. Further, molecular dynamics simulations for 1000 ps suggest that ligand interaction with the residues Asp72, Tyr70-121-334, and Phe288 of AChE, all of which fall under active site/subsite or binding pocket, might be critical for the inhibitory activity of AChE. This approach might be helpful to understand the selectivity of the given drug molecule in the treatment of Alzheimer's disease. The study provides evidence for consideration of C28H34N2O6 as a valuable small ligand molecule in treatment and prevention of AD associated disorders and further in vitro and in vivo investigations may prove its therapeutic potential.

Molecular dynamics simulations indicate an induced-fit mechanism for LSD1/CoREST-H3-histone molecular recognition

BACKGROUND

Lysine Specific Demethylase (LSD1 or KDM1A) in complex with its co-repressor protein CoREST catalyzes the demethylation of the H3 histone N-terminal tail and is currently one of the most promising epigenetic targets for drug discovery against cancer and neurodegenerative diseases. Models of non-covalent binding, such as lock and key, induced-fit, and conformational selection could help explaining the molecular mechanism of LSD1/CoREST-H3-histone association, thus guiding drug discovery and design efforts. Here, we quantify the extent to which LSD1/CoREST substrate binding is consistent with these hypothetical models using LSD1/CoREST conformational ensembles obtained through extensive explicit solvent molecular dynamics (MD) simulations.

RESULTS

We find that an induced-fit model is the most representative of LSD1/CoREST-H3-histone non-covalent binding and accounts for the local conformational changes occurring in the H3-histone binding site. We also show that conformational selection - despite in principle not ruled out by this finding - is minimal, and only relevant when global properties are considered, e.g. the nanoscale motion of the LSD1/CoREST clamp.

CONCLUSION

The induced-fit mechanism revealed by our MD simulation study will aid the inclusion of protein dynamics for the discovery and design of LSD1 inhibitors targeting the H3-histone binding region. On a general basis, our study indicates the importance of using multiple metrics or selection schemes when testing alternative hypothetical mechanistic models of non-covalent binding.

Time-dependent kinetic complexities in cholinesterase-catalyzed reactions

Cholinesterases (ChEs) display a hysteretic behavior with certain substrates and inhibitors. Kinetic cooperativity in hysteresis of ChE-catalyzed reactions is characterized by a lag or burst phase in the approach to steady state. With some substrates damped oscillations are shown to superimpose on hysteretic lags. These time dependent peculiarities are observed for both butyrylcholinesterase and acetylcholinesterase from different sources. Hysteresis in ChE-catalyzed reactions can be interpreted in terms of slow transitions between two enzyme conformers E and E'. Substrate can bind to E and/or E', both Michaelian complexes ES and Ε'S can be catalytically competent, or only one of them can make products. The formal reaction pathway depends on both the chemical structure of the substrate and the type of enzyme. In particular, damped oscillations develop when substrate exists in different, slowly interconvertible, conformational, and/or micellar forms, of which only the minor form is capable of binding and reacting with the enzyme. Biphasic pseudo-first-order progressive inhibition of ChEs by certain carbamates and organophosphates also fits with a slow equilibrium between two reactive enzyme forms. Hysteresis can be modulated by medium parameters (pH, chaotropic and kosmotropic salts, organic solvents, temperature, osmotic pressure, and hydrostatic pressure). These studies showed that water structure plays a role in hysteretic behavior of ChEs. Attempts to provide a molecular mechanism for ChE hysteresis from mutagenesis studies or crystallographic studies failed so far. In fact, several lines of evidence suggest that hysteresis is controlled by the conformation of His438, a key residue in the catalytic triad of cholinesterases. Induction time may depend on the probability of His438 to adopt the operative conformation in the catalytic triad. The functional significance of ChE hysteresis is puzzling. However, the accepted view that proteins are in equilibrium between preexisting functional and non-functional conformers, and that binding of a ligand to the functional form shifts equilibrium towards the functional conformation, suggests that slow equilibrium between two conformational states of these enzymes may have a regulatory function in damping out the response to certain ligands and irreversible inhibitors. This is particularly true for immobilized (membrane bound) enzymes where the local substrate and/or inhibitor concentrations depend on influx in crowded organellar systems, e.g. cholinergic synaptic clefts. Therefore, physiological or toxicological relevance of the hysteretic behavior and damped oscillations in ChE-catalyzed reactions and inhibition cannot be ruled out.

Implicit ligand theory: rigorous binding free energies and thermodynamic expectations from molecular docking

A rigorous formalism for estimating noncovalent binding free energies and thermodynamic expectations from calculations in which receptor configurations are sampled independently from the ligand is derived. Due to this separation, receptor configurations only need to be sampled once, facilitating the use of binding free energy calculations in virtual screening. Demonstrative calculations on a host-guest system yield good agreement with previous free energy calculations and isothermal titration calorimetry measurements. Implicit ligand theory provides guidance on how to improve existing molecular docking algorithms and insight into the concepts of induced fit and conformational selection in noncovalent macromolecular recognition.

Molecular dynamics simulation in virus research

Virus replication in the host proceeds by chains of interactions between viral and host proteins. The interactions are deeply influenced by host immune molecules and anti-viral compounds, as well as by mutations in viral proteins. To understand how these interactions proceed mechanically and how they are influenced by mutations, one needs to know the structures and dynamics of the proteins. Molecular dynamics (MD) simulation is a powerful computational method for delineating motions of proteins at an atomic-scale via theoretical and empirical principles in physical chemistry. Recent advances in the hardware and software for biomolecular simulation have rapidly improved the precision and performance of this technique. Consequently, MD simulation is quickly extending the range of applications in biology, helping to reveal unique features of protein structures that would be hard to obtain by experimental methods alone. In this review, we summarize the recent advances in MD simulations in the study of virus–host interactions and evolution, and present future perspectives on this technique.

The discovery of potential acetylcholinesterase inhibitors: a combination of pharmacophore modeling, virtual screening, and molecular docking studies

BACKGROUND

Alzheimer's disease (AD) is the most common cause of dementia characterized by progressive cognitive impairment in the elderly people. The most dramatic abnormalities are those of the cholinergic system. Acetylcholinesterase (AChE) plays a key role in the regulation of the cholinergic system, and hence, inhibition of AChE has emerged as one of the most promising strategies for the treatment of AD.

METHODS

In this study, we suggest a workflow for the identification and prioritization of potential compounds targeted against AChE. In order to elucidate the essential structural features for AChE, three-dimensional pharmacophore models were constructed using Discovery Studio 2.5.5 (DS 2.5.5) program based on a set of known AChE inhibitors.

RESULTS

The best five-features pharmacophore model, which includes one hydrogen bond donor and four hydrophobic features, was generated from a training set of 62 compounds that yielded a correlation coefficient of R = 0.851 and a high prediction of fit values for a set of 26 test molecules with a correlation of R² = 0.830. Our pharmacophore model also has a high Güner-Henry score and enrichment factor. Virtual screening performed on the NCI database obtained new inhibitors which have the potential to inhibit AChE and to protect neurons from Aβ toxicity. The hit compounds were subsequently subjected to molecular docking and evaluated by consensus scoring function, which resulted in 9 compounds with high pharmacophore fit values and predicted biological activity scores. These compounds showed interactions with important residues at the active site.

CONCLUSIONS

The information gained from this study may assist in the discovery of potential AChE inhibitors that are highly selective for its dual binding sites.

Mining electron density for functionally relevant protein polysterism in crystal structures

This review focuses on conceptual and methodological advances in our understanding and characterization of the conformational heterogeneity of proteins. Focusing on X-ray crystallography, we describe how polysterism, the interconversion of pre-existing conformational substates, has traditionally been analyzed by comparing independent crystal structures or multiple chains within a single crystal asymmetric unit. In contrast, recent studies have focused on mining electron density maps to reveal previously ‘hidden’ minor conformational substates. Functional tests of the importance of minor states suggest that evolutionary selection shapes the entire conformational landscape, including uniquely configured conformational substates, the relative distribution of these substates, and the speed at which the protein can interconvert between them. An increased focus on polysterism may shape the way protein structure and function is studied in the coming years.

Protein functional landscapes, dynamics, allostery: a tortuous path towards a universal theoretical framework

Energy landscape theories have provided a common ground for understanding the protein folding problem, which once seemed to be overwhelmingly complicated. At the same time, the native state was found to be an ensemble of interconverting states with frustration playing a more important role compared to the folding problem. The landscape of the folded protein – the native landscape – is glassier than the folding landscape; hence, a general description analogous to the folding theories is difficult to achieve. On the other hand, the native basin phase volume is much smaller, allowing a protein to fully sample its native energy landscape on the biological timescales. Current computational resources may also be used to perform this sampling for smaller proteins, to build a ‘topographical map’ of the native landscape that can be used for subsequent analysis. Several major approaches to representing this topographical map are highlighted in this review, including the construction of kinetic networks, hierarchical trees and free energy surfaces with subsequent structural and kinetic analyses. In this review, we extensively discuss the important question of choosing proper collective coordinates characterizing functional motions. In many cases, the substates on the native energy landscape, which represent different functional states, can be used to obtain variables that are well suited for building free energy surfaces and analyzing the protein's functional dynamics. Normal mode analysis can provide such variables in cases where functional motions are dictated by the molecule's architecture. Principal component analysis is a more expensive way of inferring the essential variables from the protein's motions, one that requires a long molecular dynamics simulation. Finally, the two popular models for the allosteric switching mechanism, ‘preexisting equilibrium’ and ‘induced fit’, are interpreted within the energy landscape paradigm as extreme points of a continuum of transition mechanisms. Some experimental evidence illustrating each of these two models, as well as intermediate mechanisms, is presented and discussed.

Global dynamics of proteins: bridging between structure and function

Biomolecular systems possess unique, structure-encoded dynamic properties that underlie their biological functions. Recent studies indicate that these dynamic properties are determined to a large extent by the topology of native contacts. In recent years, elastic network models used in conjunction with normal mode analyses have proven to be useful for elucidating the collective dynamics intrinsically accessible under native state conditions, including in particular the global modes of motions that are robustly defined by the overall architecture. With increasing availability of structural data for well-studied proteins in different forms (liganded, complexed, or free), there is increasing evidence in support of the correspondence between functional changes in structures observed in experiments and the global motions predicted by these coarse-grained analyses. These observed correlations suggest that computational methods may be advantageously employed for assessing functional changes in structure and allosteric mechanisms intrinsically favored by the native fold.

Acetylcholinesterase: from 3D structure to function

By rapid hydrolysis of the neurotransmitter, acetylcholine, acetylcholinesterase terminates neurotransmission at cholinergic synapses. Acetylcholinesterase is a very fast enzyme, functioning at a rate approaching that of a diffusion-controlled reaction. The powerful toxicity of organophosphate poisons is attributed primarily to their potent inhibition of acetylcholinesterase. Acetylcholinesterase inhibitors are utilized in the treatment of various neurological disorders, and are the principal drugs approved thus far by the FDA for management of Alzheimer's disease. Many organophosphates and carbamates serve as potent insecticides, by selectively inhibiting insect acetylcholinesterase. The determination of the crystal structure of Torpedo californica acetylcholinesterase permitted visualization, for the first time, at atomic resolution, of a binding pocket for acetylcholine. It also allowed identification of the active site of acetylcholinesterase, which, unexpectedly, is located at the bottom of a deep gorge lined largely by aromatic residues. The crystal structure of recombinant human acetylcholinesterase in its apo-state is similar in its overall features to that of the Torpedo enzyme; however, the unique crystal packing reveals a novel peptide sequence which blocks access to the active-site gorge.

Crystallographic snapshots of nonaged and aged conjugates of soman with acetylcholinesterase, and of a ternary complex of the aged conjugate with pralidoxime

Organophosphate compounds (OP) are potent inhibitors of acetylcholinesterases (AChEs) and can cause lethal poisoning in humans. Inhibition of AChEs by the OP soman involves phosphonylation of the catalytic serine, and subsequent dealkylation produces a form known as the "aged" enzyme. The nonaged form can be reactivated to a certain extent by nucleophiles, such as pralidoxime (2-PAM), whereas aged forms of OP-inhibited AChEs are totally resistant to reactivation. Here, we solved the X-ray crystal structures of AChE from Torpedo californica (TcAChE) conjugated with soman before and after aging. The absolute configuration of the soman stereoisomer adduct in the nonaged conjugate is P(S)C(R). A structural reorientation of the catalytic His440 side chain was observed during the aging process. Furthermore, the crystal structure of the ternary complex of the aged conjugate with 2-PAM revealed that the orientation of the oxime function does not permit nucleophilic attack on the phosphorus atom, thus providing a plausible explanation for its failure to reactivate the aged soman/AChE conjugate. Together, these three crystal structures provide an experimental basis for the design of new reactivators.

GROMOS96 43a1 performance in predicting oligosaccharide conformational ensembles within glycoproteins

In previous work [Pol-Fachin, L.; Fernandes, C. L.; Verli, H.; Carbohydr. Res.2009, 344, 491-500], we had demonstrated that GROMOS96 43a1 force field and Löwdin HF/6-31G * *-derived atomic charges, adequately represent a glycoprotein's conformational ensemble in aqueous solutions, taking as the starting geometries NMR-determined structures. Based on such data, the present work intends to evaluate the use of the main solution conformations of isolated disaccharides, to build the carbohydrate moiety of glycoproteins, for which no previous experimental information is available. The observed results suggested that the entire glycoprotein scaffold appears unable to promote major modifications in the conformational behavior of glycosidic linkages. Additionally, when compared to energy contour plots, the results support the use of solution ensembles, to refine vacuum conformations of carbohydrate databases in the assembling of glycoproteins 3D structures. Finally, such approach is applied to build a full glycosylated model for COX-1 and COX-2 enzymes.

Improving structure-based function prediction using molecular dynamics

The number of molecules with solved three-dimensional structure but unknown function is increasing rapidly. Particularly problematic are novel folds with little detectable similarity to molecules of known function. Experimental assays can determine the functions of such molecules, but are time-consuming and expensive. Computational approaches can identify potential functional sites; however, these approaches generally rely on single static structures and do not use information about dynamics. In fact, structural dynamics can enhance function prediction: we coupled molecular dynamics simulations with structure-based function prediction algorithms that identify Ca(2+) binding sites. When applied to 11 challenging proteins, both methods showed substantial improvement in performance, revealing 22 more sites in one case and 12 more in the other, with a modest increase in apparent false positives. Thus, we show that treating molecules as dynamic entities improves the performance of structure-based function prediction methods.

Conformational selection and induced fit mechanism underlie specificity in noncovalent interactions with ubiquitin

Noncovalent binding interactions between proteins are the central physicochemical phenomenon underlying biological signaling and functional control on the molecular level. Here, we perform an extensive structural analysis of a large set of bound and unbound ubiquitin conformers and study the level of residual induced fit after conformational selection in the binding process. We show that the region surrounding the binding site in ubiquitin undergoes conformational changes that are significantly more pronounced compared with the whole molecule on average. We demonstrate that these induced-fit structural adjustments are comparable in magnitude to conformational selection. Our final model of ubiquitin binding blends conformational selection with the subsequent induced fit and provides a quantitative measure of their respective contributions.

Protein flexibility: coordinate uncertainties and interpretation of structural differences

Valid interpretations of conformational movements in protein structures determined by X-ray crystallography require that the movement magnitudes exceed their uncertainty threshold. Here, it is shown that such thresholds can be obtained from the distance difference matrices (DDMs) of 1014 pairs of independently determined structures of bovine ribonuclease A and sperm whale myoglobin, with no explanations provided for reportedly minor coordinate differences. The smallest magnitudes of reportedly functional motions are just above these thresholds. Uncertainty thresholds can provide objective criteria that distinguish between true conformational changes and apparent 'noise', showing that some previous interpretations of protein coordinate changes attributed to external conditions or mutations may be doubtful or erroneous. The use of uncertainty thresholds, DDMs, the newly introduced CDDMs (contact distance difference matrices) and a novel simple rotation algorithm allows a more meaningful classification and description of protein motions, distinguishing between various rigid-fragment motions and nonrigid conformational deformations. It is also shown that half of 75 pairs of identical molecules, each from the same asymmetric crystallographic cell, exhibit coordinate differences that range from just outside the coordinate uncertainty threshold to the full magnitude of large functional movements. Thus, crystallization might often induce protein conformational changes that are comparable to those related to or induced by the protein function.

Consistent improvement of cross-docking results using binding site ensembles generated with elastic network normal modes

The representation of protein flexibility is still a challenge for the state-of-the-art flexible ligand docking protocols. In this article, we use a large and diverse benchmark to prove that is possible to improve consistently the cross-docking performance against a single receptor conformation, using an equilibrium ensemble of binding site conformers. The benchmark contained 28 proteins, and our method predicted the top-ranked near native ligand poses 20% more efficiently than using a single receptor. The multiple conformations were derived from the collective variable space defined by all heavy-atom elastic network normal modes, including backbone and side chains. We have found that the binding site displacements for best positioning of the ligand seem rather independent from the global collective motions of the protein. We also found that the number of binding site conformations needed to represent nonredundant flexibility was < 100. The ensemble of receptor conformations can be generated at our Web site at http://abagyan.scripps.edu/MRC.

Conformational changes at the tetramerization site of erythroid alpha-spectrin upon binding beta-spectrin: a spin label EPR study

We used cysteine scanning, isothermal titration calorimetry (ITC) and spin label EPR methods to study the two regions that flank the partial domain Helix C' of the N-terminal end of alpha-spectrin (residues 14-20 and residues 44-54) in the absence and presence of a model protein of the beta-spectrin C-terminal end. In the absence of beta-spectrin, residues 14-20 and 46-52 were known to be unstructured. The EPR spectral values of the inverse line width (Delta H (-1)) and of the width between the low field peak and the central peak ( aZ) of residues in part of the first unstructured region (residues 17-20) and of most residues in the second unstructured junction region (residues 46-52) changed dramatically upon association with beta-spectrin, suggesting that the two regions undergo a conformational change, becoming more rigid and likely becoming helical. ITC results showed that three of the seven residues in the junction region (residues 46-52) were very important in its association with beta-spectrin, in the following order: L49 > G46 > K48. In general, our results suggest that any mutations that affect the propensity of helical formation in the region spanning residues 17-52 in alpha-spectrin, or that affect hydrophobic clustering and/or salt-bridge stabilization of the bundled helices, would affect spectrin tetramer formation, and may lead to blood disorders.

Flexibility of aromatic residues in the active-site gorge of acetylcholinesterase: X-ray versus molecular dynamics

The high aromatic content of the deep and narrow active-site gorge of acetylcholinesterase (AChE) is a remarkable feature of this enzyme. Here, we analyze conformational flexibility of the side chains of the 14 conserved aromatic residues in the active-site gorge of Torpedo californica AChE based on the 47 three-dimensional crystal structures available for the native enzyme, and for its complexes and conjugates, and on a 20-ns molecular dynamics (MD) trajectory of the native enzyme. The degree of flexibility of these 14 aromatic side chains is diverse. Although the side-chain conformations of F330 and W279 are both very flexible, the side-chain conformations of F120, W233, W432, Y70, Y121, F288, F290 and F331 appear to be fixed. Residues located on, or adjacent to, the Omega-loop (C67-C94), namely W84, Y130, Y442, and Y334, display different flexibilities in the MD simulations and in the crystal structures. An important outcome of our study is that the majority of the side-chain conformations observed in the 47 Torpedo californica AChE crystal structures are faithfully reproduced by the MD simulation on the native enzyme. Thus, the protein can assume these conformations even in the absence of the ligand that permitted their experimental detection. These observations are pertinent to structure-based drug design.