Advertisement

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Synthetic Lethality between Gene Defects Affecting a Single Non-essential Molecular Pathway with Reversible Steps

  • Andrei Zinovyev ,

    Contributed equally to this work with: Andrei Zinovyev, Inna Kuperstein

    [email protected] (AZ); [email protected] (WDH)

    Affiliations Institut Curie, Paris, France, INSERM, U900, Paris, France, Mines ParisTech, Fontainebleau, France

    ⨯
  • Inna Kuperstein ,

    Contributed equally to this work with: Andrei Zinovyev, Inna Kuperstein

    Affiliations Institut Curie, Paris, France, INSERM, U900, Paris, France, Mines ParisTech, Fontainebleau, France

    ⨯
  • Emmanuel Barillot,

    Affiliations Institut Curie, Paris, France, INSERM, U900, Paris, France, Mines ParisTech, Fontainebleau, France

    ⨯
  • Wolf-Dietrich Heyer

    [email protected] (AZ); [email protected] (WDH)

    Affiliation Departments of Microbiology & Molecular Genetics and of Molecular & Cellular Biology, University of California, Davis, Davis, California, United States of America

    ⨯

Abstract

Systematic analysis of synthetic lethality (SL) constitutes a critical tool for systems biology to decipher molecular pathways. The most accepted mechanistic explanation of SL is that the two genes function in parallel, mutually compensatory pathways, known as between-pathway SL. However, recent genome-wide analyses in yeast identified a significant number of within-pathway negative genetic interactions. The molecular mechanisms leading to within-pathway SL are not fully understood. Here, we propose a novel mechanism leading to within-pathway SL involving two genes functioning in a single non-essential pathway. This type of SL termed within-reversible-pathway SL involves reversible pathway steps, catalyzed by different enzymes in the forward and backward directions, and kinetic trapping of a potentially toxic intermediate. Experimental data with recombinational DNA repair genes validate the concept. Mathematical modeling recapitulates the possibility of kinetic trapping and revealed the potential contributions of synthetic, dosage-lethal interactions in such a genetic system as well as the possibility of within-pathway positive masking interactions. Analysis of yeast gene interaction and pathway data suggests broad applicability of this novel concept. These observations extend the canonical interpretation of synthetic-lethal or synthetic-sick interactions with direct implications to reconstruct molecular pathways and improve therapeutic approaches to diseases such as cancer.

Author Summary

Organizing gene functions into molecular pathways is a major challenge in biology. The observation that two viable gene mutations become lethal when combined as a double mutant has been developed into a major genetic tool called synthetic lethality. The classic interpretation of synthetic lethality stipulates that the two mutations identify genes that work in parallel, mutually compensatory pathways that together perform an essential function. However, a significant number of negative interactions are caused by defects affecting a single molecular pathway. Here, we recapitulate by mathematical modeling recent experimental data that demonstrate synthetic lethality between mutations in genes acting in a single, non-essential molecular pathway. We propose a novel mechanism involving reversible pathways steps and trapping of an intermediate. The modeling also predicts that overexpression of certain genes functioning in reversible pathways will lead to synthetic lethality with gene defects in the same pathway. Our results significantly broaden the interpretation of synthetic lethal and synthetic dosage effects, which fundamentally impacts the assignment of genes to pathways. The concept of synthetic lethality has been applied to cancer therapy, and our modeling results suggest new approaches to how to target a single pathway to induce synthetic lethality in cancer cells.

Citation: Zinovyev A, Kuperstein I, Barillot E, Heyer W-D (2013) Synthetic Lethality between Gene Defects Affecting a Single Non-essential Molecular Pathway with Reversible Steps. PLoS Comput Biol 9(4): e1003016. https://doi.org/10.1371/journal.pcbi.1003016

Editor: Frederick P. Roth, University of Toronto and Mount Sinai Hospital, Canada

Received: August 21, 2012; Accepted: February 18, 2013; Published: April 4, 2013

Copyright: © 2013 Zinovyev et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: AZ, IK and EB are members of the team “Computational Systems Biology of Cancer,” Equipe labellisée par la Ligue Nationale Contre le Cancer. The study was also supported by funding to WDH from the Rothschild-Yvette Mayent Fund, NIH (CA92276, GM58015, CA154920). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Synthetic interactions between two mutations in different genes were first revealed in Drosophila by Dobzhansky in the 1940s [1]. Synthetic lethality (SL) describes that two viable single gene mutations lead to lethality (synthetic-lethal) or severely impair growth (synthetic-sick) when combined as a double mutant. This concept was implemented as a powerful research tool for molecular pathway analysis in yeast [2]–[5]. Further refinement introduced more quantitative measures of genetic epistasis [6] and lethality induced by gene overexpression in a mutant background (synthetic dosage-lethality [7]). A genetic interaction is negative or aggravating, when the combined effect of two gene defects is more severe than it is expected from a simple multiplicative model. In a positive or alleviating interaction the effect is less severe than expected. These approaches and measures are now increasingly used in mammalian cells exploiting RNA-mediated gene knockdown technologies [8], [9].

Following a proposal by Hartwell and colleagues [10], SL has been utilized as a therapeutic approach in cancer treatment employing a combination of genetic ablation (loss of tumor suppressor) and chemical inhibition. The first paradigm was set in BRCA1/2-deficient tumor cells, which are synthetic-lethal with inhibition of PolyADP-Ribose Polymerase (PARP) [11]–[13]. Small molecule PARP inhibitors are currently being evaluated in clinical trials in BRCA1- and BRCA2-deficient cancers (e.g. [14]).

