Suppressive might of a few: T follicular regulatory cells impede auto-reactivity despite being outnumbered in the germinal centres.

The selection of high-affinity B cells and the production of high-affinity antibodies are mediated by T follicular helper cells (Tfhs) within germinal centres (GCs). Therein, somatic hypermutation and selection enhance B cell affinity but risk the emergence of self-reactive B cell clones. Despite being outnumbered compared to their helper counterpart, the ablation of T follicular regulatory cells (Tfrs) results in enhanced dissemination of self-reactive antibody-secreting cells (ASCs). The specific mechanisms by which Tfrs exert their regulatory action on self-reactive B cells are largely unknown. We developed computer simulations to investigate how Tfrs regulate either selection or differentiation of B cells to prevent auto-reactivity. We observed that Tfr-induced apoptosis of self-reactive B cells during the selection phase impedes self-reactivity with physiological Tfr numbers, especially when Tfrs can access centrocyte-enriched GC areas. While this aided in selecting non-self-reactive B cells by restraining competition, higher Tfr numbers distracted non-self-reactive B cells from receiving survival signals from Tfhs. Thus, the location and number of Tfrs must be regulated to circumvent such Tfr distraction and avoid disrupting GC evolution. In contrast, when Tfrs regulate differentiation of selected centrocytes by promoting recycling to the dark zone phenotype of self-reactive GC resident pre-plasma cells (GCPCs), higher Tfr numbers were required to impede the circulation of self-reactive ASCs (s-ASCs). On the other hand, Tfr-engagement with GCPCs and subsequent apoptosis of s-ASCs can control self-reactivity with low Tfr numbers, but does not confer selection advantage to non-self-reactive B cells. The simulations predict that to restrict auto-reactivity, natural redemption of self-reactive B cells is insufficient and that Tfrs should increase the mutation probability of self-reactive B cells.

Dynamics and future prediction of LULC on Pare River basin of Arunachal Pradesh using machine learning techniques.

Anthropogenic disturbances caused by increasing population densities are a significant concern as they accelerate climate change. Thus, regular monitoring of land use/land cover (LULC) is essential to mitigate these effects. Pare River basin of Arunachala Pradesh situated in the foothills of Eastern Himalayas was selected for this study. Landsat-5 TM and Landsat-8 OLI data from 2000 (T1), 2015 (T2), and 2020 (T3) were used to prepare the LULC map. A support vector machine (SVM) classifier in the Google Earth Engine (GEE) environment was utilized for classification of LULC, while the TerrSet software environment was used for change analysis and projection using the CA-MC model. The SVM classifier produced overall all classification accuracies of 0.91, 0.85, and 0.91 with kappa values of 0.88, 0.82, and 0.89 for T1, T2, and T3, respectively. The CA-MC model, which combines Markov chain and hybrid cellular automata, was calibrated with various predictor variables, including natural, proximity, and demographic variables along with T1 and T2 LULC and validated using T3 LULC. The MLP was used for calibration, and an accuracy rate of above 0.70 was employed to generate transition potential maps (TPMs). The TPMs were used to project future LULC for 2030, 2040, and 2050. Validation analysis produced satisfactory results, with Kno, Klocation, Kquality, and Kstandard values of 0.96, 0.95, 0.95, and 0.93, respectively. Receiver operating characteristics (ROC) analysis showed an excellent area under the curve (AUC) value of 0.87. The findings of this study provide important insights to decision-makers and stakeholders in addressing the impacts of LULC changes.

Amount of antigen, T follicular helper cells and affinity of founder cells shape the diversity of germinal center B cells: A computational study.

