Chromosomal Arrangement of Phosphorelay Genes Couples Sporulation and DNA Replication.

Genes encoding proteins in a common regulatory network are frequently located close to one another on the chromosome to facilitate co-regulation or couple gene expression to growth rate. Contrasting with these observations, here, we demonstrate a functional role for the arrangement of Bacillus subtilis sporulation network genes on opposite sides of the chromosome. We show that the arrangement of two sporulation network genes, one located close to the origin and the other close to the terminus, leads to a transient gene dosage imbalance during chromosome replication. This imbalance is detected by the sporulation network to produce cell-cycle coordinated pulses of the sporulation master regulator Spo0A∼P. This pulsed response allows cells to decide between sporulation and continued vegetative growth during each cell cycle spent in starvation. The simplicity of this coordination mechanism suggests that it may be widely applicable in a variety of gene regulatory and stress-response settings. VIDEO ABSTRACT.

Slowdown of growth controls cellular differentiation.

How can changes in growth rate affect the regulatory networks behavior and the outcomes of cellular differentiation? We address this question by focusing on starvation response in sporulating Bacillus subtilis We show that the activity of sporulation master regulator Spo0A increases with decreasing cellular growth rate. Using a mathematical model of the phosphorelay-the network controlling Spo0A-we predict that this increase in Spo0A activity can be explained by the phosphorelay protein accumulation and lengthening of the period between chromosomal replication events caused by growth slowdown. As a result, only cells growing slower than a certain rate reach threshold Spo0A activity necessary for sporulation. This growth threshold model accurately predicts cell fates and explains the distribution of sporulation deferral times. We confirm our predictions experimentally and show that the concentration rather than activity of phosphorelay proteins is affected by the growth slowdown. We conclude that sensing the growth rates enables cells to indirectly detect starvation without the need for evaluating specific stress signals.

Role of Autoregulation and Relative Synthesis of Operon Partners in Alternative Sigma Factor Networks.

Despite the central role of alternative sigma factors in bacterial stress response and virulence their regulation remains incompletely understood. Here we investigate one of the best-studied examples of alternative sigma factors: the σB network that controls the general stress response of Bacillus subtilis to uncover widely relevant general design principles that describe the structure-function relationship of alternative sigma factor regulatory networks. We show that the relative stoichiometry of the synthesis rates of σB, its anti-sigma factor RsbW and the anti-anti-sigma factor RsbV plays a critical role in shaping the network behavior by forcing the σB network to function as an ultrasensitive negative feedback loop. We further demonstrate how this negative feedback regulation insulates alternative sigma factor activity from competition with the housekeeping sigma factor for RNA polymerase and allows multiple stress sigma factors to function simultaneously with little competitive interference.

Ultrasensitivity of the Bacillus subtilis sporulation decision.

Starving Bacillus subtilis cells execute a gene expression program resulting in the formation of stress-resistant spores. Sporulation master regulator, Spo0A, is activated by a phosphorelay and controls the expression of a multitude of genes, including the forespore-specific sigma factor σ(F) and the mother cell-specific sigma factor σ(E). Identification of the system-level mechanism of the sporulation decision is hindered by a lack of direct control over Spo0A activity. This limitation can be overcome by using a synthetic system in which Spo0A activation is controlled by inducing expression of phosphorelay kinase KinA. This induction results in a switch-like increase in the number of sporulating cells at a threshold of KinA. Using a combination of mathematical modeling and single-cell microscopy, we investigate the origin and physiological significance of this ultrasensitive threshold. The results indicate that the phosphorelay is unable to achieve a sufficiently fast and ultrasensitive response via its positive feedback architecture, suggesting that the sporulation decision is made downstream. In contrast, activation of σ(F) in the forespore and of σ(E) in the mother cell compartments occurs via a cascade of coherent feed-forward loops, and thereby can produce fast and ultrasensitive responses as a result of KinA induction. Unlike σ(F) activation, σ(E) activation in the mother cell compartment only occurs above the KinA threshold, resulting in completion of sporulation. Thus, ultrasensitive σ(E) activation explains the KinA threshold for sporulation induction. We therefore infer that under uncertain conditions, cells initiate sporulation but postpone making the sporulation decision to average stochastic fluctuations and to achieve a robust population response.

