Biophysical Correlates of Enhanced Immunogenicity of a Stabilized Variant of the Receptor Binding Domain of SARS-CoV-2.

The receptor binding domain (RBD) of SARS-CoV-2 is the primary target of neutralizing antibodies. We have previously reported the design and characterization of a mammalian cell expressed RBD derivative, mRBD1-3.2, that has higher thermal stability and greatly enhanced immunogenicity relative to the wild type mRBD. The protein is highly thermotolerant and immunogenic and is being explored for use in room temperature stable Covid-19 vaccine formulations. In the current study, we have investigated the folding pathway of both WT and stabilized RBD. It was found that chemical denaturation of RBD proceeds through a stable equilibrium intermediate. Thermal and chemical denaturation is reversible, as assayed by binding to the receptor ACE2. Unusually, in its native state, RBD binds to the hydrophobic probe ANS, and enhanced ANS binding is observed for the equilibrium intermediate state. Further characterization of the folding of mRBD1-3.2, both in solution and after reconstitution of lyophilized protein stored for a month at 37 °C, revealed a higher stability represented by higher Cm, faster refolding, slower unfolding, and enhanced resistance to proteolytic cleavage relative to WT. In contrast to WT RBD, the mutant showed decreased interaction with the hydrophobic moiety linoleic acid. Collectively, these data suggest that the enhanced immunogenicity results from reduced conformational fluctuations that likely enhance in vivo half-life as well as reduce the exposure of irrelevant non-neutralizing epitopes to the immune system.

Structural and Functional Characterization of Rv0792c from Mycobacterium tuberculosis: Identifying Small Molecule Inhibitor against HutC Protein.

In order to adapt in host tissues, microbial pathogens regulate their gene expression through a variety of transcription factors. Here, we have functionally characterized Rv0792c, a HutC homolog from Mycobacterium tuberculosis. In comparison to the parental strain, a strain of M. tuberculosis with a Rv0792c mutant was compromised for survival upon exposure to oxidative stress and infection in guinea pigs. RNA sequencing analysis revealed that Rv0792c regulates the expression of genes involved in stress adaptation and virulence of M. tuberculosis. Solution small-angle X-ray scattering (SAXS) data-steered model building confirmed that the C-terminal region plays a pivotal role in dimer formation. Systematic evolution of ligands by exponential enrichment (SELEX) resulted in the identification of single-strand DNA (ssDNA) aptamers that can be used as a tool to identify small-molecule inhibitors targeting Rv0792c. Using SELEX and SAXS data-based modeling, we identified residues essential for Rv0792c's aptamer binding activity. In this study, we also identified I-OMe-Tyrphostin as an inhibitor of Rv0792c's aptamer and DNA binding activity. The identified small molecule reduced the growth of intracellular M. tuberculosis in macrophages. The present study thus provides a detailed shape-function characterization of a HutC family of transcription factor from M. tuberculosis. IMPORTANCE Prokaryotes encode a large number of GntR family transcription factors that are involved in various fundamental biological processes, including stress adaptation and pathogenesis. Here, we investigated the structural and functional role of Rv0792c, a HutC homolog from M. tuberculosis. We demonstrated that Rv0792c is essential for M. tuberculosis to adapt to oxidative stress and establish disease in guinea pigs. Using a systematic evolution of ligands by exponential enrichment (SELEX) approach, we identified ssDNA aptamers from a random ssDNA library that bound to Rv0792c protein. These aptamers were thoroughly characterized using biochemical and biophysical assays. Using SAXS, we determined the structural model of Rv0792c in both the presence and absence of the aptamers. Further, using a combination of SELEX and SAXS methodologies, we identified I-OMe-Tyrphostin as a potential inhibitor of Rv0792c. Here we provide a detailed functional characterization of a transcription factor belonging to the HutC family from M. tuberculosis.

Ter-Seq: A high-throughput method to stabilize transient ternary complexes and measure associated kinetics.

Regulation of biological processes by proteins often involves the formation of transient, multimeric complexes whose characterization is mechanistically important but challenging. The bacterial toxin CcdB binds and poisons DNA Gyrase. The corresponding antitoxin CcdA extracts CcdB from its complex with Gyrase through the formation of a transient ternary complex, thus rejuvenating Gyrase. We describe a high throughput methodology called Ter-Seq to stabilize probable ternary complexes and measure associated kinetics using the CcdA-CcdB-GyrA14 ternary complex as a model system. The method involves screening a yeast surface display (YSD) saturation mutagenesis library of one partner (CcdB) for mutants that show enhanced ternary complex formation. We also isolated CcdB mutants that were either resistant or sensitive to rejuvenation, and used surface plasmon resonance (SPR) with purified proteins to validate the kinetics measured using the surface display. Positions, where CcdB mutations lead to slower rejuvenation rates, are largely involved in CcdA-binding, though there were several notable exceptions suggesting allostery. Mutations at these positions reduce the affinity towards CcdA, thereby slowing down the rejuvenation process. Mutations at GyrA14-interacting positions significantly enhanced rejuvenation rates, either due to reduced affinity or complete loss of CcdB binding to GyrA14. We examined the effect of different parameters (CcdA affinity, GyrA14 affinity, surface accessibilities, evolutionary conservation) on the rate of rejuvenation. Finally, we further validated the Ter-Seq results by monitoring the kinetics of ternary complex formation for individual CcdB mutants in solution by fluorescence resonance energy transfer (FRET) studies.

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.

Mechanistic insights into global suppressors of protein folding defects.

