Accurate positioning of functional residues with robotics-inspired computational protein design.

SignificanceComputational protein design promises to advance applications in medicine and biotechnology by creating proteins with many new and useful functions. However, new functions require the design of specific and often irregular atom-level geometries, which remains a major challenge. Here, we develop computational methods that design and predict local protein geometries with greater accuracy than existing methods. Then, as a proof of concept, we leverage these methods to design new protein conformations in the enzyme ketosteroid isomerase that change the protein's preference for a key functional residue. Our computational methods are openly accessible and can be applied to the design of other intricate geometries customized for new user-defined protein functions.

The structural basis of PTEN regulation by multi-site phosphorylation.

Phosphatase and tensin homolog (PTEN) is a phosphatidylinositol-3,4,5-triphosphate (PIP3) phospholipid phosphatase that is commonly mutated or silenced in cancer. PTEN's catalytic activity, cellular membrane localization and stability are orchestrated by a cluster of C-terminal phosphorylation (phospho-C-tail) events on Ser380, Thr382, Thr383 and Ser385, but the molecular details of this multi-faceted regulation have remained uncertain. Here we use a combination of protein semisynthesis, biochemical analysis, NMR, X-ray crystallography and computational simulations on human PTEN and its sea squirt homolog, VSP, to obtain a detailed picture of how the phospho-C-tail forms a belt around the C2 and phosphatase domains of PTEN. We also visualize a previously proposed dynamic N-terminal α-helix and show that it is key for PTEN catalysis but disordered upon phospho-C-tail interaction. This structural model provides a comprehensive framework for how C-tail phosphorylation can impact PTEN's cellular functions.

A high-affinity human TCR-like antibody detects celiac disease gluten peptide-MHC complexes and inhibits T cell activation.

Antibodies specific for peptides bound to human leukocyte antigen (HLA) molecules are valuable tools for studies of antigen presentation and may have therapeutic potential. Here, we generated human T cell receptor (TCR)-like antibodies toward the immunodominant signature gluten epitope DQ2.5-glia-α2 in celiac disease (CeD). Phage display selection combined with secondary targeted engineering was used to obtain highly specific antibodies with picomolar affinity. The crystal structure of a Fab fragment of the lead antibody 3.C11 in complex with HLA-DQ2.5:DQ2.5-glia-α2 revealed a binding geometry and interaction mode highly similar to prototypic TCRs specific for the same complex. Assessment of CeD biopsy material confirmed disease specificity and reinforced the notion that abundant plasma cells present antigen in the inflamed CeD gut. Furthermore, 3.C11 specifically inhibited activation and proliferation of gluten-specific CD4+ T cells in vitro and in HLA-DQ2.5 humanized mice, suggesting a potential for targeted intervention without compromising systemic immunity.

Directed evolution for high functional production and stability of a challenging G protein-coupled receptor.

Membrane proteins such as G protein-coupled receptors (GPCRs) carry out many fundamental biological functions, are involved in a large number of physiological responses, and are thus important drug targets. To allow detailed biophysical and structural studies, most of these important receptors have to be engineered to overcome their poor intrinsic stability and low expression levels. However, those GPCRs with especially poor properties cannot be successfully optimised even with the current technologies. Here, we present an engineering strategy, based on the combination of three previously developed directed evolution methods, to improve the properties of particularly challenging GPCRs. Application of this novel combination approach enabled the successful selection for improved and crystallisable variants of the human oxytocin receptor, a GPCR with particularly low intrinsic production levels. To analyse the selection results and, in particular, compare the mutations enriched in different hosts, we developed a Next-Generation Sequencing (NGS) strategy that combines long reads, covering the whole receptor, with exceptionally low error rates. This study thus gave insight into the evolution pressure on the same membrane protein in prokaryotes and eukaryotes. Our long-read NGS strategy provides a general methodology for the highly accurate analysis of libraries of point mutants during directed evolution.

Robustification of RosettaAntibody and Rosetta SnugDock.

In recent years, the observed antibody sequence space has grown exponentially due to advances in high-throughput sequencing of immune receptors. The rise in sequences has not been mirrored by a rise in structures, as experimental structure determination techniques have remained low-throughput. Computational modeling, however, has the potential to close the sequence-structure gap. To achieve this goal, computational methods must be robust, fast, easy to use, and accurate. Here we report on the latest advances made in RosettaAntibody and Rosetta SnugDock-methods for antibody structure prediction and antibody-antigen docking. We simplified the user interface, expanded and automated the template database, generalized the kinematics of antibody-antigen docking (which enabled modeling of single-domain antibodies) and incorporated new loop modeling techniques. To evaluate the effects of our updates on modeling accuracy, we developed rigorous tests under a new scientific benchmarking framework within Rosetta. Benchmarking revealed that more structurally similar templates could be identified in the updated database and that SnugDock broadened its applicability without losing accuracy. However, there are further advances to be made, including increasing the accuracy and speed of CDR-H3 loop modeling, before computational approaches can accurately model any antibody.