Mathematical model of a gene regulatory network reconciles effects of genetic perturbations on hematopoietic stem cell emergence.

Interlinked gene regulatory networks (GRNs) are vital for the spatial and temporal control of gene expression during development. The hematopoietic transcription factors (TFs) Scl, Gata2 and Fli1 form one such densely connected GRN which acts as a master regulator of embryonic hematopoiesis. This triad has been shown to direct the specification of the hemogenic endothelium and emergence of hematopoietic stem cells (HSCs) in response to Notch1 and Bmp4-Smad signaling. Here we employ previously published data to construct a mathematical model of this GRN network and use this model to systematically investigate the network dynamical properties. Our model uses a statistical-thermodynamic framework to describe the combinatorial regulation of gene expression and reconciles, mechanistically, several previously published but unexplained results from different genetic perturbation experiments. In particular, our results demonstrate how the interactions of Runx1, an essential hematopoietic TF, with components of the Bmp4 signaling pathway allow it to affect triad activation and acts as a key regulator of HSC emergence. We also explain why heterozygous deletion of this essential TF, Runx1, speeds up the network dynamics leading to accelerated HSC emergence. Taken together our results demonstrate that the triad, a master-level controller of definitive hematopoiesis, is an irreversible bistable switch whose dynamical properties are modulated by Runx1 and components of the Bmp4 signaling pathway.

Triggering sporulation in Bacillus subtilis with artificial two-component systems reveals the importance of proper Spo0A activation dynamics.

Sporulation initiation in Bacillus subtilis is controlled by the phosphorylated form of the master regulator Spo0A which controls transcription of a multitude of sporulation genes. In this study, we investigated the importance of temporal dynamics of phosphorylated Spo0A (Spo0A∼P) accumulation by rewiring the network controlling its phosphorylation. We showed that simultaneous induction of KinC, a kinase that can directly phosphorylate Spo0A, and Spo0A itself from separately controlled inducible promoters can efficiently trigger sporulation even under nutrient rich conditions. However, the sporulation efficiency in this artificial two-component system was significantly impaired when KinC and/or Spo0A induction was too high. Using mathematical modelling, we showed that gradual accumulation of Spo0A∼P is essential for the proper temporal order of the Spo0A regulon expression, and that reduction in sporulation efficiency results from the reversal of that order. These insights led us to identify premature repression of DivIVA as one possible explanation for the adverse effects of accelerated accumulation of Spo0A∼P on sporulation. Moreover, we found that positive feedback resulting from autoregulation of the native spo0A promoter leads to robust control of Spo0A∼P accumulation kinetics. Thus we propose that a major function of the conserved architecture of the sporulation network is controlling Spo0A activation dynamics.

Quantitative Model-Based Assessment of Multiple Sickle Cell Disease Therapeutic Approaches Alone and in Combination.

Three sickle cell disease (SCD) treatment strategies, stabilizing oxygenated hemoglobin (oxyHb), lowering 2,3-BPG, and inducing fetal hemoglobin (HbF) expression aim to prevent red blood cell (RBC) sickling by reducing tense-state sickle hemoglobin that contributes to polymer formation. Induction of 30% HbF is seen as the gold standard because 30% endogenous expression is associated with a lack of symptoms. However, the level of intervention required to achieve equivalent polymerization protection by the other strategies is uncertain, and there is little understanding of how these approaches could work in combination. We sought to develop an oxygen saturation model that could assess polymerization protection of all three approaches alone or in combination by extending the Monod-Wymann-Changeux model to include additional mechanisms. Applying the model to monotherapies suggests 51% sickle hemoglobin (HbS) occupancy with an oxyHb stabilizer or lowering RBC 2,3 BPG concentrations to 1.8 mM would produce comparable polymerization protection as 30% HbF. The model predictions are consistent with observed clinical response to the oxyHb stabilizer voxelotor and the 2,3-BPG reducer etavopivat. The model also suggests combination therapy will have added benefit in the case of dose limitations, as is the case for voxelotor, which the model predicts could be combined with 20% HbF or 2,3-BPG reduction to 3.75 mM to reach equivalent protection as 30% HbF. The proposed model represents a unified framework that is useful in supporting decisions in preclinical and early clinical development and capable of evolving with clinical experience to gain new and increasingly confident insights into treatment strategies for SCD.

