Validation of de novo designed water-soluble and transmembrane ß-barrels by in silico folding and melting.

In silico validation of de novo designed proteins with deep learning (DL)-based structure prediction algorithms has become mainstream. However, formal evidence of the relationship between a high-quality predicted model and the chance of experimental success is lacking. We used experimentally characterized de novo water-soluble and transmembrane ß-barrel designs to show that AlphaFold2 and ESMFold excel at different tasks. ESMFold can efficiently identify designs generated based on high-quality (designable) backbones. However, only AlphaFold2 can predict which sequences have the best chance of experimentally folding among similar designs. We show that ESMFold can generate high-quality structures from just a few predicted contacts and introduce a new approach based on incremental perturbation of the prediction ("in silico melting"), which can reveal differences in the presence of favorable contacts between designs. This study provides a new insight on DL-based structure prediction models explainability and on how they could be leveraged for the design of increasingly complex proteins; in particular membrane proteins which have historically lacked basic in silico validation tools.

Rationale in Custom Design of Transmembrane ß-Barrel Pores.

Biological nanopores incorporated into synthetic membranes are widely used for single-molecule analytical applications such as DNA sequencing. The ability to engineer custom membrane proteins with a pore would allow the generation of a multitude of nanopores for the sensing/sequencing of small molecules and (bio)polymers. The de novo design of transmembrane ß-barrel pores has recently enabled the generation of nanopores with custom size, shape, and properties. In this chapter, I describe the rationale of transmembrane ß-barrel design and computational methods to assemble the backbones, design sequences, and select the designs for experimental validation.

ProteinMPNN Recovers Complex Sequence Properties of Transmembrane ß-barrels.

Recent deep-learning (DL) protein design methods have been successfully applied to a range of protein design problems, including the de novo design of novel folds, protein binders, and enzymes. However, DL methods have yet to meet the challenge of de novo membrane protein (MP) and the design of complex ß-sheet folds. We performed a comprehensive benchmark of one DL protein sequence design method, ProteinMPNN, using transmembrane and water-soluble ß-barrel folds as a model, and compared the performance of ProteinMPNN to the new membrane-specific Rosetta Franklin2023 energy function. We tested the effect of input backbone refinement on ProteinMPNN performance and found that given refined and well-defined inputs, ProteinMPNN more accurately captures global sequence properties despite complex folding biophysics. It generates more diverse TMB sequences than Franklin2023 in pore-facing positions. In addition, ProteinMPNN generated TMB sequences that passed state-of-the-art in silico filters for experimental validation, suggesting that the model could be used in de novo design tasks of diverse nanopores for single-molecule sensing and sequencing. Lastly, our results indicate that the low success rate of ProteinMPNN for the design of ß-sheet proteins stems from backbone input accuracy rather than software limitations.

Treponema pallidum subsp. pallidum with an Artificially impaired TprK antigenic variation system is attenuated in the Rabbit model of syphilis.

BACKGROUND: The TprK protein of the syphilis agent, Treponema pallidum subsp. pallidum (T. pallidum), undergoes antigenic variation in seven discrete variable (V) regions via non-reciprocal segmental gene conversion. These recombination events transfer information from a repertoire of 53 silent chromosomal donor cassettes (DCs) into the single tprK expression site to continually generate TprK variants. Several lines of research developed over the last two decades support the theory that this mechanism is central to T. pallidum's ability for immune avoidance and persistence in the host. Structural and modeling data, for example, identify TprK as an integral outer membrane porin with the V regions exposed on the pathogen's surface. Furthermore, infection-induced antibodies preferentially target the V regions rather than the predicted ß-barrel scaffolding, and sequence variation abrogates the binding of antibodies elicited by antigenically different V regions. Here, we engineered a T. pallidum strain to impair its ability to vary TprK and assessed its virulence in the rabbit model of syphilis. PRINCIPAL FINDINGS: A suicide vector was transformed into the wild-type (WT) SS14 T. pallidum isolate to eliminate 96% of its tprK DCs. The resulting SS14-DCKO strain exhibited an in vitro growth rate identical to the untransformed strain, supporting that the elimination of the DCs did not affect strain viability in absence of immune pressure. In rabbits injected intradermally with the SS14-DCKO strain, generation of new TprK sequences was impaired, and the animals developed attenuated lesions with a significantly reduced treponemal burden compared to control animals. During infection, clearance of V region variants originally in the inoculum mirrored the generation of antibodies to these variants, although no new variants were generated in the SS14-DCKO strain to overcome immune pressure. Naïve rabbits that received lymph node extracts from animals infected with the SS14-DCKO strain remained uninfected. CONCLUSION: These data further support the critical role of TprK in T. pallidum virulence and persistence during infection.

