Preclinical proof of principle for orally delivered Th17 antagonist miniproteins.

Interleukin (IL)-23 and IL-17 are well-validated therapeutic targets in autoinflammatory diseases. Antibodies targeting IL-23 and IL-17 have shown clinical efficacy but are limited by high costs, safety risks, lack of sustained efficacy, and poor patient convenience as they require parenteral administration. Here, we present designed miniproteins inhibiting IL-23R and IL-17 with antibody-like, low picomolar affinities at a fraction of the molecular size. The minibinders potently block cell signaling in vitro and are extremely stable, enabling oral administration and low-cost manufacturing. The orally administered IL-23R minibinder shows efficacy better than a clinical anti-IL-23 antibody in mouse colitis and has a favorable pharmacokinetics (PK) and biodistribution profile in rats. This work demonstrates that orally administered de novo-designed minibinders can reach a therapeutic target past the gut epithelial barrier. With high potency, gut stability, and straightforward manufacturability, de novo-designed minibinders are a promising modality for oral biologics.

Elicitation of Potent Neutralizing Antibody Responses by Designed Protein Nanoparticle Vaccines for SARS-CoV-2.

A safe, effective, and scalable vaccine is needed to halt the ongoing SARS-CoV-2 pandemic. We describe the structure-based design of self-assembling protein nanoparticle immunogens that elicit potent and protective antibody responses against SARS-CoV-2 in mice. The nanoparticle vaccines display 60 SARS-CoV-2 spike receptor-binding domains (RBDs) in a highly immunogenic array and induce neutralizing antibody titers 10-fold higher than the prefusion-stabilized spike despite a 5-fold lower dose. Antibodies elicited by the RBD nanoparticles target multiple distinct epitopes, suggesting they may not be easily susceptible to escape mutations, and exhibit a lower binding:neutralizing ratio than convalescent human sera, which may minimize the risk of vaccine-associated enhanced respiratory disease. The high yield and stability of the assembled nanoparticles suggest that manufacture of the nanoparticle vaccines will be highly scalable. These results highlight the utility of robust antigen display platforms and have launched cGMP manufacturing efforts to advance the SARS-CoV-2-RBD nanoparticle vaccine into the clinic.

Induction of Potent Neutralizing Antibody Responses by a Designed Protein Nanoparticle Vaccine for Respiratory Syncytial Virus.

Respiratory syncytial virus (RSV) is a worldwide public health concern for which no vaccine is available. Elucidation of the prefusion structure of the RSV F glycoprotein and its identification as the main target of neutralizing antibodies have provided new opportunities for development of an effective vaccine. Here, we describe the structure-based design of a self-assembling protein nanoparticle presenting a prefusion-stabilized variant of the F glycoprotein trimer (DS-Cav1) in a repetitive array on the nanoparticle exterior. The two-component nature of the nanoparticle scaffold enabled the production of highly ordered, monodisperse immunogens that display DS-Cav1 at controllable density. In mice and nonhuman primates, the full-valency nanoparticle immunogen displaying 20 DS-Cav1 trimers induced neutralizing antibody responses ∼10-fold higher than trimeric DS-Cav1. These results motivate continued development of this promising nanoparticle RSV vaccine candidate and establish computationally designed two-component nanoparticles as a robust and customizable platform for structure-based vaccine design.

Principles of Protein Stability and Their Application in Computational Design.

Proteins are increasingly used in basic and applied biomedical research. Many proteins, however, are only marginally stable and can be expressed in limited amounts, thus hampering research and applications. Research has revealed the thermodynamic, cellular, and evolutionary principles and mechanisms that underlie marginal stability. With this growing understanding, computational stability design methods have advanced over the past two decades starting from methods that selectively addressed only some aspects of marginal stability. Current methods are more general and, by combining phylogenetic analysis with atomistic design, have shown drastic improvements in solubility, thermal stability, and aggregation resistance while maintaining the protein's primary molecular activity. Stability design is opening the way to rational engineering of improved enzymes, therapeutics, and vaccines and to the application of protein design methodology to large proteins and molecular activities that have proven challenging in the past.

