Coherent feedforward transcriptional regulatory motifs enhance drug resistance

Drug tolerance and persistence in bacteria, fungi and cancer cells: Role of non-genetic heterogeneity

A common feature of bacterial, fungal and cancer cell populations upon treatment is the presence of tolerant and persistent cells able to survive, and sometimes grow, even in the presence of usually inhibitory or lethal drug concentrations, driven by non-genetic differences among individual cells in a population. Here we review and compare data obtained on drug survival in bacteria, fungi and cancer cells to unravel common characteristics and cellular pathways, and to point their singularities. This comparative work also allows to cross-fertilize ideas across fields. We particularly focus on the role of gene expression variability in the emergence of cell-cell non-genetic heterogeneity because it represents a possible common basic molecular process at the origin of most persistence phenomena and could be monitored and tuned to help improve therapeutic interventions.

Dynamical model of antibiotic responses linking expression of resistance genes to metabolism explains emergence of heterogeneity during drug exposures

Antibiotic responses in bacteria are highly dynamic and heterogeneous, with sudden exposure of bacterial colonies to high drug doses resulting in the coexistence of recovered and arrested cells. The dynamics of the response is determined by regulatory circuits controlling the expression of resistance genes, which are in turn modulated by the drug's action on cell growth and metabolism. Despite advances in understanding gene regulation at the molecular level, we still lack a framework to describe how feedback mechanisms resulting from the interdependence between expression of resistance and cell metabolism can amplify naturally occurring noise and create heterogeneity at the population level. To understand how this interplay affects cell survival upon exposure, we constructed a mathematical model of the dynamics of antibiotic responses that links metabolism and regulation of gene expression, based on the tetracycline resistancetetoperon inE. coli. We use this model to interpret measurements of growth and expression of resistance in microfluidic experiments, both in single cells and in biofilms. We also implemented a stochastic model of the drug response, to show that exposure to high drug levels results in large variations of recovery times and heterogeneity at the population level. We show that stochasticity is important to determine how nutrient quality affects cell survival during exposure to high drug concentrations. A quantitative description of how microbes respond to antibiotics in dynamical environments is crucial to understand population-level behaviors such as biofilms and pathogenesis.

Evolutionary innovation through transcription factor rewiring in microbes is shaped by levels of transcription factor activity, expression, and existing connectivity

The survival of a population during environmental shifts depends on whether the rate of phenotypic adaptation keeps up with the rate of changing conditions. A common way to achieve this is via change to gene regulatory network (GRN) connections-known as rewiring-that facilitate novel interactions and innovation of transcription factors. To understand the success of rapidly adapting organisms, we therefore need to determine the rules that create and constrain opportunities for GRN rewiring. Here, using an experimental microbial model system with the soil bacterium Pseudomonas fluorescens, we reveal a hierarchy among transcription factors that are rewired to rescue lost function, with alternative rewiring pathways only unmasked after the preferred pathway is eliminated. We identify 3 key properties-high activation, high expression, and preexisting low-level affinity for novel target genes-that facilitate transcription factor innovation. Ease of acquiring these properties is constrained by preexisting GRN architecture, which was overcome in our experimental system by both targeted and global network alterations. This work reveals the key properties that determine transcription factor evolvability, and as such, the evolution of GRNs.

Quantitative systems-based prediction of antimicrobial resistance evolution

Predicting evolution is a fundamental problem in biology with practical implications for treating antimicrobial resistance, which is a complex system-level phenomenon. In this perspective article, we explore the limits of predicting antimicrobial resistance evolution, quantitatively define the predictability and repeatability of microevolutionary processes, and speculate on how these quantities vary across temporal, biological, and complexity scales. The opportunities and challenges for predicting antimicrobial resistance in the context of systems biology are also discussed. Based on recent research, we conclude that the evolution of antimicrobial resistance can be predicted using a systems biology approach integrating quantitative models with multiscale data from microbial evolution experiments.

MaxCal can infer models from coupled stochastic trajectories of gene expression and cell division

