Co-Mutations and Possible Variation Tendency of the Spike RBD and Membrane Protein in SARS-CoV-2 by Machine Learning.

Since the onset of the coronavirus disease 2019 (COVID-19) pandemic, SARS-CoV-2 variants capable of breakthrough infections have attracted global attention. These variants have significant mutations in the receptor-binding domain (RBD) of the spike protein and the membrane (M) protein, which may imply an enhanced ability to evade immune responses. In this study, an examination of co-mutations within the spike RBD and their potential correlation with mutations in the M protein was conducted. The EVmutation method was utilized to analyze the distribution of the mutations to elucidate the relationship between the mutations in the spike RBD and the alterations in the M protein. Additionally, the Sequence-to-Sequence Transformer Model (S2STM) was employed to establish mapping between the amino acid sequences of the spike RBD and M proteins, offering a novel and efficient approach for streamlined sequence analysis and the exploration of their interrelationship. Certain mutations in the spike RBD, G339D-S373P-S375F and Q493R-Q498R-Y505, are associated with a heightened propensity for inducing mutations at specific sites within the M protein, especially sites 3 and 19/63. These results shed light on the concept of mutational synergy between the spike RBD and M proteins, illuminating a potential mechanism that could be driving the evolution of SARS-CoV-2.

Deep mutational scanning of the RNase III-like domain in Trypanosoma brucei RNA editing protein KREPB4.

Kinetoplastid pathogens including Trypanosoma brucei, T. cruzi, and Leishmania species, are early diverged, eukaryotic, unicellular parasites. Functional understanding of many proteins from these pathogens has been hampered by limited sequence homology to proteins from other model organisms. Here we describe the development of a high-throughput deep mutational scanning approach in T. brucei that facilitates rapid and unbiased assessment of the impacts of many possible amino acid substitutions within a protein on cell fitness, as measured by relative cell growth. The approach leverages several molecular technologies: cells with conditional expression of a wild-type gene of interest and constitutive expression of a library of mutant variants, degron-controlled stabilization of I-SceI meganuclease to mediate highly efficient transfection of a mutant allele library, and a high-throughput sequencing readout for cell growth upon conditional knockdown of wild-type gene expression and exclusive expression of mutant variants. Using this method, we queried the effects of amino acid substitutions in the apparently non-catalytic RNase III-like domain of KREPB4 (B4), which is an essential component of the RNA Editing Catalytic Complexes (RECCs) that carry out mitochondrial RNA editing in T. brucei. We measured the impacts of thousands of B4 variants on bloodstream form cell growth and validated the most deleterious variants containing single amino acid substitutions. Crucially, there was no correlation between phenotypes and amino acid conservation, demonstrating the greater power of this method over traditional sequence homology searching to identify functional residues. The bloodstream form cell growth phenotypes were combined with structural modeling, RECC protein proximity data, and analysis of selected substitutions in procyclic form T. brucei. These analyses revealed that the B4 RNaseIII-like domain is essential for maintenance of RECC integrity and RECC protein abundances and is also involved in changes in RECCs that occur between bloodstream and procyclic form life cycle stages.

Activation of Cryptic Donor Splice Sites within the UDP-Glucuronosyltransferase (UGT)1A First-Exon Region Generates Variant Transcripts That Encode UGT1A Proteins with Truncated Aglycone-Binding Domains.