Most amino acid substitutions in a protein either lead to partial loss-of-function or are near neutral. Several studies have shown the existence of second-site mutations that can rescue defects caused by diverse loss-of-function mutations. Such global suppressor mutations are key drivers of protein evolution. However, the mechanisms responsible for such suppression remain poorly understood. To address this, we characterized multiple suppressor mutations both in isolation and in combination with inactive mutants. We examined six global suppressors of the bacterial toxin CcdB, the known M182T global suppressor of TEM-1 ß-lactamase, the N239Y global suppressor of p53-DBD and three suppressors of the SARS-CoV-2 spike Receptor Binding Domain. When coupled to inactive mutants, they promote increased in-vivo solubilities as well as regain-of-function phenotypes. In the case of CcdB, where novel suppressors were isolated, we determined the crystal structures of three such suppressors to obtain insight into the specific molecular interactions responsible for the observed effects. While most individual suppressors result in small stability enhancements relative to wildtype, which can be combined to yield significant stability increments, thermodynamic stabilisation is neither necessary nor sufficient for suppressor action. Instead, in diverse systems, we observe that individual global suppressors greatly enhance the foldability of buried site mutants, primarily through increase in refolding rate parameters measured in vitro. In the crowded intracellular environment, mutations that slow down folding likely facilitate off-pathway aggregation. We suggest that suppressor mutations that accelerate refolding can counteract this, enhancing the yield of properly folded, functional protein in vivo.

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.

Identification of stabilizing point mutations through mutagenesis of destabilized protein libraries.

Although there have been recent transformative advances in the area of protein structure prediction, prediction of point mutations that improve protein stability remains challenging. It is possible to construct and screen large mutant libraries for improved activity or ligand binding. However, reliable screens for mutants that improve protein stability do not yet exist, especially for proteins that are well folded and relatively stable. Here, we demonstrate that incorporation of a single, specific, destabilizing mutation termed parent inactivating mutation into each member of a single-site saturation mutagenesis library, followed by screening for suppressors, allows for robust and accurate identification of stabilizing mutations. We carried out fluorescence-activated cell sorting of such a yeast surface display, saturation suppressor library of the bacterial toxin CcdB, followed by deep sequencing of sorted populations. We found that multiple stabilizing mutations could be identified after a single round of sorting. In addition, multiple libraries with different parent inactivating mutations could be pooled and simultaneously screened to further enhance the accuracy of identification of stabilizing mutations. Finally, we show that individual stabilizing mutations could be combined to result in a multi-mutant that demonstrated an increase in thermal melting temperature of about 20 °C, and that displayed enhanced tolerance to high temperature exposure. We conclude that as this method is robust and employs small library sizes, it can be readily extended to other display and screening formats to rapidly isolate stabilized protein mutants.

Rapid Identification of Secondary Structure and Binding Site Residues in an Intrinsically Disordered Protein Segment.

Mycobacterium tuberculosis harbours nine toxin-antitoxin (TA) systems of the MazEF family. MazEF TA modules are of immense importance due to the perceived role of the MazF toxin in M. tuberculosis persistence and disease. The MazE antitoxin has a disordered C-terminal domain that binds the toxin, MazF and neutralizes its endoribonuclease activity. However, the structure of most MazEF TA complexes remains unsolved till date, obscuring structural and functional information about the antitoxins. We present a facile method to identify toxin binding residues on the disordered antitoxin. Charged residue scanning mutagenesis was used to screen a yeast surface displayed MazE6 antitoxin library against its purified cognate partner, the MazF6 toxin. Binding residues were deciphered by probing the relative reduction in binding to the ligand by flow cytometry. We have used this to identify putative antitoxin interface residues and local structure attained by the antitoxin upon interaction in the MazEF6 TA system and the same methodology is readily applicable to other intrinsically disordered protein regions.

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.

Mechanism of CcdA-Mediated Rejuvenation of DNA Gyrase.

Most biological processes involve formation of transient complexes where binding of a ligand allosterically modulates function. The ccd toxin-antitoxin system is involved in plasmid maintenance and bacterial persistence. The CcdA antitoxin accelerates dissociation of CcdB from its complex with DNA gyrase, binds and neutralizes CcdB, but the mechanistic details are unclear. Using a series of experimental and computational approaches, we demonstrate the formation of transient ternary and quaternary CcdA:CcdB:gyrase complexes and delineate the molecular steps involved in the rejuvenation process. Binding of region 61-72 of CcdA to CcdB induces the vital structural and dynamic changes required to facilitate dissociation from gyrase, region 50-60 enhances the dissociation process through additional allosteric effects, and segment 37-49 prevents gyrase rebinding. This study provides insights into molecular mechanisms responsible for recovery of CcdB-poisoned cells from a persister-like state. Similar methodology can be used to characterize other important transient, macromolecular complexes.

Facile measurement of protein stability and folding kinetics using a nano differential scanning fluorimeter.

With advancements in high-throughput generation of phenotypic data on mutant proteins, it has become important to individually characterize different proteins or their variants rapidly and with minimal sample consumption. We have made use of a nano differential scanning fluorimetric device, from NanoTemper technologies, to rapidly carry out isothermal chemical denaturation and measure folding/unfolding kinetics of proteins and compared these to corresponding data obtained from conventional spectrofluorimetry. We show that using sample volumes 10-50-fold lower than with conventional fluorimetric techniques, one can rapidly and accurately measure thermodynamic and kinetic stability, as well as folding/unfolding kinetics. This method also facilitates characterization of proteins that are difficult to express and purify.