Search

High-throughput deep mutational scanning (DMS) experiments have significantly impacted protein engineering, drug discovery, immunology, cancer biology, and evolutionary biology by enabling the systematic understanding of protein functions. However, the mutational space associated with proteins is astronomically large, making it overwhelming for current experimental capabilities. Therefore, alternative methods for DMS are imperative. We propose a topological deep learning (TDL) paradigm to facilitate in silico DMS. We utilize a new topological data analysis (TDA) technique based on the persistent spectral theory, also known as persistent Laplacian, to capture both topological invariants and the homotopic shape evolution of data. To validate our TDL-DMS model, we use SARS-CoV-2 datasets and show excellent accuracy and reliability for binding interface mutations. This finding is significant for SARS-CoV-2 variant forecasting and designing effective antibodies and vaccines. Our proposed model is expected to have a significant impact on drug discovery, vaccine design, precision medicine, and protein engineering., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Yeast surface display selections against mammalian cell monolayers have proven effective in isolating proteins with novel binding activity. Recent advances in this technique allow for the recovery of clones with even micromolar binding affinities. However, no efficient method has been shown for affinity-based selection in this context. This study demonstrates the effectiveness of titratable avidity reduction using dithiothreitol to achieve this goal. A series of epidermal growth factor receptor binding fibronectin domains with a range of affinities are used to quantitatively identify the number of ligands per yeast cell that yield the strongest selectivity between strong, moderate, and weak affinities. Notably, reduction of ligand display to 3,000-6,000 ligands per yeast cell of a 2 nM binder yields 16-fold better selectivity than that to a 17 nM binder. These lessons are applied to affinity maturation of an EpCAM-binding fibronectin population, yielding an enriched pool of ligands with significantly stronger affinity than that of an analogous pool sorted by standard cellular selection methods. Collectively, this study offers a facile approach for affinity selection of yeast-displayed ligands against full-length cellular targets and demonstrates the effectiveness of this method by generating EpCAM-binding ligands that are promising for further applications.

Engineered proteins provide clinically and industrially impactful molecules and utility within fundamental research, yet inefficiencies in discovering lead variants with new desired functionality, while maintaining stability, hinder progress. Improved function, which can result from a few strategic mutations, is fundamentally separate from discovering novel function, which often requires large leaps in sequence space. While a highly diverse combinatorial library covering immense sequence space would empower protein discovery, the ability to sample only a minor subset of sequence space and the typical destabilization of random mutations preclude this strategy. A balance must be reached. At library scale, compounding several destabilizing mutations renders many variants unable to properly fold and devoid of function. Broadly searching sequence space while reducing the level of destabilization may enhance evolution. We exemplify this balance with affibody, a three-helix bundle protein scaffold. Using natural ligand data sets, stability and structural computations, and deep sequencing of thousands of binding variants, a protein library was designed on a sitewise basis with a gradient of mutational levels across 29% of the protein. In direct competition of biased and uniform libraries, both with 1 × 10 9 variants, for discovery of 6 × 10 4 ligands (5 × 10 3 clusters) toward seven targets, biased amino acid frequency increased ligand discovery 13 ± 3-fold. Evolutionarily favorable amino acids, both globally and site-specifically, are further elucidated. The sitewise amino acid bias aids evolutionary discovery by reducing the level of mutant destabilization as evidenced by a midpoint of denaturation (62 ± 4 °C) 15 °C higher than that of unbiased mutants (47 ± 11 °C; p < 0.001). Sitewise diversification, identified by high-throughput evolution and rational library design, improves discovery efficiency.

This purpose of this study is to determine the efficacy of a 45-amino acid Gp2 domain, engineered to bind to epidermal growth factor receptor (EGFR), as a positron emission tomography (PET) probe of EGFR in a xenograft mouse model. The EGFR-targeted Gp2 (Gp2-EGFR) and a nonbinding control were site-specifically labeled with 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid (DOTA) chelator. Binding affinity was tested toward human EGFR and mouse EGFR. Biological activity on downstream EGFR signaling was examined in cell culture. DOTA-Gp2 molecules were labeled with 64 Cu and intravenously injected (0.6-2.3 MBq) into mice bearing EGFR high (n = 7) and EGFR low (n = 4) xenografted tumors. PET/computed tomography (CT) images were acquired at 45 min, 2 h, and 24 h. Dynamic PET (25 min) was also acquired. Tomography results were verified with gamma counting of resected tissues. Two-tailed t tests with unequal variances provided statistical comparison. DOTA-Gp2-EGFR bound strongly to human (K D = 7 ± 5 nM) and murine (K D = 29 ± 6 nM) EGFR, and nontargeted Gp2 had no detectable binding. Gp2-EGFR did not agonize EGFR nor antagonize EGF-EGFR. 64 Cu-Gp2-EGFR tracer effectively localized to EGFR high tumors at 45 min (3.2 ± 0.5%ID/g). High specificity was observed with significantly lower uptake in EGFR low tumors (0.9 ± 0.3%ID/g, p < 0.001), high tumor-to-background ratios (11 ± 6 tumor/muscle, p < 0.001). Nontargeted Gp2 tracer had low uptake in EGFR high tumors (0.5 ± 0.3%ID/g, p < 0.001). Similar data was observed at 2 h, and tumor signal was retained at 24 h (2.9 ± 0.3%ID/g). An engineered Gp2 PET imaging probe exhibited low background and target-specific EGFR high tumor uptake at 45 min, with tumor signal retained at 24 h postinjection, and compared favorably with published EGFR PET probes for alternative protein scaffolds. These beneficial in vivo characteristics, combined with thermal stability, efficient evolution, and small size of the Gp2 domain validate its use as a future class of molecular imaging agents.