To evaluate the analgesic effectiveness of bilateral erector spinae plane block versus thoracic epidural analgesia in open cardiac surgeries approached through midline sternotomy.

BACKGROUND: The efficacy of the erector spinae plane (ESP) block in mitigating postoperative pain has been shown for a range of thoracic and abdominal procedures. However, there is a paucity of literature investigating its impact on postoperative analgesia as well as its influence on weaning and subsequent recovery in comparison to thoracic epidural analgesia (TEA) in median sternotomy-based approach for open-cardiac surgeries and hence the study. METHODS: Irrespective of gender or age, 74 adult patients scheduled to undergo open cardiac surgery were enrolled and randomly allocated into two groups: the Group TEA (thoracic epidural block) and the Group ESP (bilateral Erector Spinae Plane block). The following variables were analysed prospectively and compared among the groups with regard to pain control, as determined by the VAS Scale both at rest (VASR) and during spirometry (VASS), time to extubation, quantity and frequency of rescue analgesia delivered, day of first ambulation, length of stay in the intensive care unit (ICU), and any adverse cardiac events (ACE), respiratory events (ARE), or other events, if pertinent. RESULTS: Clinical and demographic variables were similar in both groups. Both groups had overall good pain control, as determined by the VAS scale both at rest (VASR) and with spirometry (VASS) with Group ESP demonstrating superior pain regulation compared to Group TEA during the post-extubation period at 6, 9, and 12 h, respectively (P > 0.05). Although statistically insignificant, the postoperative mean rescue analgesic doses utilised in both groups were comparable, but there was a higher frequency requirement in Group TEA. The hemodynamic and respiratory profiles were comparable, except for a few arrhythmias in Group TEA. With comparable results, early recovery, fast-track extubation, and intensive care unit (ICU) stay were achieved. CONCLUSIONS: The ESP block has been found to have optimal analgesic effects during open cardiac surgery, resulting in a decreased need for additional analgesic doses and eliminating the possibility of a coagulation emergency. Consequently, it presents itself as a safer alternative to the potentially invasive thoracic epidural analgesia (TEA).

Quantitative systems pharmacology modeling provides insight into inter-mouse variability of Anti-CTLA4 response.

Clinical responses of immuno-oncology therapies are highly variable among patients. Similar response variability has been observed in syngeneic mouse models. Understanding of the variability in the mouse models may shed light on patient variability. Using a murine anti-CTLA4 antibody as a case study, we developed a quantitative systems pharmacology model to capture the molecular interactions of the antibody and relevant cellular interactions that lead to tumor cell killing. Nonlinear mixed effect modeling was incorporated to capture the inter-animal variability of tumor growth profiles in response to anti-CTLA4 treatment. The results suggested that intratumoral CD8+ T cell kinetics and tumor proliferation rate were the main drivers of the variability. In addition, simulations indicated that nonresponsive mice to anti-CTLA4 treatment could be converted to responders by increasing the number of intratumoral CD8+ T cells. The model provides a mechanistic starting point for translation of CTLA4 inhibitors from syngeneic mice to the clinic.

Modeling reveals bistability and low-pass filtering in the network module determining blood stem cell fate.

