The impact of molecular variants, crystallization conditions and the space group on ligand-protein complexes: a case study on bacterial phosphotriesterase

A bacterial phosphotriesterase was employed as an experimental paradigm to examine the effects of multiple factors, such as the molecular constructs, the ligands used during protein expression and purification, the crystallization conditions and the space group, on the visualization of molecular complexes of ligands with a target enzyme. In this case, the ligands used were organophosphates that are fragments of the nerve agents and insecticides on which the enzyme acts as a bioscavenger. 12 crystal structures of various phosphotriesterase constructs obtained by directed evolution were analyzed, with resolutions of up to 1.38 Å. Both apo forms and holo forms, complexed with the organophosphate ligands, were studied. Crystals obtained from three different crystallization conditions, crystallized in four space groups, with and without N-terminal tags, were utilized to investigate the impact of these factors on visualizing the organophosphate complexes of the enzyme. The study revealed that the tags used for protein expression can lodge in the active site and hinder ligand binding. Furthermore, the space group in which the protein crystallizes can significantly impact the visualization of bound ligands. It was also observed that the crystallization precipitants can compete with, and even preclude, ligand binding, leading to false positives or to the incorrect identification of lead drug candidates. One of the co-crystallization conditions enabled the definition of the spaces that accommodate the substituents attached to the P atom of several products of organophosphate substrates after detachment of the leaving group. The crystal structures of the complexes of phosphotriesterase with the organophosphate products reveal similar short interaction distances of the two partially charged O atoms of the P-O bonds with the exposed β-Zn2+ ion and the buried α-Zn2+ ion. This suggests that both Zn2+ ions have a role in stabilizing the transition state for substrate hydrolysis. Overall, this study provides valuable insights into the challenges and considerations involved in studying the crystal structures of ligand-protein complexes, highlighting the importance of careful experimental design and rigorous data analysis in ensuring the accuracy and reliability of the resulting phosphotriesterase-organophosphate structures.

Do "Newly Born" orphan proteins resemble "Never Born" proteins? A study using three deep learning algorithms

"Newly Born" proteins, devoid of detectable homology to any other proteins, known as orphan proteins, occur in a single species or within a taxonomically restricted gene family. They are generated by the expression of novel open reading frames, and appear throughout evolution. We were curious if three recently developed programs for predicting protein structures, namely, AlphaFold2, RoseTTAFold, and ESMFold, might be of value for comparison of such "Newly Born" proteins to random polypeptides with amino acid content similar to that of native proteins, which have been called "Never Born" proteins. The programs were used to compare the structures of two sets of "Never Born" proteins that had been expressed-Group 1, which had been shown experimentally to possess substantial secondary structure, and Group 3, which had been shown to be intrinsically disordered. Overall, although the models generated were scored as being of low quality, they nevertheless revealed some general principles. Specifically, all four members of Group 1 were predicted to be compact by all three algorithms, in agreement with the experimental data, whereas the members of Group 3 were predicted to be very extended, as would be expected for intrinsically disordered proteins, again consistent with the experimental data. These predicted differences were shown to be statistically significant by comparing their accessible surface areas. The three programs were then used to predict the structures of three orphan proteins whose crystal structures had been solved, two of which display novel folds. Surprisingly, only for the protein which did not have a novel fold, and was taxonomically restricted, rather than being a true orphan, did all three algorithms predict very similar, high-quality structures, closely resembling the crystal structure. Finally, they were used to predict the structures of seven orphan proteins with well-identified biological functions, whose 3D structures are not known. Two proteins, which were predicted to be disordered based on their sequences, are predicted by all three structure algorithms to be extended structures. The other five were predicted to be compact structures with only two exceptions in the case of AlphaFold2. All three prediction algorithms make remarkably similar and high-quality predictions for one large protein, HCO_11565, from a nematode. It is conjectured that this is due to many homologs in the taxonomically restricted family of which it is a member, and to the fact that the Dali server revealed several nonrelated proteins with similar folds. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:Proteins:3.

Design of a stable human acid-β-glucosidase: towards improved Gaucher disease therapy and mutation classification