Gene expression is inherently noisy due to small numbers of proteins and nucleic acids inside a cell. Likewise, cell division is stochastic, particularly when tracking at the level of a single cell. The two can be coupled when gene expression affects the rate of cell division. Single-cell time-lapse experiments can measure both fluctuations by simultaneously recording protein levels inside a cell and its stochastic division. These information-rich noisy trajectory data sets can be harnessed to learn about the underlying molecular and cellular details that are often not known a priori. A critical question is: How can we infer a model given data where fluctuations at two levels-gene expression and cell division-are intricately convoluted? We show the principle of maximum caliber (MaxCal)-integrated within a Bayesian framework-can be used to infer several cellular and molecular details (division rates, protein production, and degradation rates) from these coupled stochastic trajectories (CSTs). We demonstrate this proof of concept using synthetic data generated from a known model. An additional challenge in data analysis is that trajectories are often not in protein numbers, but in noisy fluorescence that depends on protein number in a probabilistic manner. We again show that MaxCal can infer important molecular and cellular rates even when data are in fluorescence, another example of CST with three confounding factors-gene expression noise, cell division noise, and fluorescence distortion-all coupled. Our approach will provide guidance to build models in synthetic biology experiments as well as general biological systems where examples of CSTs are abundant.

Systematic integration of protein-affecting mutations, gene fusions, and copy number alterations into a comprehensive somatic mutational profile

Somatic mutations occur as random genetic changes in genes through protein-affecting mutations (PAMs), gene fusions, or copy number alterations (CNAs). Mutations of different types can have a similar phenotypic effect (i.e., allelic heterogeneity) and should be integrated into a unified gene mutation profile. We developed OncoMerge to fill this niche of integrating somatic mutations to capture allelic heterogeneity, assign a function to mutations, and overcome known obstacles in cancer genetics. Application of OncoMerge to TCGA Pan-Cancer Atlas increased detection of somatically mutated genes and improved the prediction of the somatic mutation role as either activating or loss of function. Using integrated somatic mutation matrices increased the power to infer gene regulatory networks and uncovered the enrichment of switch-like feedback motifs and delay-inducing feedforward loops. These studies demonstrate that OncoMerge efficiently integrates PAMs, fusions, and CNAs and strengthens downstream analyses linking somatic mutations to cancer phenotypes.

The Extent of Antimicrobial Resistance Due to Efflux Pump Regulation

A key mechanism driving antimicrobial resistance (AMR) stems from the ability of bacteria to up-regulate efflux pumps upon exposure to drugs. The resistance gained by this up-regulation is pliable because of the tight regulation of efflux pump levels. This leads to temporary enhancement in survivability of bacteria due to higher efflux pump levels in the presence of antibiotics, which can be reversed when the cells are no longer exposed to the drug. Knowledge of the extent of resistance thus gained would inform intervention strategies aimed at mitigating AMR. Here, we combine mathematical modeling and experiments to quantify the maximum extent of resistance that efflux pump up-regulation can confer via phenotypic induction in the presence of drugs and genotypic abrogation of regulation. Our model describes the dynamics of drug transport in and out of cells coupled with the associated regulation of efflux pump levels and predicts the increase in the minimum inhibitory concentration (MIC) of drugs due to such regulation. To test the model, we measured the uptake and efflux as well as the MIC of the compound ethidium bromide (EtBr), a substrate of the efflux pump LfrA, in wild-type Mycobacterium smegmatis mc²155, as well as in two laboratory-generated strains. Our model captured the observed EtBr levels and MIC fold-changes quantitatively. Further, the model identified key parameters associated with the resulting resistance, variations in which could underlie the extent to which such resistance arises across different drug-bacteria combinations, potentially offering tunable handles to optimize interventions aimed at minimizing AMR.

Unravelling microbial efflux through mathematical modelling

AbstractMathematical modelling is a useful tool that is increasingly used in the life sciences to understand and predict the behaviour of biological systems. This review looks at how this interdisciplinary approach has advanced our understanding of microbial efflux, the process by which microbes expel harmful substances. The discussion is largely in the context of antimicrobial resistance, but applications in synthetic biology are also touched upon. The goal of this paper is to spark further fruitful collaborations between modellers and experimentalists in the efflux community and beyond.

Non-genetic resistance facilitates survival while hindering the evolution of drug resistance due to intraspecific competition

