A defined structural unit enables de novo design of small-molecule-binding proteins

Ab initio prediction of specific phospholipid complexes and membrane association of HIV-1 MPER antibodies by multi-scale simulations

A potent class of HIV-1 broadly neutralizing antibodies (bnAbs) targets the envelope glycoprotein's membrane proximal exposed region (MPER) through a proposed mechanism where hypervariable loops embed into lipid bilayers and engage headgroup moieties alongside the epitope. We address the feasibility and determinant molecular features of this mechanism using multi-scale modeling. All-atom simulations of 4E10, PGZL1, 10E8, and LN01 docked onto HIV-like membranes consistently form phospholipid complexes at key complementarity-determining region loop sites, solidifying that stable and specific lipid interactions anchor bnAbs to membrane surfaces. Ancillary protein-lipid contacts reveal surprising contributions from antibody framework regions. Coarse-grained simulations effectively capture antibodies embedding into membranes. Simulations estimating protein-membrane interaction strength for PGZL1 variants along an inferred maturation pathway show bilayer affinity is evolved and correlates with neutralization potency. The modeling demonstrated here uncovers insights into lipid participation in antibodies' recognition of membrane proteins and highlights antibody features to prioritize in vaccine design.

De novo design of transmembrane fluorescence-activating proteins

The recognition of ligands by transmembrane proteins is essential for the exchange of materials, energy and information across biological membranes. Progress has been made in the de novo design of transmembrane proteins1-6, as well as in designing water-soluble proteins to bind small molecules7-12, but de novo design of transmembrane proteins that tightly and specifically bind to small molecules remains an outstanding challenge13. Here we present the accurate design of ligand-binding transmembrane proteins by integrating deep learning and energy-based methods. We designed pre-organized ligand-binding pockets in high-quality four-helix backbones for a fluorogenic ligand, and generated a transmembrane span using gradient-guided hallucination. The designer transmembrane proteins specifically activated fluorescence of the target fluorophore with mid-nanomolar affinity, exhibiting higher brightness and quantum yield compared to those of enhanced green fluorescent protein. These proteins were highly active in the membrane fraction of live bacterial and eukaryotic cells following expression. The crystal and cryogenic electron microscopy structures of the designer protein-ligand complexes were very close to the structures of the design models. We showed that the interactions between ligands and transmembrane proteins within the membrane can be accurately designed. Our work paves the way for the creation of new functional transmembrane proteins, with a wide range of applications including imaging, ligand sensing and membrane transport.

Emergence of specific binding and catalysis from a designed generalist binding protein

The evolution of binding and catalysis played a central role in the emergence of life. While natural proteins have finely tuned affinities for their primary ligands, they also bind weakly and promiscuously to other molecules, which serve as starting points for stepwise, incremental evolution of entirely new specificities. Thus, modern proteins emerged from the joint exploration of sequence and structural space. The ability of natural proteins to bind promiscuously to small molecule fragments has been widely evaluated using methods including crystallographic fragment screening. However, this approach had not been applied to de novo proteins. Here, we apply this method to explore the promiscuity of a de novo small molecule-binding protein ABLE. As in Nature, we found ABLE was capable of forming weak complexes, which were found to be excellent starting points for evolving entirely new functions, including a binder of a turn-on fluorophore and a highly efficient and specific Kemp eliminase enzyme. This work shows how Nature and protein designers can take advantage of promiscuous binding interactions to evolve new proteins with specialized functions.

State-dependent motion of a genetically encoded fluorescent biosensor

Genetically encoded biosensors can measure biochemical properties such as small-molecule concentrations with single-cell resolution, even in vivo. Despite their utility, these sensors are "black boxes": Very little is known about the structures of their low- and high-fluorescence states or what features are required to transition between them. We used LiLac, a lactate biosensor with a quantitative fluorescence-lifetime readout, as a model system to address these questions. X-ray crystal structures and engineered high-affinity metal bridges demonstrate that LiLac exhibits a large interdomain twist motion that pulls the fluorescent protein away from a "sealed," high-lifetime state in the absence of lactate to a "cracked," low-lifetime state in its presence. Understanding the structures and dynamics of LiLac will help to think about and engineer other fluorescent biosensors.

Optimizing Solid-Phase Protein Synthesis Using CPG-2000 and a Nickel-Cleavable SNAC-tag Linker

Solid-phase chemical ligation (SPCL) is a powerful method for simplifying protein synthesis but faces challenges with orthogonal protection strategies, suboptimal solid supports, and limited cleavable linkers. In this study, we optimized SPCL by combining low-loading controlled-pore glass (CPG-2000) with a nickel-cleavable SNAC-tag linker. This system enabled the successful assembly of five peptide fragments and the efficient synthesis of a 131-amino-acid de novo protein, achieving a 25% yield.

De Novo Design of Proteins That Bind Naphthalenediimides, Powerful Photooxidants with Tunable Photophysical Properties

De novo protein design provides a framework to test our understanding of protein function and build proteins with cofactors and functions not found in nature. Here, we report the design of proteins designed to bind powerful photooxidants and the evaluation of the use of these proteins to generate diffusible small-molecule reactive species. Because excited-state dynamics are influenced by the dynamics and hydration of a photooxidant's environment, it was important to not only design a binding site but also to evaluate its dynamic properties. Thus, we used computational design in conjunction with molecular dynamics (MD) simulations to design a protein, designated NBP (NDI Binding Protein), that held a naphthalenediimide (NDI), a powerful photooxidant, in a programmable molecular environment. Solution NMR confirmed the structure of the complex. We evaluated two NDI cofactors in this de novo protein using ultrafast pump-probe spectroscopy to evaluate light-triggered intra- and intermolecular electron transfer function. Moreover, we demonstrated the utility of this platform to activate multiple molecular probes for protein labeling.

