Molecular bases for strong phenotypic effects of single synonymous codon substitutions in the E. coli ccdB toxin gene.

BACKGROUND: Single synonymous codon mutations typically have only minor or no effects on gene function. Here, we estimate the effects on cell growth of ~ 200 single synonymous codon mutations in an operonic context by mutating almost all positions of ccdB, the 101-residue long cytotoxin of the ccdAB Toxin-Antitoxin (TA) operon to most degenerate codons. Phenotypes were assayed by transforming the mutant library into CcdB sensitive and resistant E. coli strains, isolating plasmid pools, and subjecting them to deep sequencing. Since autoregulation is a hallmark of TA operons, phenotypes obtained for ccdB synonymous mutants after transformation in a RelE toxin reporter strain followed by deep sequencing provided information on the amount of CcdAB complex formed. RESULTS: Synonymous mutations in the N-terminal region involved in translation initiation showed the strongest non-neutral phenotypic effects. We observe an interplay of numerous factors, namely, location of the codon, codon usage, t-RNA abundance, formation of anti-Shine Dalgarno sequences, predicted transcript secondary structure, and evolutionary conservation in determining phenotypic effects of ccdB synonymous mutations. Incorporation of an N-terminal, hyperactive synonymous mutation, in the background of the single synonymous codon mutant library sufficiently increased translation initiation, such that mutational effects on either folding or termination of translation became more apparent. Introduction of putative pause sites not only affects the translational rate, but might also alter the folding kinetics of the protein in vivo. CONCLUSION: In summary, the study provides novel insights into diverse mechanisms by which synonymous mutations modulate gene function. This information is useful in optimizing heterologous gene expression in E. coli and understanding the molecular bases for alteration in gene expression that arise due to synonymous mutations.

Functional and Biochemical Characterization of the MazEF6 Toxin-Antitoxin System of Mycobacterium tuberculosis.

The Mycobacterium tuberculosis genome harbors nine toxin-antitoxin (TA) systems that are members of the mazEF family, unlike other prokaryotes, which have only one or two. Although the overall tertiary folds of MazF toxins are predicted to be similar, it is unclear how they recognize structurally different RNAs and antitoxins with divergent sequence specificity. Here, we have expressed and purified the individual components and complex of the MazEF6 TA system from M. tuberculosis. Size exclusion chromatography-multiangle light scattering (SEC-MALS) was performed to determine the oligomerization status of the toxin, antitoxin, and the complex in different stoichiometric ratios. The relative stabilities of the proteins were determined by nano-differential scanning fluorimetry (nano-DSF). Microscale thermophoresis (MST) and yeast surface display (YSD) were performed to measure the relative affinities between the cognate toxin-antitoxin partners. The interaction between MazEF6 complexes and cognate promoter DNA was also studied using MST. Analysis of paired-end RNA sequencing data revealed that the overexpression of MazF6 resulted in differential expression of 323 transcripts in M. tuberculosis. Network analysis was performed to identify the nodes from the top-response network. The analysis of mRNA protection ratios resulted in identification of putative MazF6 cleavage site in its native host, M. tuberculosis. IMPORTANCE M. tuberculosis harbors a large number of type II toxin-antitoxin (TA) systems, the exact roles for most of which are unclear. Prior studies have reported that overexpression of several of these type II toxins inhibits bacterial growth and contributes to the formation of drug-tolerant populations in vitro. To obtain insights into M. tuberculosis MazEF6 type II TA system function, we determined stability, oligomeric states, and binding affinities of cognate partners with each other and with their promoter operator DNA. Using RNA-seq data obtained from M. tuberculosis overexpression strains, we have identified putative MazF6 cleavage sites and targets in its native, cellular context.

Design of a highly thermotolerant, immunogenic SARS-CoV-2 spike fragment.

