A computationally designed inhibitor of an Epstein-Barr viral Bcl-2 protein induces apoptosis in infected cells.

Because apoptosis of infected cells can limit virus production and spread, some viruses have co-opted prosurvival genes from the host. This includes the Epstein-Barr virus (EBV) gene BHRF1, a homolog of human Bcl-2 proteins that block apoptosis and are associated with cancer. Computational design and experimental optimization were used to generate a novel protein called BINDI that binds BHRF1 with picomolar affinity. BINDI recognizes the hydrophobic cleft of BHRF1 in a manner similar to other Bcl-2 protein interactions but makes many additional contacts to achieve exceptional affinity and specificity. BINDI induces apoptosis in EBV-infected cancer lines, and when delivered with an antibody-targeted intracellular delivery carrier, BINDI suppressed tumor growth and extended survival in a xenograft disease model of EBV-positive human lymphoma. High-specificity-designed proteins that selectively kill target cells may provide an advantage over the toxic compounds used in current generation antibody-drug conjugates.

A tethered ligand assay to probe SARS-CoV-2:ACE2 interactions.

SignificanceIn the dynamic environment of the airways, where SARS-CoV-2 infections are initiated by binding to human host receptor ACE2, mechanical stability of the viral attachment is a crucial fitness advantage. Using single-molecule force spectroscopy techniques, we mimic the effect of coughing and sneezing, thereby testing the force stability of SARS-CoV-2 RBD:ACE2 interaction under physiological conditions. Our results reveal a higher force stability of SARS-CoV-2 binding to ACE2 compared to SARS-CoV-1, causing a possible fitness advantage. Our assay is sensitive to blocking agents preventing RBD:ACE2 bond formation. It will thus provide a powerful approach to investigate the modes of action of neutralizing antibodies and other agents designed to block RBD binding to ACE2 that are currently developed as potential COVID-19 therapeutics.

Multiple mechanisms of self-association of chemokine receptors CXCR4 and CCR5 demonstrated by deep mutagenesis.

Chemokine receptors are members of the rhodopsin-like class A GPCRs whose signaling through G proteins drives the directional movement of cells in response to a chemokine gradient. Chemokine receptors CXCR4 and CCR5 have been extensively studied due to their roles in leukocyte development and inflammation and their status as coreceptors for HIV-1 infection, among other roles. Both receptors form dimers or oligomers of unclear function. While CXCR4 has been crystallized in a dimeric arrangement, available atomic resolution structures of CCR5 are monomeric. To investigate their dimerization interfaces, we used a bimolecular fluorescence complementation (BiFC)-based screen and deep mutational scanning to find mutations that change how the receptors self-associate, either via specific oligomer assembly or alternative mechanisms of clustering in close proximity. Many disruptive mutations promoted self-associations nonspecifically, suggesting they aggregated in the membrane. A mutationally intolerant region was found on CXCR4 that matched the crystallographic dimer interface, supporting this dimeric arrangement in living cells. A mutationally intolerant region was also observed on the surface of CCR5 by transmembrane helices 3 and 4. Mutations predicted from the scan to reduce BiFC were validated and were localized in the transmembrane domains as well as the C-terminal cytoplasmic tails where they reduced lipid microdomain localization. A mutation in the dimer interface of CXCR4 had increased binding to the ligand CXCL12 and yet diminished calcium signaling. There was no change in syncytia formation with cells expressing HIV-1 Env. The data highlight that multiple mechanisms are involved in self-association of chemokine receptor chains.

Sequence basis for selectivity of ephrin-B2 ligand for Eph receptors and pathogenic henipavirus G glycoproteins.

IMPORTANCE: Ephrin-B2 (EFNB2) is a ligand for six Eph receptors in humans and regulates multiple cell developmental and signaling processes. It also functions as the cell entry receptor for Nipah virus and Hendra virus, zoonotic viruses that can cause respiratory and/or neurological symptoms in humans with high mortality. Here, we investigate the sequence basis of EFNB2 specificity for binding the Nipah virus attachment G glycoprotein over Eph receptors. We then use this information to engineer EFNB2 as a soluble decoy receptor that specifically binds the attachment glycoproteins of the Nipah virus and other related henipaviruses to neutralize infection. These findings further mechanistic understanding of protein selectivity and may facilitate the development of diagnostics or therapeutics against henipavirus infection.

TAPBPR employs a ligand-independent docking mechanism to chaperone MR1 molecules.