The human UDP-glucuronosyltransferases (UGTs) have crucial roles in metabolizing and clearing numerous small lipophilic compounds. The UGT1A locus generates nine UGT1A mRNAs, 65 spliced transcripts, and 34 circular RNAs. In this study, our analysis of published UGT-RNA capture sequencing (CaptureSeq) datasets identified novel splice junctions that predict 24 variant UGT1A transcripts derived from ligation of exon 2 to unique sequences within the UGT1A first-exon region using cryptic donor splice sites. Of these variants, seven (1A1_n1, 1A3_n3, 1A4_n4, 1A5_n1, 1A8_n2, 1A9_n2, 1A10_n7) are predicted to encode UGT1A proteins with truncated aglycone-binding domains. We assessed their expression profiles and deregulation in cancer using four RNA sequencing (RNA-Seq) datasets of paired normal and cancerous drug-metabolizing tissues from large patient cohorts. Variants were generally coexpressed with their canonical counterparts with a higher relative abundance in tumor than in normal tissues. Variants showed tissue-specific expression with high interindividual variability but overall low abundance. However, 1A8_n2 showed high abundance in normal and cancerous colorectal tissues, with levels that approached or surpassed canonical 1A8 mRNA levels in many samples. We cloned 1A8_n2 and showed expression of the predicted protein (1A8_i3) in human embryonic kidney (HEK)293T cells. Glucuronidation assays with 4-methylumbelliferone (4MU) showed that 1A8_i3 had no activity and was unable to inhibit the activity of 1A8_i1 protein. In summary, the activation of cryptic donor splice sites within the UGT1A first-exon region expands the UGT1A transcriptome and proteome. The 1A8_n2 cryptic donor splice site is highly active in colorectal tissues, representing an important cis-regulatory element that negatively regulates the function of the UGT1A8 gene through pre-mRNA splicing. SIGNIFICANT STATEMENT: The UGT1A locus generates nine canonical mRNAs, 65 alternately spliced transcripts, and 34 different circular RNAs. The present study reports a series of novel UDP-glucuronosyltransferase (UGT)1A variants resulting from use of cryptic donor splice sites in both normal and cancerous tissues, several of which are predicted to encode variant UGT1A proteins with truncated aglycone-binding domains. Of these, 1A8_n2 shows exceptionally high abundance in colorectal tissues, highlighting its potential role in the first-pass metabolism in gut through the glucuronidation pathway.

Substrate-induced condensation activates plant TIR domain proteins.

Plant nucleotide-binding leucine-rich repeat (NLR) immune receptors with an N-terminal Toll/interleukin-1 receptor (TIR) domain mediate recognition of strain-specific pathogen effectors, typically via their C-terminal ligand-sensing domains1. Effector binding enables TIR-encoded enzymatic activities that are required for TIR-NLR (TNL)-mediated immunity2,3. Many truncated TNL proteins lack effector-sensing domains but retain similar enzymatic and immune activities4,5. The mechanism underlying the activation of these TIR domain proteins remain unclear. Here we show that binding of the TIR substrates NAD+ and ATP induces phase separation of TIR domain proteins in vitro. A similar condensation occurs with a TIR domain protein expressed via its native promoter in response to pathogen inoculation in planta. The formation of TIR condensates is mediated by conserved self-association interfaces and a predicted intrinsically disordered loop region of TIRs. Mutations that disrupt TIR condensates impair the cell death activity of TIR domain proteins. Our data reveal phase separation as a mechanism for the activation of TIR domain proteins and provide insight into substrate-induced autonomous activation of TIR signalling to confer plant immunity.

Extracellular domain mutations of the EGF receptor differentially modulate high-affinity and low-affinity responses to EGF receptor ligands.

The EGF receptor is mutated in a number of cancers. In most cases, the mutations occur in the intracellular tyrosine kinase domain. However, in glioblastomas, many of the mutations are in the extracellular ligand binding domain. To determine what changes in receptor function are induced by such extracellular domain mutations, we analyzed the binding and biological response to the seven different EGF receptor ligands in three common glioblastoma mutants-R84K, A265V, and G574V. Our data indicate that all three mutations significantly increase the binding affinity of all seven ligands. In addition, the mutations increase the potency of all ligands for stimulating receptor autophosphorylation, phospholipase Cγ, Akt, and MAP kinase activity. In all mutants, the rank order of ligand potency seen at the wild-type receptor was retained, suggesting that the receptors still discriminate among the different ligands. However, the low-affinity ligands, EPR and EPG, did show larger than average enhancements of potency for stimulating Akt and MAPK but not receptor autophosphorylation and phospholipase Cγ activation. Relative to the wild-type receptor, these changes lead to an increase in the responsiveness of these mutants to physiological concentrations of ligands and an alteration in the ratio of activation of the different pathways. This may contribute to their oncogenic potential. In the context of recent findings, our data also suggest that so-called "high"-affinity biological responses arise from activation by isolated receptor dimers, whereas "low"-affinity biological responses require clustering of receptors which occurs at higher concentrations of ligand.

