Submitted by Anderson Silva (avargas@icict.fiocruz.br) on 2012-10-24T13:20:12Z No. of bitstreams: 1 chyntia_c_d_acosta_ioc_bcm_0014_2010.pdf: 6762215 bytes, checksum: eeda2289ddbc51f6797d64ab1bcdbead (MD5) / Made available in DSpace on 2012-10-24T13:20:12Z (GMT). No. of bitstreams: 1 chyntia_c_d_acosta_ioc_bcm_0014_2010.pdf: 6762215 bytes, checksum: eeda2289ddbc51f6797d64ab1bcdbead (MD5) Previous issue date: 2010 / CNPq Instituto de Investigaciones en Ciencias de la Salud-Universidad Nacional de Asunción (IICS-UNA) / Fundação Oswaldo Cruz.Instituto Oswaldo Cruz. Rio de Janeiro, RJ, Brasil / A tuberculose (TB) é uma importante causa de mortalidade, sendo Mycobacterium tuberculosis (Mtb) o agente etiológico. Globalmente a família “Latin American- Mediterranean” (LAM) é responsável por aproximadamente 15% dos casos de TB. Dentro desta família, recentemente foi descrito um novo polimorfismo caracterizado por uma deleção de 26,3 kb e designado como LSP RDRio. Esta linhagem é freqüente (30%) no Rio de Janeiro e dados preliminares sugerem que cepas pertencentes a este genótipo apresentam maior virulência. Neste contexto, é de grande importância avaliar a distribuição e transmissibilidade deste genótipo em nível regional e global. O presente estudo descritivo teve como objetivo, primeiramente, estudar a variabilidade genética entre isolados de Mtb do Paraguai, da Argentina e da Venezuela pelas técnicas de Spoligotyping e MIRU-VNTR 12loci, países onde existe pouca informação sobre a presença e natureza de genótipos de Mtb. Outro foco abordado foi o estudo da família LAM, que alberga genótipos de grande importância em nível regional. Para detectar esta família recentemente foi descrito um novo marcador genético que implica o estudo da presença do SNP Ag85C (G103A) por PCR-RFLP. Este marcador pode servir como complemento ao método de Spoligotyping, para ajudar a classificar certos spoligotypes mal definidos ou convergentes. O mesmo pode também constituir uma alternativa rápida e simples para detectar isolados LAM, o que pode ser de grande importância para países de menos recursos. Finalmente, estudamos a freqüência da linhagem RDRio ao analisarmos além da presença do LSP RDRio, outros marcadores para esta linhagem como a deleção de RD174, a presença de 2 cópias alélicas no MIRU 2 e 1 cópia no MIRU 40 e genótipo LAM. Quanto ao primeiro objetivo, foi observada uma distribuição ampla diferença nas famílias de spoligotyping, de acordo com a região estudada. A maior taxa de agrupamentos de isolados observou-se na população de isolados venezuelanos, sugerindo um processo de transmissão contínua ou introdução de algumas linhagens muito tempo atrás. Quanto ao segundo objetivo, foi observado um comportamento variado do SNP Ag85C em referência ao seu papel como marcador da família LAM. Nos isolados paraguaios e venezuelanos houve um grau significativo de concordância com os resultados de spoligotyping enquanto nos pacientes argentinos houve baixa concordância. A relação entre a presença do SNP Ag85C (G103A) e o genótipo LAM, entretanto deve ser melhor estudada. Finalmente, foi observada a presença da linhagem RDRio nos três países, com maior frequência na Venezuela. A baixa freqüência no Paraguai chama a atenção e deve ser mais bem estudada. As características genéticas da linhagem RDRio, salvo poucas exceções, estiveram presentes nos isolados dos três países estudados. Após construção de árvores tipo MST com base nos spoligotypes, observamos, como previsto, a família LAM9 como nodo central que contem a linhagem RDRio e do qual derivam os outros subtipos da família LAM. No entanto, os agrupamentos dos isolados RDRio nos MST de MIRU-VNTR devem ser melhor estudados. Concluindo, no presente estudo descrevemos os genótipos circulantes de Mtb nos três países estudados. Os resultados geraram perguntas interessantes que devem ser abordados no futuro. / Tuberculosis (TB) still remains an important cause of death and its etiological agent is Mycobacterium tuberculosis (Mtb). Globally the “Latin American-Mediterranean” (LAM) accounts for 15% of the TB cases. Within this family, recently a new Long Sequence Polymorphism has been described, characterized by a 26,3kb deletion and assigned as RDRio-LSP. Currently this lineage constitutes the predominant clade (30%) of TB cases in Rio de Janeiro, and preliminary studies suggest that strains belonging to this genotype present enhanced virulence. Consequently it is of considerable importance to evaluate the distribution and the transmissibility of this genotype regionally as well as globally. The present descriptive study had three main objectives, starting with the study of the genetic variability of the Mtb strains of three Latin-American countries, Paraguay, Argentina, and Venezuela, with little information about the presence and nature of Mtb genotypes; by using spoligotyping and MIRU-VNTR 12 loci. The other focus was on the LAM family, which regionally is of extreme importance. To detect this family recently a new genetic marker has been described, which involves the detection of the SNP Ag85C (G103A) by PCR-RFLP. This marker allows classifying ill defined spoligopatterns or convergent spoligopatterns and it constitutes a simple and fast method, thus an important alternative in low resource countries. Finally our third goal was to study the prevalence of the RDRio lineage by analyzing the RDRio-LSP as well as other markers like the deletion of RD174, 2 allelic copies in MIRU 2 and 1 in MIRU 40, as well as LAM genotype. Regarding the first objective, the population structure obtained by spoligotyping varied considerably according to the geographical location.The highest clustering rate was detected within strains from Venezuela. As for the second objective, we observed a variable behavior of the SNP Ag85C as a marker for the LAM family. Both in strains from Paraguay and Venezuela, the marker had a significant concordance when compared with spoligotyping results. Nevertheless within Argentinean strains, the opposite was observed. Thus, the relationship between the presence of the Ag85C (G103A) SNP and the LAM genotype must be further studied. As for the third objective, we observed the presence of the RDRio lineage within the three countries, but with the highest prevalence in Venezuela. The low prevalence observed within Paraguayan strains deserves more studies. The genetic characteristics of the RDRio lineage were present within the RDRio strains of the three countries, with some few exceptions. Expectedly the construction of MST trees allowed us to observe a central node of the LAM9 family that contained the RDRio lineage and of which other sub-types of the LAM family derived. The nodes observed within the MST based on MIRU, must be further analyzed. In conclusion, we managed to describe the genotypes circulating within the three studied countries. From the results obtained, interesting questions were posed that should be analyzed in the future.