Chaperones tapasin and transporter associated with antigen processing (TAP)-binding protein related (TAPBPR) associate with the major histocompatibility complex (MHC)-related protein 1 (MR1) to promote trafficking and cell surface expression. However, the binding mechanism and ligand dependency of MR1/chaperone interactions remain incompletely characterized. Here in vitro, biochemical and computational studies reveal that, unlike MHC-I, TAPBPR recognizes MR1 in a ligand-independent manner owing to the absence of major structural changes in the MR1 α2-1 helix between empty and ligand-loaded molecules. Structural characterization using paramagnetic nuclear magnetic resonance experiments combined with restrained molecular dynamics simulations reveals that TAPBPR engages conserved surfaces on MR1 to induce similar adaptations to those seen in MHC-I/TAPBPR co-crystal structures. Finally, nuclear magnetic resonance relaxation dispersion experiments using 19F-labeled diclofenac show that TAPBPR can affect the exchange kinetics of noncovalent metabolites with the MR1 groove, serving as a catalyst. Our results support a role of chaperones in stabilizing nascent MR1 molecules to enable loading of endogenous or exogenous cargo.

Engineered ACE2 decoy mitigates lung injury and death induced by SARS-CoV-2 variants.

Vaccine hesitancy and emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) escaping vaccine-induced immune responses highlight the urgency for new COVID-19 therapeutics. Engineered angiotensin-converting enzyme 2 (ACE2) proteins with augmented binding affinities for SARS-CoV-2 spike (S) protein may prove to be especially efficacious against multiple variants. Using molecular dynamics simulations and functional assays, we show that three amino acid substitutions in an engineered soluble ACE2 protein markedly augmented the affinity for the S protein of the SARS-CoV-2 WA-1/2020 isolate and multiple VOCs: B.1.1.7 (Alpha), B.1.351 (Beta), P.1 (Gamma) and B.1.617.2 (Delta). In humanized K18-hACE2 mice infected with the SARS-CoV-2 WA-1/2020 or P.1 variant, prophylactic and therapeutic injections of soluble ACE22.v2.4-IgG1 prevented lung vascular injury and edema formation, essential features of CoV-2-induced SARS, and above all improved survival. These studies demonstrate broad efficacy in vivo of an engineered ACE2 decoy against SARS-CoV-2 variants in mice and point to its therapeutic potential.

The substrate import mechanism of the human serotonin transporter.

The serotonin transporter (SERT) initiates the reuptake of extracellular serotonin in the synapse to terminate neurotransmission. The cryogenic electron microscopy structures of SERT bound to ibogaine and the physiological substrate serotonin resolved in different states have provided a glimpse of the functional conformations at atomistic resolution. However, the conformational dynamics and structural transitions to intermediate states are not fully understood. Furthermore, the molecular basis of how serotonin is recognized and transported remains unclear. In this study, we performed unbiased microsecond-long simulations of the human SERT to investigate the structural dynamics to various intermediate states and elucidated the complete substrate import pathway. Using Markov state models, we characterized a sequential order of conformational-driven ion-coupled substrate binding and transport events and calculated the free energy barriers of conformation transitions associated with the import mechanism. We find that the transition from the occluded to inward-facing state is the rate-limiting step for substrate import and that the substrate decreases the free energy barriers to achieve the inward-facing state. Our study provides insights on the molecular basis of dynamics-driven ion-substrate recognition and transport of SERT that can serve as a model for other closely related neurotransmitter transporters.

Engineered receptors for human cytomegalovirus that are orthogonal to normal human biology.

A trimeric glycoprotein complex on the surface of human cytomegalovirus (HCMV) binds to platelet-derived growth factor (PDGF) receptor α (PDGFRα) to mediate host cell recognition and fusion of the viral and cellular membranes. Soluble PDGFRα potently neutralizes HCMV in tissue culture, and its potential use as an antiviral therapeutic has the benefit that any escape mutants will likely be attenuated. However, PDGFRα binds multiple PDGF ligands in the human body as part of developmental programs in embryogenesis and continuing through adulthood. Any therapies with soluble receptor therefore come with serious efficacy and safety concerns, especially for the treatment of congenital HCMV. Soluble virus receptors that are orthogonal to human biology might resolve these concerns. This engineering problem is solved by deep mutational scanning on the D2-D3 domains of PDGFRα to identify variants that maintain interactions with the HCMV glycoprotein trimer in the presence of competing PDGF ligands. Competition by PDGFs is conformation-dependent, whereas HCMV trimer binding is independent of proper D2-D3 conformation, and many mutations at the receptor-PDGF interface are suitable for functionally separating trimer from PDGF interactions. Purified soluble PDGFRα carrying a targeted mutation succeeded in displaying wild type affinity for HCMV trimer with a simultaneous loss of PDGF binding, and neutralizes trimer-only and trimer/pentamer-expressing HCMV strains infecting fibroblasts or epithelial cells. Overall, this work makes important progress in the realization of soluble HCMV receptors for clinical application.