A variety of B cell clones seed the germinal centers, where a selection stringency expands the fitter clones to generate higher affinity antibodies. However, recent experiments suggest that germinal centers often retain a diverse set of B cell clones with a range of affinities and concurrently carry out affinity maturation. Amid a tendency to flourish germinal centers with fitter clones, how several B cell clones with differing affinities can be concurrently selected remains poorly understood. Such a permissive selection may allow non-immunodominant clones, which are often rare and of low-affinity, to somatically hypermutate and result in a broad and diverse B cell response. How the constituent elements of germinal centers, their quantity and kinetics may modulate diversity of B cells, has not been addressed well. By implementing a state-of-the-art agent-based model of germinal center, here, we study how these factors impact temporal evolution of B cell clonal diversity and its underlying balance with affinity maturation. While we find that the extent of selection stringency dictates clonal dominance, limited antigen availability on follicular dendritic cells is shown to expedite the loss of diversity of B cells as germinal centers mature. Intriguingly, the emergence of a diverse set of germinal center B cells depends on high affinity founder cells. Our analysis also reveals a substantial number of T follicular helper cells to be essential in balancing affinity maturation with clonal diversity, as a low number of T follicular helper cells impedes affinity maturation and also contracts the scope for a diverse B cell response. Our results have implications for eliciting antibody responses to non-immunodominant specificities of the pathogens by controlling the regulators of the germinal center reaction, thereby pivoting a way for vaccine development to generate broadly protective antibodies.

Tingible body macrophages arise from lymph node-resident precursors and uptake B cells by dendrites.

Antibody affinity maturation depends on the formation of germinal centers (GCs) in lymph nodes. This process generates a massive number of apoptotic B cells, which are removed by a specialized subset of phagocytes, known as tingible body macrophages (TBMs). Although defects in these cells are associated with pathological conditions, the identity of their precursors and the dynamics of dying GC B cell disposal remained unknown. Here, we demonstrate that TBMs originate from pre-existing lymph node-resident precursors that enter the lymph node follicles in a GC-dependent manner. Intravital imaging shows that TBMs are stationary cells that selectively phagocytose GC B cells via highly dynamic protrusions and accommodate the final stages of B cell apoptosis. Cell-specific depletion and chimeric mouse models revealed that GC B cells drive TBM formation from bone marrow-derived precursors stationed within lymphoid organs prior to the immune challenge. Understanding TBM dynamics and function may explain the emergence of various antibody-mediated autoimmune conditions.

BTG1 mutation yields supercompetitive B cells primed for malignant transformation.

Multicellular life requires altruistic cooperation between cells. The adaptive immune system is a notable exception, wherein germinal center B cells compete vigorously for limiting positive selection signals. Studying primary human lymphomas and developing new mouse models, we found that mutations affecting BTG1 disrupt a critical immune gatekeeper mechanism that strictly limits B cell fitness during antibody affinity maturation. This mechanism converted germinal center B cells into supercompetitors that rapidly outstrip their normal counterparts. This effect was conferred by a small shift in MYC protein induction kinetics but resulted in aggressive invasive lymphomas, which in humans are linked to dire clinical outcomes. Our findings reveal a delicate evolutionary trade-off between natural selection of B cells to provide immunity and potentially dangerous features that recall the more competitive nature of unicellular organisms.

Impact of COVID-19 on ports, multimodal logistics and transport sector in India: Responses and policy imperatives.

The COVID-19 pandemic spread rapidly in 2020 and led to full or partial lockdowns worldwide. The restrictions led to most economies contracting in the April-June quarter of 2020. As per the IMF World Economic Outlook of June 2020, the global economy was expected to contract by 4.9% in 2020, whereas the Indian economy was expected to contract by 4.5%, affecting a population of more than 1.3 billion. This led to disruption in the global supply chain network, adversely hitting all modes of transport, including ports and multimodal logistics. Based on the survey results of 98 respondents in the Indian ports, multimodal logistics and transport ("PMT") sector conducted during the lockdown in May-June 2020, an adverse impact on profitability, employment, operations and capital expenditure was identified, which is consistent with observations from previous economic crises. Additionally, we have estimated the financial impact on different categories of organisations in the industry. Based on these impacts, this study identified key factors required to support the sector and build resilience for the future through technology interventions, developing multimodal transport systems, and providing financial support and employment support for small and medium enterprises (SMEs). By providing an analysis of the impact of COVID-19 on the different categories of organisations based on data collected during and shortly after the lockdown, the study makes a novel contribution to understanding the impact of such crises on the PMT industry at a granular level and the aggregate-level impacts with guidance for policy and managers and in particular, in an emerging economy context.