Significant advancements in genome sequencing techniques and other high-throughput initiatives have resulted in the availability of complete sequences of genomes of a large number of organisms, which provide an opportunity to study detailed biological information encoded therein. Identification of functional roles of proteins can aid in comprehension of various cellular activities in an organism, which is traditionally achieved using techniques pertaining to the field of molecular biology, protein chemistry and macromolecular crystallography. The established experimental methods for protein structure and function determination, although accurate and resourceful, are laborious and time consuming. Computational analyses of sequences of gene products and exploration of evolutionary relationships can give clues on protein structure and/or function with reasonable accuracy which can be used to direct experimental studies on proteins of interest, effectively. Moreover, with growing volumes of data, there has been a growing disparity in the number of well-characterized and uncharacterized proteins, further necessitating the use of computational methods for investigating evolutionary and structure-function relationships. The remarkable progress made in the development of computational techniques (Chapter 1) has immensely contributed to the state-of-the-art biological sequence analysis and recognition of protein structure and function in a reliable manner. These methods have largely influenced the exploration of protein sequence-structure-function space. One of the relevant applications of computational approaches is in the understanding of functional make-up of human pathogens, their complex interplay with the host and implications in pathogenesis. In this thesis, sensitive profile-based search procedures have been utilized to address various aspects in the context of three pathogens- Mycobacterium tuberculosis, Plasmodium falciparum and Trypanosoma brucei, which are causative agents of potentially life- threatening diseases. The existing drugs approved for the diseases, although of immense value in controlling the disease, have several shortcomings, the most important of them being the emergence of drug resistance that render the current treatment regimens futile. Thus, the identification of practicable targets and new drugs or new combination therapies become an important necessity. Analyses on structural and functional repertoire of proteins encoded in the pathogenic genomes can provide means for rational identification of therapeutic intervention strategies. This thesis begins with the computational analyses of proteins encoded in M. tuberculosis genome. M. tuberculosis is a primary aetiological agent of tuberculosis in humans, and is o responsible for an estimated 1.5 million deaths every year. The complete genome of the pathogen was sequenced and made available more than a decade ago, which has been valuable in determination of functional roles of its gene products. Yet, functions of many M. tuberculosis proteins remain unknown. Computational prediction of protein function is an on- going process based on ever growing information made available in public databases as well as the introduction of powerful homology recognition techniques. Hence, a continuous refinement is essential to make the most of the sequence data, ensuring its accuracy and relevance. With the use of multiple sequence and structural profile-based search procedures, an enhanced structural and functional characterization of M. tuberculosis proteins, totalling to 95% of the genome was achieved (Chapter 2). Following are the key findings. o Domain definitions were obtained for a total of 3566 of 4018 proteins. Amino acid residue coverage of >70% was achieved for 2295 proteins which constitute more than half of the proteome. o Domain assignments were newly identified for 244 proteins with domain-unassigned regions. Structure prediction for these proteins corroborated all the remote homologyrelationships recognized using profile-based methods, enhancing the reliability of the predictions. o Comparison on domain compositions of proteins between M. tuberculosis and human host, revealed presence of pathogen-specific domains that are not homologous to proteins in human. Such proteins in M. tuberculosis are mainly virulence factors involved in host-pathogen interactions such as immune-dominance and aiding entry and survival in human host macrophages, hence forming attractive targets for drug discovery. o Putative structural and functional information for proteins with no recognizable domains were inferred by means of fold recognition and an iterative profile-based search against sequence database. o Attributing putative structures and functions to 955 conserved hypothetical proteins in M. tuberculosis, 137 of which are reportedly essential to the pathogen, provide a basis to re-investigate their involvement in pathogenesis and survival in the host. Proteins with no detectable homologues were recognized as M. tuberculosis H37Rv-specific, which can serve as promising drug targets. An attempt was made to identify porin-like proteins in M. tuberculosis, considering MspA porin from M. smegmatis as a template. The difficulty in recognition of putative porins in M. tuberculosis is indicative of novel outer membrane channel proteins, not characterized yet, or high representation of ion-channels, symporters and transporters to compensate for the functional role of porins. In addition, MspA-like proteins were not readily recognized in other slow-growing mycobacterial pathogens that are known to infect human host, apart from M. tuberculosis. This indicates probable acquisition of physiological adaptations, i.e. absence of porins, to confer drug-resistance, in the course of their co-evolution with human hosts. Evolutionary relationships recognized between sequence (Pfam) and structural (SCOP) families aided in association of potential structures and/or functions for 55 uncharacterized Pfam domains recognized in M. tuberculosis. Such associations deliver useful insights into the structure and function of a protein housing the uncharacterized domain. The functional inferences drawn for M. tuberculosis proteins based on the predictions can provide valuable basis for experimental endeavours in understanding mechanisms of pathogenesis and can significantly impact anti-tubercular drug discovery programmes. An interesting outcome benefitted from the exercise of exploring relationships between Pfam and SCOP families, was the identification of evolutionary relationship between a Pfam domain of unknown function DUF2652 and class III nucleotidyl cyclases. A detailed investigation was undertaken to assess this relationship (Chapter 3). Nucleotidyl cyclases synthesize cyclic nucleotides which are critical second messengers in signalling pathways. The DUF2652 family predominantly comprises of bacterial proteins belonging to three lineages- Actinobacteria, Bacteroidetes and Proteobacteria. Thus, recognition of evolutionary relationship between these bacterial proteins and nucleotide cyclases is of particular interest due to the indispensability of cyclic nucleotides in regulation of varied biological activities in bacteria. Use of fold recognition program suggested presence of nucleotide cyclase-characteristic topological motif (βααββαβ) in all the members of the DUF2652 family. Detailed analyses on structural and functional features of the uncharacterized set of bacterial proteins corresponding to 50 bacterial genomes, using profile- based alignments, revealed presence of key features typical of nucleotidyl cyclases, including metal-binding aspartates, substrate-specifying residues and transition-state stabilizing residues. Depending on the features, 20 proteins of Actinobacteria lineage, predominantly mycobacteria, of unknown structure and function were identified as putative nucleotide cyclases, 23 proteins of Bacteroidetes lineage were associated with guanylyl cyclases, while 8 uncharacterized proteins of Proteobacteria were recognized as nucleotide cyclase-like proteins (7 adenylyl and one guanylyl cyclase). Sequence similarity-based clustering of the predicted nucleotide cyclase-like proteins with established nucleotide cyclases indicated the apparent evolutionarily distinctness of the subfamily of class III nucleotidyl cyclases predicted. Furthermore, analysis of evolutionarily conserved gene clusters of the predicted nucleotide cyclase-like proteins indicated functional associations that support the predictions on their participation in cellular signalling events. The inferences made can be experimentally investigated further to ascertain the involvement of the uncharacterized bacterial proteins in signalling pathways, which can help in understanding the pathobiology of pathogenic species of interest. The next objective was the recognition of biologically relevant protein-protein interactions across M. tuberculosis and human host (Chapter 4). M. tuberculosis is well known for its ability to successfully co-evolve with human host in terms of establishing infection, survival and persistence. The current knowledge on the mechanisms of host invasion, immune evasion and persistence in the host environment can be attributed, and is limited, to the experimental studies pursued by numerous groups. Chapter 4 presents an approach for computational identification of biologically feasible protein-protein interactions across M. tuberculosis and human host. The approach utilizes crystal structures of intra-organism protein-protein complexes which are transient in nature. Identification of homologues of host and pathogen proteins in the database of known protein-protein interactions, formed the initial step, followed by identification of conserved interfacial patch and integration of information on tissue-specific expression of human proteins and subcellular localization of human and M. tuberculosis proteins. In addition, appropriate filters were used to extract biologically feasible host-pathogen protein-protein interactions. This resulted in recognition of 386 interactions potentially mediated by 59 M. tuberculosis proteins and 90 human proteins. A predominance of host-pathogen interactions (193 protein-protein interactions) brought about by M. tuberculosis proteins participating in cell wall processes, was observed, which is in concurrence with the experimental studies on immuno-modulatory activities brought about by such proteins. These set of mycobacterial proteins were predicted to interact with diverse set of host proteins such as those involved in ubiquitin conjugation pathways, metabolic pathways, signalling pathways, regulation of cell proliferation, transport, apoptosis and autophagy. The predictions have the potential to complement experimental observations at the molecular level. Details on couple of interesting cases are presented in the chapter, one of which is the probable mechanism of immune evasion adopted by M. tuberculosis to inhibit lysozyme activity in macrophages, and second is the mechanism of nutrient uptake from host. The set of M. tuberculosis proteins predicted to mediate interactions with host proteins have the potential to warrant an experimental follow-up on probable mechanisms of pathogenesis and also serve as attractive targets for chemotherapeutic interventions. proteins known to participate in P. falciparum metabolism. Pathway holes, where evidence on metabolic step exists but the catalysing enzyme is not known, have also been addressed in the study, several of which have been suggested to play an important role in growth and development of the parasite during its intra-erythrocytic stages in human host. A subsequent objective was the recognition P. falciparum proteins potentially capable of remodelling erythrocytes to suit their niche (Chapter 7). Exploitative mechanisms are brought about by the parasite to remodel erythrocytes for growth and survival during intra-erythrocytic stages of its life-cycle, the understanding of which is limited to experimental studies. To achieve physicochemically viable protein-protein interactions potentially mediated by proteins of human erythrocytes and P. falciparum proteins, a structure-influenced protocol, similar to the one demonstrated in Chapter 4, was employed. Information on subcellular localization and protein expression is crucial especially for parasites like P. falciparum, which reside in One of the major shortcomings with current treatment regimen for tuberculosis is the emergence of multidrug (MDR) and extensively drug-resistant (XDR) strains that render first-line and second-line drug treatments futile. This entails a need to explore target space in M. tuberculosis as well as explore the potential of existing drugs for repurposing against tuberculosis. A drug repurposing strategy i.e. exploring within-target-family selectivity of small molecules, has been implemented (Chapter 5) to contribute towards time and cost-saving anti-tubercular drug development efforts. With the use of profile-based search procedures, evolutionary relationships between targets (other than proteins of M. tuberculosis) of FDA-approved drugs and M. tuberculosis proteins were investigated. A key filter to exclude drugs capable of acting on human proteins substantially reduced the chances of obtaining anti-targets. Thus, total of 130 FDA-approved drugs were recognized that can be repurposed against 78 M. tuberculosis proteins, belonging to the functional categories- intermediary metabolism and respiration, information pathways, cell wall and cell processes and lipid metabolism. The catalogue of structure and function of M. tuberculosis proteins and their involvement in host-pathogen protein-protein interactions compiled from chapters 2 and 4 served as a guiding tool to explore the functional importance of targets identified. Many of the potential targets identified have been experimentally shown to be essential for growth and survival of the pathogen earlier, thus gaining importance in terms of pharmaceutical relevance. Polypharmacological drugs or drugs capable of acting of multiple targets were also identified (92 drugs) in the study. These drugs have the potential to stand tolerance against development of drug resistance in the pathogen. Comparative sequence and structure-based analysis of M. tuberculosis proteins homologous to known targets yielded credible inferences on putative binding sites of FDA-approved drugs in potential targets. Instances where information on binding sites could not be readily inferred from known targets, potentially druggable sites have been predicted. Comparison with earlier experimental studies that report anti-tubercular potential of several approved drugs enhanced the credibility of 74 of 130 FDA-approved drugs that can be readily prioritized for clinical studies. An additional exercise was pursued to identify prospective anti-tubercular agents by means of structural comparison between ChEMBL compounds and 130 FDA-approved drugs. Only those compounds were retained that showed considerably high structural similarity with approved drugs. Such compounds with minor changes in terms of physicochemical properties provide a basis for exploration of compounds that may exhibit higher affinities to bind to M. tuberculosis targets. The set of approved drugs recognized as repurpose-able candidates against tuberculosis, in concert with the structurally similar compounds, can significantly impact anti-tubercular drug development and drug discovery. The next part of the thesis focuses on Plasmodium falciparum, an obligate intracellular protozoan parasite responsible for malaria. The parasite genome features unusual characteristics including abundance of low complexity regions and pronounced sequence divergence that render protein structure and function recognition difficult. The parasite also manifests remarkable plasticity in its metabolic organization throughout its developmental stages in two hosts-human and mosquito; thus obtaining an exhaustive list of metabolic proteins in the parasite gains importance. Considering the utility of multiple sensitive profile-based search approaches in enhanced annotation of M. tuberculosis genome, a similar exercise was employed to recognize potential metabolic proteins in P. falciparum (Chapter 6). A total of 172 metabolic proteins were identified as participants of 78 metabolic pathways, over and above 609heterogeneous environmental conditions at different stages in their lifecycle. Inclusion of such data aided in extraction of 208 biologically relevant protein-protein interactions potentially mediated by 59 P. falciparum proteins and 30 erythrocyte proteins. Host-parasite protein-protein interactions were predicted pertaining to several major strategies spanning intra-erythrocytic stages in P. falciparum pathogenesis including- gaining entry into the host erythrocytes (category: RBC invasion, protease), redirecting parasitic proteins to erythrocyte membrane (category: protein traffic), modulating erythrocyte machinery (category: rosette formation, putative adhesin, chaperone, kinase), evading immunity (category: immune evasion) and eventually egress (category: merozoite egress) to infect other uninfected erythrocytes. Elaborate means to analyse and evaluate the functional viability of a predicted interaction in terms of geometrical packing at the interfacial region, electrostatic complementarity of the interacting surfaces and interaction energies is also demonstrated. The protein-protein interactions, thus predicted between human erythrocytes and P. falciparum, have the potential to provide a useful basis in understanding probable mechanisms of pathogenesis, and indeed in pinning down attractive targets for antimalarial drug discovery. The emergence of drug resistance against all known antimalarial agents, currently in use, necessitates discovery and development of either new antimalarial agents or unexplored combination of drugs that may not only reduce mortality and morbidity of malaria, but also reduce the risk of resistance to antimalarial drugs. In an attempt to contribute towards the same, Chapter 8 explores the established concept of within-target-family selectivity of small molecules to recognize antimalarial potential of the approved drugs. Eighty six FDA-approved drugs, predominantly constituted by antibacterial agents, were identified as feasible candidates for repurposing against 90 P. falciparum proteins. Most of the potential parasite targets identified are known to participate in housekeeping machinery, protein biosynthesis, metabolic pathways and cell growth and differentiation, and thus are pharmaceutically relevant. During intra-erythrocytic growth of P. falciparum, the parasite resides within the erythrocyte, within a protective encasing, known as parasitophorous vacuole. Hence a drug, intended to target a parasite protein residing in an organelle, must be sufficiently hydrophilic or hydrophobic to be able to permeate cell membranes and reach its site of activity. On the basis of lipophilicity of the drugs, a physical property determined experimentally, 57 of 86 FDA-approved drugs were recognized as feasible candidates for use against P. falciparum during the course of blood-stages of infection, which can be prioritized for antimalarial drug development programmes. The final section of the thesis focuses on the protozoan parasite Trypanosoma brucei, a causative agent of African sleeping sickness (Chapter 9). This disease is endemic to sub-Saharan regions of Africa. Despite the availability of completely sequenced genome of T. brucei, structure and function for about 50% of the proteins encoded in the genome remain unknown. Absence of prophylactic chemotherapy and vaccine, compounded with emergence of drug-resistance renders anti-trypanosomal drug discovery challenging. Thus, considering the utility of frameworks established in earlier chapters for recognition of protein structure, function and drug-targets, similar steps were undertaken to understand functional repertoire of the parasite and use drug repurposing methods to accelerate anti-trypanosomal drug discovery efforts. Structures and functions were reliably recognized for 70% of the gene products (5894) encoded in T. brucei genome, with the use of multiple profile-based search procedures, coupled with information on presence of transmembrane domains and signal peptide cleavage sites. Consequently, a total of 282 uncharacterized T. brucei proteins could be newly coined as potential metabolic proteins. Integration of information on stage-specific expression profiles with Trypanosoma-specific and T-.brucei-specific proteins identified in the study, aided in pinning down potential attractive targets. Additionally, exploration of evolutionary relationships between targets of FDA-approved drugs and T. brucei proteins, 68 FDA-approved drugs were predicted as repurpose-able candidates against 42 potential T. brucei targets which primarily include proteins involved in regulatory processes and metabolism. Several targets predicted are reportedly essential in assisting the parasite to switch between differentiation forms (bloodstream and procyclic) in the course of its lifecycle. These targets are of high therapeutic relevance, hence the corresponding drug-target associations provide a useful resource for experimental endeavours. In summary, this thesis presents computational analyses on three pathogenic genomes in terms of enhancing the understanding of functional repertoire of the pathogens, addressing metabolic pathway holes, exploring probable mechanisms of pathogenesis brought about by potential host-pathogen protein-protein interactions, and identifying feasible FDA-approved drug candidates to repurpose against the pathogens. The studies are pursued primarily by taking advantage of powerful homology-detection techniques and the ever-growing biological information made available in public databases. Indeed, the inferences drawn for the three pathogenic genomes serve an excellent resource for an experimental follow-up. The set of protocols presented in the thesis are highly generic in nature, as demonstrated for three pathogens, and can be utilized for genome-wide analyses on many other pathogens of interest. The supplemental data associated with the chapters is provided in a compact disc attached with this thesis.