Sequence-dependent scale for translocon-mediated insertion of interfacial helices in membranes

Biological membranes consist of a lipid bilayer studded with integral and peripheral membrane proteins. Most α-helical membrane proteins require protein-conducting insertases known as translocons to assist in their membrane insertion and folding. While the sequence-dependent propensities for a helix to either translocate through the translocon or insert into the membrane have been codified into numerical hydrophobicity scales, the corresponding propensity to partition into the membrane interface remains unrevealed. By engineering diagnostic glycosylation sites around test peptide sequences inserted into a host protein, we devised a system that can differentiate between water-soluble, surface-bound, and transmembrane (TM) states of the sequence based on its glycosylation pattern. Using this system, we determined the sequence-dependent propensities for transfer from the translocon to a TM, interfacial, or extramembrane space and compared these propensities with the corresponding probability distributions determined from the sequences and structures of experimentally determined proteins.

Confinement and Catalysis within De Novo Designed Peptide Barrels

De novo protein design has advanced such that many peptide assemblies and protein structures can be generated predictably and quickly. The drive now is to bring functions to these structures, for example, small-molecule binding and catalysis. The formidable challenge of binding and orienting multiple small molecules to direct chemistry is particularly important for paving the way to new functionalities. To address this, here we describe the design, characterization, and application of small-molecule:peptide ternary complexes in aqueous solution. This uses α-helical barrel (αHB) peptide assemblies, which comprise 5 or more α helices arranged around central channels. These channels are solvent accessible, and their internal dimensions and chemistries can be altered predictably. Thus, αHBs are analogous to "molecular flasks" made in supramolecular, polymer, and materials chemistry. Using Förster resonance energy transfer as a readout, we demonstrate that specific αHBs can accept two different organic dyes, 1,6-diphenyl-1,3,5-hexatriene and Nile red, in close proximity. In addition, two anthracene molecules can be accommodated within an αHB to promote anthracene photodimerization. However, not all ternary complexes are productive, either in energy transfer or photodimerization, illustrating the control that can be exerted by judicious choice and design of the αHB.

Transition transferases prime bacterial capsule polymerization

Capsules are long-chain carbohydrate polymers that envelop the surfaces of many bacteria, protecting them from host immune responses. Capsule biosynthesis enzymes are potential drug targets and valuable biotechnological tools for generating vaccine antigens. Despite their importance, it remains unknown how structurally variable capsule polymers of Gram-negative pathogens are linked to the conserved glycolipid anchoring these virulence factors to the bacterial membrane. Using Actinobacillus pleuropneumoniae as an example, we demonstrate that CpsA and CpsC generate a poly(glycerol-3-phosphate) linker to connect the glycolipid with capsules containing poly(galactosylglycerol-phosphate) backbones. We reconstruct the entire capsule biosynthesis pathway in A. pleuropneumoniae serotypes 3 and 7, solve the X-ray crystal structure of the capsule polymerase CpsD, identify its tetratricopeptide repeat domain as essential for elongating poly(glycerol-3-phosphate) and show that CpsA and CpsC stimulate CpsD to produce longer polymers. We identify the CpsA and CpsC product as a wall teichoic acid homolog, demonstrating similarity between the biosynthesis of Gram-positive wall teichoic acid and Gram-negative capsules.

Target-conditioned diffusion generates potent TNFR superfamily antagonists and agonists

Despite progress in designing protein-binding proteins, the shape matching of designs to targets is lower than in many native protein complexes, and design efforts have failed for the tumor necrosis factor receptor 1 (TNFR1) and other protein targets with relatively flat and polar surfaces. We hypothesized that free diffusion from random noise could generate shape-matched binders for challenging targets and tested this approach on TNFR1. We obtain designs with low picomolar affinity whose specificity can be completely switched to other family members using partial diffusion. Designs function as antagonists or as superagonists when presented at higher valency for OX40 and 4-1BB. The ability to design high-affinity and high-specificity antagonists and agonists for pharmacologically important targets in silico presages a coming era in protein design in which binders are made by computation rather than immunization or random screening approaches.

ProteinReDiff: Complex-based ligand-binding proteins redesign by equivariant diffusion-based generative models

Proteins, serving as the fundamental architects of biological processes, interact with ligands to perform a myriad of functions essential for life. Designing functional ligand-binding proteins is pivotal for advancing drug development and enhancing therapeutic efficacy. In this study, we introduce ProteinReDiff, an diffusion framework targeting the redesign of ligand-binding proteins. Using equivariant diffusion-based generative models, ProteinReDiff enables the creation of high-affinity ligand-binding proteins without the need for detailed structural information, leveraging instead the potential of initial protein sequences and ligand SMILES strings. Our evaluations across sequence diversity, structural preservation, and ligand binding affinity underscore ProteinReDiff's potential to advance computational drug discovery and protein engineering.

De novo protein design with a denoising diffusion network independent of pretrained structure prediction models

The recent success of RFdiffusion, a method for protein structure design with a denoising diffusion probabilistic model, has relied on fine-tuning the RoseTTAFold structure prediction network for protein backbone denoising. Here, we introduce SCUBA-diffusion (SCUBA-D), a protein backbone denoising diffusion probabilistic model freshly trained by considering co-diffusion of sequence representation to enhance model regularization and adversarial losses to minimize data-out-of-distribution errors. While matching the performance of the pretrained RoseTTAFold-based RFdiffusion in generating experimentally realizable protein structures, SCUBA-D readily generates protein structures with not-yet-observed overall folds that are different from those predictable with RoseTTAFold. The accuracy of SCUBA-D was confirmed by the X-ray structures of 16 designed proteins and a protein complex, and by experiments validating designed heme-binding proteins and Ras-binding proteins. Our work shows that deep generative models of images or texts can be fruitfully extended to complex physical objects like protein structures by addressing outstanding issues such as the data-out-of-distribution errors.

