Structure-sequence features based prediction of phosphosites of serine/threonine protein kinases of Mycobacterium tuberculosis.

Elucidation of signaling events in a pathogen is potentially important to tackle the infection caused by it. Such events mediated by protein phosphorylation play important roles in infection, and therefore, to predict the phosphosites and substrates of the serine/threonine protein kinases, we have developed a Machine learning-based approach for Mycobacterium tuberculosis serine/threonine protein kinases using kinase-peptide structure-sequence data. This approach utilizes features derived from kinase three-dimensional-structure environment and known phosphosite sequences to generate support vector machine (SVM)-based kinase-specific predictions of phosphosites of serine/threonine protein kinases (STPKs) with no or scarce data of their substrates. SVM outperformed the four machine learning algorithms we tried (random forest, logistic regression, SVM, and k-nearest neighbors) with an area under the curve receiver-operating characteristic value of 0.88 on the independent testing dataset and a 10-fold cross-validation accuracy of ~81.6% for the final model. Our predicted phosphosites of M. tuberculosis STPKs form a useful resource for experimental biologists enabling elucidation of STPK mediated posttranslational regulation of important cellular processes.

Metabolomics Studies To Decipher Stress Responses in Mycobacterium smegmatis Point to a Putative Pathway of Methylated Amine Biosynthesis.

Mycobacterium smegmatis, the saprophytic soil mycobacterium, is routinely used as a surrogate system to study the human pathogen Mycobacterium tuberculosis It has also been reported as an opportunistic pathogen in immunocompromised hosts. In addition, it can exist in several ecological setups, thereby suggesting its capacity to adapt to a variety of environmental cues. In this study, we employed untargeted proton nuclear magnetic resonance (1H-NMR)-based metabolomics to identify metabolites and metabolic pathways critical for early adaptive responses to acidic stress, oxidative stress, and nutrient starvation in Mycobacterium smegmatis We identified 31, 20, and 46 metabolites that showed significant changes in levels in response to acidic, oxidative, and nutrient starvation stresses, respectively. Pathway analyses showed significant perturbations in purine-pyrimidine, amino-acid, nicotinate-nicotinamide, and energy metabolism pathways. Besides these, differential levels of intermediary metabolites involved in α-glucan biosynthesis pathway were observed. We also detected high levels of organic osmolytes, methylamine, and betaine during nutrient starvation and oxidative stress. Further, tracing the differential levels of these osmolytes through computational search tools, gene expression studies (using reverse transcription-PCR [RT-PCR]), and enzyme assays, we detected the presence of a putative pathway of biosynthesis of betaine, methylamine, and dimethylamine previously unreported in Mycobacterium smegmatisIMPORTANCE Alterations in metabolite levels provide fast and direct means to regulate enzymatic reactions and, therefore, metabolic pathways. This study documents, for the first time, the metabolic changes that occur in Mycobacterium smegmatis as a response to three stresses, namely, acidic stress, oxidative stress, and nutrient starvation. These stresses are also faced by intracellular mycobacteria during infection and therefore may be extended to frame therapeutic interventions for pathogenic mycobacteria. In addition to the purine-pyrimidine, amino acid, nicotinate-nicotinamide, and energy metabolism pathways that were found to be affected in response to different stresses, a novel putative methylamine biosynthesis pathway was identified to be present in Mycobacterium smegmatis.

Using structural knowledge in the protein data bank to inform the search for potential host-microbe protein interactions in sequence space: application to Mycobacterium tuberculosis.

BACKGROUND: A comprehensive map of the human-M. tuberculosis (MTB) protein interactome would help fill the gaps in our understanding of the disease, and computational prediction can aid and complement experimental studies towards this end. Several sequence-based in silico approaches tap the existing data on experimentally validated protein-protein interactions (PPIs); these PPIs serve as templates from which novel interactions between pathogen and host are inferred. Such comparative approaches typically make use of local sequence alignment, which, in the absence of structural details about the interfaces mediating the template interactions, could lead to incorrect inferences, particularly when multi-domain proteins are involved. RESULTS: We propose leveraging the domain-domain interaction (DDI) information in PDB complexes to score and prioritize candidate PPIs between host and pathogen proteomes based on targeted sequence-level comparisons. Our method picks out a small set of human-MTB protein pairs as candidates for physical interactions, and the use of functional meta-data suggests that some of them could contribute to the in vivo molecular cross-talk between pathogen and host that regulates the course of the infection. Further, we present numerical data for Pfam domain families that highlights interaction specificity on the domain level. Not every instance of a pair of domains, for which interaction evidence has been found in a few instances (i.e. structures), is likely to functionally interact. Our sorting approach scores candidates according to how "distant" they are in sequence space from known examples of DDIs (templates). Thus, it provides a natural way to deal with the heterogeneity in domain-level interactions. CONCLUSIONS: Our method represents a more informed application of local alignment to the sequence-based search for potential human-microbial interactions that uses available PPI data as a prior. Our approach is somewhat limited in its sensitivity by the restricted size and diversity of the template dataset, but, given the rapid accumulation of solved protein complex structures, its scope and utility are expected to keep steadily improving.

