Exactly Solvable Model for a Deconfined Quantum Critical Point in 1D.

We construct an exactly solvable lattice model for a deconfined quantum critical point (DQCP) in (1+1) dimensions. This DQCP occurs in an unusual setting, namely, at the edge of a (2+1) dimensional bosonic symmetry protected topological (SPT) phase with Z_{2}×Z_{2} symmetry. The DQCP describes a transition between two gapped edges that break different Z_{2} subgroups of the full Z_{2}×Z_{2} symmetry. Our construction is based on an exact mapping between the SPT edge theory and a Z_{4} spin chain. This mapping reveals that DQCPs in this system are directly related to ordinary Z_{4} symmetry breaking critical points.

Volume-Law to Area-Law Entanglement Transition in a Nonunitary Periodic Gaussian Circuit.

We consider Gaussian quantum circuits that alternate unitary gates and postselected weak measurements, with spatial translation symmetry and time p eriodicity. We show analytically that such models can host different kinds of measurement-induced phase transitions detected by entanglement entropy, by mapping the unitary gates and weak measurements onto Möbius transformations. We demonstrate the existence of a log-law to area-law transition, as well as a volume-law to area-law transition at a finite measurement amplitude. For the latter, we compute the critical exponent ν for the Hartley, von Neumann and Rényi entropies exactly.

Temporal encoding of bacterial identity and traits in growth dynamics.

In biology, it is often critical to determine the identity of an organism and phenotypic traits of interest. Whole-genome sequencing can be useful for this but has limited power for trait prediction. However, we can take advantage of the inherent information content of phenotypes to bypass these limitations. We demonstrate, in clinical and environmental bacterial isolates, that growth dynamics in standardized conditions can differentiate between genotypes, even among strains from the same species. We find that for pairs of isolates, there is little correlation between genetic distance, according to phylogenetic analysis, and phenotypic distance, as determined by growth dynamics. This absence of correlation underscores the challenge in using genomics to infer phenotypes and vice versa. Bypassing this complexity, we show that growth dynamics alone can robustly predict antibiotic responses. These findings are a foundation for a method to identify traits not easily traced to a genetic mechanism.

Prophage Hunter: an integrative hunting tool for active prophages.

Identifying active prophages is critical for studying coevolution of phage and bacteria, investigating phage physiology and biochemistry, and engineering designer phages for diverse applications. We present Prophage Hunter, a tool aimed at hunting for active prophages from whole genome assembly of bacteria. Combining sequence similarity-based matching and genetic features-based machine learning classification, we developed a novel scoring system that exhibits higher accuracy than current tools in predicting active prophages on the validation datasets. The option of skipping similarity matching is also available so that there's higher chance for novel phages to be discovered. Prophage Hunter provides a one-stop web service to extract prophage genomes from bacterial genomes, evaluate the activity of the prophages, identify phylogenetically related phages, and annotate the function of phage proteins. Prophage Hunter is freely available at https://pro-hunter.bgi.com/.

Metabolic division of labor in microbial systems.

Metabolic pathways are often engineered in single microbial populations. However, the introduction of heterologous circuits into the host can create a substantial metabolic burden that limits the overall productivity of the system. This limitation could be overcome by metabolic division of labor (DOL), whereby distinct populations perform different steps in a metabolic pathway, reducing the burden each population will experience. While conceptually appealing, the conditions when DOL is advantageous have not been rigorously established. Here, we have analyzed 24 common architectures of metabolic pathways in which DOL can be implemented. Our analysis reveals general criteria defining the conditions that favor DOL, accounting for the burden or benefit of the pathway activity on the host populations as well as the transport and turnover of enzymes and intermediate metabolites. These criteria can help guide engineering of metabolic pathways and have implications for understanding evolution of natural microbial communities.

Nonlinear response to cancer nanotherapy due to macrophage interactions revealed by mathematical modeling and evaluated in a murine model via CRISPR-modulated macrophage polarization.

Tumor-associated macrophages (TAMs) have been shown to both aid and hinder tumor growth, with patient outcomes potentially hinging on the proportion of M1, pro-inflammatory/growth-inhibiting, to M2, growth-supporting, phenotypes. Strategies to stimulate tumor regression by promoting polarization to M1 are a novel approach that harnesses the immune system to enhance therapeutic outcomes, including chemotherapy. We recently found that nanotherapy with mesoporous particles loaded with albumin-bound paclitaxel (MSV-nab-PTX) promotes macrophage polarization towards M1 in breast cancer liver metastases (BCLM). However, it remains unclear to what extent tumor regression can be maximized based on modulation of the macrophage phenotype, especially for poorly perfused tumors such as BCLM. Here, for the first time, a CRISPR system is employed to permanently modulate macrophage polarization in a controlled in vitro setting. This enables the design of 3D co-culture experiments mimicking the BCLM hypovascularized environment with various ratios of polarized macrophages. We implement a mathematical framework to evaluate nanoparticle-mediated chemotherapy in conjunction with TAM polarization. The response is predicted to be not linearly dependent on the M1:M2 ratio. To investigate this phenomenon, the response is simulated via the model for a variety of M1:M2 ratios. The modeling indicates that polarization to an all-M1 population may be less effective than a combination of both M1 and M2. Experimental results with the CRISPR system confirm this model-driven hypothesis. Altogether, this study indicates that response to nanoparticle-mediated chemotherapy targeting poorly perfused tumors may benefit from a fine-tuned M1:M2 ratio that maintains both phenotypes in the tumor microenvironment during treatment.