Acid-β-glucosidase (GCase, EC3.2.1.45), the lysosomal enzyme which hydrolyzes the simple glycosphingolipid, glucosylceramide (GlcCer), is encoded by the GBA1 gene. Biallelic mutations in GBA1 cause the human inherited metabolic disorder, Gaucher disease (GD), in which GlcCer accumulates, while heterozygous GBA1 mutations are the highest genetic risk factor for Parkinson's disease (PD). Recombinant GCase (e.g., Cerezyme^® ) is produced for use in enzyme replacement therapy for GD and is largely successful in relieving disease symptoms, except for the neurological symptoms observed in a subset of patients. As a first step toward developing an alternative to the recombinant human enzymes used to treat GD, we applied the PROSS stability-design algorithm to generate GCase variants with enhanced stability. One of the designs, containing 55 mutations compared to wild-type human GCase, exhibits improved secretion and thermal stability. Furthermore, the design has higher enzymatic activity than the clinically used human enzyme when incorporated into an AAV vector, resulting in a larger decrease in the accumulation of lipid substrates in cultured cells. Based on stability-design calculations, we also developed a machine learning-based approach to distinguish benign from deleterious (i.e., disease-causing) GBA1 mutations. This approach gave remarkably accurate predictions of the enzymatic activity of single-nucleotide polymorphisms in the GBA1 gene that are not currently associated with GD or PD. This latter approach could be applied to other diseases to determine risk factors in patients carrying rare mutations.

A practical guide to teaching with Proteopedia

Proteopedia (proteopedia.org) is an open resource to explore the structure-function relationship of proteins and other biomolecules. This guide provides practical advice on how to incorporate Proteopedia into teaching the structure and function of proteins and other biomolecules. For 11 activities, we discuss desired outcomes, setting expectations, preparing students for the tasks, using resources within Proteopedia, and evaluating student work. We point out features of Proteopedia that make it especially suitable for teaching and give examples of how to avoid common pitfalls.

Torpedo californica acetylcholinesterase is stabilized by binding of a divalent metal ion to a novel and versatile 4D motif

Stabilization of Torpedo californica acetylcholinesterase by the divalent cations Ca⁺², Mg⁺², and Mn⁺² was investigated. All three substantially protect the enzyme from thermal inactivation. Electron paramagnetic resonance revealed one high‐affinity binding site for Mn⁺² and several much weaker sites. Differential scanning calorimetry showed a single irreversible thermal transition. All three cations raise both the temperature of the transition and the activation energy, with the transition becoming more cooperative. The crystal structures of the Ca⁺² and Mg⁺² complexes with Torpedo acetylcholinesterase were solved. A principal binding site was identified. In both cases, it consists of four aspartates (a 4D motif), within which the divalent ion is embedded, together with several water molecules. It makes direct contact with two of the aspartates, and indirect contact, via waters, with the other two. The 4D motif has been identified in 31 acetylcholinesterase sequences and 28 butyrylcholinesterase sequences. Zebrafish acetylcholinesterase also contains the 4D motif; it, too, is stabilized by divalent metal ions. The ASSAM server retrieved 200 other proteins that display the 4D motif, in many of which it is occupied by a divalent cation. It is a very versatile motif, since, even though tightly conserved in terms of RMSD values, it can contain from one to as many as three divalent metal ions, together with a variable number of waters. This novel motif, which binds primarily divalent metal ions, is shared by a broad repertoire of proteins. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:Protein_Science:3.

PDB‐ID(s): 7B38, 7B8E and 7B2W;

IceBear: an intuitive and versatile web application for research-data tracking from crystallization experiment to PDB deposition

The web-based IceBear software is a versatile tool to monitor the results of crystallization experiments and is designed to facilitate supervisor and student communications. It also records and tracks all relevant information from crystallization setup to PDB deposition in protein crystallography projects. Fully automated data collection is now possible at several synchrotrons, which means that the number of samples tested at the synchrotron is currently increasing rapidly. Therefore, the protein crystallography research communities at the University of Oulu, Weizmann Institute of Science and Diamond Light Source have joined forces to automate the uploading of sample metadata to the synchrotron. In IceBear, each crystal selected for data collection is given a unique sample name and a crystal page is generated. Subsequently, the metadata required for data collection are uploaded directly to the ISPyB synchrotron database by a shipment module, and for each sample a link to the relevant ISPyB page is stored. IceBear allows notes to be made for each sample during cryocooling treatment and during data collection, as well as in later steps of the structure determination. Protocols are also available to aid the recycling of pins, pucks and dewars when the dewar returns from the synchrotron. The IceBear database is organized around projects, and project members can easily access the crystallization and diffraction metadata for each sample, as well as any additional information that has been provided via the notes. The crystal page for each sample connects the crystallization, diffraction and structural information by providing links to the IceBear drop-viewer page and to the ISPyB data-collection page, as well as to the structure deposited in the Protein Data Bank.

