VWD domain stabilization by autocatalytic Asp-Pro cleavage.

Domains known as von Willebrand factor type D (VWD) are found in extracellular and cell-surface proteins including von Willebrand factor, mucins, and various signaling molecules and receptors. Many VWD domains have a glycine-aspartate-proline-histidine (GDPH) amino-acid sequence motif, which is hydrolytically cleaved post-translationally between the aspartate (Asp) and proline (Pro). The Fc IgG binding protein (FCGBP), found in intestinal mucus secretions and other extracellular environments, contains 13 VWD domains, 11 of which have a GDPH cleavage site. In this study, we investigated the structural and biophysical consequences of Asp-Pro peptide cleavage in a representative FCGBP VWD domain. We found that endogenous Asp-Pro cleavage increases the resistance of the domain to exogenous proteolytic degradation. Tertiary structural interactions made by the newly generated chain termini, as revealed by a crystal structure of an FCGBP segment containing the VWD domain, may explain this observation. Notably, the Gly-Asp peptide bond, upstream of the cleavage site, assumed the cis configuration in the structure. In addition to these local features of the cleavage site, a global organizational difference was seen when comparing the FCGBP segment structure with the numerous other structures containing the same set of domains. Together, these data illuminate the outcome of GDPH cleavage and demonstrate the plasticity of proteins with VWD domains, which may contribute to their evolution for function in a dynamic extracellular environment.

PGR5-LIKE PHOTOSYNTHETIC PHENOTYPE1A redox states alleviate photoinhibition during changes in light intensity.

Plants have evolved photosynthetic regulatory mechanisms to maintain homeostasis in response to light changes during diurnal transitions and those caused by passing clouds or by wind. One such adaptation directs photosynthetic electron flow to a cyclic pathway to alleviate excess energy surges. Here, we assign a function to regulatory cysteines of PGR5-like protein 1A (PGRL1A), a constituent of the PROTON GRADIENT REGULATION5 (PGR5)-dependent cyclic electron flow (CEF) pathway. During step increases from darkness to low light intensity in Arabidopsis (Arabidopsis thaliana), the intermolecular disulfide of the PGRL1A 59-kDa complex was reduced transiently within seconds to the 28-kDa form. In contrast, step increases from darkness to high light stimulated a stable, partially reduced redox state in PGRL1A. Mutations of 2 cysteines in PGRL1A, Cys82 and Cys183, resulted in a constitutively pseudo-reduced state. The mutant displayed higher proton motive force (PMF) and nonphotochemical quenching (NPQ) than the wild type (WT) and showed altered donor and acceptor dynamic flow around PSI. These changes were found to correspond with the redox state of PGRL1A. Continuous light regimes did not affect mutant growth compared to the WT. However, under fluctuating regimes of high light, the mutant showed better growth than the WT. In contrast, in fluctuating regimes of low light, the mutant displayed a growth penalty that can be attributed to constant stimulation of CEF under low light. Treatment with photosynthetic inhibitors indicated that PGRL1A redox state control depends on the penultimate Fd redox state. Our results showed that redox state changes in PGRL1A are crucial to optimize photosynthesis.

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.

Diversity of CysD domains in gel-forming mucins.

CysD domains are disulfide-rich modules embedded within long O-glycosylated regions of mucin glycoproteins. CysD domains are thought to mediate intermolecular adhesion during the intracellular bioassembly of mucin polymers and perhaps also after secretion in extracellular mucus hydrogels. The human genome encodes 18 CysD domains distributed across three different mucins. To date, experimental structural information is available only for the first CysD domain (CysD1) of the intestinal mucin MUC2, which is one of the most divergent of the CysDs. To provide experimental data on a CysD that is representative of a larger branch of the fold family, we determined the crystal structure of the seventh CysD domain (CysD7) from MUC5AC, a mucin found in the respiratory tract and stomach. The MUC5AC CysD7 structure revealed a single calcium-binding site, contrasting with the two sites in MUC2 CysD1. The MUC5AC CysD7 structure also contained an additional α-helix absent from MUC2 CysD1, with potential functional implications for intermolecular interactions. Lastly, the experimental structure emphasized the flexibility of the loop analogous to the main adhesion loop of MUC2 CysD1, suggesting that both sequence divergence and physical plasticity in this region may contribute to the adaptation of mucin CysD domains.