A designed Zn 2+ sensor domain transmits binding information to transmembrane histidine kinases

Generating stimulus-responsive, allosteric signaling de novo is a significant challenge in protein design. In natural systems like bacterial histidine kinases (HKs), signal transduction occurs when ligand binding initiates a signal that is amplified across biological membranes over long distances to induce large-scale rearrangements and phosphorylation relays. Here, we ask whether our understanding of protein design and multi-domain, intramolecular signaling has progressed sufficiently to enable engineering of a HK with tunable de novo components. We generated de novo metal-binding sensor domains and substituted them for the native sensor domain of a transmembrane HK, affording chimeras that transduce signals initiated from a de novo sensor. Signaling depended on the designed sensor's stability and the interdomain linker's phase and length. These results show the usefulness of de novo design to elucidate biochemical mechanisms and principles for design of new signaling systems.

The future of machine learning for small-molecule drug discovery will be driven by data

Many studies have prophesied that the integration of machine learning techniques into small-molecule therapeutics development will help to deliver a true leap forward in drug discovery. However, increasingly advanced algorithms and novel architectures have not always yielded substantial improvements in results. In this Perspective, we propose that a greater focus on the data for training and benchmarking these models is more likely to drive future improvement, and explore avenues for future research and strategies to address these data challenges.

Efficient Generation of Protein Pockets with PocketGen

Designing protein-binding proteins is critical for drug discovery. However, the AI-based design of such proteins is challenging due to the complexity of ligand-protein interactions, the flexibility of ligand molecules and amino acid side chains, and sequence-structure dependencies. We introduce PocketGen, a deep generative model that simultaneously produces both the residue sequence and atomic structure of the protein regions where ligand interactions occur. PocketGen ensures consistency between sequence and structure by using a graph transformer for structural encoding and a sequence refinement module based on a protein language model. The bilevel graph transformer captures interactions at multiple scales, including atom, residue, and ligand levels. To enhance sequence refinement, PocketGen integrates a structural adapter into the protein language model, ensuring that structure-based predictions align with sequence-based predictions. PocketGen can generate high-fidelity protein pockets with superior binding affinity and structural validity. It operates ten times faster than physics-based methods and achieves a 95% success rate, defined as the percentage of generated pockets with higher binding affinity than reference pockets. Additionally, it attains an amino acid recovery rate exceeding 64%.

Rationalizing Diverse Binding Mechanisms to the Same Protein Fold: Insights for Ligand Recognition and Biosensor Design

The engineering of novel protein-ligand binding interactions, particularly for complex drug-like molecules, is an unsolved problem, which could enable many practical applications of protein biosensors. In this work, we analyzed two engineered biosensors, derived from the plant hormone sensor PYR1, to recognize either the agrochemical mandipropamid or the synthetic cannabinoid WIN55,212-2. Using a combination of quantitative deep mutational scanning experiments and molecular dynamics simulations, we demonstrated that mutations at common positions can promote protein-ligand shape complementarity and revealed prominent differences in the electrostatic networks needed to complement diverse ligands. MD simulations indicate that both PYR1 protein-ligand complexes bind a single conformer of their target ligand that is close to the lowest free-energy conformer. Computational design using a fixed conformer and rigid body orientation led to new WIN55,212-2 sensors with nanomolar limits of detection. This work reveals mechanisms by which the versatile PYR1 biosensor scaffold can bind diverse ligands. This work also provides computational methods to sample realistic ligand conformers and rigid body alignments that simplify the computational design of biosensors for novel ligands of interest.

Opportunities and challenges in design and optimization of protein function

The field of protein design has made remarkable progress over the past decade. Historically, the low reliability of purely structure-based design methods limited their application, but recent strategies that combine structure-based and sequence-based calculations, as well as machine learning tools, have dramatically improved protein engineering and design. In this Review, we discuss how these methods have enabled the design of increasingly complex structures and therapeutically relevant activities. Additionally, protein optimization methods have improved the stability and activity of complex eukaryotic proteins. Thanks to their increased reliability, computational design methods have been applied to improve therapeutics and enzymes for green chemistry and have generated vaccine antigens, antivirals and drug-delivery nano-vehicles. Moreover, the high success of design methods reflects an increased understanding of basic rules that govern the relationships among protein sequence, structure and function. However, de novo design is still limited mostly to α-helix bundles, restricting its potential to generate sophisticated enzymes and diverse protein and small-molecule binders. Designing complex protein structures is a challenging but necessary next step if we are to realize our objective of generating new-to-nature activities.

Binding and sensing diverse small molecules using shape-complementary pseudocycles

We describe an approach for designing high-affinity small molecule-binding proteins poised for downstream sensing. We use deep learning-generated pseudocycles with repeating structural units surrounding central binding pockets with widely varying shapes that depend on the geometry and number of the repeat units. We dock small molecules of interest into the most shape complementary of these pseudocycles, design the interaction surfaces for high binding affinity, and experimentally screen to identify designs with the highest affinity. We obtain binders to four diverse molecules, including the polar and flexible methotrexate and thyroxine. Taking advantage of the modular repeat structure and central binding pockets, we construct chemically induced dimerization systems and low-noise nanopore sensors by splitting designs into domains that reassemble upon ligand addition.

Generalized biomolecular modeling and design with RoseTTAFold All-Atom