A community proposal to integrate structural bioinformatics activities in ELIXIR (3D-Bioinfo Community)

Structural bioinformatics provides the scientific methods and tools to analyse, archive, validate, and present the biomolecular structure data generated by the structural biology community. It also provides an important link with the genomics community, as structural bioinformaticians also use the extensive sequence data to predict protein structures and their functional sites. A very broad and active community of structural bioinformaticians exists across Europe, and 3D-Bioinfo will establish formal platforms to address their needs and better integrate their activities and initiatives. Our mission will be to strengthen the ties with the structural biology research communities in Europe covering life sciences, as well as chemistry and physics and to bridge the gap between these researchers in order to fully realize the potential of structural bioinformatics. Our Community will also undertake dedicated educational, training and outreach efforts to facilitate this, bringing new insights and thus facilitating the development of much needed innovative applications e.g. for human health, drug and protein design. Our combined efforts will be of critical importance to keep the European research efforts competitive in this respect.

Here we highlight the major European contributions to the field of structural bioinformatics, the most pressing challenges remaining and how Europe-wide interactions, enabled by ELIXIR and its platforms, will help in addressing these challenges and in coordinating structural bioinformatics resources across Europe. In particular, we present recent activities and future plans to consolidate an ELIXIR 3D-Bioinfo Community in structural bioinformatics and propose means to develop better links across the community. These include building new consortia, organising workshops to establish data standards and seeking community agreement on benchmark data sets and strategies. We also highlight existing and planned collaborations with other ELIXIR Communities and other European infrastructures, such as the structural biology community supported by Instruct-ERIC, with whom we have synergies and overlapping common interests.

A Second Look at the Crystal Structures of Drosophila melanogaster Acetylcholinesterase in Complex with Tacrine Derivatives Provides Insights Concerning Catalytic Intermediates and the Design of Specific Insecticides

Over recent decades, crystallographic software for data processing and structure refinement has improved dramatically, resulting in more accurate and detailed crystal structures. It is, therefore, sometimes valuable to have a second look at "old" diffraction data, especially when earlier interpretation of the electron density maps was rather difficult. Here, we present updated crystal structures of Drosophila melanogaster acetylcholinesterase (DmAChE) originally published in [Harel et al., Prot Sci (2000) 9:1063-1072], which reveal features previously unnoticed. Thus, previously unmodeled density in the native active site can be interpreted as stable acetylation of the catalytic serine. Similarly, a strong density in the DmAChE/ZA complex originally attributed to a sulfate ion is better interpreted as a small molecule that is covalently bound. This small molecule can be modeled as either a propionate or a glycinate. The complex is reminiscent of the carboxylate butyrylcholinesterase complexes observed in crystal structures of human butyrylcholinesterases from various sources, and demonstrates the remarkable ability of cholinesterases to stabilize covalent complexes with carboxylates. A very strong peak of density (10 σ) at covalent distance from the Cβ of the catalytic serine is present in the DmAChE/ZAI complex. This can be undoubtedly attributed to an iodine atom, suggesting an unanticipated iodo/hydroxyl exchange between Ser238 and the inhibitor, possibly driven by the intense X-ray irradiation. Finally, the binding of tacrine-derived inhibitors, such as ZA (1DX4) or the iodinated analog, ZAI (1QON) results in the appearance of an open channel that connects the base of the active-site gorge to the solvent. This channel, which arises due to the absence of the conserved tyrosine present in vertebrate cholinesterases, could be exploited to design inhibitors specific to insect cholinesterases. The present study demonstrates that updated processing of older diffraction images, and the re-refinement of older diffraction data, can produce valuable information that could not be detected in the original analysis, and strongly supports the preservation of the diffraction images in public data banks.

Polyproline-rich peptides associated with Torpedo californica acetylcholinesterase tetramers