Combinatorial regulation of gene expression is ubiquitous in eukaryotes with multiple inputs converging on regulatory control elements. The dynamic properties of these elements determine the functionality of genetic networks regulating differentiation and development. Here we propose a method to quantitatively characterize the regulatory output of distant enhancers with a biophysical approach that recursively determines free energies of protein-protein and protein-DNA interactions from experimental analysis of transcriptional reporter libraries. We apply this method to model the Scl-Gata2-Fli1 triad-a network module important for cell fate specification of hematopoietic stem cells. We show that this triad module is inherently bistable with irreversible transitions in response to physiologically relevant signals such as Notch, Bmp4 and Gata1 and we use the model to predict the sensitivity of the network to mutations. We also show that the triad acts as a low-pass filter by switching between steady states only in response to signals that persist for longer than a minimum duration threshold. We have found that the auto-regulation loops connecting the slow-degrading Scl to Gata2 and Fli1 are crucial for this low-pass filtering property. Taken together our analysis not only reveals new insights into hematopoietic stem cell regulatory network functionality but also provides a novel and widely applicable strategy to incorporate experimental measurements into dynamical network models.

A Mechanistic Site-Of-Action Model: A Tool for Informing Right Target, Right Compound, And Right Dose for Therapeutic Antagonistic Antibody Programs.

Quantitative modeling is increasingly utilized in the drug discovery and development process, from the initial stages of target selection, through clinical studies. The modeling can provide guidance on three major questions-is this the right target, what are the right compound properties, and what is the right dose for moving the best possible candidate forward. In this manuscript, we present a site-of-action modeling framework which we apply to monoclonal antibodies against soluble targets. We give a comprehensive overview of how we construct the model and how we parametrize it and include several examples of how to apply this framework for answering the questions postulated above. The utilities and limitations of this approach are discussed.

One size does not fit all: navigating the multi-dimensional space to optimize T-cell engaging protein therapeutics.

T-cell engaging biologics is a class of novel and promising immune-oncology compounds that leverage the immune system to eradicate cancer. Here, we compared and contrasted a bispecific diabody-Fc format, which displays a relatively short antigen-binding arm distance, with our bispecific IgG platform. By generating diverse panels of antigen-expressing cells where B cell maturation antigen is either tethered to the cell membrane or located to the juxtamembrane region and masked by elongated structural spacer units, we presented a systematic approach to investigate the role of antigen epitope location and molecular formats in immunological synapse formation and cytotoxicity. We demonstrated that diabody-Fc is more potent for antigen epitopes located in the membrane distal region, while bispecific IgG is more efficient for membrane-proximal epitopes. Additionally, we explored other parameters, including receptor density, antigen-binding affinity, and kinetics. Our results show that molecular format and antigen epitope location, which jointly determine the intermembrane distance between target cells and T cells, allow decoupling of cytotoxicity and cytokine release, while antigen-binding affinities appear to be positively correlated with both readouts. Our work offers new insight that could potentially lead to a wider therapeutic window for T-cell engaging biologics in general.

Preclinical Evaluation of 89Zr-Df-IAB22M2C PET as an Imaging Biomarker for the Development of the GUCY2C-CD3 Bispecific PF-07062119 as a T Cell Engaging Therapy.