Protein design with a machine-learned potential about backbone designability.

Proteins are fundamental molecules that mediate diverse biological processes, and protein design can shed light on the molecular mechanisms underlying their biological functions. Huang and colleagues have developed a sequence-independent statistical model for de novo protein design using neural networks (NNs) to learn the distribution of backbone structures with minimal side-chain information.

In silico evolution of autoinhibitory domains for a PD-L1 antagonist using deep learning models.

There has been considerable progress in the development of computational methods for designing protein-protein interactions, but engineering high-affinity binders without extensive screening and maturation remains challenging. Here, we test a protein design pipeline that uses iterative rounds of deep learning (DL)-based structure prediction (AlphaFold2) and sequence optimization (ProteinMPNN) to design autoinhibitory domains (AiDs) for a PD-L1 antagonist. With the goal of creating an anticancer agent that is inactive until reaching the tumor environment, we sought to create autoinhibited (or masked) forms of the PD-L1 antagonist that can be unmasked by tumor-enriched proteases. Twenty-three de novo designed AiDs, varying in length and topology, were fused to the antagonist with a protease-sensitive linker, and binding to PD-L1 was measured with and without protease treatment. Nine of the fusion proteins demonstrated conditional binding to PD-L1, and the top-performing AiDs were selected for further characterization as single-domain proteins. Without any experimental affinity maturation, four of the AiDs bind to the PD-L1 antagonist with equilibrium dissociation constants (KDs) below 150 nM, with the lowest KD equal to 0.9 nM. Our study demonstrates that DL-based protein modeling can be used to rapidly generate high-affinity protein binders.

Decoding CRISPR-Cas PAM recognition with UniDesign.

The critical first step in Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-associated (CRISPR-Cas) protein-mediated gene editing is recognizing a preferred protospacer adjacent motif (PAM) on target DNAs by the protein's PAM-interacting amino acids (PIAAs). Thus, accurate computational modeling of PAM recognition is useful in assisting CRISPR-Cas engineering to relax or tighten PAM requirements for subsequent applications. Here, we describe a universal computational protein design framework (UniDesign) for designing protein-nucleic acid interactions. As a proof of concept, we applied UniDesign to decode the PAM-PIAA interactions for eight Cas9 and two Cas12a proteins. We show that, given native PIAAs, the UniDesign-predicted PAMs are largely identical to the natural PAMs of all Cas proteins. In turn, given natural PAMs, the computationally redesigned PIAA residues largely recapitulated the native PIAAs (74% and 86% in terms of identity and similarity, respectively). These results demonstrate that UniDesign faithfully captures the mutual preference between natural PAMs and native PIAAs, suggesting it is a useful tool for engineering CRISPR-Cas and other nucleic acid-interacting proteins. UniDesign is open-sourced at https://github.com/tommyhuangthu/UniDesign.

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.

Click, Compute, Create: A Review of Web-based Tools for Enzyme Engineering.

Enzyme engineering, though pivotal across various biotechnological domains, is often plagued by its time-consuming and labor-intensive nature. This review aims to offer an overview of supportive in silico methodologies for this demanding endeavor. Starting from methods to predict protein structures, to classification of their activity and even the discovery of new enzymes we continue with describing tools used to increase thermostability and production yields of selected targets. Subsequently, we discuss computational methods to modulate both, the activity as well as selectivity of enzymes. Last, we present recent approaches based on cutting-edge machine learning methods to redesign enzymes. With exception of the last chapter, there is a strong focus on methods easily accessible via web-interfaces or simple Python-scripts, therefore readily useable for a diverse and broad community.

Recent advances in de novo computational design and redesign of intrinsically disordered proteins and intrinsically disordered protein regions.