Acetylcholinesterase (AChE) terminates cholinergic neurotransmission by hydrolyzing acetylcholine. The collagen-tailed AChE tetramer is a product of 2 genes, ACHE and ColQ. The AChE tetramer consists of 4 identical AChE subunits and one polyproline-rich peptide, whose function is to hold the 4 AChE subunits together. Our goal was to determine the amino acid sequence of the polyproline-rich peptide(s) in Torpedo californica AChE (TcAChE) tetramers to aid in the analysis of images that will be acquired by cryo-EM. Collagen-tailed AChE was solubilized from Torpedo californica electric organ, converted to 300 kDa tetramers by digestion with trypsin, and purified by affinity chromatography. Polyproline-rich peptides were released by denaturing the TcAChE tetramers in a boiling water bath, and reducing disulfide bonds with dithiothreitol. Carbamidomethylated peptides were separated from TcAChE protein on a spin filter before they were analyzed by liquid chromatography tandem mass spectrometry on a high resolution Orbitrap Fusion Lumos mass spectrometer. Of the 64 identified collagen-tail (ColQ) peptides, 60 were from the polyproline-rich region near the N-terminus of ColQ. The most abundant proline-rich peptides were SVNKCCLLTPPPPPMFPPPFFTETNILQE, at 40% of total mass-spectral signal intensity, and SVNKCCLLTPPPPPMFPPPFFTETNILQEVDLNNLPLEIKPTEPSCK, at 27% of total intensity. The high abundance of these 2 peptides makes them candidates for the principal form of the polyproline-rich peptide in the trypsin-treated TcAChE tetramers.

Enzyme Evolution: An Epistatic Ratchet versus a Smooth Reversible Transition

Evolutionary trajectories are deemed largely irreversible. In a newly diverged protein, reversion of mutations that led to the functional switch typically results in loss of both the new and the ancestral functions. Nonetheless, evolutionary transitions where reversions are viable have also been described. The structural and mechanistic causes of reversion compatibility versus incompatibility therefore remain unclear. We examined two laboratory evolution trajectories of mammalian paraoxonase-1, a lactonase with promiscuous organophosphate hydrolase (OPH) activity. Both trajectories began with the same active-site mutant, His115Trp, which lost the native lactonase activity and acquired higher OPH activity. A neo-functionalization trajectory amplified the promiscuous OPH activity, whereas the re-functionalization trajectory restored the native activity, thus generating a new lactonase that lacks His115. The His115 revertants of these trajectories indicated opposite trends. Revertants of the neo-functionalization trajectory lost both the evolved OPH and the original lactonase activity. Revertants of the trajectory that restored the original lactonase function were, however, fully active. Crystal structures and molecular simulations show that in the newly diverged OPH, the reverted His115 and other catalytic residues are displaced, thus causing loss of both the original and the new activity. In contrast, in the re-functionalization trajectory, reversion compatibility of the original lactonase activity derives from mechanistic versatility whereby multiple residues can fulfill the same task. This versatility enables unique sequence-reversible compositions that are inaccessible when the active site was repurposed toward a new function.

An intrinsically disordered proteins community for ELIXIR

Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) are now recognised as major determinants in cellular regulation. This white paper presents a roadmap for future e-infrastructure developments in the field of IDP research within the ELIXIR framework. The goal of these developments is to drive the creation of high-quality tools and resources to support the identification, analysis and functional characterisation of IDPs. The roadmap is the result of a workshop titled “An intrinsically disordered protein user community proposal for ELIXIR” held at the University of Padua. The workshop, and further consultation with the members of the wider IDP community, identified the key priority areas for the roadmap including the development of standards for data annotation, storage and dissemination; integration of IDP data into the ELIXIR Core Data Resources; and the creation of benchmarking criteria for IDP-related software. Here, we discuss these areas of priority, how they can be implemented in cooperation with the ELIXIR platforms, and their connections to existing ELIXIR Communities and international consortia. The article provides a preliminary blueprint for an IDP Community in ELIXIR and is an appeal to identify and involve new stakeholders.

The four-helix bundle in cholinesterase dimers: Structural and energetic determinants of stability

The crystal structures of truncated forms of cholinesterases provide good models for assessing the role of non-covalent interactions in dimer assembly in the absence of cross-linking disulfide bonds. These structures identify the four-helix bundle that serves as the interface for formation of acetylcholinesterase and butyrylcholinesterase dimers. Here we performed a theoretical comparison of the structural and energetic factors governing dimerization. This included identification of inter-subunit and intra-subunit hydrogen bonds and hydrophobic interactions, evaluation of solvent-accessible surfaces, and estimation of electrostatic contributions to dimerization. To reveal the contribution to dimerization of individual amino acids within the contact area, free energy perturbation alanine screening was performed. Markov state modelling shows that the loop between the α13 and α14 helices in BChE is unstable, and occupies 4 macro-states. The order of magnitude of mean first passage times between these macrostates is ~10^-8 s. Replica exchange molecular dynamics umbrella sampling calculations revealed that the free energy of human BChE dimerization is -15.5 kcal/mol, while that for human AChE is -26.4 kcal/mol. Thus, the C-terminally truncated human butyrylcholinesterase dimer is substantially less stable than that of human acetylcholinesterase. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:CHEMBIOINT:1.

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.