Rising rates of resistance to antimicrobial drugs threaten the effective treatment of infections across the globe. Drug resistance has been established to emerge from non-genetic mechanisms as well as from genetic mechanisms. However, it is still unclear how non-genetic resistance affects the evolution of genetic drug resistance. We develop deterministic and stochastic population models that incorporate resource competition to quantitatively investigate the transition from non-genetic to genetic resistance during the exposure to static and cidal drugs. We find that non-genetic resistance facilitates the survival of cell populations during drug treatment while hindering the development of genetic resistance due to competition between the non-genetically and genetically resistant subpopulations. Non-genetic resistance in the presence of subpopulation competition increases the fixation times of drug resistance mutations, while increasing the probability of mutation before population extinction during cidal drug treatment. Intense intraspecific competition during drug treatment leads to extinction of susceptible and non-genetically resistant subpopulations. Alternating between drug and no drug conditions results in oscillatory population dynamics, increased resistance mutation fixation timescales, and reduced population survival. These findings advance our fundamental understanding of the evolution of resistance and may guide novel treatment strategies for patients with drug-resistant infections.

Mathematical Modelling Highlights the Potential for Genetic Manipulation as an Adjuvant to Counter Efflux-Mediated MDR in Salmonella

Bacteria have developed resistance to antibiotics by various mechanisms, notable amongst these is the use of permeation barriers and the expulsion of antibiotics via efflux pumps. The resistance-nodulation-division (RND) family of efflux pumps is found in Gram-negative bacteria and a major contributor to multidrug resistance (MDR). In particular, Salmonella encodes five RND efflux pump systems: AcrAB, AcrAD, AcrEF, MdsAB and MdtAB which have different substrate ranges including many antibiotics. We produce a spatial partial differential equation (PDE) model governing the diffusion and efflux of antibiotic in Salmonella, via these RND efflux pumps. Using parameter fitting techniques on experimental data, we are able to establish the behaviour of multiple wild-type and efflux mutant Salmonella strains, which enables us to produce efflux profiles for each individual efflux pump system. By combining the model with a gene regulatory network (GRN) model of efflux regulation, we simulate how the bacteria respond to their environment. Finally, performing a parameter sensitivity analysis, we look into various different targets to inhibit the efflux pumps. The model provides an in silico framework with which to test these potential adjuvants to counter MDR.

UPMaBoSS: A Novel Framework for Dynamic Cell Population Modeling

Mathematical modeling aims at understanding the effects of biological perturbations, suggesting ways to intervene and to reestablish proper cell functioning in diseases such as cancer or in autoimmune disorders. This is a difficult task for obvious reasons: the level of details needed to describe the intra-cellular processes involved, the numerous interactions between cells and cell types, and the complex dynamical properties of such populations where cells die, divide and interact constantly, to cite a few. Another important difficulty comes from the spatial distribution of these cells, their diffusion and motility. All of these aspects cannot be easily resolved in a unique mathematical model or with a unique formalism. To cope with some of these issues, we introduce here a novel framework, UPMaBoSS (for Update Population MaBoSS), dedicated to modeling dynamic populations of interacting cells. We rely on the preexisting tool MaBoSS, which enables probabilistic simulations of cellular networks. A novel software layer is added to account for cell interactions and population dynamics, but without considering the spatial dimension. This modeling approach can be seen as an intermediate step towards more complex spatial descriptions. We illustrate our methodology by means of a case study dealing with TNF-induced cell death. Interestingly, the simulation of cell population dynamics with UPMaBoSS reveals a mechanism of resistance triggered by TNF treatment. Relatively easy to encode, UPMaBoSS simulations require only moderate computational power and execution time. To ease the reproduction of simulations, we provide several Jupyter notebooks that can be accessed within the CoLoMoTo Docker image, which contains all software and models used for this study.

Does transcriptional heterogeneity facilitate the development of genetic drug resistance?

Non-genetic forms of antimicrobial (drug) resistance can result from cell-to-cell variability that is not encoded in the genetic material. Data from recent studies also suggest that non-genetic mechanisms can facilitate the development of genetic drug resistance. We speculate on how the interplay between non-genetic and genetic mechanisms may affect microbial adaptation and evolution during drug treatment. We argue that cellular heterogeneity arising from fluctuations in gene expression, epigenetic modifications, as well as genetic changes contribute to drug resistance at different timescales, and that the interplay between these mechanisms enhance pathogen resistance. Accordingly, developing a better understanding of the role of non-genetic mechanisms in drug resistance and how they interact with genetic mechanisms will enhance our ability to combat antimicrobial resistance. Also see the video abstract here: https://youtu.be/aefGpdh-bgU.