In the early 2000s, the concept of "unstructured biology" has emerged to be an important field in protein science by generating various new research directions. Many novel strategies and methods have been developed that are focused on effectively identifying/predicting intrinsically disordered proteins (IDPs) and intrinsically disordered protein regions (IDPRs), identifying their potential functions, disorder based drug design etc. Due to the range of functions of IDPs/IDPRs and their involvement in various debilitating diseases they are of contemporary interest to the scientific community. Recent researches are focused on designing/redesigning specific IDPs/IDPRs de novo. These de novo design/redesigns of IDPs/IDPRs are carried out by altering compositional biases and specific sequence patterning parameters. The main focus of these researches is to influence specific molecular functions, phase behavior, cellular phenotypes etc. In this review, we first provide the differences of natively folded and natively unfolded or IDPs with respect to their potential energy landscapes. Here, we provide current understandings on the different computational design strategies and methods that have been utilized in de novo design and redesigns of IDPs and IDPRs. Finally, we conclude the review by discussing the challenges that have been faced during the computational design/design attempts of IDPs/IDPRs.

Transferrin receptor targeting by de novo sheet extension.

The de novo design of polar protein-protein interactions is challenging because of the thermodynamic cost of stripping water away from the polar groups. Here, we describe a general approach for designing proteins which complement exposed polar backbone groups at the edge of beta sheets with geometrically matched beta strands. We used this approach to computationally design small proteins that bind to an exposed beta sheet on the human transferrin receptor (hTfR), which shuttles interacting proteins across the blood-brain barrier (BBB), opening up avenues for drug delivery into the brain. We describe a design which binds hTfR with a 20 nM Kd, is hyperstable, and crosses an in vitro microfluidic organ-on-a-chip model of the human BBB. Our design approach provides a general strategy for creating binders to protein targets with exposed surface beta edge strands.

De novo design and Rosetta-based assessment of high-affinity antibody variable regions (Fv) against the SARS-CoV-2 spike receptor binding domain (RBD).

The continued emergence of new SARS-CoV-2 variants has accentuated the growing need for fast and reliable methods for the design of potentially neutralizing antibodies (Abs) to counter immune evasion by the virus. Here, we report on the de novo computational design of high-affinity Ab variable regions (Fv) through the recombination of VDJ genes targeting the most solvent-exposed hACE2-binding residues of the SARS-CoV-2 spike receptor binding domain (RBD) protein using the software tool OptMAVEn-2.0. Subsequently, we carried out computational affinity maturation of the designed variable regions through amino acid substitutions for improved binding with the target epitope. Immunogenicity of designs was restricted by preferring designs that match sequences from a 9-mer library of "human Abs" based on a human string content score. We generated 106 different antibody designs and reported in detail on the top five that trade-off the greatest computational binding affinity for the RBD with human string content scores. We further describe computational evaluation of the top five designs produced by OptMAVEn-2.0 using a Rosetta-based approach. We used Rosetta SnugDock for local docking of the designs to evaluate their potential to bind the spike RBD and performed "forward folding" with DeepAb to assess their potential to fold into the designed structures. Ultimately, our results identified one designed Ab variable region, P1.D1, as a particularly promising candidate for experimental testing. This effort puts forth a computational workflow for the de novo design and evaluation of Abs that can quickly be adapted to target spike epitopes of emerging SARS-CoV-2 variants or other antigenic targets.

Computational design of closely related proteins that adopt two well-defined but structurally divergent folds.

The plasticity of naturally occurring protein structures, which can change shape considerably in response to changes in environmental conditions, is critical to biological function. While computational methods have been used for de novo design of proteins that fold to a single state with a deep free-energy minimum [P.-S. Huang, S. E. Boyken, D. Baker, Nature 537, 320-327 (2016)], and to reengineer natural proteins to alter their dynamics [J. A. Davey, A. M. Damry, N. K. Goto, R. A. Chica, Nat. Chem. Biol. 13, 1280-1285 (2017)] or fold [P. A. Alexander, Y. He, Y. Chen, J. Orban, P. N. Bryan, Proc. Natl. Acad. Sci. U.S.A. 106, 21149-21154 (2009)], the de novo design of closely related sequences which adopt well-defined but structurally divergent structures remains an outstanding challenge. We designed closely related sequences (over 94% identity) that can adopt two very different homotrimeric helical bundle conformations-one short (∼66 Šheight) and the other long (∼100 Šheight)-reminiscent of the conformational transition of viral fusion proteins. Crystallographic and NMR spectroscopic characterization shows that both the short- and long-state sequences fold as designed. We sought to design bistable sequences for which both states are accessible, and obtained a single designed protein sequence that populates either the short state or the long state depending on the measurement conditions. The design of sequences which are poised to adopt two very different conformations sets the stage for creating large-scale conformational switches between structurally divergent forms.