Trp207 regulation of voltage-dependent activation of human Hv1 proton channel.

In voltage-gated Na+ and K+ channels, the hydrophobicity of noncharged residues in the S4 helix has been shown to regulate the S4 movement underlying the process of voltage-sensing domain (VSD) activation. In voltage-gated proton channel Hv1, there is a bulky noncharged tryptophan residue located at the S4 transmembrane segment. This tryptophan remains entirely conserved across all Hv1 members but is not seen in other voltage-gated ion channels, indicating that the tryptophan contributes different roles in VSD activation. The conserved tryptophan of human voltage-gated proton channel Hv1 is Trp207 (W207). Here, we showed that W207 modifies human Hv1 voltage-dependent activation, and small residues replacement at position 207 strongly perturbs Hv1 channel opening and closing, and the size of the side chain instead of the hydrophobic group of W207 regulates the transition between closed and open states of the channel. We conclude that the large side chain of tryptophan controls the energy barrier during the Hv1 VSD transition.

Missense mutations in the central domains of cardiac myosin binding protein-C and their potential contribution to hypertrophic cardiomyopathy.

Myosin binding protein-C (MyBP-C) is a multidomain protein that regulates muscle contraction. Mutations in MYBPC3, the gene encoding for the cardiac variant (henceforth called cMyBP-C), are amongst the most frequent causes of hypertrophic cardiomyopathy. Most mutations lead to a truncated version of cMyBP-C, which is most likely unstable. However, missense mutations have also been reported, which tend to cluster in the central domains of the cMyBP-C molecule. This suggests that these central domains are more than just a passive spacer between the better characterized N- and C-terminal domains. Here, we investigated the potential impact of four different missense mutations, E542Q, G596R, N755K, and R820Q, which are spread over the domains C3 to C6, on the function of MyBP-C on both the isolated protein level and in cardiomyocytes in vitro. Effect on domain stability, interaction with thin filaments, binding to myosin, and subcellular localization behavior were assessed. Our studies show that these missense mutations result in slightly different phenotypes at the molecular level, which are mutation specific. The expected functional readout of each mutation provides a valid explanation for why cMyBP-C fails to work as a brake in the regulation of muscle contraction, which eventually results in a hypertrophic cardiomyopathy phenotype. We conclude that missense mutations in cMyBP-C must be evaluated in context of their domain localization, their effect on interaction with thin filaments and myosin, and their effect on protein stability to explain how they lead to disease.

A regulatory role for the unstructured C-terminal domain of the CtBP transcriptional corepressor.

The C-terminal binding protein (CtBP) is a transcriptional corepressor that plays critical roles in development, tumorigenesis, and cell fate. CtBP proteins are structurally similar to alpha hydroxyacid dehydrogenases and feature a prominent intrinsically disordered region in the C terminus. In the mammalian system, CtBP proteins lacking the C-terminal domain (CTD) are able to function as transcriptional regulators and oligomerize, putting into question the significance of this unstructured domain for gene regulation. Yet, the presence of an unstructured CTD of ∼100 residues, including some short motifs, is conserved across Bilateria, indicating the importance of maintaining this domain over evolutionary time. To uncover the significance of the CtBP CTD, we functionally tested naturally occurring Drosophila isoforms of CtBP that possess or lack the CTD, namely CtBP(L) and CtBP(S). We used the CRISPRi system to recruit dCas9-CtBP(L) and dCas9-CtBP(S) to endogenous promoters to directly compare their transcriptional impacts in vivo. Interestingly, CtBP(S) was able to significantly repress transcription of the Mpp6 promoter, while CtBP(L) was much weaker, suggesting that the long CTD may modulate CtBP's repression activity. In contrast, in cell culture, the isoforms behaved similarly on a transfected Mpp6 reporter gene. The context-specific differences in activity of these two developmentally regulated isoforms suggests that the CTD may help provide a spectrum of repression activity suitable for developmental programs.