Molecular determinants of chaperone interactions on MHC-I for folding and antigen repertoire selection.

The interplay between a highly polymorphic set of MHC-I alleles and molecular chaperones shapes the repertoire of peptide antigens displayed on the cell surface for T cell surveillance. Here, we demonstrate that the molecular chaperone TAP-binding protein related (TAPBPR) associates with a broad range of partially folded MHC-I species inside the cell. Bimolecular fluorescence complementation and deep mutational scanning reveal that TAPBPR recognition is polarized toward the α2 domain of the peptide-binding groove, and depends on the formation of a conserved MHC-I disulfide epitope in the α2 domain. Conversely, thermodynamic measurements of TAPBPR binding for a representative set of properly conformed, peptide-loaded molecules suggest a narrower MHC-I specificity range. Using solution NMR, we find that the extent of dynamics at "hotspot" surfaces confers TAPBPR recognition of a sparsely populated MHC-I state attained through a global conformational change. Consistently, restriction of MHC-I groove plasticity through the introduction of a disulfide bond between the α1/α2 helices abrogates TAPBPR binding, both in solution and on a cellular membrane, while intracellular binding is tolerant of many destabilizing MHC-I substitutions. Our data support parallel TAPBPR functions of 1) chaperoning unstable MHC-I molecules with broad allele-specificity at early stages of their folding process, and 2) editing the peptide cargo of properly conformed MHC-I molecules en route to the surface, which demonstrates a narrower specificity. Our results suggest that TAPBPR exploits localized structural adaptations, both near and distant to the peptide-binding groove, to selectively recognize discrete conformational states sampled by MHC-I alleles, toward editing the repertoire of displayed antigens.

ACE2-based decoy receptors for SARS coronavirus 2.

SARS coronavirus 2 is neutralized by proteins that block receptor-binding sites on spikes that project from the viral envelope. In particular, substantial research investment has advanced monoclonal antibody therapies to the clinic where they have shown partial efficacy in reducing viral burden and hospitalization. An alternative is to use the host entry receptor, angiotensin-converting enzyme 2 (ACE2), as a soluble decoy that broadly blocks SARS-associated coronaviruses with limited potential for viral escape. Here, we summarize efforts to engineer higher affinity variants of soluble ACE2 that rival the potency of affinity-matured antibodies. Strategies have also been used to increase the valency of ACE2 decoys for avid spike interactions and to improve pharmacokinetics via IgG fusions. Finally, the intrinsic catalytic activity of ACE2 for the turnover of the vasoconstrictor angiotensin II may directly address COVID-19 symptoms and protect against lung and cardiovascular injury, conferring dual mechanisms of action unachievable by monoclonal antibodies. Soluble ACE2 derivatives therefore have the potential to be next generation therapeutics for addressing the immediate needs of the current pandemic and possible future outbreaks.

Computationally Designed ACE2 Decoy Receptor Binds SARS-CoV-2 Spike (S) Protein with Tight Nanomolar Affinity.

Even with the availability of vaccines, therapeutic options for COVID-19 still remain highly desirable, especially in hospitalized patients with moderate or severe disease. Soluble ACE2 (sACE2) is a promising therapeutic candidate that neutralizes SARS CoV-2 infection by acting as a decoy. Using computational mutagenesis, we designed a number of sACE2 derivatives carrying three to four mutations. The top-predicted sACE2 decoy based on the in silico mutagenesis scan was subjected to molecular dynamics and free-energy calculations for further validation. After illuminating the mechanism of increased binding for our designed sACE2 derivative, the design was verified experimentally by flow cytometry and BLI-binding experiments. The computationally designed sACE2 decoy (ACE2-FFWF) bound the receptor-binding domain of SARS-CoV-2 tightly with low nanomolar affinity and ninefold affinity enhancement over the wild type. Furthermore, cell surface expression was slightly greater than wild-type ACE2, suggesting that the design is well-folded and stable. Having an arsenal of high-affinity sACE2 derivatives will help to buffer against the emergence of SARS CoV-2 variants. Here, we show that computational methods have become sufficiently accurate for the design of therapeutics for current and future viral pandemics.