The canonical interpretation of SL stipulates two mutually compensatory, parallel pathways capable of performing the same essential function [2]–[4]. Thus, disrupting a single pathway is viable, while disrupting both pathways is lethal. This concept of between-pathway synthetic lethality (bpSL) (Figure 1A) led to the creation of computational approaches aiming at reconstructing interaction networks from pair-wise gene deletions or siRNA-induced gene knock-down screens in yeast and mammals [15]–[17]. However, recent genome-wide genetic interaction data revealed multiple negative interactions between mutations affecting the same molecular pathway or complex [3], [5],[16],[18]–[21]. For example, it was estimated that ∼9% [19] and in another study 14% [17] of all negative genetic interaction clusters belong to the same biological pathway. Several mechanistic models were suggested to explain within-pathway SL (wpSL) [6], [16], [20], [22]. The deletion of a gene might lead only to a partial degradation of an essential pathway which might be tolerable, whereas the double mutation leads to complete pathway degradation and lethality (Figure 1B1). This is especially relevant for the interpretation of siRNA-based screens where the efficiency of a particular gene knock down is uncertain. A second possible mechanism suggests that steps in an essential pathway are internally redundant (Figure 1B2). Lastly, two mutations may cumulatively degrade an essential protein complex, whereas they are individually viable (Figure 1B3). This mechanism is consistent with the observations that molecular complexes are frequently characterized by the dominance of negative over positive genetic interactions between their components [18]. wpSL interactions between defects in components of a single protein complex are highly enriched for complexes with an essential component [16], [22]. It was estimated that the contribution from within-complex interactions to the total number of within-pathway negative interactions does not exceed 7% [19]. Common to these mechanistic explanations of wpSL is that they involve either an essential pathway or an essential protein complex.

thumbnail
Download:
Figure 1. Schematic representation of mechanisms of (A) between-pathways Synthetic Lethality, (B) within-pathway Synthetic Lethality, and (C) within-reversible-pathway Synthetic Lethality.

S: substrate, I: intermediate, P: product, C1–4: components of a protein complex. Red crosses define single mutations. For details see text.

https://doi.org/10.1371/journal.pcbi.1003016.g001

Here, we highlight a novel scenario of wpSL involving two components of a non-essential pathway. The view of molecular pathways as unidirectional, linear reaction cascades is too simplistic. Pathway steps can be reversible which leads to forward and backward propagation of molecular events along the pathway increasing robustness and fidelity of the overall process [23]–[28]. Affecting both forward and reverse steps of the pathway by abrogating the corresponding enzymes creates scenarios in which the pathway flow can be trapped in an intermediate state that may be toxic to the cell or deprive the cell of a limiting resource (Figure 1C). This can create a genetic scenario we define as within-reversible-pathway synthetic lethality (wrpSL), which is the subject of this study.

Here, we study bpSL and wrpSL scenarios using mathematical modeling to better understand the system properties of these genetic relationships. We present a simplified model of the pathway applicable for its formal analytical study and performed in silico simulations for bpSL and wrpSL as well as synthetic dosage effects. Our main experimentally confirmed examples of wrpSL are in the homologous recombination DNA repair pathway. Homologous recombination (HR) is an important mechanism to maintain genome integrity [29] (Section S1 and Figure S1 in for more discussion). Analysis of yeast gene interaction and pathway data suggests broad applicability of this novel concept.

Results

Homologous recombination and other DNA repair pathways are a hotspot for negative genetic interactions

In order to assess the importance of within-pathway negative interactions we ranked all pathways from the KEGG database [30] according to their normalized proportion of negative interactions [5] within each KEGG pathway (Figure 2). This analysis confirms the previous conclusion [19] that only a minority of within-pathway negative interactions can be explained by negative interactions within a complex (Figure 1B3). In our analysis only 12% of all negative interactions were of this type (compared to 7% in [19]). Interestingly, HR ranks at the top with 27 within-pathway negative interactions between 20 KEGG pathway components (Figure 2), of which only a single one affects components of the same protein complex.

thumbnail
Download:
Figure 2. Normalized proportion of negative genetic interactions within one KEGG pathway.

The proportion is related to the number of all possible pairwise interactions S(S−1)/2, where S is the size of the KEGG pathway. Negative interactions between components of a molecular complex were excluded (compare the fourth and the fifth columns). We counted the normalized number of negative genetic interactions within one pathway using recent data on the genome-wide screening of genetic interactions in yeast from [5], using the most stringent filter on the epistasis measure ε. Definitions of yeast signaling and metabolic pathways were taken from KEGG database [30]. Only KEGG pathways with a normalized proportion of ≥1% are shown. DNA repair KEGG pathways are highlighted in bold.

https://doi.org/10.1371/journal.pcbi.1003016.g002

Reversibility of the pathway steps in homologous recombination pathway

Recent studies show that individual reaction steps in HR are reversible [23]–[25] (Figure S1 in Text S1). The Rad51-ssDNA filament is a key intermediate in HR, as it performs the signature reactions of homology search and DNA strand invasion. The formation of this filament is catalyzed by specific co-factors (see Section S1 and Figure S1 in Text S1). The helicase Srs2 specifically targets the Rad51-ssDNA filament for disruption to reverse filament formation [31], [32], [33]. The reversibility of the Rad51-ssDNA filament sets a new paradigm and draws attention to additional reversible steps and their mechanisms in HR, other DNA repair processes, and unrelated molecular pathways.

Simplest abstract toy model of DNA repair pathway with reversible steps and a toxic intermediate

We derived the simplest linear mathematical model of a main DNA repair pathway with reversible steps and a toxic intermediate, and a compensatory pathway, which can recapitulate bpSL and wrpSL (Figure 3). Each state transition is catalyzed by an abstract enzyme, which may correspond to several biological entities (compare Figure S1 in Text S1 with Figure 3). In wrpSL trapping of toxic intermediate I is caused by defects in the first backward reaction (I→S, R1, k−1) and the second forward reaction (I→P, F2, k2). The reversibility of the second step (reaction P→I, R1) is not essential for wrpSL to occur, but might be important for quantitative pathway characteristics. Introduction of a final irreversible step (Figure S2 in Text S1) would result in a kinetic proofreading mechanism [34] (see Figures S2, S3 in Text S1 and discussion there). Such a mechanism increases the robustness of DNA repair, as it avoids a futile P↔I cycle. However, in this simplest model we eliminated the final irreversible step to allow us analyzing the most essential features of wrpSL (Section S2 in Text S1). Figure 4 explores conditions for various cellular fates (Normal Robust, NR: no single knockout leads to lethality (Figure S4 in Text S1); Normal Fragile, NF: single knockout can lead to lethality (Figure S4 in Text S1); Compensated, C: repair is performed by compensatory pathway; death due to DNA Damage, DD: steady state probability of DNA damage >50%; and Death due to Toxic intermediates, DT: steady state probability of toxic intermediate >50%). Figure 5 visualizes parametric conditions (see Section S2B in Text S1 for discussion).