Testing and isolation to prevent overloaded healthcare facilities and reduce death rates in the SARS-CoV-2 pandemic in Italy.

Background: During the first wave of COVID-19, hospital and intensive care unit beds got overwhelmed in Italy leading to an increased death burden. Based on data from Italian regions, we disentangled the impact of various factors contributing to the bottleneck situation of healthcare facilities, not well addressed in classical SEIR-like models. A particular emphasis was set on the undetected fraction (dark figure), on the dynamically changing hospital capacity, and on different testing, contact tracing, quarantine strategies. Methods: We first estimated the dark figure for different Italian regions. Using parameter estimates from literature and, alternatively, with parameters derived from a fit to the initial phase of COVID-19 spread, the model was optimized to fit data (infected, hospitalized, ICU, dead) published by the Italian Civil Protection. Results: We show that testing influenced the infection dynamics by isolation of newly detected cases and subsequent interruption of infection chains. The time-varying reproduction number (R t) in high testing regions decreased to <1 earlier compared to the low testing regions. While an early test and isolate (TI) scenario resulted in up to ~31% peak reduction of hospital occupancy, the late TI scenario resulted in an overwhelmed healthcare system. Conclusions: An early TI strategy would have decreased the overall hospital usage drastically and, hence, death toll (∼34% reduction in Lombardia) and could have mitigated the lack of healthcare facilities in the course of the pandemic, but it would not have kept the hospitalization amount within the pre-pandemic hospital limit.

Inland heat waves (IHWs) and associated impacts on hydro-biology of aquatic ecosystems in lower Ganga basin, India.

The present study assessed the occurrence and impact of heat waves on the ecology of two ecosystems namely Bhomra wetland and Ganga River stretch, India. The water samples collected from these ecosystems were analyzed for estimating the hydrological and biological variables during heat wave. The inland heat index (IHI) was derived from the climatic variables, relative humidity and temperature. The study indicated the predominant and periodic occurrence of inland heat waves (IHW) with indices ranging from 34.8 to 42.8 °C and 35.9 to 43.5 °C at the Bhomra and Ganga River stretch respectively during the summer months (March-June). The first two components of the principal component analysis of physico-chemical parameters and heat index explained 45.6% and 59% of the variation in the Bhomra and Ganga River stretch respectively. PCA showed a similar pattern in variation of IHWs and dissolved oxygen, nutrients, hardness and alkalinity, but a distinct pattern with conductivity and TDS in the wetland. IHW exhibited a similar pattern of variation with TDS, conductivity, dissolved oxygen, pH and hardness and distinct pattern with alkalinity, phosphate and nitrate in the river stretch. The first two components of PCA of IHI and plankton abundance explained 89% of the variation and IHI had a similar pattern of variation with the abundance of diatoms and a diverse pattern of variation with blue-green and green algae in the studied ecosystems which might affect the food availability of the associated fishes. The study suggests that IHW influences the water quality and primary producers and also summarizes the impact of IHW on ecosystem services and necessitates mitigation measures.

Sepsis induces interleukin 6, gp130/JAK2/STAT3, and muscle wasting.