Deep-learning methods have revolutionized protein structure prediction and design but are presently limited to protein-only systems. We describe RoseTTAFold All-Atom (RFAA), which combines a residue-based representation of amino acids and DNA bases with an atomic representation of all other groups to model assemblies that contain proteins, nucleic acids, small molecules, metals, and covalent modifications, given their sequences and chemical structures. By fine-tuning on denoising tasks, we developed RFdiffusion All-Atom (RFdiffusionAA), which builds protein structures around small molecules. Starting from random distributions of amino acid residues surrounding target small molecules, we designed and experimentally validated, through crystallography and binding measurements, proteins that bind the cardiac disease therapeutic digoxigenin, the enzymatic cofactor heme, and the light-harvesting molecule bilin.

De novo design of drug-binding proteins with predictable binding energy and specificity

The de novo design of small molecule-binding proteins has seen exciting recent progress; however, high-affinity binding and tunable specificity typically require laborious screening and optimization after computational design. We developed a computational procedure to design a protein that recognizes a common pharmacophore in a series of poly(ADP-ribose) polymerase-1 inhibitors. One of three designed proteins bound different inhibitors with affinities ranging from <5 nM to low micromolar. X-ray crystal structures confirmed the accuracy of the designed protein-drug interactions. Molecular dynamics simulations informed the role of water in binding. Binding free energy calculations performed directly on the designed models were in excellent agreement with the experimentally measured affinities. We conclude that de novo design of high-affinity small molecule-binding proteins with tuned interaction energies is feasible entirely from computation.

De novo design of cavity-containing proteins with a backbone-centered neural network energy function

The design of small-molecule-binding proteins requires protein backbones that contain cavities. Previous design efforts were based on naturally occurring cavity-containing backbone architectures. Here, we designed diverse cavity-containing backbones without predefined architectures by introducing tailored restraints into the backbone sampling driven by SCUBA (Side Chain-Unknown Backbone Arrangement), a neural network statistical energy function. For 521 out of 5816 designs, the root-mean-square deviations (RMSDs) of the Cα atoms for the AlphaFold2-predicted structures and our designed structures are within 2.0 Å. We experimentally tested 10 designed proteins and determined the crystal structures of two of them. One closely agrees with the designed model, while the other forms a domain-swapped dimer, where the partial structures are in agreement with the designed structures. Our results indicate that data-driven methods such as SCUBA hold great potential for designing de novo proteins with tailored small-molecule-binding function.

Rationalizing diverse binding mechanisms to the same protein fold: insights for ligand recognition and biosensor design

The engineering of novel protein-ligand binding interactions, particularly for complex drug-like molecules, is an unsolved problem which could enable many practical applications of protein biosensors. In this work, we analyzed two engineer ed biosensors, derived from the plant hormone sensor PYR1, to recognize either the agrochemical mandipropamid or the synthetic cannabinoid WIN55,212-2. Using a combination of quantitative deep mutational scanning experiments and molecular dynamics simulations, we demonstrated that mutations at common positions can promote protein-ligand shape complementarity and revealed prominent differences in the electrostatic networks needed to complement diverse ligands. MD simulations indicate that both PYR1 protein-ligand complexes bind a single conformer of their target ligand that is close to the lowest free energy conformer. Computational design using a fixed conformer and rigid body orientation led to new WIN55,212-2 sensors with nanomolar limits of detection. This work reveals mechanisms by which the versatile PYR1 biosensor scaffold can bind diverse ligands. This work also provides computational methods to sample realistic ligand conformers and rigid body alignments that simplify the computational design of biosensors for novel ligands of interest.

Discovery of lupus nephritis targeted inhibitors based on De novo molecular design: comprehensive application of vinardo scoring, ADMET analysis, and molecular dynamics simulation

Lupus Nephritis (LN) is an autoimmune disease affecting the kidneys, and conventional drug studies have limitations due to its imprecise and complex pathogenesis. Therefore, the aim of this study was to design a novel Lupus Nephritis-targeted drug with good clinical due potential, high potency and selectivity by computer-assisted approach.NIK belongs to the serine/threonine protein kinase, which is gaining attention as a drug target for Lupus Nephritis. we used bioinformatics, homology modelling and sequence comparison analysis, small molecule ab initio design, ADMET analysis, molecular docking, molecular dynamics simulation, and MM/PBSA analysis to design and explore the selectivity and efficiency of a novel Lupus Nephritis-targeting drug, ClImYnib, and a classical NIK inhibitor, NIK SMI1. We used bioinformatics techniques to determine the correlation between lupus nephritis and the NF-κB signaling pathway. De novo drugs design was used to create a NIK-targeted inhibitor, ClImYnib, with lower toxicity, after which we used molecular dynamics to simulate NIK SMI1 against ClImYnib, and the simulation results showed that ClImYnib had better selectivity and efficiency. Our research delves into the molecular mechanism of protein ligands, and we have designed and validated an excellent NIK inhibitor using multiple computational simulation methods. More importantly, it provides an idea of target designing small molecules.

Deep Confident Steps to New Pockets: Strategies for Docking Generalization

Accurate blind docking has the potential to lead to new biological breakthroughs, but for this promise to be realized, docking methods must generalize well across the proteome. Existing benchmarks, however, fail to rigorously assess generalizability. Therefore, we develop DockGen, a new benchmark based on the ligand-binding domains of proteins, and we show that existing machine learning-based docking models have very weak generalization abilities. We carefully analyze the scaling laws of ML-based docking and show that, by scaling data and model size, as well as integrating synthetic data strategies, we are able to significantly increase the generalization capacity and set new state-of-the-art performance across benchmarks. Further, we propose Confidence Bootstrapping, a new training paradigm that solely relies on the interaction between diffusion and confidence models and exploits the multi-resolution generation process of diffusion models. We demonstrate that Confidence Bootstrapping significantly improves the ability of ML-based docking methods to dock to unseen protein classes, edging closer to accurate and generalizable blind docking methods.