Identification and characterization of calreticulin as a novel plasminogen receptor.

Calreticulin (CRT) was originally identified as a key calcium-binding protein of the endoplasmic reticulum. Subsequently, CRT was shown to possess multiple intracellular functions, including roles in calcium homeostasis and protein folding. Recently, several extracellular functions have been identified for CRT, including roles in cancer cell invasion and phagocytosis of apoptotic and cancer cells by macrophages. In the current report, we uncover a novel function for extracellular CRT and report that CRT functions as a plasminogen-binding receptor that regulates the conversion of plasminogen to plasmin. We show that human recombinant or bovine tissue-derived CRT dramatically stimulated the conversion of plasminogen to plasmin by tissue plasminogen activator or urokinase-type plasminogen activator. Surface plasmon resonance analysis revealed that CRT-bound plasminogen (KD = 1.8 µM) with moderate affinity. Plasminogen binding and activation by CRT were inhibited by ε-aminocaproic acid, suggesting that an internal lysine residue of CRT interacts with plasminogen. We subsequently show that clinically relevant CRT variants (lacking four or eight lysines in carboxyl-terminal region) exhibited decreased plasminogen activation. Furthermore, CRT-deficient fibroblasts generated 90% less plasmin and CRT-depleted MDA MB 231 cells also demonstrated a significant reduction in plasmin generation. Moreover, treatment of fibroblasts with mitoxantrone dramatically stimulated plasmin generation by WT but not CRT-deficient fibroblasts. Our results suggest that CRT is an important cellular plasminogen regulatory protein. Given that CRT can empower cells with plasmin proteolytic activity, this discovery may provide new mechanistic insight into the established role of CRT in cancer.

The heparin-binding hemagglutinin protein of Mycobacterium tuberculosis is a nucleoid-associated protein.

Nucleoid-associated proteins (NAPs) regulate multiple cellular processes such as gene expression, virulence, and dormancy throughout bacterial species. NAPs help in the survival and adaptation of Mycobacterium tuberculosis (Mtb) within the host. Fourteen NAPs have been identified in Escherichia coli; however, only seven NAPs are documented in Mtb. Given its complex lifestyle, it is reasonable to assume that Mtb would encode for more NAPs. Using bioinformatics tools and biochemical experiments, we have identified the heparin-binding hemagglutinin (HbhA) protein of Mtb as a novel sequence-independent DNA-binding protein which has previously been characterized as an adhesion molecule required for extrapulmonary dissemination. Deleting the carboxy-terminal domain of HbhA resulted in a complete loss of its DNA-binding activity. Atomic force microscopy showed HbhA-mediated architectural modulations in the DNA, which may play a regulatory role in transcription and genome organization. Our results showed that HbhA colocalizes with the nucleoid region of Mtb. Transcriptomics analyses of a hbhA KO strain revealed that it regulates the expression of ∼36% of total and ∼29% of essential genes. Deletion of hbhA resulted in the upregulation of ∼73% of all differentially expressed genes, belonging to multiple pathways suggesting it to be a global repressor. The results show that HbhA is a nonessential NAP regulating gene expression globally and acting as a plausible transcriptional repressor.

Cyclase-associated protein interacts with actin filament barbed ends to promote depolymerization and formin displacement.

Cyclase-associated protein (CAP) has emerged as a central player in cellular actin turnover, but its molecular mechanisms of action are not yet fully understood. Recent studies revealed that the N terminus of CAP interacts with the pointed ends of actin filaments to accelerate depolymerization in conjunction with cofilin. Here, we use in vitro microfluidics-assisted TIRF microscopy to show that the C terminus of CAP promotes depolymerization at the opposite (barbed) ends of actin filaments. In the absence of actin monomers, full-length mouse CAP1 and C-terminal halves of CAP1 (C-CAP1) and CAP2 (C-CAP2) accelerate barbed end depolymerization. Using mutagenesis and structural modeling, we show that these activities are mediated by the WH2 and CARP domains of CAP. In addition, we observe that CAP collaborates with profilin to accelerate barbed end depolymerization and that these effects depend on their direct interaction, providing the first known example of CAP-profilin collaborative effects in regulating actin. In the presence of actin monomers, CAP1 attenuates barbed end growth and promotes formin dissociation. Overall, these findings demonstrate that CAP uses distinct domains and mechanisms to interact with opposite ends of actin filaments and drive turnover. Further, they contribute to the emerging view of actin barbed ends as sites of dynamic molecular regulation, where numerous proteins compete and cooperate with each other to tune polymer dynamics, similar to the rich complexity seen at microtubule ends.

