Structure-based design of combinatorial mutagenesis libraries

Problem of Domain/Building Block Preservation in the Evolution of Biological Macromolecules and Evolutionary Computation

Structurally and functionally isolated domains in biological macromolecular evolution, both natural and artificial, are largely similar to "schemata", building blocks (BBs), in evolutionary computation (EC). The problem of preserving in subsequent evolutionary searches the already found domains / BBs is well known and quite relevant in biology as well as in EC. Both biology and EC are seeing parallel and independent development of several approaches to identifying and preserving previously identified domains / BBs. First, we notice the similarity of DNA shuffling methods in synthetic biology and multi-parent recombination algorithms in EC. Furthermore, approaches to computer identification of domains in proteins that are being developed in biology can be aligned with BB identification methods in EC. Finally, approaches to chimeric protein libraries optimization in biology can be compared to evolutionary search methods based on probabilistic models in EC. We propose to validate the prospects of mutual exchange of ideas and transfer of algorithms and approaches between evolutionary systems biology and EC in these three principal directions. A crucial aim of this transfer is the design of new advanced experimental techniques capable of solving more complex problems of in vitro evolution.

Functional Deimmunization of Botulinum Neurotoxin Protease Domain via Computationally Driven Library Design and Ultrahigh-Throughput Screening

Botulinum neurotoxin serotype A (BoNT/A) is a widely used cosmetic agent that also has diverse therapeutic applications; however, adverse antidrug immune responses and associated loss of efficacy have been reported in clinical uses. Here, we describe computational design and ultrahigh-throughput screening of a massive BoNT/A light-chain (BoNT/A-LC) library optimized for reduced T cell epitope content and thereby dampened immunogenicity. We developed a functional assay based on bacterial co-expression of BoNT/A-LC library members with a Förster resonance energy transfer (FRET) sensor for BoNT/A-LC enzymatic activity, and we employed high-speed fluorescence-activated cell sorting (FACS) to identify numerous computationally designed variants having wild-type-like enzyme kinetics. Many of these variants exhibited decreased immunogenicity in humanized HLA transgenic mice and manifested in vivo paralytic activity when incorporated into full-length toxin. One variant achieved near-wild-type paralytic potency and a 300% reduction in antidrug antibody response in vivo. Thus, we have achieved a striking level of BoNT/A-LC functional deimmunization by combining computational library design and ultrahigh-throughput screening. This strategy holds promise for deimmunizing other biologics with complex superstructures and mechanisms of action.

Deimmunization of protein therapeutics - Recent advances in experimental and computational epitope prediction and deletion

Biotherapeutics, and antimicrobial proteins in particular, are of increasing interest for human medicine. An important challenge in the development of such therapeutics is their potential immunogenicity, which can induce production of anti-drug-antibodies, resulting in altered pharmacokinetics, reduced efficacy, and potentially severe anaphylactic or hypersensitivity reactions. For this reason, the development and application of effective deimmunization methods for protein drugs is of utmost importance. Deimmunization may be achieved by unspecific shielding approaches, which include PEGylation, fusion to polypeptides (e.g., XTEN or PAS), reductive methylation, glycosylation, and polysialylation. Alternatively, the identification of epitopes for T cells or B cells and their subsequent deletion through site-directed mutagenesis represent promising deimmunization strategies and can be accomplished through either experimental or computational approaches. This review highlights the most recent advances and current challenges in the deimmunization of protein therapeutics, with a special focus on computational epitope prediction and deletion tools.

Globally deimmunized lysostaphin evades human immune surveillance and enables highly efficacious repeat dosing

There is a critical need for novel therapies to treat methicillin-resistant Staphylococcus aureus (MRSA) and other drug-resistant pathogens, and lysins are among the vanguard of innovative antibiotics under development. Unfortunately, lysins' own microbial origins can elicit detrimental antidrug antibodies (ADAs) that undermine efficacy and threaten patient safety. To create an enhanced anti-MRSA lysin, a novel variant of lysostaphin was engineered by T cell epitope deletion. This "deimmunized" lysostaphin dampened human T cell activation, mitigated ADA responses in human HLA transgenic mice, and enabled safe and efficacious repeated dosing during a 6-week longitudinal infection study. Furthermore, the deimmunized lysostaphin evaded established anti-wild-type immunity, thereby providing significant anti-MRSA protection for animals that were immune experienced to the wild-type enzyme. Last, the enzyme synergized with daptomycin to clear a stringent model of MRSA endocarditis. By mitigating T cell-driven antidrug immunity, deimmunized lysostaphin may enable safe, repeated dosing to treat refractory MRSA infections.

CodonAdjust: a software for in silico design of a mutagenesis library with specific amino acid profiles