Molecular co-assembled strategy tuning protein conformation for cartilage regeneration

The assembly of oligopeptide and polypeptide molecules can reconstruct various ordered advanced structures through intermolecular interactions to achieve protein-like biofunction. Here, we develop a "molecular velcro"-inspired peptide and gelatin co-assembly strategy, in which amphiphilic supramolecular tripeptides are attached to the molecular chain of gelatin methacryloyl via intra-/intermolecular interactions. We perform molecular docking and dynamics simulations to demonstrate the feasibility of this strategy and reveal the advanced structural transition of the co-assembled hydrogel, which brings more ordered β-sheet content and 10-fold or more compressive strength improvement. We conduct transcriptome analysis to reveal the role of co-assembled hydrogel in promoting cell proliferation and chondrogenic differentiation. Subcutaneous implantation evaluation confirms considerably reduced inflammatory responses and immunogenicity in comparison with type I collagen. We demonstrate that bone mesenchymal stem cells-laden co-assembled hydrogel can be stably fixed in rabbit knee joint defects by photocuring, which significantly facilitates hyaline cartilage regeneration after three months. This co-assembly strategy provides an approach for developing cartilage regenerative biomaterials.

Sparks of function by de novo protein design

Information in proteins flows from sequence to structure to function, with each step causally driven by the preceding one. Protein design is founded on inverting this process: specify a desired function, design a structure executing this function, and find a sequence that folds into this structure. This 'central dogma' underlies nearly all de novo protein-design efforts. Our ability to accomplish these tasks depends on our understanding of protein folding and function and our ability to capture this understanding in computational methods. In recent years, deep learning-derived approaches for efficient and accurate structure modeling and enrichment of successful designs have enabled progression beyond the design of protein structures and towards the design of functional proteins. We examine these advances in the broader context of classical de novo protein design and consider implications for future challenges to come, including fundamental capabilities such as sequence and structure co-design and conformational control considering flexibility, and functional objectives such as antibody and enzyme design.

De novo protein design-From new structures to programmable functions

Methods from artificial intelligence (AI) trained on large datasets of sequences and structures can now "write" proteins with new shapes and molecular functions de novo, without starting from proteins found in nature. In this Perspective, I will discuss the state of the field of de novo protein design at the juncture of physics-based modeling approaches and AI. New protein folds and higher-order assemblies can be designed with considerable experimental success rates, and difficult problems requiring tunable control over protein conformations and precise shape complementarity for molecular recognition are coming into reach. Emerging approaches incorporate engineering principles-tunability, controllability, and modularity-into the design process from the beginning. Exciting frontiers lie in deconstructing cellular functions with de novo proteins and, conversely, constructing synthetic cellular signaling from the ground up. As methods improve, many more challenges are unsolved.

Protein sequence design on given backbones with deep learning

Deep learning methods for protein sequence design focus on modeling and sampling the many- dimensional distribution of amino acid sequences conditioned on the backbone structure. To produce physically foldable sequences, inter-residue couplings need to be considered properly. These couplings are treated explicitly in iterative methods or autoregressive methods. Non-autoregressive models treating these couplings implicitly are computationally more efficient, but still await tests by wet experiment. Currently, sequence design methods are evaluated mainly using native sequence recovery rate and native sequence perplexity. These metrics can be complemented by sequence-structure compatibility metrics obtained from energy calculation or structure prediction. However, existing computational metrics have important limitations that may render the generalization of computational test results to performance in real applications unwarranted. Validation of design methods by wet experiments should be encouraged.

De novo design of drug-binding proteins with predictable binding energy and specificity

The de novo design of small-molecule-binding proteins has seen exciting recent progress; however, the ability to achieve exquisite affinity for binding small molecules while tuning specificity has not yet been demonstrated directly from computation. Here, we develop a computational procedure that results in the highest affinity binders to date with predetermined relative affinities, targeting a series of PARP1 inhibitors. Two of four designed proteins bound with affinities ranging from < 5 nM to low μM, in a predictable manner. X-ray crystal structures confirmed the accuracy of the designed protein-drug interactions. Molecular dynamics simulations informed the role of water in binding. Binding free-energy calculations performed directly on the designed models are in excellent agreement with the experimentally measured affinities, suggesting that the de novo design of small-molecule-binding proteins with tuned interaction energies is now feasible entirely from computation. We expect these methods to open many opportunities in biomedicine, including rapid sensor development, antidote design, and drug delivery vehicles.

De novo design of diverse small molecule binders and sensors using Shape Complementary Pseudocycles

A general method for designing proteins to bind and sense any small molecule of interest would be widely useful. Due to the small number of atoms to interact with, binding to small molecules with high affinity requires highly shape complementary pockets, and transducing binding events into signals is challenging. Here we describe an integrated deep learning and energy based approach for designing high shape complementarity binders to small molecules that are poised for downstream sensing applications. We employ deep learning generated psuedocycles with repeating structural units surrounding central pockets; depending on the geometry of the structural unit and repeat number, these pockets span wide ranges of sizes and shapes. For a small molecule target of interest, we extensively sample high shape complementarity pseudocycles to generate large numbers of customized potential binding pockets; the ligand binding poses and the interacting interfaces are then optimized for high affinity binding. We computationally design binders to four diverse molecules, including for the first time polar flexible molecules such as methotrexate and thyroxine, which are expressed at high levels and have nanomolar affinities straight out of the computer. Co-crystal structures are nearly identical to the design models. Taking advantage of the modular repeating structure of pseudocycles and central location of the binding pockets, we constructed low noise nanopore sensors and chemically induced dimerization systems by splitting the binders into domains which assemble into the original pseudocycle pocket upon target molecule addition.

