Computational Protein Design Research Articles

Reaching New Heights in Genetic Code Manipulation with High Throughput Screening.

The chemical and physical properties of proteins are limited by the 20 canonical amino acids. Genetic code manipulation allows for the incorporation of noncanonical amino acids (ncAAs) that enhance or alter protein functionality. This review explores advances in the three main strategies for introducing ncAAs into biosynthesized proteins, focusing on the role of high throughput screening in these advancements. The first section discusses engineering aminoacyl-tRNA synthetases (aaRSs) and tRNAs, emphasizing how novel selection methods improve characteristics including ncAA incorporation efficiency and selectivity. The second section examines high-throughput techniques for improving protein translation machinery, enabling accommodation of alternative genetic codes. This includes opportunities to enhance ncAA incorporation through engineering cellular components unrelated to translation. The final section highlights various discovery platforms for high-throughput screening of ncAA-containing proteins, showcasing innovative binding ligands and enzymes that are challenging to create with only canonical amino acids. These advances have led to promising drug leads and biocatalysts. Overall, the ability to discover unexpected functionalities through high-throughput methods significantly influences ncAA incorporation and its applications. Future innovations in experimental techniques, along with advancements in computational protein design and machine learning, are poised to further elevate this field.

Challenges and Solutions for Leave-One-Out Biosensor Design in the Context of a Rugged Fitness Landscape

The leave-one-out (LOO) green fluorescent protein (GFP) approach to biosensor design combines computational protein design with split protein reconstitution. LOO-GFPs reversibly fold and gain fluorescence upon encountering the target peptide, which can be redefined by computational design of the LOO site. Such an approach can be used to create reusable biosensors for the early detection of emerging biological threats. Enlightening biophysical inferences for nine LOO-GFP biosensor libraries are presented, with target sequences from dengue, influenza, or HIV, replacing beta strands 7, 8, or 11. An initially low hit rate was traced to components of the energy function, manifesting in the over-rewarding of over-tight side chain packing. Also, screening by colony picking required a low library complexity, but designing a biosensor against a peptide of at least 12 residues requires a high-complexity library. This double-bind was solved using a “piecemeal” iterative design strategy. Also, designed LOO-GFPs fluoresced in the unbound state due to unwanted dimerization, but this was solved by fusing a fully functional prototype LOO-GFP to a fiber-forming protein, Drosophila ultrabithorax, creating a biosensor fiber. One influenza hemagglutinin biosensor is characterized here in detail, showing a shifted excitation/emission spectrum, a micromolar affinity for the target peptide, and an unexpected photo-switching ability.

Design of intrinsically disordered protein variants with diverse structural properties.

Intrinsically disordered proteins (IDPs) perform a broad range of functions in biology, suggesting that the ability to design IDPs could help expand the repertoire of proteins with novel functions. Computational design of IDPs with specific conformational properties has, however, been difficult because of their substantial dynamics and structural complexity. We describe a general algorithm for designing IDPs with specific structural properties. We demonstrate the power of the algorithm by generating variants of naturally occurring IDPs that differ in compaction, long-range contacts, and propensity to phase separate. We experimentally tested and validated our designs and analyzed the sequence features that determine conformations. We show how our results are captured by a machine learning model, enabling us to speed up the algorithm. Our work expands the toolbox for computational protein design and will facilitate the design of proteins whose functions exploit the many properties afforded by protein disorder.

Extending computational protein design to intrinsically disordered proteins.

Advances in the accuracy and throughput of molecular simulations usher in a new era in the structural biology of disordered proteins.

Harnessing Deep Learning Methods for Voltage-Gated Ion Channel Drug Discovery.

Voltage-gated ion channels (VGICs) are pivotal in regulating electrical activity in excitable cells and are critical pharmaceutical targets for treating many diseases including cardiac arrhythmia and neuropathic pain. Despite their significance, challenges such as achieving target selectivity persist in VGIC drug development. Recent progress in deep learning, particularly diffusion models, has enabled the computational design of protein binders for any clinically relevant protein based solely on its structure. These developments coincide with a surge in experimental structural data for VGICs, providing a rich foundation for computational design efforts. This review explores the recent advancements in computational protein design using deep learning and diffusion methods, focusing on their application in designing protein binders to modulate VGIC activity. We discuss the potential use of these methods to computationally design protein binders targeting different regions of VGICs, including the pore domain, voltage-sensing domains, and interface with auxiliary subunits. We provide a comprehensive overview of the different design scenarios, discuss key structural considerations, and address the practical challenges in developing VGIC-targeting protein binders. By exploring these innovative computational methods, we aim to provide a framework for developing novel strategies that could significantly advance VGIC pharmacology and lead to the discovery of effective and safe therapeutics.