Sequence tolerance of immunoglobulin variable domain framework regions to noncanonical intradomain disulfide linkages.

Most immunoglobulin (Ig) domains bear only a single highly conserved canonical intradomain, inter-ß-sheet disulfide linkage formed between Cys23-Cys104, and incorporation of rare noncanonical disulfide linkages at other locations can enhance Ig domain stability. Here, we exhaustively surveyed the sequence tolerance of Ig variable (V) domain framework regions (FRs) to noncanonical disulfide linkages. Starting from a destabilized VH domain lacking a Cys23-Cys104 disulfide linkage, we generated and screened phage-displayed libraries of engineered VHs, bearing all possible pairwise combinations of Cys residues in neighboring ß-strands of the Ig fold FRs. This approach identified seven novel Cys pairs in VH FRs (Cys4-Cys25, Cys4-Cys118, Cys5-Cys120, Cys6-Cys119, Cys22-Cys88, Cys24-Cys86, and Cys45-Cys100; the international ImMunoGeneTics information system numbering), whose presence rescued domain folding and stability. Introduction of a subset of these noncanonical disulfide linkages (three intra-ß-sheet: Cys4-Cys25, Cys22-Cys88, and Cys24-Cys86, and one inter-ß-sheet: Cys6-Cys119) into a diverse panel of VH, VL, and VHH domains enhanced their thermostability and protease resistance without significantly impacting expression, solubility, or binding to cognate antigens. None of the noncanonical disulfide linkages identified were present in the natural human VH repertoire. These data reveal an unexpected permissiveness of Ig V domains to noncanonical disulfide linkages at diverse locations in FRs, absent in the human repertoire, whose presence is compatible with antigen recognition and improves domain stability. Our work represents the most complete assessment to date of the role of engineered noncanonical disulfide bonding within FRs in Ig V domain structure and function.

A Molecular Analysis of the Aminopeptidase P-Related Domain of PID-5 from Caenorhabditis elegans.

A novel protein, PID-5, has been shown to be a requirement for germline immortality and has recently been implicated in RNA-induced epigenetic silencing in the Caenorhabditis elegans embryo. Importantly, it has been shown to contain both an eTudor and aminopeptidase P-related domain. However, the silencing mechanism has not yet been fully characterised. In this study, bioinformatic tools were used to compare pre-existing aminopeptidase P molecular structures to the AlphaFold2-predicted aminopeptidase P-related domain of PID-5 (PID-5 APP-RD). Structural homology, metal composition, inhibitor-bonding interactions, and the potential for dimerisation were critically assessed through computational techniques, including structural superimposition and protein-ligand docking. Results from this research suggest that the metallopeptidase-like domain shares high structural homology with known aminopeptidase P enzymes and possesses the canonical 'pita-bread fold'. However, the absence of conserved metal-coordinating residues indicates that only a single Zn2+ may be bound at the active site. The PID-5 APP-RD may form transient interactions with a known aminopeptidase P inhibitor and may therefore recognise substrates in a comparable way to the known structures. However, loss of key catalytic residues suggests the domain will be inactive. Further evidence suggests that heterodimerisation with C. elegans aminopeptidase P is feasible and therefore PID-5 is predicted to regulate proteolytic cleavage in the silencing pathway. PID-5 may interact with PID-2 to bring aminopeptidase P activity to the Z-granule, where it could influence WAGO-4 activity to ensure the balanced production of 22G-RNA signals for transgenerational silencing. Targeted experiments into APPs implicated in malaria and cancer are required in order to build upon the biological and therapeutic significance of this research.