Treponema pallidum subsp. pallidum with an Artificially Impaired TprK Antigenic Variation System is Attenuated in the Rabbit Model of Syphilis.

Background: The TprK protein of the syphilis agent, Treponema pallidum subsp. pallidum ( T. pallidum ), undergoes antigenic variation in seven discrete variable (V) regions via non-reciprocal segmental gene conversion. These recombination events transfer information from a repertoire of 53 silent chromosomal donor cassettes (DCs) into the single tprK expression site to continually generate TprK variants. Several lines of research developed over the last two decades support the theory that this mechanism is central to T. pallidum 's ability for immune avoidance and persistence in the host. Structural and modeling data, for example, identify TprK as an integral outer membrane porin with the V regions exposed on the pathogen's surface. Furthermore, infection-induced antibodies preferentially target the V regions rather than the predicted ß-barrel scaffolding, and sequence variation abrogates the binding of antibodies elicited by antigenically different V regions. Here, we engineered a T. pallidum strain to impair its ability to vary TprK and assessed its virulence in the rabbit model of syphilis. Principal findings: A suicide vector was transformed into the wild-type (WT) SS14 T. pallidum isolate to eliminate 96% of its tprK DCs. The resulting SS14-DC KO strain exhibited an in vitro growth rate identical to the untransformed strain, supporting that the elimination of the DCs did not affect strain viability in absence of immune pressure. In rabbits injected intradermally with the SS14-DC KO strain, generation of new TprK sequences was impaired, and the animals developed attenuated lesions with a significantly reduced treponemal burden compared to control animals. During infection, clearance of V region variants originally in the inoculum mirrored the generation of antibodies to these variants, although no new variants were generated in the SS14-DC KO strain to overcome immune pressure. Naïve rabbits that received lymph node extracts from animals infected with the SS14-DC KO strain remained uninfected. Conclusion: These data further support the critical role of TprK in T. pallidum virulence and persistence during infection. Author Summary: Syphilis is still endemic in low- and middle-income countries, and it has been resurgent in high-income nations, including the U.S., for years. In endemic areas, there is still significant morbidity and mortality associated with this disease, particularly when its causative agent, the spirochete Treponema pallidum subsp . pallidum ( T. pallidum ) infects the fetus during pregnancy. Improving our understanding of syphilis pathogenesis and T. pallidum biology could help investigators devise better control strategies for this serious infection. Now that tools to genetically manipulate this pathogen are available, we can engineer T. pallidum strains lacking specific genes or genomic regions known (or believed) to be associated with virulence. This approach can shed light on the role of the ablated genes or sequences in disease development using loss-of-function strains. Here, we derived a knockout (KO) T. pallidum mutant (SS14-DC KO ) impaired in its ability to undergo antigenic variation of TprK, a protein that has long been hypothesized to be central in evasion of the host immune response and pathogen persistence during infection. When compared to the WT isolate, which is still capable of antigenic variation, the SS14-DC KO strain is significantly attenuated in its ability to proliferate and to induce early disease manifestations in infected rabbits. Our results further support the importance of TprK antigenic variation in syphilis pathogenesis and pathogen persistence.

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.