BACKGROUND: Sepsis and inflammation can cause intensive care unit-acquired weakness (ICUAW). Increased interleukin-6 (IL-6) plasma levels are a risk factor for ICUAW. IL-6 signalling involves the glycoprotein 130 (gp130) receptor and the JAK/STAT-pathway, but its role in sepsis-induced muscle wasting is uncertain. In a clinical observational study, we found that the IL-6 target gene, SOCS3, was increased in skeletal muscle of ICUAW patients indicative for JAK/STAT-pathway activation. We tested the hypothesis that the IL-6/gp130-pathway mediates ICUAW muscle atrophy. METHODS: We sequenced RNA (RNAseq) from tibialis anterior (TA) muscle of cecal ligation and puncture-operated (CLP) and sham-operated wildtype (WT) mice. The effects of the IL-6/gp130/JAK2/STAT3-pathway were investigated by analysing the atrophy phenotype, gene expression, and protein contents of C2C12 myotubes. Mice lacking Il6st, encoding gp130, in myocytes (cKO) and WT controls, as well as mice treated with the JAK2 inhibitor AG490 or vehicle were exposed to CLP or sham surgery for 24 or 96 h. RESULTS: Analyses of differentially expressed genes in RNAseq (≥2-log2-fold change, P < 0.01) revealed an activation of IL-6-signalling and JAK/STAT-signalling pathways in muscle of septic mice, which occurred after 24 h and lasted at least for 96 h during sepsis. IL-6 treatment of C2C12 myotubes induced STAT3 phosphorylation (three-fold, P < 0.01) and Socs3 mRNA expression (3.1-fold, P < 0.01) and caused myotube atrophy. Knockdown of Il6st diminished IL-6-induced STAT3 phosphorylation (-30.0%; P < 0.01), Socs3 mRNA expression, and myotube atrophy. JAK2 (- 29.0%; P < 0.01) or STAT3 inhibition (-38.7%; P < 0.05) decreased IL-6-induced Socs3 mRNA expression. Treatment with either inhibitor attenuated myotube atrophy in response to IL-6. CLP-operated septic mice showed an increased STAT3 phosphorylation and Socs3 mRNA expression in TA muscle, which was reduced in septic Il6st-cKO mice by 67.8% (P < 0.05) and 85.6% (P < 0.001), respectively. CLP caused a loss of TA muscle weight, which was attenuated in Il6st-cKO mice (WT: -22.3%, P < 0.001, cKO: -13.5%, P < 0.001; WT vs. cKO P < 0.001). While loss of Il6st resulted in a reduction of MuRF1 protein contents, Atrogin-1 remained unchanged between septic WT and cKO mice. mRNA expression of Trim63/MuRF1 and Fbxo32/Atrogin-1 were unaltered between CLP-treated WT and cKO mice. AG490 treatment reduced STAT3 phosphorylation (-22.2%, P < 0.05) and attenuated TA muscle atrophy in septic mice (29.6% relative reduction of muscle weight loss, P < 0.05). The reduction in muscle atrophy was accompanied by a reduction in Fbxo32/Atrogin-1-mRNA (-81.3%, P < 0.05) and Trim63/MuRF1-mRNA expression (-77.6%, P < 0.05) and protein content. CONCLUSIONS: IL-6 via the gp130/JAK2/STAT3-pathway mediates sepsis-induced muscle atrophy possibly contributing to ICUAW.

High SARS-CoV-2 seroprevalence in children and adults in the Austrian ski resort of Ischgl.

BACKGROUND: In early March 2020, a SARS-CoV-2 outbreak in the ski resort Ischgl in Austria initiated the spread of SARS-CoV-2 throughout Austria and Northern Europe. METHODS: Between April 21st and 27th 2020, a cross-sectional epidemiologic study targeting the full population of Ischgl (n = 1867), of which 79% could be included (n = 1473, incl. 214 children), was performed. For each individual, the study involved a SARS-CoV-2 PCR, antibody testing and structured questionnaires. A mathematical model was used to help understand the influence of the determined seroprevalence on virus transmission. RESULTS: The seroprevalence was 42.4% (95% confidence interval (CI) 39.8-44.7). Individuals under 18 showed a significantly lower seroprevalence of 27.1% (95% CI 21.3-33.6) than adults (45%; 95% CI 42.2-47.7; OR of 0.455, 95% CI 0.356-0.682, p < 0.001). Of the seropositive individuals, 83.7% had not been diagnosed to have had SARS-CoV-2 infection previously. The clinical course was generally mild. Over the previous two months, two COVID-19-related deaths had been recorded, corresponding to an infection fatality rate of 0.25% (95% CI 0.03-0.91). Only 8 (0.5 %) individuals were newly diagnosed to be infected with SARS-CoV-2 during this study. CONCLUSIONS: Ischgl was hit early and hard by SARS-CoV-2 leading to a high local seroprevalence of 42.4%, which was lower in individuals below the age of 18 than in adults. Mathematical modeling suggests that a drastic decline of newly infected individuals in Ischgl by the end of April occurred due to the dual impact from the non-pharmacological interventions and a high immunization of the Ischgl population.