The Versatile Biocatalyst of Cytochrome P450 CYP102A1: Structure, Function, and Engineering.

Wild-type cytochrome P450 CYP102A1 from Bacillus megaterium is a highly efficient monooxygenase for the oxidation of long-chain fatty acids. The unique features of CYP102A1, such as high catalytic activity, expression yield, regio- and stereoselectivity, and self-sufficiency in electron transfer as a fusion protein, afford the requirements for an ideal biocatalyst. In the past three decades, remarkable progress has been made in engineering CYP102A1 for applications in drug discovery, biosynthesis, and biotechnology. The repertoire of engineered CYP102A1 variants has grown tremendously, whereas the substrate repertoire is avalanched to encompass alkanes, alkenes, aromatics, organic solvents, pharmaceuticals, drugs, and many more. In this article, we highlight the major advances in the past five years in our understanding of the structure and function of CYP102A1 and the methodologies used to engineer CYP102A1 for novel applications. The objective is to provide a succinct review of the latest developments with reference to the body of CYP102A1-related literature.

Recent advances in de novo protein design: Principles, methods, and applications.

The computational de novo protein design is increasingly applied to address a number of key challenges in biomedicine and biological engineering. Successes in expanding applications are driven by advances in design principles and methods over several decades. Here, we review recent innovations in major aspects of the de novo protein design and include how these advances were informed by principles of protein architecture and interactions derived from the wealth of structures in the Protein Data Bank. We describe developments in de novo generation of designable backbone structures, optimization of sequences, design scoring functions, and the design of the function. The advances not only highlight design goals reachable now but also point to the challenges and opportunities for the future of the field.

Modularity-based parallel protein design algorithm with an implementation using shared memory programming.

Given a target protein structure, the prime objective of protein design is to find amino acid sequences that will fold/acquire to the given three-dimensional structure. The protein design problem belongs to the non-deterministic polynomial-time-hard class as sequence search space increases exponentially with protein length. To ensure better search space exploration and faster convergence, we propose a protein modularity-based parallel protein design algorithm. The modular architecture of the protein structure is exploited by considering an intermediate structural organization between secondary structure and domain defined as protein unit (PU). Here, we have incorporated a divide-and-conquer approach where a protein is split into PUs and each PU region is explored in a parallel fashion. It has been further analyzed that our shared memory implementation of modularity-based parallel sequence search leads to better search space exploration compared to the case of traditional full protein design. Sequence-based analysis on design sequences depicts an average of 39.7% sequence similarity on the benchmark data set. Structure-based comparison of the modeled structures of the design protein with the target structure exhibited an average root-mean-square deviation of 1.17 Å and an average template modeling score of 0.89. The selected modeled structures of the design protein sequences are validated using 100 ns molecular dynamics simulations where 80% of the proteins have shown better or similar stability to the respective target proteins. Our study informs that our modularity-based protein design algorithm can be extended to protein interaction design as well.

Advances in the Computational Design of Small-Molecule-Controlled Protein-Based Circuits for Synthetic Biology.