Bebtelovimab-bound SARS-CoV-2 RBD mutants: resistance profiling and validation with escape mutations, clinical results, and viral genome sequences.

The dynamic evolution of SARS-CoV-2 variants necessitates ongoing advancements in therapeutic strategies. Despite the promise of monoclonal antibody (mAb) therapies like bebtelovimab, concerns persist regarding resistance mutations, particularly single-to-multipoint mutations in the receptor-binding domain (RBD). Our study addresses this by employing interface-guided computational protein design to predict potential bebtelovimab-resistance mutations. Through extensive physicochemical analysis, mutational preferences, precision-recall metrics, protein-protein docking, and energetic analyses, combined with all-atom, and coarse-grained molecular dynamics (MD) simulations, we elucidated the structural-dynamics-binding features of the bebtelovimab-RBD complexes. Identification of susceptible RBD residues under positive selection pressure, coupled with validation against bebtelovimab-escape mutations, clinically reported resistance mutations, and viral genomic sequences enhances the translational significance of our findings and contributes to a better understanding of the resistance mechanisms of SARS-CoV-2.

Computationally designed mRNA-launched protein nanoparticle vaccines.

Both protein nanoparticle and mRNA vaccines were clinically de-risked during the COVID-19 pandemic1-6. These vaccine modalities have complementary strengths: antigen display on protein nanoparticles can enhance the magnitude, quality, and durability of antibody responses7-10, while mRNA vaccines can be rapidly manufactured11 and elicit antigen-specific CD4 and CD8 T cells12,13. Here we leverage a computationally designed icosahedral protein nanoparticle that was redesigned for optimal secretion from eukaryotic cells14 to develop an mRNA-launched nanoparticle vaccine for SARS-CoV-2. The nanoparticle, which displays 60 copies of a stabilized variant of the Wuhan-Hu-1 Spike receptor binding domain (RBD)15, formed monodisperse, antigenically intact assemblies upon secretion from transfected cells. An mRNA vaccine encoding the secreted RBD nanoparticle elicited 5- to 28-fold higher levels of neutralizing antibodies than an mRNA vaccine encoding membrane-anchored Spike, induced higher levels of CD8 T cells than the same immunogen when delivered as an adjuvanted protein nanoparticle, and protected mice from vaccine-matched and -mismatched SARS-CoV-2 challenge. Our data establish that delivering protein nanoparticle immunogens via mRNA vaccines can combine the benefits of each modality and, more broadly, highlight the utility of computational protein design in genetic immunization strategies.

Massively parallel free energy calculations for in silico affinity maturation of designed miniproteins.

Computational protein design efforts continue to make remarkable advances, yet the discovery of high-affinity binders typically requires large-scale experimental screening of site-saturated mutant (SSM) libraries. Here, we explore how massively parallel free energy methods can be used for in silico affinity maturation of de novo designed binding proteins. Using an expanded ensemble (EE) approach, we perform exhaustive relative binding free energy calculations for SSM variants of three miniproteins designed to bind influenza A H1 hemagglutinin by Chevalier et al. (2017). We compare our predictions to experimental ΔΔ G values inferred from a Bayesian analysis of the high-throughput sequencing data, and to state-of-the-art predictions made using the Flex ddG Rosetta protocol. A systematic comparison reveals prediction accuracies around 2 kcal/mol, and identifies net charge changes, large numbers of alchemical atoms, and slow side chain conformational dynamics as key contributors to the uncertainty of the EE predictions. Flex ddG predictions are more accurate on average, but highly conservative. In contrast, EE predictions can better classify stabilizing and destabilizing mutations. We also explored the ability of SSM scans to rationalize known affinity-matured variants containing multiple mutations, which are non-additive due to epistatic effects. Simple electrostatic models fail to explain non-additivity, but observed mutations are found at positions with higher Shannon entropies. Overall, this work suggests that simulation-based free energy methods can provide predictive information for in silico affinity maturation of designed miniproteins, with many feasible improvements to the efficiency and accuracy within reach.

An integrative approach to protein sequence design through multiobjective optimization.

