Deep sequencing methods for protein engineering and design

An integrated technology for quantitative wide mutational scanning of human antibody Fab libraries

Antibodies are engineerable quantities in medicine. Learning antibody molecular recognition would enable the in silico design of high affinity binders against nearly any proteinaceous surface. Yet, publicly available experiment antibody sequence-binding datasets may not contain the mutagenic, antigenic, or antibody sequence diversity necessary for deep learning approaches to capture molecular recognition. In part, this is because limited experimental platforms exist for assessing quantitative and simultaneous sequence-function relationships for multiple antibodies. Here we present MAGMA-seq, an integrated technology that combines multiple antigens and multiple antibodies and determines quantitative biophysical parameters using deep sequencing. We demonstrate MAGMA-seq on two pooled libraries comprising mutants of nine different human antibodies spanning light chain gene usage, CDR H3 length, and antigenic targets. We demonstrate the comprehensive mapping of potential antibody development pathways, sequence-binding relationships for multiple antibodies simultaneously, and identification of paratope sequence determinants for binding recognition for broadly neutralizing antibodies (bnAbs). MAGMA-seq enables rapid and scalable antibody engineering of multiple lead candidates because it can measure binding for mutants of many given parental antibodies in a single experiment.

Killer yeasts: expanding frontiers in the age of synthetic biology

Killer yeasts secrete protein toxins that are selectively lethal to other yeast and filamentous fungi. These exhibit exceptional genetic and functional diversity, and have several biotechnological applications. However, despite decades of research, several limitations hinder their widespread adoption. In this perspective we contend that technical advances in synthetic biology present an unprecedented opportunity to unlock the full potential of yeast killer systems across a spectrum of applications. By leveraging these new technologies, engineered killer toxins may emerge as a pivotal new tool to address antifungal resistance and food security. Finally, we speculate on the biotechnological potential of re-engineering host double-stranded (ds) RNA mycoviruses, from which many toxins derive, as a safe and noninfectious system to produce designer RNA.

Is Novelty Predictable?

Machine learning-based design has gained traction in the sciences, most notably in the design of small molecules, materials, and proteins, with societal applications ranging from drug development and plastic degradation to carbon sequestration. When designing objects to achieve novel property values with machine learning, one faces a fundamental challenge: how to push past the frontier of current knowledge, distilled from the training data into the model, in a manner that rationally controls the risk of failure. If one trusts learned models too much in extrapolation, one is likely to design rubbish. In contrast, if one does not extrapolate, one cannot find novelty. Herein, we ponder how one might strike a useful balance between these two extremes. We focus in particular on designing proteins with novel property values, although much of our discussion is relevant to machine learning-based design more broadly.

An integrated technology for quantitative wide mutational scanning of human antibody Fab libraries

Antibodies are engineerable quantities in medicine. Learning antibody molecular recognition would enable the in silico design of high affinity binders against nearly any proteinaceous surface. Yet, publicly available experiment antibody sequence-binding datasets may not contain the mutagenic, antigenic, or antibody sequence diversity necessary for deep learning approaches to capture molecular recognition. In part, this is because limited experimental platforms exist for assessing quantitative and simultaneous sequence-function relationships for multiple antibodies. Here we present MAGMA-seq, an integrated technology that combines multiple antigens and multiple antibodies and determines quantitative biophysical parameters using deep sequencing. We demonstrate MAGMA-seq on two pooled libraries comprising mutants of ten different human antibodies spanning light chain gene usage, CDR H3 length, and antigenic targets. We demonstrate the comprehensive mapping of potential antibody development pathways, sequence-binding relationships for multiple antibodies simultaneously, and identification of paratope sequence determinants for binding recognition for broadly neutralizing antibodies (bnAbs). MAGMA-seq enables rapid and scalable antibody engineering of multiple lead candidates because it can measure binding for mutants of many given parental antibodies in a single experiment.

Advances in ligand-specific biosensing for structurally similar molecules

The specificity of biological systems makes it possible to develop biosensors targeting specific metabolites, toxins, and pollutants in complex medical or environmental samples without interference from structurally similar compounds. For the last two decades, great efforts have been devoted to creating proteins or nucleic acids with novel properties through synthetic biology strategies. Beyond augmenting biocatalytic activity, expanding target substrate scopes, and enhancing enzymes' enantioselectivity and stability, an increasing research area is the enhancement of molecular specificity for genetically encoded biosensors. Here, we summarize recent advances in the development of highly specific biosensor systems and their essential applications. First, we describe the rational design principles required to create libraries containing potential mutants with less promiscuity or better specificity. Next, we review the emerging high-throughput screening techniques to engineer biosensing specificity for the desired target. Finally, we examine the computer-aided evaluation and prediction methods to facilitate the construction of ligand-specific biosensors.

MBE: model-based enrichment estimation and prediction for differential sequencing data

Characterizing differences in sequences between two conditions, such as with and without drug exposure, using high-throughput sequencing data is a prevalent problem involving quantifying changes in sequence abundances, and predicting such differences for unobserved sequences. A key shortcoming of current approaches is their extremely limited ability to share information across related but non-identical reads. Consequently, they cannot use sequencing data effectively, nor be directly applied in many settings of interest. We introduce model-based enrichment (MBE) to overcome this shortcoming. We evaluate MBE using both simulated and real data. Overall, MBE improves accuracy compared to current differential analysis methods.

Position-Specific Enrichment Ratio Matrix scores predict antibody variant properties from deep sequencing data

MOTIVATION

