Response of Saccharomyces cerevisiae to severe osmotic stress: evidence for a novel activation mechanism of the HOG MAP kinase pathway

CsPOM1, a DYRK Family Kinase, Plays Diverse Roles in Fungal Development, Virulence, and Stress Tolerance in the Anthracnose Pathogen Colletotrichum scovillei

Colletotrichum scovillei is the major anthracnose fungus of sweet pepper and chili pepper (Capsicum annuum L.), causing significant losses in the yield and quality of the pepper fruits. Molecular mechanisms governing development and pathogenicity have been widely studied in many foliar fungal pathogens, but the information on fruit diseases is still limited. In this study, we determined the functional roles of the dual-specificity tyrosine phosphorylation-regulated kinase CsPOM1 in C. scovillei. Knockout mutant for CsPOM1 gene was obtained via homology-dependent gene replacement. The ΔCspom1 mutant exhibited a reduction in vegetative growth on osmotic stress, surface hydrophobicity, and conidiation compared with wild-type. Conidia of the ΔCspom1 mutant were already two-celled before inoculation on an induction surface, indicating that CsPOM1 negatively regulates conidial cell division. The ΔCspom1 mutant, similar to wild-type, formed appressoria on the plant surface, but was significantly reduced on hydrophobic coverslips, probably due to a defect in the recognition of surface hydrophobicity. Treatment of conidia with cutin monomers restored appressorium formation on hydrophobic coverslips in the ΔCspom1 mutant. On pepper fruits, the ΔCspom1 mutant exhibited delayed penetration and invasive growth, leading to significantly reduced virulence. Collectively, the results showed that CsPOM1 is important for stress tolerance, conidiation, surface hydrophobicity, appressorium formation, and virulence in C. scovillei.

The high osmolarity glycerol (HOG) pathway in fungi†

While fungi are widely occupying nature, many species are responsible for devastating mycosis in humans. Such niche diversity explains how quick fungal adaptation is necessary to endow the capacity of withstanding fluctuating environments and to cope with host-imposed conditions. Among all the molecular mechanisms evolved by fungi, the most studied one is the activation of the phosphorelay signalling pathways, of which the high osmolarity glycerol (HOG) pathway constitutes one of the key molecular apparatus underpinning fungal adaptation and virulence. In this review, we summarize the seminal knowledge of the HOG pathway with its more recent developments. We specifically described the HOG-mediated stress adaptation, with a particular focus on osmotic and oxidative stress, and point out some lags in our understanding of its involvement in the virulence of pathogenic species including, the medically important fungi Candida albicans, Cryptococcus neoformans, and Aspergillus fumigatus, compared to the model yeast Saccharomyces cerevisiae. Finally, we also highlighted some possible applications of the HOG pathway modifications to improve the fungal-based production of natural products in the industry.

Molecular cloning and sequence analysis of a mitogen-activated protein kinase gene in the Antarctic yeast Rhodotorula mucilaginosa AN5

BACKGROUND

The mitogen-activated protein kinase (MAPK) cascades play important roles in various signaling transduction networks of biotic and abiotic stress responses. However, MAPK signaling pathways in cold-active yeast Rhodotorula mucilaginosa have not been reported comprehensively.

METHODS AND RESULTS

In the present study, MAPK gene (RmMAPK) was first cloned and characterized from Antarctic sea ice yeast R. mucilaginosa AN5. The full length of the RmMAPK gene is 1086 bp and encodes a 361 amino acids protein with a predicted molecular mass of 40.9 kDa and a pI of 5.25. The RmMAPK contains 11 MAPK conserved subdomains and the phosphorylation motif TGY located in the activation loop of the kinase. Quantitative real-time PCR and western blot assay revealed that the expression and phosphorylation level of RmMAPK up-regulated rapidly and significantly when yeast cells were subjected to low temperature (4 °C), high salinity (120‰ NaCl) and heavy metal (2 mmol/L CuCl₂).

CONCLUSIONS

All data suggested that the MAPK cascades might act as a key function in response to extreme stresses, such as low temperature, high salinity and heavy metal.

Histidine Kinase Sln1 and cAMP/PKA Signaling Pathways Antagonistically Regulate Sporisorium scitamineum Mating and Virulence via Transcription Factor Prf1

Many prokaryotes and eukaryotes utilize two-component signaling pathways to counter environmental stress and regulate virulence genes associated with infection. In this study, we identified and characterized a conserved histidine kinase (SsSln1), which is the sensor of the two-component system of Sln1-Ypd1-Ssk1 in Sporisorium scitamineum. SsSln1 null mutant exhibited enhanced mating and virulence capabilities in S. scitamineum, which is opposite to what has been reported in Candida albicans. Further investigations revealed that the deletion of SsSLN1 enhanced SsHog1 phosphorylation and nuclear localization and thus promoted S. scitamineum mating. Interestingly, SsSln1 and cAMP/PKA signaling pathways antagonistically regulated the transcription of pheromone-responsive transcription factor SsPrf1, for regulating S. scitamineum mating and virulence. In short, the study depicts a novel mechanism in which the cross-talk between SsSln1 and cAMP/PKA pathways antagonistically regulates mating and virulence by balancing the transcription of the SsPRF1 gene in S. scitamineum.

A Target-Based In Vivo Test System to Identify Novel Fungicides with Mode of Action in the HOG Pathway

Resistance management plays a key role in modern plant protection. There is a growing need to identify new fungicide targets and new modes of action. In this context, it is also mandatory to find new compounds acting on successful target locations. For the latter, so-called target-site-specific test systems emerged to search for inhibitors. Most of them are based on in vitro assays, in which interaction between a compound and a purified target protein is demonstrated. Consequently, getting essential information about potentially toxic effects in the living cell or in the whole organism is not possible. Thus, we present a fluorescent-labelled mutant strain of the rice blast fungus Magnaporthe oryzae as a rapid tool for fluorescence-based identification and visualization of fungicides in vivo with the mode of action in the high osmolarity glycerol (HOG)-signaling pathway. The HOG pathway represents an excellent target for antifungal agents such as the phenylpyrrole fungicides, since almost no relevant resistances have occurred to date, despite 30 years of extensive usage of this fungicide class.

Osmolyte Signatures for the Protection of Aspergillus sydowii Cells under Halophilic Conditions and Osmotic Shock