With recent methodological advances in the field of computational protein design, in particular those based on deep learning, there is an increasing need for frameworks that allow for coherent, direct integration of different models and objective functions into the generative design process. Here we demonstrate how evolutionary multiobjective optimization techniques can be adapted to provide such an approach. With the established Non-dominated Sorting Genetic Algorithm II (NSGA-II) as the optimization framework, we use AlphaFold2 and ProteinMPNN confidence metrics to define the objective space, and a mutation operator composed of ESM-1v and ProteinMPNN to rank and then redesign the least favorable positions. Using the two-state design problem of the foldswitching protein RfaH as an in-depth case study, and PapD and calmodulin as examples of higher-dimensional design problems, we show that the evolutionary multiobjective optimization approach leads to significant reduction in the bias and variance in RfaH native sequence recovery, compared to a direct application of ProteinMPNN. We suggest that this improvement is due to three factors: (i) the use of an informative mutation operator that accelerates the sequence space exploration, (ii) the parallel, iterative design process inherent to the genetic algorithm that improves upon the ProteinMPNN autoregressive sequence decoding scheme, and (iii) the explicit approximation of the Pareto front that leads to optimal design candidates representing diverse tradeoff conditions. We anticipate this approach to be readily adaptable to different models and broadly relevant for protein design tasks with complex specifications.

How to Sample Dozens of Substitutions per Site with λ Dynamics.

Alchemical free energy methods are useful in computer-aided drug design and computational protein design because they provide rigorous statistical mechanics-based estimates of free energy differences from molecular dynamics simulations. λ dynamics is a free energy method with the ability to characterize combinatorial chemical spaces spanning thousands of related systems within a single simulation, which gives it a distinct advantage over other alchemical free energy methods that are mostly limited to pairwise comparisons. Recently developed methods have improved the scalability of λ dynamics to perturbations at many sites; however, the size of chemical space that can be explored at each individual site has previously been limited to fewer than ten substituents. As the number of substituents increases, the volume of alchemical space corresponding to nonphysical alchemical intermediates grows exponentially relative to the size corresponding to the physical states of interest. Beyond nine substituents, λ dynamics simulations become lost in an alchemical morass of intermediate states. In this work, we introduce new biasing potentials that circumvent excessive sampling of intermediate states by favoring sampling of physical end points relative to alchemical intermediates. Additionally, we present a more scalable adaptive landscape flattening algorithm for these larger alchemical spaces. Finally, we show that this potential enables more efficient sampling in both protein and drug design test systems with up to 24 substituents per site, enabling, for the first time, simultaneous simulation of all 20 amino acids.

Computational design of de novo bioenergetic membrane proteins.

The major energy-producing reactions of biochemistry occur at biological membranes. Computational protein design now provides the opportunity to elucidate the underlying principles of these processes and to construct bioenergetic pathways on our own terms. Here, we review recent achievements in this endeavour of 'synthetic bioenergetics', with a particular focus on new enabling tools that facilitate the computational design of biocompatible de novo integral membrane proteins. We use recent examples to showcase some of the key computational approaches in current use and highlight that the overall philosophy of 'surface-swapping' - the replacement of solvent-facing residues with amino acids bearing lipid-soluble hydrophobic sidechains - is a promising avenue in membrane protein design. We conclude by highlighting outstanding design challenges and the emerging role of AI in sequence design and structure ideation.

Data hazards in synthetic biology.

Data science is playing an increasingly important role in the design and analysis of engineered biology. This has been fueled by the development of high-throughput methods like massively parallel reporter assays, data-rich microscopy techniques, computational protein structure prediction and design, and the development of whole-cell models able to generate huge volumes of data. Although the ability to apply data-centric analyses in these contexts is appealing and increasingly simple to do, it comes with potential risks. For example, how might biases in the underlying data affect the validity of a result and what might the environmental impact of large-scale data analyses be? Here, we present a community-developed framework for assessing data hazards to help address these concerns and demonstrate its application to two synthetic biology case studies. We show the diversity of considerations that arise in common types of bioengineering projects and provide some guidelines and mitigating steps. Understanding potential issues and dangers when working with data and proactively addressing them will be essential for ensuring the appropriate use of emerging data-intensive AI methods and help increase the trustworthiness of their applications in synthetic biology.

Rationally seeded computational protein design of ɑ-helical barrels