Made available in DSpace on 2014-07-29T15:30:42Z (GMT). No. of bitstreams: 1 DISSSERTACAO CRISTINA.pdf: 3986017 bytes, checksum: add734135c93af474b22c106a3a51d41 (MD5) Previous issue date: 2007-02-23 / Tuberculosis (TB) caused by the Mycobacterium tuberculosis is responsible for more than 2 million deaths annually in the world. Although one third of the world population carries the bacilli, only 5% of the infected people develop active disease. 80% of TB cases are concentrated among 21 countries and Brazil is one of them. Due to the flaws existent inherent to the TB diagnostic, many works try to discover TB antigens to be used on an ELISA assay to detect active TB in underdeveloped countries, endemic for other mycobacterias such as leprosy (L). Several researches are identifying bacilli proteins obtained under nutritional culture stress, in order to mimic the intracellular milieu. It has been selected the immunodominant ones as a disease marker. The objective of this work was to characterize the humoral immune response among TB patients to rMPT-51 and rGlcB using an ELISA indirect test. It was included voluntaries that were HIV negative and the ones without pregnancy, chronic disease or under immunosuppressant treatment. Serum IgM and IgG against MPT-51 and GlcB recombinant antigen from 49 patients with active tuberculosis were measured by indirect ELISA, paired by sex and age with: healthy PPD negative individuals (controls) and lepromatous leprosy patients (LL). Patients with TB (0.810±0.319) showed higher levels of IgM against rMPT-51 than both LL individuals (0.454±0.195) and control (0.448±0.162) with statistical significant, p= 0.001 and p<0.001) respectively. These test showed 96.9% specificity and 67.3% sensitivity. Conversely, tuberculosis, controls, and LL individuals showed lower of MPT-51 IgG levels, which could not be distinguished among the groups. rGlcB antigen was able to distinguished TB patients from controls for IgM levels (specificity, 95.9 % and sensitivity, 8.2% ) and IgG levels (specificity, 99% and sensitivity, 18.2%). In order to evaluate the profiles of IgM and IgG against rMPT-51 and rGlcB before and after chemotherapy, the sera from 11 patients was collected and paired according to the treatment status (before or after). IgM and IgG against MPT-51 remained with the same profile levels before and after the treatment. The levels of serum IgG against rGlcB clearly diminished after the chemotherapy (p<0.01). Our results suggest that serum IgM levels against recombinant MPT-51 could to be useful to distinguish between active TB, controls and LL individuals. In addition, after TB treatment the IgG response to rGlcB diminished suggesting that it could be used to follow up of the TB treatment / A tuberculose (TB) causada pelo Mycobacterium tuberculosis, resulta em mais de dois milhões de óbitos anualmente. Segundo dados da OMS cerca de 30% da população mundial está infectada com o bacilo e 5% desses infectados desenvolvem a doença ativa. O Brasil juntamente com outros 21 países albergam 80% de todos os casos de TB no mundo. Devido às falhas existentes nos diagnósticos atuais, muitos estudos tentam descobrir antígenos do Mycobacterium tuberculosis que podem ser usados no ensaio ELISA, um teste de baixo custo. Esse teste seria de grande valia na identificação da tuberculose, principalmente em países em desenvolvimento endêmicos também para outras micobactérias como M. leprae, causador da hanseníase. Várias pesquisas têm identificado componentes do bacilo, como as proteínas secretadas pelas micobactérias em cultura, sob condições de estresse nutricional, mimetizando a situação vivida pelo bacilo no meio intracelular. Dessas proteínas são selecionadas aquelas imunodominantes que podem ser usadas como marcadores da doença. O objetivo desse trabalho foi avaliar o reconhecimento de duas dessas proteínas: os antígenos recombinantes MPT-51 e GlcB do M. tuberculosis, por anticorpos séricos da classe IgM e IgG de pacientes com tuberculose ativa, pelo método imunoenzimático indireto (ELISA). Foram adotados como critério de inclusão aqueles indivíduos de qualquer grupo que fossem HIV negativos, sem doenças crônicas ou uso de medicamentos imunossupressores e mulheres não gestantes. Quarenta e nove pacientes com tuberculose ativa foram selecionados e comparados com os grupos: controles saudáveis PPD não reatores e pacientes hansenianos portadores da forma Virchoviana, pareados por sexo e idade. Os pacientes com TB (0,810±0,319) mostraram maiores concentrações de IgM anti-MPT-51 que os seus respectivos controles: pacientes com hanseníase 0,454±0,195) e controles saudáveis (0,448±0,162), com diferença estatística, p=0,001 e p<0,001 respectivamente. Os ensaios de ELISA nas dosagens de IgM e IgG anti-MPT-51 mostraram especificidade de 96,9% e 98,0% e sensibilidade 67,3% e 4,1% respectivamente. Para o antígeno GlcB, os ensaios de ELISA na dosagem de IgM e IgG mostraram especificidade de 95,9% e 99% e sensibilidade 8,2% e 18,2%. Onze pacientes com TB foram monitorados durante o tratamento, com realização de dosagens dos níveis de anticorpos IgM e IgG específicos ao rMPT-51 e rGlcB antes e após a terapia. Os níveis de anticorpos IgM e IgG anti-MPT-51 antes e após a terapia não sofreram alterações significativas. Entretanto, os níveis de IgG anti-GlcB diminuíram após a terapia (p<0,01). Nossos resultados sugerem que o rMPT-51 pode ser usado como marcador da tuberculose, quando mensuradas as concentrações de IgM específicos, pois foram capazes de discriminar pacientes TB, de controles e hansenianos. Apesar do ELISA ter demonstrado baixa sensibilidade quando o antígeno rGlcB foi utilizado, os níveis séricos dos anticorpos da classe IgG diminuíram após o tratamento da tuberculose (p<0,01) sugerindo que essa técnica poderia ser utilizada para o acompanhamento da terapêutica