GroEL2 of Mycobacterium tuberculosis Reveals the Importance of Structural Pliability in Chaperonin Function.

UNLABELLED: Intracellular protein folding is mediated by molecular chaperones, the best studied among which are the chaperonins GroEL and GroES. Conformational changes and allosteric transitions between different metastable states are hallmarks of the chaperonin mechanism. These conformational transitions between three structural domains of GroEL are anchored at two hinges. Although hinges are known to be critical for mediating the communication between different domains of GroEL, the relative importance of hinges on GroEL oligomeric assembly, ATPase activity, conformational changes, and functional activity is not fully characterized. We have exploited the inability of Mycobacterium tuberculosis GroEL2 to functionally complement an Escherichia coli groEL mutant to address the importance of hinge residues in the GroEL mechanism. Various chimeras of M. tuberculosis GroEL2 and E. coli GroEL allowed us to understand the role of hinges and dissect the consequences of oligomerization and substrate binding capability on conformational transitions. The present study explains the concomitant conformational changes observed with GroEL hinge variants and is best supported by the normal mode analysis. IMPORTANCE: Conformational changes and allosteric transitions are hallmarks of the chaperonin mechanism. We have exploited the inability of M. tuberculosis GroEL2 to functionally complement a strain of E. coli in which groEL expression is repressed to address the importance of hinges. The significance of conservation at the hinge regions stands out as a prominent feature of the GroEL mechanism in binding to GroES and substrate polypeptides. The hinge residues play a significant role in the chaperonin activity in vivo and in vitro.

Delineation of key regulatory elements identifies points of vulnerability in the mitogen-activated signaling network.

Drug development efforts against cancer are often hampered by the complex properties of signaling networks. Here we combined the results of an RNAi screen targeting the cellular signaling machinery, with graph theoretical analysis to extract the core modules that process both mitogenic and oncogenic signals to drive cell cycle progression. These modules encapsulated mechanisms for coordinating seamless transition of cells through the individual cell cycle stages and, importantly, were functionally conserved across different cancer cell types. Further analysis also enabled extraction of the core signaling axes that progressively guide commitment of cells to the division cycle. Importantly, pharmacological targeting of the least redundant nodes in these axes yielded a synergistic disruption of the cell cycle in a tissue-type-independent manner. Thus, the core elements that regulate temporally distinct stages of the cell cycle provide attractive targets for the development of multi-module-based chemotherapeutic strategies.

Confronting Tuberculosis: A Synthetic Quinoline-Isonicotinic Acid Hydrazide Hybrid Compound as a Potent Lead Molecule Against Mycobacterium tuberculosis.

The current tuberculosis (TB) treatment is challenged by a complex first-line treatment for drug-sensitive (DS) TB. Additionally, the prevalence of multidrug (MDR)- and extensively drug (XDR)-resistant TB necessitates the search for new drug prototypes. We synthesized and screened 30 hybrid compounds containing aminopyridine and 2-chloro-3-formyl quinoline to arrive at a compound with potent antimycobacterial activity, UH-NIP-16. Subsequently, antimycobacterial activity against DS and MDR Mycobacterium tuberculosis (M.tb) strains were performed. It demonstrated an MIC50 value of 1.86 ± 0.21 µM for laboratory pathogenic M.tb strain H37Rv and 3.045 ± 0.813 µM for a clinical M.tb strain CDC1551. UH-NIP-16 also decreased the MIC50 values of streptomycin, isoniazid, ethambutol, and bedaquiline to about 45, 55, 68, and 76%, respectively, when used in combination, potentiating their activities. The molecule was active against a clinical MDR M.tb strain. Cytotoxicity on PBMCs from healthy donors and on human cell lines was found to be negligible. Further, blind docking of UH-NIP-16 using Auto Dock Vina and MGL tools onto diverse M.tb proteins showed high binding affinities with multiple M.tb proteins, the top five targets being metabolically critical proteins CelA1, DevS, MmaA4, lysine acetyltransferase, and immunity factor for tuberculosis necrotizing toxin. These bindings were confirmed by fluorescence spectroscopy using a representative protein, MmaA4. Envisaging that a pathogen will have a lower probability of developing resistance to a hybrid molecule with multiple targets, we propose that UH-NIP-16 can be further developed as a lead molecule with the bacteriostatic potential against M.tb, both alone and in combination with first-line drugs.