thumbnail
Download:
Figure 3. Abstract representation of the recombinational DNA repair pathway.

The scheme represents a simplified version of the pathway depicted in Figure S1 in Text S1 using the same abbreviations. Dynamic states S, I, P represent DNA damage substrate (S), toxic intermediate (I), and the product of repair (P). F1 (e.g. Rad51), F2 (e.g. Rad54), R1 (e.g. Srs2), R2 (e.g. Mph1) are enzymes in the main pathway, and F3-EC represent enzymes in the compensatory pathway. ks signify the kinetic rates of the model state transition steps. The two types of synthetic lethality (SL) are indicated: Red Xs – classic scenario of between-pathway SL (e.g. BRCA1 or BRCA2 mutant and PARP inhibition); Green Xs – within-reversible-pathway SL between two mutations in genes acting in a single non-essential pathway leading to the accumulation of a toxic intermediate (e.g. srs2 rad54 double mutant). For more discussion see also Section S1 in Text S1.

https://doi.org/10.1371/journal.pcbi.1003016.g003

thumbnail
Download:
Figure 4. Pathway steady states for various combinations of parameters in the mathematical model of DNA repair with reversible steps and a toxic intermediate.

The classification of pathway steady states depends on the values of three control parameters, representing the ratios of some kinetic rates of the model (for more details see Section S2A in Text S1): The diagrams visualize the values of the control parameters where thickness of the solid arrows shows the relative value of the corresponding kinetic rate and the dashed arrows represent kinetic rates which values are irrelevant for a given scenario. Color coding shows partial orderings of the parameters important for a given pathway state (thickness of edges of the same color should be compared but not between color.

https://doi.org/10.1371/journal.pcbi.1003016.g004

thumbnail
Download:
Figure 5. Visualization of the steady states for the toy model in dependence of various parameters.

The phase diagrams (r2×r1 plane) show the qualitative behavior for small, intermediate and large values of the model parameters. Red color corresponds to the state for which Ps>0.5 (DNA is repaired with probability >50%). If Ps<0.5 then the color is chosen as green if the probability of trapping in the intermediate state (I) is bigger than the propability of the initial unrepaired state (S), and as blue in the opposite case. The case k1 = 0 (F1↓) is represented separately on the right. In this case, another parameter is used instead of r1 for the phase plane. The case k−1 = 0 (R1↓) is treated separately on the top and only the r2 value is varied. The 14 model simulations listed in Figure 6 are shown by the circled numbers in the position of the chosen parameters.

https://doi.org/10.1371/journal.pcbi.1003016.g005

Examples of numerical simulations of the toy model

Using analytical study and numerical simulations with some characteristic choices of kinetic rate values, we explored the dynamical behaviors of the simplest model (see Figure 6). Here, we discuss the qualitative results and interpretations, while the more formal derivation of these statements is found in Section S3 in Text S1. To illustrate the static and dynamic properties of the toy model, we selected two typical positions (Figure 5 #1, #2) corresponding to NR and NF pathway states, respectively. From these “normal” conditions we simulated a number of single and double mutant/overexpression conditions as shown in Figure 5 (see also Figure S4 in Text S1).

thumbnail
Download:
Figure 6. Modeling possible scenarios of single and double synthetic lethal mutations.

Pathway steady states depicted as Normal states (N) as Normal Robust (NR) and Normal Fragile (NF); Mutant states (M) as compensated state dependent on compensatory pathway EC(C); death from DNA damage (DD) and death from toxicity (DT). Dynamic plots show prediction of model for evolution of substrate (S), intermediate (I) and product (P) amounts over time corresponding to the choice of kinetic parameters shown on the Model diagram. (X-axis)-time, (Y-axis)-substances level. F1, F2, R1, R2, and EC refer to the enzymes catalyzing the two forward and two backward reactions as well as the compensatory pathway, respectively (see Figure 3). (↓)-complete knock-down or mutational loss of function; (↑)-over-expression. A deletion mutant was simulated by setting the corresponding kinetic rate constant to zero. An overexpression mutant was simulated by setting the corresponding kinetic rate constant sufficiently high to have a qualitative effect onto the steady state or the dynamics in the simulations. Some double mutants are not shown due to their triviality (such as F1↓F2↓) or difficulties with interpretation (such as F1↓R1↓). For more modeling scenarios see Figure S5 in Text S1.

https://doi.org/10.1371/journal.pcbi.1003016.g006

Normal pathway vs. irreversible pathway: reversibility does not affect efficiency, but is essential for pathway robustness.

We assume that in the NR and NF states under normal conditions forward steps of the pathways are faster than the backward steps. In this case the presence of reversible steps does not significantly affect the numerical solutions of the model equations, i.e. the efficiency of DNA repair (Figure 5 cases 1–3; Figure 6, rows 1–3; Figure S4 in Text S1). However, the presence of reversible steps is important from the robustness point of view (Figure 5 cases 3, 7; Figure 6, rows 3, 7; Figure S4 in Text S1). The gain in robustness through backward reactions is also illustrated by the fact that the difference between the NR and NF states is entirely defined by the ratio of the two backward reactions (rb). The model recapitulates the fact that single enzyme deficiencies are not capable to completely block the overall DNA repair network. The F2 knock-out in NF is the only exception. When the first backward reaction is relatively slow , the compensatory pathway cannot rescue the lethal phenotype (Figure 5 case 8; Figure 7, row 8; Figure S4 in Text S1). This underlines the importance of efficient backward reactions.

thumbnail
Download:
Figure 7. Modeling possible scenarios of single and double synthetic lethal mutations (continued).

Pathway steady states depicted as Normal states (N) as Normal Robust (NR) and Normal Fragile (NF); Mutant states (M) as compensated state dependent on compensatory pathway EC(C); death from DNA damage (DD) and death from toxicity (DT). Dynamic plots show prediction of model for evolution of substrate (S), intermediate (I) and product (P) amounts over time corresponding to the choice of kinetic parameters shown on the Model diagram. (X-axis)-time, (Y-axis)-substances level. F1, F2, R1, R2, and EC refer to the enzymes catalyzing the two forward and two backward reactions as well as the compensatory pathway, respectively (see Figure 3). (↓)-complete knock-down or mutational loss of function; (↑)-over-expression. A deletion mutant was simulated by setting the corresponding kinetic rate constant to zero. An overexpression mutant was simulated by setting the corresponding kinetic rate constant sufficiently high to have a qualitative effect onto the steady state or the dynamics in the simulations. Some double mutants are not shown due to their triviality (such as F1↓F2↓) or difficulties with interpretation (such as F1↓R1↓). For more modeling scenarios see Figure S5 in Text S1.

https://doi.org/10.1371/journal.pcbi.1003016.g007

Between-pathway Synthetic Lethality (bpSL): Disruption of main and compensatory pathways.

The model recapitulates canonical bpSL by putting two kinetic rates to zero simultaneously: k1 = k3 = 0 which causes cell death from accumulation of DNA damage (DD) (Figure 5 case 4; Figure 6, row 4; Figure S4 in Text S1). Alternatively, in the case of k2 = k3 = 0 the pathway state will be DT, resulting in cell death from accumulation of toxic intermediates.

Within-reversible-pathway Synthetic Lethality (wrpSL): Trapping a toxic intermediate.

The wrpSL case is modeled by putting . As it can be seen (Figure 5 case 5; Figure 6, row 5; Figure S4 in Text S1), this creates a kinetic trap in the toxic intermediate state I (DT). The model recapitulates the situation observed in the srs2 rad54 double mutant [23], [35], [36]. Modeling predicts that mammalian cells deficient in RAD54 or its close paralog RAD54B should also be sensitive to inhibition of enzymes, such as FBH1, FANCJ, or RECQ5 [37]–[39] (Figure 5 cases 5, 8; Figure 6, row 5; Figure 7, row 8; Figure S4 in Text S1), that are postulated or shown to dissociate the RAD51-ssDNA filament (Figures S1, S2 in Text S1).

Aspects of reversibility.

Blocking the pathway upstream by eliminating the F1 reaction in the normal pathway states (NR and NF) does not lead to cell death. However, this affects the kinetics of the process, because repair is carried out by the compensatory pathway (state C), which is slower than the principal one (Figure 5 case 6; Figure 6, row 6; Figure S4 in Text S1). Therefore, the time of repair increases and the cell fate depends on the efficiency of the compensatory pathway, which might be cell type-specific. The situation may result in cell death if the kinetics of the process are too slow or the efficiency of DNA damage repair by the compensatory pathway is below an essential threshold. This scenario may represent cancer cells with certain DNA repair defects, which can be compensated by alternative pathways. In such cases, therapeutic targeting of the compensatory pathways may result in fast accumulation of DNA damage, selectively killing the cancer cells, representing the classic bpSL concept.

When the pathway is blocked downstream at F2 (Figure 5 case 7; Figure 6, row 7; Figure S4 in Text S1) the time of repair significantly increases, as the pathway is only partially compensated. This scenario can be rescued by reducing k1 or increasing k−1 (state C). The scenario recapitulates RAD54-deficient cells, where sub-lethal levels of toxic intermediates and unrepaired DNA increase genetic instability [36], [40]. This last scenario of a downstream block (F2 knockout: Figure 5 case 7; Figure 6, row 7; Figure S4 in Text S1) can be rescued by overexpression of R1 (Figure 5 case 9; Figure 7, row 9; Figure S4 in Text S1). The modeling predicts that this will accelerate the backward reactions, favor the compensatory route and increase repair of damaged DNA by the compensatory pathway (state C). The modeling also highlights the benefit of reversibility. In Figure 7 (row 8), a decrease in reversibility (10-fold reduction in R1) leads to significantly less repaired product in mutants that block the downstream step (F2). Complete inhibition of the R1 step leads to wrpSL in cells lacking F2, such as RAD54B-deficient cells (Figure 5 case 8; Figure 7, row 8; Figure S4 in Text S1, see above). The modeling suggests that cells deficient in the downstream forward reaction (F2, e.g. RAD54B−/−) are potentially under selective pressure to strengthen the upstream backward reaction (R1, e.g. overexpression of FBH1, FANCJ, or RECQ5).

Activating mutations or overexpression can create within-reversible-pathway-synthetic dosage lethality.

The computational model allows making predictions in less intuitive cases of model perturbations. The lethal phenotype by R2 overexpression (Figure 5 case 10; Figure 7, row 10; Figure S4 in Text S1) can be partially rescued by eliminating F1 (Figure 5 case 12; Figure 7, row 12; Figure S4 in Text S1), as repair is carried out by the compensatory pathway (state C). The modeling illustrates how a defect in an early HR gene such as RAD51, BRCA2, RAD52, or the RAD51 paralogs would provide a selective advantage in cells overexpressing enzymes of the R2 backward reaction (yeast Mph1, human RTEL1 or FANCM).

F1 overexpression amplifies the effect of an F2 mutation, creating synthetic-lethal dosage and leading to wrpSL due to the accumulation of toxic intermediates (state DT) (Figure 5 case 13; Figure 7, row 13; Figure S4 in Text S1). This suggests that inhibition of F2 (for example RAD54/RAD54B) would be particularly effective to sensitize cancer cells overexpressing F1 components (e.g. RAD51 paralogs, RAD51, RAD52, BRAC2 [41]–[45]), a non-intuitive insight derived from the mathematical modeling.

Blocking the compensatory pathway and overexpressing R2 also creates synthetic lethal dosage interaction, causing accumulation of toxic intermediates due to the increase in backward reactions kinetics (state DT) (Figure 5 case 14; Figure 7, row 14; Figure S4 in Text S1). However, if the final irreversible step of DNA repair will be taken into consideration (P→I2a/b→EP, Figure 3), slow conversion of all favorable DNA configurations from the futile I↔P cycle into the completely repaired DNA state will rescue cells (similar to scenario depicted in Figure 7, row 10). This example is interesting, because it demonstrates that interference with compensatory pathways by drugs might have a therapeutic importance if the backward steps of the main pathway are considered. Finally, the computational model predicts that a similar trap in the intermediate state can be created by eliminating R1 and overexpressing R2 (Figure 5 case 3; Figure S4 in Text S1; Figure S5, row 3 in Text S1).

Discussion

Synthetic lethality/sickness and synthetic dosage lethality are important genetic tools to assign individual gene functions into molecular pathways [2]–[4], [7]–[9]. The canonical interpretation for two mutants found to be synthetically lethal or sick stipulates that the encoded gene products function in different parallel pathways that can mutually compensate (bpSL) [2]–[4], [7]–[9], [15]–[17]. However, computational analysis of genetic interaction data combined with protein interaction data revealed multiple negative interactions between mutations affecting functions in the same molecular pathway or complex (wpSL) [3], [5], [16], [18]–[21]. Several mechanisms of wpSL have been proposed (Figure 1B), and they all involve either essential pathways or essential protein complexes. In extension of this fundamental concept of wpSL, there are several cases of SL between mutants encoding proteins acting in HR, a pathway that is not essential in yeast [35], [46]–[48]. We term this novel genetic interaction within-reversible-pathway Synthetic Lethality (wrpSL; Figure 1C) and provide a novel mechanistic explanation for wpSL, which can create SL within a non-essential pathway or between hypomorphic mutations in an essential pathway that is different from a model invoking sequential pathway degradation by accumulation of partial defects of successive steps (Figure 1B1).

Here, we explore by mathematical modeling the system properties of wrpSL. The modeling must make assumptions about the system properties (state transition rates, relative pathway efficiencies, etc.) and identifies several conditions to be met for wrpSL. 1) Reversibility of pathway steps. In fact, only the first pathway step must be reversible, whereas reversibility of the second pathway steps enables additional genetic scenarios. 2) Possibility of kinetic trapping of an intermediate state of the pathway when both the backward and forward reactions are compromised. The trapping per se can be detrimental due to blockage of cell signaling, sequestering an essential compound, or toxicity. We have assumed lethal toxicity in our model. 3) The possibility of rescue by a parallel compensatory pathway may not be strictly required, but highlights the applicability of this concept to non-essential pathways.