The eubacterial genome is maintained in a negatively supercoiled state which facilitates its compaction and storage in a small cellular space. Genome supercoiling can potentially influence various DNA transaction processes such as DNA replication, transcription, recombination, chromosome segregation and gene expression. Alterations in the genome supercoiling have global impact on the gene expression and cell growth. Inside the cell, the genome supercoiling is maintained judiciously by DNA topoisomerases to optimize DNA transaction processes. These enzymes solve the problems associated with the DNA topology by cutting and rejoining the DNA. Due to their essential cellular functions and global regulatory roles, DNA topoisomerases are fascinating candidates for the study of the effect of topology perturbation on a global scale. Genus Mycobacterium includes a large number of species including the well-studied Mycobacterium smegmatis (Msm) as well as various pathogens–Mycobacterium leprae, Mycobacterium abscessus and Mycobacterium tuberculosis (Mtb), the last one being the causative agent of the deadly disease Tuberculosis (TB), which claims millions of lives worldwide annually. The organism combats various stresses and alterations in its environment during the pathogenesis and virulence. During such adaptation, various metabolic pathways and transcriptional networks are reconfigured. Considering their global regulatory role, DNA topoisomerases and genome supercoiling may have an influence on the mycobacterial survival and adaptation. Biochemical studies from our laboratory have revealed several distinctive characteristics of mycobacterial DNA gyrase and topoisomerase I. DNA gyrase has been shown to be a strong decatenase apart from its characteristic supercoiling activity. Similarly, the mycobacterial topoisomerase I exhibits several distinct features such as the ability to bind both single- as well as double-stranded DNA, site specific DNA binding and absence of Zn2+ fingers required for DNA relaxation activity in other Type I enzymes. Although, efforts have been made to understand the biochemistry and mechanism of mycobacterial topoisomerases, in vivo significance and regulatory roles remain to be explored. The present study is aimed at understanding the mechanism, in vivo functions, regulation and genome wide distribution of mycobacterial topoisomerases. Chapter 1 of the thesis provides introduction on DNA topology, genome supercoiling and DNA topoisomerases. The importance of genome supercoiling and its regulatory roles has been discussed. Further, the regulation of topoisomerase activity and the role in the virulence gene regulation is described. Finally, a brief overview of Mtb genome, disease epidemiology, and pathogenesis is presented along with the description of the work on mycobacterial topoisomerases. In Chapter 2, the studies are directed to understand the DNA relaxation mechanism of mycobacterial Type IA topoisomerase which lack Zn2+ fingers. The N-terminal domain (NTD) of the Type IA topoisomerases harbor DNA cleavage and religation activities, but the carboxyl terminal domain (CTD) is highly diverse. Most of these enzymes contain a varied number of Zn2+ finger motifs in the CTD. The Zn2+ finger motifs were found to be essential in Escherichia coli TopoI but dispensable in the Thermotoga maritima enzyme. Although, the CTD of mycobacterial TopoI lacks Zn2+ fingers, it is indispensable for the DNA relaxation activity of the enzyme. The divergent CTD harbors three stretches of basic amino acids needed for the strand passage step of the reaction as demonstrated by a new assay. It is elucidated that the basic amino acids constitute an independent DNA-binding site apart from the NTD and assist the simultaneous binding of two molecules of DNA to the enzyme, as required during the strand passage step of the catalysis. It is hypothesized that the loss of Zn2+ fingers from the mycobacterial TopoI could be associated with Zn2+ export and homeostasis. In Chapter 3, the studies have been carried out to understand the regulation of mycobacterial TopoI. Identification of Transcription Start Site (TSS) suggested the presence of multiple promoters which were found to be sensitive to genome supercoiling. The promoter activity was found to be specific to mycobacteria as the promoter(s) did not show activity in E. coli. Analysis of the putative promoter elements suggested the non-optimal spacing of the putative -35 and -10 promoter elements indicating the involvement of supercoiling for the optimal alignment during the transcription. Moreover, upon genome relaxation, the occupancy of RNA polymerase was decreased on the promoter region of topoI gene implicating the role of DNA topology in the Supercoiling Sensitive Transcription (SST) of TopoI gene from mycobacteria. The involvement of intrinsic promoter elements in such regulation has been proposed. In Chapter 4, the importance of TopoI for the Mtb growth and survival has been validated. Mtb contains only one Type IA topoisomerase (Rv3646c), a sole DNA relaxase in the cell, and hence a candidate drug target. To validate the essentiality of Mtb topoisomerase I for bacterial growth and survival, conditionally regulated strain of topoI in Mtb was generated. The conditional knockdown mutant exhibited delayed growth on agar plate and in liquid culture the growth was drastically impaired when TopoI expression was suppressed. Additionally, novobiocin and isoniazid showed enhanced inhibitory potential against the conditional mutant. Analysis of the nucleoid revealed its altered architecture upon TopoI depletion. These studies establish the essentiality of TopoI for the Mtb growth and open up new avenues for targeting the enzyme. In Chapter 5, the influence of perturbation of TopoI activity on the Msm growth and physiology has been studied. Notably, Msm contains an additional DNA relaxation enzyme– an atypical Type II topoisomerase TopoNM. The TopoI depleted strain exhibited slow growth and drastic change in phenotypic characters. Moreover, the genome architecture was disturbed upon depletion of TopoI. Further, the proteomic and transcript analysis indicated the altered expression of the genes involved in central metabolic pathways and core DNA transaction processes in the mutant. The study suggests the importance of TopoI in the maintenance of cellular phenotype and growth characteristics of fast growing mycobacteria having additional topoisomerases. In Chapter 6, the ChIP-Seq method is used to decipher the genome wide distribution of the DNA gyrase, topoisomerase I (TopoI) and RNA polymerase (RNAP). Analysis of the ChIP-Seq data revealed the genome wide distribution of topoisomerases along with RNAP. Importantly, the signals of topoisomerases and RNAP was found to be co-localized on the genome suggesting their functional association in the twin supercoiled domain model, originally proposed by J. C. Wang. Closer inspection of the occupancy profile of topoisomerases and RNAP on transcription units (TUs) revealed their co-existence validating the topoisomerases occupancy within the twin supercoiled domains. On the genomic scale, the distribution of topoisomerases was found to be more at the ori domains compared to the ter domain which appeared to be an attribute of higher torsional stress at ori. The reappearance of gyrase binding at the ter domain (and the lack of it in the ter domain of E. coli) suggests a role for Mtb gyrase in the decatenation of the daughter chromosomes at the end of replication. The eubacterial genome is maintained in a negatively supercoiled state which facilitates its compaction and storage in a small cellular space. Genome supercoiling can potentially influence various DNA transaction processes such as DNA replication, transcription, recombination, chromosome segregation and gene expression. Alterations in the genome supercoiling have global impact on the gene expression and cell growth. Inside the cell, the genome supercoiling is maintained judiciously by DNA topoisomerases to optimize DNA transaction processes. These enzymes solve the problems associated with the DNA topology by cutting and rejoining the DNA. Due to their essential cellular functions and global regulatory roles, DNA topoisomerases are fascinating candidates for the study of the effect of topology perturbation on a global scale. Genus Mycobacterium includes a large number of species including the well-studied Mycobacterium smegmatis (Msm) as well as various pathogens–Mycobacterium leprae, Mycobacterium abscessus and Mycobacterium tuberculosis (Mtb), the last one being the causative agent of the deadly disease Tuberculosis (TB), which claims millions of lives worldwide annually. The organism combats various stresses and alterations in its environment during the pathogenesis and virulence. During such adaptation, various metabolic pathways and transcriptional networks are reconfigured. Considering their global regulatory role, DNA topoisomerases and genome supercoiling may have an influence on the mycobacterial survival and adaptation. Biochemical studies from our laboratory have revealed several distinctive characteristics of mycobacterial DNA gyrase and topoisomerase I. DNA gyrase has been shown to be a strong decatenase apart from its characteristic supercoiling activity. Similarly, the mycobacterial topoisomerase I exhibits several distinct features such as the ability to bind both single- as well as double-stranded DNA, site specific DNA binding and absence of Zn2+ fingers required for DNA relaxation activity in other Type I enzymes. Although, efforts have been made to understand the biochemistry and mechanism of mycobacterial topoisomerases, in vivo significance and regulatory roles remain to be explored. The present study is aimed at understanding the mechanism, in vivo functions, regulation and genome wide distribution of mycobacterial topoisomerases. Chapter 1 of the thesis provides introduction on DNA topology, genome supercoiling and DNA topoisomerases. The importance of genome supercoiling and its regulatory roles has been discussed. Further, the regulation of topoisomerase activity and the role in the virulence gene regulation is described. Finally, a brief overview of Mtb genome, disease epidemiology, and pathogenesis is presented along with the description of the work on mycobacterial topoisomerases. In Chapter 2, the studies are directed to understand the DNA relaxation mechanism of mycobacterial Type IA topoisomerase which lack Zn2+ fingers. The N-terminal domain (NTD) of the Type IA topoisomerases harbor DNA cleavage and religation activities, but the carboxyl terminal domain (CTD) is highly diverse. Most of these enzymes contain a varied number of Zn2+ finger motifs in the CTD. The Zn2+ finger motifs were found to be essential in Escherichia coli TopoI but dispensable in the Thermotoga maritima enzyme. Although, the CTD of mycobacterial TopoI lacks Zn2+ fingers, it is indispensable for the DNA relaxation activity of the enzyme. The divergent CTD harbors three stretches of basic amino acids needed for the strand passage step of the reaction as demonstrated by a new assay. It is elucidated that the basic amino acids constitute an independent DNA-binding site apart from the NTD and assist the simultaneous binding of two molecules of DNA to the enzyme, as required during the strand passage step of the catalysis. It is hypothesized that the loss of Zn2+ fingers from the mycobacterial TopoI could be associated with Zn2+ export and homeostasis. In Chapter 3, the studies have been carried out to understand the regulation of mycobacterial TopoI. Identification of Transcription Start Site (TSS) suggested the presence of multiple promoters which were found to be sensitive to genome supercoiling. The promoter activity was found to be specific to mycobacteria as the promoter(s) did not show activity in E. coli. Analysis of the putative promoter elements suggested the non-optimal spacing of the putative -35 and -10 promoter elements indicating the involvement of supercoiling for the optimal alignment during the transcription. Moreover, upon genome relaxation, the occupancy of RNA polymerase was decreased on the promoter region of topoI gene implicating the role of DNA topology in the Supercoiling Sensitive Transcription (SST) of TopoI gene from mycobacteria. The involvement of intrinsic promoter elements in such regulation has been proposed. In Chapter 4, the importance of TopoI for the Mtb growth and survival has been validated. Mtb contains only one Type IA topoisomerase (Rv3646c), a sole DNA relaxase in the cell, and hence a candidate drug target. To validate the essentiality of Mtb topoisomerase I for bacterial growth and survival, conditionally regulated strain of topoI in Mtb was generated. The conditional knockdown mutant exhibited delayed growth on agar plate and in liquid culture the growth was drastically impaired when TopoI expression was suppressed. Additionally, novobiocin and isoniazid showed enhanced inhibitory potential against the conditional mutant. Analysis of the nucleoid revealed its altered architecture upon TopoI depletion. These studies establish the essentiality of TopoI for the Mtb growth and open up new avenues for targeting the enzyme. In Chapter 5, the influence of perturbation of TopoI activity on the Msm growth and physiology has been studied. Notably, Msm contains an additional DNA relaxation enzyme– an atypical Type II topoisomerase TopoNM. The TopoI depleted strain exhibited slow growth and drastic change in phenotypic characters. Moreover, the genome architecture was disturbed upon depletion of TopoI. Further, the proteomic and transcript analysis indicated the altered expression of the genes involved in central metabolic pathways and core DNA transaction processes in the mutant. The study suggests the importance of TopoI in the maintenance of cellular phenotype and growth characteristics of fast growing mycobacteria having additional topoisomerases. In Chapter 6, the ChIP-Seq method is used to decipher the genome wide distribution of the DNA gyrase, topoisomerase I (TopoI) and RNA polymerase (RNAP). Analysis of the ChIP-Seq data revealed the genome wide distribution of topoisomerases along with RNAP. Importantly, the signals of topoisomerases and RNAP was found to be co-localized on the genome suggesting their functional association in the twin supercoiled domain model, originally proposed by J. C. Wang. Closer inspection of the occupancy profile of topoisomerases and RNAP on transcription units (TUs) revealed their co-existence validating the topoisomerases occupancy within the twin supercoiled domains. On the genomic scale, the distribution of topoisomerases was found to be more at the ori domains compared to the ter domain which appeared to be an attribute of higher torsional stress at ori. The reappearance of gyrase binding at the ter domain (and the lack of it in the ter domain of E. coli) suggests a role for Mtb gyrase in the decatenation of the daughter chromosomes at the end of replication.