Structural architecture of a dimeric class C GPCR based on co-trafficking of sweet taste receptor subunits.

Class C G protein-coupled receptors (GPCRs) are obligatory dimers that are particularly important for neuronal responses to endogenous and environmental stimuli. Ligand recognition through large extracellular domains leads to the reorganization of transmembrane regions to activate G protein signaling. Although structures of individual domains are known, the complete architecture of a class C GPCR and the mechanism of interdomain coupling during receptor activation are unclear. By screening a mutagenesis library of the human class C sweet taste receptor subunit T1R2, we enhanced surface expression and identified a dibasic intracellular retention motif that modulates surface expression and co-trafficking with its heterodimeric partner T1R3. Using a highly expressed T1R2 variant, dimerization sites along the entire subunit within all the structural domains were identified by a comprehensive mutational scan for co-trafficking with T1R3 in human cells. The data further reveal that the C terminus of the extracellular cysteine-rich domain needs to be properly folded for T1R3 dimerization and co-trafficking, but not for surface expression of T1R2 alone. These results guided the modeling of the T1R2-T1R3 dimer in living cells, which predicts a twisted arrangement of domains around the central axis, and a continuous folded structure between transmembrane domain loops and the cysteine-rich domains. These insights have implications for how conformational changes between domains are coupled within class C GPCRs.

Conformational Engineering of HIV-1 Env Based on Mutational Tolerance in the CD4 and PG16 Bound States.

HIV-1 infection is initiated by viral Env engaging the host receptor CD4, triggering Env to transition from a "closed" to "open" conformation during the early events of virus-cell membrane fusion. To understand how Env sequence accommodates this conformational change, mutational landscapes decoupled from virus replication were determined for Env from BaL (clade B) and DU422 (clade C) isolates interacting with CD4 or antibody PG16 that preferentially recognizes closed trimers. Sequence features uniquely important to each bound state were identified, including glycosylation and binding sites. Notably, the Env apical domain and trimerization interface are under selective pressure for PG16 binding. Based on this key observation, mutations were found that increase presentation of quaternary epitopes associated with properly conformed trimers when Env is expressed at the plasma membrane. Many mutations reduce electrostatic repulsion at the Env apex and increase PG16 recognition of Env sequences from clades A and B. Other mutations increase hydrophobic packing at the gp120 inner-outer domain interface and were broadly applicable for engineering Env from diverse strains spanning tiers 1, 2, and 3 across clades A, B, C, and BC recombinants. Core mutations predicted to introduce steric strain in the open state show markedly reduced CD4 interactions. Finally, we demonstrate how our methodology can be adapted to interrogate interactions between membrane-associated Env and the matrix domain of Gag. These findings and methods may assist vaccine design.IMPORTANCE HIV-1 Env is dynamic and undergoes large conformational changes that drive fusion of virus and host cell membranes. Three Env proteins in a trimer contact each other at their apical tips to form a closed conformation that presents epitopes recognized by broadly neutralizing antibodies. The apical tips separate, among other changes, to form an open conformation that binds tightly to host receptors. Understanding how Env sequence facilitates these structural changes can inform the biophysical mechanism and aid immunogen design. Using deep mutational scans decoupled from virus replication, we report mutational landscapes for Env from two strains interacting with conformation-dependent binding proteins. Residues in the Env trimer interface and apical domains are preferentially conserved in the closed conformation, and conformational diversity is facilitated by electrostatic repulsion and an underpacked core between domains. Specific mutations are described that enhance presentation of the trimeric closed conformation across diverse HIV-1 strains.

Deep mutagenesis in the study of COVID-19: a technical overview for the proteomics community.

INTRODUCTION: The spike (S) of SARS coronavirus 2 (SARS-CoV-2) engages angiotensin-converting enzyme 2 (ACE2) on a host cell to trigger viral-cell membrane fusion and infection. The extracellular region of ACE2 can be administered as a soluble decoy to compete for binding sites on the receptor-binding domain (RBD) of S, but it has only moderate affinity and efficacy. The RBD, which is targeted by neutralizing antibodies, may also change and adapt through mutation as SARS-CoV-2 becomes endemic, posing challenges for therapeutic and vaccine development. AREAS COVERED: Deep mutagenesis is a Big Data approach to characterizing sequence variants. A deep mutational scan of ACE2 expressed on human cells identified mutations that increase S affinity and guided the engineering of a potent and broad soluble receptor decoy. A deep mutational scan of the RBD displayed on the surface of yeast has revealed residues tolerant of mutational changes that may act as a source for drug resistance and antigenic drift. EXPERT OPINION: Deep mutagenesis requires a selection of diverse sequence variants; an in vitro evolution experiment that is tracked with next-generation sequencing. The choice of expression system, diversity of the variant library and selection strategy have important consequences for data quality and interpretation.