Membrane composition and lipid to protein ratio modulate amyloid kinetics of yeast prion protein.

Understanding of prion aggregation in a membrane environment may help to ameliorate neurodegenerative complications caused by the amyloid forms of prions. Here, we investigated the membrane binding-induced aggregation of yeast prion protein Sup35. Using the combination of fluorescence correlation spectroscopy (FCS) at single molecule resolution and other biophysical studies, we establish that lipid composition and lipid/protein ratio are key modulators of the aggregation kinetics of Sup35. In the presence of a zwitterionic membrane (DMPC), Sup35 exhibited novel biphasic aggregation kinetics at lipid/protein ratios ranging between 20 : 1 and 70 : 1 (termed here as the optimum lipid concentration, OLC). In ratios below (low lipid concentration, LLC) and above (ELC, excess lipid concentration) that range, the aggregation was found to be monophasic. In contrast, in the presence of negatively charged membranes, we did not observe any bi-phasic aggregation kinetics in the entire range of protein to lipid ratios. Our results provide a mechanistic description of the role that membrane concentration/composition-modulated aggregation may play in neurodegenerative diseases.

Development of the reproduction number from coronavirus SARS-CoV-2 case data in Germany and implications for political measures.

BACKGROUND: SARS-CoV-2 has induced a worldwide pandemic and subsequent non-pharmaceutical interventions (NPIs) to control the spread of the virus. As in many countries, the SARS-CoV-2 pandemic in Germany has led to a consecutive roll-out of different NPIs. As these NPIs have (largely unknown) adverse effects, targeting them precisely and monitoring their effectiveness are essential. We developed a compartmental infection dynamics model with specific features of SARS-CoV-2 that allows daily estimation of a time-varying reproduction number and published this information openly since the beginning of April 2020. Here, we present the transmission dynamics in Germany over time to understand the effect of NPIs and allow adaptive forecasts of the epidemic progression. METHODS: We used a data-driven estimation of the evolution of the reproduction number for viral spreading in Germany as well as in all its federal states using our model. Using parameter estimates from literature and, alternatively, with parameters derived from a fit to the initial phase of COVID-19 spread in different regions of Italy, the model was optimized to fit data from the Robert Koch Institute. RESULTS: The time-varying reproduction number (Rt) in Germany decreased to <1 in early April 2020, 2-3 weeks after the implementation of NPIs. Partial release of NPIs both nationally and on federal state level correlated with moderate increases in Rt until August 2020. Implications of state-specific Rt on other states and on national level are characterized. Retrospective evaluation of the model shows excellent agreement with the data and usage of inpatient facilities well within the healthcare limit. While short-term predictions may work for a few weeks, long-term projections are complicated by unpredictable structural changes. CONCLUSIONS: The estimated fraction of immunized population by August 2020 warns of a renewed outbreak upon release of measures. A low detection rate prolongs the delay reaching a low case incidence number upon release, showing the importance of an effective testing-quarantine strategy. We show that real-time monitoring of transmission dynamics is important to evaluate the extent of the outbreak, short-term projections for the burden on the healthcare system, and their response to policy changes.

Genomic Landscape and Targeted Treatment of Gallbladder Cancer: Results of a First Ongoing Prospective Study.