Comparing NMR and X-ray protein structure: Lindemann-like parameters and NMR disorder

Disordered protein chains and segments are fast becoming a major pathway for our understanding of biological function, especially in more evolved species. However, the standard definition of disordered residues: the inability to constrain them in X-ray derived structures, is not easily applied to NMR derived structures. We carry out a statistical comparison between proteins whose structure was resolved using NMR and using X-ray protocols. We start by establishing a connection between these two protocols for obtaining protein structure. We find a close statistical correspondence between NMR and X-ray structures if fluctuations inherent to the NMR protocol are taken into account. Intuitively this tends to lend support to the validity of both NMR and X-ray protocols in deriving biomolecular models that correspond to in vivo conditions. We then establish Lindemann-like parameters for NMR derived structures and examine what order/disorder cutoffs for these parameters are most consistent with X-ray data and how consistent are they. Finally, we find critical value of [Formula: see text] for the best correspondence between X-ray and NMR derived order/disorder assignment, judged by maximizing the Matthews correlation, and a critical value [Formula: see text] if a balance between false positive and false negative prediction is sought. We examine a few non-conforming cases, and examine the origin of the structure derived in X-ray. This study could help in assigning meaningful disorder from NMR experiments.

Recent developments in structural studies on acetylcholinesterase

This review focuses on several recent developments concerning structure-function relationships in vertebrate acetylcholinesterase. These include studies on high-resolution structures of human acetylcholinesterase and its complexes; the first crystal structure of a snake venom acetylcholinesterase, in which open and closed states of the 'back door' are visualized; a powerful algorithm for redesigning proteins for enhanced expression in prokaryotic systems, as applied to human acetylcholinesterase, which has hitherto been an intractable target; in situ implementation of 'click chemistry' in crystalline acetylcholinesterase, which yields novel insights into the steric and dynamic changes involved in the reaction within the active-site gorge; and a study that demonstrates the effect of crystallization conditions on ligand alignment within a protein complex, in this case the methylene blue-Torpedo californica acetylcholinesterase complex, which highlights the relevance of the precipitant employed to structure-based drug design. This is an article for the special issue XVth International Symposium on Cholinergic Mechanisms.

Automated Structure- and Sequence-Based Design of Proteins for High Bacterial Expression and Stability

Upon heterologous overexpression, many proteins misfold or aggregate, thus resulting in low functional yields. Human acetylcholinesterase (hAChE), an enzyme mediating synaptic transmission, is a typical case of a human protein that necessitates mammalian systems to obtain functional expression. We developed a computational strategy and designed an AChE variant bearing 51 mutations that improved core packing, surface polarity, and backbone rigidity. This variant expressed at ∼2,000-fold higher levels in E. coli compared to wild-type hAChE and exhibited 20°C higher thermostability with no change in enzymatic properties or in the active-site configuration as determined by crystallography. To demonstrate broad utility, we similarly designed four other human and bacterial proteins. Testing at most three designs per protein, we obtained enhanced stability and/or higher yields of soluble and active protein in E. coli. Our algorithm requires only a 3D structure and several dozen sequences of naturally occurring homologs, and is available at http://pross.weizmann.ac.il.

•

A new computational method is used to stabilize five recalcitrant proteins

•

Designed variants show higher expression and stability with unmodified function

•

A designed human acetylcholinesterase variant expresses solubly in bacteria

•

The method is fully automated and implemented on a webserver

In vitro evaluation of the catalytic activity of paraoxonases and phosphotriesterases predicts the enzyme circulatory levels required for in vivo protection against organophosphate intoxications