Specific quinone reductase 2 inhibitors reduce metabolic burden and reverse Alzheimer's disease phenotype in mice.

Biological aging can be described as accumulative, prolonged metabolic stress and is the major risk factor for cognitive decline and Alzheimer's disease (AD). Recently, we identified and described a quinone reductase 2 (QR2) pathway in the brain, in which QR2 acts as a removable memory constraint and metabolic buffer within neurons. QR2 becomes overexpressed with age, and it is possibly a novel contributing factor to age-related metabolic stress and cognitive deficit. We found that, in human cells, genetic removal of QR2 produced a shift in the proteome opposing that found in AD brains while simultaneously reducing oxidative stress. We therefore created highly specific QR2 inhibitors (QR2is) to enable evaluation of chronic QR2 inhibition as a means to reduce biological age-related metabolic stress and cognitive decline. QR2is replicated results obtained by genetic removal of QR2, while local QR2i microinjection improved hippocampal and cortical-dependent learning in rats and mice. Continuous consumption of QR2is in drinking water improved cognition and reduced pathology in the brains of AD-model mice (5xFAD), with a noticeable between-sex effect on treatment duration. These results demonstrate the importance of QR2 activity and pathway function in the healthy and neurodegenerative brain and what we believe to be the great therapeutic potential of QR2is as first-in-class drugs.

Phylogeny and biogeography of the algal DMS-releasing enzyme in the global ocean.

Phytoplankton produce the volatile dimethyl sulfide (DMS), an important infochemical mediating microbial interactions, which is also emitted to the atmosphere and affecting the global climate. Albeit the enzymatic source for DMS in eukaryotes was elucidated, namely a DMSP lyase (DL) called Alma1, we still lack basic knowledge regarding its taxonomic distribution. We defined unique sequence motifs which enable the identification of DL homologs (DLHs) in model systems and environmental populations. We used these motifs to predict DLHs in diverse algae by analyzing hundreds of genomic and transcriptomic sequences from model systems under stress conditions and from environmental samples. Our findings show that the DL enzyme is more taxonomically widespread than previously thought, as it is encoded by known algal taxa as haptophytes and dinoflagellates, but also by chlorophytes, pelagophytes and diatoms, which were conventionally considered to lack the DL enzyme. By exploring the Tara Oceans database, we showed that DLHs are widespread across the oceans and are predominantly expressed by dinoflagellates. Certain dinoflagellate DLHs were differentially expressed between the euphotic and mesopelagic zones, suggesting a functional specialization and an involvement in the metabolic plasticity of mixotrophic dinoflagellates. In specific regions as the Southern Ocean, DLH expression by haptophytes and diatoms was correlated with environmental drivers such as nutrient availability. The expanded repertoire of putative DL enzymes from diverse microbial origins and geographic niches suggests new potential players in the marine sulfur cycle and provides a foundation to study the cellular function of the DL enzyme in marine microbes.

gUMI-BEAR, a modular, unsupervised population barcoding method to track variants and evolution at high resolution.

Cellular lineage tracking provides a means to observe population makeup at the clonal level, allowing exploration of heterogeneity, evolutionary and developmental processes and individual clones' relative fitness. It has thus contributed significantly to understanding microbial evolution, organ differentiation and cancer heterogeneity, among others. Its use, however, is limited because existing methods are highly specific, expensive, labour-intensive, and, critically, do not allow the repetition of experiments. To address these issues, we developed gUMI-BEAR (genomic Unique Molecular Identifier Barcoded Enriched Associated Regions), a modular, cost-effective method for tracking populations at high resolution. We first demonstrate the system's application and resolution by applying it to track tens of thousands of Saccharomyces cerevisiae lineages growing together under varying environmental conditions applied across multiple generations, revealing fitness differences and lineage-specific adaptations. Then, we demonstrate how gUMI-BEAR can be used to perform parallel screening of a huge number of randomly generated variants of the Hsp82 gene. We further show how our method allows isolation of variants, even if their frequency in the population is low, thus enabling unsupervised identification of modifications that lead to a behaviour of interest.

Stable Mammalian Serum Albumins Designed for Bacterial Expression.