Computational protein design is advancing rapidly. Here we describe efficient routes starting from validated parallel and antiparallel peptide assemblies to design two families of α-helical barrel proteins with central channels that bind small molecules. Computational designs are seeded by the sequences and structures of defined de novo oligomeric barrel-forming peptides, and adjacent helices are connected by loop building. For targets with antiparallel helices, short loops are sufficient. However, targets with parallel helices require longer connectors; namely, an outer layer of helix–turn–helix–turn–helix motifs that are packed onto the barrels. Throughout these computational pipelines, residues that define open states of the barrels are maintained. This minimizes sequence sampling, accelerating the design process. For each of six targets, just two to six synthetic genes are made for expression in Escherichia coli. On average, 70% of these genes express to give soluble monomeric proteins that are fully characterized, including high-resolution structures for most targets that match the design models with high accuracy.

Computational Protein Design - Where it goes?

Proteins have been playing a critical role in the regulation of diverse biological processes related to human life. With the increasing demand, functional proteins are sparse in this immense sequence space. Therefore, protein design has become an important task in various fields, including medicine, food, energy, materials, etc. Directed evolution has recently led to significant achievements. Molecular modification of proteins through directed evolution technology has significantly advanced the fields of enzyme engineering, metabolic engineering, medicine, and beyond. However, it is impossible to identify desirable sequences from a large number of synthetic sequences alone. As a result, computational methods, including data-driven machine learning and physics-based molecular modeling, have been introduced to protein engineering to produce more functional proteins. This review focuses on recent advances in computational protein design, highlighting the applicability of different approaches as well as their limitations.

De novo design of pH-responsive self-assembling helical protein filaments

Biological evolution has led to precise and dynamic nanostructures that reconfigure in response to pH and other environmental conditions. However, designing micrometre-scale protein nanostructures that are environmentally responsive remains a challenge. Here we describe the de novo design of pH-responsive protein filaments built from subunits containing six or nine buried histidine residues that assemble into micrometre-scale, well-ordered fibres at neutral pH. The cryogenic electron microscopy structure of an optimized design is nearly identical to the computational design model for both the subunit internal geometry and the subunit packing into the fibre. Electron, fluorescent and atomic force microscopy characterization reveal a sharp and reversible transition from assembled to disassembled fibres over 0.3 pH units, and rapid fibre disassembly in less than 1 s following a drop in pH. The midpoint of the transition can be tuned by modulating buried histidine-containing hydrogen bond networks. Computational protein design thus provides a route to creating unbound nanomaterials that rapidly respond to small pH changes.

Machine learning-aided design and screening of an emergent protein function in synthetic cells

Recently, utilization of Machine Learning (ML) has led to astonishing progress in computational protein design, bringing into reach the targeted engineering of proteins for industrial and biomedical applications. However, the design of proteins for emergent functions of core relevance to cells, such as the ability to spatiotemporally self-organize and thereby structure the cellular space, is still extremely challenging. While on the generative side conditional generative models and multi-state design are on the rise, for emergent functions there is a lack of tailored screening methods as typically needed in a protein design project, both computational and experimental. Here we describe a proof-of-principle of how such screening, in silico and in vitro, can be achieved for ML-generated variants of a protein that forms intracellular spatiotemporal patterns. For computational screening we use a structure-based divide-and-conquer approach to find the most promising candidates, while for the subsequent in vitro screening we use synthetic cell-mimics as established by Bottom-Up Synthetic Biology. We then show that the best screened candidate can indeed completely substitute the wildtype gene in Escherichia coli. These results raise great hopes for the next level of synthetic biology, where ML-designed synthetic proteins will be used to engineer cellular functions.

Machine learning for functional protein design.

Recent breakthroughs in AI coupled with the rapid accumulation of protein sequence and structure data have radically transformed computational protein design. New methods promise to escape the constraints of natural and laboratory evolution, accelerating the generation of proteins for applications in biotechnology and medicine. To make sense of the exploding diversity of machine learning approaches, we introduce a unifying framework that classifies models on the basis of their use of three core data modalities: sequences, structures and functional labels. We discuss the new capabilities and outstanding challenges for the practical design of enzymes, antibodies, vaccines, nanomachines and more. We then highlight trends shaping the future of this field, from large-scale assays to more robust benchmarks, multimodal foundation models, enhanced sampling strategies and laboratory automation.

DexDesign: an OSPREY-based algorithm for designing de novo D-peptide inhibitors.