The three-dimensional structures of biological macromolecules and molecular assemblies are becoming increasingly important with the changing methodologies of drug discovery. The structures aid in understanding of protein function at the molecular level: be it a macromolecular assembly, a cytosolic enzyme or an intermembrane receptor molecule. X-ray crystallography is the most powerful technique to obtain the three-dimensional structures of such molecules at or near atomic resolution. With such a wide-spread importance, crystallography is an integral part of structural biology and also of the current drug discovery programs. The present thesis mainly deals with application of the crystallographic techniques for understanding the structure and function of adenylosuccinate lyase (ASL) from bacterial pathogens Salmonella typhimurium and Mycobacterium tuberculosis as well as its non-pathogenic counterpart Mycobacterium smegmatis. Studies were also carried out to understand the structure-function relationship of the protease in the plant virus Sesbania Mosaic Virus (SeMV). The thesis has been divided into six chapters. The first chapter contains an introduction to nucleotide synthesis and ASL superfamily of enzymes known as the aspartase/fumarase superfamily based on the published literature. Chapter 2 provides the details of the techniques used for the investigations presented in this thesis. Chapters 3-5 deal with the structural and functional studies carried out on ASL from the three bacterial organisms. Chapter 6 deals with the simulation studies carried out on SeMV protease. Mechanism and importance of nucleotide synthesis is introduced in Chapter 1, with special emphasis on purine de novo and salvage pathways. ASL is introduced as an important enzyme for purine synthesis. Its superfamily, the aspartase/fumarase superfamily of enzymes is described in detail with respect to its structure, function and pathophysiology. Objectives of the present study are outlined towards the end of the chapter. The experimental and computational techniques utilized during the course of my research are described in Chapter 2. These techniques include gene cloning, protein expression and purification, kinetic and biophysical characterization of proteins, crystallization, X-ray diffraction, data collection and processing, structure solution, refinement, model building, validation and structural analysis, phylogenetic studies, molecular docking and molecular dynamic simulation studies. Adenylosuccinate lyase is an important enzyme participating in purine biosynthesis. With the emergence of drug resistant variants of various pathogens, ASL has been recognized as a drug target against microbial infections. Chapter 3 deals with the structural and functional characterization of ASL from Salmonella typhimurium. Two constructs of the StASL gene were cloned and expressed leading to the purification of truncated (residues 1-366) and full-length (residues 1-456) polypeptides. Crystallization of the two polypeptides resulted in three independent structures. The full-length structure was very similar to the E. coli ASL structure consistent with 95% amino acid sequence identity between the two polypeptides. However, the truncated structures showed large distortions, especially of the active site residues, accounting for the catalytic inactivity of the truncated polypeptide in spite of retaining all residues considered important for function. The full-length ASL was catalytically active. A unique feature observed in StASL, not reported in other ASLs, was its allosteric regulation by the substrate. Kinetic studies also revealed hysteretic behavior of the enzyme. The electron density map of the full-length structure showed two novel densities on the molecular 2-fold axis into each of which a molecule of cadavarine could be fitted. Docking studies revealed a ligand-binding site at the inter-subunit interface between the two observed densities which might represent a potential allosteric site. Combining the structural and kinetic results, a possible morpheein model of allosteric regulation of StASL was hypothesized. Chapter 4 deals with the crystallographic and kinetic investigations on ASL from Mycobacterium smegmatis and Mycobacterium tuberculosis. MsASL and MtbASL were cloned, purified and crystallized. The X-ray crystal structure of MsASL was determined at 2.16 Å resolution. It is the first report of an apo-ASL structure with a partially ordered active site C3 loop. Diffracting crystals of MtbASL could not be obtained and a model for its structure was derived using MsASL as a template. Most of the active site residues were found to be conserved with the exception of Ser 148 and Gly 319 of MsASL. Ser 148 is structurally equivalent to a threonine in most other ASLs. Gly 319 is replaced by an arginine residue in most ASLs. The two enzymes were catalytically much less active when compared to ASLs from other organisms. Arg319Gly substitution and reduced flexibility of the C3 loop might account for the low catalytic activity of mycobacterial ASLs. The low activity is consistent with the slow growth rate of Mycobacteria, their high GC containing genomes as well as with their dependence on other salvage pathways for the supply of purine nucleotides. Chapter 5 deals with the identification of the catalytic residues important for ASL catalysis in view of the earlier conflicting reports on the identity of these residues. pH-dependent kinetic studies were performed on full-length StASL. The theory behind these studies is also described in this chapter. Two residues with pKa values of 6.6 and 7.7 were identified as essential for the enzymatic activity. These results were interpreted along with structural comparison of MsASL and other superfamily enzymes with ordered C3 loops. They suggest that His 149 and either Lys 285 or Ser 279 of MsASL are the residues most likely to function as the catalytic acid and base, respectively. The final Chapter 6 of the thesis deals with the structural and dynamic studies carried out on Sesbania mosaic virus (SeMV) protease. The chapter begins with a general introduction to viruses, followed by a brief summary of SeMV. The goal of this study is to understand the interactions between the protease and VPg at a structural level using the information available from biochemical studies. Crystallographic studies initiated for the mutant H275APro and Y315APro were unsuccessful due to the insolubility of the proteins. Co-crystallization or soaking experiments of wild type protease with cognate peptides were unsuccessful due to the inability of the enzyme to bind to its substrates in the absence of VPg. Higher resolution structure of wild type protease did not yield any new insights when compared to the earlier reported structure determined at a lower resolution. In the absence of structural insights, molecular dynamic simulations were carried out on wild type protease structure and in silico generated mutants using GROMACS package. The studies showed the importance of flipping of residue Phe 301 and opening-closing of the loop region corresponding to residues 301-308 for the catalytic mechanism. The thesis concludes with Future perspectives of the various studies carried out on ASL and SeMV protease. The atomic coordinates determined from the work presented in this thesis have been deposited in the PDB and the assigned PDB codes are reported in the respective chapters. Publications cited in the thesis are listed in the Bibliography section.

This thesis is concerned with the first ever X-ray crystallographic and complimentary solution studies on mycobacterial lectins. Lectins, described as multivalent carbohydrate binding proteins of non-immune origin, are found in all kingdoms of life. As explained in the introductory chapter, those from plants and animals are the best characterized in terms of structure and function. Although not that extensive, important studies have been carried out on viral, fungal and parasite lectins as well. Bacterial lectins studied so far can be classified in to fimbrial, surface and secretory (or toxic). Applications of lectins include blood typing, cell separation and purification of glycoconjugates, mitogenic stimulation of lymphocytes, mapping of neuronal pathways and drug targeting and delivery. The work reported in the thesis lies at the intersection of two major long range programs in this laboratory, one on lectins and the other on mycobacterial proteins. Three putative lectins Rv1419 and Rv2813 from M. tuberculosis and MSMEG_3662 from M. smegmatis were chosen for exploratory studies on the basis of preliminary genomic searches. Exploratory studies on Rv1419, Rv2813 and MSMEG_3662 are described in the second chapter. MSMEG_3662 contains two domains, a LysM domain and a lectin domain (MSL) connected by a long polypeptide chain. The two M. tuberculosis proteins, full length MSMEG_3662 and MSL were cloned, expressed, purified and characterized. Rv2813 did not show any appreciable agglutination activity. It showed ATPase activity. Clearly the protein was not a lectin. Rv1419, full length MSMEG_3662 and MSL exhibited lectin characteristics. Among them, Rv1419 and MSL could be crystallized. Preliminary X-ray diffraction studies on them were carried out. Rv1419 could be successfully expressed only once. However, that was enough for the determination of crystal structure and the glycan array analysis of the lectin (Chapter 3). The monomeric lectin has a β-trefoil fold. It has high affinity for LacNAc and its Neu5Ac derivatives. Modeling studies using complexes of homologous structures, led to the identification of two carbohydrate binding sites on the lectins. Sequence comparisons of Rv1419 with homologous proteins with known structures and phylogenetic analysis involving them provide interesting insights into the relationship among trefoil lectins from different sources. X-ray crystal structure analysis of MSL and its complexes with mannose and methyl-α-mannose, the first comprehensive effort of its kind on a mycobacterial lectin, reveals a structure very similar to β-prism II fold lectins from plant sources, but with extensive unprecedented domain swapping in dimer formation (Chapter 4). The two subunits in a dimer often show small differences in structure, but the two domains, not always related by 2-fold symmetry, have the same structure. Each domain carries three sugar-binding sites, similar to those in plant lectins, one on each Greek key motif. The occurrence of β-prism II fold lectins in bacteria, with characteristics similar to those from plants, indicates that this family of lectins is of ancient origin and had evolved into a mature system before bacteria and plants diverged. In plants, the number of binding sites per domain varies between one and three, whereas the number is two in the recently reported lectin domains from Pseudomonas putida and Pseudomonas aeruginosa. An analysis of the sequences of the lectins and the lectin domains shows that the level of sequence similarity among the three Greek keys in each domain has a correlation with the number of binding sites in it. Furthermore, sequence conservation among the lectins from different species is the highest for that Greek key which carries a binding site in all of them. Thus, it would appear that carbohydrate binding influences the course of the evolution of the lectin. LysM domains have been recognized in bacteria and eukaryotes as carbohydrate-binding protein modules, but the mechanism of their binding to chitooligosaccharides is underexplored. Binding of a full length MSMEG_3662 containing LysM and lectin (MSL) domains to chitooligosaccharides has been studied using isothermal titration calorimetry and fluorescence titration (Chapter 5). This investigation demonstrates the presence of two binding sites of non-identical affinities per dimeric MSL-LysM molecule. Affinity of the molecule for chitooligosaccharides correlates with the length of the carbohydrate chain. Its binding to chitooligosaccharides is characterized by negative cooperativity in the interactions of the two domains. Apparently, the flexibility of the long linker that connects the LysM and MSL domains plays a facilitating role in this recognition. The LysM domain in MSL-LysM, like other bacterial domains but unlike plant LysM domains, recognizes equally well peptidoglycan fragments as well as chitin polymers. Interestingly, in the present case two LysM domains are enough for binding to peptidoglycan in contrast to the three reportedly required by the LysM domains of Bacillus subtilis and Lactococcus lactis. Also, the affinity of MSL-LysM for chitooligosaccharides is higher than that of LysM-chitooligosaccharide interactions reported so far. A part of the work presented in this thesis has been reported in the following publications: • Patra D, Mishra P, Surolia A, Vijayan M. 2014. Structure, interactions and evolutionary implications of a domain-swapped lectin dimer from Mycobacterium smegmatis. Glycobiology, 24:956-965. • Patra D, Sharma A, Chandran D, Vijayan M. 2011. Cloning, expression, purification, crystallization and preliminary X-ray studies of the mannose-binding lectin domain of MSMEG_3662 from Mycobacterium smegmatis. Acta Crystallogr Sect F Struct Biol Cryst Commun, 67:596-599. • Patra D, Srikalaivani R, Misra A, Singh DD, Selvaraj M, Vijayan M. 2010. Cloning, expression, purification, crystallization and preliminary X-ray studies of a secreted lectin (Rv1419) from Mycobacterium tuberculosis. Acta Crystallogr Sect F Struct Biol Cryst Commun, 66:1662-1665.