Catalytic scavengers of organophosphates (OPs) are considered very promising antidote candidates for preventing the adverse effects of OP intoxication as stand alone treatments. This study aimed at correlating the in-vivo catalytic efficiency ((k_cat/K_M)[Enzyme]_pl), established prior to the OP challenge, with the severity of symptoms and survival rates of intoxicated animals. The major objective was to apply a theoretical approach to estimate a lower limit for (k_cat/K_M)[Enzyme]_pl that will be adequate for establishing the desired k_cat/K_M value and plasma concentration of efficacious catalytic bioscavengers. Published data sets by our group and others, from in vivo protection experiments executed in the absence of any supportive medicine, were analyzed. The k_cat/K_M values of eight OP hydrolyzing enzymes and their plasma concentrations in four species exposed to OPs via s.c., i.m. and oral gavage, were analyzed. Our results show that regardless of the OP type and the animal species employed, sign-free animals were observed following bioscavenger treatment provided the theoretically estimated time period required to detoxify 96% of the OP (t_96%) in vivo was ≤10 s. This, for example, can be achieved by an enzyme with k_cat/K_M = 5 × 10⁷ M^-1 min^-1 and a plasma concentration of 0.4 μM ((k_cat/K_M)[Enzyme]_pl = 20 min^-1). Experiments in which animals were intoxicated by i.v. OP injections did not always conform to this rule, and in some cases resulted in high mortality rates. We suggest that in vivo evaluation of catalytic scavengers should avoid the unrealistic bolus i.v. route of OP exposure.

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.

Catalytic stimulation by restrained active-site floppiness--the case of high density lipoprotein-bound serum paraoxonase-1

Despite the abundance of membrane-associated enzymes, the mechanism by which membrane binding stabilizes these enzymes and stimulates their catalysis remains largely unknown. Serum paraoxonase-1 (PON1) is a lipophilic lactonase whose stability and enzymatic activity are dramatically stimulated when associated with high-density lipoprotein (HDL) particles. Our mutational and structural analyses, combined with empirical valence bond simulations, reveal a network of hydrogen bonds that connect HDL binding residues with Asn168--a key catalytic residue residing >15Å from the HDL contacting interface. This network ensures precise alignment of N168, which, in turn, ligates PON1's catalytic calcium and aligns the lactone substrate for catalysis. HDL binding restrains the overall motion of the active site and particularly of N168, thus reducing the catalytic activation energy barrier. We demonstrate herein that disturbance of this network, even at its most far-reaching periphery, undermines PON1's activity. Membrane binding thus immobilizes long-range interactions via second- and third-shell residues that reduce the active site's floppiness and pre-organize the catalytic residues. Although this network is critical for efficient catalysis, as demonstrated here, unraveling these long-rage interaction networks is challenging, let alone their implementation in artificial enzyme design.

What's in a name? Why these proteins are intrinsically disordered: Why these proteins are intrinsically disordered

"What's in a name? That which we call a rose By any other name would smell as sweet." From "Romeo and Juliet", William Shakespeare (1594) This article opens a series of publications on disambiguation of the basic terms used in the field of intrinsically disordered proteins. We start from the beginning, namely from the explanation of what the expression "intrinsically disordered protein" actually means and why this particular term has been chosen as the common denominator for this class of proteins characterized by broad structural, dynamic and functional characteristics.

Catalytic metal ion rearrangements underline promiscuity and evolvability of a metalloenzyme

Although largely deemed as structurally conserved, catalytic metal ion sites can rearrange, thereby contributing to enzyme evolvability. Here, we show that in paraoxonase-1, a lipo-lactonase, catalytic promiscuity and divergence into an organophosphate hydrolase are correlated with an alternative mode of the catalytic Ca(2+). We describe the crystal structures of active-site mutants bearing mutations at position 115. The histidine at this position acts as a base to activate the lactone-hydrolyzing water molecule. Mutations to Trp or Gln indeed diminish paraoxonase-1's lactonase activity; however, the promiscuous organophosphate hydrolase activity is enhanced. The structures reveal a 1.8-Å upward displacement towards the enzyme's surface of the catalytic Ca(2+) in the His115 mutants and configurational changes in the ligating side chains and water molecules, relative to the wild-type enzyme. Biochemical analysis and molecular dynamics simulations suggest that this alternative, upward metal mode mediates the promiscuous hydrolysis of organophosphates. The upward Ca(2+) mode observed in the His115 mutants also appears to mediate the wild type's paraoxonase activity. However, whereas the upward mode dominates in the Trp115 mutant, it is scarcely populated in wild type. Thus, the plasticity of active-site metal ions may permit alternative, latent, promiscuous activities and also provide the basis for the divergence of new enzymatic functions.

Evolved stereoselective hydrolases for broad-spectrum G-type nerve agent detoxification