With over 270 unique occurrences in the human genome, peptide-recognizing PDZ domains play a central role in modulating polarization, signaling, and trafficking pathways. Mutations in PDZ domains lead to diseases such as cancer and cystic fibrosis, making PDZ domains attractive targets for therapeutic intervention. D-peptide inhibitors offer unique advantages as therapeutics, including increased metabolic stability and low immunogenicity. Here, we introduce DexDesign, a novel OSPREY-based algorithm for computationally designing de novo D-peptide inhibitors. DexDesign leverages three novel techniques that are broadly applicable to computational protein design: the Minimum Flexible Set, K*-based Mutational Scan, and Inverse Alanine Scan. We apply these techniques and DexDesign to generate novel D-peptide inhibitors of two biomedically important PDZ domain targets: CAL and MAST2. We introduce a framework for analyzing de novo peptides-evaluation along a replication/restitution axis-and apply it to the DexDesign-generated D-peptides. Notably, the peptides we generated are predicted to bind their targets tighter than their targets' endogenous ligands, validating the peptides' potential as lead inhibitors. We also provide an implementation of DexDesign in the free and open source computational protein design software OSPREY.

Computationally designed sensors detect endogenous Ras activity and signaling effectors at subcellular resolution.

The utility of genetically encoded biosensors for sensing the activity of signaling proteins has been hampered by a lack of strategies for matching sensor sensitivity to the physiological concentration range of the target. Here we used computational protein design to generate intracellular sensors of Ras activity (LOCKR-based Sensor for Ras activity (Ras-LOCKR-S)) and proximity labelers of the Ras signaling environment (LOCKR-based, Ras activity-dependent Proximity Labeler (Ras-LOCKR-PL)). These tools allow the detection of endogenous Ras activity and labeling of the surrounding environment at subcellular resolution. Using these sensors in human cancer cell lines, we identified Ras-interacting proteins in oncogenic EML4-Alk granules and found that Src-Associated in Mitosis 68-kDa (SAM68) protein specifically enhances Ras activity in the granules. The ability to subcellularly localize endogenous Ras activity should deepen our understanding of Ras function in health and disease and may suggest potential therapeutic strategies.

Interface design of SARS-CoV-2 symmetrical nsp7 dimer and machine learning-guided nsp7 sequence prediction reveals physicochemical properties and hotspots for nsp7 stability, adaptation, and therapeutic design.

The COVID-19 pandemic, driven by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), necessitates a profound understanding of the virus and its lifecycle. As an RNA virus with high mutation rates, SARS-CoV-2 exhibits genetic variability leading to the emergence of variants with potential implications. Among its key proteins, the RNA-dependent RNA polymerase (RdRp) is pivotal for viral replication. Notably, RdRp forms dimers via non-structural protein (nsp) subunits, particularly nsp7, crucial for efficient viral RNA copying. Similar to the main protease (Mpro) of SARS-CoV-2, there is a possibility that the nsp7 might also undergo mutational selection events to generate more stable and adaptable versions of nsp7 dimer during virus evolution. However, efforts to obtain such cohesive and comprehensive information are lacking. To address this, we performed this study focused on deciphering the molecular intricacies of nsp7 dimerization using a multifaceted approach. Leveraging computational protein design (CPD), machine learning (ML), AlphaFold v2.0-based structural analysis, and several related computational approaches, we aimed to identify critical residues and mutations influencing nsp7 dimer stability and adaptation. Our methodology involved identifying potential hotspot residues within the dimeric nsp7 interface using an interface-based CPD approach. Through Rosetta-based symmetrical protein design, we designed and modulated nsp7 dimerization, considering selected interface residues. Analysis of physicochemical features revealed acceptable structural changes and several structural and residue-specific insights emphasizing the intricate nature of such protein-protein complexes. Our ML models, particularly the random forest regressor (RFR), accurately predicted binding affinities and ML-guided sequence predictions corroborated CPD findings, elucidating potential nsp7 mutations and their impact on binding affinity. Validation against clinical sequencing data demonstrated the predictive accuracy of our approach. Moreover, AlphaFold v2.0 structural analyses validated optimal dimeric configurations of affinity-enhancing designs, affirming methodological precision. Affinity-enhancing designs exhibited favourable energetics and higher binding affinity as compared to their counterparts. The obtained physicochemical properties, molecular interactions, and sequence predictions advance our understanding of SARS-CoV-2 evolution and inform potential avenues for therapeutic intervention against COVID-19.