Deep sequencing of antibody and related protein libraries after phage or yeast-surface display sorting is widely used to identify variants with increased affinity, specificity, and/or improvements in key biophysical properties. Conventional approaches for identifying optimal variants typically use the frequencies of observation in enriched libraries or the corresponding enrichment ratios. However, these approaches disregard the vast majority of deep sequencing data and often fail to identify the best variants in the libraries.

RESULTS

Here, we present a method, Position-Specific Enrichment Ratio Matrix (PSERM) scoring, that uses entire deep sequencing datasets from pre- and post-selections to score each observed protein variant. The PSERM scores are the sum of the site-specific enrichment ratios observed at each mutated position. We find that PSERM scores are much more reproducible and correlate more strongly with experimentally measured properties than frequencies or enrichment ratios, including for multiple antibody properties (affinity and non-specific binding) for a clinical-stage antibody (emibetuzumab). We expect that this method will be broadly applicable to diverse protein engineering campaigns.

AVAILABILITY AND IMPLEMENTATION

All deep sequencing datasets and code to perform the analyses presented within are available via https://github.com/Tessier-Lab-UMich/PSERM_paper.

Position-Specific Enrichment Ratio Matrix scores predict antibody variant properties from deep sequencing data

Motivation

Deep sequencing of antibody and related protein libraries after phage or yeast-surface display sorting is widely used to identify variants with increased affinity, specificity and/or improvements in key biophysical properties. Conventional approaches for identifying optimal variants typically use the frequencies of observation in enriched libraries or the corresponding enrichment ratios. However, these approaches disregard the vast majority of deep sequencing data and often fail to identify the best variants in the libraries.

Results

Here, we present a method, Position-Specific Enrichment Ratio Matrix (PSERM) scoring, that uses entire deep sequencing datasets from pre- and post-selections to score each observed protein variant. The PSERM scores are the sum of the site-specific enrichment ratios observed at each mutated position. We find that PSERM scores are much more reproducible and correlate more strongly with experimentally measured properties than frequencies or enrichment ratios, including for multiple antibody properties (affinity and non-specific binding) for a clinical-stage antibody (emibetuzumab). We expect that this method will be broadly applicable to diverse protein engineering campaigns.

Availability

All deep sequencing datasets and code to do the analyses presented within are available via GitHub.

Contact

Peter Tessier, ptessier@umich.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Designed loop extension followed by combinatorial screening confers high specificity to a broad matrix metalloproteinaseinhibitor

Matrix metalloproteinases (MMPs) are key drivers of various diseases, including cancer. Development of selective probes and drugs capable of selectively inhibiting the individual members of the large MMP family remains a persistent challenge. The inhibitory N-terminal domain of tissue inhibitor of metalloproteinases-2 (N-TIMP2), a natural broad MMP inhibitor, can provide a scaffold for protein engineering to create more selective MMP inhibitors. Here, we pursued a unique approach harnessing both computational design and combinatorial screening to confer high binding specificity toward a target MMP in preference to an anti-target MMP. We designed a loop extension of N-TIMP2 to allow new interactions with the non-conserved MMP surface and generated an efficient focused library for yeast surface display, which was then screened for high binding to the target MMP-14 and low binding to anti-target MMP-3. Deep sequencing analysis identified the most promising variants, which were expressed, purified, and tested for selectivity of inhibition. Our best N-TIMP2 variant exhibited 29 pM binding affinity to MMP-14 and 2.4 µM affinity to MMP-3, revealing 7500-fold greater specificity than WT N-TIMP2. High-confidence structural models were obtained by including NGS data in the AlphaFold multiple sequence alignment. The modeling together with experimental mutagenesis validated our design predictions, demonstrating that the loop extension packs tightly against non-conserved residues on MMP-14 and clashes with MMP-3. This study demonstrates how introduction of loop extensions in a manner guided by target protein conservation data and loop design can offer an attractive strategy to achieve specificity in design of protein ligands.

Data-Driven Synthetic Cell Factories Development for Industrial Biomanufacturing

Revolutionary breakthroughs in artificial intelligence (AI) and machine learning (ML) have had a profound impact on a wide range of scientific disciplines, including the development of artificial cell factories for biomanufacturing. In this paper, we review the latest studies on the application of data-driven methods for the design of new proteins, pathways, and strains. We first briefly introduce the various types of data and databases relevant to industrial biomanufacturing, which are the basis for data-driven research. Different types of algorithms, including traditional ML and more recent deep learning methods, are also presented. We then demonstrate how these data-based approaches can be applied to address various issues in cell factory development using examples from recent studies, including the prediction of protein function, improvement of metabolic models, and estimation of missing kinetic parameters, design of non-natural biosynthesis pathways, and pathway optimization. In the last section, we discuss the current limitations of these data-driven approaches and propose that data-driven methods should be integrated with mechanistic models to complement each other and facilitate the development of synthetic strains for industrial biomanufacturing.

Identification of stabilizing point mutations through mutagenesis of destabilized protein libraries

Although there have been recent transformative advances in the area of protein structure prediction, prediction of point mutations that improve protein stability remains challenging. It is possible to construct and screen large mutant libraries for improved activity or ligand binding. However, reliable screens for mutants that improve protein stability do not yet exist, especially for proteins that are well folded and relatively stable. Here, we demonstrate that incorporation of a single, specific, destabilizing mutation termed parent inactivating mutation into each member of a single-site saturation mutagenesis library, followed by screening for suppressors, allows for robust and accurate identification of stabilizing mutations. We carried out fluorescence-activated cell sorting of such a yeast surface display, saturation suppressor library of the bacterial toxin CcdB, followed by deep sequencing of sorted populations. We found that multiple stabilizing mutations could be identified after a single round of sorting. In addition, multiple libraries with different parent inactivating mutations could be pooled and simultaneously screened to further enhance the accuracy of identification of stabilizing mutations. Finally, we show that individual stabilizing mutations could be combined to result in a multi-mutant that demonstrated an increase in thermal melting temperature of about 20 °C, and that displayed enhanced tolerance to high temperature exposure. We conclude that as this method is robust and employs small library sizes, it can be readily extended to other display and screening formats to rapidly isolate stabilized protein mutants.

