Deep mutational scanning of hepatitis B virus reveals a mechanism for cis-preferential reverse transcription.

Hepatitis B virus (HBV) is a small double-stranded DNA virus that chronically infects 296 million people. Over half of its compact genome encodes proteins in two overlapping reading frames, and during evolution, multiple selective pressures can act on shared nucleotides. This study combines an RNA-based HBV cell culture system with deep mutational scanning (DMS) to uncouple cis- and trans-acting sequence requirements in the HBV genome. The results support a leaky ribosome scanning model for polymerase translation, provide a fitness map of the HBV polymerase at single-nucleotide resolution, and identify conserved prolines adjacent to the HBV polymerase termination codon that stall ribosomes. Further experiments indicated that stalled ribosomes tether the nascent polymerase to its template RNA, ensuring cis-preferential RNA packaging and reverse transcription of the HBV genome.

Structured States of Disordered Proteins from Genomic Sequences.

Protein flexibility ranges from simple hinge movements to functional disorder. Around half of all human proteins contain apparently disordered regions with little 3D or functional information, and many of these proteins are associated with disease. Building on the evolutionary couplings approach previously successful in predicting 3D states of ordered proteins and RNA, we developed a method to predict the potential for ordered states for all apparently disordered proteins with sufficiently rich evolutionary information. The approach is highly accurate (79%) for residue interactions as tested in more than 60 known disordered regions captured in a bound or specific condition. Assessing the potential for structure of more than 1,000 apparently disordered regions of human proteins reveals a continuum of structural order with at least 50% with clear propensity for three- or two-dimensional states. Co-evolutionary constraints reveal hitherto unseen structures of functional importance in apparently disordered proteins.

3D RNA and Functional Interactions from Evolutionary Couplings.

Non-coding RNAs are ubiquitous, but the discovery of new RNA gene sequences far outpaces the research on the structure and functional interactions of these RNA gene sequences. We mine the evolutionary sequence record to derive precise information about the function and structure of RNAs and RNA-protein complexes. As in protein structure prediction, we use maximum entropy global probability models of sequence co-variation to infer evolutionarily constrained nucleotide-nucleotide interactions within RNA molecules and nucleotide-amino acid interactions in RNA-protein complexes. The predicted contacts allow all-atom blinded 3D structure prediction at good accuracy for several known RNA structures and RNA-protein complexes. For unknown structures, we predict contacts in 160 non-coding RNA families. Beyond 3D structure prediction, evolutionary couplings help identify important functional interactions-e.g., at switch points in riboswitches and at a complex nucleation site in HIV. Aided by increasing sequence accumulation, evolutionary coupling analysis can accelerate the discovery of functional interactions and 3D structures involving RNA.

The SpoVA membrane complex is required for dipicolinic acid import during sporulation and export during germination.

In response to starvation, endospore-forming bacteria differentiate into stress-resistant spores that can remain dormant for years yet rapidly germinate and resume growth in response to nutrients. The small molecule dipicolinic acid (DPA) plays a central role in both the stress resistance of the dormant spore and its exit from dormancy during germination. The spoVA locus is required for DPA import during sporulation and has been implicated in its export during germination, but the molecular bases are unclear. Here, we define the minimal set of proteins encoded in the Bacillus subtilis spoVA operon required for DPA import and demonstrate that these proteins form a membrane complex. Structural modeling of these components combined with mutagenesis and in vivo analysis reveal that the C and Eb subunits form a membrane channel, while the D subunit functions as a cytoplasmic plug. We show that point mutations that impair the interactions between D and the C-Eb membrane complex reduce the efficiency of DPA import during sporulation and reciprocally accelerate DPA release during germination. Our data support a model in which DPA transport into spores involves cycles of unplugging and then replugging the C-Eb membrane channel, while nutrient detection during germination triggers DPA release by unplugging it.

Coordination of bacterial cell wall and outer membrane biosynthesis.

Gram-negative bacteria surround their cytoplasmic membrane with a peptidoglycan (PG) cell wall and an outer membrane (OM) with an outer leaflet composed of lipopolysaccharide (LPS)1. This complex envelope presents a formidable barrier to drug entry and is a major determinant of the intrinsic antibiotic resistance of these organisms2. The biogenesis pathways that build the surface are also targets of many of our most effective antibacterial therapies3. Understanding the molecular mechanisms underlying the assembly of the Gram-negative envelope therefore promises to aid the development of new treatments effective against the growing problem of drug-resistant infections. Although the individual pathways for PG and OM synthesis and assembly are well characterized, almost nothing is known about how the biogenesis of these essential surface layers is coordinated. Here we report the discovery of a regulatory interaction between the committed enzymes for the PG and LPS synthesis pathways in the Gram-negative pathogen Pseudomonas aeruginosa. We show that the PG synthesis enzyme MurA interacts directly and specifically with the LPS synthesis enzyme LpxC. Moreover, MurA was shown to stimulate LpxC activity in cells and in a purified system. Our results support a model in which the assembly of the PG and OM layers in many proteobacterial species is coordinated by linking the activities of the committed enzymes in their respective synthesis pathways.