Albumin is the most abundant protein in the blood serum of mammals and has essential carrier and physiological roles. Albumins are also used in a wide variety of molecular and cellular experiments and in the cultivated meat industry. Despite their importance, however, albumins are challenging for heterologous expression in microbial hosts, likely due to 17 conserved intramolecular disulfide bonds. Therefore, albumins used in research and biotechnological applications either derive from animal serum, despite severe ethical and reproducibility concerns, or from recombinant expression in yeast or rice. We use the PROSS algorithm to stabilize human and bovine serum albumins, finding that all are highly expressed in E. coli. Design accuracy is verified by crystallographic analysis of a human albumin variant with 16 mutations. This albumin variant exhibits ligand binding properties similar to those of the wild type. Remarkably, a design with 73 mutations relative to human albumin exhibits over 40 °C improved stability and is stable beyond the boiling point of water. Our results suggest that proteins with many disulfide bridges have the potential to exhibit extreme stability when subjected to design. The designed albumins may be used to make economical, reproducible, and animal-free reagents for molecular and cell biology. They also open the way to high-throughput screening to study and enhance albumin carrier properties.

Allosteric communication between ligand binding domains modulates substrate inhibition in adenylate kinase.

Enzymes play a vital role in life processes; they control chemical reactions and allow functional cycles to be synchronized. Many enzymes harness large-scale motions of their domains to achieve tremendous catalytic prowess and high selectivity for specific substrates. One outstanding example is provided by the three-domain enzyme adenylate kinase (AK), which catalyzes phosphotransfer between ATP to AMP. Here we study the phenomenon of substrate inhibition by AMP and its correlation with domain motions. Using single-molecule FRET spectroscopy, we show that AMP does not block access to the ATP binding site, neither by competitive binding to the ATP cognate site nor by directly closing the LID domain. Instead, inhibitory concentrations of AMP lead to a faster and more cooperative domain closure by ATP, leading in turn to an increased population of the closed state. The effect of AMP binding can be modulated through mutations throughout the structure of the enzyme, as shown by the screening of an extensive AK mutant library. The mutation of multiple conserved residues reduces substrate inhibition, suggesting that substrate inhibition is an evolutionary well conserved feature in AK. Combining these insights, we developed a model that explains the complex activity of AK, particularly substrate inhibition, based on the experimentally observed opening and closing rates. Notably, the model indicates that the catalytic power is affected by the microsecond balance between the open and closed states of the enzyme. Our findings highlight the crucial role of protein motions in enzymatic activity.

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.

The identification and characterization of an oxalyl-CoA synthetase from grass pea (Lathyrus sativus L.).

Oxalic acid is a small metabolite found in many plants. It serves as protection from herbivores, a chelator of metal ions, a regulator of calcium levels, and additional tasks. However, it is also a strong di-carboxylic acid that can compromise plant viability by reducing cellular pH. Several metabolic pathways have evolved to control oxalate levels in plants by enzymatic degradation. Among them is the pathway that utilizes oxalyl-CoA synthetase (OCS, EC 6.2.1.8) and ATP to convert oxalate to oxalyl-CoA. Oxalyl-CoA can then be degraded to CO2 or utilized as a precursor for the synthesis of other compounds. In grass pea (Lathyrus sativus L.), a grain legume grown in Asia and Africa for human and animal consumption, the neurotoxic compound ß-N-oxalyl-l-α,ß-diaminopropionic acid (ß-ODAP) is synthesized from oxalyl-CoA and l-α,ß-diaminopropionic acid (l-DAPA). Here, we report on the identification and characterization of oxalyl CoA-synthetase from grass pea (LsOCS). The gene encoding LsOCS was amplified from grass pea, and then expressed and purified from E. coli cells as an untagged, monomeric protein of 56 kDa. Its catalytic efficiency with oxalate, K oxalate M = 71.5 ± 13.3 µM, V max = 8.2 ± 0.8 µmole min-1 mg-1, was similar to that of OCS homologs from Arabidopsis thaliana (AtAAE3) and Medicago truncatula (MtAAE3). The enzyme was crystalized in complex with AMP and is the first OCS whose structure was determined in the thioester-forming conformation. Finally, we propose that substituting LsOCS with an oxalate oxidase or decarboxylase may reduce the levels of ß-ODAP in grass pea.