Aspergillus sydowii is a moderate halophile fungus extensively studied for its biotechnological potential and halophile responses, which has also been reported as a coral reef pathogen. In a recent publication, the transcriptomic analysis of this fungus, when growing on wheat straw, showed that genes related to cell wall modification and cation transporters were upregulated under hypersaline conditions but not under 0.5 M NaCl, the optimal salinity for growth in this strain. This led us to study osmolyte accumulation as a mechanism to withstand moderate salinity. In this work, we show that A. sydowii accumulates trehalose, arabitol, mannitol, and glycerol with different temporal dynamics, which depend on whether the fungus is exposed to hypo- or hyperosmotic stress. The transcripts coding for enzymes responsible for polyalcohol synthesis were regulated in a stress-dependent manner. Interestingly, A. sydowii contains three homologs (Hog1, Hog2 and MpkC) of the Hog1 MAPK, the master regulator of hyperosmotic stress response in S. cerevisiae and other fungi. We show a differential regulation of these MAPKs under different salinity conditions, including sustained basal Hog1/Hog2 phosphorylation levels in the absence of NaCl or in the presence of 2.0 M NaCl, in contrast to what is observed in S. cerevisiae. These findings indicate that halophilic fungi such as A. sydowii utilize different osmoadaptation mechanisms to hypersaline conditions.

MAP kinase Hog1 mediates a cytochrome P450 oxidoreductase to promote the Sporisorium scitamineum cell survival under oxidative stress

The MAP kinase high osmolarity glycerol 1 (Hog1) plays a central role in responding to external oxidative stress in budding yeast Saccchromyces cerevisiae. However, the downstream responsive elements regulated by Hog1 remain poorly understood. In this study, we report that a Sporisorium scitamineum orthologue of Hog1, named as SsHog1, induced transcriptional expression of a putative cytochrome P450 oxidoreductase encoding gene SsCPR1, to antagonize oxidative stress. We found that upon exposure to hydrogen peroxide (H₂ O₂ ), SsHog1 underwent strikingly phosphorylation, which was proved to be critical for transcriptional induction of SsCPR1. Loss of SsCPR1 led to hypersensitive to oxidative stress similar as the sshog1Δ mutant did, but was resistant to osmotic stress, which is different from the sshog1Δ mutant. On the other hand, overexpression of SsCPR1 in the sshog1Δ mutant could partially restore its ability of oxidative stress tolerance, which indicated that the Hog1 MAP kinase regulates the oxidative stress response specifically through cytochrome P450 (SsCpr1) pathway. Overall, our findings highlight a novel MAPK signalling pathway mediated by Hog1 in regulation of the oxidative stress response via the cytochrome P450 system, which plays an important role in host-fungus interaction.

Contribution of the mitogen-activated protein kinase Hog1 to the halotolerance of the marine yeast Debaryomyces hansenii

Halotolerant species are adapted to dealing continually with hyperosmotic environments, having evolved strategies that are uncommon in other organisms. The HOG pathway is the master system that regulates the cellular adaptation under these conditions; nevertheless, apart from the importance of Debaryomyces hansenii as an organism representative of the halotolerant class, its HOG1 pathway has been poorly studied, due to the difficulty of applying conventional recombinant DNA technology. Here we describe for the first time the phenotypic characterisation of a null HOG1 mutant of D. hansenii. Dhhog1Δ strain was found moderately resistant to 1 M NaCl and sensitive to higher concentrations. Under hyperosmotic shock, DhHog1 fully upregulated transcription of DhSTL1 and partially upregulated that of DhGPD1. High osmotic stress lead to long-term inner glycerol accumulation that was partially dependent on DhHog1. These observations indicated that the HOG pathway is required for survival under high external osmolarity but dispensable under low and mid-osmotic conditions. It was also found that DhHog1 can regulate response to alkali stress during hyperosmotic conditions and that it plays a role in oxidative and endoplasmic reticulum stress. Taken together, these results provide new insight into the contribution of this MAPK in halotolerance of this yeast.

Osmostress enhances activating phosphorylation of Hog1 MAP kinase by mono-phosphorylated Pbs2 MAP2K

The MAP kinase (MAPK) Hog1 is the central regulator of osmoadaptation in yeast. When cells are exposed to high osmolarity, the functionally redundant Sho1 and Sln1 osmosensors, respectively, activate the Ste11‐Pbs2‐Hog1 MAPK cascade and the Ssk2/Ssk22‐Pbs2‐Hog1 MAPK cascade. In a canonical MAPK cascade, a MAPK kinase kinase (MAP3K) activates a MAPK kinase (MAP2K) by phosphorylating two conserved Ser/Thr residues in the activation loop. Here, we report that the MAP3K Ste11 phosphorylates only one activating phosphorylation site (Thr‐518) in Pbs2, whereas the MAP3Ks Ssk2/Ssk22 can phosphorylate both Ser‐514 and Thr‐518 under optimal osmostress conditions. Mono‐phosphorylated Pbs2 cannot phosphorylate Hog1 unless the reaction between Pbs2 and Hog1 is enhanced by osmostress. The lack of the osmotic enhancement of the Pbs2‐Hog1 reaction suppresses Hog1 activation by basal MAP3K activities and prevents pheromone‐to‐Hog1 crosstalk in the absence of osmostress. We also report that the rapid‐and‐transient Hog1 activation kinetics at mildly high osmolarities and the slow and prolonged activation kinetics at severely high osmolarities are both caused by a common feedback mechanism.

Visualizing fungicide action: an in vivo tool for rapid validation of fungicides with target location HOG pathway

BACKGROUND

The mitogen-activated protein kinase MoHog1p was fused with a green fluorescent protein (GFP) in the filamentous fungus Magnaporthe oryzae. The MoHOG1::GFP mutant was found to be an excellent tool visualizing in vivo fungicide-dependent translocation of MoHog1p into the nucleus. Validation of pathway specificity was achieved by generating fluorescence-labelled MoHog1p in the ΔMohik1 'loss of function' mutant strain.

RESULTS

GFP-labelled MoHog1p expressed in the wildtype and in ΔMohik1 demonstrates that fludioxonil is acting on the HOG pathway and even more precisely that fungicide action is dependent on the group III histidine kinase MoHik1p. GFP-tagged MoHog1p translocated into the nucleus upon fungicide treatment in the MoHOG1::GFP mutant within seconds, but did not do so in the ΔMohik1/HOG1::GFP mutant.

CONCLUSION