PURPOSE: A sensitive and specific imaging biomarker to monitor immune activation and quantify pharmacodynamic responses would be useful for development of immunomodulating anti-cancer agents. PF-07062119 is a T cell engaging bispecific antibody that binds to CD3 and guanylyl cyclase C, a protein that is over-expressed by colorectal cancers. Here, we used 89Zr-Df-IAB22M2C (89Zr-Df-Crefmirlimab), a human CD8-specific minibody to monitor CD8+ T cell infiltration into tumors by positron emission tomography. We investigated the ability of 89Zr-Df-IAB22M2C to track anti-tumor activity induced by PF-07062119 in a human CRC adoptive transfer mouse model (with injected activated/expanded human T cells), as well as the correlation of tumor radiotracer uptake with CD8+ immunohistochemical staining. PROCEDURES: NOD SCID gamma mice bearing human CRC LS1034 tumors were treated with four different doses of PF-07062119, or a non-targeted CD3 BsAb control, and imaged with 89Zr-Df-IAB22M2C PET at days 4 and 9. Following PET/CT imaging, mice were euthanized and dissected for ex vivo distribution analysis of 89Zr-Df-IAB22M2C in tissues on days 4 and 9, with additional data collected on day 6 (supplementary). Data were analyzed and reported as standard uptake value and %ID/g for in vivo imaging and ex vivo tissue distribution. In addition, tumor tissues were evaluated by immunohistochemistry for CD8+ T cells. RESULTS: The results demonstrated substantial mean uptake of 89Zr-Df-IAB22M2C (%ID/g) in PF-07062119-treated tumors, with significant increases in comparison to non-targeted BsAb-treated controls, as well as PF-07062119 dose-dependent responses over time of treatment. A moderate correlation was observed between tumor tissue radioactivity uptake and CD8+ cell density, demonstrating the value of the imaging agent for non-invasive assessment of intra-tumoral CD8+ T cells and the mechanism of action for PF-07062119. CONCLUSION: Immune-imaging technologies for quantitative cellular measures would be a valuable biomarker in immunotherapeutic clinical development. We demonstrated a qualification of 89Zr-IAB22M2C PET to evaluate PD responses (mice) to a novel immunotherapeutic.

Discovery and optimization of a novel anti-GUCY2c x CD3 bispecific antibody for the treatment of solid tumors.

We report here the discovery and optimization of a novel T cell retargeting anti-GUCY2C x anti-CD3ε bispecific antibody for the treatment of solid tumors. Using a combination of hybridoma, phage display and rational design protein engineering, we have developed a fully humanized and manufacturable CD3 bispecific antibody that demonstrates favorable pharmacokinetic properties and potent in vivo efficacy. Anti-GUCY2C and anti-CD3ε antibodies derived from mouse hybridomas were first humanized into well-behaved human variable region frameworks with full retention of binding and T-cell mediated cytotoxic activity. To address potential manufacturability concerns, multiple approaches were taken in parallel to optimize and de-risk the two antibody variable regions. These approaches included structure-guided rational mutagenesis and phage display-based optimization, focusing on improving stability, reducing polyreactivity and self-association potential, removing chemical liabilities and proteolytic cleavage sites, and de-risking immunogenicity. Employing rapid library construction methods as well as automated phage display and high-throughput protein production workflows enabled efficient generation of an optimized bispecific antibody with desirable manufacturability properties, high stability, and low nonspecific binding. Proteolytic cleavage and deamidation in complementarity-determining regions were also successfully addressed. Collectively, these improvements translated to a molecule with potent single-agent in vivo efficacy in a tumor cell line adoptive transfer model and a cynomolgus monkey pharmacokinetic profile (half-life>4.5 days) suitable for clinical development. Clinical evaluation of PF-07062119 is ongoing.

PF-07059013: A Noncovalent Modulator of Hemoglobin for Treatment of Sickle Cell Disease.

Sickle cell disease (SCD) is a genetic disorder caused by a single point mutation (ß6 Glu → Val) on the ß-chain of adult hemoglobin (HbA) that results in sickled hemoglobin (HbS). In the deoxygenated state, polymerization of HbS leads to sickling of red blood cells (RBC). Several downstream consequences of polymerization and RBC sickling include vaso-occlusion, hemolytic anemia, and stroke. We report the design of a noncovalent modulator of HbS, clinical candidate PF-07059013 (23). The seminal hit molecule was discovered by virtual screening and confirmed through a series of biochemical and biophysical studies. After a significant optimization effort, we arrived at 23, a compound that specifically binds to Hb with nanomolar affinity and displays strong partitioning into RBCs. In a 2-week multiple dose study using Townes SCD mice, 23 showed a 37.8% (±9.0%) reduction in sickling compared to vehicle treated mice. 23 (PF-07059013) has advanced to phase 1 clinical trials.

A Novel GUCY2C-CD3 T-Cell Engaging Bispecific Construct (PF-07062119) for the Treatment of Gastrointestinal Cancers.