Mycobacterium tuberculosis (M.tb), the causative agent of tuberculosis (TB), has remained the largest killer among infectious diseases for over a century. The increasing emergence of drug resistant varieties such as the multidrug resistant (MDR) and extremely drug resistant (XDR) strains are only increasing the global burden of the disease. Available statistics indicate that nearly one-third of the world’s population is infected, where the bacteria remains in the latent state but can reactivate into an actively growing stage to cause disease when the individual is immunocompromised. It is thus immensely important to rethink newer strategies for containing and combating the spread of this disease. Extraction of iron from the host cell is one of the many factors that enable the bacterium to survive in the harsh environments of the host macrophages and promote tuberculosis. Host–pathogen interactions can be interpreted as the battle of two systems, each aiming to overcome the other. From the host’s perspective, iron is essential for diverse processes such as oxygen transport, repression, detoxification and DNA synthesis. Infact, during infection, both the host and the pathogen are known to fight for the available iron, thereby influencing the outcome of the infection. It is of no surprise therefore, that many studies have investigated several components of the iron regulatory machinery of M.tb and the host. However, very few attempts have been made to study the interactions between these components and how such interactions lead to a better adapted phenotype. Such studies require exploration at multiple levels of structural and functional complexity, thereby necessitating the use of a multiscale approach. Systems biology adopts an integrated approach to study and understand the function of biological systems. It involves building large scale models based on individual biochemical interactions, followed by model validation and predictions of the system’s response to perturbations, such as a gene knock-out or exposure to drug. In multiscale modeling, an approach employed in this thesis, a particular biological phenomenon is studied at different spatiotemporal levels. Studying responses at multiple scales provides a broader picture of the communications that occur between a host and pathogen. Moreover, such an analysis also provides valuable insights into how perturbation at a particular level can elicit responses at another level and help in the identification of crucial inter-level communications that can possibly be hindered or activated for a desired physiological outcome. The broad objectives of this thesis was to obtain a comprehensive in silico understanding of mycobacterial iron homeostasis and metabolism, the influence of iron on host-pathogen interactions, identification of key players that mediate such interactions, determination of the molecular consequences of inhibiting the key players and finally the global response of M.tb to altered iron concentration. Perturbation of iron homeostasis holds a strong therapeutic potential, given its essentiality in both the host and the pathogen. Understanding the workings of iron metabolism and regulation in M.tb has been a main objective, so as to ultimately obtain insights about specific therapeutic strategies that capitalize on the criticality of iron concentration. An in-depth study of iron metabolism and regulation is performed at different levels of temporal and spatial scales using diverse methods, each appropriate to investigate biological events associated with the different scales. The specific investigations carried out in the thesis are as follows, a) Reconstruction of a host-pathogen interaction (HPI) model, with focus on iron homeostasis. This study represented the inter-cellular level analysis and was crucial for the identification of key players that mediate communication between the host and pathogen. Additionally, the model also provided a mathematical framework to study the effect of perturbations and gene knock-outs. b) Understanding the influence of iron on IdeR, an iron-responsive transcription factor, also identified as a key player in the HPI model. The study was carried out at the molecular level to identify atomistic details of how IdeR senses iron and the resulting structural modifications, which finally enables IdeR-DNA interaction. The study enabled identification of residues for the functioning of IdeR. c) Genome scale identification of genes that are regulated by IdeR to obtain an overview of the various biological processes affected by changing iron concentrations and IdeR mutation in M.tb. d) To understand the direct and indirect influences of iron and IdeR on the M.tb proteome using large scale protein-protein interaction network. The study enabled identification of highest differentially regulated genes and altered activity of the different biological processes under differing iron concentrations and regulation. e) Systems level analysis of the M.tb metabolome to investigate the metabolic re-adjustments undertaken by M.tb to adapt to altered iron concentration and regulation. The conceptual details and the background of each of the methods used to study the specific aims are provided in the Methodology chapter (Chapter 2). Construction of the host-pathogen interaction (HPI) model and the insights obtained from this study are presented in Chapter 3. A rule based HPI model was built with a focus on the iron regulatory mechanisms in both the host and pathogen. The model consisted of 194 rules, of which 4 rules represented interactions between the host and pathogen. The model not only represented an overview of iron metabolism but also allowed prediction of critical interaction that had the potential to form bottleneck in the system so as to control bacterial proliferation. Infact, model simulation led to the identification of 5 bottlenecks or chokepoints in the system, which if perturbed, could successfully interfere with the host-pathogen dynamics in favour of the host. The model also provided a framework to test perturbation strategies based on the bottlenecks. The study also established the importance of an iron responsive transcription factor, IdeR for regulating iron concentration in the pathogen and mediating host-pathogen interactions. Additionally, the importance of mycobactin and transferrin as key molecular players, involved in host-pathogen dynamics was also determined. The model provided a mathematical framework to test TB pathogenesis and provided significant insights about key molecular players and perturbation strategies that can be used to enhance therapeutic strategies. Given the importance of IdeR in HPI, its molecular mechanism of activation and dimerization was explored in Chapter 4. The main objective of the study was to explore the structural details of IdeR and its iron sensing capacity at the molecular level. A combination of molecular dynamics and protein structure network (PSN) were used to analyse IdeR monomers and dimers in the presence and absence of iron. PSNs used in this thesis are based on non-covalent interactions between sidechain atoms and are quite efficient in identifying iron induced subtle conformational variations. The study distinctly indicated the role of iron in IdeR stability. Further, it was observed that IdeR monomers can take up two major conformations, the ‘open’ and ‘close’ conformation with the iron bound structure preferring the ‘close’ conformation. Major structural changes, such as the N-terminal folding and increased propensity for dimerization were observed upon iron binding. Interestingly, careful analysis of structure suggests a role of these structural modifications towards DNA binding and has been tested in the next chapter. Overall, the results clearly highlight the influence of iron on IdeR activation and dimerization. The predisposition of IdeR to bind to DNA in the presence of metal is clearly visible even when the simulations are performed solely on protein molecules. However, to confirm the conjectures proposed in this chapter and to obtain the atomistic details of IdeR-DNA interactions, the IdeR-DNA complex was investigated. Chapter 5 focuses on the mechanistic details of IdeR-DNA interactions and the influence of iron on the same. IdeR is known to bind to a specific stretch of DNA, known as the ‘iron-box’ motif to form a dimer-of-dimer complex. Molecular dynamics followed by protein-DNA bipartite network analysis was performed on a set of four IdeR-DNA complexes to obtain a molecular level understanding of IdeR-DNA interactions. A striking observation was the dissociation of IdeR-DNA complex in the absence of iron, undoubtedly establishing the importance of iron for IdeR-DNA binding. At the residue level, hydrogen bond and non-covalent interactions clearly established the importance of N-terminal residues for DNA binding, thereby confirming the conjecture put forth in the previous chapter. An important aspect studied in this chapter is the allosteric nature of IdeR-DNA binding. Recent years have witnessed a paradigm shift in the understanding of allostery. Unlike the classical definition of allostery that was based on static structures, the newer definition is based on the conformational ensemble as represented by the shift in the energy landscape of the protein. The allosteric nature of IdeR-DNA complex was probed using simulated trajectories and indeed they suggest iron to be an allosteric regulator of the protein. Finally, based on the known experimental data and observations presented in Chapters 4 and 5, a multi-step model of IdeR activation and DNA binding has been proposed. In chapter 6, a global perspective of IdeR regulation in M.tb was obtained. This was important to gain insights about the influences of iron and its regulation at the M.tb cellular level. A genome scale identification of all possible IdeR targets based on the presence of ‘iron-box’ motif in the promoter region of the genes was carried out. An interesting aspect of this study was the use of energetic information from previous molecular dynamics study as an input for generation of the motif. A total of 255 such IdeR targets were identified and converted into an IdeR target network (IdeRnet). Along with IdeRnet, an unbiased systems level protein-protein interaction network was also generated. To study the response of the pathogen to external perturbations, iron-specific gene expression data was integrated into the network as node weights and edge weights. Analysis of IdeRnet provides interesting associations between fatty acid metabolism and IdeR regulations. Specific genes such as fadD32, DesA3 or lppW have been found to be affected by IdeR mutation. While IdeRnet discusses the direct associations, the global level responses are monitored by analysing pathways for the flow of information in the protein-protein interaction network (PPInet). Comparisons of the PPInets under conditions such as altering iron concentrations and lack of iron homeostasis led to the identification of the ‘top-most’ active paths under the different conditions. The study clearly suggests a halt in the protein synthesis machinery and decreased energy consumption under iron scarcity and an uninhibited consumption of energy when iron homeostasis is perturbed. In the final chapter (Chapter 7), flux balance analyses has been used to investigate the influence of iron on M.tb metabolism. The importance of iron for metabolic enzymes has already been established in the previous chapter. Additionally, M.tb is known to produce siderophores, an important metabolite that requires amino acids as its precursors, for iron extraction. All this, together highlighted the importance of iron and its regulation of M.tb metabolism. Flux balance analysis has been used previously to study the metabolic alterations that occur in an organism under different conditions. For this study, iron specific gene expression data was also incorporated into the model as reaction bounds and the flux values so obtained were compared in different environmental conditions. The study provided valuable insights into the metabolic adjustments taken up by M.tb under iron stress conditions and correlates well with the responses observed from the interactome as well as experimental observations. Most significantly, changes were observed in the energy preferences of the cell. For instance, it was noted that while the wild type strain of M.tb prefers synthesis of ATP via glycolysis, the IdeR mutant strain preferred oxidative phosphorylation. The picture becomes clearer when one accounts for the uncontrolled utilization of energy and rapid activation of protein synthesis machinery in the IdeR mutant strain. Biological systems are inherently multiscale in nature and therefore for a successful drug target regime, analysis of the genome to the phenome, which captures interactions at multiple levels, is essential. In this thesis, a detailed understanding of iron homeostasis and regulation in M.tb at multiple levels has been attempted. More importantly, insights obtained from one level, formed questions in the next level. The study was initiated at the inter-cellular level, where the influence of iron on HPI was modeled and analysed. From this study, IdeR, an iron-responsive transcription factor was identified as a key player that had the potential to alter host-pathogen interactions in the favour of the host. For a complete understanding of how IdeR regulates iron homeostasis, it was imperative to obtain a molecular level insight of its mechanism of action. Finally, the various aspects of IdeR regulation were investigated at the cellular level by analysing direct and indirect influences of IdeR on M.tb proteome and metabolome. The study suggests certain therapeutic interventions, such as 1) reduction in the concentration of free transferrin various, 2) mutations at the N-terminal sites of IdeR, 3) regulation of proteins involved in production of mycolic acids by iron and 4) perturbation of altering energy sources, which capitalize on iron and should be investigated in detail. In summary, the consequences of iron on TB infection were studied by threading different levels. This is based on the belief that most biological functions involve multiple spatio-temporal levels with frequent cross talks between the different levels, thereby making such multiscale approaches very useful.