Learning from prepandemic data to forecast viral escape.

Effective pandemic preparedness relies on anticipating viral mutations that are able to evade host immune responses to facilitate vaccine and therapeutic design. However, current strategies for viral evolution prediction are not available early in a pandemic-experimental approaches require host polyclonal antibodies to test against1-16, and existing computational methods draw heavily from current strain prevalence to make reliable predictions of variants of concern17-19. To address this, we developed EVEscape, a generalizable modular framework that combines fitness predictions from a deep learning model of historical sequences with biophysical and structural information. EVEscape quantifies the viral escape potential of mutations at scale and has the advantage of being applicable before surveillance sequencing, experimental scans or three-dimensional structures of antibody complexes are available. We demonstrate that EVEscape, trained on sequences available before 2020, is as accurate as high-throughput experimental scans at anticipating pandemic variation for SARS-CoV-2 and is generalizable to other viruses including influenza, HIV and understudied viruses with pandemic potential such as Lassa and Nipah. We provide continually revised escape scores for all current strains of SARS-CoV-2 and predict probable further mutations to forecast emerging strains as a tool for continuing vaccine development ( evescape.org ).

scPerturb: harmonized single-cell perturbation data.

Analysis across a growing number of single-cell perturbation datasets is hampered by poor data interoperability. To facilitate development and benchmarking of computational methods, we collect a set of 44 publicly available single-cell perturbation-response datasets with molecular readouts, including transcriptomics, proteomics and epigenomics. We apply uniform quality control pipelines and harmonize feature annotations. The resulting information resource, scPerturb, enables development and testing of computational methods, and facilitates comparison and integration across datasets. We describe energy statistics (E-statistics) for quantification of perturbation effects and significance testing, and demonstrate E-distance as a general distance measure between sets of single-cell expression profiles. We illustrate the application of E-statistics for quantifying similarity and efficacy of perturbations. The perturbation-response datasets and E-statistics computation software are publicly available at scperturb.org. This work provides an information resource for researchers working with single-cell perturbation data and recommendations for experimental design, including optimal cell counts and read depth.

Three-dimensional structures of membrane proteins from genomic sequencing.

We show that amino acid covariation in proteins, extracted from the evolutionary sequence record, can be used to fold transmembrane proteins. We use this technique to predict previously unknown 3D structures for 11 transmembrane proteins (with up to 14 helices) from their sequences alone. The prediction method (EVfold_membrane) applies a maximum entropy approach to infer evolutionary covariation in pairs of sequence positions within a protein family and then generates all-atom models with the derived pairwise distance constraints. We benchmark the approach with blinded de novo computation of known transmembrane protein structures from 23 families, demonstrating unprecedented accuracy of the method for large transmembrane proteins. We show how the method can predict oligomerization, functional sites, and conformational changes in transmembrane proteins. With the rapid rise in large-scale sequencing, more accurate and more comprehensive information on evolutionary constraints can be decoded from genetic variation, greatly expanding the repertoire of transmembrane proteins amenable to modeling by this method.

Disease variant prediction with deep generative models of evolutionary data.

Quantifying the pathogenicity of protein variants in human disease-related genes would have a marked effect on clinical decisions, yet the overwhelming majority (over 98%) of these variants still have unknown consequences1-3. In principle, computational methods could support the large-scale interpretation of genetic variants. However, state-of-the-art methods4-10 have relied on training machine learning models on known disease labels. As these labels are sparse, biased and of variable quality, the resulting models have been considered insufficiently reliable11. Here we propose an approach that leverages deep generative models to predict variant pathogenicity without relying on labels. By modelling the distribution of sequence variation across organisms, we implicitly capture constraints on the protein sequences that maintain fitness. Our model EVE (evolutionary model of variant effect) not only outperforms computational approaches that rely on labelled data but also performs on par with, if not better than, predictions from high-throughput experiments, which are increasingly used as evidence for variant classification12-16. We predict the pathogenicity of more than 36 million variants across 3,219 disease genes and provide evidence for the classification of more than 256,000 variants of unknown significance. Our work suggests that models of evolutionary information can provide valuable independent evidence for variant interpretation that will be widely useful in research and clinical settings.

A Systematic p53 Mutation Library Links Differential Functional Impact to Cancer Mutation Pattern and Evolutionary Conservation.