Mega-scale experimental analysis of protein folding stability in biology and design.

Advances in DNA sequencing and machine learning are providing insights into protein sequences and structures on an enormous scale1. However, the energetics driving folding are invisible in these structures and remain largely unknown2. The hidden thermodynamics of folding can drive disease3,4, shape protein evolution5-7 and guide protein engineering8-10, and new approaches are needed to reveal these thermodynamics for every sequence and structure. Here we present cDNA display proteolysis, a method for measuring thermodynamic folding stability for up to 900,000 protein domains in a one-week experiment. From 1.8 million measurements in total, we curated a set of around 776,000 high-quality folding stabilities covering all single amino acid variants and selected double mutants of 331 natural and 148 de novo designed protein domains 40-72 amino acids in length. Using this extensive dataset, we quantified (1) environmental factors influencing amino acid fitness, (2) thermodynamic couplings (including unexpected interactions) between protein sites, and (3) the global divergence between evolutionary amino acid usage and protein folding stability. We also examined how our approach could identify stability determinants in designed proteins and evaluate design methods. The cDNA display proteolysis method is fast, accurate and uniquely scalable, and promises to reveal the quantitative rules for how amino acid sequences encode folding stability.

FREEDA: An automated computational pipeline guides experimental testing of protein innovation.

Cell biologists typically focus on conserved regions of a protein, overlooking innovations that can shape its function over evolutionary time. Computational analyses can reveal potential innovations by detecting statistical signatures of positive selection that lead to rapid accumulation of beneficial mutations. However, these approaches are not easily accessible to non-specialists, limiting their use in cell biology. Here, we present an automated computational pipeline FREEDA that provides a simple graphical user interface requiring only a gene name; integrates widely used molecular evolution tools to detect positive selection in rodents, primates, carnivores, birds, and flies; and maps results onto protein structures predicted by AlphaFold. Applying FREEDA to >100 centromere proteins, we find statistical evidence of positive selection within loops and turns of ancient domains, suggesting innovation of essential functions. As a proof-of-principle experiment, we show innovation in centromere binding of mouse CENP-O. Overall, we provide an accessible computational tool to guide cell biology research and apply it to experimentally demonstrate functional innovation.

The evolution of plant NLR immune receptors and downstream signal components.

Along with the emergence of green plants on this planet one billion years ago, the nucleotide binding site leucine-rich repeat (NLR) gene family originated and diverged into at least three subclasses. Two of them, with either characterized N-terminal toll/interleukin-1 receptor (TIR) or coiled-coil (CC) domain, serve as major types of immune receptor of effector-triggered immunity (ETI) in plants, whereas the one having a N-terminal Resistance to powdery mildew8 (RPW8) domain, functions as signal transfer component to them. In this review, we briefly summarized the history of identification of diverse NLR subclasses across Viridiplantae lineages during the establishment of NLR category, and highlighted recent advances on the evolution of NLR genes and several key downstream signal components under the background of ecological adaption.

Large-scale mapping and mutagenesis of human transcriptional effector domains.

Human gene expression is regulated by more than 2,000 transcription factors and chromatin regulators1,2. Effector domains within these proteins can activate or repress transcription. However, for many of these regulators we do not know what type of effector domains they contain, their location in the protein, their activation and repression strengths, and the sequences that are necessary for their functions. Here, we systematically measure the effector activity of more than 100,000 protein fragments tiling across most chromatin regulators and transcription factors in human cells (2,047 proteins). By testing the effect they have when recruited at reporter genes, we annotate 374 activation domains and 715 repression domains, roughly 80% of which are new and have not been previously annotated3-5. Rational mutagenesis and deletion scans across all the effector domains reveal aromatic and/or leucine residues interspersed with acidic, proline, serine and/or glutamine residues are necessary for activation domain activity. Furthermore, most repression domain sequences contain sites for small ubiquitin-like modifier (SUMO)ylation, short interaction motifs for recruiting corepressors or are structured binding domains for recruiting other repressive proteins. We discover bifunctional domains that can both activate and repress, some of which dynamically split a cell population into high- and low-expression subpopulations. Our systematic annotation and characterization of effector domains provide a rich resource for understanding the function of human transcription factors and chromatin regulators, engineering compact tools for controlling gene expression and refining predictive models of effector domain function.