Processing Oscillatory Signals by Incoherent Feedforward Loops.

From the timing of amoeba development to the maintenance of stem cell pluripotency, many biological signaling pathways exhibit the ability to differentiate between pulsatile and sustained signals in the regulation of downstream gene expression. While the networks underlying this signal decoding are diverse, many are built around a common motif, the incoherent feedforward loop (IFFL), where an input simultaneously activates an output and an inhibitor of the output. With appropriate parameters, this motif can exhibit temporal adaptation, where the system is desensitized to a sustained input. This property serves as the foundation for distinguishing input signals with varying temporal profiles. Here, we use quantitative modeling to examine another property of IFFLs-the ability to process oscillatory signals. Our results indicate that the system's ability to translate pulsatile dynamics is limited by two constraints. The kinetics of the IFFL components dictate the input range for which the network is able to decode pulsatile dynamics. In addition, a match between the network parameters and input signal characteristics is required for optimal "counting". We elucidate one potential mechanism by which information processing occurs in natural networks, and our work has implications in the design of synthetic gene circuits for this purpose.

Chlorfenapyr Crystal Polymorphism and Insecticidal Activity.

Four crystalline polymorphs of the proinsecticide chlorfenapyr [4-bromo-2-(4-chlorophenyl)-1-ethoxymethyl-5-trifluoromethyl-1H-pyrrole-3-carbonitrile] have been identified and characterized by polarized light optical microscopy, differential scanning calorimetry, Raman spectroscopy, X-ray diffraction, and electron diffraction. Three of the four structures were considered polytypic. Chlorfenapyr polymorphs show similar lethality against fruit flies (Drosophila melanogaster) and mosquitoes (Anopheles quadrimaculatus) with the least stable polymorph showing slightly higher lethality. Similar activities may be expected to be consistent with structural similarities. Knockdown kinetics, however, depend on an internal metabolic activating step, which further complicates polymorph-dependent bioavailability.

Applications of Model-Based Target Pharmacology Assessment in Defining Drug Design and DMPK Strategies: GSK Experiences.

Many critical decisions faced in early discovery require a thorough understanding of the dynamic behavior of pharmacological pathways following target engagement. From fundamental decisions on the optimal target to pursue and the ultimate drug product profile (combination of modality, potency, and compound properties) expected to elicit the desired clinical outcome to tactical program decisions such as what chemical series to pursue, what chemical properties require optimization, and what compounds to synthesize and progress, all demand detailed consideration of pharmacodynamics. Model-based target pharmacology assessment (mTPA) is a computational approach centered around large-scale virtual exploration of pharmacokinetic and pharmacodynamic models built early in discovery to guide these decisions. The present work summarizes several examples (use cases) from programs at GlaxoSmithKline that demonstrate the utility of mTPA throughout the drug discovery lifecycle.

Quantitative immunopeptidomics reveals a tumor stroma-specific target for T cell therapy.

T cell receptor (TCR)-based immunotherapy has emerged as a promising therapeutic approach for the treatment of patients with solid cancers. Identifying peptide-human leukocyte antigen (pHLA) complexes highly presented on tumors and rarely expressed on healthy tissue in combination with high-affinity TCRs that when introduced into T cells can redirect T cells to eliminate tumor but not healthy tissue is a key requirement for safe and efficacious TCR-based therapies. To discover promising shared tumor antigens that could be targeted via TCR-based adoptive T cell therapy, we employed population-scale immunopeptidomics using quantitative mass spectrometry across ~1500 tumor and normal tissue samples. We identified an HLA-A*02:01-restricted pan-cancer epitope within the collagen type VI α-3 (COL6A3) gene that is highly presented on tumor stroma across multiple solid cancers due to a tumor-specific alternative splicing event that rarely occurs outside the tumor microenvironment. T cells expressing natural COL6A3-specific TCRs demonstrated only modest activity against cells presenting high copy numbers of COL6A3 pHLAs. One of these TCRs was affinity-enhanced, enabling transduced T cells to specifically eliminate tumors in vivo that expressed similar copy numbers of pHLAs as primary tumor specimens. The enhanced TCR variants exhibited a favorable safety profile with no detectable off-target reactivity, paving the way to initiate clinical trials using COL6A3-specific TCRs to target an array of solid tumors.

Caspase Activation and Caspase-Mediated Cleavage of APP Is Associated with Amyloid ß-Protein-Induced Synapse Loss in Alzheimer's Disease.