Small-molecule binding and sensing with a designed protein family

Despite transformative advances in protein design with deep learning, the design of small-molecule-binding proteins and sensors for arbitrary ligands remains a grand challenge. Here we combine deep learning and physics-based methods to generate a family of proteins with diverse and designable pocket geometries, which we employ to computationally design binders for six chemically and structurally distinct small-molecule targets. Biophysical characterization of the designed binders revealed nanomolar to low micromolar binding affinities and atomic-level design accuracy. The bound ligands are exposed at one edge of the binding pocket, enabling the de novo design of chemically induced dimerization (CID) systems; we take advantage of this to create a biosensor with nanomolar sensitivity for cortisol. Our approach provides a general method to design proteins that bind and sense small molecules for a wide range of analytical, environmental, and biomedical applications.

Hallucination of closed repeat proteins containing central pockets

In pseudocyclic proteins, such as TIM barrels, β barrels, and some helical transmembrane channels, a single subunit is repeated in a cyclic pattern, giving rise to a central cavity that can serve as a pocket for ligand binding or enzymatic activity. Inspired by these proteins, we devised a deep-learning-based approach to broadly exploring the space of closed repeat proteins starting from only a specification of the repeat number and length. Biophysical data for 38 structurally diverse pseudocyclic designs produced in Escherichia coli are consistent with the design models, and the three crystal structures we were able to obtain are very close to the designed structures. Docking studies suggest the diversity of folds and central pockets provide effective starting points for designing small-molecule binders and enzymes.

An expandable, modular de novo protein platform for precision redox engineering

The electron-conducting circuitry of life represents an as-yet untapped resource of exquisite, nanoscale biomolecular engineering. Here, we report the characterization and structure of a de novo diheme "maquette" protein, 4D2, which we subsequently use to create an expanded, modular platform for heme protein design. A well-folded monoheme variant was created by computational redesign, which was then utilized for the experimental validation of continuum electrostatic redox potential calculations. This demonstrates how fundamental biophysical properties can be predicted and fine-tuned. 4D2 was then extended into a tetraheme helical bundle, representing a 7 nm molecular wire. Despite a molecular weight of only 24 kDa, electron cryomicroscopy illustrated a remarkable level of detail, indicating the positioning of the secondary structure and the heme cofactors. This robust, expressible, highly thermostable and readily designable modular platform presents a valuable resource for redox protein design and the future construction of artificial electron-conducting circuitry.

Transient water wires mediate selective proton transport in designed channel proteins

Selective proton transport through proteins is essential for forming and using proton gradients in cells. Protons are conducted along hydrogen-bonded 'wires' of water molecules and polar side chains, which, somewhat surprisingly, are often interrupted by dry apolar stretches in the conduction pathways, inferred from static protein structures. Here we hypothesize that protons are conducted through such dry spots by forming transient water wires, often highly correlated with the presence of the excess protons in the water wire. To test this hypothesis, we performed molecular dynamics simulations to design transmembrane channels with stable water pockets interspersed by apolar segments capable of forming flickering water wires. The minimalist designed channels conduct protons at rates similar to viral proton channels, and they are at least 106-fold more selective for H+ over Na+. These studies inform the mechanisms of biological proton conduction and the principles for engineering proton-conductive materials.

Design and Selection of Heterodimerizing Helical Hairpins for Synthetic Biology

Synthetic biology applications would benefit from protein modules of reduced complexity that function orthogonally to cellular components. As many subcellular processes depend on peptide-protein or protein-protein interactions, de novo designed polypeptides that can bring together other proteins controllably are particularly useful. Thanks to established sequence-to-structure relationships, helical bundles provide good starting points for such designs. Typically, however, such designs are tested in vitro and function in cells is not guaranteed. Here, we describe the design, characterization, and application of de novo helical hairpins that heterodimerize to form 4-helix bundles in cells. Starting from a rationally designed homodimer, we construct a library of helical hairpins and identify complementary pairs using bimolecular fluorescence complementation in E. coli. We characterize some of the pairs using biophysics and X-ray crystallography to confirm heterodimeric 4-helix bundles. Finally, we demonstrate the function of an exemplar pair in regulating transcription in both E. coli and mammalian cells.

Computational structural-based GPCR optimization for user-defined ligand: Implications for the development of biosensors

Organisms have developed effective mechanisms to sense the external environment. Human-designed biosensors exploit this natural optimization, where different biological machinery have been adapted to detect the presence of user-defined molecules. Specifically, the pheromone pathway in the model organism Saccharomyces cerevisiae represents a suitable candidate as a synthetic signaling system. Indeed, it expresses just one G-Protein Coupled Receptor (GPCR), Ste2, able to recognize pheromone and initiate the expression of pheromone-dependent genes. To date, the standard procedure to engineer this system relies on the substitution of the yeast GPCR with another one and on the modification of the yeast G-protein to bind the inserted receptor. Here, we propose an innovative computational procedure, based on geometrical and chemical optimization of protein binding pockets, to select the amino acid substitutions required to make the native yeast GPCR able to recognize a user-defined ligand. This procedure would allow the yeast to recognize a wide range of ligands, without a-priori knowledge about a GPCR recognizing them or the corresponding G protein. We used Monte Carlo simulations to design on Ste2 a binding pocket able to recognize epinephrine, selected as a test ligand. We validated Ste2 mutants via molecular docking and molecular dynamics. We verified that the amino acid substitutions we identified make Ste2 able to accommodate and remain firmly bound to epinephrine. Our results indicate that we sampled efficiently the huge space of possible mutants, proposing such a strategy as a promising starting point for the development of a new kind of S.cerevisiae-based biosensors.