Mapping Interaction Sites on Human Chemokine Receptors by Deep Mutational Scanning.

Chemokine receptors CXCR4 and CCR5 regulate WBC trafficking and are engaged by the HIV-1 envelope glycoprotein gp120 during infection. We combine a selection of human CXCR4 and CCR5 libraries comprising nearly all of ∼7000 single amino acid substitutions with deep sequencing to define sequence-activity landscapes for surface expression and ligand interactions. After consideration of sequence constraints for surface expression, known interaction sites with HIV-1-blocking Abs were appropriately identified as conserved residues following library sorting for Ab binding, validating the use of deep mutational scanning to map functional interaction sites in G protein-coupled receptors. Chemokine CXCL12 was found to interact with residues extending asymmetrically into the CXCR4 ligand-binding cavity, similar to the binding surface of CXCR4 recognized by an antagonistic viral chemokine previously observed crystallographically. CXCR4 mutations distal from the chemokine binding site were identified that enhance chemokine recognition. This included disruptive mutations in the G protein-coupling site that diminished calcium mobilization, as well as conservative mutations to a membrane-exposed site (CXCR4 residues H792.45 and W1614.50) that increased ligand binding without loss of signaling. Compared with CXCR4-CXCL12 interactions, CCR5 residues conserved for gp120 (HIV-1 BaL strain) interactions map to a more expansive surface, mimicking how the cognate chemokine CCL5 makes contacts across the entire CCR5 binding cavity. Acidic substitutions in the CCR5 N terminus and extracellular loops enhanced gp120 binding. This study demonstrates how comprehensive mutational scanning can define functional interaction sites on receptors, and novel mutations that enhance receptor activities can be found simultaneously.

Structural and energetic basis of folded-protein transport by the FimD usher.

Type 1 pili, produced by uropathogenic Escherichia coli, are multisubunit fibres crucial in recognition of and adhesion to host tissues. During pilus biogenesis, subunits are recruited to an outer membrane assembly platform, the FimD usher, which catalyses their polymerization and mediates pilus secretion. The recent determination of the crystal structure of an initiation complex provided insight into the initiation step of pilus biogenesis resulting in pore activation, but very little is known about the elongation steps that follow. Here, to address this question, we determine the structure of an elongation complex in which the tip complex assembly composed of FimC, FimF, FimG and FimH passes through FimD. This structure demonstrates the conformational changes required to prevent backsliding of the nascent pilus through the FimD pore and also reveals unexpected properties of the usher pore. We show that the circular binding interface between the pore lumen and the folded substrate participates in transport by defining a low-energy pathway along which the nascent pilus polymer is guided during secretion.

Reduced axonal surface expression and phosphoinositide sensitivity in Kv7 channels disrupts their function to inhibit neuronal excitability in Kcnq2 epileptic encephalopathy.

Neuronal Kv7/KCNQ channels are voltage-gated potassium channels composed of Kv7.2/KCNQ2 and Kv7.3/KCNQ3 subunits. Enriched at the axonal membrane, they potently suppress neuronal excitability. De novo and inherited dominant mutations in Kv7.2 cause early onset epileptic encephalopathy characterized by drug resistant seizures and profound psychomotor delay. However, their precise pathogenic mechanisms remain elusive. Here, we investigated selected epileptic encephalopathy causing mutations in calmodulin (CaM)-binding helices A and B of Kv7.2. We discovered that R333W, K526N, and R532W mutations located peripheral to CaM contact sites decreased axonal surface expression of heteromeric channels although only R333W mutation reduced CaM binding to Kv7.2. These mutations also altered gating modulation by phosphatidylinositol 4,5-bisphosphate (PIP2), revealing novel PIP2 binding residues. While these mutations disrupted Kv7 function to suppress excitability, hyperexcitability was observed in neurons expressing Kv7.2-R532W that displayed severe impairment in voltage-dependent activation. The M518 V mutation at the CaM contact site in helix B caused most defects in Kv7 channels by severely reducing their CaM binding, K+ currents, and axonal surface expression. Interestingly, the M518 V mutation induced ubiquitination and accelerated proteasome-dependent degradation of Kv7.2, whereas the presence of Kv7.3 blocked this degradation. Furthermore, expression of Kv7.2-M518V increased neuronal death. Together, our results demonstrate that epileptic encephalopathy mutations in helices A and B of Kv7.2 cause abnormal Kv7 expression and function by disrupting Kv7.2 binding to CaM and/or modulation by PIP2. We propose that such multiple Kv7 channel defects could exert more severe impacts on neuronal excitability and health, and thus serve as pathogenic mechanisms underlying Kcnq2 epileptic encephalopathy.