Time dependent asymptotic analysis of the gene regulatory network of the AcrAB-TolC efflux pump system in gram-negative bacteria

Efflux pumps are a mechanism of intrinsic and evolved resistance in bacteria. If an efflux pump can expel an antibiotic so that its concentration within the cell is below a killing threshold the bacteria are resistant to the antibiotic. Efflux pumps may be specific or they may pump various different substances. This is why many efflux pumps confer multi drug resistance (MDR). In particular over expression of the AcrAB−TolC efflux pump system confers MDR in both Salmonella and Escherichia coli. We consider the complex gene regulation network that controls expression of genes central to controlling the efflux associated genes acrAB and acrEF in Salmonella. We present the first mathematical model of this gene regulatory network in the form of a system of ordinary differential equations. Using a time dependent asymptotic analysis, we examine in detail the behaviour of the efflux system on various different timescales. Asymptotic approximations of the steady states provide an analytical comparison of targets for efflux inhibition.

Optimal transition paths of phenotypic switching in a non-Markovian self-regulation gene expression

Gene expression is a complex biochemical process involving multiple reaction steps, creating molecular memory because the probability of waiting time between consecutive reaction steps no longer follows exponential distributions. What effect the molecular memory has on metastable states in gene expression remains not fully understood. Here, we study transition paths of switching between bistable states for a non-Markovian model of gene expression equipped with a self-regulation. Employing the large deviation theory for this model, we analyze the optimal transition paths of switching between bistable states in gene expression, interestingly finding that dynamic behaviors in gene expression along the optimal transition paths significantly depend on the molecular memory. Moreover, we discover that the molecular memory can prolong the time of switching between bistable states in gene expression along the optimal transition paths. Our results imply that the molecular memory may be an unneglectable factor to affect switching between metastable states in gene expression.

Advancing Antimicrobial Resistance Research Through Quantitative Modeling and Synthetic Biology

Antimicrobial resistance (AMR) is an emerging global health crisis that is undermining advances in modern medicine and, if unmitigated, threatens to kill 10 million people per year worldwide by 2050. Research over the last decade has demonstrated that the differences between genetically identical cells in the same environment can lead to drug resistance. Fluctuations in gene expression, modulated by gene regulatory networks, can lead to non-genetic heterogeneity that results in the fractional killing of microbial populations causing drug therapies to fail; this non-genetic drug resistance can enhance the probability of acquiring genetic drug resistance mutations. Mathematical models of gene networks can elucidate general principles underlying drug resistance, predict the evolution of resistance, and guide drug resistance experiments in the laboratory. Cells genetically engineered to carry synthetic gene networks regulating drug resistance genes allow for controlled, quantitative experiments on the role of non-genetic heterogeneity in the development of drug resistance. In this perspective article, we emphasize the contributions that mathematical, computational, and synthetic gene network models play in advancing our understanding of AMR to discover effective therapies against drug-resistant infections.

Model genotype-phenotype mappings and the algorithmic structure of evolution

Cancers are complex dynamic systems that undergo evolution and selection. Personalized medicine approaches in the clinic increasingly rely on predictions of tumour response to one or more therapies; these predictions are complicated by the inevitable evolution of the tumour. Despite enormous amounts of data on the mutational status of cancers and numerous therapies developed in recent decades to target these mutations, many of these treatments fail after a time due to the development of resistance in the tumour. The emergence of these resistant phenotypes is not easily predicted from genomic data, since the relationship between genotypes and phenotypes, termed the genotype-phenotype (GP) mapping, is neither injective nor functional. We present a review of models of this mapping within a generalized evolutionary framework that takes into account the relation between genotype, phenotype, environment and fitness. Different modelling approaches are described and compared, and many evolutionary results are shown to be conserved across studies despite using different underlying model systems. In addition, several areas for future work that remain understudied are identified, including plasticity and bet-hedging. The GP-mapping provides a pathway for understanding the potential routes of evolution taken by cancers, which will be necessary knowledge for improving personalized therapies.