Regulation of transcription in prokaryotes is primarily governed at the transcription initiation step. This feature has been extensively characterized in model prokaryotes notably Escherichia coli and Bacillus subtilis. Transcription initiation was initially thought to be governed primarily by initiation factors that recruit the RNA polymerase (RNAP) enzyme to initiate expression of given gene. Recent studies reveal multiple mechanisms at play including additional protein factors that can modulate gene expression. Nonetheless, understanding transcription factors is key to rationalize the nuanced changes in prokaryotic gene expression in response to diverse environmental stimuli. This is particularly relevant in the case of the human pathogen, Mycobacterium tuberculosis, especially due to the ability of this bacterium to survive in the host, often for several decades prior to the onset of the disease. Transcription initiation factors, also called σ factors in prokaryotes, are diverse in size and sensory/regulatory mechanisms. Indeed, the number of alternate σ factors vary substantially from six in E. coli to more than 118 in Plesiocystis pacifica. The large number of alternative σ factors has been suggested to be correlated with the diversity of micro-environments experienced by a bacterial cell. Studies on several prokaryotic σ factors reveal common features in these proteins that was not evident earlier due to poor sequence conservation. A central theme that emerges from these studies is that a minimalistic architecture of two domains can recognize promoter DNA and recruit the RNAP enzyme to initiate transcription. Additional domains are required when certain promoter elements are missing or to enable a specific, context dependent regulatory mechanism. The work reported in this thesis was influenced by previous studies in this laboratory and elsewhere on M. tuberculosis σ factors. While these studies revealed multiple features of transcription initiation, several aspects of this mechanism, including some classes of σ factors remain to be examined. The focus of this study was to examine an under-explored sub-group of σ factors, classified as the ECF41 sub-group. This sub-group has an additional domain at the Carboxy-terminus that has been hypothesised to influence σ factor activity. Towards this goal, M. tuberculosis σJ was examined. Previous studies suggested a role for this σ factor in modulating the response to hydrogen peroxide stress. An intriguing feature based on sequence analysis was that neither did this extra-cytoplasmic function σ factor have an anti-σ factor that can respond to oxidative stress nor was it directly associated with a mechanism to sense oxidative stress. The specific goal of the research described here was to understand the structural and mechanistic features that govern σJ activity. This thesis is organized as follows- The first chapter provides a brief introduction to prokaryotic transcription and regulatory mechanisms that govern this process. This chapter also has the literature necessary to phrase the problem in characterizing this family of proteins with particular reference to the unique physiology of Mycobacterium tuberculosis. A summary of the previous work is provided in this chapter to place the current study in context of previous studies and highlight the lacunae in our understanding of the transcription mechanism in M. tuberculosis. Chapter two describes the structural characterization of M. tuberculosis σJ by single-crystal X-ray diffraction. The poor sequence similarity of σJ to known σ factors precluded efforts to obtain phase information by molecular replacement methods. Here we also describe the steps that were essential to obtain diffraction quality crystals and the subsequent steps to account for pseudo-merohedral twinning, an imperfection that could have potentially been a limitation for structure determination. The crystal structure of σJ provide an example of successful phase determination with data collected on near-perfectly twinned crystals using single-wavelength anomalous dispersion. Chapter three describes computational efforts to understand the regulatory mechanisms of M. tuberculosis σJ. Classical Molecular Dynamics (MD) simulations were performed to understand the role of a C-terminal SnoaL_2 domain in this transcription factor. The MD simulations suggest that the C-terminal SnoaL_2 domain limits inter-domain movements between σJ2 (the pribnow box binding domain) and σJ4 (the -35 promoter element binding domain) and confers a compact three domain organization to this protein. The biochemical and functional characterization of M. tuberculosis σJ is described in chapter four. This includes in vitro studies on σJ and cognate promoter DNA interactions performed using Surface Plasmon Resonance (SPR) and Electrophoretic Mobility Shift Assays (EMSA). The ex vivo reporter based experiments to examine the effect of SnoaL_2 domain on σJ activity are also described. Spectroscopic studies on σJ interactions with a small molecule limonene-1,2-epoxide suggested a potential novel role for the SnoaL_2 domain in σJ. Chapter five summarizes the work on M. tuberculosis σJ reported in this thesis. We note that this study opens up a new perspective to understand σ factors. In particular, M. tuberculosis σJ suggests that the domain organization is likely to be retained in ECF41 sub-group of σ factors. This study also hints at broader implications in the distinction between one-component systems and transcription factors. Bioinformatic analysis suggest that observations similar to that noted in M. tuberculosis σJ are likely to be more widespread across diverse phyla than currently acknowledged. This thesis has three annexures. Annexure-I summarizes experimental details of the work performed on the M. tuberculosis σ/anti-σ factor complex σH/RshA. Annexure-II summarizes experimental details and strategies that could not be incorporated in the main body of this thesis. Annexure-III describes a short project performed on a bi-domain protein tyrosine phosphatase PTP99A.