Differential Structural Features of Two Mutant ADAR1p150 Zα Domains Associated with Aicardi-Goutières Syndrome.

The Zα domain of ADARp150 is critical for proper Z-RNA substrate binding and is a key factor in the type-I interferon response pathway. Two point-mutations in this domain (N173S and P193A), which cause neurodegenerative disorders, are linked to decreased A-to-I editing in disease models. To understand this phenomenon at the molecular level, we biophysically and structurally characterized these two mutated domains, revealing that they bind Z-RNA with a decreased affinity. Less efficient binding to Z-RNA can be explained by structural changes in beta-wing, part of the Z-RNA-protein interface, and alteration of conformational dynamics of the proteins.

Structural basis for transcription factor ZBTB7A recognition of DNA and effects of ZBTB7A somatic mutations that occur in human acute myeloid leukemia.

ZBTB7A belongs to a small family of transcription factors having three members in humans (7A, 7B, and 7C). They share a BTB/POZ protein interaction domain at the amino end and a zinc-finger DNA-binding domain at the carboxyl end. They control the transcription of a wide range of genes, having varied functions in hematopoiesis, oncogenesis, and metabolism (in particular glycolysis). ZBTB7A-binding profiles at gene promoters contain a consensus G(a/c)CCC motif, followed by a CCCC sequence in some instances. Structural and mutational investigations suggest that DNA-specific contacts with the four-finger tandem array of ZBTB7A are formed sequentially, initiated from ZF1-ZF2 binding to G(a/c)CCC before spreading to ZF3-ZF4, which bind the DNA backbone and the 3' CCCC sequence, respectively. Here, we studied some mutations found in t(8;21)-positive acute myeloid leukemia patients that occur within the ZBTB7A DNA-binding domain. We determined that these mutations generally impair ZBTB7A DNA binding, with the most severe disruptions resulting from mutations in ZF1 and ZF2, and the least from a frameshift mutation in ZF3 that results in partial mislocalization. Information provided here on ZBTB7A-DNA interactions is likely applicable to ZBTB7B/C, which have overlapping functions with ZBTB7A in controlling primary metabolism.

Imprinted SARS-CoV-2 humoral immunity induces convergent Omicron RBD evolution.

Continuous evolution of Omicron has led to a rapid and simultaneous emergence of numerous variants that display growth advantages over BA.5 (ref. 1). Despite their divergent evolutionary courses, mutations on their receptor-binding domain (RBD) converge on several hotspots. The driving force and destination of such sudden convergent evolution and its effect on humoral immunity remain unclear. Here we demonstrate that these convergent mutations can cause evasion of neutralizing antibody drugs and convalescent plasma, including those from BA.5 breakthrough infection, while maintaining sufficient ACE2-binding capability. BQ.1.1.10 (BQ.1.1 + Y144del), BA.4.6.3, XBB and CH.1.1 are the most antibody-evasive strains tested. To delineate the origin of the convergent evolution, we determined the escape mutation profiles and neutralization activity of monoclonal antibodies isolated from individuals who had BA.2 and BA.5 breakthrough infections2,3. Owing to humoral immune imprinting, BA.2 and especially BA.5 breakthrough infection reduced the diversity of the neutralizing antibody binding sites and increased proportions of non-neutralizing antibody clones, which, in turn, focused humoral immune pressure and promoted convergent evolution in the RBD. Moreover, we show that the convergent RBD mutations could be accurately inferred by deep mutational scanning profiles4,5, and the evolution trends of BA.2.75 and BA.5 subvariants could be well foreseen through constructed convergent pseudovirus mutants. These results suggest that current herd immunity and BA.5 vaccine boosters may not efficiently prevent the infection of Omicron convergent variants.