Structural Basis of RNA Recognition by Mycobacterium tuberculosis MazF Homologues

The MazEF toxin-antitoxin (TA) system consists of the antitoxin MazE and the toxin MazF. MazF is a sequence-specific endoribonuclease that upon activation causes cellular growth arrest and increass the level of persisters. Moreover, MazF-induced cells are in a quasi-dormant state that cells remain metabolically active while stop dividing. The quasi-dormancy is similar to the nonreplicating state of M. tuberculosis during latent tuberculosis, thus suggesting the role of mazEF in M. tuberculosis dormancy and persistence. M. tuberculosis has nine mazEF TA modules, each with different RNA cleavage specificities and implicated in selective gene expression during stress conditions. To date only the Bacillus subtilis MazF-RNA complex structure has been determined. As M. tuberculosis MazF homologues recognize distinct RNA sequences, their molecular mechanisms of substrate specificity remain unclear. By taking advantage of X-ray crystallography, we have determined structures of two M. tuberculosis MazF-RNA complexes, MazF-mt1 (Rv2801c) and MazF-mt3 (Rv1991c) in complex with an uncleavable RNA substrate. These structures have provided the molecular basis of sequence-specific RNA recognition and cleavage by MazF toxins.

Both MazF-mt1-RNA and MazF-mt3-RNA complexes showed similar structural organization with one molecule of RNA bound to a MazF-mt1 or MazF-mt3 dimer and occupying the same pocket within the MazF dimer interface. Similar to B. subtilis MazF-RNA complex, MazF-mt1 and MazF-mt3 displayed a conserved active site architecture, where two highly conserved residues, Arg and Thr, form hydrogen bonds with the scissile phosphate group in the cleavage site of the bound RNA. The MazF-mt1-RNA complex also showed specific interactions with its three-base RNA recognition element. Compared with the B. subtilis MazF-RNA complex, our structures showed that residues involved in sequence-specific recognition of target RNA vary between the MazF homologues, therefore explaining the molecular basis for their different RNA recognition sequences. In addition, local conformational changes of the loops in the RNA binding site of MazF-mt1 appear to play a role in MazF targeting different RNA lengths and sequences. In contrast, the MazF-mt3-RNA complex is in a non-optimal RNA binding state with a symmetry-related MazF-mt3 molecule found to make interactions with the bound RNA in the crystal. The crystal-packing interactions were further examined by isothermal titration calorimetry (ITC) studies on selected MazF-mt3 mutants. Our attempts to utilize a MazF-mt3 mutant bearing mutations involved in crystal contacts all crystallized with few nucleotides, which are still found to interact with a symmetry mate. However, these different crystal forms revealed the conformational flexibility of loops in the RNA binding interface of MazF-mt3, suggesting their role in RNA binding and recognition, which will require further studies on additional MazF-mt3-RNA complex interactions.

In conclusion, the structures of the MazF-mt1-RNA and MazF-mt3-RNA complexes provide the first structural information on any M. tuberculosis MazF homologues. Supplemented with structure-guided mutational studies on MazF toxicity in vivo, this study has addressed the structural basis of different RNA cleavage specificities among MazF homologues. Our work will guide future studies on the function of other M. tuberculosis MazF and MazE-MazF homologues, and will help delineate their physiological roles in M. tuberculosis stress responses and pathogenesis.

RNA Virus Evolution: a Cross-scale, Disease Dynamic Perspective

RNA viruses are an important cause of global morbidity and mortality. The rapid evolutionary rates of RNA virus pathogens, caused by high replication rates and error-prone polymerases, can make the pathogens difficult to control. RNA viruses can undergo immune escape within their hosts and develop resistance to the treatment and vaccines we design to fight them. Understanding the spread and evolution of RNA pathogens is essential for reducing human suffering. In this dissertation, I make use of the rapid evolutionary rate of viral pathogens to answer several questions about how RNA viruses spread and evolve. To address each of the questions, I link mathematical techniques for modeling viral population dynamics with phylogenetic and coalescent techniques for analyzing and modeling viral genetic sequences and evolution. The first project uses multi-scale mechanistic modeling to show that decreases in viral substitution rates over the course of an acute infection, combined with the timing of infectious hosts transmitting new infections to susceptible individuals, can account for discrepancies in viral substitution rates in different host populations. The second project combines coalescent models with within-host mathematical models to identify driving evolutionary forces in chronic hepatitis C virus infection. The third project compares the effects of intrinsic and extrinsic viral transmission rate variation on viral phylogenies.