PURPOSE: Gastrointestinal cancers remain areas of high unmet need despite advances in targeted and immunotherapies. Here, we demonstrate potent, tumor-selective efficacy with PF-07062119, a T-cell engaging CD3 bispecific targeting tumors expressing Guanylyl Cyclase C (GUCY2C), which is expressed widely across colorectal cancer and other gastrointestinal malignancies. In addition, to address immune evasion mechanisms, we explore combinations with immune checkpoint blockade agents and with antiangiogenesis therapy. EXPERIMENTAL DESIGN: PF-07062119 activity was evaluated in vitro in multiple tumor cell lines, and in vivo in established subcutaneous and orthotopic human colorectal cancer xenograft tumors with adoptive transfer of human T cells. Efficacy was also evaluated in mouse syngeneic tumors using human CD3ε transgenic mice. IHC and mass cytometry were performed to demonstrate drug biodistribution, recruitment of activated T cells, and to identify markers of immune evasion. Combination studies were performed with anti-PD-1/PD-L1 and anti-VEGF antibodies. Toxicity and pharmacokinetic studies were done in cynomolgus macaque. RESULTS: We demonstrate that GUCY2C-positive tumors can be targeted with an anti-GUCY2C/anti-CD3ε bispecific, with selective drug biodistribution to tumors. PF-07062119 showed potent T-cell-mediated in vitro activity and in vivo efficacy in multiple colorectal cancer human xenograft tumor models, including KRAS- and BRAF-mutant tumors, as well as in the immunocompetent mouse syngeneic tumor model. PF-07062119 activity was further enhanced when combined with anti-PD-1/PD-L1 treatment or in combination with antiangiogenic therapy. Toxicity studies in cynomolgus indicated a monitorable and manageable toxicity profile. CONCLUSIONS: These data highlight the potential for PF-07062119 to demonstrate efficacy and improve patient outcomes in colorectal cancer and other gastrointestinal malignancies.

Importance of Feedback and Feedforward Loops to Adaptive Immune Response Modeling.

The human adaptive immune system is a very complex network of different types of cells, cytokines, and signaling molecules. This complex network makes it difficult to understand the system level regulations. To properly explain the immune system, it is necessary to explicitly investigate the presence of different feedback and feedforward loops (FFLs) and their crosstalks. Considering that these loops increase the complexity of the system, the mathematical modeling has been proved to be an important tool to explain such complex biological systems. This review focuses on these regulatory loops and discusses their importance on systems modeling of the immune system.

Functional requirements of cellular differentiation: lessons from Bacillus subtilis.

Successful execution of differentiation programs requires cells to assess multitudes of internal and external cues and respond with appropriate gene expression programs. Here, we review how Bacillus subtilis sporulation network deals with these tasks focusing on the lessons generalizable to other systems. With feedforward loops controlling both production and activation of downstream transcriptional regulators, cells achieve ultrasensitive threshold-like responses. The arrangement of sporulation network genes on the chromosome and transcriptional feedback loops allow coordination of sporulation decision with DNA-replication. Furthermore, to assess the starvation conditions without sensing specific metabolites, cells respond to changes in their growth rates with increased activity of sporulation master regulator. These design features of the sporulation network enable cells to robustly decide between vegetative growth and sporulation.

Bistable responses in bacterial genetic networks: designs and dynamical consequences.

A key property of living cells is their ability to react to stimuli with specific biochemical responses. These responses can be understood through the dynamics of underlying biochemical and genetic networks. Evolutionary design principles have been well studied in networks that display graded responses, with a continuous relationship between input signal and system output. Alternatively, biochemical networks can exhibit bistable responses so that over a range of signals the network possesses two stable steady states. In this review, we discuss several conceptual examples illustrating network designs that can result in a bistable response of the biochemical network. Next, we examine manifestations of these designs in bacterial master-regulatory genetic circuits. In particular, we discuss mechanisms and dynamic consequences of bistability in three circuits: two-component systems, sigma-factor networks, and a multistep phosphorelay. Analyzing these examples allows us to expand our knowledge of evolutionary design principles networks with bistable responses.