Sculpting conducting nanopore size and shape through de novo protein design.

Transmembrane ß-barrels (TMBs) are widely used for single molecule DNA and RNA sequencing and have considerable potential for a broad range of sensing and sequencing applications. Current engineering approaches for nanopore sensors are limited to naturally occurring channels such as CsgG, which have evolved to carry out functions very different from sensing, and hence provide sub-optimal starting points. In contrast, de novo protein design can in principle create an unlimited number of new nanopores with any desired properties. Here we describe a general approach to the design of transmembrane ß-barrel pores with different diameter and pore geometry. NMR and crystallographic characterization shows that the designs are stably folded with structures close to the design models. We report the first examples of de novo designed TMBs with 10, 12 and 14 stranded ß-barrels. The designs have distinct conductances that correlate with their pore diameter, ranging from 110 pS (~0.5 nm pore diameter) to 430 pS (~1.1 nm pore diameter), and can be converted into sensitive small-molecule sensors with high signal to noise ratio. The capability to generate on demand ß-barrel pores of defined geometry opens up fundamentally new opportunities for custom engineering of sequencing and sensing technologies.

Design and optimization of enzymatic activity in a de novo ß-barrel scaffold.

While native scaffolds offer a large diversity of shapes and topologies for enzyme engineering, their often unpredictable behavior in response to sequence modification makes de novo generated scaffolds an exciting alternative. Here we explore the customization of the backbone and sequence of a de novo designed eight stranded ß-barrel protein to create catalysts for a retro-aldolase model reaction. We show that active and specific catalysts can be designed in this fold and use directed evolution to further optimize activity and stereoselectivity. Our results support previous suggestions that different folds have different inherent amenability to evolution and this property could account, in part, for the distribution of natural enzymes among different folds.

Evolutionary adaptation of the protein folding pathway for secretability.

Secretory preproteins of the Sec pathway are targeted post-translationally and cross cellular membranes through translocases. During cytoplasmic transit, mature domains remain non-folded for translocase recognition/translocation. After translocation and signal peptide cleavage, mature domains fold to native states in the bacterial periplasm or traffic further. We sought the structural basis for delayed mature domain folding and how signal peptides regulate it. We compared how evolution diversified a periplasmic peptidyl-prolyl isomerase PpiA mature domain from its structural cytoplasmic PpiB twin. Global and local hydrogen-deuterium exchange mass spectrometry showed that PpiA is a slower folder. We defined at near-residue resolution hierarchical folding initiated by similar foldons in the twins, at different order and rates. PpiA folding is delayed by less hydrophobic native contacts, frustrated residues and a ß-turn in the earliest foldon and by signal peptide-mediated disruption of foldon hierarchy. When selected PpiA residues and/or its signal peptide were grafted onto PpiB, they converted it into a slow folder with enhanced in vivo secretion. These structural adaptations in a secretory protein facilitate trafficking.

B-Cell Epitope Mapping of TprC and TprD Variants of Treponema pallidum Subspecies Informs Vaccine Development for Human Treponematoses.

Several recent studies have focused on the identification, functional analysis, and structural characterization of outer membrane proteins (OMPs) of Treponema pallidum (Tp). The Tp species encompasses the highly related pallidum, pertenue, and endemicum subspecies of this pathogen, known to be the causative agents of syphilis, yaws, and bejel, respectively. These studies highlighted the importance of identifying surface-exposed OMP regions and the identification of B-cell epitopes that could be protective and used in vaccine development efforts. We previously reported that the TprC and TprD OMPs of Tp are predicted to contain external loops scattered throughout the entire length of the proteins, several of which show a low degree of sequence variability among strains and subspecies. In this study, these models were corroborated using AlphaFold2, a state-of-the-art protein structure modeling software. Here, we identified B-cell epitopes across the full-length TprC and TprD variants using the Geysan pepscan mapping approach with antisera from rabbits infected with syphilis, yaws, and bejel strains and from animals immunized with refolded recombinant TprC proteins from three syphilis strains. Our results show that the humoral response is primarily directed to sequences predicted to be on surface-exposed loops of TprC and TprD proteins, and that the magnitude of the humoral response to individual epitopes differs among animals infected with various syphilis strains and Tp subspecies. Rather than exhibiting strain-specificity, antisera showed various degrees of cross-reactivity with variant sequences from other strains. The data support the further exploration of TprC and TprD as vaccine candidates.