The mathematical model is validated by the experimentally observed recombination-dependent SL of the srs2 rad54 double mutant in budding yeast [35] (Figure 3, Figure S1 in Text S1, Figure 6, row 5). Srs2-defective cells are unable to reverse Rad51-ssDNA filaments. These Rad51-ssDNA filaments represent toxic intermediates that accumulate in the cell due to kinetic trapping and interfere with cell viability. The key functions of the Rad54 protein are to assist in DNA strand invasion and allowing DNA synthesis off the invading 3′-end [36]. Hence, in the srs2 rad54 double mutant Rad51-ssDNA filaments and/or D-loops may accumulate forming a toxic intermediate that leads to cell death (Figure S1 in Text S1; Figure 3 green pair and Figure 6, row 5). This interpretation is supported by the observation that lethality in this double mutant is suppressed by a defect in Rad51-ssDNA filament formation (mutations in RAD51, RAD55, RAD57, or RAD52) [49] (see Figure S1 in Text S1), what has been termed recombination-dependent lethality. Preventing Rad51-ssDNA filament formation allows bypass of recombination by alternative means of DNA repair (for DSBs: Nonhomologous endjoining or single-strand annealing; for gaps: Translesion synthesis or fork regression [23]; see Figure S1 in Text S1). The recombination-dependent lethality of srs2 rad54 is not unique and is also found in additional double mutants in recombinational repair genes including the double mutants mph1 mus81, mph1 mms4, srs2 sgs1 and sgs1 (or top3 or rmi1) and mus81 (or mms4) which likely reflect additional examples of wrpSL possibly involving different toxic intermediates [35], [46]–[53]. As discussed in detail in Section S1 in Text S1, the synthetic lethalities involving sgs1 are more complex, because of the multiple roles of Sgs1-Top3-Rmi1 in HR, and could be caused also by other mechanisms of SL.