Background Prognosis of gallbladder cancer (GBC) has not changed in the past 20 years. Comprehensive genomic profiling (CGP) carries potential to determine the actionability for multiple targets, including ERBB2 , ERBB3 , MET , ROSI , FGFR , and PIK3 . This study evaluates the role of CGP and targeted therapies. Methods This is a multicenter, prospective, single-arm study. All consecutive patients of unresectable and/or metastatic GBC of age ≥18 years were enrolled. Hybrid capture-based CGP was performed by Foundation Medicine CDx. All patients received first-line chemotherapy with gemcitabine-cisplatin regimen. Patients with ERBB2/3 amplification received trastuzumab with capecitabine or nab-paclitaxel, and patients with MET amplification were treated with crizotinib. For ERBB2/3 mutations, lapatinib plus capecitabine regimen was used. Results Fifty patients were studied with a median age of 56 years (range 26-83) and a male-to-female ratio of 1:1.6. ERBB2 and ERBB3 amplification was seen in 9 (18%) and 2 (4%) patients, respectively. Four patients with ERBB2 amplification received trastuzumab and/or lapatinib, showed partial response, and maintained response beyond 12 weeks. One patient had mixed response, whereas two patients progressed on trastuzumab and lapatinib. Three patients with ERBB3 mutations showed response to lapatinib-capecitabine. One patient with MET amplification responded to crizotinib for 4 weeks. PIK3 mutations were present in 14% of cases and were independent of ERBB aberrations. Conclusion GBC is enriched in 28% of patients with ERBB2 and ERBB3 amplifications and/or mutations. Responses are seen with lapatinib in concurrent ERBB2 mutation and amplification. ERBB3 mutation showed response to lapatinib. MET and PIK3 are new findings in GBC, which may be targeted.

Osmolytes ameliorate the effects of stress in the absence of the heat shock protein Hsp104 in Saccharomyces cerevisiae.

Aggregation of the prion protein has strong implications in the human prion disease. Sup35p is a yeast prion, and has been used as a model protein to study the disease mechanism. We have studied the pattern of Sup35p aggregation inside live yeast cells under stress, by using confocal microscopy, fluorescence activated cell sorting and western blotting. Heat shock proteins are a family of proteins that are produced by yeast cells in response to exposure to stressful conditions. Many of the proteins behave as chaperones to combat stress-induced protein misfolding and aggregation. In spite of this, yeast also produce small molecules called osmolytes during stress. In our work, we tried to find the reason as to why yeast produce osmolytes and showed that the osmolytes are paramount to ameliorate the long-term effects of lethal stress in Saccharomyces cerevisiae, either in the presence or absence of Hsp104p.

Lineage space and the propensity of bacterial cells to undergo growth transitions.

The molecular makeup of the offspring of a dividing cell gradually becomes phenotypically decorrelated from the parent cell by noise and regulatory mechanisms that amplify phenotypic heterogeneity. Such regulatory mechanisms form networks that contain thresholds between phenotypes. Populations of cells can be poised near the threshold so that a subset of the population probabilistically undergoes the phenotypic transition. We sought to characterize the diversity of bacterial populations around a growth-modulating threshold via analysis of the effect of non-genetic inheritance, similar to conditions that create antibiotic-tolerant persister cells and other examples of bet hedging. Using simulations and experimental lineage data in Escherichia coli, we present evidence that regulation of growth amplifies the dependence of growth arrest on cellular lineage, causing clusters of related cells undergo growth arrest in certain conditions. Our simulations predict that lineage correlations and the sensitivity of growth to changes in toxin levels coincide in a critical regime. Below the critical regime, the sizes of related growth arrested clusters are distributed exponentially, while in the critical regime clusters sizes are more likely to become large. Furthermore, phenotypic diversity can be nearly as high as possible near the critical regime, but for most parameter values it falls far below the theoretical limit. We conclude that lineage information is indispensable for understanding regulation of cellular growth.