Experimental evolution reveals a general role for the methyltransferase Hmt1 in noise buffering

Cell-to-cell heterogeneity within an isogenic population has been observed in prokaryotic and eukaryotic cells. Such heterogeneity often manifests at the level of individual protein abundance and may have evolutionary benefits, especially for organisms in fluctuating environments. Although general features and the origins of cellular noise have been revealed, details of the molecular pathways underlying noise regulation remain elusive. Here, we used experimental evolution of Saccharomyces cerevisiae to select for mutations that increase reporter protein noise. By combining bulk segregant analysis and CRISPR/Cas9-based reconstitution, we identified the methyltransferase Hmt1 as a general regulator of noise buffering. Hmt1 methylation activity is critical for the evolved phenotype, and we also show that two of the Hmt1 methylation targets can suppress noise. Hmt1 functions as an environmental sensor to adjust noise levels in response to environmental cues. Moreover, Hmt1-mediated noise buffering is conserved in an evolutionarily distant yeast species, suggesting broad significance of noise regulation.

Experimental evolution in yeast reveals that the methyltransferase Hmt1 functions as a mediator connecting environmental stimuli to cellular noise; Hmt1-mediated noise buffering is conserved in an evolutionarily distant yeast.

Mathematical modeling and computational prediction of cancer drug resistance

Diverse forms of resistance to anticancer drugs can lead to the failure of chemotherapy. Drug resistance is one of the most intractable issues for successfully treating cancer in current clinical practice. Effective clinical approaches that could counter drug resistance by restoring the sensitivity of tumors to the targeted agents are urgently needed. As numerous experimental results on resistance mechanisms have been obtained and a mass of high-throughput data has been accumulated, mathematical modeling and computational predictions using systematic and quantitative approaches have become increasingly important, as they can potentially provide deeper insights into resistance mechanisms, generate novel hypotheses or suggest promising treatment strategies for future testing. In this review, we first briefly summarize the current progress of experimentally revealed resistance mechanisms of targeted therapy, including genetic mechanisms, epigenetic mechanisms, posttranslational mechanisms, cellular mechanisms, microenvironmental mechanisms and pharmacokinetic mechanisms. Subsequently, we list several currently available databases and Web-based tools related to drug sensitivity and resistance. Then, we focus primarily on introducing some state-of-the-art computational methods used in drug resistance studies, including mechanism-based mathematical modeling approaches (e.g. molecular dynamics simulation, kinetic model of molecular networks, ordinary differential equation model of cellular dynamics, stochastic model, partial differential equation model, agent-based model, pharmacokinetic–pharmacodynamic model, etc.) and data-driven prediction methods (e.g. omics data-based conventional screening approach for node biomarkers, static network approach for edge biomarkers and module biomarkers, dynamic network approach for dynamic network biomarkers and dynamic module network biomarkers, etc.). Finally, we discuss several further questions and future directions for the use of computational methods for studying drug resistance, including inferring drug-induced signaling networks, multiscale modeling, drug combinations and precision medicine.

Identification of co-evolving temporal networks

BACKGROUND

Biological networks describes the mechanisms which govern cellular functions. Temporal networks show how these networks evolve over time. Studying the temporal progression of network topologies is of utmost importance since it uncovers how a network evolves and how it resists to external stimuli and internal variations. Two temporal networks have co-evolving subnetworks if the evolving topologies of these subnetworks remain similar to each other as the network topology evolves over a period of time. In this paper, we consider the problem of identifying co-evolving subnetworks given a pair of temporal networks, which aim to capture the evolution of molecules and their interactions over time. Although this problem shares some characteristics of the well-known network alignment problems, it differs from existing network alignment formulations as it seeks a mapping of the two network topologies that is invariant to temporal evolution of the given networks. This is a computationally challenging problem as it requires capturing not only similar topologies between two networks but also their similar evolution patterns.

RESULTS

We present an efficient algorithm, Tempo, for solving identifying co-evolving subnetworks with two given temporal networks. We formally prove the correctness of our method. We experimentally demonstrate that Tempo scales efficiently with the size of network as well as the number of time points, and generates statistically significant alignments-even when evolution rates of given networks are high. Our results on a human aging dataset demonstrate that Tempo identifies novel genes contributing to the progression of Alzheimer's, Huntington's and Type II diabetes, while existing methods fail to do so.