Atligator Web: A Graphical User Interface for Analysis and Design of Protein-Peptide Interactions

A key functionality of proteins is based on their ability to form interactions with other proteins or peptides. These interactions are neither easily described nor fully understood, which is why the design of specific protein-protein interaction interfaces remains a challenge for protein engineering. We recently developed the software ATLIGATOR to extract common interaction patterns between different types of amino acids and store them in a database. The tool enables the user to better understand frequent interaction patterns and find groups of interactions. Furthermore, frequent motifs can be directly transferred from the database to a user-defined scaffold as a starting point for the engineering of new binding capabilities. Since three-dimensional visualization is a crucial part of ATLIGATOR, we created ATLIGATOR web-a web server offering an intuitive graphical user interface (GUI) available at https://atligator.uni-bayreuth.de. This new interface empowers users to apply ATLIGATOR by providing easy access with having all parts directly connected. Moreover, we extended the web by a design functionality so that, overall, ATLIGATOR web facilitates the use of ATLIGATOR with a more intuitive UI and advanced design options.

Transcription factor-based biosensors for screening and dynamic regulation

Advances in synthetic biology and genetic engineering are bringing into the spotlight a wide range of bio-based applications that demand better sensing and control of biological behaviours. Transcription factor (TF)-based biosensors are promising tools that can be used to detect several types of chemical compounds and elicit a response according to the desired application. However, the wider use of this type of device is still hindered by several challenges, which can be addressed by increasing the current metabolite-activated transcription factor knowledge base, developing better methods to identify new transcription factors, and improving the overall workflow for the design of novel biosensor circuits. These improvements are particularly important in the bioproduction field, where researchers need better biosensor-based approaches for screening production-strains and precise dynamic regulation strategies. In this work, we summarize what is currently known about transcription factor-based biosensors, discuss recent experimental and computational approaches targeted at their modification and improvement, and suggest possible future research directions based on two applications: bioproduction screening and dynamic regulation of genetic circuits.

Rapid biosensor development using plant hormone receptors as reprogrammable scaffolds

A general method to generate biosensors for user-defined molecules could provide detection tools for a wide range of biological applications. Here, we describe an approach for the rapid engineering of biosensors using PYR1 (Pyrabactin Resistance 1), a plant abscisic acid (ABA) receptor with a malleable ligand-binding pocket and a requirement for ligand-induced heterodimerization, which facilitates the construction of sense-response functions. We applied this platform to evolve 21 sensors with nanomolar to micromolar sensitivities for a range of small molecules, including structurally diverse natural and synthetic cannabinoids and several organophosphates. X-ray crystallography analysis revealed the mechanistic basis for new ligand recognition by an evolved cannabinoid receptor. We demonstrate that PYR1-derived receptors are readily ported to various ligand-responsive outputs, including enzyme-linked immunosorbent assay (ELISA)-like assays, luminescence by protein-fragment complementation and transcriptional circuits, all with picomolar to nanomolar sensitivity. PYR1 provides a scaffold for rapidly evolving new biosensors for diverse sense-response applications.

A Structural Analysis of Proteinaceous Nanotube Cavities and Their Applications in Nanotechnology

Protein nanotubes offer unique properties to the materials science field that allow them to fulfill various functions in drug delivery, biosensors and energy storage. Protein nanotubes are chemically diverse, modular, biodegradable and nontoxic. Furthermore, although the initial design or repurposing of such nanotubes is highly complex, the field has matured to understand underlying chemical and physical properties to a point where applications are successfully being developed. An important feature of a nanotube is its ability to bind ligands via its internal cavities. As ligands of interest vary in size, shape and chemical properties, cavities have to be able to accommodate very specific features. As such, understanding cavities on a structural level is essential for their effective application. The objective of this review is to present the chemical and physical diversity of protein nanotube cavities and highlight their potential applications in materials science, specifically in biotechnology.

ATLIGATOR: editing protein interactions with an atlas-based approach

MOTIVATION

Recognition of specific molecules by proteins is a fundamental cellular mechanism and relevant for many applications. Being able to modify binding is a key interest and can be achieved by repurposing established interaction motifs. We were specifically interested in a methodology for the design of peptide binding modules. By leveraging interaction data from known protein structures, we plan to accelerate the design of novel protein or peptide binders.

RESULTS

We developed ATLIGATOR-a computational method to support the analysis and design of a protein's interaction with a single side chain. Our program enables the building of interaction atlases based on structures from the PDB. From these atlases pocket definitions are extracted that can be searched for frequent interactions. These searches can reveal similarities in unrelated proteins as we show here for one example. Such frequent interactions can then be grafted onto a new protein scaffold as a starting point of the design process. The ATLIGATOR tool is made accessible through a python API as well as a CLI with python scripts.

AVAILABILITY AND IMPLEMENTATION

