Advances in preclinical TCR characterization: leveraging cell avidity to identify functional TCRs.

T-cell therapy has emerged as an effective approach for treating viral infections and cancers. However, a significant challenge is the selection of T-cell receptors (TCRs) that exhibit the desired functionality. Conventionally in vitro techniques, such as peptide sensitivity measurements and cytotoxicity assays, provide valuable insights into TCR potency but are labor-intensive. In contrast, measuring ligand binding properties (z-Movi technology) could provide an accelerated processing while showing robust correlations with T-cell functions. In this study, we assessed whether cell avidity can predict functionality also in the context of TCR-engineered T cells. To this end, we developed a flexible system for TCR re-expression by generating a Jurkat-derived T cell clone lacking TCR and CD3 expression through CRISPR-Cas9-mediated TRBC knockout. The knockin of a transgenic TCR into the TRAC locus restored TCR/CD3 expression, allowing for CD3-based purification of TCR-engineered T cells. Subsequently, we characterized these engineered cell lines by functional readouts, and assessment of binding properties through the z-Movi technology. Our findings revealed a strong correlation between the cell avidities and functional sensitivities of Jurkat TCR-T cells. Altogether, by integrating cell avidity measurements with our versatile T cell engineering platform, we established an accelerated system for enhancing the in vitro selection of clinically relevant TCRs.

Active T1 transitions in cellular networks.

In amorphous solids as in tissues, neighbor exchanges can relax local stresses and allow the material to flow. In this paper, we use an anisotropic vertex model to study T1 rearrangements in polygonal cellular networks. We consider two different physical realizations of the active anisotropic stresses: (i) anisotropic bond tension and (ii) anisotropic cell stress. Interestingly, the two types of active stress lead to patterns of relative orientation of T1 transitions and cell elongation that are different. Our work suggests that these two realizations of anisotropic active stresses can be observed in vivo. We describe and explain these results through the lens of a continuum description of the tissue as an anisotropic active material. We furthermore discuss the energetics of the dynamic tissue and express the energy balance in terms of internal elastic energy, mechanical work, chemical work and heat. This allows us to define active T1 transitions that can perform mechanical work while consuming chemical energy.

eGFRD in all dimensions.

Biochemical reactions often occur at low copy numbers but at once in crowded and diverse environments. Space and stochasticity therefore play an essential role in biochemical networks. Spatial-stochastic simulations have become a prominent tool for understanding how stochasticity at the microscopic level influences the macroscopic behavior of such systems. While particle-based models guarantee the level of detail necessary to accurately describe the microscopic dynamics at very low copy numbers, the algorithms used to simulate them typically imply trade-offs between computational efficiency and biochemical accuracy. eGFRD (enhanced Green's Function Reaction Dynamics) is an exact algorithm that evades such trade-offs by partitioning the N-particle system into M ≤ N analytically tractable one- and two-particle systems; the analytical solutions (Green's functions) then are used to implement an event-driven particle-based scheme that allows particles to make large jumps in time and space while retaining access to their state variables at arbitrary simulation times. Here we present "eGFRD2," a new eGFRD version that implements the principle of eGFRD in all dimensions, thus enabling efficient particle-based simulation of biochemical reaction-diffusion processes in the 3D cytoplasm, on 2D planes representing membranes, and on 1D elongated cylinders representative of, e.g., cytoskeletal tracks or DNA; in 1D, it also incorporates convective motion used to model active transport. We find that, for low particle densities, eGFRD2 is up to 6 orders of magnitude faster than conventional Brownian dynamics. We exemplify the capabilities of eGFRD2 by simulating an idealized model of Pom1 gradient formation, which involves 3D diffusion, active transport on microtubules, and autophosphorylation on the membrane, confirming recent experimental and theoretical results on this system to hold under genuinely stochastic conditions.

Discrete gene replication events drive coupling between the cell cycle and circadian clocks.

Many organisms possess both a cell cycle to control DNA replication and a circadian clock to anticipate changes between day and night. In some cases, these two rhythmic systems are known to be coupled by specific, cross-regulatory interactions. Here, we use mathematical modeling to show that, additionally, the cell cycle generically influences circadian clocks in a nonspecific fashion: The regular, discrete jumps in gene-copy number arising from DNA replication during the cell cycle cause a periodic driving of the circadian clock, which can dramatically alter its behavior and impair its function. A clock built on negative transcriptional feedback either phase-locks to the cell cycle, so that the clock period tracks the cell division time, or exhibits erratic behavior. We argue that the cyanobacterium Synechococcus elongatus has evolved two features that protect its clock from such disturbances, both of which are needed to fully insulate it from the cell cycle and give it its observed robustness: a phosphorylation-based protein modification oscillator, together with its accompanying push-pull read-out circuit that responds primarily to the ratios of different phosphoform concentrations, makes the clock less susceptible to perturbations in protein synthesis; the presence of multiple, asynchronously replicating copies of the same chromosome diminishes the effect of replicating any single copy of a gene.