CONCLUSIONS

Studying temporal networks in general and human aging specifically using Tempo enables us to identify age related genes from non age related genes successfully. More importantly, Tempo takes the network alignment problem one huge step forward by moving beyond the classical static network models.

Multiscale effects of heating and cooling on genes and gene networks

Fluctuating environments such as changes in ambient temperature represent a fundamental challenge to life. Cells must protect gene networks that protect them from such stresses, making it difficult to understand how temperature affects gene network function in general. Here, we focus on single genes and small synthetic network modules to reveal four key effects of nonoptimal temperatures at different biological scales: (i) a cell fate choice between arrest and resistance, (ii) slower growth rates, (iii) Arrhenius reaction rates, and (iv) protein structure changes. We develop a multiscale computational modeling framework that captures and predicts all of these effects. These findings promote our understanding of how temperature affects living systems and enables more robust cellular engineering for real-world applications.

Most organisms must cope with temperature changes. This involves genes and gene networks both as subjects and agents of cellular protection, creating difficulties in understanding. Here, we study how heating and cooling affect expression of single genes and synthetic gene circuits in Saccharomyces cerevisiae. We discovered that nonoptimal temperatures induce a cell fate choice between stress resistance and growth arrest. This creates dramatic gene expression bimodality in isogenic cell populations, as arrest abolishes gene expression. Multiscale models incorporating population dynamics, temperature-dependent growth rates, and Arrhenius scaling of reaction rates captured the effects of cooling, but not those of heating in resistant cells. Molecular-dynamics simulations revealed how heating alters the conformational dynamics of the TetR repressor, fully explaining the experimental observations. Overall, nonoptimal temperatures induce a cell fate decision and corrupt gene and gene network function in computationally predictable ways, which may aid future applications of engineered microbes in nonstandard temperatures.

Modeling of signaling crosstalk-mediated drug resistance and its implications on drug combination

The efficacy of pharmacological perturbation to the signaling transduction network depends on the network topology. However, whether and how signaling dynamics mediated by crosstalk contributes to the drug resistance are not fully understood and remain to be systematically explored. In this study, motivated by a realistic signaling network linked by crosstalk between EGF/EGFR/Ras/MEK/ERK pathway and HGF/HGFR/PI3K/AKT pathway, we develop kinetic models for several small networks with typical crosstalk modules to investigate the role of the architecture of crosstalk in inducing drug resistance. Our results demonstrate that crosstalk inhibition diminishes the response of signaling output to the external stimuli. Moreover, we show that signaling crosstalk affects the relative sensitivity of drugs, and some types of crosstalk modules that could yield resistance to the targeted drugs were identified. Furthermore, we quantitatively evaluate the relative efficacy and synergism of drug combinations. For the modules that are resistant to the targeted drug, we identify drug targets that can not only increase the relative drug efficacy but also act synergistically. In addition, we analyze the role of the strength of crosstalk in switching a module between drug-sensitive and drug-resistant. Our study provides mechanistic insights into the signaling crosstalk-mediated mechanisms of drug resistance and provides implications for the design of synergistic drug combinations to reduce drug resistance.

Negative Regulation Gene Circuits for Efflux Pump Control

Synthetic biologists aim to design biological systems for a variety of industrial and medical applications, ranging from biofuel to drug production. Synthetic gene circuits regulating efflux pump protein expression can achieve this by driving desired substrates such as biofuels, pharmaceuticals, or other chemicals out of the cell in a precisely controlled manner. However, efflux pumps may introduce implicit negative feedback by pumping out intracellular inducer molecules that control gene circuits, which then can alter gene circuit function. Therefore, synthetic gene circuits must be carefully designed and constructed for precise efflux control. Here, we provide protocols for quantitatively modeling and building synthetic gene constructs for efflux pump regulation.

Stochasticity in the Genotype-Phenotype Map: Implications for the Robustness and Persistence of Bet-Hedging

