Structural Annotation of Mycobacterium tuberculosis Proteome

Of the similar to 4000 ORFs identified through the genome sequence of Mycobacterium tuberculosis (TB) H37Rv, experimentally determined structures are available for 312. Since knowledge of protein structures is essential to obtain a high-resolution understanding of the underlying biology, we seek to obtain a structural annotation for the genome, using computational methods. Structural models were obtained and validated for similar to 2877 ORFs, covering similar to 70% of the genome. Functional annotation of each protein was based on fold-based functional assignments and a novel binding site based ligand association. New algorithms for binding site detection and genome scale binding site comparison at the structural level, recently reported from the laboratory, were utilized. Besides these, the annotation covers detection of various sequence and sub-structural motifs and quaternary structure predictions based on the corresponding templates. The study provides an opportunity to obtain a global perspective of the fold distribution in the genome. The annotation indicates that cellular metabolism can be achieved with only 219 folds. New insights about the folds that predominate in the genome, as well as the fold-combinations that make up multi-domain proteins are also obtained. 1728 binding pockets have been associated with ligands through binding site identification and sub-structure similarity analyses. The resource (http://proline.physics.iisc.ernet.in/Tbstructuralannotation), being one of the first to be based on structure-derived functional annotations at a genome scale, is expected to be useful for better understanding of TB and for application in drug discovery. The reported annotation pipeline is fairly generic and can be applied to other genomes as well.

UDP-glucose 4, 6-dehydratase Activity Plays an Important Role in Maintaining Cell Wall Integrity and Virulence of Candida albicans

Candida albicans, a human fungal pathogen, undergoes morphogenetic changes that are associated with virulence. We report here that GAL102 in C. albicans encodes a homolog of dTDP-glucose 4,6-dehydratase, an enzyme that affects cell wall properties as well as virulence of many pathogenic bacteria. We found that GAL102 deletion leads to greater sensitivity to antifungal drugs and cell wall destabilizing agents like Calcofluor white and Congo red. The mutant also formed biofilms consisting mainly of hyphal cells that show less turgor. The NMR analysis of cell wall mannans of gal102 deletion strain revealed that a major constituent of mannan is missing and the phosphomannan component known to affect virulence is greatly reduced. We also observed that there was a substantial reduction in the expression of genes involved in biofilm formation but increase in the expression of genes encoding glycosylphosphatidylinositol-anchored proteins in the mutant. These, along with altered mannosylation of cell wall proteins together might be responsible for multiple phenotypes displayed by the mutant. Finally, the mutant was unable to grow in the presence of resident peritoneal macrophages and elicited a weak pro-inflammatory cytokine response in vitro. Similarly, this mutant elicited a poor serum pro-inflammatory cytokine response as judged by IFN gamma and TNF alpha levels and showed reduced virulence in a mouse model of systemic candidiasis. Importantly, an Ala substitution for a conserved Lys residue in the active site motif YXXXK, that abrogates the enzyme activity also showed reduced virulence and increased filamentation similar to the gal102 deletion strain. Since inactivating the enzyme encoded by GAL102 makes the cells sensitive to antifungal drugs and reduces its virulence, it can serve as a potential drug target in combination therapies for C. albicans and related pathogens.

Chimeras of Escherichia coli and Mycobacterium tuberculosis Single-Stranded DNA Binding Proteins: Characterization and Function in Escherichia coli

Single stranded DNA binding proteins (SSBs) are vital for the survival of organisms. Studies on SSBs from the prototype, Escherichia coli (EcoSSB) and, an important human pathogen, Mycobacterium tuberculosis (MtuSSB) had shown that despite significant variations in their quaternary structures, the DNA binding and oligomerization properties of the two are similar. Here, we used the X-ray crystal structure data of the two SSBs to design a series of chimeric proteins (m beta 1, m beta 1'beta 2, m beta 1-beta 5, m beta 1-beta 6 and m beta 4-beta 5) by transplanting beta 1, beta 1'beta 2, beta 1-beta 5, beta 1-beta 6 and beta 4-beta 5 regions, respectively of the N-terminal (DNA binding) domain of MtuSSB for the corresponding sequences in EcoSSB. In addition, m beta 1'beta 2(ESWR) SSB was generated by mutating the MtuSSB specific `PRIY' sequence in the beta 2 strand of m beta 1'beta 2 SSB to EcoSSB specific `ESWR' sequence. Biochemical characterization revealed that except for m beta 1 SSB, all chimeras and a control construct lacking the C-terminal domain (Delta C SSB) bound DNA in modes corresponding to limited and unlimited modes of binding. However, the DNA on MtuSSB may follow a different path than the EcoSSB. Structural probing by protease digestion revealed that unlike other SSBs used, m beta 1 SSB was also hypersensitive to chymotrypsin treatment. Further, to check for their biological activities, we developed a sensitive assay, and observed that m beta 1-beta 6, MtuSSB, m beta 1'beta 2 and m beta 1-beta 5 SSBs complemented E. coli Delta ssb in a dose dependent manner. Complementation by the m beta 1-beta 5 SSB was poor. In contrast, m beta 1'beta 2(ESWR) SSB complemented E. coli as well as EcoSSB. The inefficiently functioning SSBs resulted in an elongated cell/filamentation phenotype of E. coli. Taken together, our observations suggest that specific interactions within the DNA binding domain of the homotetrameric SSBs are crucial for their biological function.

Retention Time Variability as a Mechanism for Animal Mediated Long-Distance Dispersal

Long-distance dispersal (LDD) events, although rare for most plant species, can strongly influence population and community dynamics. Animals function as a key biotic vector of seeds and thus, a mechanistic and quantitative understanding of how individual animal behaviors scale to dispersal patterns at different spatial scales is a question of critical importance from both basic and applied perspectives. Using a diffusion-theory based analytical approach for a wide range of animal movement and seed transportation patterns, we show that the scale (a measure of local dispersal) of the seed dispersal kernel increases with the organisms' rate of movement and mean seed retention time. We reveal that variations in seed retention time is a key determinant of various measures of LDD such as kurtosis (or shape) of the kernel, thinkness of tails and the absolute number of seeds falling beyond a threshold distance. Using empirical data sets of frugivores, we illustrate the importance of variability in retention times for predicting the key disperser species that influence LDD. Our study makes testable predictions linking animal movement behaviors and gut retention times to dispersal patterns and, more generally, highlights the potential importance of animal behavioral variability for the LDD of seeds.

Mathematical Model of Viral Kinetics In Vitro Estimates the Number of E2-CD81 Complexes Necessary for Hepatitis C Virus Entry

Interaction between the hepatitis C virus (HCV) envelope protein E2 and the host receptor CD81 is essential for HCV entry into target cells. The number of E2-CD81 complexes necessary for HCV entry has remained difficult to estimate experimentally. Using the recently developed cell culture systems that allow persistent HCV infection in vitro, the dependence of HCV entry and kinetics on CD81 expression has been measured. We reasoned that analysis of the latter experiments using a mathematical model of viral kinetics may yield estimates of the number of E2-CD81 complexes necessary for HCV entry. Here, we constructed a mathematical model of HCV viral kinetics in vitro, in which we accounted explicitly for the dependence of HCV entry on CD81 expression. Model predictions of viral kinetics are in quantitative agreement with experimental observations. Specifically, our model predicts triphasic viral kinetics in vitro, where the first phase is characterized by cell proliferation, the second by the infection of susceptible cells and the third by the growth of cells refractory to infection. By fitting model predictions to the above data, we were able to estimate the threshold number of E2-CD81 complexes necessary for HCV entry into human hepatoma-derived cells. We found that depending on the E2-CD81 binding affinity, between 1 and 13 E2-CD81 complexes are necessary for HCV entry. With this estimate, our model captured data from independent experiments that employed different HCV clones and cells with distinct CD81 expression levels, indicating that the estimate is robust. Our study thus quantifies the molecular requirements of HCV entry and suggests guidelines for intervention strategies that target the E2-CD81 interaction. Further, our model presents a framework for quantitative analyses of cell culture studies now extensively employed to investigate HCV infection.

Curcumin-Arteether Combination Therapy of Plasmodium berghei-Infected Mice Prevents Recrudescence Through Immunomodulation

Earlier studies in this laboratory have shown the potential of artemisinin-curcumin combination therapy in experimental malaria. In a parasite recrudescence model in mice infected with Plasmodium berghei (ANKA), a single dose of alpha, beta-arteether (ART) with three oral doses of curcumin prevented recrudescence, providing almost 95% protection. The parasites were completely cleared in blood with ART-alone (AE) or ART+curcumin (AC) treatments in the short-term, although the clearance was faster in the latter case involving increased ROS generation. But, parasites in liver and spleen were not cleared in AE or AC treatments, perhaps, serving as a reservoir for recrudescence. Parasitemia in blood reached up to 60% in AE-treated mice during the recrudescence phase, leading to death of animals. A transient increase of up to 2-3% parasitemia was observed in AC-treatment, leading to protection and reversal of splenomegaly. A striking increase in spleen mRNA levels for TLR2, IL-10 and IgG-subclass antibodies but a decrease in those for INF gamma and IL-12 was observed in AC-treatment. There was a striking increase in IL-10 and IgG subclass antibody levels but a decrease in INF gamma levels in sera leading to protection against recrudescence. AC-treatment failed to protect against recrudescence in TLR2(-/-) and IL-10(-/-) animals. IL-10 injection to AE-treated wild type mice and AC-treated TLR22/2 mice was able to prolong survival. Blood from the recrudescence phase in AE-treatment, but not from AC-treatment, was able to reinfect and kill naive animals. Sera from the recrudescence phase of AC-treated animals reacted with several parasite proteins compared to that from AE-treated animals. It is proposed that activation of TLR2-mediated innate immune response leading to enhanced IL-10 production and generation of anti-parasite antibodies contribute to protective immunity in AC-treated mice. These results indicate a potential for curcumin-based combination therapy to be tested for prevention of recrudescence in falciparum and relapse in vivax malaria.

A Unique Modification of the Eukaryotic Initiation Factor 5A Shows the Presence of the Complete Hypusine Pathway in Leishmania donovani

Deoxyhypusine hydroxylase (DOHH) catalyzes the final step in the post-translational synthesis of an unusual amino acid hypusine (N-(sic)-(4-amino-2-hydroxybutyl) lysine), which is present on only one cellular protein, eukaryotic initiation factor 5A (eIF5A). We present here the molecular and structural basis of the function of DOHH from the protozoan parasite, Leishmania donovani, which causes visceral leishmaniasis. The L. donovani DOHH gene is 981 bp and encodes a putative polypeptide of 326 amino acids. DOHH is a HEAT-repeat protein with eight tandem repeats of alpha-helical pairs. Four conserved histidine-glutamate sequences have been identified that may act as metal coordination sites. A similar to 42 kDa recombinant protein with a His-tag was obtained by heterologous expression of DOHH in Escherichia coli. Purified recombinant DOHH effectively catalyzed the hydroxylation of the intermediate, eIF5A-deoxyhypusine (eIF5A-Dhp), in vitro. L. donovani DOHH (LdDOHH) showed similar to 40.6% sequence identity with its human homolog. The alignment of L. donovani DOHH with the human homolog shows that there are two significant insertions in the former, corresponding to the alignment positions 159-162 (four amino acid residues) and 174-183 (ten amino acid residues) which are present in the variable loop connecting the N- and C-terminal halves of the protein, the latter being present near the substrate binding site. Deletion of the ten-amino-acid-long insertion decreased LdDOHH activity to 14% of the wild type recombinant LdDOHH. Metal chelators like ciclopirox olamine (CPX) and mimosine significantly inhibited the growth of L. donovani and DOHH activity in vitro. These inhibitors were more effective against the parasite enzyme than the human enzyme. This report, for the first time, confirms the presence of a complete hypusine pathway in a kinetoplastid unlike eubacteria and archaea. The structural differences between the L. donovani DOHH and the human homolog may be exploited for structure based design of selective inhibitors against the parasite.

Sesbania Mosaic Virus (SeMV) Infectious Clone: Possible Mechanism of 3 ` and 5 ` End Repair and Role of Polyprotein Processing in Viral Replication

Sesbania mosaic virus (SeMV) is a positive stranded RNA virus belonging to the genus Sobemovirus. Construction of an infectious clone is an essential step for deciphering the virus gene functions in vivo. Using Agrobacterium based transient expression system we show that SeMV icDNA is infectious on Sesbania grandiflora and Cyamopsis tetragonoloba plants. The efficiency of icDNA infection was found to be significantly high on Cyamopsis plants when compared to that on Sesbania grandiflora. The coat protein could be detected within 6 days post infiltration in the infiltrated leaves. Different species of viral RNA (double stranded and single stranded genomic and subgenomic RNA) could be detected upon northern analysis, suggesting that complete replication had taken place. Based on the analysis of the sequences at the genomic termini of progeny RNA from SeMV icDNA infiltrated leaves and those of its 3' and 5' terminal deletion mutants, we propose a possible mechanism for 3' and 5' end repair in vivo. Mutation of the cleavage sites in the polyproteins encoded by ORF 2 resulted in complete loss of infection by the icDNA, suggesting the importance of correct polyprotein processing at all the four cleavage sites for viral replication. Complementation analysis suggested that ORF 2 gene products can act in trans. However, the trans acting ability of ORF 2 gene products was abolished upon deletion of the N-terminal hydrophobic domain of polyprotein 2a and 2ab, suggesting that these products necessarily function at the replication site, where they are anchored to membranes.

Extent of Structural Asymmetry in Homodimeric Proteins: Prevalence and Relevance

Most homodimeric proteins have symmetric structure. Although symmetry is known to confer structural and functional advantage, asymmetric organization is also observed. Using a non-redundant dataset of 223 high-resolution crystal structures of biologically relevant homodimers, we address questions on the prevalence and significance of asymmetry. We used two measures to quantify global and interface asymmetry, and assess the correlation of several molecular and structural parameters with asymmetry. We have identified rare cases (11/223) of biologically relevant homodimers with pronounced global asymmetry. Asymmetry serves as a means to bring about 2:1 binding between the homodimer and another molecule; it also enables cellular signalling arising from asymmetric macromolecular ligands such as DNA. Analysis of these cases reveals two possible mechanisms by which possible infinite array formation is prevented. In case of homodimers associating via non-topologically equivalent surfaces in their tertiary structures, ligand-dependent mechanisms are used. For stable dimers binding via large surfaces, ligand-dependent structural change regulates polymerisation/depolymerisation; for unstable dimers binding via smaller surfaces that are not evolutionarily well conserved, dimerisation occurs only in the presence of the ligand. In case of homodimers associating via interaction surfaces with parts of the surfaces topologically equivalent in the tertiary structures, steric hindrance serves as the preventive mechanism of infinite array. We also find that homodimers exhibiting grossly symmetric organization rarely exhibit either perfect local symmetry or high local asymmetry. Binding of small ligands at the interface does not cause any significant variation in interface asymmetry. However, identification of biologically relevant interface asymmetry in grossly symmetric homodimers is confounded by the presence of similar small magnitude changes caused due to artefacts of crystallisation. Our study provides new insights regarding accommodation of asymmetry in homodimers.

Structural and Mutational Studies on Substrate Specificity and Catalysis of Salmonella typhimurium D-Cysteine Desulfhydrase

Salmonella typhimurium DCyD (StDCyD) is a fold type II pyridoxal 5' phosphate (PLP)-dependent enzyme that catalyzes the degradation of D-Cys to H2S and pyruvate. It also efficiently degrades beta-chloro-D-alanine (beta CDA). D-Ser is a poor substrate while the enzyme is inactive with respect to L-Ser and 1-amino-1-carboxy cyclopropane (ACC). Here, we report the X-ray crystal structures of StDCyD and of crystals obtained in the presence of D-Cys, beta CDA, ACC, D-Ser, L-Ser, D-cycloserine (DCS) and L-cycloserine (LCS) at resolutions ranging from 1.7 to 2.6 angstrom. The polypeptide fold of StDCyD consisting of a small domain (residues 48-161) and a large domain (residues 1-47 and 162-328) resembles other fold type II PLP dependent enzymes. The structures obtained in the presence of D-Cys and beta CDA show the product, pyruvate, bound at a site 4.0-6.0 angstrom away from the active site. ACC forms an external aldimine complex while D- and L-Ser bind non-covalently suggesting that the reaction with these ligands is arrested at C alpha proton abstraction and transimination steps, respectively. In the active site of StDCyD cocrystallized with DCS or LCS, electron density for a pyridoxamine phosphate (PMP) was observed. Crystals soaked in cocktail containing these ligands show density for PLP-cycloserine. Spectroscopic observations also suggest formation of PMP by the hydrolysis of cycloserines. Mutational studies suggest that Ser78 and Gln77 are key determinants of enzyme specificity and the phenolate of Tyr287 is responsible for C alpha proton abstraction from D-Cys. Based on these studies, a probable mechanism for the degradation of D-Cys by StDCyD is proposed.

Viral Kinetics Suggests a Reconciliation of the Disparate Observations of the Modulation of Claudin-1 Expression on Cells Exposed to Hepatitis C Virus

The tight junction protein claudin-1 (CLDN1) is necessary for hepatitis C virus (HCV) entry into target cells. Recent studies have made disparate observations of the modulation of the expression of CLDN1 on cells following infection by HCV. In one study, the mean CLDN1 expression on cells exposed to HCV declined, whereas in another study HCV infected cells showed increased CLDN1 expression compared to uninfected cells. Consequently, the role of HCV in modulating CLDN1 expression, and hence the frequency of cellular superinfection, remains unclear. Here, we present a possible reconciliation of these disparate observations. We hypothesized that viral kinetics and not necessarily HCV-induced receptor modulation underlies these disparate observations. To test this hypothesis, we constructed a mathematical model of viral kinetics in vitro that mimicked the above experiments. Model predictions provided good fits to the observed evolution of the distribution of CLDN1 expression on cells following exposure to HCV. Cells with higher CLDN1 expression were preferentially infected and outgrown by cells with lower CLDN1 expression, resulting in a decline of the mean CLDN1 expression with time. At the same time, because the susceptibility of cells to infection increased with CLDN1 expression, infected cells tended to have higher CLDN1 expression on average than uninfected cells. Our study thus presents an explanation of the disparate observations of CLDN1 expression following HCV infection and points to the importance of considering viral kinetics in future studies of receptor expression on cells exposed to HCV.

Binding of Gemini Bisbenzimidazole Drugs with Human Telomeric G-Quadruplex Dimers: Effect of the Spacer in the Design of Potent Telomerase Inhibitors

The study of anticancer agents that act via stabilization of telomeric G-quadruplex DNA (G4DNA) is important because such agents often inhibit telomerase activity. Several types of G4DNA binding ligands are known. In these studies, the target structures often involve a single G4 DNA unit formed by short DNA telomeric sequences. However, the 3'-terminal single-stranded human telomeric DNA can form higher-order structures by clustering consecutive quadruplex units (dimers or nmers). Herein, we present new synthetic gemini (twin) bisbenzimidazole ligands, in which the oligo-oxyethylene spacers join the two bisbenzimidazole units for the recognition of both monomeric and dimeric G4DNA, derived from d(T2AG3)4 and d(T2AG3) 8 human telomeric DNA, respectively. The spacer between the two bisbenzimidazoles in the geminis plays a critical role in the G4DNA stability. We report here (i) synthesis of new effective gemini anticancer agents that are selectively more toxic towards the cancer cells than the corresponding normal cells; (ii) formation and characterization of G4DNA dimers in solution as well as computational construction of the dimeric G4DNA structures. The gemini ligands direct the folding of the single-stranded DNA into an unusually stable parallel-stranded G4DNA when it was formed in presence of the ligands in KCl solution and the gemini ligands show spacer length dependent potent telomerase inhibition properties.

A Novel Property of DNA - As a Bioflotation Reagent in Mineral Processing

Environmental concerns regarding the use of certain chemicals in the froth flotation of minerals have led investigators to explore biological entities as potential substitutes for the reagents in vogue. Despite the fact that several microorganisms have been used for the separation of a variety of mineral systems, a detailed characterization of the biochemical molecules involved therein has not been reported so far. In this investigation, the selective flotation of sphalerite from a sphalerite-galena mineral mixture has been achieved using the cellular components of Bacillus species. The key constituent primarily responsible for the flotation of sphalerite has been identified as DNA, which functions as a bio-collector. Furthermore, using reconstitution studies, the obligatory need for the presence of non-DNA components as bio-depressants for galena has been demonstrated. A probable model involving these entities in the selective flotation of sphalerite from the mineral mixture has been discussed.

Molecular Epidemiology of HIV-1 Subtypes in India: Origin and Evolutionary History of the Predominant Subtype C

Background: India has the third largest HIV-1 epidemic with 2.4 million infected individuals. Molecular epidemiological analysis has identified the predominance of HIV-1 subtype C (HIV-1C). However, the previous reports have been limited by sample size, and uneven geographical distribution. The introduction of HIV-1C in India remains uncertain due to this lack of structured studies. To fill the gap, we characterised the distribution pattern of HIV-1 subtypes in India based on data collection from nationwide clinical cohorts between 2007 and 2011. We also reconstructed the time to the most recent common ancestor (tMRCA) of the predominant HIV-1C strains. Methodology/Principal Findings: Blood samples were collected from 168 HIV-1 seropositive subjects from 7 different states. HIV-1 subtypes were determined using two or three genes, gag, pol, and env using several methods. Bayesian coalescent-based approach was used to reconstruct the time of introduction and population growth patterns of the Indian HIV-1C. For the first time, a high prevalence (10%) of unique recombinant forms (BC and A1C) was observed when two or three genes were used instead of one gene (p<0.01; p = 0.02, respectively). The tMRCA of Indian HIV-1C was estimated using the three viral genes, ranged from 1967 (gag) to 1974 (env). Pol-gene analysis was considered to provide the most reliable estimate 1971, (95% CI: 1965-1976)]. The population growth pattern revealed an initial slow growth phase in the mid-1970s, an exponential phase through the 1980s, and a stationary phase since the early 1990s. Conclusions/Significance: The Indian HIV-1C epidemic originated around 40 years ago from a single or few genetically related African lineages, and since then largely evolved independently. The effective population size in the country has been broadly stable since the 1990s. The evolving viral epidemic, as indicated by the increase of recombinant strains, warrants a need for continued molecular surveillance to guide efficient disease intervention strategies.

Insights into the Fold Organization of TIM Barrel from Interaction Energy Based Structure Networks

There are many well-known examples of proteins with low sequence similarity, adopting the same structural fold. This aspect of sequence-structure relationship has been extensively studied both experimentally and theoretically, however with limited success. Most of the studies consider remote homology or ``sequence conservation'' as the basis for their understanding. Recently ``interaction energy'' based network formalism (Protein Energy Networks (PENs)) was developed to understand the determinants of protein structures. In this paper we have used these PENs to investigate the common non-covalent interactions and their collective features which stabilize the TIM barrel fold. We have also developed a method of aligning PENs in order to understand the spatial conservation of interactions in the fold. We have identified key common interactions responsible for the conservation of the TIM fold, despite high sequence dissimilarity. For instance, the central beta barrel of the TIM fold is stabilized by long-range high energy electrostatic interactions and low-energy contiguous vdW interactions in certain families. The other interfaces like the helix-sheet or the helix-helix seem to be devoid of any high energy conserved interactions. Conserved interactions in the loop regions around the catalytic site of the TIM fold have also been identified, pointing out their significance in both structural and functional evolution. Based on these investigations, we have developed a novel network based phylogenetic analysis for remote homologues, which can perform better than sequence based phylogeny. Such an analysis is more meaningful from both structural and functional evolutionary perspective. We believe that the information obtained through the ``interaction conservation'' viewpoint and the subsequently developed method of structure network alignment, can shed new light in the fields of fold organization and de novo computational protein design.