In silico evolution of nucleic acid-binding proteins from a nonfunctional scaffold

Directed evolution emulates the process of natural selection to produce proteins with improved or altered functions. These approaches have proven to be very powerful but are technically challenging and particularly time and resource intensive. To bypass these limitations, we constructed a system to perform the entire process of directed evolution in silico. We employed iterative computational cycles of mutation and evaluation to predict mutations that confer high-affinity binding activities for DNA and RNA to an initial de novo designed protein with no inherent function. Beneficial mutations revealed modes of nucleic acid recognition not previously observed in natural proteins, highlighting the ability of computational directed evolution to access new molecular functions. Furthermore, the process by which new functions were obtained closely resembles natural evolution and can provide insights into the contributions of mutation rate, population size and selective pressure on functionalization of macromolecules in nature.

Engineering of the Ligand Specificity of Transcriptional Regulator XylS by Deep Mutational Scanning

Deep mutational scanning is a method for protein engineering. Here, we applied it to alter the ligand specificity of the transcriptional regulator XylS from Pseudomonas putida to recognize p-toluic acid instead of the native ligand m-toluic acid. For this purpose, we used an antibiotic resistance gene-based dual screening system, which was constructed for the directed evolution of XylS toward the above-mentioned ligand specificity. We constructed a xylS mutant library in which each codon for the amino acid residue of the putative ligand-binding domain (residues 1-213, except 7th residue) was randomized to generate all possible single amino acid-substituted XylS variants and introduced it into Escherichia coli harboring the selection plasmid for the screening system. The cells were cultured in the presence of appropriate antibiotics and m-toluic acid or p-toluic acid, and the frequency of each mutation present in the library was examined using a next-generation sequencer before and after cultivation. Heatmaps showing the enrichment score of each XylS variant were obtained. By searching for a p-toluic-acid-specific heatmap pattern, we focused on G71 and H77. Analysis of the ligand specificities of G71- or H77-substituted XylS variants revealed that several G71-substituted XylS variants responded specifically to p-toluic acid. Thus, the 71st residue was found to be an unprecedented residue that is important for switching ligand specificity. Our study demonstrated the usefulness of deep mutational scanning in engineering the ligand specificity of a transcriptional regulator without structural information. We also discussed the advantages and disadvantages of deep mutational scanning compared with directed evolution.

Machine Learning-driven Protein Library Design: A Path Toward Smarter Libraries

Proteins are small yet valuable biomolecules that play a versatile role in therapeutics and diagnostics. The intricate sequence-structure-function paradigm in the realm of proteins opens the possibility for directly mapping amino acid sequence to function. However, the rugged nature of the protein fitness landscape and an astronomical number of possible mutations even for small proteins make navigating this system a daunting task. Moreover, the scarcity of functional proteins and the ease with which deleterious mutations are introduced, due to complex epistatic relationships, compound the existing challenges. This highlights the need for auxiliary tools in current techniques such as rational design and directed evolution. To that end, the state-of-the-art machine learning can offer time and cost efficiency in finding high fitness proteins, circumventing unnecessary wet-lab experiments. In the context of improving library design, machine learning provides valuable insights via its unique features such as high adaptation to complex systems, multi-tasking, and parallelism, and the ability to capture hidden trends in input data. Finally, both the advancements in computational resources and the rapidly increasing number of sequences in protein databases will allow more promising and detailed insights delivered from machine learning to protein library design. In this chapter, fundamental concepts and a method for machine learning-driven library design leveraging deep sequencing datasets will be discussed. We elaborate on (1) basic knowledge about machine learning algorithms, (2) the benefit of machine learning in library design, and (3) methodology for implementing machine learning in library design.

Yeast Surface Display: New Opportunities for a Time-Tested Protein Engineering System

Yeast surface display has proven to be a powerful tool for the discovery of antibodies and other novel binding proteins and for engineering the affinity and selectivity of existing proteins for their targets. In the decades since the first demonstrations of the approach, the range of yeast display applications has greatly expanded to include many different protein targets and has grown to encompass methods for rapid protein characterization. Here, we briefly summarize the development of yeast display methodologies and highlight several selected examples of recent applications to timely and challenging protein engineering and characterization problems.

Heterogeneity of the GFP fitness landscape and data-driven protein design

Studies of protein fitness landscapes reveal biophysical constraints guiding protein evolution and empower prediction of functional proteins. However, generalisation of these findings is limited due to scarceness of systematic data on fitness landscapes of proteins with a defined evolutionary relationship. We characterized the fitness peaks of four orthologous fluorescent proteins with a broad range of sequence divergence. While two of the four studied fitness peaks were sharp, the other two were considerably flatter, being almost entirely free of epistatic interactions. Mutationally robust proteins, characterized by a flat fitness peak, were not optimal templates for machine-learning-driven protein design - instead, predictions were more accurate for fragile proteins with epistatic landscapes. Our work paves insights for practical application of fitness landscape heterogeneity in protein engineering.

Oligo Pools as an Affordable Source of Synthetic DNA for Cost-Effective Library Construction in Protein- and Metabolic Pathway Engineering