The TP53 gene is frequently mutated in human cancer. Research has focused predominantly on six major "hotspot" codons, which account for only ∼30% of cancer-associated p53 mutations. To comprehensively characterize the consequences of the p53 mutation spectrum, we created a synthetically designed library and measured the functional impact of ∼10,000 DNA-binding domain (DBD) p53 variants in human cells in culture and in vivo. Our results highlight the differential outcome of distinct p53 mutations in human patients and elucidate the selective pressure driving p53 conservation throughout evolution. Furthermore, while loss of anti-proliferative functionality largely correlates with the occurrence of cancer-associated p53 mutations, we observe that selective gain-of-function may further favor particular mutants in vivo. Finally, when combined with additional acquired p53 mutations, seemingly neutral TP53 SNPs may modulate phenotypic outcome and, presumably, tumor progression.

The relaxin receptor RXFP1 signals through a mechanism of autoinhibition.

The relaxin family peptide receptor 1 (RXFP1) is the receptor for relaxin-2, an important regulator of reproductive and cardiovascular physiology. RXFP1 is a multi-domain G protein-coupled receptor (GPCR) with an ectodomain consisting of a low-density lipoprotein receptor class A (LDLa) module and leucine-rich repeats. The mechanism of RXFP1 signal transduction is clearly distinct from that of other GPCRs, but remains very poorly understood. In the present study, we determine the cryo-electron microscopy structure of active-state human RXFP1, bound to a single-chain version of the endogenous agonist relaxin-2 and the heterotrimeric Gs protein. Evolutionary coupling analysis and structure-guided functional experiments reveal that RXFP1 signals through a mechanism of autoinhibition. Our results explain how an unusual GPCR family functions, providing a path to rational drug development targeting the relaxin receptors.

Structure of the peptidoglycan polymerase RodA resolved by evolutionary coupling analysis.

The shape, elongation, division and sporulation (SEDS) proteins are a large family of ubiquitous and essential transmembrane enzymes with critical roles in bacterial cell wall biology. The exact function of SEDS proteins was for a long time poorly understood, but recent work has revealed that the prototypical SEDS family member RodA is a peptidoglycan polymerase-a role previously attributed exclusively to members of the penicillin-binding protein family. This discovery has made RodA and other SEDS proteins promising targets for the development of next-generation antibiotics. However, little is known regarding the molecular basis of SEDS activity, and no structural data are available for RodA or any homologue thereof. Here we report the crystal structure of Thermus thermophilus RodA at a resolution of 2.9 Å, determined using evolutionary covariance-based fold prediction to enable molecular replacement. The structure reveals a ten-pass transmembrane fold with large extracellular loops, one of which is partially disordered. The protein contains a highly conserved cavity in the transmembrane domain, reminiscent of ligand-binding sites in transmembrane receptors. Mutagenesis experiments in Bacillus subtilis and Escherichia coli show that perturbation of this cavity abolishes RodA function both in vitro and in vivo, indicating that this cavity is catalytically essential. These results provide a framework for understanding bacterial cell wall synthesis and SEDS protein function.

An ANXA11 P93S variant dysregulates TDP­43 and causes corticobasal syndrome.

INTRODUCTION: Variants of uncertain significance (VUS) surged with affordable genetic testing, posing challenges for determining pathogenicity. We examine the pathogenicity of a novel VUS P93S in Annexin A11 (ANXA11) - an amyotrophic lateral sclerosis/frontotemporal dementia-associated gene - in a corticobasal syndrome kindred. Established ANXA11 mutations cause ANXA11 aggregation, altered lysosomal-RNA granule co-trafficking, and transactive response DNA binding protein of 43 kDa (TDP-43) mis-localization. METHODS: We described the clinical presentation and explored the phenotypic diversity of ANXA11 variants. P93S's effect on ANXA11 function and TDP-43 biology was characterized in induced pluripotent stem cell-derived neurons alongside multiomic neuronal and microglial profiling. RESULTS: ANXA11 mutations were linked to corticobasal syndrome cases. P93S led to decreased lysosome colocalization, neuritic RNA, and nuclear TDP-43 with cryptic exon expression. Multiomic microglial signatures implicated immune dysregulation and interferon signaling pathways. DISCUSSION: This study establishes ANXA11 P93S pathogenicity, broadens the phenotypic spectrum of ANXA11 mutations, underscores neuronal and microglial dysfunction in ANXA11 pathophysiology, and demonstrates the potential of cellular models to determine variant pathogenicity. HIGHLIGHTS: ANXA11 P93S is a pathogenic variant. Corticobasal syndrome is part of the ANXA11 phenotypic spectrum. Hybridization chain reaction fluorescence in situ hybridization (HCR FISH) is a new tool for the detection of cryptic exons due to TDP-43-related loss of splicing regulation. Microglial ANXA11 and related immune pathways are important drivers of disease. Cellular models are powerful tools for adjudicating variants of uncertain significance.