Amyloid ß-protein (Aß) toxicity is hypothesized to play a seminal role in Alzheimer's disease (AD) pathogenesis. However, it remains unclear how Aß causes synaptic dysfunction and synapse loss. We hypothesize that one mechanism of Aß-induced synaptic injury is related to the cleavage of amyloid ß precursor protein (APP) at position D664 by caspases that release the putatively cytotoxic C31 peptide. In organotypic slice cultures derived from mice with a knock-in mutation in the APP gene (APP D664A) to inhibit caspase cleavage, Aß-induced synaptic injury is markedly reduced in two models of Aß toxicity. Loss of dendritic spines is also attenuated in mice treated with caspase inhibitors. Importantly, the time-dependent dendritic spine loss is correlated with localized activation of caspase-3 but is absent in APP D664A cultures. We propose that the APP cytosolic domain plays an essential role in Aß-induced synaptic damage in the injury pathway mediated by localized caspase activation.

Massive computational acceleration by using neural networks to emulate mechanism-based biological models.

For many biological applications, exploration of the massive parametric space of a mechanism-based model can impose a prohibitive computational demand. To overcome this limitation, we present a framework to improve computational efficiency by orders of magnitude. The key concept is to train a neural network using a limited number of simulations generated by a mechanistic model. This number is small enough such that the simulations can be completed in a short time frame but large enough to enable reliable training. The trained neural network can then be used to explore a much larger parametric space. We demonstrate this notion by training neural networks to predict pattern formation and stochastic gene expression. We further demonstrate that using an ensemble of neural networks enables the self-contained evaluation of the quality of each prediction. Our work can be a platform for fast parametric space screening of biological models with user defined objectives.

Cyclin D/CDK4/6 activity controls G1 length in mammalian cells.

The length of the G1 phase in the cell cycle shows significant variability in different cell types and tissue types. To gain insights into the control of G1 length, we generated an E2F activity reporter that captures free E2F activity after dissociation from Rb sequestration and followed its kinetics of activation at the single-cell level, in real time. Our results demonstrate that its activity is precisely coordinated with S phase progression. Quantitative analysis indicates that there is a pre-S phase delay between E2F transcriptional dynamic and activity dynamics. This delay is variable among different cell types and is strongly modulated by the cyclin D/CDK4/6 complex activity through Rb phosphorylation. Our findings suggest that the main function of this complex is to regulate the appropriate timing of G1 length.

Expression level is a key determinant of E2F1-mediated cell fate.

The Rb/E2F network has a critical role in regulating cell cycle progression and cell fate decisions. It is dysfunctional in virtually all human cancers, because of genetic lesions that cause overexpression of activators, inactivation of repressors, or both. Paradoxically, the downstream target of this network, E2F1, is rarely strongly overexpressed in cancer. E2F1 can induce both proliferation and apoptosis but the factors governing these critical cell fate decisions remain unclear. Previous studies have focused on qualitative mechanisms such as differential cofactors, posttranslational modification or state of other signaling pathways as modifiers of the cell fate decisions downstream of E2F1 activation. In contrast, the importance of the expression levels of E2F1 itself in dictating the downstream phenotypes has not been rigorously studied, partly due to the limited resolution of traditional population-level measurements. Here, through single-cell quantitative analysis, we demonstrate that E2F1 expression levels have a critical role in determining the fate of individual cells. Low levels of exogenous E2F1 promote proliferation, moderate levels induce G1, G2 and mitotic cell cycle arrest, and very high levels promote apoptosis. These multiple anti-proliferative mechanisms result in a strong selection pressure leading to rapid elimination of E2F1-overexpressing cells from the population. RNA-sequencing and RT-PCR revealed that low levels of E2F1 are sufficient to induce numerous cell cycle-promoting genes, intermediate levels induce growth arrest genes (i.e., p18, p19 and p27), whereas higher levels are necessary to induce key apoptotic E2F1 targets APAF1, PUMA, HRK and BIM. Finally, treatment of a lung cancer cell line with a proteasome inhibitor, MLN2238, resulted in an E2F1-dependent mitotic arrest and apoptosis, confirming the role of endogenous E2F1 levels in these phenotypes. The strong anti-proliferative activity of moderately overexpressed E2F1 in multiple cancer types suggests that targeting E2F1 for upregulation may represent an attractive therapeutic strategy in cancer.

Addressing biological uncertainties in engineering gene circuits.

Synthetic biology has grown tremendously over the past fifteen years. It represents a new strategy to develop biological understanding and holds great promise for diverse practical applications. Engineering of a gene circuit typically involves computational design of the circuit, selection of circuit components, and test and optimization of circuit functions. A fundamental challenge in this process is the predictable control of circuit function due to multiple layers of biological uncertainties. These uncertainties can arise from different sources. We categorize these uncertainties into incomplete quantification of parts, interactions between heterologous components and the host, or stochastic dynamics of chemical reactions and outline potential design strategies to minimize or exploit them.