Nongenetic variation in phenotypes, or bet-hedging, has been observed as a driver of drug resistance in both bacterial infections and cancers. Here, we study how bet-hedging emerges in genotype-phenotype (GP) mapping through a simple interaction model: a molecular switch. We use simple chemical reaction networks to implement stochastic switches that map gene products to phenotypes, and investigate the impact of structurally distinct mappings on the evolution of phenotypic heterogeneity. Bet-hedging naturally emerges within this model, and is robust to evolutionary loss through mutations to both the expression of individual genes, and to the network itself. This robustness explains an apparent paradox of bet-hedging-why does it persist in environments where natural selection necessarily acts to remove it? The structure of the underlying molecular mechanism, itself subject to selection, can slow the evolutionary loss of bet-hedging to ensure a survival mechanism against environmental catastrophes even when they are rare. Critically, these properties, taken together, have profound implications for the use of treatment-holidays to combat bet-hedging-driven resistant disease, as the efficacy of breaks from treatment will ultimately be determined by the structure of the GP mapping.

Identification of large disjoint motifs in biological networks

BACKGROUND

Biological networks provide great potential to understand how cells function. Network motifs, frequent topological patterns, are key structures through which biological networks operate. Finding motifs in biological networks remains to be computationally challenging task as the size of the motif and the underlying network grow. Often, different copies of a given motif topology in a network share nodes or edges. Counting such overlapping copies introduces significant problems in motif identification.

RESULTS

In this paper, we develop a scalable algorithm for finding network motifs. Unlike most of the existing studies, our algorithm counts independent copies of each motif topology. We introduce a set of small patterns and prove that we can construct any larger pattern by joining those patterns iteratively. By iteratively joining already identified motifs with those patterns, our algorithm avoids (i) constructing topologies which do not exist in the target network (ii) repeatedly counting the frequency of the motifs generated in subsequent iterations. Our experiments on real and synthetic networks demonstrate that our method is significantly faster and more accurate than the existing methods including SUBDUE and FSG.

CONCLUSIONS

We conclude that our method for finding network motifs is scalable and computationally feasible for large motif sizes and a broad range of networks with different sizes and densities. We proved that any motif with four or more edges can be constructed as a join of the small patterns.

Efflux Pump Control Alters Synthetic Gene Circuit Function

Synthetic biology aims to design new biological systems for predefined purposes, such as the controlled secretion of biofuels, pharmaceuticals, or other chemicals. Synthetic gene circuits regulating an efflux pump from the ATP-binding cassette (ABC) protein family could achieve this. However, ABC efflux pumps can also drive out intracellular inducer molecules that control the gene circuits. This will introduce an implicit feedback that could alter gene circuit function in ways that are poorly understood. Here, we used two synthetic gene circuits inducible by tetracycline family molecules to regulate the expression of a yeast ABC pump (Pdr5p) that pumps out the inducer. Pdr5p altered the dose-responses of the original gene circuits substantially in Saccharomyces cerevisiae. While one aspect of the change could be attributed to the efflux pumping function of Pdr5p, another aspect remained unexplained. Quantitative modeling indicated that reduced regulator gene expression in addition to efflux pump function could fully explain the altered dose-responses. These predictions were validated experimentally. Overall, we highlight how efflux pumps can alter gene circuit dynamics and demonstrate the utility of mathematical modeling in understanding synthetic gene circuit function in new circumstances.

Noise and low-level dynamics can coordinate multicomponent bet hedging mechanisms

To counter future uncertainty, cells can stochastically express stress response mechanisms to diversify their population and hedge against stress. This approach allows a small subset of the population to survive without the prohibitive cost of constantly expressing resistance machinery at the population level. However, expression of multiple genes in concert is often needed to ensure survival, requiring coordination of infrequent events across many downstream targets. This raises the question of how cells orchestrate the timing of multiple rare events without adding cost. To investigate this, we used a stochastic model to study regulation of downstream target genes by a transcription factor. We compared several upstream regulator profiles, including constant expression, pulsatile dynamics, and noisy expression. We found that pulsatile dynamics and noise are sufficient to coordinate expression of multiple downstream genes. Notably, this is true even when fluctuations in the upstream regulator are far below the dissociation constants of the regulated genes, as with infrequently activated genes. As an example, we simulated the dynamics of the multiple antibiotic resistance activator (MarA) and 40 diverse downstream genes it regulates, determining that low-level dynamics in MarA are sufficient to coordinate expression of resistance mechanisms. We also demonstrated that noise can play a similar coordinating role. Importantly, we found that these benefits are present without a corresponding increase in the population-level cost. Therefore, our model suggests that low-level dynamics or noise in a transcription factor can coordinate expression of multiple stress response mechanisms by engaging them simultaneously without adding to the overall cost.