Here, we developed a rapid in vivo tool for fluorescent-based validation of fungicides targeting the HOG-signaling pathway. Furthermore, using the fluorescent mutants generated in this study, we are able to visualize that fungicide action is dependent on the histidine kinase MoHik1p but operates in a different mechanism of pathway activation compared to osmotic stress. © 2018 Society of Chemical Industry.

Sho1 and Msb2 Play Complementary but Distinct Roles in Stress Responses, Sexual Differentiation, and Pathogenicity of Cryptococcus neoformans

The high-osmolarity glycerol response (HOG) pathway is pivotal in environmental stress response, differentiation, and virulence of Cryptococcus neoformans, which causes fatal meningoencephalitis. A putative membrane sensor protein, Sho1, has been postulated to regulate HOG pathway, but its regulatory mechanism remains elusive. In this study, we characterized the function of Sho1 with relation to the HOG pathway in C. neoformans. Sho1 played minor roles in osmoresistance, thermotolerance, and maintenance of membrane integrity mainly in a HOG-independent manner. However, it was dispensable for cryostress resistance, primarily mediated through the HOG pathway. A mucin-like transmembrane (TM) protein, Msb2, which interacts with Sho1 in Saccharomyces cerevisiae, was identified in C. neoformans, but found not to interact with Sho1. MSB2 codeletion with SHO1 further decreased osmoresistance and membrane integrity, but not thermotolerance, of sho1Δ mutant, indicating that both factors play to some level redundant but also discrete roles in C. neoformans. Sho1 and Msb2 played redundant roles in promoting the filamentous growth in sexual differentiation in a Cpk1-independent manner, in contrast to the inhibitory effect of the HOG pathway in the process. However, both factors contributed independently to Cpk1 phosphorylation during vegetative growth and endoplasmic reticulum (ER) stress response. Finally, Sho1 and Msb2 play distinct but complementary roles in the pulmonary virulence of C. neoformans. Overall, Sho1 and Msb2 play complementary but distinct roles in stress response, differentiation, and pathogenicity of C. neoformans.

The Lsm1-7/Pat1 complex binds to stress-activated mRNAs and modulates the response to hyperosmotic shock

RNA-binding proteins (RBPs) establish the cellular fate of a transcript, but an understanding of these processes has been limited by a lack of identified specific interactions between RNA and protein molecules. Using MS2 RNA tagging, we have purified proteins associated with individual mRNA species induced by osmotic stress, STL1 and GPD1. We found members of the Lsm1-7/Pat1 RBP complex to preferentially bind these mRNAs, relative to the non-stress induced mRNAs, HYP2 and ASH1. To assess the functional importance, we mutated components of the Lsm1-7/Pat1 RBP complex and analyzed the impact on expression of osmostress gene products. We observed a defect in global translation inhibition under osmotic stress in pat1 and lsm1 mutants, which correlated with an abnormally high association of both non-stress and stress-induced mRNAs to translationally active polysomes. Additionally, for stress-induced proteins normally triggered only by moderate or high osmostress, in the mutants the protein levels rose high already at weak hyperosmosis. Analysis of ribosome passage on mRNAs through co-translational decay from the 5’ end (5P-Seq) showed increased ribosome accumulation in lsm1 and pat1 mutants upstream of the start codon. This effect was particularly strong for mRNAs induced under osmostress. Thus, our results indicate that, in addition to its role in degradation, the Lsm1-7/Pat1 complex acts as a selective translational repressor, having stronger effect over the translation initiation of heavily expressed mRNAs. Binding of the Lsm1-7/Pat1p complex to osmostress-induced mRNAs mitigates their translation, suppressing it in conditions of weak or no stress, and avoiding a hyperresponse when triggered.

When confronted with external physical or chemical stress, cells respond by increasing the mRNA output of a small number of genes required for stress survival, while shutting down the majority of other genes. Moreover, each mRNA is regulated under stress to either enhance or diminish its translation into proteins. The overall purpose is for the cell to optimize gene expression for survival and recovery during rapidly changing conditions. Much of this regulation is mediated by RNA-binding proteins. We have isolated proteins binding to specific mRNAs induced by stress, to investigate how they affect the stress response. We found members of one protein complex to be bound to stress-induced mRNAs. When mutants lacking these proteins were exposed to stress, ribosomes were more engaged with translating mRNAs than in the wild-type. In the mutants, it was also possible to trigger expression of stress proteins with only minimal stress levels. Tracing the passage of ribosomes over mRNAs, we saw that ribosomes accumulated around the start codon in the mutants. These findings indicate that the protein complex is required to moderate the stress response and prevent it from overreacting, which would be harmful for the cell.

Activation of the Hog1 MAPK by the Ssk2/Ssk22 MAP3Ks, in the absence of the osmosensors, is not sufficient to trigger osmostress adaptation in Saccharomyces cerevisiae

Yeast cells respond to hyperosmotic stress by activating the high-osmolarity glycerol (HOG) pathway, which consists of two branches, Hkr1/Msb2-Sho1 and Sln1, which trigger phosphorylation and nuclear internalization of the Hog1 mitogen-activated protein kinase. In the nucleus, Hog1 regulates gene transcription and cell cycle progression, which allows the cell to respond and adapt to hyperosmotic conditions. This study demonstrates that the uncoupling of the known sensors of both branches of the pathway at the level of Ssk1 and Ste11 impairs cell growth in hyperosmotic medium. However, under these conditions, Hog1 was still phosphorylated and internalized into the nucleus, suggesting the existence of an alternative Hog1 activation mechanism. In the ssk1ste11 mutant, phosphorylated Hog1 failed to associate with chromatin and to activate transcription of canonical hyperosmolarity-responsive genes. Accordingly, Hog1 also failed to induce glycerol production at the levels of a wild-type strain. Inactivation of the Ptp2 phosphatase moderately rescued growth impairment of the ssk1ste11 mutant under hyperosmotic conditions, indicating that downregulation of the HOG pathway only partially explains the phenotypes displayed by the ssk1ste11 mutant. Cell cycle defects were also observed in response to stress when Hog1 was phosphorylated in the ssk1ste11 mutant. Taken together, these observations indicate that Hog1 phosphorylation by noncanonical upstream mechanisms is not sufficient to trigger a protective response to hyperosmotic stress.

A Slowed Cell Cycle Stabilizes the Budding Yeast Genome