Principles and Methods in Computational Membrane Protein Design.

After decades of progress in computational protein design, the design of proteins folding and functioning in lipid membranes appears today as the next frontier. Some notable successes in the de novo design of simplified model membrane protein systems have helped articulate fundamental principles of protein folding, architecture and interaction in the hydrophobic lipid environment. These principles are reviewed here, together with the computational methods and approaches that were used to identify them. We provide an overview of the methodological innovations in the generation of new protein structures and functions and in the development of membrane-specific energy functions. We highlight the opportunities offered by new machine learning approaches applied to protein design, and by new experimental characterization techniques applied to membrane proteins. Although membrane protein design is in its infancy, it appears more reachable than previously thought.

De novo design of transmembrane ß barrels.

Transmembrane ß-barrel proteins (TMBs) are of great interest for single-molecule analytical technologies because they can spontaneously fold and insert into membranes and form stable pores, but the range of pore properties that can be achieved by repurposing natural TMBs is limited. We leverage the power of de novo computational design coupled with a "hypothesis, design, and test" approach to determine TMB design principles, notably, the importance of negative design to slow ß-sheet assembly. We design new eight-stranded TMBs, with no homology to known TMBs, that insert and fold reversibly into synthetic lipid membranes and have nuclear magnetic resonance and x-ray crystal structures very similar to the computational models. These advances should enable the custom design of pores for a wide range of applications.

Incorporation of sensing modalities into de novo designed fluorescence-activating proteins.

Through the efforts of many groups, a wide range of fluorescent protein reporters and sensors based on green fluorescent protein and its relatives have been engineered in recent years. Here we explore the incorporation of sensing modalities into de novo designed fluorescence-activating proteins, called mini-fluorescence-activating proteins (mFAPs), that bind and stabilize the fluorescent cis-planar state of the fluorogenic compound DFHBI. We show through further design that the fluorescence intensity and specificity of mFAPs for different chromophores can be tuned, and the fluorescence made sensitive to pH and Ca2+ for real-time fluorescence reporting. Bipartite split mFAPs enable real-time monitoring of protein-protein association and (unlike widely used split GFP reporter systems) are fully reversible, allowing direct readout of association and dissociation events. The relative ease with which sensing modalities can be incorporated and advantages in smaller size and photostability make de novo designed fluorescence-activating proteins attractive candidates for optical sensor engineering.

Improving the Efficiency of Ligand-Binding Protein Design with Molecular Dynamics Simulations.

Custom-designed ligand-binding proteins represent a promising class of macromolecules with exciting applications toward the design of new enzymes or the engineering of antibodies and small-molecule recruited proteins for therapeutic interventions. However, several challenges remain in designing a protein sequence such that the binding site organization results in high affinity interaction with a bound ligand. Here, we study the dynamics of explicitly solvated designed proteins through all-atom molecular dynamics (MD) simulations to gain insight into the causes that lead to the low affinity or instability of most of these designs, despite the prediction of their success by the computational design methodology. Simulations ranging from 500 to 1000 ns per replicate were conducted on 37 designed protein variants encompassing two distinct folds and a range of ligand affinities, resulting in more than 180 µs of combined sampling. The simulations provide retrospective insights into the properties affecting ligand affinity that can prove useful in guiding further steps of design optimization. Features indicate that entropic components are particularly important for affinity, which are not easily incorporated in the empirical models often used in design protocols. Additionally, we demonstrate that the application of machine learning approaches built upon the output from the simulations can help discriminate between successful and failed binders, such that MD could act as a screening step in protein design, resulting in a more efficient process.