A compact vocabulary of paratope-epitope interactions enables predictability of antibody-antigen binding.

Antibody-antigen binding relies on the specific interaction of amino acids at the paratope-epitope interface. The predictability of antibody-antigen binding is a prerequisite for de novo antibody and (neo-)epitope design. A fundamental premise for the predictability of antibody-antigen binding is the existence of paratope-epitope interaction motifs that are universally shared among antibody-antigen structures. In a dataset of non-redundant antibody-antigen structures, we identify structural interaction motifs, which together compose a commonly shared structure-based vocabulary of paratope-epitope interactions. We show that this vocabulary enables the machine learnability of antibody-antigen binding on the paratope-epitope level using generative machine learning. The vocabulary (1) is compact, less than 104 motifs; (2) distinct from non-immune protein-protein interactions; and (3) mediates specific oligo- and polyreactive interactions between paratope-epitope pairs. Our work leverages combined structure- and sequence-based learning to demonstrate that machine-learning-driven predictive paratope and epitope engineering is feasible.

An expanded benchmark for antibody-antigen docking and affinity prediction reveals insights into antibody recognition determinants.

Accurate predictive modeling of antibody-antigen complex structures and structure-based antibody design remain major challenges in computational biology, with implications for biotherapeutics, immunity, and vaccines. Through a systematic search for high-resolution structures of antibody-antigen complexes and unbound antibody and antigen structures, in conjunction with identification of experimentally determined binding affinities, we have assembled a non-redundant set of test cases for antibody-antigen docking and affinity prediction. This benchmark more than doubles the number of antibody-antigen complexes and corresponding affinities available in our previous benchmarks, providing an unprecedented view of the determinants of antibody recognition and insights into molecular flexibility. Initial assessments of docking and affinity prediction tools highlight the challenges posed by this diverse set of cases, which includes camelid nanobodies, therapeutic monoclonal antibodies, and broadly neutralizing antibodies targeting viral glycoproteins. This dataset will enable development of advanced predictive modeling and design methods for this therapeutically relevant class of protein-protein interactions.

Macromolecular modeling and design in Rosetta: recent methods and frameworks.

The Rosetta software for macromolecular modeling, docking and design is extensively used in laboratories worldwide. During two decades of development by a community of laboratories at more than 60 institutions, Rosetta has been continuously refactored and extended. Its advantages are its performance and interoperability between broad modeling capabilities. Here we review tools developed in the last 5 years, including over 80 methods. We discuss improvements to the score function, user interfaces and usability. Rosetta is available at http://www.rosettacommons.org.

Novel sampling strategies and a coarse-grained score function for docking homomers, flexible heteromers, and oligosaccharides using Rosetta in CAPRI rounds 37-45.

Critical Assessment of PRediction of Interactions (CAPRI) rounds 37 through 45 introduced larger complexes, new macromolecules, and multistage assemblies. For these rounds, we used and expanded docking methods in Rosetta to model 23 target complexes. We successfully predicted 14 target complexes and recognized and refined near-native models generated by other groups for two further targets. Notably, for targets T110 and T136, we achieved the closest prediction of any CAPRI participant. We created several innovative approaches during these rounds. Since round 39 (target 122), we have used the new RosettaDock 4.0, which has a revamped coarse-grained energy function and the ability to perform conformer selection during docking with hundreds of pregenerated protein backbones. Ten of the complexes had some degree of symmetry in their interactions, so we tested Rosetta SymDock, realized its shortcomings, and developed the next-generation symmetric docking protocol, SymDock2, which includes docking of multiple backbones and induced-fit refinement. Since the last CAPRI assessment, we also developed methods for modeling and designing carbohydrates in Rosetta, and we used them to successfully model oligosaccharide-protein complexes in round 41. Although the results were broadly encouraging, they also highlighted the pressing need to invest in (a) flexible docking algorithms with the ability to model loop and linker motions and in (b) new sampling and scoring methods for oligosaccharide-protein interactions.

Toward the computational design of protein crystals with improved resolution.