WITHDRAWN: Curcumin shifts RAS-induced pro-proliferative MEK/ERK-signaling toward pro-apoptotic p38MAPK/JNK1-signaling, triggering p53 activation and apoptosis.

This article has been withdrawn by the authors. A mistake was made during the preparation of Fig 1C, NKE panel. The Western blot data shown for p-ERK1/2 and actin are not from this set, but rather a similar set of data from a different experiment. The authors apologize to the readers.

Role of functionality in two-component signal transduction: a stochastic study.

We present a stochastic formalism for signal transduction processes in a bacterial two-component system. Using elementary mass action kinetics, the proposed model takes care of signal transduction in terms of a phosphotransfer mechanism between the cognate partners of a two-component system, viz., the sensor kinase and the response regulator. Based on the difference in functionality of the sensor kinase, the noisy phosphotransfer mechanism has been studied for monofunctional and bifunctional two-component systems using the formalism of the linear noise approximation. Steady-state analysis of both models quantifies different physically realizable quantities, e.g., the variance, the Fano factor (variance/mean), and mutual information. The resultant data reveal that both systems reliably transfer information of extracellular environment under low external stimulus and in a high-kinase-and-phosphatase regime. We extend our analysis further by studying the role of the two-component system in downstream gene regulation.

Analysis of DevR regulated genes in Mycobacterium tuberculosis.

The DevRS two component system of Mycobacterium tuberculosis is responsible for its dormancy in host and becomes operative under hypoxic condition. It is experimentally known that phosphorylated DevR controls the expression of several downstream genes in a complex manner. In the present work we propose a theoretical model to show role of binding sites in DevR mediated gene expression. Individual and collective role of binding sites in regulating DevR mediated gene expression has been shown via modeling. Objective of the present work is twofold. First, to describe qualitatively the temporal dynamics of wild type genes and their known mutants. Based on these results we propose that DevR controlled gene expression follows a specific pattern which is efficient in describing other DevR mediated gene expression. Second, to analyze behavior of the system from information theoretical point of view. Using the tools of information theory we have calculated molecular efficiency of the system and have shown that it is close to the maximum limit of isothermal efficiency.

Quantification of noise in bifunctionality-induced post-translational modification.

We present a generic analytical scheme for the quantification of fluctuations due to bifunctionality-induced signal transduction within the members of a bacterial two-component system. The proposed model takes into account post-translational modifications in terms of elementary phosphotransfer kinetics. Sources of fluctuations due to autophosphorylation, kinase, and phosphatase activity of the sensor kinase have been considered in the model via Langevin equations, which are then solved within the framework of linear noise approximation. The resultant analytical expression of phosphorylated response regulators are then used to quantify the noise profile of biologically motivated single and branched pathways. Enhancement and reduction of noise in terms of extra phosphate outflux and influx, respectively, have been analyzed for the branched system. Furthermore, the role of fluctuations of the network output in the regulation of a promoter with random activation-deactivation dynamics has been analyzed.

Positive feedback and temperature mediated molecular switch controls differential gene regulation in Bordetella pertussis.

Based on the phosphorelay kinetics operative within BvgAS two component system we propose a mathematical framework for signal transduction and gene regulation of phenotypic phases in Bordetella pertussis. The proposed model identifies a novel mechanism of transcriptional interference between two promoters present in the bvg locus. To understand the system behavior under elevated temperature, the developed model has been studied in two different ways. First, a quasi-steady state analysis has been carried out for the two component system, comprising of sensor BvgS and response regulator BvgA. The quasi-steady state analysis reveals temperature induced sharp molecular switch, leading to amplification in the output of BvgA. Accumulation of a large pool of BvgA thus results into differential regulation of the downstream genes, including the gene encoding toxin. Numerical integration of the full network kinetics is then carried out to explore time dependent behavior of different system components, that qualitatively capture the essential features of experimental results performed in vivo. Furthermore, the developed model has been utilized to study mutants that are impaired in their ability to phosphorylate the transcription factor, BvgA, of the signaling network.