Mycobacterium tuberculosis, the etiological agent of tuberculosis, has adapted with the host environment and evolved to survive in harsh conditions in the host. The pathogen has successfully evolved strategies not only to evade the host immune system but also to thrive within the host cells. Upon infection, the pathogen is either cleared due to the host immune response, or it survives and causes active tuberculosis (TB) infection. In a number of cases however, the pathogen is neither killed nor does it actively proliferate, but it remains dormant in the host until the environment becomes favorable. This dormant state of pathogen is responsible for latent TB infection (LTBI). WHO reports indicated that as much as a third of the whole world’s population is exposed to the pathogen, of which a significant proportion could be latently infected (WHO report, 2015). These individuals do not show symptoms of active TB infection and hence are difficult to detect. The latent TB infected (LTBI) individuals serve as a reservoir for the pathogen, which can lead to epidemics when the conditions change. Hence, it is necessary to understand the host -pathogen interactions during LTBI, as this might provide clues to developing new strategies to detect and curb a latent infection. Host-pathogen interactions are multifaceted, in which both species attempt to recognize and respond to each other, all of these through specific molecules making distinct interactions with the other species. The outcome of the infection is thus decided by a complex set of host-pathogen interactions. The complexity arises since a large number of molecular components are involved, also multiplicity of interactions among these components and due to several feedback, feed forwards or other regulatory or influential loops within the system. The complexity of biological systems makes modeling and simulation an essential and critical part of systems– level studies. Systems biology studies provide an integrated framework to analyze and understand the function of biological systems. This work addresses some of these issues with an unbiased systems-level analysis so as to identify and understand the important global changes both in the host and in the pathogen during LTBI. The broad objectives of the work was to identify the key processes that vary in the host during latent infection, the set of metabolic reactions in the host which can be modulated to control the reactivation of infection, global adaptation in Mycobacterium tuberculosis (Mtb) and then to utilize this knowledge to identify strategies for tackling latent infection. A review of literature of the current understanding of latency from the pathogen and the host perspective is described in chapter 1. From this, it is clear that most available studies have focused on the role of individual molecules and individual biological processes such as granuloma formation, toll-like receptor signaling, T cell responses as well as cytokine signaling, in either initiating or maintaining a latent infection, but there is no report till date about whether and how these processes are connected with each other. While transcriptome based studies have identified lists of differentially expressed genes in LTBI as compared to healthy controls, no further understanding is currently available for many of them, regarding the processes they may be involved in and what interactions they make, which may be important for understanding LTBI. The first part of the work is a systematic meta-analysis of genome-scale protein interaction networks rendered condition-specific with transcriptome data of patients with LTBI, which has provided a global unbiased picture of the transcriptional and metabolic variations in the host and in the pathogen during the latent infection. To start with, publicly available gene expression data related to LTBI, active TB and healthy controls were considered. In all, 183 datasets summing up to 105 LTBI, 41 active TB and 37 healthy control samples were analyzed. (Chapter 2). Standard analysis of the transcriptome profiles of these datasets indicated that there was zero overlap among them and that not a single gene was seen in common among all datasets for the same condition. An extensive human protein-protein interaction network was constructed using information available from multiple resources that comprehensively contained structural or physical interactions and genetic interactions or functional influences. Nodes in this network represented individual proteins and edges represented interactions between pairs of nodes. The identity of each node and the nature of interaction of each edge along with the type of evidence that was used as the basis for drawing the edge, was collated for the network. The gene expression data was integrated into the human protein-protein interaction (PPI) network for each condition, which essentially had weighted nodes and directed edges, specific to that condition, from which specific comparative networks were derived. The highest ranked perturbations in LTBI were identified through a network mining protocol previously established in the laboratory. This involved computing all versus all shortest paths on the comparative network, scoring the paths based on connectedness and various centrality measures of the nodes and the edges and finally ranking the paths based on the cumulative path scores. Intriguingly, the top-ranked set of perturbations were found to form a connected sub-network by themselves, referred to as a top perturbed sub-network (top-net), indicating that they were functionally linked or perhaps even orchestrated in some sense. Th17 signaling appears to be dominant. About 40 genes were identified in the unique set of LTBI condition as compared to the active TB condition, and these genes showed enrichment for processes such as apoptosis, cell cycle as well as natural killer cell mediated toxicity. Construction and analysis of a miRNA network indicated that 32 of these have strong associations with miRNA explaining the role of the latter in controlling LTBI. 3 other genes from the top-net are already established drug targets for different diseases with known drugs associated with them, which are BCL2, HSP90AA1 and NR3C1. These 3 proteins can be explored further as drug targets in LTBI whose manipulation using existing drugs may result in inhibiting the underlying biological process and thereby result in disturbing the state of latency. As a second objective, global variations in the host transcriptome were identified during ascorbic acid induced dormancy (Chapter 3). Ascorbic acid or Vitamin C is a nutrient supplement required in the diet. This organic compound has a known antioxidant property, as it is known to scavenge the free radicals. In a recent study, Taneja et al, demonstrated that Vitamin C could induce dormancy in Mtb. On similar lines, experiments were done in THP-1 cells infected with Mtb to determine the host responses during ascorbic acid (AA) induced dormancy. The raw gene expression data was provided by our collaborator Prof. Jaya Tyagi that included 0 hour, 4 days and 6 days time points with infection and vitamin C versus infection alone or vitamin C alone as controls. The transcriptome data was normalized and integrated into the human PPI network as described for the meta-analyses. It was experimentally determined that ascorbic acid induces dormancy in 4 days post infection. The top-ranked paths of perturbation were analyzed and compared for three different conditions: (i) uninfected condition, (ii) AA treated and infected condition, and (iii) AA, isoniazid and infected condition. The dormant pathogen is known to be drug-tolerant and thus as a marker for the state of dormancy, the lack of effect of isoniazid is also monitored in the infected host cells. The analysis revealed that there were some broad similarities as compared to LTBI from patient samples but AA induced dormancy in cell lines stood out a separate group indicating that there were significant differences such as involving Interferon Induced Transmembrane Proteins (IFITMs), vacuolar ATPase as well as GDF15, which belongs to TGF-beta signaling pathway. The highest ranked perturbed paths contained genes involved in innate immune responses of which ISG15, IFITMs, HLAs and ATPases emerge as the most altered in the dormant condition. CCR7 emerges as a key discriminator, which is subdued in the latent samples but highly induced in infection conditions. Pathway-based analysis of different conditions showed that oxidative stress, glutathione metabolism, proteasome degradation as well as type II interferon signaling are significantly up-regulated in AA induced dormancy. The dormant bacteria reside in the host cells and are known to modulate the host metabolism for their own benefit. So, the third objective was to understand the metabolic variations in the host during LTBI (Chapter 4). A genome-scale metabolic (GSM) model of alveolar macrophage was used in this study. The metabolic model contains information of the reactions, metabolites and the genes encoding enzymes that catalyze a particular reaction. Flux balance analysis (FBA), a constraint-based metabolic modeling method, is used for analyzing the alterations in the metabolism under different infection conditions. In order to mimic the physiological condition, gene expression data was used for constraining the bounds of the reactions in the model. Two different expression studies were used for analysis: GSE25534 (from Chapter 2) and ascorbic acid induced dormancy (Chapter 3). The analysis was carried out for latent TB versus healthy control and latent TB versus active TB to identify the most altered metabolic processes in LTBI. Differences in fluxes between the two conditions were calculated. A new classification scheme was devised to categorize the reactions on the basis of flux differences. In this chapter, higher fluxes in LTBI condition were identified for reactions involved in transport of small metabolites as well as amino acids. Solute carrier proteins responsible for the transport of the metabolites were identified and their biological significance is discussed. Reduced glutathione (GSH), arachidonic acid, prostaglandins, pantothenate were identified as important metabolites in LTBI condition and their physiological role has been described. Sub-system analysis for different conditions shows differential regulation for arachidonic acid metabolism, fatty acid metabolism, folate metabolism, pyruvate metabolism, glutathione metabolism, ROS detoxification, triacylglycerol synthesis and transport as well as tryptophan metabolism. From the study, transporter proteins and reactions altered during LTBI were identified, which again provide clues for understanding the molecular basis of establishing a latent infection. Mycoabcterium tuberculosis is known to undergo dormancy during stress conditions. In this chapter, the main objective was to identify the global variations in the dormant Mtb (Chapter 5). To carry out the analysis, the Mtb PPI network was constructed using information from available resources. Gene expression data of two different dormancy models, Wayne growth model and multiple-stress model, were used for the study. To identify the key players involved in reversal of dormancy, the transcriptome data of reaeration condition was also used. In this study, the Max-flow algorithm was implemented to identify the feasible paths or flows in different condition. The flows with higher scores indicate that more information is traversed by the path, and hence is important for the study. From the analysis of Wayne growth model (hypoxia model), important transcriptional regulators such as SigB, SigE, SigH, regulators in the two-component system such as MprA, MtrA, PhoP, RegX3 and TrcR were identified in stress condition. Multiple-stress model studied the growth of bacteria in low oxygen concentration, high carbon dioxide levels, low pH and nutrient starvation. The gene expression data was integrated in the Mtb PPI network and implementation of Max-flow algorithm showed that MprA, part of the MprA-MprB two-component system, is involved in the regulation of persistent condition. WhiB1 also features in the paths of dormant condition and its role in persistence can be explored. In reaeration model, WhiB1 and WhiB4 are present in the top flows of this condition indicating that the redox state is perturbed in the pathogen and the interactions of these proteins are important to understand the reversal of dormant condition. From the study, Rv2034, Rv2035, HigA, Rv1989, Rv1990 and Rv0837 proteins belonging to toxin-antitoxin systems were also identified in the dormant bacteria, indicating their role in adaptation during stress condition. The role of Rv2034 has been studied in persistence, but the function of other proteins can be analyzed to provide new testable hypotheses about the role of these proteins in dormancy. Thus, the flows or paths perturbed during dormancy were identified in this study. To get a better understanding of the metabolic network active in mycobacteria under different conditions, experiments were performed in Mycobacterium smegmatis MC2 155. The non-pathogenic strain of genus Mycobacteria, Mycobacterium smegmatis, is used as a surrogate to carry out molecular biology studies of Mtb. Mycobacterium smegmatis MC2 155 (Msm) is the commonly used laboratory strain for experimental purpose. In order to obtain a clear understanding of how comparable are the metabolic networks between the virulent M. tuberculosis H37Rv and the model system Msm, the latter model is first studied systematically. In Chapter 6, first the functional annotation of the Msm genome was carried out and the genes were categorized into different Tuberculist classes based on homology with the Mtb genome. A high-throughput growth characterization was carried out to characterize the strain systematically in terms of different carbon, nitrogen or other sources that promoted growth and thus served as nutrients and those that did not, together yielding a genome-phenome correlation in Msm. Gene expression was measured and used for explaining the observed phenotypic behavior of the organism. Together with the genome sequence, the transcriptome and phenome analysis, a set of about 257 different metabolic pathways were identified to be feasible in wild-type Msm. About 284 different carbon, nitrogen source and nutrient supplements were tested in this experiment and 167 of them supported growth of Msm. This indicates that the compounds enter the cells and are metabolized efficiently, thus yielding similar phenotypes. The expressed genes and metabolites supporting growth were mapped to the metabolic network of Msm, thus helping in the identification of feasible metabolic routes in Msm. A comparative study between Msm and Mtb revealed that these organisms share similarity in the nutrient sources that are utilized for growth. The study provides experimental proof to identify the feasible metabolic routes in Msm, and this can be used for understanding the metabolic capability in the two organisms under different conditions providing a basis to understand adaptations during dormancy. In the last part of the work presented in this thesis, the metabolic shift in the pathogen was studied using a genome-scale metabolic model of Mtb (Chapter 7). The model contains information of the reactions, metabolites and genes involved in the reactions. Flux balance analysis (FBA) was carried out by integrating normalized gene expression data (Wayne model and multiple-stress model transcriptome considered in Chapter 5) to identify the set of reactions, which have a higher flux in the dormant condition as compared to the control replicating condition. Glutamate metabolism along with propionyl CoA metabolism emerge as major up-regulated processes in dormant Mtb. Next, with an objective of identifying essential genes in dormant Mtb, a systematic in silico single gene knock-out analysis was carried out where each gene and it's associated reaction was knocked out of the model, one at a time and the ability of the model to reach its objective function assessed. About 168 common genes in Wayne model and multiple-stress model were identified as important in Mtb after the knockout analysis. Essentiality is in essence a systems property and requires to be probed through multiple angles. Towards this, essential genes were identified in Mtb using a multi-level multi-scale systems biology approach. About 283 genes were identified as essential on the basis of combined analysis of transcriptome data, FBA, network analysis and phyletic retention studies in Mtb. 168 genes identified as important in dormant Mtb were compared with 283 essential genes and about 91 genes were found to be essential. Finally, among the set of essential genes, those that satisfy other criteria for a drug target were analyzed using the list of high-confidence drug targets of Mtb available in the laboratory along with their associated drug or drug-like molecules. 38 out of the 168 important genes in Mtb were found to have one or more drugs associated with them from the DrugBank database. Colchicin-Rv1655, Raloxifene-Rv1653, Bexarotene-Rv3804, Rosiglitazone-Rv3804 are top-scoring drug-target pairs that can be explored for killing dormant bacilli. The study has thus been useful in identifying important proteins, reactions and drug targets in dormant Mtb. In summary, the thesis presents a comprehensive systems-level understanding of various aspects of host responses and pathogen adaptation during latent TB infection. Key host and pathogen factors involved in LTBI are identified that serve as useful pointers for deriving strategies for tackling a latent infection.

No description available.