ScaffoldSeq is software designed for the numerous applications-including directed evolution analysis-in which a user generates a population of DNA sequences encoding for partially diverse proteins with related functions and would like to characterize the single site and pairwise amino acid frequencies across the population. A common scenario for enzyme maturation, antibody screening, and alternative scaffold engineering involves naïve and evolved populations that contain diversified regions, varying in both sequence and length, within a conserved framework. Analyzing the diversified regions of such populations is facilitated by high-throughput sequencing platforms; however, length variability within these regions (e.g., antibody CDRs) encumbers the alignment process. To overcome this challenge, the ScaffoldSeq algorithm takes advantage of conserved framework sequences to quickly identify diverse regions. Beyond this, unintended biases in sequence frequency are generated throughout the experimental workflow required to evolve and isolate clones of interest prior to DNA sequencing. ScaffoldSeq software uniquely handles this issue by providing tools to quantify and remove background sequences, cluster similar protein families, and dampen the impact of dominant clones. The software produces graphical and tabular summaries for each region of interest, allowing users to evaluate diversity in a site-specific manner as well as identify epistatic pairwise interactions. The code and detailed information are freely available at http://research.cems.umn.edu/hackel. Proteins 2016; 84:869-874. © 2016 Wiley Periodicals, Inc., (© 2016 Wiley Periodicals, Inc.)

Discovering new binding function via a combinatorial library in small protein scaffolds requires balance between appropriate mutations to introduce favorable intermolecular interactions while maintaining intramolecular integrity. Sitewise constraints exist in a non-spatial gradient from diverse to conserved in evolved antibody repertoires; yet non-antibody scaffolds generally do not implement this strategy in combinatorial libraries. Despite the fact that biased amino acid distributions, typically elevated in tyrosine, serine, and glycine, have gained wider use in synthetic scaffolds, these distributions are still predominantly applied uniformly to diversified sites. While select sites in fibronectin domains and DARPins have shown benefit from sitewise designs, they have not been deeply evaluated. Inspired by this disparity between diversity distributions in natural libraries and synthetic scaffold libraries, we hypothesized that binders resulting from discovery and evolution would exhibit a non-spatial, sitewise gradient of amino acid diversity. To identify sitewise diversities consistent with efficient evolution in the context of a hydrophilic fibronectin domain, >105 binders to six targets were evolved and sequenced. Evolutionarily favorable amino acid distributions at 25 sites reveal Shannon entropies (range: 0.3-3.9; median: 2.1; standard deviation: 1.1) supporting the diversity gradient hypothesis. Sitewise constraints in evolved sequences are consistent with complementarity, stability, and consensus biases. Implementation of sitewise constrained diversity enables direct selection of nanomolar affinity binders validating an efficient strategy to balance inter- and intra-molecular interaction demands at each site.

Small protein ligands can provide superior physiological distribution compared with antibodies, and improved stability, production, and specific conjugation. Systematic evaluation of the PDB identified a scaffold to push the limits of small size and robust evolution of stable, high-affinity ligands: 45-residue T7 phage gene 2 protein (Gp2) contains an α helix opposite a β sheet with two adjacent loops amenable to mutation. De novo ligand discovery from 10(8) mutants and directed evolution toward four targets yielded target-specific binders with affinities as strong as 200 ± 100 pM, Tms from 65 °C ± 3 °C to 80°C ± 1 °C, and retained activity after thermal denaturation. For cancer targeting, a Gp2 domain for epidermal growth factor receptor was evolved with 18 ± 8 nM affinity, receptor-specific binding, and high thermal stability with refolding. The efficiency of evolving new binding function and the size, affinity, specificity, and stability of evolved domains render Gp2 a uniquely effective ligand scaffold., (Copyright © 2015 Elsevier Ltd. All rights reserved.)