Occurrence and Function of Hoogsteen Base Pairs in Nucleic Acids

Nucleic acids (DNA and RNA) play essential roles in the central dogma of biology for the storage and transfer of genetic information. The unique chemical and conformational structures of nucleic acids – the double helix composed of complementary Watson-Crick base pairs, provide the structural basis to carry out their biological functions. DNA double helix can dynamically accommodate Watson-Crick and Hoogsteen base-pairing, in which the purine base is flipped by ~180° degrees to adopt syn rather than anti conformation as in Watson-Crick base pairs. There is growing evidence that Hoogsteen base pairs play important roles in DNA replication, recognition, damage or mispair accommodation and repair. Here, we constructed a database for existing Hoogsteen base pairs in DNA duplexes by a structure-based survey from the Protein Data Bank, and structural analyses based on the resulted Hoogsteen structures revealed that Hoogsteen base pairs occur in a wide variety of biological contexts and can induce DNA kinking towards the major groove. As there were documented difficulties in modeling Hoogsteen or Watson-Crick by crystallography, we collaborated with the Richardsons’ lab and identified potential Hoogsteen base pairs that were mis-modeled as Watson-Crick base pairs which suggested that Hoogsteen can be more prevalent than it was thought to be. We developed solution NMR method combined with the site-specific isotope labeling to characterize the formation of, or conformational exchange with Hoogsteen base pairs in large DNA-protein complexes under solution conditions, in the absence of the crystal packing force. We showed that there are enhanced chemical exchange, potentially between Watson-Crick and Hoogsteen, at a sharp kink site in the complex formed by DNA and the Integration Host Factor protein. In stark contrast to B-form DNA, we found that Hoogsteen base pairs are strongly disfavored in A-form RNA duplex. Chemical modifications N1-methyl adenosine and N1-methyl guanosine that block Watson-Crick base-pairing, can be absorbed as Hoogsteen base pairs in DNA, but rather potently destabilized A-form RNA and caused helix melting. The intrinsic instability of Hoogsteen base pairs in A-form RNA endows the N1-methylation as a functioning post-transcriptional modification that was known to facilitate RNA folding, translation and potentially play roles in the epitranscriptome. On the other hand, the dynamic property of DNA that can accommodate Hoogsteen base pairs could be critical to maintaining the genome stability.

A New Method for Modeling Free Surface Flows and Fluid-structure Interaction with Ocean Applications

The computational modeling of ocean waves and ocean-faring devices poses numerous challenges. Among these are the need to stably and accurately represent both the fluid-fluid interface between water and air as well as the fluid-structure interfaces arising between solid devices and one or more fluids. As techniques are developed to stably and accurately balance the interactions between fluid and structural solvers at these boundaries, a similarly pressing challenge is the development of algorithms that are massively scalable and capable of performing large-scale three-dimensional simulations on reasonable time scales. This dissertation introduces two separate methods for approaching this problem, with the first focusing on the development of sophisticated fluid-fluid interface representations and the second focusing primarily on scalability and extensibility to higher-order methods.

We begin by introducing the narrow-band gradient-augmented level set method (GALSM) for incompressible multiphase Navier-Stokes flow. This is the first use of the high-order GALSM for a fluid flow application, and its reliability and accuracy in modeling ocean environments is tested extensively. The method demonstrates numerous advantages over the traditional level set method, among these a heightened conservation of fluid volume and the representation of subgrid structures.

Next, we present a finite-volume algorithm for solving the incompressible Euler equations in two and three dimensions in the presence of a flow-driven free surface and a dynamic rigid body. In this development, the chief concerns are efficiency, scalability, and extensibility (to higher-order and truly conservative methods). These priorities informed a number of important choices: The air phase is substituted by a pressure boundary condition in order to greatly reduce the size of the computational domain, a cut-cell finite-volume approach is chosen in order to minimize fluid volume loss and open the door to higher-order methods, and adaptive mesh refinement (AMR) is employed to focus computational effort and make large-scale 3D simulations possible. This algorithm is shown to produce robust and accurate results that are well-suited for the study of ocean waves and the development of wave energy conversion (WEC) devices.