Source code can be downloaded at github (https://www.github.com/Hoecker-Lab/atligator), installed from PyPI ('atligator') and is implemented in Python 3.

From peptides to proteins: coiled-coil tetramers to single-chain 4-helix bundles

The design of completely synthetic proteins from first principles-de novo protein design-is challenging. This is because, despite recent advances in computational protein-structure prediction and design, we do not understand fully the sequence-to-structure relationships for protein folding, assembly, and stabilization. Antiparallel 4-helix bundles are amongst the most studied scaffolds for de novo protein design. We set out to re-examine this target, and to determine clear sequence-to-structure relationships, or design rules, for the structure. Our aim was to determine a common and robust sequence background for designing multiple de novo 4-helix bundles. In turn, this could be used in chemical and synthetic biology to direct protein-protein interactions and as scaffolds for functional protein design. Our approach starts by analyzing known antiparallel 4-helix coiled-coil structures to deduce design rules. In terms of the heptad repeat, abcdefg -i.e., the sequence signature of many helical bundles-the key features that we identify are: a = Leu, d = Ile, e = Ala, g = Gln, and the use of complementary charged residues at b and c. Next, we implement these rules in the rational design of synthetic peptides to form antiparallel homo- and heterotetramers. Finally, we use the sequence of the homotetramer to derive in one step a single-chain 4-helix-bundle protein for recombinant production in E. coli. All of the assembled designs are confirmed in aqueous solution using biophysical methods, and ultimately by determining high-resolution X-ray crystal structures. Our route from peptides to proteins provides an understanding of the role of each residue in each design.

Computational design of a sensitive, selective phase-changing sensor protein for the VX nerve agent

The VX nerve agent is one of the deadliest chemical warfare agents. Specific, sensitive, real-time detection methods for this neurotoxin have not been reported. The creation of proteins that use biological recognition to fulfill these requirements using directed evolution or library screening methods has been hampered because its toxicity makes laboratory experimentation extraordinarily expensive. A pair of VX-binding proteins were designed using a supercharged scaffold that couples a large-scale phase change from unstructured to folded upon ligand binding, enabling fully internal binding sites that present the maximum surface area possible for high affinity and specificity in target recognition. Binding site residues were chosen using a new distributed evolutionary algorithm implementation in protCAD. Both designs detect VX at parts per billion concentrations with high specificity. Computational design of fully buried molecular recognition sites, in combination with supercharged phase-changing chassis proteins, enables the ready development of a new generation of small-molecule biosensors.

Protein design via deep learning

Proteins with desired functions and properties are important in fields like nanotechnology and biomedicine. De novo protein design enables the production of previously unseen proteins from the ground up and is believed as a key point for handling real social challenges. Recent introduction of deep learning into design methods exhibits a transformative influence and is expected to represent a promising and exciting future direction. In this review, we retrospect the major aspects of current advances in deep-learning-based design procedures and illustrate their novelty in comparison with conventional knowledge-based approaches through noticeable cases. We not only describe deep learning developments in structure-based protein design and direct sequence design, but also highlight recent applications of deep reinforcement learning in protein design. The future perspectives on design goals, challenges and opportunities are also comprehensively discussed.

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.

Accelerating therapeutic protein design

Protein structures provide for defined microenvironments that can support complex pharmacological functions, otherwise unachievable by small molecules. The advent of therapeutic proteins has thus greatly broadened the range of manageable disorders. Leveraging the knowledge and recent advances in de novo protein design methods has the prospect of revolutionizing how protein drugs are discovered and developed. This review lays out the main challenges facing therapeutic proteins discovery and development, and how present and future advancements of protein design can accelerate the protein drug pipelines.

The stability and dynamics of computationally designed proteins

Protein stability, dynamics and function are intricately linked. Accordingly, protein designers leverage dynamics in their designs and gain insight to their successes and failures by analyzing their proteins' dynamics. Molecular dynamics (MD) simulations are a powerful computational tool for quantifying both local and global protein dynamics. This review highlights studies where MD simulations were applied to characterize the stability and dynamics of designed proteins and where dynamics were incorporated into computational protein design. First, we discuss the structural basis underlying the extreme stability and thermostability frequently observed in computationally designed proteins. Next, we discuss examples of designed proteins, where dynamics were not explicitly accounted for in the design process, whose coordinated motions or active site dynamics, as observed by MD simulation, enhanced or detracted from their function. Many protein functions depend on sizeable or subtle conformational changes, so we finally discuss the computational design of proteins to perform a specific function that requires consideration of motion by multi-state design.

DEPACT and PACMatch: A Workflow of Designing De Novo Protein Pockets to Bind Small Molecules

Engineering of new functional proteins such as enzymes and biosensors involves the design of new protein pockets for the specific binding of small molecules. Here, we report a workflow composed of two new computational methods to execute this task. The DEPACT (Design Pocket as a Cluster based on Templates) method is a data-driven approach to design and evaluate small-molecule-binding pockets as isolated clusters, while the PACMatch method is a computational approach to match pocket residues in a cluster model to positions on given protein scaffolds. Using DEPACT and its scoring function, pocket clusters of natural-pocket-like chemical compositions and protein-ligand interaction strength can be designed. DEPACT can design pocket clusters containing water- or metal-ion-mediated protein-ligand interactions. While being able to efficiently treat relatively large pocket cluster models (e.g., of around 10 pocket residues), PACMatch outperforms previous methods in test cases of recovering the native positions of pocket residues in natural enzyme-substrate complexes.

Transferability of Geometric Patterns from Protein Self-Interactions to Protein-Ligand Interactions

There is significant interest in developing machine learning methods to model protein-ligand interactions but a scarcity of experimentally resolved protein-ligand structures to learn from. Protein self-contacts are a much larger source of structural data that could be leveraged, but currently it is not well understood how this data source differs from the target domain. Here, we characterize the 3D geometric patterns of protein self-contacts as probability distributions. We then present a flexible statistical framework to assess the transferability of these patterns to protein-ligand contacts. We observe that the level of transferability from protein self-contacts to protein-ligand contacts depends on contact type, with many contact types exhibiting high transferability. We then demonstrate the potential of leveraging information from these geometric patterns to aid in ligand pose-selection problems in protein-ligand docking. We publicly release our extracted data on geometric interaction patterns to enable further exploration of this problem.