Further modeling revealed additional genetic conditions including overexpression of specific pathway enzymes that are predicted to lead to wrpSL (Figure 6). The mathematical modeling also reveals the importance of reversible pathway steps, which are validated by genetic and biochemical experiments in yeast [23]–[25]. First, the existence of reversible pathway steps does not affect normal pathway progression (Figure 6, rows 1–3). Second, reversible pathway steps allow much more efficient and timely use of compensatory pathways (Figure 6, row 6). Third, reversible pathways coupled with compensatory pathways avoid lethality of single mutations (Figure 6, row 7). The existence of reversible intermediates in HR, and possibly other molecular pathways, has been proposed to increase the robustness of the overall DNA repair system [23]–[25], and here we provide quantitative modeling evidence and formal analysis of this assertion.

An important question is how general wrpSL might be or whether it is an idiosyncrasy of the recombinational repair pathway. Even if wrpSL were restricted to HR, this concept provides significant potential application in anti-cancer therapy. However, there is considerable evidence that many molecular pathways include reversible steps catalyzed by different enzymes in the forward and backward directions (see Figure 8). Any of those processes can be theoretically trapped into one of their intermediate states if two regulators of forward and backward steps are inactive. In these cases, the accumulating intermediate might be toxic, block proper signal propagation or prevent resource recycling. Focusing on three examples of reversible protein modifications (phosphorylation by Cdc5/dephosphorylation by Cdc14, sumoylation by Slx5–Slx8/desumoylation by Ulp1, Nup60, ubiquitylation by Rad6–Rad18/deubiquitylation by Bre5, Ubp3 or degradation dependent on Doa1, Rpn6; see Figure S6 in Text S1 for details), we found ample evidence in published genetic interaction data that are consistent with the wrpSL mechanism. These examples have not been fully explored, but are consistent with the wrpSL concept and amenable to test specific predictions.