The crystal structure of Mycobacterium tuberculosis NrdH at 0.87 Å suggests a possible mode of its activity.

Members of the NrdH family of redox proteins, which consists of small glutaredoxin-like proteins with thioredoxin-like activity, serve as the reducing partners of class Ib ribonucleotide reductases. Here, we report the crystal structure of NrdH from Mycobacterium tuberculosis, refined to a crystallographic R factor of 14.02% (Rfree = 15.53%) at 0.87 Å resolution. The tertiary structure of M. tuberculosis NrdH has a typical thioredoxin fold as expected. The extremely high resolution of the structure allows us to dissect the functionality of the protein in great depth. Structural superimposition of M. tuberculosis NrdH and thioredoxin reductase over the Escherichia coli thioredoxin reductase-thioredoxin complex suggests the ability of NrdH to accept electrons from M. tuberculosis thioredoxin reductase. This raises the important question of why glutaredoxins are unable to accept electrons from thioredoxin reductases and why thioredoxins are unable to reduce ribonucleotide reductases. Furthermore, forms of NrdH from other organisms have been shown to be a specific reductant of class Ib ribonucleotide reductases. We attempt to explain this substrate specificity by modeling the C-terminal peptide of a ribunucleotide subunit, NrdE, in the active site of NrdH using the already available Grx-NrdA-Cter-peptide structure. Statistical coupling analysis of NrdH, glutaredoxins, and thioredoxins reveals different sets of co-evolving contiguous clusters of amino acid residues, which might explain the differences in the biochemical properties of these structurally similar yet functionally distinct subclasses of proteins.

The role of thioredoxin proteins in Mycobacterium tuberculosis probed by proteome-wide target profiling.

Mycobacterium tuberculosis encounters diverse microenvironments, including oxidative assault (ROS and RNS), when it attempts to establish itself within its human host. Therefore, redox sensory and regulation processes are assumed significant importance, as these are essential processes for M. tuberculosis to survive under these hostile conditions. M. tuberculosis contains thioredoxin system to maintain redox homeostasis, which establish a balance between the thiol/dithiol couple. Still very less is known about it. In the present study, we attempted to capture the targets of all the M. tuberculosis thioredoxin proteins (viz., TrxB and TrxC) and a thioredoxin-like protein, NrdH, under aerobic and hypoxic conditions by performing thioredoxin trapping chromatography followed by mass spectrometry. We found that TrxC captured the maximum number of targets in both the physiological conditions and most of the targets of TrxB and NrdH showing overlap with targets of TrxC, indicating that TrxC acts as main thioredoxin. Further the PANTHER classification system provides involvement of targets in various metabolic processes and Gene Ontology analysis suggests that glutamine biosynthetic process and Fe-S cluster biosynthesis are the most enriched processes in the target list of TrxC and TrxB respectively. Also, we suggest that the thioredoxin system might play an important role under hypoxia by targeting those proteins which are responsible to sense and maintain hypoxic conditions. Furthermore, our studies establish a link between TrxB and iron-sulfur cluster biogenesis in M. tuberculosis. Ultimately, these findings open a new direction to target the thioredoxin system for screening new anti-mycobacterial drug targets.

Staufen-2 functions as a cofactor for enhanced Rev-mediated nucleocytoplasmic trafficking of HIV-1 genomic RNA via the CRM1 pathway.