Deep Mutational Scans as a Guide to Engineering High Affinity T Cell Receptor Interactions with Peptide-bound Major Histocompatibility Complex.

Proteins are often engineered to have higher affinity for their ligands to achieve therapeutic benefit. For example, many studies have used phage or yeast display libraries of mutants within complementarity-determining regions to affinity mature antibodies and T cell receptors (TCRs). However, these approaches do not allow rapid assessment or evolution across the entire interface. By combining directed evolution with deep sequencing, it is now possible to generate sequence fitness landscapes that survey the impact of every amino acid substitution across the entire protein-protein interface. Here we used the results of deep mutational scans of a TCR-peptide-MHC interaction to guide mutational strategies. The approach yielded stable TCRs with affinity increases of >200-fold. The substitutions with the greatest enrichments based on the deep sequencing were validated to have higher affinity and could be combined to yield additional improvements. We also conducted in silico binding analyses for every substitution to compare them with the fitness landscape. Computational modeling did not effectively predict the impacts of mutations distal to the interface and did not account for yeast display results that depended on combinations of affinity and protein stability. However, computation accurately predicted affinity changes for mutations within or near the interface, highlighting the complementary strengths of computational modeling and yeast surface display coupled with deep mutational scanning for engineering high affinity TCRs.

Intracellular delivery system for antibody-Peptide drug conjugates.

Antibodies armed with biologic drugs could greatly expand the therapeutic potential of antibody-drug conjugates for cancer therapy, broadening their application to disease targets currently limited by intracellular delivery barriers. Additional selectivity and new therapeutic approaches could be realized with intracellular protein drugs that more specifically target dysregulated pathways in hematologic cancers and other malignancies. A multifunctional polymeric delivery system for enhanced cytosolic delivery of protein drugs has been developed that incorporates endosomal-releasing activity, antibody targeting, and a biocompatible long-chain ethylene glycol component for optimized safety, pharmacokinetics, and tumor biodistribution. The pH-responsive polymeric micelle carrier, with an internalizing anti-CD22 monoclonal targeting antibody, effectively delivered a proapoptotic Bcl-2 interacting mediator (BIM) peptide drug that suppressed tumor growth for the duration of treatment and prolonged survival in a xenograft mouse model of human B-cell lymphoma. Antitumor drug activity was correlated with a mechanistic induction of the Bcl-2 pathway biomarker cleaved caspase-3 and a marked decrease in the Ki-67 proliferation biomarker. Broadening the intracellular target space by more effective delivery of protein/peptide drugs could expand the repertoire of antibody-drug conjugates to currently undruggable disease-specific targets and permit tailored drug strategies to stratified subpopulations and personalized medicines.

Deep mutational scanning reveals sequence to function constraints for SWEET family transporters.

Protein science is entering a transformative phase enabled by deep mutational scans that provide an unbiased view of the residue level interactions that mediate function. However, it has yet to be extensively used to characterize the mutational and evolutionary landscapes of plant proteins. Here, we apply the method to explore sequence-function relationships within the sugar transporter AtSWEET13. DMS results describe how mutational interrogation throughout different regions of the protein affects AtSWEET13 abundance and transport function. Our results identify novel transport-enhancing mutations that are validated using the FRET sensor assays. Extending DMS results to phylogenetic analyses reveal the role of transmembrane helix 4 (TM4) which makes the SWEET family transporters distinct from prokaryotic SemiSWEETs. We show that transmembrane helix 4 is intolerant to motif swapping with other clade-specific SWEET TM4 compositions, despite accommodating single point-mutations towards aromatic and charged polar amino acids. We further show that the transfer learning approaches based on physics and ML based In silico variant prediction tools have limited utility for engineering plant proteins as they were unable to reproduce our experimental results. We conclude that DMS can produce datasets which, when combined with the right predictive computational frameworks, can direct plant engineering efforts through derivative phenotype selection and evolutionary insights.