During cell division, aberrant DNA structures are detected by regulators called checkpoints that slow division to allow error correction. In addition to checkpoint-induced delay, it is widely assumed, though rarely shown, that merely slowing the cell cycle might allow more time for error detection and correction, thus resulting in a more stable genome. Fidelity by a slowed cell cycle might be independent of checkpoints. Here we tested the hypothesis that a slowed cell cycle stabilizes the genome, independent of checkpoints, in the budding yeast Saccharomyces cerevisiae We were led to this hypothesis when we identified a gene (ERV14, an ER cargo membrane protein) that when mutated, unexpectedly stabilized the genome, as measured by three different chromosome assays. After extensive studies of pathways rendered dysfunctional in erv14 mutant cells, we are led to the inference that no particular pathway is involved in stabilization, but rather the slowed cell cycle induced by erv14 stabilized the genome. We then demonstrated that, in genetic mutations and chemical treatments unrelated to ERV14, a slowed cell cycle indeed correlates with a more stable genome, even in checkpoint-proficient cells. Data suggest a delay in G2/M may commonly stabilize the genome. We conclude that chromosome errors are more rarely made or are more readily corrected when the cell cycle is slowed (even ∼15 min longer in an ∼100-min cell cycle). And, some chromosome errors may not signal checkpoint-mediated responses, or do not sufficiently signal to allow correction, and their correction benefits from this "time checkpoint."

Transcriptome analysis of Sporisorium scitamineum reveals critical environmental signals for fungal sexual mating and filamentous growth

Background

Sporisorium scitamineum causes the sugarcane smut disease, one of the most serious constraints to global sugarcane production. S. scitamineum possesses a sexual mating system composed of two mating-type loci, a and b locus. We previously identified and deleted the b locus in S. scitamineum, and found that the resultant SsΔMAT-1b mutant was defective in mating and pathogenicity.

Results

To further understand the function of b-mating locus, we carried out transcriptome analysis by comparing the transcripts of the mutant strain SsΔMAT-1b, from which the SsbE1 and SsbW1 homeodomain transcription factors have previously been deleted, with those from the wild-type MAT-1 strain. Also the transcripts from SsΔMAT-1b X MAT-2 were compared with those from wild-type MAT-1 X MAT-2 mating. A total of 209 genes were up-regulated (p < 0.05) in the SsΔMAT-1b mutant, compared to the wild-type MAT-1 strain, while 148 genes down-regulated (p < 0.05). In the mixture, 120 genes were up-regulated (p < 0.05) in SsΔMAT-1b X MAT-2, which failed to mate, compared to the wild-type MAT-1 X MAT-2 mating, and 271 genes down-regulated (p < 0.05). By comparing the up- and down-regulated genes in these two sets, it was found that 15 up-regulated and 37 down-regulated genes were common in non-mating haploid and mating mixture, which indeed could be genes regulated by b-locus. Furthermore, GO and KEGG enrichment analysis suggested that carbon metabolism pathway and stress response mediated by Hog1 MAPK signaling pathway were altered in the non-mating sets.

Conclusions

Experimental validation results indicate that the bE/bW heterodimeric transcriptional factor, encoded by the b-locus, could regulate S. scitamineum sexual mating and/or filamentous growth via modulating glucose metabolism and Hog1-mediating oxidative response.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2691-5) contains supplementary material, which is available to authorized users.

A Thermodynamic Model of Monovalent Cation Homeostasis in the Yeast Saccharomyces cerevisiae

Cationic and heavy metal toxicity is involved in a substantial number of diseases in mammals and crop plants. Therefore, the understanding of tightly regulated transporter activities, as well as conceiving the interplay of regulatory mechanisms, is of substantial interest. A generalized thermodynamic description is developed for the complex interplay of the plasma membrane ion transporters, membrane potential and the consumption of energy for maintaining and restoring specific intracellular cation concentrations. This concept is applied to the homeostasis of cation concentrations in the yeast cells of S. cerevisiae. The thermodynamic approach allows to model passive ion fluxes driven by the electrochemical potential differences, but also primary or secondary active transport processes driven by the inter- play of different ions (symport, antiport) or by ATP consumption (ATPases). The model-confronted with experimental data-reproduces the experimentally observed potassium and proton fluxes induced by the external stimuli KCl and glucose. The estimated phenomenological constants combine kinetic parameters and transport coefficients. These are in good agreement with the biological understanding of the transporters thus providing a better understanding of the control exerted by the coupled fluxes. The model predicts the flux of additional ion species, like e.g. chloride, as a potential candidate for counterbalancing positive charges. Furthermore, the effect of a second KCl stimulus is simulated, predicting a reduced cellular response for cells that were first exposed to a high KCl stimulus compared to cells pretreated with a mild KCl stimulus. By describing the generalized forces that are responsible for a given flow, the model provides information and suggestions for new experiments. Furthermore, it can be extended to other systems such as e.g. Candida albicans, or selected plant cells.

Influence of heat shock and osmotic stresses on the growth and viability of Saccharomyces cerevisiae SUBSC01

Background

With a preceding scrutiny of bacterial cellular responses against heat shock and oxidative stresses, current research further investigated such impact on yeast cell. Present study attempted to observe the influence of high temperature (44–46 °C) on the growth and budding pattern of Saccharomyces cerevisiae SUBSC01. Effect of elevated sugar concentrations as another stress stimulant was also observed. Cell growth was measured through the estimation of the optical density at 600 nm (OD₆₀₀) and by the enumeration of colony forming units on the agar plates up to 450 min.

Results

Subsequent transformation in the yeast morphology and the cellular arrangement were noticed. A delayed and lengthy lag phase was observed when yeast strain was grown at 30, 37, and 40 °C, while at 32.5 °C, optimal growth pattern was noticed. Cells were found to lose culturability completely at 46 °C whereby cells without the cytoplasmic contents were also observed under the light microscope. Thus the critical growth temperature was recorded as 45 °C which was the highest temperature at which S. cerevisiae SUBSC01 could grow. However, a complete growth retardation was observed at 45 °C with the high concentrations of dextrose (0.36 g/l) and sucrose (0.18 g/l). Notably, yeast budding was found at 44 and 45 °C up to 270 min of incubation, which was further noticed to be suppressed at 46 °C.

Conclusions

Present study revealed that the optimal and the critical growth temperatures of S. cerevisiae SUBSC01 were 32.5 and 45 °C, respectively; and also projected on the inhibitory concentrations of sugars on yeast growth at that temperature.