Nucleocytoplasmic shuttling of viral elements, supported by several host factors, is essential for the replication of the human immunodeficiency virus (HIV). HIV-1 uses a nuclear RNA export pathway mediated by viral protein Rev to transport its Rev response element (RRE)-containing partially spliced and unspliced transcripts aided by the host nuclear RNA export protein CRM1. The factor(s) interacting with the CRM1-Rev complex are potential antiretroviral target(s) and could serve as a retroviral model system to study nuclear export machinery adapted by these viruses. We earlier reported that cellular Staufen-2 interacts with Rev, facilitating viral-RNA export. Here, we identified the formation of a complex between Staufen-2, CRM1 and Rev. Molecular docking and simulations mapped the interacting residues in the RNA-binding Domain 4 of Staufen-2 as R336 and R337, which were experimentally verified to be critical for interactions among Staufen-2, CRM1 and Rev by mutational analysis. Staufen-2 mutants defective in interaction with CRM1 or Rev failed to supplement the Rev-RNA export activity and viral production, demonstrating the importance of these interactions. Rev-dependent reporter assays and proviral DNA-construct transfection-based studies in Staufen-2 knockout cells in the presence of leptomycin-B (LMB) revealed a significant reduction in CRM1-mediated Rev-dependent RNA export with decreased virus production as compared to Staufen-2 knockout background or LMB treatment alone, suggesting the relevance of these interactions in augmenting RNA export activity of Rev. Our observations provide further insights into the mechanistic intricacies of unspliced viral-RNA export to the cytoplasm and support the notion that abrogating such interactions can reduce HIV-1 proliferation.

The PPE18 of Mycobacterium tuberculosis interacts with TLR2 and activates IL-10 induction in macrophage.

The pathophysiological functions of proline-glutamic acid (PE)/proline-proline-glutamic acid (PPE) family of proteins of Mycobacterium tuberculosis are not well understood. In this study, we demonstrate that one of the PPE proteins, PPE18 can stimulate macrophages to secrete IL-10, known to favor a Th2 type response. The recombinant PPE18 was found to specifically interact with the TLR2 leading to an early and sustained activation of p38 MAPK, which is critical for IL-10 induction. In silico docking analyses and mutation experiments indicate that PPE18 specifically interacts with the leucine rich repeat 11 approximately 15 domain of TLR2 and the site of interaction is different from that of a synthetic lipopeptide Pam(3)CSK(4) known to activate predominantly ERK 1/2. When PMA-differentiated THP-1 macrophages were infected with a mutant Mycobacterium tuberculosis strain lacking the PPE18, produced poorer levels of IL-10 as compared with those infected with the wild-type strain. In contrast, an M. smegmatis strain overexpressing the PPE18 induced higher levels of IL-10 in infected macrophages. Our data indicate that the PPE18 protein may trigger an anti-inflammatory response by inducing IL-10 production.

A novel nucleoid-associated protein of Mycobacterium tuberculosis is a sequence homolog of GroEL.

The Mycobacterium tuberculosis genome sequence reveals remarkable absence of many nucleoid-associated proteins (NAPs), such as HNS, Hfq or DPS. In order to characterize the nucleoids of M. tuberculosis, we have attempted to identify NAPs, and report an interesting finding that a chaperonin-homolog, GroEL1, is nucleoid associated. We report that M. tuberculosis GroEL1 binds DNA with low specificity but high affinity, suggesting that it might have naturally evolved to bind DNA. We are able to demonstrate that GroEL1 can effectively function as a DNA-protecting agent against DNase I or hydroxyl-radicals. Moreover, Atomic Force Microscopic studies reveal that GroEL1 can condense a large DNA into a compact structure. We also provide in vivo evidences that include presence of GroEL1 in purified nucleoids, in vivo crosslinking followed by Southern hybridizations and immunofluorescence imaging in M. tuberculosis confirming that GroEL1: DNA interactions occur in natural biological settings. These findings therefore reveal that M. tuberculosis GroEL1 has evolved to be associated with nucleoids.

Chaperonin Abundance Enhances Bacterial Fitness.

The ability of chaperonins to buffer mutations that affect protein folding pathways suggests that their abundance should be evolutionarily advantageous. Here, we investigate the effect of chaperonin overproduction on cellular fitness in Escherichia coli. We demonstrate that chaperonin abundance confers 1) an ability to tolerate higher temperatures, 2) improved cellular fitness, and 3) enhanced folding of metabolic enzymes, which is expected to lead to enhanced energy harvesting potential.