Synthetic biology approaches living systems with an engineering perspective and promises to deliver solutions to global challenges in healthcare and sustainability. A critical component is the design of biomolecular circuits with programmable input-output behaviors. Such circuits typically rely on a sensor module that recognizes molecular inputs, which is coupled to a functional output via protein-level circuits or regulating the expression of a target gene. While gene expression outputs can be customized relatively easily by exchanging the target genes, sensing new inputs is a major limitation. There is a limited repertoire of sensors found in nature, and there are often difficulties with interfacing them with engineered circuits. Computational protein design could be a key enabling technology to address these challenges, as it allows for the engineering of modular and tunable sensors that can be tailored to the circuit's application. In this article, we review recent computational approaches to design protein-based sensors for small-molecule inputs with particular focus on those based on the widely used Rosetta software suite. Furthermore, we review mechanisms that have been harnessed to couple ligand inputs to functional outputs. Based on recent literature, we illustrate how the combination of protein design and synthetic biology enables new sensors for diverse applications ranging from biomedicine to metabolic engineering. We conclude with a perspective on how strategies to address frontiers in protein design and cellular circuit design may enable the next generation of sense-response networks, which may increasingly be assembled from de novo components to display diverse and engineerable input-output behaviors.

Iterated local search with partition crossover for computational protein design.

Structure-based computational protein design (CPD) refers to the problem of finding a sequence of amino acids which folds into a specific desired protein structure, and possibly fulfills some targeted biochemical properties. Recent studies point out the particularly rugged CPD energy landscape, suggesting that local search optimization methods should be designed and tuned to easily escape local minima attraction basins. In this article, we analyze the performance and search dynamics of an iterated local search (ILS) algorithm enhanced with partition crossover. Our algorithm, PILS, quickly finds local minima and escapes their basins of attraction by solution perturbation. Additionally, the partition crossover operator exploits the structure of the residue interaction graph in order to efficiently mix solutions and find new unexplored basins. Our results on a benchmark of 30 proteins of various topology and size show that PILS consistently finds lower energy solutions compared to Rosetta fixbb and a classic ILS, and that the corresponding sequences are mostly closer to the native.

Novel variants of engineered water soluble mu opioid receptors with extensive mutations and removal of cysteines.

We have shown that water-soluble variants of the human mu opioid receptor (wsMOR) containing a reduced number of hydrophobic residues at the lipid-facing residues of the transmembrane (TM) helices can be expressed in E. coli. In this study, we tested the consequences of increasing the number of mutations on the surface of the transmembrane domain on the receptor's aqueous solubility and ligand binding properties, along with mutation of 11 cysteine residues regardless of their solvent exposure value and location in the protein. We computationally engineered 10 different variants of MOR, and tested four of them for expression in E. coli. We found that all four variants were successfully expressed and could be purified in high quantities. The variants have alpha helical structural content similar to that of the native MOR, and they also display binding affinities for the MOR antagonist (naltrexone) similar to the wsMOR variants we engineered previously that contained many fewer mutations. Furthermore, for these full-length variants, the helical content remains unchanged over a wide range of pH values (pH 6 ~ 9). This study demonstrates the flexibility and robustness of the water-soluble MOR variants with respect to additional designed mutations in the TM domain and changes in pH, whereupon the protein's structural integrity and its ligand binding affinity are maintained. These variants of the full-length MOR with less hydrophobic surface residues and less cysteines can be obtained in large amounts from expression in E. coli and can serve as novel tools to investigate structure-function relationships of the receptor.

Perturbing the energy landscape for improved packing during computational protein design.

The FastDesign protocol in the molecular modeling program Rosetta iterates between sequence optimization and structure refinement to stabilize de novo designed protein structures and complexes. FastDesign has been used previously to design novel protein folds and assemblies with important applications in research and medicine. To promote sampling of alternative conformations and sequences, FastDesign includes stages where the energy landscape is smoothened by reducing repulsive forces. Here, we discover that this process disfavors larger amino acids in the protein core because the protein compresses in the early stages of refinement. By testing alternative ramping strategies for the repulsive weight, we arrive at a scheme that produces lower energy designs with more native-like sequence composition in the protein core. We further validate the protocol by designing and experimentally characterizing over 4000 proteins and show that the new protocol produces higher stability proteins.