Identification and characterization of the key enzyme in the biosynthesis of the neurotoxin ß-ODAP in grass pea.

Grass pea (Lathyrus sativus L.) is a grain legume commonly grown in Asia and Africa for food and forage. It is a highly nutritious and robust crop, capable of surviving both droughts and floods. However, it produces a neurotoxic compound, ß-N-oxalyl-L-α,ß-diaminopropionic acid (ß-ODAP), which can cause a severe neurological disorder when consumed as a primary diet component. While the catalytic activity associated with ß-ODAP formation was demonstrated more than 50 years ago, the enzyme responsible for this activity has not been identified. Here, we report on the identity, activity, 3D structure, and phylogenesis of this enzyme-ß-ODAP synthase (BOS). We show that BOS belongs to the benzylalcohol O-acetyltransferase, anthocyanin O-hydroxycinnamoyltransferase, anthranilate N-hydroxycinnamoyl/benzoyltransferase, deacetylvindoline 4-O-acetyltransferase superfamily of acyltransferases and is structurally similar to hydroxycinnamoyl transferase. Using molecular docking, we propose a mechanism for its catalytic activity, and using heterologous expression in tobacco leaves (Nicotiana benthamiana), we demonstrate that expression of BOS in the presence of its substrates is sufficient for ß-ODAP production in vivo. The identification of BOS may pave the way toward engineering ß-ODAP-free grass pea cultivars, which are safe for human and animal consumption.

Computationally designed dual-color MRI reporters for noninvasive imaging of transgene expression.

Imaging of gene-expression patterns in live animals is difficult to achieve with fluorescent proteins because tissues are opaque to visible light. Imaging of transgene expression with magnetic resonance imaging (MRI), which penetrates to deep tissues, has been limited by single reporter visualization capabilities. Moreover, the low-throughput capacity of MRI limits large-scale mutagenesis strategies to improve existing reporters. Here we develop an MRI system, called GeneREFORM, comprising orthogonal reporters for two-color imaging of transgene expression in deep tissues. Starting from two promiscuous deoxyribonucleoside kinases, we computationally designed highly active, orthogonal enzymes ('reporter genes') that specifically phosphorylate two MRI-detectable synthetic deoxyribonucleosides ('reporter probes'). Systemically administered reporter probes exclusively accumulate in cells expressing the designed reporter genes, and their distribution is displayed as pseudo-colored MRI maps based on dynamic proton exchange for noninvasive visualization of transgene expression. We envision that future extensions of GeneREFORM will pave the way to multiplexed deep-tissue mapping of gene expression in live animals.

SARS-CoV-2 variant prediction and antiviral drug design are enabled by RBD in vitro evolution.

SARS-CoV-2 variants of interest and concern will continue to emerge for the duration of the COVID-19 pandemic. To map mutations in the receptor-binding domain (RBD) of the spike protein that affect binding to angiotensin-converting enzyme 2 (ACE2), the receptor for SARS-CoV-2, we applied in vitro evolution to affinity-mature the RBD. Multiple rounds of random mutagenic libraries of the RBD were sorted against decreasing concentrations of ACE2, resulting in the selection of higher affinity RBD binders. We found that mutations present in more transmissible viruses (S477N, E484K and N501Y) were preferentially selected in our high-throughput screen. Evolved RBD mutants include prominently the amino acid substitutions found in the RBDs of B.1.620, B.1.1.7 (Alpha), B1.351 (Beta) and P.1 (Gamma) variants. Moreover, the incidence of RBD mutations in the population as presented in the GISAID database (April 2021) is positively correlated with increased binding affinity to ACE2. Further in vitro evolution increased binding by 1,000-fold and identified mutations that may be more infectious if they evolve in the circulating viral population, for example, Q498R is epistatic to N501Y. We show that our high-affinity variant RBD-62 can be used as a drug to inhibit infection with SARS-CoV-2 and variants Alpha, Beta and Gamma in vitro. In a model of SARS-CoV-2 challenge in hamster, RBD-62 significantly reduced clinical disease when administered before or after infection. A 2.9 Å cryo-electron microscopy structure of the high-affinity complex of RBD-62 and ACE2, including all rapidly spreading mutations, provides a structural basis for future drug and vaccine development and for in silico evaluation of known antibodies.