Structural aspects of lesional and non-lesional skin microbiota reveal key community changes in leprosy patients from India.

Although skin is the primary affected organ in Leprosy, the role of the skin microbiome in its pathogenesis is not well understood. Recent reports have shown that skin of leprosy patients (LP) harbours perturbed microbiota which grants inflammation and disease progression. Herein, we present the results of nested Polymerase Chain Reaction-Denaturing Gradient Gel Electrophoresis (PCR-DGGE) which was initially performed for investigating the diversity of bacterial communities from lesional skin (LS) and non-lesional skin (NLS) sites of LP (n = 11). Further, we performed comprehensive analysis of 16S rRNA profiles corresponding to skin samples from participants (n = 90) located in two geographical locations i.e. Hyderabad and Miraj in India. The genus Staphylococcus was observed to be one of the representative bacteria characterizing healthy controls (HC; n = 30), which in contrast was underrepresented in skin microbiota of LP. Taxa affiliated to phyla Firmicutes and Proteobacteria were found to be signatures of HC and LS, respectively. Observed diversity level changes, shifts in core microbiota, and community network structure support the evident dysbiosis in normal skin microbiota due to leprosy. Insights obtained indicate the need for exploring skin microbiota modulation as a potential therapeutic option for leprosy.

Mapping conformational transitions in cyclic AMP receptor protein: crystal structure and normal-mode analysis of Mycobacterium tuberculosis apo-cAMP receptor protein.

Cyclic AMP (cAMP) receptor protein, which acts as the sensor of cAMP levels in cells, is a well-studied transcription factor that is best known for allosteric changes effected by the binding of cAMP. Although genetic and biochemical data on the protein are available from several sources, structural information about the cAMP-free protein has been lacking. Therefore, the precise atomic events that take place upon binding of cAMP, leading to conformational changes in the protein and its activation to bind DNA, have been elusive. In this work we solved the cAMP-free crystal structure of the Mycobacterium tuberculosis homolog of cAMP receptor protein at 2.9 A resolution, and carried out normal-mode analysis to map conformational transitions among its various conformational states. In our structure, the cAMP-binding domain holds onto the DNA-binding domain via strong hydrophobic interactions, thereby freezing the latter in a conformation that is not competent to bind DNA. The two domains release each other in the presence of cAMP, making the DNA-binding domain more flexible and allowing it to bind its cognate DNA via an induced-fit mechanism. The structure of the cAMP-free protein and results of the normal-mode analysis therefore highlight an elegant mechanism of the allosteric changes effected by the binding of cAMP.

CoRNeA: A Pipeline to Decrypt the Inter-Protein Interfaces from Amino Acid Sequence Information.

Decrypting the interface residues of the protein complexes provides insight into the functions of the proteins and, hence, the overall cellular machinery. Computational methods have been devised in the past to predict the interface residues using amino acid sequence information, but all these methods have been majorly applied to predict for prokaryotic protein complexes. Since the composition and rate of evolution of the primary sequence is different between prokaryotes and eukaryotes, it is important to develop a method specifically for eukaryotic complexes. Here, we report a new hybrid pipeline for predicting the protein-protein interaction interfaces in a pairwise manner from the amino acid sequence information of the interacting proteins. It is based on the framework of Co-evolution, machine learning (Random Forest), and Network Analysis named CoRNeA trained specifically on eukaryotic protein complexes. We use Co-evolution, physicochemical properties, and contact potential as major group of features to train the Random Forest classifier. We also incorporate the intra-contact information of the individual proteins to eliminate false positives from the predictions keeping in mind that the amino acid sequence of a protein also holds information for its own folding and not only the interface propensities. Our prediction on example datasets shows that CoRNeA not only enhances the prediction of true interface residues but also reduces false positive rates significantly.

A novel function of Mycobacterium tuberculosis chaperonin paralog GroEL1 in copper homeostasis.

Among the two GroEL paralogs in Mycobacterium tuberculosis, GroEL1 and GroEL2, GroEL1 has a characteristic histidine-rich C terminus. Since histidine richness is likely to be involved in metal binding, we attempted to decipher the role of GroEL1 in chelating metals and the consequence on M. tuberculosis physiology. Isothermal titration calorimetry showed that GroEL1 binds copper and other metals. Mycobacterial viability assay, redox balance, and DNA protection assay concluded that GroEL1 protects from copper stress in vitro. Solution X-ray scattering and constrained modeling of GroEL1 -/+ copper ions showed reorientation of the apical domain as seen in functional assembly. We conclude that the duplication of chaperonin genes in M. tuberculosis might have led to their evolutionary divergence and consequent functional divergence of chaperonins.