A preferred strategy for preventing nerve agents intoxication is catalytic scavenging by enzymes that hydrolyze them before they reach their targets. Using directed evolution, we simultaneously enhanced the activity of a previously described serum paraoxonase 1 (PON1) variant for hydrolysis of the toxic S(P) isomers of the most threatening G-type nerve agents. The evolved variants show ≤340-fold increased rates and catalytic efficiencies of 0.2-5 × 10(7) M(-1) min(-1). Our selection for prevention of acetylcholinesterase inhibition also resulted in the complete reversion of PON1's stereospecificity, from an enantiomeric ratio (E) < 6.3 × 10(-4) in favor of the R(P) isomer of a cyclosarin analog in wild-type PON1, to E > 2,500 for the S(P) isomer in an evolved variant. Given their ability to hydrolyze G-agents, these evolved variants may serve as broad-range G-agent prophylactics.

Structural and functional characterization of the interaction of the photosensitizing probe methylene blue with Torpedo californica acetylcholinesterase

The photosensitizer, methylene blue (MB), generates singlet oxygen that irreversibly inhibits Torpedo californica acetylcholinesterase (TcAChE). In the dark, it inhibits reversibly. Binding is accompanied by a bathochromic absorption shift, used to demonstrate displacement by other acetylcholinesterase inhibitors interacting with the catalytic "anionic" subsite (CAS), the peripheral "anionic" subsite (PAS), or bridging them. MB is a noncompetitive inhibitor of TcAChE, competing with reversible inhibitors directed at both "anionic" subsites, but a single site is involved in inhibition. MB also quenches TcAChE's intrinsic fluorescence. It binds to TcAChE covalently inhibited by a small organophosphate (OP), but not an OP containing a bulky pyrene. Differential scanning calorimetry shows an ~8° increase in the denaturation temperature of the MB/TcAChE complex relative to native TcAChE, and a less than twofold increase in cooperativity of the transition. The crystal structure reveals a single MB stacked against Trp279 in the PAS, oriented down the gorge toward the CAS; it is plausible that irreversible inhibition is associated with photooxidation of this residue and others within the active-site gorge. The kinetic and spectroscopic data showing that inhibitors binding at the CAS can impede binding of MB are reconciled by docking studies showing that the conformation adopted by Phe330, midway down the gorge, in the MB/TcAChE crystal structure, precludes simultaneous binding of a second MB at the CAS. Conversely, binding of ligands at the CAS dislodges MB from its preferred locus at the PAS. The data presented demonstrate that TcAChE is a valuable model for understanding the molecular basis of local photooxidative damage.

Backbone and side chain NMR assignments for the intrinsically disordered cytoplasmic domain of human neuroligin-3

Neuroligins act as heterophilic adhesion molecules at neuronal synapses. Their cytoplasmic domains interact with synaptic scaffolding proteins, and have been shown to be intrinsically disordered. Here we report the backbone and side chain (1)H, (13)C and (15)N resonance assignments for the cytoplasmic domain of human neuroligin 3.

Structure of recombinant human carboxylesterase 1 isolated from whole cabbage looper larvae

The use of whole insect larvae as a source of recombinant proteins offers a more cost-effective method of producing large quantities of human proteins than conventional cell-culture approaches. Human carboxylesterase 1 has been produced in and isolated from whole Trichoplusia ni larvae. The recombinant protein was crystallized and its structure was solved to 2.2 resolution. The results indicate that the larvae-produced enzyme is essentially identical to that isolated from cultured Sf21 cells, supporting the use of this expression system to produce recombinant enzymes for crystallization studies.

Catalytic versatility and backups in enzyme active sites: the case of serum paraoxonase 1

The origins of enzyme specificity are well established. However, the molecular details underlying the ability of a single active site to promiscuously bind different substrates and catalyze different reactions remain largely unknown. To better understand the molecular basis of enzyme promiscuity, we studied the mammalian serum paraoxonase 1 (PON1) whose native substrates are lipophilic lactones. We describe the crystal structures of PON1 at a catalytically relevant pH and of its complex with a lactone analogue. The various PON1 structures and the analysis of active-site mutants guided the generation of docking models of the various substrates and their reaction intermediates. The models suggest that promiscuity is driven by coincidental overlaps between the reactive intermediate for the native lactonase reaction and the ground and/or intermediate states of the promiscuous reactions. This overlap is also enabled by different active-site conformations: the lactonase activity utilizes one active-site conformation whereas the promiscuous phosphotriesterase activity utilizes another. The hydrolysis of phosphotriesters, and of the aromatic lactone dihydrocoumarin, is also driven by an alternative catalytic mode that uses only a subset of the active-site residues utilized for lactone hydrolysis. Indeed, PON1's active site shows a remarkable level of networking and versatility whereby multiple residues share the same task and individual active-site residues perform multiple tasks (e.g., binding the catalytic calcium and activating the hydrolytic water). Overall, the coexistence of multiple conformations and alternative catalytic modes within the same active site underlines PON1's promiscuity and evolutionary potential.