Rapid generation of potent antibodies by autonomous hypermutation in yeast.

The predominant approach for antibody generation remains animal immunization, which can yield exceptionally selective and potent antibody clones owing to the powerful evolutionary process of somatic hypermutation. However, animal immunization is inherently slow, not always accessible and poorly compatible with many antigens. Here, we describe 'autonomous hypermutation yeast surface display' (AHEAD), a synthetic recombinant antibody generation technology that imitates somatic hypermutation inside engineered yeast. By encoding antibody fragments on an error-prone orthogonal DNA replication system, surface-displayed antibody repertoires continuously mutate through simple cycles of yeast culturing and enrichment for antigen binding to produce high-affinity clones in as little as two weeks. We applied AHEAD to generate potent nanobodies against the SARS-CoV-2 S glycoprotein, a G-protein-coupled receptor and other targets, offering a template for streamlined antibody generation at large.

Interaction specificity of clustered protocadherins inferred from sequence covariation and structural analysis.

Clustered protocadherins, a large family of paralogous proteins that play important roles in neuronal development, provide an important case study of interaction specificity in a large eukaryotic protein family. A mammalian genome has more than 50 clustered protocadherin isoforms, which have remarkable homophilic specificity for interactions between cellular surfaces. A large antiparallel dimer interface formed by the first 4 extracellular cadherin (EC) domains controls this interaction. To understand how specificity is achieved between the numerous paralogs, we used a combination of structural and computational approaches. Molecular dynamics simulations revealed that individual EC interactions are weak and undergo binding and unbinding events, but together they form a stable complex through polyvalency. Strongly evolutionarily coupled residue pairs interacted more frequently in our simulations, suggesting that sequence coevolution can inform the frequency of interaction and biochemical nature of a residue interaction. With these simulations and sequence coevolution, we generated a statistical model of interaction energy for the clustered protocadherin family that measures the contributions of all amino acid pairs at the interface. Our interaction energy model assesses specificity for all possible pairs of isoforms, recapitulating known pairings and predicting the effects of experimental changes in isoform specificity that are consistent with literature results. Our results show that sequence coevolution can be used to understand specificity determinants in a protein family and prioritize interface amino acid substitutions to reprogram specific protein-protein interactions.

Deep generative models of genetic variation capture the effects of mutations.

The functions of proteins and RNAs are defined by the collective interactions of many residues, and yet most statistical models of biological sequences consider sites nearly independently. Recent approaches have demonstrated benefits of including interactions to capture pairwise covariation, but leave higher-order dependencies out of reach. Here we show how it is possible to capture higher-order, context-dependent constraints in biological sequences via latent variable models with nonlinear dependencies. We found that DeepSequence ( https://github.com/debbiemarkslab/DeepSequence ), a probabilistic model for sequence families, predicted the effects of mutations across a variety of deep mutational scanning experiments substantially better than existing methods based on the same evolutionary data. The model, learned in an unsupervised manner solely on the basis of sequence information, is grounded with biologically motivated priors, reveals the latent organization of sequence families, and can be used to explore new parts of sequence space.

Evidence that regulation of intramembrane proteolysis is mediated by substrate gating during sporulation in Bacillus subtilis.

During the morphological process of sporulation in Bacillus subtilis two adjacent daughter cells (called the mother cell and forespore) follow different programs of gene expression that are linked to each other by signal transduction pathways. At a late stage in development, a signaling pathway emanating from the forespore triggers the proteolytic activation of the mother cell transcription factor σK. Cleavage of pro-σK to its mature and active form is catalyzed by the intramembrane cleaving metalloprotease SpoIVFB (B), a Site-2 Protease (S2P) family member. B is held inactive by two mother-cell membrane proteins SpoIVFA (A) and BofA. Activation of pro-σK processing requires a site-1 signaling protease SpoIVB (IVB) that is secreted from the forespore into the space between the two cells. IVB cleaves the extracellular domain of A but how this cleavage activates intramembrane proteolysis has remained unclear. Structural studies of the Methanocaldococcus jannaschii S2P homolog identified closed (substrate-occluded) and open (substrate-accessible) conformations of the protease, but the biological relevance of these conformations has not been established. Here, using co-immunoprecipitation and fluorescence microscopy, we show that stable association between the membrane-embedded protease and its substrate requires IVB signaling. We further show that the cytoplasmic cystathionine-ß-synthase (CBS) domain of the B protease is not critical for this interaction or for pro-σK processing, suggesting the IVB-dependent interaction site is in the membrane protease domain. Finally, we provide evidence that the B protease domain adopts both open and closed conformations in vivo. Collectively, our data support a substrate-gating model in which IVB-dependent cleavage of A on one side of the membrane triggers a conformational change in the membrane-embedded protease from a closed to an open state allowing pro-σK access to the caged interior of the protease.