Conformation-Specific Inhibitory Anti-MMP-7 Monoclonal Antibody Sensitizes Pancreatic Ductal Adenocarcinoma Cells to Chemotherapeutic Cell Kill.

Matrix metalloproteases (MMPs) undergo post-translational modifications including pro-domain shedding. The activated forms of these enzymes are effective drug targets, but generating potent biological inhibitors against them remains challenging. We report the generation of anti-MMP-7 inhibitory monoclonal antibody (GSM-192), using an alternating immunization strategy with an active site mimicry antigen and the activated enzyme. Our protocol yielded highly selective anti-MMP-7 monoclonal antibody, which specifically inhibits MMP-7's enzyme activity with high affinity (IC50 = 132 ± 10 nM). The atomic model of the MMP-7-GSM-192 Fab complex exhibited antibody binding to unique epitopes at the rim of the enzyme active site, sterically preventing entry of substrates into the catalytic cleft. In human PDAC biopsies, tissue staining with GSM-192 showed characteristic spatial distribution of activated MMP-7. Treatment with GSM-192 in vitro induced apoptosis via stabilization of cell surface Fas ligand and retarded cell migration. Co-treatment with GSM-192 and chemotherapeutics, gemcitabine and oxaliplatin elicited a synergistic effect. Our data illustrate the advantage of precisely targeting catalytic MMP-7 mediated disease specific activity.

A Molecular Lid Mechanism of K+ Channel Blocker Action Revealed by a Cone Peptide.

Many venomous organisms carry in their arsenal short polypeptides that block K+ channels in a highly selective manner. These toxins may compete with the permeating ions directly via a "plug" mechanism or indirectly via a "pore-collapse" mechanism. An alternative "lid" mechanism was proposed but remained poorly defined. Here we study the Drosophila Shaker channel block by Conkunitzin-S1 and Conkunitzin-C3, two highly similar toxins derived from cone venom. Despite their similarity, the two peptides exhibited differences in their binding poses and biophysical assays, implying discrete action modes. We show that while Conkunitzin-S1 binds tightly to the channel turret and acts via a "pore-collapse" mechanism, Conkunitzin-C3 does not contact this region. Instead, Conk-C3 uses a non-conserved Arg to divert the permeant ions and trap them in off-axis cryptic sites above the SF, a mechanism we term a "molecular-lid". Our study provides an atomic description of the "lid" K+ blocking mode and offers valuable insights for the design of therapeutics based on venom peptides.

Community-Wide Experimental Evaluation of the PROSS Stability-Design Method.

Recent years have seen a dramatic improvement in protein-design methodology. Nevertheless, most methods demand expert intervention, limiting their widespread adoption. By contrast, the PROSS algorithm for improving protein stability and heterologous expression levels has been successfully applied to a range of challenging enzymes and binding proteins. Here, we benchmark the application of PROSS as a stand-alone tool for protein scientists with no or limited experience in modeling. Twelve laboratories from the Protein Production and Purification Partnership in Europe (P4EU) challenged the PROSS algorithm with 14 unrelated protein targets without support from the PROSS developers. For each target, up to six designs were evaluated for expression levels and in some cases, for thermal stability and activity. In nine targets, designs exhibited increased heterologous expression levels either in prokaryotic and/or eukaryotic expression systems under experimental conditions that were tailored for each target protein. Furthermore, we observed increased thermal stability in nine of ten tested targets. In two prime examples, the human Stem Cell Factor (hSCF) and human Cadherin-Like Domain (CLD12) from the RET receptor, the wild type proteins were not expressible as soluble proteins in E. coli, yet the PROSS designs exhibited high expression levels in E. coli and HEK293 cells, respectively, and improved thermal stability. We conclude that PROSS may improve stability and expressibility in diverse cases, and that improvement typically requires target-specific expression conditions. This study demonstrates the strengths of community-wide efforts to probe the generality of new methods and recommends areas for future research to advance practically useful algorithms for protein science.

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.

Correction: Optimizing antibody affinity and stability by the automated design of the variable light-heavy chain interfaces.

[This corrects the article DOI: 10.1371/journal.pcbi.1007207.].

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.