In protein engineering, generation of mutagenesis libraries is a key step to study the functions of mutants. To generate mutants with a desired composition of amino acids (AAs), a codon consisting of a mixture of nucleotides is widely applied. Several computational methods have been proposed to calculate a codon nucleotide composition for generating a given amino acid profile based on mathematical optimization. However, these previous methods need to manually tune weights of amino acids in objective functions, which are time-consuming and, more importantly, lack publicly available software implementations. Here, we develop CodonAdjust, a software to adjust a codon nucleotide composition for mimicking a given amino acid profile. We propose different options of CodonAdjust, which provide various customizations in practical scenarios such as setting a guaranteeing threshold for the frequencies of amino acids without any manual tasks. We demonstrate the capability of CodonAdjust in the experiments on the complementarity-determining regions (CDRs) of antibodies and T-cell receptors (TCRs) as well as millions of amino acid profiles from Pfam. These results suggest that CodonAdjust is a productive software for codon design and may accelerate library generation. CodonAdjust is freely available at https://github.com/tiffany-nguyen/CodonAdjust. Paper edited by Dr. Jeffery Saven, Board Member for PEDS.

Pareto Optimization of Combinatorial Mutagenesis Libraries

In order to increase the hit rate of discovering diverse, beneficial protein variants via high-throughput screening, we have developed a computational method to optimize combinatorial mutagenesis libraries for overall enrichment in two distinct properties of interest. Given scoring functions for evaluating individual variants, POCoM (Pareto Optimal Combinatorial Mutagenesis) scores entire libraries in terms of averages over their constituent members, and designs optimal libraries as sets of mutations whose combinations make the best trade-offs between average scores. This represents the first general-purpose method to directly design combinatorial libraries for multiple objectives characterizing their constituent members. Despite being rigorous in mapping out the Pareto frontier, it is also very fast even for very large libraries (e.g., designing 30 mutation, billion-member libraries in only hours). We here instantiate POCoM with scores based on a target's protein structure and its homologs' sequences, enabling the design of libraries containing variants balancing these two important yet quite different types of information. We demonstrate POCoM's generality and power in case study applications to green fluorescent protein, cytochrome P450, and β-lactamase. Analysis of the POCoM library designs provides insights into the trade-offs between structure- and sequence-based scores, as well as the impacts of experimental constraints on library designs. POCoM libraries incorporate mutations that have previously been found favorable experimentally, while diversifying the contexts in which these mutations are situated and maintaining overall variant quality.

A machine learning approach for reliable prediction of amino acid interactions and its application in the directed evolution of enantioselective enzymes

Directed evolution is an important research activity in synthetic biology and biotechnology. Numerous reports describe the application of tedious mutation/screening cycles for the improvement of proteins. Recently, knowledge-based approaches have facilitated the prediction of protein properties and the identification of improved mutants. However, epistatic phenomena constitute an obstacle which can impair the predictions in protein engineering. We present an innovative sequence-activity relationship (innov'SAR) methodology based on digital signal processing combining wet-lab experimentation and computational protein design. In our machine learning approach, a predictive model is developed to find the resulting property of the protein when the n single point mutations are permuted (2n combinations). The originality of our approach is that only sequence information and the fitness of mutants measured in the wet-lab are needed to build models. We illustrate the application of the approach in the case of improving the enantioselectivity of an epoxide hydrolase from Aspergillus niger. n = 9 single point mutants of the enzyme were experimentally assessed for their enantioselectivity and used as a learning dataset to build a model. Based on combinations of the 9 single point mutations (29), the enantioselectivity of these 512 variants were predicted, and candidates were experimentally checked: better mutants with higher enantioselectivity were indeed found.

Application of fourier transform and proteochemometrics principles to protein engineering

Background

Connecting the dots between the protein sequence and its function is of fundamental interest for protein engineers. In-silico methods are useful in this quest especially when structural information is not available. In this study we propose a mutant library screening tool called iSAR (innovative Sequence Activity Relationship) that relies on the physicochemical properties of the amino acids, digital signal processing and partial least squares regression to uncover these sequence-function correlations.

Results

We show that the digitalized representation of the protein sequence in the form of a Fourier spectrum can be used as an efficient descriptor to model the sequence-activity relationship of proteins. The iSAR methodology that we have developed identifies high fitness mutants from mutant libraries relying on physicochemical properties of the amino acids, digital signal processing and regression techniques. iSAR correlates variations caused by mutations in spectra with biological activity/fitness. It takes into account the impact of mutations on the whole spectrum and does not focus on local fitness alone. The utility of the method is illustrated on 4 datasets: cytochrome P450 for thermostability, TNF-alpha for binding affinity, GLP-2 for potency and enterotoxins for thermostability. The choice of the datasets has been made such as to illustrate the ability of the method to perform when limited training data is available and also when novel mutations appear in the test set, that have not been featured in the training set.

Conclusion

The combination of Fast Fourier Transform and Partial Least Squares regression is efficient in capturing the effects of mutations on the function of the protein. iSAR is a fast algorithm which can be implemented with limited computational resources and can make effective predictions even if the training set is limited in size.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2407-8) contains supplementary material, which is available to authorized users.

Enriching Peptide Libraries for Binding Affinity and Specificity Through Computationally Directed Library Design