Virtually all SARS-CoV-2 vaccines currently in clinical testing are stored in a refrigerated or frozen state prior to use. This is a major impediment to deployment in resource-poor settings. Furthermore, several of them use viral vectors or mRNA. In contrast to protein subunit vaccines, there is limited manufacturing expertise for these nucleic-acid-based modalities, especially in the developing world. Neutralizing antibodies, the clearest known correlate of protection against SARS-CoV-2, are primarily directed against the receptor-binding domain (RBD) of the viral spike protein, suggesting that a suitable RBD construct might serve as a more accessible vaccine ingredient. We describe a monomeric, glycan-engineered RBD protein fragment that is expressed at a purified yield of 214 mg/l in unoptimized, mammalian cell culture and, in contrast to a stabilized spike ectodomain, is tolerant of exposure to temperatures as high as 100 °C when lyophilized, up to 70 °C in solution and stable for over 4 weeks at 37 °C. In prime:boost guinea pig immunizations, when formulated with the MF59-like adjuvant AddaVax, the RBD derivative elicited neutralizing antibodies with an endpoint geometric mean titer of ∼415 against replicative virus, comparing favorably with several vaccine formulations currently in the clinic. These features of high yield, extreme thermotolerance, and satisfactory immunogenicity suggest that such RBD subunit vaccine formulations hold great promise to combat COVID-19.

Protein model discrimination attempts using mutational sensitivity, predicted secondary structure, and model quality information.

Structure prediction methods often generate a large number of models for a target sequence. Even if the correct fold for the target sequence is sampled in this dataset, it is difficult to distinguish it from other decoy structures. An attempt to solve this problem using experimental mutational sensitivity data for the CcdB protein was described previously by exploiting the correlation of residue depth with mutational sensitivity (r ~ 0.6). We now show that such a correlation extends to four other proteins with localized active sites, and for which saturation mutagenesis datasets exist. We also examine whether incorporation of predicted secondary structure information and the DOPE model quality assessment score, in addition to mutational sensitivity, improves the accuracy of model discrimination using a decoy dataset of 163 targets from CASP. Although most CASP models would have been subjected to model quality assessment prior to submission, we find that the DOPE score makes a substantial contribution to the observed improvement. We therefore also applied the approach to CcdB and four other proteins for which reliable experimental mutational data exist and observe that inclusion of experimental mutational data results in a small qualitative improvement in model discrimination relative to that seen with just the DOPE score. This is largely because of our limited ability to quantitatively predict effects of point mutations on in vivo protein activity. Further improvements in the methodology are required to facilitate improved utilization of single mutant data.

Combining cysteine scanning with chemical labeling to map protein-protein interactions and infer bound structure in an intrinsically disordered region.

The Mycobacterium tuberculosis genome harbours nine toxin-antitoxin (TA) systems of the mazEF family. These consist of two proteins, a toxin and an antitoxin, encoded in an operon. While the toxin has a conserved fold, the antitoxins are structurally diverse and the toxin binding region is typically intrinsically disordered before binding. We describe high throughput methodology for accurate mapping of interfacial residues and apply it to three MazEF complexes. The method involves screening one partner protein against a panel of chemically masked single cysteine mutants of its interacting partner, displayed on the surface of yeast cells. Such libraries have much lower diversity than those generated by saturation mutagenesis, simplifying library generation and data analysis. Further, because of the steric bulk of the masking reagent, labeling of virtually all exposed epitope residues should result in loss of binding, and buried residues are inaccessible to the labeling reagent. The binding residues are deciphered by probing the loss of binding to the labeled cognate partner by flow cytometry. Using this methodology, we have identified the interfacial residues for MazEF3, MazEF6 and MazEF9 TA systems of M. tuberculosis. In the case of MazEF9, where a crystal structure was available, there was excellent agreement between our predictions and the crystal structure, superior to those with AlphaFold2. We also report detailed biophysical characterization of the MazEF3 and MazEF9 TA systems and measured the relative affinities between cognate and non-cognate toxin-antitoxin partners in order to probe possible cross-talk between these systems.