The construction of custom libraries is critical for rational protein engineering and directed evolution. Array‐synthesized oligo pools of thousands of user‐defined sequences (up to ∼350 bases in length) have emerged as a low‐cost commercially available source of DNA. These pools cost ≤10 % (depending on error rate and length) of other commercial sources of custom DNA, and this significant cost difference can determine whether an enzyme engineering project can be realized on a given research budget. However, while being cheap, oligo pools do suffer from a low concentration of individual oligos and relatively high error rates. Several powerful techniques that specifically make use of oligo pools have been developed and proven valuable or even essential for next‐generation protein and pathway engineering strategies, such as sequence‐function mapping, enzyme minimization, or de‐novo design. Here we consolidate the knowledge on these techniques and their applications to facilitate the use of oligo pools within the protein engineering community.

Climbing Up and Down Binding Landscapes through Deep Mutational Scanning of Three Homologous Protein-Protein Complexes

Protein-protein interactions (PPIs) have evolved to display binding affinities that can support their function. As such, cognate and noncognate PPIs could be highly similar structurally but exhibit huge differences in binding affinities. To understand this phenomenon, we study three homologous protease-inhibitor PPIs that span 9 orders of magnitude in binding affinity. Using state-of-the-art methodology that combines protein randomization, affinity sorting, deep sequencing, and data normalization, we report quantitative binding landscapes consisting of ΔΔG_bind values for the three PPIs, gleaned from tens of thousands of single and double mutations. We show that binding landscapes of the three complexes are strikingly different and depend on the PPI evolutionary optimality. We observe different patterns of couplings between mutations for the three PPIs with negative and positive epistasis appearing most frequently at hot-spot and cold-spot positions, respectively. The evolutionary trends observed here are likely to be universal to other biological complexes in the cell.

PacBio sequencing output increased through uniform and directional fivefold concatenation

Advances in sequencing technology have allowed researchers to sequence DNA with greater ease and at decreasing costs. Main developments have focused on either sequencing many short sequences or fewer large sequences. Methods for sequencing mid-sized sequences of 600-5,000 bp are currently less efficient. For example, the PacBio Sequel I system yields ~ 100,000-300,000 reads with an accuracy per base pair of 90-99%. We sought to sequence several DNA populations of ~ 870 bp in length with a sequencing accuracy of 99% and to the greatest depth possible. We optimised a simple, robust method to concatenate genes of ~ 870 bp five times and then sequenced the resulting DNA of ~ 5,000 bp by PacBioSMRT long-read sequencing. Our method improved upon previously published concatenation attempts, leading to a greater sequencing depth, high-quality reads and limited sample preparation at little expense. We applied this efficient concatenation protocol to sequence nine DNA populations from a protein engineering study. The improved method is accompanied by a simple and user-friendly analysis pipeline, DeCatCounter, to sequence medium-length sequences efficiently at one-fifth of the cost.

Applications of Deep Mutational Scanning in Virology

Several recently developed high-throughput techniques have changed the field of molecular virology. For example, proteomics studies reveal complete interactomes of a viral protein, genome-wide CRISPR knockout and activation screens probe the importance of every single human gene in aiding or fighting a virus, and ChIP-seq experiments reveal genome-wide epigenetic changes in response to infection. Deep mutational scanning is a relatively novel form of protein science which allows the in-depth functional analysis of every nucleotide within a viral gene or genome, revealing regions of importance, flexibility, and mutational potential. In this review, we discuss the application of this technique to RNA viruses including members of the Flaviviridae family, Influenza A Virus and Severe Acute Respiratory Syndrome Coronavirus 2. We also briefly discuss the reverse genetics systems which allow for analysis of viral replication cycles, next-generation sequencing technologies and the bioinformatics tools that facilitate this research.

Substitution of a Surface-Exposed Residue Involved in an Allosteric Network Enhances Tryptophan Synthase Function in Cells

Networks of noncovalent amino acid interactions propagate allosteric signals throughout proteins. Tryptophan synthase (TS) is an allosterically controlled bienzyme in which the indole product of the alpha subunit (αTS) is transferred through a 25 Å hydrophobic tunnel to the active site of the beta subunit (βTS). Previous nuclear magnetic resonance and molecular dynamics simulations identified allosteric networks in αTS important for its function. We show here that substitution of a distant, surface-exposed network residue in αTS enhances tryptophan production, not by activating αTS function, but through dynamically controlling the opening of the indole channel and stimulating βTS activity. While stimulation is modest, the substitution also enhances cell growth in a tryptophan-auxotrophic strain of Escherichia coli compared to complementation with wild-type αTS, emphasizing the biological importance of the network. Surface-exposed networks provide new opportunities in allosteric drug design and protein engineering, and hint at potential information conduits through which the functions of a metabolon or even larger proteome might be coordinated and regulated.

Deep Mutational Scanning of Viral Glycoproteins and Their Host Receptors

Deep mutational scanning or deep mutagenesis is a powerful tool for understanding the sequence diversity available to viruses for adaptation in a laboratory setting. It generally involves tracking an in vitro selection of protein sequence variants with deep sequencing to map mutational effects based on changes in sequence abundance. Coupled with any of a number of selection strategies, deep mutagenesis can explore the mutational diversity available to viral glycoproteins, which mediate critical roles in cell entry and are exposed to the humoral arm of the host immune response. Mutational landscapes of viral glycoproteins for host cell attachment and membrane fusion reveal extensive epistasis and potential escape mutations to neutralizing antibodies or other therapeutics, as well as aiding in the design of optimized immunogens for eliciting broadly protective immunity. While less explored, deep mutational scans of host receptors further assist in understanding virus-host protein interactions. Critical residues on the host receptors for engaging with viral spikes are readily identified and may help with structural modeling. Furthermore, mutations may be found for engineering soluble decoy receptors as neutralizing agents that specifically bind viral targets with tight affinity and limited potential for viral escape. By untangling the complexities of how sequence contributes to viral glycoprotein and host receptor interactions, deep mutational scanning is impacting ideas and strategies at multiple levels for combatting circulating and emergent virus strains.