Ineffective Phosphorylation of Mitogen-Activated Protein Kinase Hog1p in Response to High Osmotic Stress in the Yeast Kluyveromyces lactis

When treated with a hyperosmotic stimulus, Kluyveromyces lactis cells respond by activating the mitogen-activated protein kinase (MAPK) K. lactis Hog1 (KlHog1) protein via two conserved branches, SLN1 and SHO1. Mutants affected in only one branch can cope with external hyperosmolarity by activating KlHog1p by phosphorylation, except for single ΔKlste11 and ΔKlste50 mutants, which showed high sensitivity to osmotic stress, even though the other branch (SLN1) was intact. Inactivation of both branches by deletion of KlSHO1 and KlSSK2 also produced sensitivity to high salt. Interestingly, we have observed that in ΔKlste11 and ΔKlsho1 ΔKlssk2 mutants, which exhibit sensitivity to hyperosmotic stress, and contrary to what would be expected, KlHog1p becomes phosphorylated. Additionally, in mutants lacking both MAPK kinase kinases (MAPKKKs) present in K. lactis (KlSte11p and KlSsk2p), the hyperosmotic stress induced the phosphorylation and nuclear internalization of KlHog1p, but it failed to induce the transcriptional expression of KlSTL1 and the cell was unable to grow in high-osmolarity medium. KlHog1p phosphorylation via the canonical HOG pathway or in mutants where the SHO1 and SLN1 branches have been inactivated requires not only the presence of KlPbs2p but also its kinase activity. This indicates that when the SHO1 and SLN1 branches are inactivated, high-osmotic-stress conditions activate an independent input that yields active KlPbs2p, which, in turn, renders KlHog1p phosphorylation ineffective. Finally, we found that KlSte11p can alleviate the sensitivity to hyperosmotic stress displayed by a ΔKlsho1 ΔKlssk2 mutant when it is anchored to the plasma membrane by adding the KlSho1p transmembrane segments, indicating that this chimeric protein can substitute for KlSho1p and KlSsk2p.

Coordinated gene regulation in the initial phase of salt stress adaptation

Stress triggers complex transcriptional responses, which include both gene activation and repression. We used time-resolved reporter assays in living yeast cells to gain insights into the coordination of positive and negative control of gene expression upon salt stress. We found that the repression of "housekeeping" genes coincides with the transient activation of defense genes and that the timing of this expression pattern depends on the severity of the stress. Moreover, we identified mutants that caused an alteration in the kinetics of this transcriptional control. Loss of function of the vacuolar H(+)-ATPase (vma1) or a defect in the biosynthesis of the osmolyte glycerol (gpd1) caused a prolonged repression of housekeeping genes and a delay in gene activation at inducible loci. Both mutants have a defect in the relocation of RNA polymerase II complexes at stress defense genes. Accordingly salt-activated transcription is delayed and less efficient upon partially respiratory growth conditions in which glycerol production is significantly reduced. Furthermore, the loss of Hog1 MAP kinase function aggravates the loss of RNA polymerase II from housekeeping loci, which apparently do not accumulate at inducible genes. Additionally the Def1 RNA polymerase II degradation factor, but not a high pool of nuclear polymerase II complexes, is needed for efficient stress-induced gene activation. The data presented here indicate that the finely tuned transcriptional control upon salt stress is dependent on physiological functions of the cell, such as the intracellular ion balance, the protective accumulation of osmolyte molecules, and the RNA polymerase II turnover.

Pathway connectivity and signaling coordination in the yeast stress-activated signaling network

Stressed cells coordinate a multi-faceted response spanning many levels of physiology. Yet knowledge of the complete stress-activated regulatory network as well as design principles for signal integration remains incomplete. We developed an experimental and computational approach to integrate available protein interaction data with gene fitness contributions, mutant transcriptome profiles, and phospho-proteome changes in cells responding to salt stress, to infer the salt-responsive signaling network in yeast. The inferred subnetwork presented many novel predictions by implicating new regulators, uncovering unrecognized crosstalk between known pathways, and pointing to previously unknown ‘hubs’ of signal integration. We exploited these predictions to show that Cdc14 phosphatase is a central hub in the network and that modification of RNA polymerase II coordinates induction of stress-defense genes with reduction of growth-related transcripts. We find that the orthologous human network is enriched for cancer-causing genes, underscoring the importance of the subnetwork's predictions in understanding stress biology.

Analysis of osmoadaptation system in budding yeast suggests that regulated degradation of glycerol synthesis enzyme is key to near-perfect adaptation

In order to maintain its turgor pressure at a desired homeostatic level, budding yeast, Saccharomyces cerevisiae responds to the external variation of the osmotic pressure by varying its internal osmotic pressure through regulation of synthesis and transport of the intracellular glycerol. Hog1PP (dually phosphorylated Hog1), a final effector in the signalling pathway of the hyper osmotic stress, regulates the glycerol synthesis both at transcriptional and non-transcriptional stages. It is known that for a step-change in salt concentration leading to moderate osmotic shock, Hog1PP activity shows a transient response before it returns to the vicinity of pre-stimulus level. It is believed that an integrating process in a negative feedback loop can be a design strategy to yield such an adaptive response. Several negative feedback loops have been identified in the osmoadaptation system in yeast. However, the precise location of the integrating process in the osmoadaptation system which includes signalling, gene regulation, metabolism and biophysical modules is unclear. To address this issue, we developed a reduced model which captures various experimental observations of the osmoadaptation behaviour of wild type and mutant strains. Dynamic simulations and steady state analysis suggested that known information about the osmoadaptation system of budding yeast does not necessarily give a perfect integrating process through the known feedback loops of Hog1PP. On the other hand, regulation of glycerol synthesising enzyme degradation can result in a near integrating process leading to a near-perfect adaptation.

Osmostress-induced cell volume loss delays yeast Hog1 signaling by limiting diffusion processes and by Hog1-specific effects