Prediction of Function Determining and Buried Residues Through Analysis of Saturation Mutagenesis Datasets.

Mutational scanning can be used to probe effects of large numbers of point mutations on protein function. Positions affected by mutation are primarily at either buried or at exposed residues directly involved in function, hereafter designated as active-site residues. In the absence of prior structural information, it has not been easy to distinguish between these two categories of residues. We curated and analyzed a set of twelve published deep mutational scanning datasets. The analysis revealed differential patterns of mutational sensitivity and substitution preferences at buried and exposed positions. Prediction of buried-sites solely from the mutational sensitivity data was facilitated by incorporating predicted sequence-based accessibility values. For active-site residues we observed mean sensitivity, specificity and accuracy of 61, 90 and 88% respectively. For buried residues the corresponding figures were 59, 90 and 84% while for exposed non active-site residues these were 98, 44 and 82% respectively. We also identified positions which did not follow these general trends and might require further experimental re-validation. This analysis highlights the ability of deep mutational scans to provide important structural and functional insights, even in the absence of three-dimensional structures determined using conventional structure determination techniques, and also discuss some limitations of the methodology.

Prediction of Residue-specific Contributions to Binding and Thermal Stability Using Yeast Surface Display.

Accurate prediction of residue burial as well as quantitative prediction of residue-specific contributions to protein stability and activity is challenging, especially in the absence of experimental structural information. This is important for prediction and understanding of disease causing mutations, and for protein stabilization and design. Using yeast surface display of a saturation mutagenesis library of the bacterial toxin CcdB, we probe the relationship between ligand binding and expression level of displayed protein, with in vivo solubility in E. coli and in vitro thermal stability. We find that both the stability and solubility correlate well with the total amount of active protein on the yeast cell surface but not with total amount of expressed protein. We coupled FACS and deep sequencing to reconstruct the binding and expression mean fluorescent intensity of each mutant. The reconstructed mean fluorescence intensity (MFIseq) was used to differentiate between buried site, exposed non active-site and exposed active-site positions with high accuracy. The MFIseq was also used as a criterion to identify destabilized as well as stabilized mutants in the library, and to predict the melting temperatures of destabilized mutants. These predictions were experimentally validated and were more accurate than those of various computational predictors. The approach was extended to successfully identify buried and active-site residues in the receptor binding domain of the spike protein of SARS-CoV-2, suggesting it has general applicability.

VapC21 Toxin Contributes to Drug-Tolerance and Interacts With Non-cognate VapB32 Antitoxin in Mycobacterium tuberculosis.

The prokaryotic ubiquitous Toxin-antitoxin (TA) modules encodes for a stable toxin and an unstable antitoxin. VapBC subfamily is the most abundant Type II TA system in M. tuberculosis genome. However, the exact physiological role for most of these Type II TA systems are still unknown. Here, we have comprehensively characterized the VapBC21 TA locus from M. tuberculosis. The overexpression of VapC21 inhibited mycobacterial growth in a bacteriostatic manner and as expected, growth inhibition was abrogated upon co-expression of the cognate antitoxin, VapB21. We observed that the deletion of vapC21 had no noticeable influence on the in vitro and in vivo growth of M. tuberculosis. Using co-expression and biophysical studies, we observed that in addition to VapB21, VapC21 is also able to interact with non-cognate antitoxin, VapB32. The strength of interaction varied between the cognate and non-cognate TA pairs. The overexpression of VapC21 resulted in differential expression of approximately 435 transcripts in M. tuberculosis. The transcriptional profiles obtained upon ectopic expression of VapC21 was similar to those reported in M. tuberculosis upon exposure to stress conditions such as nutrient starvation and enduring hypoxic response. Further, VapC21 overexpression also led to increased expression of WhiB7 regulon and bacterial tolerance to aminoglycosides and ethambutol. Taken together, these results indicate that a complex network of interactions exists between non-cognate TA pairs and VapC21 contributes to drug tolerance in vitro.