Computational design and experimental optimization of protein binders with prospects for biomedical applications

Protein-based binders have become increasingly more attractive candidates for drug and imaging agent development. Such binders could be evolved from a number of different scaffolds, including antibodies, natural protein effectors and unrelated small protein domains of different geometries. While both computational and experimental approaches could be utilized for protein binder engineering, in this review we focus on various computational approaches for protein binder design and demonstrate how experimental selection could be applied to subsequently optimize computationally-designed molecules. Recent studies report a number of designed protein binders with pM affinities and high specificities for their targets. These binders usually characterized with high stability, solubility, and low production cost. Such attractive molecules are bound to become more common in various biotechnological and biomedical applications in the near future.

Accurate Broadcasting of Substrate Fitness for Lactazole Biosynthetic Pathway from Reactivity-Profiling mRNA Display

We report a method for the high-throughput reactivity profiling of genetically encoded libraries as a tool to study substrate fitness landscapes for RiPP (ribosomally synthesized and post-translationally modified peptide) biosynthetic enzymes. This method allowed us to rapidly analyze the substrate preferences of the lactazole biosynthetic pathway using a saturation mutagenesis mRNA display library of lactazole precursor peptides. We demonstrate that the assay produces accurate and reproducible in vitro data, enabling the quantification of reaction yields with temporal resolution. Our results recapitulate the previously established knowledge on lactazole biosynthesis and expand it by identifying the extent of substrate promiscuity exhibited by the enzymes. This work lays a foundation for the construction and screening of mRNA display-based combinatorial thiopeptide libraries for the discovery of lactazole-inspired thiopeptides with de novo designed biological activities.

Inferring Protein Sequence-Function Relationships with Large-Scale Positive-Unlabeled Learning

Machine learning can infer how protein sequence maps to function without requiring a detailed understanding of the underlying physical or biological mechanisms. It is challenging to apply existing supervised learning frameworks to large-scale experimental data generated by deep mutational scanning (DMS) and related methods. DMS data often contain high-dimensional and correlated sequence variables, experimental sampling error and bias, and the presence of missing data. Notably, most DMS data do not contain examples of negative sequences, making it challenging to directly estimate how sequence affects function. Here, we develop a positive-unlabeled (PU) learning framework to infer sequence-function relationships from large-scale DMS data. Our PU learning method displays excellent predictive performance across ten large-scale sequence-function datasets, representing proteins of different folds, functions, and library types. The estimated parameters pinpoint key residues that dictate protein structure and function. Finally, we apply our statistical sequence-function model to design highly stabilized enzymes.

Structure-Guided Molecular Engineering of a Vascular Endothelial Growth Factor Antagonist to Treat Retinal Diseases

Background

Ocular neovascularization is a hallmark of retinal diseases including neovascular age-related macular degeneration and diabetic retinopathy, two leading causes of blindness in adults. Neovascularization is driven by the interaction of soluble vascular endothelial growth factor (VEGF) ligands with transmembrane VEGF receptors (VEGFR), and inhibition of the VEGF pathway has shown tremendous clinical promise. However, anti-VEGF therapies require invasive intravitreal injections at frequent intervals and high doses, and many patients show incomplete responses to current drugs due to the lack of sustained VEGF signaling suppression.

Methods

We synthesized insights from structural biology with molecular engineering technologies to engineer an anti-VEGF antagonist protein. Starting from the clinically approved decoy receptor protein aflibercept, we strategically designed a yeast-displayed mutagenic library of variants and isolated clones with superior VEGF affinity compared to the clinical drug. Our lead engineered protein was expressed in the choroidal space of rat eyes via nonviral gene delivery.

Results

Using a structure-informed directed evolution approach, we identified multiple promising anti-VEGF antagonist proteins with improved target affinity. Improvements were primarily mediated through reduction in dissociation rate, and structurally significant convergent sequence mutations were identified. Nonviral gene transfer of our engineered antagonist protein demonstrated robust and durable expression in the choroid of treated rats one month post-injection.

Conclusions

We engineered a novel anti-VEGF protein as a new weapon against retinal diseases and demonstrated safe and noninvasive ocular delivery in rats. Furthermore, our structure-guided design approach presents a general strategy for discovery of targeted protein drugs for a vast array of applications.

Construction and analysis of the protein-protein interaction network for the olfactory system of the silkworm Bombyx mori

Olfaction plays an essential role in feeding and information exchange in insects. Previous studies on the olfaction of silkworms have provided a wealth of information about genes and proteins, yet, most studies have only focused on a single gene or protein related to the insect's olfaction. The aim of the current study is to determine key proteins in the olfactory system of the silkworm, and further understand protein-protein interactions (PPIs) in the olfactory system of Lepidoptera. To achieve this goal, we integrated Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and network analyses. Furthermore, we selected 585 olfactory-related proteins and constructed a (PPI) network for the olfactory system of the silkworm. Network analysis led to the identification of several key proteins, including GSTz1, LOC733095, BGIBMGA002169-TA, BGIBMGA010939-TA, GSTs2, GSTd2, Or-2, and BGIBMGA013255-TA. A comprehensive evaluation of the proteins showed that glutathione S-transferases (GSTs) had the highest ranking. GSTs also had the highest enrichment levels in GO and KEGG. In conclusion, our analysis showed that key nodes in the biological network had a significant impact on the network, and the key proteins identified via network analysis could serve as new research targets to determine their functions in olfaction.