Signal transmission progresses via a series of transient protein-protein interactions and protein movements, which require diffusion within a cell packed with different molecules. Yeast Hog1, the effector protein kinase of the High Osmolarity Glycerol pathway, translocates transiently from the cytosol to the nucleus during adaptation to high external osmolarity. We followed the dynamics of osmostress-induced cell volume loss and Hog1 nuclear accumulation upon exposure of cells to different NaCl concentrations. While Hog1 nuclear accumulation peaked within five minutes following mild osmotic shock it was delayed up to six-fold under severe stress. The timing of Hog1 nuclear accumulation correlated with the degree of cell volume loss and the cells capacity to recover. Also the nuclear translocation of Msn2, the transcription factor of the general stress response pathway, is delayed upon severe osmotic stress suggesting a general phenomenon. We show by direct measurements that the general diffusion rate of Hog1 in the cytoplasm as well as its rate of nuclear transport are dramatically reduced following severe volume reduction. However, neither Hog1 phosphorylation nor Msn2 nuclear translocation were as much delayed as Hog1 nuclear translocation. Our data provide direct evidence that signaling slows down during cell volume compression, probably as a consequence of molecular crowding. Hence one purpose of osmotic adaptation is to restore optimal diffusion rates for biochemical and cell biological processes. In addition, there may be mechanisms slowing down especially Hog1 nuclear translocation under severe stress in order to prioritize Hog1 cytosolic targets.

Initiation of the transcriptional response to hyperosmotic shock correlates with the potential for volume recovery

The control of activity and localization of transcription factors is critical for appropriate transcriptional responses. In eukaryotes, signal transduction components such as mitogen-activated protein kinase (MAPK) shuttle into the nucleus to activate transcription. It is not known in detail how different amounts of nuclear MAPK over time affect the transcriptional response. In the present study, we aimed to address this issue by studying the high osmolarity glycerol (HOG) system in Saccharomyces cerevisiae. We employed a conditional osmotic system, which changes the period of the MAPK Hog1 signal independent of the initial stress level. We determined the dynamics of the Hog1 nuclear localization and cell volume by single-cell analysis in well-controlled microfluidics systems and compared the responses with the global transcriptional output of cell populations. We discovered that the onset of the initial transcriptional response correlates with the potential of cells for rapid adaptation; cells that are capable of recovering quickly initiate the transcriptional responses immediately, whereas cells that require longer time to adapt also respond later. This is reflected by Hog1 nuclear localization, Hog1 promoter association and the transcriptional response, but not Hog1 phosphorylation, suggesting that a presently uncharacterized rapid adaptive mechanism precedes the Hog1 nuclear response. Furthermore, we found that the period of Hog1 nuclear residence affects the amplitude of the transcriptional response rather than the spectrum of responsive genes.

Severe osmotic compression triggers a slowdown of intracellular signaling, which can be explained by molecular crowding

Regulation of the cellular volume is fundamental for cell survival and function. Deviations from equilibrium trigger dedicated signaling and transcriptional responses that mediate water homeostasis and volume recovery. Cells are densely packed with proteins, and molecular crowding may play an important role in cellular processes. Indeed, increasing molecular crowding has been shown to modify the kinetics of biochemical reactions in vitro; however, the effects of molecular crowding in living cells are mostly unexplored. Here, we report that, in yeast, a sudden reduction in cellular volume, induced by severe osmotic stress, slows down the dynamics of several signaling cascades, including the stress-response pathways required for osmotic adaptation. We show that increasing osmotic compression decreases protein mobility and can eventually lead to a dramatic stalling of several unrelated signaling and cellular processes. The rate of these cellular processes decreased exponentially with protein density when approaching stalling osmotic compression. This suggests that, under compression, the cytoplasm behaves as a soft colloid undergoing a glass transition. Our results shed light on the physical mechanisms that force cells to cope with volume fluctuations to maintain an optimal protein density compatible with cellular functions.

Ssk1p-independent activation of Ssk2p plays an important role in the osmotic stress response in Saccharomyces cerevisiae: alternative activation of Ssk2p in osmotic stress

In Saccharomyces cerevisiae, external high osmolarity activates the HOG MAPK pathway, which controls various aspects of osmoregulation. MAPKKK Ssk2 is activated by Ssk1 in the SLN1 branch of the osmoregulatory HOG MAPK pathway under hyperosmotic stress. We observed that Ssk2 can be activated independent of Ssk1 upon osmotic shock by an unidentified mechanism. The domain for the Ssk1p-independent activation was identified to be located between the amino acids 177∼240. This region might be involved in the binding of an unknown regulator to Ssk2 which in turn activates Ssk2p without Ssk1p under hyperosmotic stress. The osmotic stress response through the Ssk1p-independent Ssk2p activation is strong, although its duration is short compared with the Ssk1p-dependent activation. The alternative Ssk2p activation is also important for the salt resistance.

Boolean model of yeast apoptosis as a tool to study yeast and human apoptotic regulations

Programmed cell death (PCD) is an essential cellular mechanism that is evolutionary conserved, mediated through various pathways and acts by integrating different stimuli. Many diseases such as neurodegenerative diseases and cancers are found to be caused by, or associated with, regulations in the cell death pathways. Yeast Saccharomyces cerevisiae, is a unicellular eukaryotic organism that shares with human cells components and pathways of the PCD and is therefore used as a model organism. Boolean modeling is becoming promising approach to capture qualitative behavior and describe essential properties of such complex networks. Here we present large literature-based and to our knowledge first Boolean model that combines pathways leading to apoptosis (a type of PCD) in yeast. Analysis of the yeast model confirmed experimental findings of anti-apoptotic role of Bir1p and pro-apoptotic role of Stm1p and revealed activation of the stress protein kinase Hog proposing the maximal level of activation upon heat stress. In addition we extended the yeast model and created an in silico humanized yeast in which human pro- and anti-apoptotic regulators Bcl-2 family and Valosin-contain protein (VCP) are included in the model. We showed that accumulation of Bax in silico humanized yeast shows apoptotic markers and that VCP is essential target of Akt Signaling. The presented Boolean model provides comprehensive description of yeast apoptosis network behavior. Extended model of humanized yeast gives new insights of how complex human disease like neurodegeneration can initially be tested.

Molecular cloning, characterization and differential expression of DRK1 in Sporothrix schenckii