In vitro detoxification of cyclosarin in human blood pre-incubated ex vivo with recombinant serum paraoxonases

An ex vivo protocol was developed to assay the antidotal capacity of rePON1 variants to protect endogenous acetylcholinesterase and butyrylcholinesterase in human whole blood against OP nerve agents. This protocol permitted us to address the relationship between blood rePON1 concentrations, their kinetic parameters, and the level of protection conferred by rePON1 on the cholinesterases in human blood, following a challenge with cyclosarin (GF). The experimental data thus obtained were in good agreement with the predicted percent residual activities of blood cholinesterases calculated on the basis of the rate constants for inhibition of human acetylcholinesterase and butyrylcholinesterase by GF, the concentration of the particular rePON1 variant, and its k(cat)/K(m) value for GF. This protocol thus provides a rapid and reliable ex vivo screening tool for identification of rePON1 bioscavenger candidates suitable for protection of humans against organophosphorus-based toxicants. The results also permitted the refinement of a mathematical model for estimating the efficacious dose of rePON1s variants required for prophylaxis in humans.

Backdoor opening mechanism in acetylcholinesterase based on X-ray crystallography and molecular dynamics simulations

The transient opening of a backdoor in the active-site wall of acetylcholinesterase, one of nature's most rapid enzymes, has been suggested to contribute to the efficient traffic of substrates and products. A crystal structure of Torpedo californica acetylcholinesterase in complex with the peripheral-site inhibitor aflatoxin is now presented, in which a tyrosine at the bottom of the active-site gorge rotates to create a 3.4-Å wide exit channel. Molecular dynamics simulations show that the opening can be further enlarged by movement of Trp84. The crystallographic and molecular dynamics simulation data thus point to the interface between Tyr442 and Trp84 as the key element of a backdoor, whose opening permits rapid clearance of catalysis products from the active site. Furthermore, the crystal structure presented provides a novel template for rational design of inhibitors and reactivators, including anti-Alzheimer drugs and antidotes against organophosphate poisoning.

Context-dependent resistance to proteolysis of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs), also known as intrinsically unstructured proteins (IUPs), lack a well-defined 3D structure in vitro and, in some cases, also in vivo. Here, we discuss the question of proteolytic sensitivity of IDPs, with a view to better explaining their in vivo characteristics. After an initial assessment of the status of IDPs in vivo, we briefly survey the intracellular proteolytic systems. Subsequently, we discuss the evidence for IDPs being inherently sensitive to proteolysis. Such sensitivity would not, however, result in enhanced degradation if the protease-sensitive sites were sequestered. Accordingly, IDP access to and degradation by the proteasome, the major proteolytic complex within eukaryotic cells, are discussed in detail. The emerging picture appears to be that IDPs are inherently sensitive to proteasomal degradation along the lines of the "degradation by default" model. However, available data sets of intracellular protein half-lives suggest that intrinsic disorder does not imply a significantly shorter half-life. We assess the power of available systemic half-life measurements, but also discuss possible mechanisms that could protect IDPs from intracellular degradation. Finally, we discuss the relevance of the proteolytic sensitivity of IDPs to their function and evolution.

Structure of estradiol metal chelate and estrogen receptor complex: the basis for designing a new class of selective estrogen receptor modulators

Selective estrogen receptor modulators, such as 17β-estradiol derivatives bound to metal complexes, have been synthesized as targeted probes for the diagnosis and treatment of breast cancer. Here, we report the detailed 3D structure of estrogen receptor α ligand-binding domain (ERα-LBD) bound with a novel estradiol-derived metal complex, estradiol-pyridine tetra acetate europium(III), at 2.6 Å resolution. This structure provides important information pertinent to the design of novel functional ERα targeted probes for clinical applications.

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

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