Measurements drive progress in directed evolution for precise engineering of biological systems

Precise engineering of biological systems requires quantitative, high-throughput measurements, exemplified by progress in directed evolution. New approaches allow high-throughput measurements of phenotypes and their corresponding genotypes. When integrated into directed evolution, these quantitative approaches enable the precise engineering of biological function. At the same time, the increasingly routine availability of large, high-quality data sets supports the integration of machine learning with directed evolution. Together, these advances herald striking capabilities for engineering biology.

Discovery, Design, and Structural Characterization of Alkane-Producing Enzymes across the Ferritin-like Superfamily

To complement established rational and evolutionary protein design approaches, significant efforts are being made to utilize computational modeling and the diversity of naturally occurring protein sequences. Here, we combine structural biology, genomic mining, and computational modeling to identify structural features critical to aldehyde deformylating oxygenases (ADOs), an enzyme family that has significant implications in synthetic biology and chemoenzymatic synthesis. Through these efforts, we discovered latent ADO-like function across the ferritin-like superfamily in various species of Bacteria and Archaea. We created a machine learning model that uses protein structural features to discriminate ADO-like activity. Computational enzyme design tools were then utilized to introduce ADO-like activity into the small subunit of Escherichia coli class I ribonucleotide reductase. The integrated approach of genomic mining, structural biology, molecular modeling, and machine learning has the potential to be utilized for rapid discovery and modulation of functions across enzyme families.

Practically useful protein-design methods combining phylogenetic and atomistic calculations

Our ability to design new or improved biomolecular activities depends on understanding the sequence-function relationships in proteins. The large size and fold complexity of most proteins, however, obscure these relationships, and protein-optimization methods continue to rely on laborious experimental iterations. Recently, a deeper understanding of the roles of stability-threshold effects and biomolecular epistasis in proteins has led to the development of hybrid methods that combine phylogenetic analysis with atomistic design calculations. These methods enable reliable and even single-step optimization of protein stability, expressibility, and activity in proteins that were considered outside the scope of computational design. Furthermore, ancestral-sequence reconstruction produces insights on missing links in the evolution of enzymes and binders that may be used in protein design. Through the combination of phylogenetic and atomistic calculations, the long-standing goal of general computational methods that can be universally applied to study and optimize proteins finally seems within reach.

Quantitative mapping of binding specificity landscapes for homologous targets by using a high-throughput method

To facilitate investigations of protein-protein interactions (PPIs), we developed a novel platform for quantitative mapping of protein binding specificity landscapes, which combines the multi-target screening of a mutagenesis library into high- and low-affinity populations with sophisticated next-generation sequencing analysis. Importantly, this method generates accurate models to predict affinity and specificity values for any mutation within a protein complex, and requires only a few experimental binding affinity measurements using purified proteins for calibration. We demonstrated the utility of the approach by mapping quantitative landscapes for interactions between the N-terminal domain of the tissue inhibitor of metalloproteinase 2 (N-TIMP2) and three matrix metalloproteinases (MMPs) having homologous structures but different affinities (MMP-1, MMP-3, and MMP-14). The binding landscapes for N-TIMP2/MMP-1 and N-TIMP2/MMP-3 showed the PPIs to be almost fully optimized, with most single mutations giving a loss of affinity. In contrast, the non-optimized PPI for N-TIMP2/MMP-14 was reflected in a wide range of binding affinities, where single mutations exhibited a far more attenuated effect on the PPI. Our new platform reliably and comprehensively identified not only hot- and cold-spot residues, but also specificity-switch mutations that shape target affinity and specificity. Thus, our approach provides a methodology giving an unprecedentedly rich quantitative analysis of the binding specificity landscape, which will broaden the understanding of the mechanisms and evolutionary origins of specific PPIs and facilitate the rational design of specific inhibitors for structurally similar target proteins.

Variant effect predictions capture some aspects of deep mutational scanning experiments

BACKGROUND

Deep mutational scanning (DMS) studies exploit the mutational landscape of sequence variation by systematically and comprehensively assaying the effect of single amino acid variants (SAVs; also referred to as missense mutations, or non-synonymous Single Nucleotide Variants - missense SNVs or nsSNVs) for particular proteins. We assembled SAV annotations from 22 different DMS experiments and normalized the effect scores to evaluate variant effect prediction methods. Three trained on traditional variant effect data (PolyPhen-2, SIFT, SNAP2), a regression method optimized on DMS data (Envision), and a naïve prediction using conservation information from homologs.

RESULTS

On a set of 32,981 SAVs, all methods captured some aspects of the experimental effect scores, albeit not the same. Traditional methods such as SNAP2 correlated slightly more with measurements and better classified binary states (effect or neutral). Envision appeared to better estimate the precise degree of effect. Most surprising was that the simple naïve conservation approach using PSI-BLAST in many cases outperformed other methods. All methods captured beneficial effects (gain-of-function) significantly worse than deleterious (loss-of-function). For the few proteins with multiple independent experimental measurements, experiments differed substantially, but agreed more with each other than with predictions.

CONCLUSIONS

DMS provides a new powerful experimental means of understanding the dynamics of the protein sequence space. As always, promising new beginnings have to overcome challenges. While our results demonstrated that DMS will be crucial to improve variant effect prediction methods, data diversity hindered simplification and generalization.

Active Site-Induced Evolutionary Constraints Follow Fold Polarity Principles in Soluble Globular Enzymes