A thermodynamically consistent model of the post-translational Kai circadian clock.

The principal pacemaker of the circadian clock of the cyanobacterium S. elongatus is a protein phosphorylation cycle consisting of three proteins, KaiA, KaiB and KaiC. KaiC forms a homohexamer, with each monomer consisting of two domains, CI and CII. Both domains can bind and hydrolyze ATP, but only the CII domain can be phosphorylated, at two residues, in a well-defined sequence. While this system has been studied extensively, how the clock is driven thermodynamically has remained elusive. Inspired by recent experimental observations and building on ideas from previous mathematical models, we present a new, thermodynamically consistent, statistical-mechanical model of the clock. At its heart are two main ideas: i) ATP hydrolysis in the CI domain provides the thermodynamic driving force for the clock, switching KaiC between an active conformational state in which its phosphorylation level tends to rise and an inactive one in which it tends to fall; ii) phosphorylation of the CII domain provides the timer for the hydrolysis in the CI domain. The model also naturally explains how KaiA, by acting as a nucleotide exchange factor, can stimulate phosphorylation of KaiC, and how the differential affinity of KaiA for the different KaiC phosphoforms generates the characteristic temporal order of KaiC phosphorylation. As the phosphorylation level in the CII domain rises, the release of ADP from CI slows down, making the inactive conformational state of KaiC more stable. In the inactive state, KaiC binds KaiB, which not only stabilizes this state further, but also leads to the sequestration of KaiA, and hence to KaiC dephosphorylation. Using a dedicated kinetic Monte Carlo algorithm, which makes it possible to efficiently simulate this system consisting of more than a billion reactions, we show that the model can describe a wealth of experimental data.

Period Robustness and Entrainability of the Kai System to Changing Nucleotide Concentrations.

Circadian clocks must be able to entrain to time-varying signals to keep their oscillations in phase with the day-night rhythm. On the other hand, they must also exhibit input compensation: their period must remain approximately one day in different constant environments. The posttranslational oscillator of the Kai system can be entrained by transient or oscillatory changes in the ATP fraction, yet is insensitive to constant changes in this fraction. We study in three different models of this system how these two seemingly conflicting criteria are met. We find that one of these (our recently published Paijmans model) exhibits the best tradeoff between input compensation and entrainability: on the footing of equal phase-response curves, it exhibits the strongest input compensation. Performing stochastic simulations at the level of individual hexamers allows us to identify a new, to our knowledge, mechanism, which is employed by the Paijmans model to achieve input compensation: at lower ATP fraction, the individual hexamers make a shorter cycle in the phosphorylation state space, which compensates for the slower pace at which they traverse the cycle.

The Berg-Purcell limit revisited.

Biological systems often have to measure extremely low concentrations of chemicals with high precision. When dealing with such small numbers of molecules, the inevitable randomness of physical transport processes and binding reactions will limit the precision with which measurements can be made. An important question is what the lower bound on the noise would be in such measurements. Using the theory of diffusion-influenced reactions, we derive an analytical expression for the precision of concentration estimates that are obtained by monitoring the state of a receptor to which a diffusing ligand can bind. The variance in the estimate consists of two terms, one resulting from the intrinsic binding kinetics and the other from the diffusive arrival of ligand at the receptor. The latter term is identical to the fundamental limit derived by Berg and Purcell (Biophys. J., 1977), but disagrees with a more recent expression by Bialek and Setayeshgar. Comparing the theoretical predictions against results from particle-based simulations confirms the accuracy of the resulting expression and reaffirms the fundamental limit established by Berg and Purcell.

Active shape programming drives Drosophila wing disc eversion.

How complex 3D tissue shape emerges during animal development remains an important open question in biology and biophysics. Here, we discover a mechanism for 3D epithelial shape change based on active, in-plane cellular events that is analogous to inanimate "shape programmable" materials, which undergo blueprinted 3D shape transformations from in-plane gradients of spontaneous strains. We study eversion of the Drosophila wing disc pouch, when the epithelium transforms from a dome into a curved fold, quantifying 3D tissue shape changes and mapping spatial patterns of cellular behaviors on the evolving geometry using cellular topology. Using a physical model inspired by shape programming, we find that active cell rearrangements are the major contributor to pouch eversion and validate this conclusion using a knockdown of MyoVI, which reduces rearrangements and disrupts morphogenesis. This work shows that shape programming is a mechanism for animal tissue morphogenesis and suggests that patterns in nature could present design strategies for shape-programmable materials.

Nonlinear rheology of cellular networks.