The dimorphism of Sporothrix schenckii (S. schenckii) reflects a developmental switch in morphology and lifestyle that is necessary for virulence. DRK1, a hybrid histidine kinase, functions as a global regulator of dimorphism and virulence in Blastomyces dermatitidis (B. dermatitidis) and Histoplasma capsulatum (H. capsulatum). The partial cDNA sequence of DRK1 of S. schenckii, designated SsDRK1, was obtained using degenerate primers based on the conserved domain of the DRK1 of other fungi. The complete cDNA sequence of SsDRK1 was obtained by 5' and 3' RACE. The full-length cDNA is 4743 bp in size and has an open reading frame (ORF) of 4071 bp, encoding 1356 amino acid residues. The predicted molecular mass of SsDRK1 is 147.3 kDa with an estimated theoretical isoelectric point of 5.46. The deduced amino acid sequence of SsDRK1 shows 65% identity to that of B. dermatitidis. The SsDRK1 was predicted to be a soluble histidine kinase and to contain three parts: sensor domain, linker domain and functional domain. Quantitative real-time RT-PCR revealed that SsDRK1 was more highly expressed in the yeast stage compared with that in the mycelial stage, which indicated that the SsDRK1 may be involved in the dimorphic switch in S. schenckii.

Low temperature highlights the functional role of the cell wall integrity pathway in the regulation of growth in Saccharomyces cerevisiae

Unlike other stresses, the physiological significance and molecular mechanisms involved in the yeast cold response are largely unknown. In the present study, we show that the CWI (cell wall integrity) pathway plays an important role in the growth of Saccharomyces cerevisiae at low temperatures. Cells lacking the Wsc1p (wall integrity and stress response component 1) membrane sensor or the MAPKs (mitogen-activated protein kinases) Bck1p (bypass of C kinase 1), Mkk (Mapk kinase) 1p/Mkk2p or Slt2p (suppressor of lyt2) exhibited cold sensitivity. However, there was no evidence of either a cold-provoked perturbation of the cell wall or a differential cold expression program mediated by Slt2p. The results of the present study suggest that Slt2p is activated by different inputs in response to nutrient signals and mediates growth control through TORC1 (target of rapamycin 1 complex)-Sch9p (suppressor of cdc25) and PKA (protein kinase A) at low temperatures. We found that absence of TOR1 (target of rapamycin 1) causes cold sensitivity, whereas a ras2Δ mutant shows increased cold growth. Lack of Sch9p alleviates the phenotype of slt2Δ and bck1Δ mutant cells, as well as attenuation of PKA activity by overexpression of BCY1 (bypass of cyclase mutations 1). Interestingly, swi4Δ mutant cells display cold sensitivity, but the phenotype is neither mediated by the Slt2p-regulated induction of Swi4p (switching deficient 4)-responsive promoters nor influenced by osmotic stabilization. Hence, cold signalling through the CWI pathway has distinct features and might mediate still unknown effectors and targets.

An integrated pathway system modeling of Saccharomyces cerevisiae HOG pathway: a Petri net based approach

Biochemical networks comprise many diverse components and interactions between them. It has intracellular signaling, metabolic and gene regulatory pathways which are highly integrated and whose responses are elicited by extracellular actions. Previous modeling techniques mostly consider each pathway independently without focusing on the interrelation of these which actually functions as a single system. In this paper, we propose an approach of modeling an integrated pathway using an event-driven modeling tool, i.e., Petri nets (PNs). PNs have the ability to simulate the dynamics of the system with high levels of accuracy. The integrated set of signaling, regulatory and metabolic reactions involved in Saccharomyces cerevisiae's HOG pathway has been collected from the literature. The kinetic parameter values have been used for transition firings. The dynamics of the system has been simulated and the concentrations of major biological species over time have been observed. The phenotypic characteristics of the integrated system have been investigated under two conditions, viz., under the absence and presence of osmotic pressure. The results have been validated favorably with the existing experimental results. We have also compared our study with the study of idFBA (Lee et al., PLoS Comput Biol 4:e1000-e1086, 2008) and pointed out the differences between both studies. We have simulated and monitored concentrations of multiple biological entities over time and also incorporated feedback inhibition by Ptp2 which has not been included in the idFBA study. We have concluded that our study is the first to the best of our knowledge to model signaling, metabolic and regulatory events in an integrated form through PN model framework. This study is useful in computational simulation of system dynamics for integrated pathways as there are growing evidences that the malfunctioning of the interplay among these pathways is associated with disease.

Proteomic analysis of Mrr1p- and Tac1p-associated differential protein expression in azole-resistant clinical isolates of Candida albicans

Azole resistance in Candida albicans is frequently caused by the overexpression of multi-drug efflux pump genes MDR1, CDR1, and CDR2 due to gain-of-function mutations in the zinc cluster transcription factors Mrr1p and Tac1p. In this study, we performed a comparative proteomic analysis to identify proteins whose expression level is influenced by these transcription factors. Both 2-DE and PMF were used to examine the expression profiles of six pairs of matched C. albicans isolates carrying gain-of-function mutations in either MRR1 or TAC1 resulting in the overexpression of either MDR1 or CDR1 and CDR2. Using this approach, 17 differentially expressed proteins were identified in the MDR1-overexpressing isolates, while 14 were identified in the isolates that overexpress CDR1 and CDR2. Furthermore, we found that the expression of many of these proteins was increased in a wild-type strain of C. albicans after the introduction of a gain-of-function allele of MRR1 or TAC1. Moreover, disruption of MRR1 and TAC1 in isolates carrying gain-of-function mutations resulted in decreased expression of these proteins, confirming their regulation by Mrr1p or Tac1p. Several proteins involved in heat shock and carbohydrate metabolism were differentially expressed in all clinical isolate sets, but these proteins were not dependent upon either Tac1p or Mrr1p.

Characterization of the adaptive response and growth upon hyperosmotic shock in Saccharomyces cerevisiae

Molecular and physiological details of osmoadaptation in yeast Saccharomyces cerevisiae are well characterized. It is well known that a cell, upon osmotic shock, delays its growth, produces a compatible solute like glycerol in yeast to maintain the osmotic equilibrium. Many genes are regulated by the hyperosmolarity glycerol (HOG) singling pathway, some of which in turn control the carbon flux in the glycolytic pathway for glycerol synthesis and reduced growth. The whole process of survival of cells under hyperosmotic stress is controlled at multiple levels in signaling and metabolic pathways. To better understand the multi-level regulations in yeast to osmotic shock, a mathematical model is formulated which integrates the growth and the osmoadaptation process. The model included the HOG pathway which consists of Sho1 and Sln1 signaling branches, gene regulation, metabolism and cell growth on glucose and ethanol. Experiments were performed to characterize the effect of various concentrations of salt on the wild-type and mutant strains. The model was able to successfully predict the experimental observations for both the wild-type and mutant strains. Further, the model was used to analyze the effects of various regulatory mechanisms prevalent in the signaling and metabolic pathways which are essential in achieving optimum growth in a saline medium. The analysis demonstrated the relevance of the combined effects of regulation at several points in the signaling and metabolic pathways including activation of GPD1 and GPD2, inhibition of PYK and PDC1, closure of the Fps1 channel, volume effect on the glucose uptake rate, downregulation of ethanol synthesis and upregulation of ALD6 for acetate synthesis. The analysis demonstrated that these combined effects orchestrated the phenomena of adaptation to osmotic stress in yeast.