A recent analysis of evolutionary rates in >500 globular soluble enzymes revealed pervasive conservation gradients toward catalytic residues. By looking at amino acid preference profiles rather than evolutionary rates in the same data set, we quantified the effects of active sites on site-specific constraints for physicochemical traits. We found that conservation gradients respond to constraints for polarity, hydrophobicity, flexibility, rigidity and structure in ways consistent with fold polarity principles; while sites far from active sites seem to experience no physicochemical constraint, rather being highly variable and favoring amino acids of low metabolic cost. Globally, our results highlight that amino acid variation contains finer information about protein structure than usually regarded in evolutionary models, and that this information is retrievable automatically with simple fits. We propose that analyses of the kind presented here incorporated into models of protein evolution should allow for better description of the physical chemistry that underlies molecular evolution.

Generating quantitative binding landscapes through fractional binding selections combined with deep sequencing and data normalization

Quantifying the effects of various mutations on binding free energy is crucial for understanding the evolution of protein-protein interactions and would greatly facilitate protein engineering studies. Yet, measuring changes in binding free energy (ΔΔG_bind) remains a tedious task that requires expression of each mutant, its purification, and affinity measurements. We developed an attractive approach that allows us to quantify ΔΔG_bind for thousands of protein mutants in one experiment. Our protocol combines protein randomization, Yeast Surface Display technology, deep sequencing, and a few experimental ΔΔG_bind data points on purified proteins to generate ΔΔG_bind values for the remaining numerous mutants of the same protein complex. Using this methodology, we comprehensively map the single-mutant binding landscape of one of the highest-affinity interaction between BPTI and Bovine Trypsin (BT). We show that ΔΔG_bind for this interaction could be quantified with high accuracy over the range of 12 kcal mol⁻¹ displayed by various BPTI single mutants.

Quantifying the effect of mutations on binding free energy is important to understand protein-protein interaction (PPI). Here the authors develop a method based on yeast display and next-generation sequencing to generate quantitative binding landscapes for any PPI regardless of their Kd value.

Protein tolerance to random circular permutation correlates with thermostability and local energetics of residue-residue contacts

Adenylate kinase (AK) orthologs with a range of thermostabilities were subjected to random circular permutation, and deep mutational scanning was used to evaluate where new protein termini were nondisruptive to activity. The fraction of circularly permuted variants that retained function in each library correlated with AK thermostability. In addition, analysis of the positional tolerance to new termini, which increase local conformational flexibility, showed that bonds were either functionally sensitive to cleavage across all homologs, differentially sensitive, or uniformly tolerant. The mobile AMP-binding domain, which displays the highest calculated contact energies, presented the greatest tolerance to new termini across all AKs. In contrast, retention of function in the lid and core domains was more dependent upon AK melting temperature. These results show that family permutation profiling identifies primary structure that has been selected by evolution for dynamics that are critical to activity within an enzyme family. These findings also illustrate how deep mutational scanning can be applied to protein homologs in parallel to differentiate how topology, stability, and local energetics govern mutational tolerance.

Current Trends in Protein Engineering: Updates and Progress

Proteins are one of the most important and resourceful biomolecules that find applications in health, industry, medicine, research, and biotechnology. Given its tremendous relevance, protein engineering has emerged as significant biotechnological intervention in this area. Strategic utilization of protein engineering methods and approaches has enabled better enzymatic properties, better stability, increased catalytic activity and most importantly, interesting and wide range applicability of proteins. In fact, the commercialization of engineered proteins have manifested in economically beneficial and viable solutions for industry and healthcare sector. Protein engineering has also evolved to become a powerful tool contributing significantly to the developments in both synthetic biology and metabolic engineering. The present review revisits the current trends in protein engineering approaches such as rational design, directed evolution, de novo design, computational approaches etc. and encompasses the recent progresses made in this field over the last few years. The review also throws light on advanced or futuristic protein engineering aspects, which are being explored for design and development of novel proteins with improved properties or advanced applications.

Advances in protein structure prediction and design

The prediction of protein three-dimensional structure from amino acid sequence has been a grand challenge problem in computational biophysics for decades, owing to its intrinsic scientific interest and also to the many potential applications for robust protein structure prediction algorithms, from genome interpretation to protein function prediction. More recently, the inverse problem - designing an amino acid sequence that will fold into a specified three-dimensional structure - has attracted growing attention as a potential route to the rational engineering of proteins with functions useful in biotechnology and medicine. Methods for the prediction and design of protein structures have advanced dramatically in the past decade. Increases in computing power and the rapid growth in protein sequence and structure databases have fuelled the development of new data-intensive and computationally demanding approaches for structure prediction. New algorithms for designing protein folds and protein-protein interfaces have been used to engineer novel high-order assemblies and to design from scratch fluorescent proteins with novel or enhanced properties, as well as signalling proteins with therapeutic potential. In this Review, we describe current approaches for protein structure prediction and design and highlight a selection of the successful applications they have enabled.

Tools and systems for evolutionary engineering of biomolecules and microorganisms

Evolutionary approaches have been providing solutions to various bioengineering challenges in an efficient manner. In addition to traditional adaptive laboratory evolution and directed evolution, recent advances in synthetic biology and fluidic systems have opened a new era of evolutionary engineering. Synthetic genetic circuits have been created to control mutagenesis and enable screening of various phenotypes, particularly metabolite production. Fluidic systems can be used for high-throughput screening and multiplexed continuous cultivation of microorganisms. Moreover, continuous directed evolution has been achieved by combining all the steps of evolutionary engineering. Overall, modern tools and systems for evolutionary engineering can be used to establish the artificial equivalent to natural evolution for various research applications.

Adaptive Landscapes in the Age of Synthetic Biology