Morphogenesis depends crucially on the complex rheological properties of cell tissues and on their ability to maintain mechanical integrity while rearranging at long times. In this paper, we study the rheology of polygonal cellular networks described by a vertex model in the presence of fluctuations. We use a triangulation method to decompose shear into cell shape changes and cell rearrangements. Considering the steady-state stress under constant shear, we observe nonlinear shear-thinning behavior at all magnitudes of the fluctuations, and an even stronger nonlinear regime at lower values of the fluctuations. We successfully capture this nonlinear rheology by a mean-field model that describes the tissue in terms of cell elongation and cell rearrangements. We furthermore introduce anisotropic active stresses in the vertex model and analyze their effect on rheology. We include this anisotropy in the mean-field model and show that it recapitulates the behavior observed in the simulations. Our work clarifies how tissue rheology is related to stochastic cell rearrangements and provides a simple biophysical model to describe biological tissues. Further, it highlights the importance of nonlinearities when discussing tissue mechanics.

Epithelial Viscoelasticity Is Regulated by Mechanosensitive E-cadherin Turnover.

Studying how epithelia respond to mechanical stresses is key to understanding tissue shape changes during morphogenesis. Here, we study the viscoelastic properties of the Drosophila wing epithelium during pupal morphogenesis by quantifying mechanical stress and cell shape as a function of time. We find a delay of 8 h between maximal tissue stress and maximal cell elongation, indicating a viscoelastic deformation of the tissue. We show that this viscoelastic behavior emerges from the mechanosensitivity of endocytic E-cadherin turnover. The increase in E-cadherin turnover in response to stress is mediated by mechanosensitive relocalization of the E-cadherin binding protein p120-catenin (p120) from cell junctions to cytoplasm. Mechanosensitivity of E-cadherin turnover is lost in p120 mutant wings, where E-cadherin turnover is constitutively high. In this mutant, the relationship between mechanical stress and stress-dependent cell dynamics is altered. Cells in p120 mutant deform and undergo cell rearrangements oriented along the stress axis more rapidly in response to mechanical stress. These changes imply a lower viscosity of wing epithelium. Taken together, our findings reveal that p120-dependent mechanosensitive E-cadherin turnover regulates viscoelastic behavior of epithelial tissues.

Robustness of synthetic oscillators in growing and dividing cells.

Synthetic biology sets out to implement new functions in cells, and to develop a deeper understanding of biological design principles. Elowitz and Leibler [Nature (London) 403, 335 (2000)NATUAS0028-083610.1038/35002125] showed that by rational design of the reaction network, and using existing biological components, they could create a network that exhibits periodic gene expression, dubbed the repressilator. More recently, Stricker et al. [Nature (London) 456, 516 (2008)NATUAS0028-083610.1038/nature07389] presented another synthetic oscillator, called the dual-feedback oscillator, which is more stable. Detailed studies have been carried out to determine how the stability of these oscillators is affected by the intrinsic noise of the interactions between the components and the stochastic expression of their genes. However, as all biological oscillators reside in growing and dividing cells, an important question is how these oscillators are perturbed by the cell cycle. In previous work we showed that the periodic doubling of the gene copy numbers due to DNA replication can couple not only natural, circadian oscillators to the cell cycle [Paijmans et al., Proc. Natl. Acad. Sci. (USA) 113, 4063 (2016)PNASA60027-842410.1073/pnas.1507291113], but also these synthetic oscillators. Here we expand this study. We find that the strength of the locking between oscillators depends not only on the positions of the genes on the chromosome, but also on the noise in the timing of gene replication: noise tends to weaken the coupling. Yet, even in the limit of high levels of noise in the replication times of the genes, both synthetic oscillators show clear signatures of locking to the cell cycle. This work enhances our understanding of the design of robust biological oscillators inside growing and diving cells.

Lower bound on the precision of transcriptional regulation and why facilitated diffusion can reduce noise in gene expression.

The diffusive arrival of transcription factors at the promoter sites on DNA sets a lower bound on how accurately a cell can regulate its protein levels. Using results from the literature on diffusion-influenced reactions, we derive an analytical expression for the lower bound on the precision of transcriptional regulation. In our theory, transcription factors can perform multiple rounds of one-dimensional (1D) diffusion along the DNA and 3D diffusion in the cytoplasm before binding to the promoter. Comparing our expression for the lower bound on the precision against results from Green's function reaction dynamics simulations shows that the theory is highly accurate under biologically relevant conditions. Our results demonstrate that, to an excellent approximation, the promoter switches between the transcription-factor bound and unbound state in a Markovian fashion. This remains true even in the presence of sliding, i.e., with 1D diffusion along the DNA. This has two important implications: (1) Minimizing the noise in the promoter state is equivalent to minimizing the search time of transcription factors for their promoters; (2) the complicated dynamics of 3D diffusion in the cytoplasm and 1D diffusion along the DNA can be captured in a well-stirred model by renormalizing the promoter association and dissociation rates, making it possible to efficiently simulate the promoter dynamics using Gillespie simulations. Based on the recent experimental observation that sliding can speed up the promoter search by a factor of 4, our theory predicts that sliding can enhance the precision of transcriptional regulation by a factor of 2.