Balancing Specificity and Promiscuity in Enzyme Evolution: Multidimensional Activity Transitions in the Alkaline Phosphatase Superfamily.

Highly proficient, promiscuous enzymes can be springboards for functional evolution, able to avoid loss of function during adaptation by their capacity to promote multiple reactions. We employ a systematic comparative study of structure, sequence, and substrate specificity to track the evolution of specificity and reactivity between promiscuous members of clades of the alkaline phosphatase (AP) superfamily. Construction of a phylogenetic tree of protein sequences maps out the likely transition zone between arylsulfatases (ASs) and phosphonate monoester hydrolases (PMHs). Kinetic analysis shows that all enzymes characterized have four chemically distinct phospho- and sulfoesterase activities, with rate accelerations ranging from 1011- to 1017-fold for their primary and 109- to 1012-fold for their promiscuous reactions, suggesting that catalytic promiscuity is widespread in the AP-superfamily. This functional characterization and crystallography reveal a novel class of ASs that is so similar in sequence to known PMHs that it had not been recognized as having diverged in function. Based on analysis of snapshots of catalytic promiscuity "in transition", we develop possible models that would allow functional evolution and determine scenarios for trade-off between multiple activities. For the new ASs, we observe largely invariant substrate specificity that would facilitate the transition from ASs to PMHs via trade-off-free molecular exaptation, that is, evolution without initial loss of primary activity and specificity toward the original substrate. This ability to bypass low activity generalists provides a molecular solution to avoid adaptive conflict.

De novo design of a fluorescence-activating ß-barrel.

The regular arrangements of ß-strands around a central axis in ß-barrels and of α-helices in coiled coils contrast with the irregular tertiary structures of most globular proteins, and have fascinated structural biologists since they were first discovered. Simple parametric models have been used to design a wide range of α-helical coiled-coil structures, but to date there has been no success with ß-barrels. Here we show that accurate de novo design of ß-barrels requires considerable symmetry-breaking to achieve continuous hydrogen-bond connectivity and eliminate backbone strain. We then build ensembles of ß-barrel backbone models with cavity shapes that match the fluorogenic compound DFHBI, and use a hierarchical grid-based search method to simultaneously optimize the rigid-body placement of DFHBI in these cavities and the identities of the surrounding amino acids to achieve high shape and chemical complementarity. The designs have high structural accuracy and bind and fluorescently activate DFHBI in vitro and in Escherichia coli, yeast and mammalian cells. This de novo design of small-molecule binding activity, using backbones custom-built to bind the ligand, should enable the design of increasingly sophisticated ligand-binding proteins, sensors and catalysts that are not limited by the backbone geometries available in known protein structures.

Escherichia coli D-malate dehydrogenase, a generalist enzyme active in the leucine biosynthesis pathway.

The enzymes of the ß-decarboxylating dehydrogenase superfamily catalyze the oxidative decarboxylation of D-malate-based substrates with various specificities. Here, we show that, in addition to its natural function affording bacterial growth on D-malate as a carbon source, the D-malate dehydrogenase of Escherichia coli (EcDmlA) naturally expressed from its chromosomal gene is capable of complementing leucine auxotrophy in a leuB(-) strain lacking the paralogous isopropylmalate dehydrogenase enzyme. To our knowledge, this is the first example of an enzyme that contributes with a physiologically relevant level of activity to two distinct pathways of the core metabolism while expressed from its chromosomal locus. EcDmlA features relatively high catalytic activity on at least three different substrates (L(+)-tartrate, D-malate, and 3-isopropylmalate). Because of these properties both in vivo and in vitro, EcDmlA may be defined as a generalist enzyme. Phylogenetic analysis highlights an ancient origin of DmlA, indicating that the enzyme has maintained its generalist character throughout evolution. We discuss the implication of these findings for protein evolution.