For nearly a century adaptive landscapes have provided overviews of the evolutionary process and yet they remain metaphors. We redefine adaptive landscapes in terms of biological processes rather than descriptive phenomenology. We focus on the underlying mechanisms that generate emergent properties such as epistasis, dominance, trade-offs and adaptive peaks. We illustrate the utility of landscapes in predicting the course of adaptation and the distribution of fitness effects. We abandon aged arguments concerning landscape ruggedness in favor of empirically determining landscape architecture. In so doing, we transform the landscape metaphor into a scientific framework within which causal hypotheses can be tested.

Biophysical and Mechanistic Models for Disease-Causing Protein Variants

The rapid decrease in DNA sequencing cost is revolutionizing medicine and science. In medicine, genome sequencing has revealed millions of missense variants that change protein sequences, yet we only understand the molecular and phenotypic consequences of a small fraction. Within protein science, high-throughput deep mutational scanning experiments enable us to probe thousands of variants in a single, multiplexed experiment. We review efforts that bring together these topics via experimental and computational approaches to determine the consequences of missense variants in proteins. We focus on the role of changes in protein stability as a driver for disease, and how experiments, biophysical models, and computation are providing a framework for understanding and predicting how changes in protein sequence affect cellular protein stability.

Evolution on the Biophysical Fitness Landscape of an RNA Virus

Viral evolutionary pathways are determined by the fitness landscape, which maps viral genotype to fitness. However, a quantitative description of the landscape and the evolutionary forces on it remain elusive. Here, we apply a biophysical fitness model based on capsid folding stability and antibody binding affinity to predict the evolutionary pathway of norovirus escaping a neutralizing antibody. The model is validated by experimental evolution in bulk culture and in a drop-based microfluidics that propagates millions of independent small viral subpopulations. We demonstrate that along the axis of binding affinity, selection for escape variants and drift due to random mutations have the same direction, an atypical case in evolution. However, along folding stability, selection and drift are opposing forces whose balance is tuned by viral population size. Our results demonstrate that predictable epistatic tradeoffs between molecular traits of viral proteins shape viral evolution.

Plasticity and Constraints of tRNA Aminoacylation Define Directed Evolution of Aminoacyl-tRNA Synthetases

Genetic incorporation of noncanonical amino acids (ncAAs) has become a powerful tool to enhance existing functions or introduce new ones into proteins through expanded chemistry. This technology relies on the process of nonsense suppression, which is made possible by directing aminoacyl-tRNA synthetases (aaRSs) to attach an ncAA onto a cognate suppressor tRNA. However, different mechanisms govern aaRS specificity toward its natural amino acid (AA) substrate and hinder the engineering of aaRSs for applications beyond the incorporation of a single l-α-AA. Directed evolution of aaRSs therefore faces two interlinked challenges: the removal of the affinity for cognate AA and improvement of ncAA acylation. Here we review aspects of AA recognition that directly influence the feasibility and success of aaRS engineering toward d- and β-AAs incorporation into proteins in vivo. Emerging directed evolution methods are described and evaluated on the basis of aaRS active site plasticity and its inherent constraints.

A Yeast System for Discovering Optogenetic Inhibitors of Eukaryotic Translation Initiation

The precise spatiotemporal regulation of protein synthesis is essential for many complex biological processes such as memory formation, embryonic development, and tumor formation. Current methods used to study protein synthesis offer only a limited degree of spatiotemporal control. Optogenetic methods, in contrast, offer the prospect of controlling protein synthesis noninvasively within minutes and with a spatial scale as small as a single synapse. Here, we present a hybrid yeast system where growth depends on the activity of human eukaryotic initiation factor 4E (eIF4E) that is suitable for screening optogenetic designs for the down-regulation of protein synthesis. We used this system to screen a diverse initial panel of 15 constructs designed to couple a light switchable domain (PYP, RsLOV, AsLOV, Dronpa) to 4EBP2 (eukaryotic initiation factor 4E binding protein 2), a native inhibitor of translation initiation. We identified cLIPS1 (circularly permuted LOV inhibitor of protein synthesis 1), a fusion of a segment of 4EBP2 and a circularly permuted version of the LOV2 domain from Avena sativa, as a photoactivated inhibitor of translation. Adapting the screen for higher throughput, we tested small libraries of cLIPS1 variants and found cLIPS2, a construct with an improved degree of optical control. We show that these constructs can both inhibit translation in yeast harboring a human eIF4E in vivo, and bind human eIF4E in vitro in a light-dependent manner. This hybrid yeast system thus provides a convenient way for discovering optogenetic constructs that can regulate human eIF4E-dependent translation initiation in a mechanistically defined manner.

An Automated Data-Driven Pipeline for Improving Heterologous Enzyme Expression

Enzymes are the ultimate entities responsible for chemical transformations in natural and engineered biosynthetic pathways. However, many natural enzymes suffer from suboptimal functional expression due to poor intrinsic protein stability. Further, stability enhancing mutations often come at the cost of impaired function. Here we demonstrate an automated protein engineering strategy for stabilizing enzymes while retaining catalytic function using deep mutational scanning coupled to multiple-filter based screening and combinatorial mutagenesis. We validated this strategy by improving the functional expression of a Type III polyketide synthase from the Atropa belladonna biosynthetic pathway for tropane alkaloids. The best variant had a total of 8 mutations with over 25-fold improved activity over wild-type in E. coli cell lysates, an improved melting temperature of 11.5 ± 0.6 °C, and only minimal reduction in catalytic efficiency. We show that the multiple-filter approach maintains acceptable sensitivity with homology modeling structures up to 4 Å RMS. Our results highlight an automated protein engineering tool for improving the stability and solubility of difficult to express enzymes, which has impact for biotechnological applications.