Substantial advances have been made in the computational design of protein interfaces over the last 20 years. However, the interfaces targeted by design have typically been stable and high-affinity. Here, we report the development of a generic computational design method to stabilize the weak interactions at crystallographic interfaces. Initially, we analyzed structures reported in the Protein Data Bank to determine whether crystals with more stable interfaces result in higher resolution structures. We found that for 22 variants of a single protein crystallized by a single individual, the Rosetta-calculated `crystal score' correlates with the reported diffraction resolution. We next developed and tested a computational design protocol, seeking to identify point mutations that would improve resolution in a highly stable variant of staphylococcal nuclease (SNase). Using a protocol based on fixed protein backbones, only one of the 11 initial designs crystallized, indicating modeling inaccuracies and forcing us to re-evaluate our strategy. To compensate for slight changes in the local backbone and side-chain environment, we subsequently designed on an ensemble of minimally perturbed protein backbones. Using this strategy, four of the seven designed proteins crystallized. By collecting diffraction data from multiple crystals per design and solving crystal structures, we found that the designed crystals improved the resolution modestly and in unpredictable ways, including altering the crystal space group. Post hoc, in silico analysis of the three observed space groups for SNase showed that the native space group was the lowest scoring for four of six variants (including the wild type), but that resolution did not correlate with crystal score, as it did in the preliminary results. Collectively, our results show that calculated crystal scores can correlate with reported resolution, but that the correlation is absent when the problem is inverted. This outcome suggests that more comprehensive modeling of the crystallographic state is necessary to design high-resolution protein crystals from poorly diffracting crystals.

Caulobacter crescentus Hfq structure reveals a conserved mechanism of RNA annealing regulation.

We have solved the X-ray crystal structure of the RNA chaperone protein Hfq from the alpha-proteobacterium Caulobacter crescentus to 2.15-Å resolution, resolving the conserved core of the protein and the entire C-terminal domain (CTD). The structure reveals that the CTD of neighboring hexamers pack in crystal contacts, and that the acidic residues at the C-terminal tip of the protein interact with positive residues on the rim of Hfq, as has been recently proposed for a mechanism of modulating RNA binding. De novo computational models predict a similar docking of the acidic tip residues against the core of Hfq. We also show that C. crescentus Hfq has sRNA binding and RNA annealing activities and is capable of facilitating the annealing of certain Escherichia coli sRNA:mRNA pairs in vivo. Finally, we describe how the Hfq CTD and its acidic tip residues provide a mechanism to modulate annealing activity and substrate specificity in various bacteria.

Non-H3 CDR template selection in antibody modeling through machine learning.

Antibodies are proteins generated by the adaptive immune system to recognize and counteract a plethora of pathogens through specific binding. This adaptive binding is mediated by structural diversity in the six complementary determining region (CDR) loops (H1, H2, H3, L1, L2 and L3), which also makes accurate structural modeling of CDRs challenging. Both homology and de novo modeling approaches have been used; to date, the former has achieved greater accuracy for the non-H3 loops. The homology modeling of non-H3 CDRs is more accurate because non-H3 CDR loops of the same length and type can be grouped into a few structural clusters. Most antibody-modeling suites utilize homology modeling for the non-H3 CDRs, differing only in the alignment algorithm and how/if they utilize structural clusters. While RosettaAntibody and SAbPred do not explicitly assign query CDR sequences to clusters, two other approaches, PIGS and Kotai Antibody Builder, utilize sequence-based rules to assign CDR sequences to clusters. While the manually curated sequence rules can identify better structural templates, because their curation requires extensive literature search and human effort, they lag behind the deposition of new antibody structures and are infrequently updated. In this study, we propose a machine learning approach (Gradient Boosting Machine [GBM]) to learn the structural clusters of non-H3 CDRs from sequence alone. The GBM method simplifies feature selection and can easily integrate new data, compared to manual sequence rule curation. We compare the classification results using the GBM method to that of RosettaAntibody in a 3-repeat 10-fold cross-validation (CV) scheme on the cluster-annotated antibody database PyIgClassify and we observe an improvement in the classification accuracy of the concerned loops from 84.5% ± 0.24% to 88.16% ± 0.056%. The GBM models reduce the errors in specific cluster membership misclassifications when the involved clusters have relatively abundant data. Based on the factors identified, we suggest methods that can enrich structural classes with sparse data to further improve prediction accuracy in future studies.

Plasma Cells Are the Most Abundant Gluten Peptide MHC-expressing Cells in Inflamed Intestinal Tissues From Patients With Celiac Disease.

BACKGROUND & AIMS: Development of celiac disease is believed to involve the transglutaminase-dependent response of CD4+ T cells toward deamidated gluten peptides in the intestinal mucosa of individuals with specific HLA-DQ haplotypes. We investigated the antigen presentation process during this mucosal immune response. METHODS: We generated monoclonal antibodies (mAbs) specific for the peptide-MHC (pMHC) complex of HLA-DQ2.5 and the immunodominant gluten epitope DQ2.5-glia-α1a using phage display. We used these mAbs to assess gluten peptide presentation and phenotypes of presenting cells by flow cytometry and enzyme-linked immune absorbent spot (ELISPOT) in freshly prepared single-cell suspensions from intestinal biopsies from 40 patients with celiac disease (35 untreated and 5 on a gluten-free diet) as well as 18 subjects with confirmed noninflamed gut mucosa (controls, 12 presumed healthy, 5 undergoing pancreatoduodenectomy, and 1 with potential celiac disease). RESULTS: Using the mAbs, we detected MHC complexes on cells from intestinal biopsies from patients with celiac disease who consume gluten, but not from patients on gluten-free diets. We found B cells and plasma cells to be the most abundant cells that present DQ2.5-glia-α1a in the inflamed mucosa. We identified a subset of plasma cells that expresses B-cell receptors (BCR) specific for gluten peptides or the autoantigen transglutaminase 2 (TG2). Expression of MHC class II (MHCII) was not restricted to these specific plasma cells in patients with celiac disease but was observed in an average 30% of gut plasma cells from patients and controls. CONCLUSIONS: A population of plasma cells from intestinal biopsies of patients with celiac disease express MHCII; this is the most abundant cell type presenting the immunodominant gluten peptide DQ2.5-glia-α1a in the tissues from these patients. These results indicate that plasma cells in the gut can function as antigen-presenting cells and might promote and maintain intestinal inflammation in patients with celiac disease or other inflammatory disorders.

Repertoire Analysis of Antibody CDR-H3 Loops Suggests Affinity Maturation Does Not Typically Result in Rigidification.

Antibodies can rapidly evolve in specific response to antigens. Affinity maturation drives this evolution through cycles of mutation and selection leading to enhanced antibody specificity and affinity. Elucidating the biophysical mechanisms that underlie affinity maturation is fundamental to understanding B-cell immunity. An emergent hypothesis is that affinity maturation reduces the conformational flexibility of the antibody's antigen-binding paratope to minimize entropic losses incurred upon binding. In recent years, computational and experimental approaches have tested this hypothesis on a small number of antibodies, often observing a decrease in the flexibility of the complementarity determining region (CDR) loops that typically comprise the paratope and in particular the CDR-H3 loop, which contributes a plurality of antigen contacts. However, there were a few exceptions and previous studies were limited to a small handful of cases. Here, we determined the structural flexibility of the CDR-H3 loop for thousands of recent homology models of the human peripheral blood cell antibody repertoire using rigidity theory. We found no clear delineation in the flexibility of naïve and antigen-experienced antibodies. To account for possible sources of error, we additionally analyzed hundreds of human and mouse antibodies in the Protein Data Bank through both rigidity theory and B-factor analysis. By both metrics, we observed only a slight decrease in the CDR-H3 loop flexibility when comparing affinity matured antibodies to naïve antibodies, and the decrease was not as drastic as previously reported. Further analysis, incorporating molecular dynamics simulations, revealed a spectrum of changes in flexibility. Our results suggest that rigidification may be just one of many biophysical mechanisms for increasing affinity.

Targeting the CoREST complex with dual histone deacetylase and demethylase inhibitors.

Here we report corin, a synthetic hybrid agent derived from the class I HDAC inhibitor (entinostat) and an LSD1 inhibitor (tranylcypromine analog). Enzymologic analysis reveals that corin potently targets the CoREST complex and shows more sustained inhibition of CoREST complex HDAC activity compared with entinostat. Cell-based experiments demonstrate that corin exhibits a superior anti-proliferative profile against several melanoma lines and cutaneous squamous cell carcinoma lines compared to its parent monofunctional inhibitors but is less toxic to melanocytes and keratinocytes. CoREST knockdown, gene expression, and ChIP studies suggest that corin's favorable pharmacologic effects may rely on an intact CoREST complex. Corin was also effective in slowing tumor growth in a melanoma mouse xenograft model. These studies highlight the promise of a new class of two-pronged hybrid agents that may show preferential targeting of particular epigenetic regulatory complexes and offer unique therapeutic opportunities.

Acidic C-terminal domains autoregulate the RNA chaperone Hfq.

The RNA chaperone Hfq is an Sm protein that facilitates base pairing between bacterial small RNAs (sRNAs) and mRNAs involved in stress response and pathogenesis. Hfq possesses an intrinsically disordered C-terminal domain (CTD) that may tune the function of the Sm domain in different organisms. In Escherichia coli, the Hfq CTD increases kinetic competition between sRNAs and recycles Hfq from the sRNA-mRNA duplex. Here, de novo Rosetta modeling and competitive binding experiments show that the acidic tip of the E. coli Hfq CTD transiently binds the basic Sm core residues necessary for RNA annealing. The CTD tip competes against non-specific RNA binding, facilitates dsRNA release, and prevents indiscriminate DNA aggregation, suggesting that this acidic peptide mimics nucleic acid to auto-regulate RNA binding to the Sm ring. The mechanism of CTD auto-inhibition predicts the chaperone function of Hfq in bacterial genera and illuminates how Sm proteins may evolve new functions.

The Rosetta All-Atom Energy Function for Macromolecular Modeling and Design.

Over the past decade, the Rosetta biomolecular modeling suite has informed diverse biological questions and engineering challenges ranging from interpretation of low-resolution structural data to design of nanomaterials, protein therapeutics, and vaccines. Central to Rosetta's success is the energy function: a model parametrized from small-molecule and X-ray crystal structure data used to approximate the energy associated with each biomolecule conformation. This paper describes the mathematical models and physical concepts that underlie the latest Rosetta energy function, called the Rosetta Energy Function 2015 (REF15). Applying these concepts, we explain how to use Rosetta energies to identify and analyze the features of biomolecular models. Finally, we discuss the latest advances in the energy function that extend its capabilities from soluble proteins to also include membrane proteins, peptides containing noncanonical amino acids, small molecules, carbohydrates, nucleic acids, and other macromolecules.

Modeling and docking of antibody structures with Rosetta.

We describe Rosetta-based computational protocols for predicting the 3D structure of an antibody from sequence (RosettaAntibody) and then docking the antibody to protein antigens (SnugDock). Antibody modeling leverages canonical loop conformations to graft large segments from experimentally determined structures, as well as offering (i) energetic calculations to minimize loops, (ii) docking methodology to refine the VL-VH relative orientation and (iii) de novo prediction of the elusive complementarity determining region (CDR) H3 loop. To alleviate model uncertainty, antibody-antigen docking resamples CDR loop conformations and can use multiple models to represent an ensemble of conformations for the antibody, the antigen or both. These protocols can be run fully automated via the ROSIE web server (http://rosie.rosettacommons.org/) or manually on a computer with user control of individual steps. For best results, the protocol requires roughly 1,000 CPU-hours for antibody modeling and 250 CPU-hours for antibody-antigen docking. Tasks can be completed in under a day by using public supercomputers.

Modeling oblong proteins and water-mediated interfaces with RosettaDock in CAPRI rounds 28-35.

The 28th-35th rounds of the Critical Assessment of PRotein Interactions (CAPRI) served as a practical benchmark for our RosettaDock protein-protein docking protocols, highlighting strengths and weaknesses of the approach. We achieved acceptable or better quality models in three out of 11 targets. For the two α-repeat protein-green fluorescent protein (αrep-GFP) complexes, we used a novel ellipsoidal partial-global docking method (Ellipsoidal Dock) to generate models with 2.2 Å/1.5 Å interface RMSD, capturing 49%/42% of the native contacts, for the 7-/5-repeat αrep complexes. For the DNase-immunity protein complex, we used a new predictor of hydrogen-bonding networks, HBNet with Bridging Waters, to place individual water models at the complex interface; models were generated with 1.8 Å interface RMSD and 12% native water contacts recovered. The targets for which RosettaDock failed to create an acceptable model were typically difficult in general, as six had no acceptable models submitted by any CAPRI predictor. The UCH-L5-RPN13 and UCH-L5-INO80G de-ubiquitinating enzyme-inhibitor complexes comprised inhibitors undergoing significant structural changes upon binding, with the partners being highly interwoven in the docked complexes. Our failure to predict the nucleosome-enzyme complex in Target 95 was largely due to tight constraints we placed on our model based on sparse biochemical data suggesting two specific cross-interface interactions, preventing the correct structure from being sampled. While RosettaDock's three successes show that it is a state-of-the-art docking method, the difficulties with highly flexible and multi-domain complexes highlight the need for better flexible docking and domain-assembly methods. Proteins 2017; 85:479-486. © 2016 Wiley Periodicals, Inc.