Effect and evolution of gene expression noise on the fitness landscape

Gene expression is a stochastic process that affects cellular and population fitness. Noise in gene expression can enhance fitness by increasing cell to cell variability as well as the time cells spend in favorable expression states. Using a stochastic model of gene expression together with a fitness function that incorporates the costs and benefits of gene expression in a stressful environment, we show that the fitness landscape is shaped by gene expression noise in more complex ways than previously anticipated. We find that mutations modulating the properties of expression noise enable cell populations to optimize their position on the fitness landscape. Additionally, we find that low levels of expression noise evolve under conditions where the fitness benefits of expression exceed the fitness costs, and that high levels of expression noise evolve when the expression costs exceed the fitness benefits. The results presented in this study expand our understanding of the interplay between stochastic gene expression and fitness in selective environments.

Role of relaxation time scale in noisy signal transduction

Intra-cellular fluctuations, mainly triggered by gene expression, are an inevitable phenomenon observed in living cells. It influences generation of phenotypic diversity in genetically identical cells. Such variation of cellular components is beneficial in some contexts but detrimental in others. To quantify the fluctuations in a gene product, we undertake an analytical scheme for studying few naturally abundant linear as well as branched chain network motifs. We solve the Langevin equations associated with each motif under the purview of linear noise approximation and derive the expressions for Fano factor and mutual information in close analytical form. Both quantifiable expressions exclusively depend on the relaxation time (decay rate constant) and steady state population of the network components. We investigate the effect of relaxation time constraints on Fano factor and mutual information to indentify a time scale domain where a network can recognize the fluctuations associated with the input signal more reliably. We also show how input population affects both quantities. We extend our calculation to long chain linear motif and show that with increasing chain length, the Fano factor value increases but the mutual information processing capability decreases. In this type of motif, the intermediate components act as a noise filter that tune up input fluctuations and maintain optimum fluctuations in the output. For branched chain motifs, both quantities vary within a large scale due to their network architecture and facilitate survival of living system in diverse environmental conditions.

Bacterial temporal dynamics enable optimal design of antibiotic treatment

There is a critical need to better use existing antibiotics due to the urgent threat of antibiotic resistant bacteria coupled with the reduced effort in developing new antibiotics. β-lactam antibiotics represent one of the most commonly used classes of antibiotics to treat a broad spectrum of Gram-positive and -negative bacterial pathogens. However, the rise of extended spectrum β-lactamase (ESBL) producing bacteria has limited the use of β-lactams. Due to the concern of complex drug responses, many β-lactams are typically ruled out if ESBL-producing pathogens are detected, even if these pathogens test as susceptible to some β-lactams. Using quantitative modeling, we show that β-lactams could still effectively treat pathogens producing low or moderate levels of ESBLs when administered properly. We further develop a metric to guide the design of a dosing protocol to optimize treatment efficiency for any antibiotic-pathogen combination. Ultimately, optimized dosing protocols could allow reintroduction of a repertoire of first-line antibiotics with improved treatment outcomes and preserve last-resort antibiotics.

Antibiotic resistance is a growing problem that the World Health Organization describes as “one of the top three threats to global health.” To date, bacteria have developed resistance to all antibiotics used in clinical settings. Unfortunately, the evolution of antibiotic resistant bacteria is accelerating, as antibiotics continue to be misused and overused. As the antibiotic pipeline is drying up, it becomes increasingly critical to utilize the antibiotics already on the market more effectively. The key to designing better regimens lies in the ability to predict how bacteria will respond to a particular antibiotic treatment. Because of this, we need a simple metric that characterizes this pathogen-antibiotic interaction that can be easily measured and used to design dosing protocols that will effectively clear an infection. To help guide the design of effective protocols, we use quantitative modeling to develop a metric that is easy to measure and quantifies the pathogen-antibiotic interaction. Through optimized antibiotic regimens, our strategy could extend the use of first-line antibiotics, improve treatment outcome, and preserve last-resort antibiotics.