thumbnail
Download:
Figure 8. Examples of molecular processes with alternative pathways and potential to kinetic trap.

A) Multiple post-translational protein modifications (phosphorylation followed by ubiquitination; acetylation followed by phosphorylation; sumotargeted ubiquitination, etc.). M1 and M2 represent two types of post-translational modifications. MOD1–4 represent enzymes catalyzing the reactions. Kinetic trapping of an intermediate modification can drastically disturb the balance between signaling pathways (e.g. #23 in Figure 2) B) Protein folding control by chaperones. Protein folding is controlled by chaperones (#18 in Figure 2) and may generate partially unfolded proteins as toxic intermediates which are subject to degradation [55]. Regulation of protein folding homeostasis is essential for protein pool control [56]; C) Lack of balance between production and detoxification of ROS leads to significant increases or drop in ROS levels that can be detrimental for cell signaling [57]–[59] (#19 in Figure 2). D) Coordinated sumoylation-desumoylation is important for proper signal propagation [60]–[62]. E) Glycan biosynthesis and protein glycosylation depend on the availability of common carrier dolichol phosphate (P-DOL). Correct recycling of P-DOL is important for sustaining the pool and utilization of this carrier in other glycans biosynthesis pathways. Kinetic trapping can consume the pool of P-DOL and perturb cell signaling [63], [64] (#2, 3 and 26 in Figure 2) F) Protein ubiquitylation is not only a tagging signal for degradation, but also involved in signaling. The correct tuning between two functions of ubiquitylation depends on the type and the length of ubiquitin (UB) polymer transferred to the protein and at the ubiquitylation site [65]. Monoubiquitylated proteins participate in signal transduction [66], [67], whereas K48-polyubiquitylated proteins are redirected to the proteasome for proteolysis. Thus, the balance between protein homeostasis and ubiquitin-dependent signaling is essential [68].

https://doi.org/10.1371/journal.pcbi.1003016.g008

In summary, genetic and biochemical data strongly support the significance of the wrpSL mechanism in HR, and existing data are consistent with the notion that wrpSL could be a general, widely applicable type of genetic interaction. This may refine our understanding of relationships between gene products and will help to improve pathway reconstruction. In particular, our mathematical modeling provides a conceptual framework for guiding systematic exploitation of mutations and changes in the expression profiles of HR genes and potentially genes of other pathways to induce SL.

Materials and Methods

Model formalism and model solution

The simplest mathematical model of Figure 3 was converted into a set of linear ordinary differential equations using the standard chemical kinetics formalism. The steady state model properties were analyzed analytically and exemplified with numerical simulations. Classification of the pathway states according to the extreme (large or small) values of the control parameters and the corresponding asymptotic solutions follow the methodology of the asymptotology of reaction networks [54].

Numerical simulations

All numerical simulations were performed using SBTOOLBOX package for Matlab (Section S4 in Text S1).

Supporting Information

Text S1.

Text S1 contains Figure S1–S6, Section S1–S4, and Table S1.

https://doi.org/10.1371/journal.pcbi.1003016.s001

(PDF)

Acknowledgments

We thank J. Hall, M.-H. Stern, T. Popova, L. Calzone, R. Janke, S. Ceballos, C. Fasching, S. Kowalczykowski, M.Lichten and M. Savageau for valuable comments on the manuscript. WDH thanks the Institut Curie and especially Dr. Alain Nicolas for the hospitality during his sabbatical stay.

Author Contributions

Conceived and designed the experiments: AZ IK WDH. Performed the experiments: AZ. Analyzed the data: AZ IK EB WDH. Wrote the paper: AZ IK EB WDH.

References

  1. 1. Dobzhansky T (1946) Genetics of natural populations. XIII. Recombination and variability in populations of Drosophila pseudoobscura. Genetics 31: 269–290.
  2. 2. Collins SR, Miller KM, Maas NL, Roguev A, Fillingham J, et al. (2007) Functional dissection of protein complexes involved in yeast chromosome biology using a genetic interaction map. Nature 446: 806–810.
  3. 3. Bandyopadhyay S, Mehta M, Kuo D, Sung MK, Chuang R, et al. (2010) Rewiring of Genetic Networks in Response to DNA Damage. Science 330: 1385–1389.
  4. 4. Tong AHY, Evangelista M, Parsons AB, Xu H, Bader GD, et al. (2001) Systematic genetic analysis with ordered arrays of yeast deletion mutants. Science 294: 2363–2368.
  5. 5. Costanzo M, Baryshnikova A, Bellay J, Kim Y, Spear ED, et al. (2010) The genetic landscape of a cell. Science 327: 425–431.
  6. 6. Boone C, Bussey H, Andrews BJ (2007) Exploring genetic interactions and networks with yeast. Nature Reviews Genetics 8: 437–449.
  7. 7. Kroll ES, Hyland KM, Hieter P, Li JJ (1996) Establishing genetic interactions by a synthetic dosage lethality phenotype. Genetics 143: 95–102.
  8. 8. Silva JM, Marran K, Parker JS, Silva J, Golding M, et al. (2008) Profiling essential genes in human mammary cells by multiplex RNAi screening. Science 319: 617–620.
  9. 9. Iorns E, Lord CJ, Turner N, Ashworth A (2007) Utilizing RNA interference to enhance cancer drug discovery. Nature Reviews Drug Discovery 6: 556–568.
  10. 10. Hartwell LH, Szankasi P, Roberts CJ, Murray AW, Friend SH (1997) Integrating genetic approaches into the discovery of anticancer drugs. Science 278: 1064–1068.
  11. 11. Farmer H, McCabe N, Lord CJ, Tutt ANJ, Johnson DA, et al. (2005) Targeting the DNA repair defect in BRCA mutant cells as a therapeutic strategy. Nature 434: 917–921.
  12. 12. Bryant HE, Schultz N, Thomas HD, Parker KM, Flower D, et al. (2005) Specific killing of BRCA2-deficient tumours with inhibitors of poly(ADP-ribose) polymerase. Nature 434: 913–917.
  13. 13. Helleday T (2011) The underlying mechanism for the PARP and BRCA synthetic lethality: Clearing up the misunderstandings. Molecular Oncology 5: 387–393.
  14. 14. Fong PC, Boss DS, Yap TA, Tutt A, Wu P, et al. (2009) Inhibition of poly(ADP-ribose) polymerase in tumors from BRCA mutation carriers. New England Journal of Medicine 361: 123–134.
  15. 15. Mizuarai S, Kotani H (2010) Synthetic lethal interactions for the development of cancer therapeutics: biological and methodological advancements. Human Genetics 128: 567–575.
  16. 16. Baryshnikova A, Costanzo M, Kim Y, Ding H, Koh J, et al. (2010) Quantitative analysis of fitness and genetic interactions in yeast on a genome scale. Nature Methods 7: 1017–1024.
  17. 17. Kelley R, Ideker T (2005) Systematic interpretation of genetic interactions using protein networks. Nature Biotechnology 23: 561–566.
  18. 18. Michaut M, Baryshnikova A, Costanzo M, Myers CL, Andrews BJ, et al. (2011) Protein complexes are central in the yeast genetic landscape. PLoS Comput Biol 7: e1001092.
  19. 19. Bellay J, Atluri G, Sing TL, Toufighi K, Costanzo M, et al. (2011) Putting genetic interactions in context through a global modular decomposition. Genome Research 21: 1375–1387.
  20. 20. Costanzo M, Baryshnikova A, Myers CL, Andrews B, Boone C (2011) Charting the genetic interaction map of a cell. Current Opinion in Biotechnology 22: 66–74.
  21. 21. Sopko R, Huang D, Preston N, Chua G, Papp B, et al. (2006) Mapping pathways and phenotypes by systematic gene overexpression. Molecular Cell 21: 319–330.
  22. 22. Bandyopadhyay S, Kelley R, Krogan NJ, Ideker T (2008) Functional maps of protein complexes from quantitative genetic interaction data. PLoS Comput Biol 4: e1000065.
  23. 23. Heyer WD, Ehmsen KT, Liu J (2010) Regulation of homologous recombination in eukaryotes. Annual Reviews Genetics 44: 113–139.
  24. 24. Kanaar R, Wyman C, Rothstein R (2008) Quality control of DNA break metabolism: in the ‘end’, it's a good thing. EMBO Journal 27: 581–588.
  25. 25. Symington LS, Heyer WD (2006) Some disassembly required: role of DNA translocases in the disruption of recombination intermediates and dead-end complexes. Genes and Development 20: 2479–2486.
  26. 26. Alon U (2003) Biological Networks: The tinkerer as an engineer. Science 301: 1866–1867.
  27. 27. Savageau M (1971) Parameter sensitivity as a criterion for evaluating and comparing the performance of biochemical systems. Nature 229: 542–544.
  28. 28. Kitano H (2004) Biological robustness. Nature Reviews Genetics 5: 826–837.
  29. 29. Li X, Heyer WD (2008) Homologous recombination in DNA repair and DNA damage tolerance. Cell Research 18: 99–113.
  30. 30. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M (2012) KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Research 40: D109–114.
  31. 31. Veaute X, Jeusset J, Soustelle C, Kowalczykowski SC, Le Cam E, et al. (2003) The Srs2 helicase prevents recombination by disrupting Rad51 nucleoprotein filaments. Nature 423: 309–312.
  32. 32. Krejci L, Van Komen S, Li Y, Villemain J, Reddy MS, et al. (2003) DNA helicase Srs2 disrupts the Rad51 presynaptic filament. Nature 423: 305–309.
  33. 33. Liu J, Renault L, Veaute X, Fabre F, Stahlberg H, et al. (2011) Rad51 paralogs Rad55-Rad57 balance the anti-recombinase function of Srs2 in Rad51 pre-synaptic filament formation. Nature 479: 245–248.
  34. 34. Hopfield JJ (1974) Kintic proofreading: a new mechanism for reducing errors in biosynthetic processes requiting high specifcity. Proceedings National Academy of Sciences USA 71: 4135–4139.
  35. 35. Palladino F, Klein HL (1992) Analysis of mitotic and meiotic defects in Saccharomyces cerevisiae SRS2 DNA helicase mutants. Genetics 132: 23–37.
  36. 36. Heyer WD, Li X, Rolfsmeier M, Zhang XP (2006) Rad54: the Swiss Army knife of homologous recombination? Nucleic Acids Research 34: 4115–4125.
  37. 37. Fugger K, Mistrik M, Danielsen JR, Dinant C, Falck J, et al. (2009) Human Fbh1 helicase contributes to genome maintenance via pro- and anti-recombinase activities. Journal of Cell Biology 186: 655–663.
  38. 38. Sommers JA, Rawtani N, Gupta R, Bugreev DV, Mazin AV, et al. (2009) FANCJ uses its motor ATPase to destabilize protein-DNA complexes, unwind triplexes, and inhibit RAD51 strand exchange. Journal of Biological Chemistry 284: 7502–7514.
  39. 39. Schwendener S, Raynard S, Paliwal S, Cheng A, Kanagaraj R, et al. (2010) Physical Interaction of RECQ5 Helicase with RAD51 Facilitates Its Anti-recombinase Activity. Journal of Biological Chemistry 285: 15739–15745.
  40. 40. McManus KJ, Barrett IJ, Nouhi Y, Hieter P (2009) Specific synthetic lethal killing of RAD54B-deficient human colorectal cancer cells by FEN1 silencing. Proceedings National Academy of Sciences USA 106: 3276–3281.
  41. 41. Klein HL (2008) The consequences of Rad51 overexpression for normal and tumor cells. DNA Repair 7: 686–693.
  42. 42. Feng ZH, Scott SP, Bussen W, Sharma GG, Guo G, et al. (2010) Rad52 inactivation is synthetically lethal with BRCA2 deficiency. Proceedings National Academy of Sciences USA 108: 686–691.
  43. 43. Slupianek A, Schmutte C, Tombline G, Nieborowska-Skorska M, Hoser G, et al. (2001) BCR/ABL regulates mammalian RecA homologs, resulting in drug resistance. Molecular Cell 8: 795–806.
  44. 44. Maacke H, Opitz S, Jost K, Hamdorf W, Henning W, et al. (2000) Over-expression of wild-type Rad51 correlates with histological grading of invasive ductal breast cancer. International Journal of Cancer 88: 907–913.
  45. 45. Christodoulopoulos G, Malapetsa A, Schipper H, Golub E, Radding C, et al. (1999) Chlorambucil induction of HsRad51 in B-cell chronic lymphocytic leukemia. Clinical Cancer Research 5: 2178–2184.
  46. 46. Fabre F, Chan A, Heyer WD, Gangloff S (2002) Alternate pathways involving Sgs1/Top3, Mus81/Mms4, and Srs2 prevent formation of toxic recombination intermediates from single-stranded gaps created by DNA replication. Proceedings National Academy of Sciences USA 99: 16887–16892.
  47. 47. Gangloff S, Soustelle C, Fabre F (2000) Homologous recombination is responsible for cell death in the absence of the Sgs1 and Srs2 helicases. Nature Genetics 25: 192–194.
  48. 48. Maftahi M, Hope JC, Delgado-Cruzata L, Han CS, Freyer GA (2002) The severe slow growth of Delta srs2 Delta rqh1 in Schizosaccharomyces pombe is suppressed by loss of recombination and checkpoint genes. Nucleic Acids Research 30: 4781–4792.
  49. 49. Schild D (1995) Suppression of a new allele of the yeast RAD52 gene by overexpression of RAD51, mutations in srs2 and ccr4, or mating-type heterozygosity. Genetics 140: 115–127.
  50. 50. Trowbridge K, McKim K, Brill SJ, Sekelsky J (2007) Synthetic lethality of Drosophila in the absence of the MUS81 endonuclease and the DmBlm helicase is associated with elevated apoptosis. Genetics 176: 1993–2001.
  51. 51. De Muyt A, Jessop L, Kolar E, Sourirajan A, Chen JH, et al. (2012) BLM Helicase Ortholog Sgs1 Is a Central Regulator of Meiotic Recombination Intermediate Metabolism. Molecular Cell 46: 43–53.
  52. 52. Zakharyevich K, Tang SM, Ma YM, Hunter N (2012) Delineation of Joint Molecule Resolution Pathways in Meiosis Identifies a Crossover-Specific Resolvase. Cell 149: 334–347.
  53. 53. St Onge RP, Mani R, Oh J, Proctor M, Fung E, et al. (2007) Systematic pathway analysis using high-resolution fitness profiling of combinatorial gene deletions. Nature Genetics 39: 199–206.
  54. 54. Gorban A, Radulescu O, Zinovyev A (2010) Asymptotology of chemical reactions networks.. Chemical Engineering Sciences 65: 2310–2324.
  55. 55. Hartl FU, Bracher A, Hayer-Hartl M (2011) Molecular chaperones in protein folding and proteostasis. Nature 475: 324–332.
  56. 56. Voisine C, Pedersen JS, Morimoto RI (2010) Chaperone networks: tipping the balance in protein folding diseases. Neurobiological Diseases 40: 12–20.
  57. 57. Bartz RR, Piantadosi CA (2010) Clinical review: oxygen as a signaling molecule. Critical Care 14: 234.
  58. 58. Hamanaka RB, Chandel NS (2012) Mitochondrial reactive oxygen species regulate cellular signaling and dictate biological outcomes. Trends in Biochemical Sciences 35: 505–513.
  59. 59. Karihtala P, Soini Y (2007) Reactive oxygen species and antioxidant mechanisms in human tissues and their relation to malignancies. APMIS 115: 81–103.
  60. 60. Gareau JR, Lima CD (2010) The SUMO pathway: emerging mechanisms that shape specificity, conjugation and recognition. Nature Reviews Molecular Cell Biology 11: 861–871.
  61. 61. Geiss-Friedlander R, Melchior F (2007) Concepts in sumoylation: a decade on. Nature Reviews Molecular Cell Biology 8: 947–956.
  62. 62. Wilkinson KA, Henley JM (2010) Mechanisms, regulation and consequences of protein SUMOylation. Biochemical Journal 428: 133–145.
  63. 63. Moremen KW, Tiemeyer M, Nairn AV (2012) Vertebrate protein glycosylation: diversity, synthesis and function. Nature Reviews Molecular Cell Biology 13: 448–462.
  64. 64. Schwarz F, Aebi M (2011) Mechanisms and principles of N-linked protein glycosylation. Current Opinion Structural Biology 21: 576–582.
  65. 65. Kulathu Y, Komander D (2012) Atypical ubiquitylation - the unexplored world of polyubiquitin beyond Lys48 and Lys63 linkages. Nature Reviews Molecular Cell Biology 13: 508–523.
  66. 66. Sigismund S, Polo S, Di Fiore PP (2004) Signaling through monoubiquitination. Curr Top Microbiology Immunology 286: 149–185.
  67. 67. Vertegaal AC (2011) Uncovering ubiquitin and ubiquitin-like signaling networks. Chemical Reviews 111: 7923–7940.
  68. 68. Chen X, Zhang B, Fischer JA (2002) A specific protein substrate for a deubiquitinating enzyme: Liquid facets is the substrate of Fat facets. Genes and Development 16: 289–294.