Garde2020-Minimal model describing metabolic oscillations in Bacillus subtilis biofilms
Ontology highlight
ABSTRACT:
SUBMITTER: 
Ahmad Zyoud
PROVIDER: BIOMD0000000932 | BioModels | 2024-09-02
REPOSITORIES: BioModels
Dataset's files
Items per page: 5 1 - 4 of 4 |
Publications
Differential equation-based minimal model describing metabolic oscillations in <i>Bacillus subtilis</i> biofilms.
Garde Ravindra R Ibrahim Bashar B Kovács Ákos T ÁT Schuster Stefan S
Royal Society open science 20200205 2
Biofilms offer an excellent example of ecological interaction among bacteria. Temporal and spatial oscillations in biofilms are an emerging topic. In this paper, we describe the metabolic oscillations in <i>Bacillus subtilis</i> biofilms by applying the smallest theoretical chemical reaction system showing Hopf bifurcation proposed by Wilhelm and Heinrich in 1995. The system involves three differential equations and a single bilinear term. We specifically select parameters that are suitable for ...[more]
Similar Datasets
Project description:RNAseq analysis of budding yeast respiratory oscillations was performed at high temporal resolution (4 min, 24 samples), covering 2.5 cycles of the short period oscillation of the IFO 0233 strain of budding yeast. This and related systems were previously analyzed by microarrays and RNAseq, but the expression of noncoding transcripts has not been analyzed. The oscillations were recorded by online measurements and the system underwent a multi-generational transient, where periods first gradually decreased (from 0.7 h to 0.6 h), and then a bifurcation of system dynamics occurred. After this bifurcation periods were longer but the oscillation was unstable. The RNAseq samples were taken in the 2.5 cycles before this bifurcation occurred, i.e., the data provides a snapshot of a transient state of the respiratory oscillations.
2022-07-11 | E-MTAB-11901 | biostudies-arrayexpress
Project description:During commitment of a multipotent stem or progenitor cell to a particular lineage, a large number of genes alter their expression in a coordinated manner orchestrated by the gene regulatory network (GRN). The constraints imposed by the GRN govern how cells move in the high-dimensional gene expression state space and can be understood as a dynamical system in which phenotypic cell states (cell types) are attractors that stabilize the cell-type characteristic gene expression pattern against molecular noise. Despite insights from various theoretical models, it remains elusive how multipotent cells, when committing to a specific lineage, exit their attractor and enter a new distinct attractor. Here we show, using single-cell resolution monitoring of transcript patterns by qPCR that commitment of multipotent blood progenitor cells to either the erythroid or the myeloid lineage is preceded by a destabilization of the progenitorsâ?? attractor state and a slowing-down of relaxation of cells from outlier states, indicating a critical state transition (â??tipping pointâ??). The high-dimensionality of the system (many genes) and availability of individual trajectories of a large ensemble of systems (many cells) affords a novel signature for critical transition which can be predicted from theory: Decrease of correlation between cells and concomitant increase of correlation between genes as the cell population approaches the tipping point. Consistent with a destabilizing bifurcation that simultaneously opens access to the erythroid and myeloid attractors, differentiation signal for either lineage caused some cells to commit to the â??wrongâ?? fate; moreover providing conflicting signals resulted in a delayed decision at the bifurcation point that however was ultimately resolved by commitment to one fate. These results suggest that the theoretical framework of â??early-warning signsâ?? and critical transitions can be applied to ensembles of high-dimensional systems, offering a formal tool for analyzing single-cell omics data beyond current descriptive computational pattern recognition. Mouse blood progenitor cells (EML cell line) was exposed to EPO, IL-3/GM-CSF or a mixture of both cytokines and gene expression change was measured in sorted subpopulations wrt Sca1 progenitor surface marker expression. In total, there was 4 conditions (including control), three time points (including d0) and 20 samples (10 samples in duplicates) were analyzed. Two independent experiments were performed for each condition. The untreated progenitor cell population was used as control.
2016-09-01 | E-GEOD-70405 | biostudies-arrayexpress
Project description:Ravindra Garde, Bashar Ibrahim & Stefan Schuster. Extending the minimal model of metabolic oscillations in Bacillus subtilis biofilms. Scientific Reports 10, 1 (2020).
Biofilms are composed of microorganisms attached to a solid surface or floating on top of a liquid surface. They pose challenges in the field of medicine but can also have useful applications in industry. Regulation of biofilm growth is complex and still largely elusive. Oscillations are thought to be advantageous for biofilms to cope with nutrient starvation and chemical attacks. Recently, a minimal mathematical model has been employed to describe the oscillations in Bacillus subtilis biofilms. In this paper, we investigate four different modifications to that minimal model in order to better understand the oscillations in biofilms. Our first modification is towards making a gradient of metabolites from the center of the biofilm to the periphery. We find that it does not improve the model and is therefore, unnecessary. We then use realistic Michaelis-Menten kinetics to replace the highly simple mass-action kinetics for one of the reactions. Further, we use reversible reactions to mimic the diffusion in biofilms. As the final modification, we check the combined effect of using Michaelis-Menten kinetics and reversible reactions on the model behavior. We find that these two modifications alone or in combination improve the description of the biological scenario.
2024-09-02 | BIOMD0000001053 | BioModels
Project description:The bacterial cell wall has been a celebrated target for antibiotics and holds real promise as a target for the discovery of new chemical matter to surmount pervasive multi-drug resistance among pathogenic bacteria. While the walls of Gram-negative bacteria are composed primarily of peptidoglycan, those of Gram-positives are more substantial and contain, in addition, large amounts of the polymer teichoic acid, covalently attached to peptidoglycan. Wall teichoic acids are a diverse group of phosphate-rich, extracellular polysaccharides that have been largely regarded as ancillary cell surface components. Recently, wall teichoic acid was shown to be essential to the proper rod-shaped cell morphology of the prototype Gram-positive bacterium Bacillus subtilis and an important virulence factor for the human pathogen Staphylococcus aureus. Thus wall teichoic acid synthesis is an intriguing target for the development of new cell wall-active antibiotics. Nevertheless, recent studies have shown that the dispensability of genes encoding teichoic acid biosynthetic enzymes in both B. subtilis and S. aureus is paradoxical and complex. Here, we report here on the discovery of a promoter (PywaC), which is sensitive to lesions in teichoic acid synthesis. Using this promoter we developed a luminescent, cell-based, reporter system to take a chemical-genetic approach to understanding the complexity of wall teichoic acid biogenesis using a large collection of antibiotics of well characterized biological activity. Our results reveal surprising interactions among undecaprenol, peptidoglycan and teichoic acid biosynthesis that help explain the complexity of teichoic acid gene dispensability. Furthermore, the new reporter assay represents an exciting avenue for the discovery of novel antibacterial molecules that impinge broadly on Gram-positive bacterial cell wall biogenesis. Keywords: comparison between depleted and repleted tagD mutant
2008-03-01 | GSE7373 | GEO
Project description:During commitment of a multipotent stem or progenitor cell to a particular lineage, a large number of genes alter their expression in a coordinated manner orchestrated by the gene regulatory network (GRN). The constraints imposed by the GRN govern how cells move in the high-dimensional gene expression state space and can be understood as a dynamical system in which phenotypic cell states (cell types) are attractors that stabilize the cell-type characteristic gene expression pattern against molecular noise. Despite insights from various theoretical models, it remains elusive how multipotent cells, when committing to a specific lineage, exit their attractor and enter a new distinct attractor. Here we show, using single-cell resolution monitoring of transcript patterns by qPCR that commitment of multipotent blood progenitor cells to either the erythroid or the myeloid lineage is preceded by a destabilization of the progenitors’ attractor state and a slowing-down of relaxation of cells from outlier states, indicating a critical state transition (“tipping point”). The high-dimensionality of the system (many genes) and availability of individual trajectories of a large ensemble of systems (many cells) affords a novel signature for critical transition which can be predicted from theory: Decrease of correlation between cells and concomitant increase of correlation between genes as the cell population approaches the tipping point. Consistent with a destabilizing bifurcation that simultaneously opens access to the erythroid and myeloid attractors, differentiation signal for either lineage caused some cells to commit to the “wrong” fate; moreover providing conflicting signals resulted in a delayed decision at the bifurcation point that however was ultimately resolved by commitment to one fate. These results suggest that the theoretical framework of “early-warning signs” and critical transitions can be applied to ensembles of high-dimensional systems, offering a formal tool for analyzing single-cell omics data beyond current descriptive computational pattern recognition.
2016-09-01 | GSE70405 | GEO
Project description:It is widely accepted that repair of double-strand breaks in bacteria that either sporulate or that undergo extended periods of stationary phase relies not only on homologous recombination but also on a minimal nonhomologous end joining (NHEJ) system consisting of a dedicated multifunctional ATP-dependent DNA Ligase D (LigD) and the DNA-end-binding protein Ku. Bacillus subtilis is one of the bacterial members with a NHEJ system that contributes to genome stability during the stationary phase and germination of spores, having been characterized exclusively in vivo. Here, the in vitro analysis of the functional properties of the purified B. subtilis LigD (BsuLigD) and Ku (BsuKu) proteins is presented. The results show that the essential biochemical signatures exhibited by BsuLigD agree with its proposed function in NHEJ: i) inherent polymerization activity showing preferential insertion of NMPs, ii) specific recognition of the phosphate group at the downstream 5' end, iii) intrinsic ligase activity, iv) ability to promote realignments of the template and primer strands during elongation of mispaired 3' ends, and v) it is recruited to DNA by BsuKu that stimulates the inherent polymerization and ligase activities of the enzyme allowing it to deal with and to hold different and unstable DNA realignments.
| S-EPMC3656841 | biostudies-literature
Project description:An 8,883-bp mini-pXO1 plasmid containing a replicon from Bacillus anthracis pXO1 (181.6 kb) was identified by making large deletions in the original plasmid using a newly developed Cre-loxP system. Portions of the truncated mini-pXO1 were cloned into an Escherichia coli vector unable to replicate in B. anthracis. A 5.95-kb region encompassing three putative genes was identified as the minimal pXO1 fragment required for replication of the resulting recombinant shuttle plasmid (named pMR) in B. anthracis. Deletion analysis showed that the only genes essential for replication were the pXO1-14 and pXO1-16 genes, which are transcribed in opposite directions and encode predicted proteins of 66.5 and 67.1 kDa, respectively. The ORF14 protein contains a helix-turn-helix motif, while the ORF16 upstream region contains attributes of a theta-replicating plasmid origin of replication (Ori), namely, an exclusively A+T-containing segment, five 9-bp direct repeats, an inverted repeat, and a sigma(A)-dependent promoter for the putative replication initiator Rep protein (ORF16). Spontaneous mutations generated in the ORF14, ORF16, and Ori regions of pMR during PCR amplification produced a temperature-sensitive plasmid that is unable to replicate in B. anthracis at 37 degrees C. The efficacy of transformation of plasmid-free B. anthracis Ames and Sterne strains by the original pMR was approximately 10(3) CFU/microg, while Bacillus cereus strains 569 and ATCC 10987 were transformed with efficiencies of 10(4) and 10(2) CFU/microg, respectively. Around 95% of B. anthracis cells retained pMR after one round of sporulation and germination.
| S-EPMC2725576 | biostudies-literature
Project description:The p53 transcription factor is a regulator of key cellular processes including DNA repair, cell cycle arrest, and apoptosis. In this theoretical study, we investigate how the complex circuitry of the p53 network allows for stochastic yet unambiguous cell fate decision-making. The proposed Markov chain model consists of the regulatory core and two subordinated bistable modules responsible for cell cycle arrest and apoptosis. The regulatory core is controlled by two negative feedback loops (regulated by Mdm2 and Wip1) responsible for oscillations, and two antagonistic positive feedback loops (regulated by phosphatases Wip1 and PTEN) responsible for bistability. By means of bifurcation analysis of the deterministic approximation we capture the recurrent solutions (i.e., steady states and limit cycles) that delineate temporal responses of the stochastic system. Direct switching from the limit-cycle oscillations to the "apoptotic" steady state is enabled by the existence of a subcritical Neimark-Sacker bifurcation in which the limit cycle loses its stability by merging with an unstable invariant torus. Our analysis provides an explanation why cancer cell lines known to have vastly diverse expression levels of Wip1 and PTEN exhibit a broad spectrum of responses to DNA damage: from a fast transition to a high level of p53 killer (a p53 phosphoform which promotes commitment to apoptosis) in cells characterized by high PTEN and low Wip1 levels to long-lasting p53 level oscillations in cells having PTEN promoter methylated (as in, e.g., MCF-7 cell line).
2024-09-02 | BIOMD0000000943 | BioModels
Project description:Wolf2000 - Cellular interaction on glycolytic
oscillations in yeast
A two-cell model of glycolysis.
This model is described in the article:
Effect of cellular
interaction on glycolytic oscillations in yeast: a theoretical
investigation.
Wolf J, Heinrich R.
Biochem. J. 2000 Jan; 345 Pt 2:
321-334
Abstract:
On the basis of a detailed model of yeast glycolysis, the
effect of intercellular dynamics is analysed theoretically. The
model includes the main steps of anaerobic glycolysis, and the
production of ethanol and glycerol. Transmembrane diffusion of
acetaldehyde is included, since it has been hypothesized that
this substance mediates the interaction. Depending on the
kinetic parameter, the single-cell model shows both stationary
and oscillatory behaviour. This agrees with experimental data
with respect to metabolite concentrations and phase shifts. The
inclusion of intercellular coupling leads to a variety of
dynamical modes, such as synchronous oscillations, and
different kinds of asynchronous behavior. These oscillations
can co-exist, leading to bi- and tri-rhythmicity. The
corresponding parameter regions have been identified by a
bifurcation analysis. The oscillatory dynamics of synchronized
cell populations are investigated by calculating the phase
responses to acetaldehyde pulses. Simulations are performed
with respect to the synchronization of two subpopulations that
are oscillating out of phase before mixing. The effect of the
various process on synchronization is characterized
quantitatively. While continuous exchange of acetaldehyde might
synchronize the oscillations for appropriate sets of parameter
values, the calculated synchronization time is longer than that
observed experimentally. It is concluded either that addition
to the transmembrane exchange of acetaldehyde, other processes
may contribute to intercellular coupling, or that intracellular
regulator feedback plays a role in the acceleration of the
synchronization. for appropriate sets of parameter values, the
calculated synchronization time is longer than that observed
experimentally. It is concluded either that addition to the
transmembrane exchange of acetaldehyde, other processes may
contribute to intercellular coupling, or that intracellular
regulator feedback plays a role in the acceleration of the
synchronization.
This model is hosted on
BioModels Database
and identified by:
BIOMD0000000691.
To cite BioModels Database, please use:
Chelliah V et al. BioModels: ten-year
anniversary. Nucl. Acids Res. 2015, 43(Database
issue):D542-8.
To the extent possible under law, all copyright and related or
neighbouring rights to this encoded model have been dedicated to
the public domain worldwide. Please refer to
CC0
Public Domain Dedication for more information.
2024-09-02 | BIOMD0000000691 | BioModels
Project description:Glucocorticoids act through the glucocorticoid receptor (GR), a ligand activated transcription factor, to modulate gene expression profiles in target cells and tissues. The principal physiological glucocorticoid in humans, cortisol, is released in a pulsatile fashion from the adrenal gland, which results in oscillations in blood concentrations, with a period of 1 - 2 hours. The biological consequence of fluctuating cortisol concentrations has not been explored. To identify the role of cortisol concentration oscillations for target cell response a flow-through cell culture system was developed. Transcription profiling was performed in the absence and presence of cortisol. Cells were exposed to constant cortisol levels of 100ng/ml or 200 ng/ml or oscillating cortisol concentrations.
2010-07-10 | E-MEXP-2447 | biostudies-arrayexpress