Facilitated oligomerization of mycobacterial GroEL: evidence for phosphorylation-mediated oligomerization.

The distinctive feature of the GroES-GroEL chaperonin system in mediating protein folding lies in its ability to exist in a tetradecameric state, form a central cavity, and encapsulate the substrate via the GroES lid. However, recombinant GroELs of Mycobacterium tuberculosis are unable to act as effective molecular chaperones when expressed in Escherichia coli. We demonstrate here that the inability of M. tuberculosis GroEL1 to act as a functional chaperone in E. coli can be alleviated by facilitated oligomerization. The results of directed evolution involving random DNA shuffling of the genes encoding M. tuberculosis GroEL homologues followed by selection for functional entities suggested that the loss of chaperoning ability of the recombinant mycobacterial GroEL1 and GroEL2 in E. coli might be due to their inability to form canonical tetradecamers. This was confirmed by the results of domain-swapping experiments that generated M. tuberculosis-E. coli chimeras bearing mutually exchanged equatorial domains, which revealed that E. coli GroEL loses its chaperonin activity due to alteration of its oligomerization capabilities and vice versa for M. tuberculosis GroEL1. Furthermore, studying the oligomerization status of native GroEL1 from cell lysates of M. tuberculosis revealed that it exists in multiple oligomeric forms, including single-ring and double-ring variants. Immunochemical and mass spectrometric studies of the native M. tuberculosis GroEL1 revealed that the tetradecameric form is phosphorylated on serine-393, while the heptameric form is not, indicating that the switch between the single- and double-ring variants is mediated by phosphorylation.

Mycobacterium tuberculosis heat shock protein 60 modulates immune response to PPD by manipulating the surface expression of TLR2 on macrophages.

The T-helper (Th) 1 T-cell response to purified protein derivative (PPD) is known to be suppressed in tuberculosis patients which favours intracellular survival of the bacilli. We demonstrate that the Mycobacterium tuberculosis heat shock protein 60 (Mtbhsp60) plays an important role to skew the anti-PPD T-cell response towards the Th2 type when macrophages were used as antigen presenting cells. We found that the PPD-induced IL-12 p40 was downregulated in macrophages by Mtbhsp60. The Mtbhsp60 preferentially induced Toll-like receptor (TLR) 2 without affecting TLR4 expression on macrophages. Interaction of Mtbhsp60 with TLR2 resulted in significant suppression of nuclear c-rel and consequently IL-12 p40 levels in PPD-activated macrophages. Our findings reveal a unique role of the Mtbhsp60 favouring development of Th2 type response by upregulating surface expression of TLR2 on macrophages which could be a survival strategy adopted by the bacilli.

Dynamic changes in protein functional linkage networks revealed by integration with gene expression data.

Response of cells to changing environmental conditions is governed by the dynamics of intricate biomolecular interactions. It may be reasonable to assume, proteins being the dominant macromolecules that carry out routine cellular functions, that understanding the dynamics of protein:protein interactions might yield useful insights into the cellular responses. The large-scale protein interaction data sets are, however, unable to capture the changes in the profile of protein:protein interactions. In order to understand how these interactions change dynamically, we have constructed conditional protein linkages for Escherichia coli by integrating functional linkages and gene expression information. As a case study, we have chosen to analyze UV exposure in wild-type and SOS deficient E. coli at 20 minutes post irradiation. The conditional networks exhibit similar topological properties. Although the global topological properties of the networks are similar, many subtle local changes are observed, which are suggestive of the cellular response to the perturbations. Some such changes correspond to differences in the path lengths among the nodes of carbohydrate metabolism correlating with its loss in efficiency in the UV treated cells. Similarly, expression of hubs under unique conditions reflects the importance of these genes. Various centrality measures applied to the networks indicate increased importance for replication, repair, and other stress proteins for the cells under UV treatment, as anticipated. We thus propose a novel approach for studying an organism at the systems level by integrating genome-wide functional linkages and the gene expression data.