Peptide reagents with high affinity or specificity for their target protein interaction partner are of utility for many important applications. Optimization of peptide binding by screening large libraries is a proven and powerful approach. Libraries designed to be enriched in peptide sequences that are predicted to have desired affinity or specificity characteristics are more likely to yield success than random mutagenesis. We present a library optimization method in which the choice of amino acids to encode at each peptide position can be guided by available experimental data or structure-based predictions. We discuss how to use analysis of predicted library performance to inform rounds of library design. Finally, we include protocols for more complex library design procedures that consider the chemical diversity of the amino acids at each peptide position and optimize a library score based on a user-specified input model.

Computationally optimized deimmunization libraries yield highly mutated enzymes with low immunogenicity and enhanced activity

Therapeutic proteins of wide-ranging function hold great promise for treating disease, but immune surveillance of these macromolecules can drive an antidrug immune response that compromises efficacy and even undermines safety. To eliminate widespread T-cell epitopes in any biotherapeutic and thereby mitigate this key source of detrimental immune recognition, we developed a Pareto optimal deimmunization library design algorithm that optimizes protein libraries to account for the simultaneous effects of combinations of mutations on both molecular function and epitope content. Active variants identified by high-throughput screening are thus inherently likely to be deimmunized. Functional screening of an optimized 10-site library (1,536 variants) of P99 β-lactamase (P99βL), a component of ADEPT cancer therapies, revealed that the population possessed high overall fitness, and comprehensive analysis of peptide-MHC II immunoreactivity showed the population possessed lower average immunogenic potential than the wild-type enzyme. Although similar functional screening of an optimized 30-site library (2.15 × 10⁹ variants) revealed reduced population-wide fitness, numerous individual variants were found to have activity and stability better than the wild type despite bearing 13 or more deimmunizing mutations per enzyme. The immunogenic potential of one highly active and stable 14-mutation variant was assessed further using ex vivo cellular immunoassays, and the variant was found to silence T-cell activation in seven of the eight blood donors who responded strongly to wild-type P99βL. In summary, our multiobjective library-design process readily identified large and mutually compatible sets of epitope-deleting mutations and produced highly active but aggressively deimmunized constructs in only one round of library screening.

EpiSweep: Computationally Driven Reengineering of Therapeutic Proteins to Reduce Immunogenicity While Maintaining Function

Therapeutic proteins are yielding ever more advanced and efficacious new drugs, but the biological origins of these highly effective therapeutics render them subject to immune surveillance within the patient's body. When recognized by the immune system as a foreign agent, protein drugs elicit a coordinated response that can manifest a range of clinical complications including rapid drug clearance, loss of functionality and efficacy, delayed infusion-like allergic reactions, more serious anaphylactic shock, and even induced auto-immunity. It is thus often necessary to deimmunize an exogenous protein in order to enable its clinical application; critically, the deimmunization process must also maintain the desired therapeutic activity.To meet the growing need for effective, efficient, and broadly applicable protein deimmunization technologies, we have developed the EpiSweep suite of protein design algorithms. EpiSweep seamlessly integrates computational prediction of immunogenic T cell epitopes with sequence- or structure-based assessment of the impacts of mutations on protein stability and function, in order to select combinations of mutations that make Pareto optimal trade-offs between the competing goals of low immunogenicity and high-level function. The methods are applicable both to the design of individual functionally deimmunized variants as well as the design of combinatorial libraries enriched in functionally deimmunized variants. After validating EpiSweep in a series of retrospective case studies providing comparisons to conventional approaches to T cell epitope deletion, we have experimentally demonstrated it to be highly effective in prospective application to deimmunization of a number of different therapeutic candidates. We conclude that our broadly applicable computational protein design algorithms guide the engineer towards the most promising deimmunized therapeutic candidates, and thereby have the potential to accelerate development of new protein drugs by shortening time frames and improving hit rates.

Depletion of T cell epitopes in lysostaphin mitigates anti-drug antibody response and enhances antibacterial efficacy in vivo

The enzyme lysostaphin possesses potent anti-staphylococcal activity and represents a promising antibacterial drug candidate; however, its immunogenicity poses a barrier to clinical translation. Here, structure-based biomolecular design enabled widespread depletion of lysostaphin DRB1(∗)0401 restricted T cell epitopes, and resulting deimmunized variants exhibited striking reductions in anti-drug antibody responses upon administration to humanized HLA-transgenic mice. This reduced immunogenicity translated into improved efficacy in the form of protection against repeated challenges with methicillin-resistant Staphylococcus aureus (MRSA). In contrast, while wild-type lysostaphin was efficacious against the initial MRSA infection, it failed to clear subsequent bacterial challenges that were coincident with escalating anti-drug antibody titers. These results extend the existing deimmunization literature, in which reduced immunogenicity and retained efficacy are assessed independently of each other. By correlating in vivo efficacy with longitudinal measures of anti-drug antibody development, we provide the first direct evidence that T cell epitope depletion manifests enhanced biotherapeutic efficacy.