Mathematical models of cell factories: moving towards the core of industrial biotechnology

Industrial biotechnology involves the utilization of cell factories for the production of fuels and chemicals. Traditionally, the development of highly productive microbial strains has relied on random mutagenesis and screening. The development of predictive mathematical models provides a new paradigm for the rational design of cell factories. Instead of selecting among a set of strains resulting from random mutagenesis, mathematical models allow the researchers to predict in silico the outcomes of different genetic manipulations and engineer new strains by performing gene deletions or additions leading to a higher productivity of the desired chemicals. In this review we aim to summarize the main modelling approaches of biological processes and illustrate the particular applications that they have found in the field of industrial microbiology.

The HOG pathway dictates the short-term translational response after hyperosmotic shock

Cellular responses to environmental changes occur on different levels. We investigated the translational response of yeast cells after mild hyperosmotic shock by isolating mRNA associated with multiple ribosomes (polysomes) followed by array analysis. Globally, recruitment of preexisting mRNAs to ribosomes (translational response) is faster than the transcriptional response. Specific functional groups of mRNAs are recruited to ribosomes without any corresponding increase in total mRNA. Among mRNAs under strong translational up-regulation upon shock, transcripts encoding membrane-bound proteins including hexose transporters were enriched. Similarly, numerous mRNAs encoding cytoplasmic ribosomal proteins run counter to the overall trend of down-regulation and are instead translationally mobilized late in the response. Surprisingly, certain transcriptionally induced mRNAs were excluded from ribosomal association after shock. Importantly, we verify, using constructs with intact 5' and 3' untranslated regions, that the observed changes in polysomal mRNA are reflected in protein levels, including cases with only translational up-regulation. Interestingly, the translational regulation of the most highly osmostress-regulated mRNAs was more strongly dependent on the stress-activated protein kinases Hog1 and Rck2 than the transcriptional regulation. Our results show the importance of translational control for fine tuning of the adaptive responses.

Defects in the secretory pathway and high Ca2+ induce multiple P-bodies

mRNA is sequestered and turned over in cytoplasmic processing bodies (PBs), which are induced by various cellular stresses. Unexpectedly, in Saccharomyces cerevisiae, mutants of the small GTPase Arf1 and various secretory pathway mutants induced a significant increase in PB number, compared with PB induction by starvation or oxidative stress. Exposure of wild-type cells to osmotic stress or high extracellular Ca(2+) mimicked this increase in PB number. Conversely, intracellular Ca(2+)-depletion strongly reduced PB formation in the secretory mutants. In contrast to PB induction through starvation or osmotic stress, PB formation in secretory mutants and by Ca(2+) required the PB components Pat1 and Scd6, and calmodulin, indicating that different stressors act through distinct pathways. Consistent with this hypothesis, when stresses were combined, PB number did not correlate with the strength of the translational block, but rather with the type of stress encountered. Interestingly, independent of the stressor, PBs appear as spheres of approximately 40-100 nm connected to the endoplasmic reticulum (ER), consistent with the idea that translation and silencing/degradation occur in a spatially coordinated manner at the ER. We propose that PB assembly in response to stress occurs at the ER and depends on intracellular signals that regulate PB number.

Interactions among HAMP domain repeats act as an osmosensing molecular switch in group III hybrid histidine kinases from fungi

The members of group III hybrid histidine kinases (HHK) are ubiquitous in fungi. Group III HHK have been implicated to function as osmosensors in the high osmolarity glycerol (HOG) pathway that is essential for fungal survival under high osmolarity stress. Recent literature suggests that group III HHK are also involved in conidia formation, virulence in several filamentous fungi, and are an excellent molecular target for antifungal agents. Thus, group III HHK constitute a very important group of sensor kinases. Structurally, group III HHK are distinct from Sln1p, the osmosensing HHK that regulates the HOG pathway in Saccharomyces cerevisiae. Group III HHK lack any transmembrane domain and typically contain HAMP domain repeats at the N terminus. Until now, it is not clear how group III HHK function as an osmosensor to regulate the HOG pathway. To investigate this, we undertook molecular characterization of DhNIK1, an ortholog from osmotolerant yeast Debaryomyces hansenii. We show here that DhNIK1 could complement sln1 mutation in S. cerevisiae thereby confirming its role as a bona fide osmosensor. We further investigated the role of HAMP domains by deleting them systematically. Our results clearly indicate that the HAMP4 domain is crucial for osmosensing by DhNik1p. Most importantly, we also show that the alternative interaction among the HAMP domains regulates the activity of DhNik1p like an "on-off switch" and thus provides, for the first time, an insight into the molecular mechanism of osmosensing by this group of HHKs.

Toxicity testing in the 21st century: a vision and a strategy

With the release of the landmark report Toxicity Testing in the 21st Century: A Vision and a Strategy, the U.S. National Academy of Sciences, in 2007, precipitated a major change in the way toxicity testing is conducted. It envisions increased efficiency in toxicity testing and decreased animal usage by transitioning from current expensive and lengthy in vivo testing with qualitative endpoints to in vitro toxicity pathway assays on human cells or cell lines using robotic high-throughput screening with mechanistic quantitative parameters. Risk assessment in the exposed human population would focus on avoiding significant perturbations in these toxicity pathways. Computational systems biology models would be implemented to determine the dose-response models of perturbations of pathway function. Extrapolation of in vitro results to in vivo human blood and tissue concentrations would be based on pharmacokinetic models for the given exposure condition. This practice would enhance human relevance of test results, and would cover several test agents, compared to traditional toxicological testing strategies. As all the tools that are necessary to implement the vision are currently available or in an advanced stage of development, the key prerequisites to achieving this paradigm shift are a commitment to change in the scientific community, which could be facilitated by a broad discussion of the vision, and obtaining necessary resources to enhance current knowledge of pathway perturbations and pathway assays in humans and to implement computational systems biology models. Implementation of these strategies would result in a new toxicity testing paradigm firmly based on human biology.