Gene regulatory logic for reading the Sonic Hedgehog signaling gradient in the vertebrate neural tube.

Secreted signals, known as morphogens, provide the positional information that organizes gene expression and cellular differentiation in many developing tissues. In the vertebrate neural tube, Sonic Hedgehog (Shh) acts as a morphogen to control the pattern of neuronal subtype specification. Using an in vivo reporter of Shh signaling, mouse genetics, and systems modeling, we show that a spatially and temporally changing gradient of Shh signaling is interpreted by the regulatory logic of a downstream transcriptional network. The design of the network, which links three transcription factors to Shh signaling, is responsible for differential spatial and temporal gene expression. In addition, the network renders cells insensitive to fluctuations in signaling and confers hysteresis--memory of the signal. Our findings reveal that morphogen interpretation is an emergent property of the architecture of a transcriptional network that provides robustness and reliability to tissue patterning.

'Neighbourhood watch' model: embryonic epiblast cells assess positional information in relation to their neighbours.

In many developing and regenerating systems, tissue pattern is established through gradients of informative morphogens, but we know little about how cells interpret these. Using experimental manipulation of early chick embryos, including misexpression of an inducer (VG1 or ACTIVIN) and an inhibitor (BMP4), we test two alternative models for their ability to explain how the site of primitive streak formation is positioned relative to the rest of the embryo. In one model, cells read morphogen concentrations cell-autonomously. In the other, cells sense changes in morphogen status relative to their neighbourhood. We find that only the latter model can account for the experimental results, including some counter-intuitive predictions. This mechanism (which we name the 'neighbourhood watch' model) illuminates the classic 'French Flag Problem' and how positional information is interpreted by a sheet of cells in a large developing system.

Amino acid sensor conserved from bacteria to humans.

SignificanceAmino acids are the building blocks of life and important signaling molecules. Despite their common structure, no universal mechanism for amino acid recognition by cellular receptors is currently known. We discovered a simple motif, which binds amino acids in various receptor proteins from all major life-forms. In humans, this motif is found in subunits of calcium channels that are implicated in pain and neurodevelopmental disorders. Our findings suggest that γ-aminobutyric acid-derived drugs bind to the same motif in human proteins that binds natural ligands in bacterial receptors, thus enabling future improvement of important drugs.

Effect of abrocitinib on skin biomarkers in patients with moderate-to-severe atopic dermatitis.

BACKGROUND: This is the first report on the effects of abrocitinib, a Janus kinase 1-selective inhibitor, on the expression of skin biomarkers in patients with moderate-to-severe atopic dermatitis (AD). METHODS: JADE MOA (NCT03915496) was a double-blind Phase 2a trial. Adults were randomly assigned 1:1:1 to receive monotherapy with once-daily abrocitinib 200 mg, abrocitinib 100 mg, or placebo for 12 weeks. The primary endpoint was change from baseline in markers of inflammation (matrix metalloproteinase [MMP]-12), epidermal hyperplasia (keratin-16 [KRT16]), T-helper 2 (Th2) immune response (C-C motif chemokine ligand [CCL]17, CCL18, and CCL26), and Th22 immune response (S100 calcium binding protein A8, A9, and A12 [S100A8, S100A9, and S100A12]) in skin through 12 weeks. RESULTS: A total of 46 patients received abrocitinib 200 mg (n = 14), abrocitinib 100 mg (n = 16), or placebo (n = 16). Abrocitinib improved AD clinical signs and reduced itch. Gene expression of MMP-12, KRT16, S100A8, S100A9, and S100A12 was significantly decreased from baseline with abrocitinib 200 mg (at Weeks 2, 4, and 12) and abrocitinib 100 mg (at Weeks 4 and 12) in a dose-dependent manner. Abrocitinib 200 mg resulted in significant decreases from baseline in CCL17 expression at Week 12 and CCL18 expression at Weeks 2, 4, and 12; no significant decreases were observed for CCL26. CONCLUSIONS: Alongside improvements in clinical signs and symptoms of AD, 12 weeks of abrocitinib treatment resulted in downregulation of genes associated with inflammation, epidermal hyperplasia, and Th2 and Th22 immune responses in the skin of patients with moderate-to-severe AD.

Structure modeling hints at a granular organization of the Golgi ribbon.

BACKGROUND: In vertebrate cells, the Golgi functional subunits, mini-stacks, are linked into a tri-dimensional network. How this "ribbon" architecture relates to Golgi functions remains unclear. Are all connections between mini-stacks equal? Is the local structure of the ribbon of functional importance? These are difficult questions to address, without a quantifiable readout of the output of ribbon-embedded mini-stacks. Endothelial cells produce secretory granules, the Weibel-Palade bodies (WPB), whose von Willebrand Factor (VWF) cargo is central to hemostasis. The Golgi apparatus controls WPB size at both mini-stack and ribbon levels. Mini-stack dimensions delimit the size of VWF "boluses" whilst the ribbon architecture allows their linear co-packaging, thereby generating WPBs of different lengths. This Golgi/WPB size relationship suits mathematical analysis. RESULTS: WPB lengths were quantized as multiples of the bolus size and mathematical modeling simulated the effects of different Golgi ribbon organizations on WPB size, to be compared with the ground truth of experimental data. An initial simple model, with the Golgi as a single long ribbon composed of linearly interlinked mini-stacks, was refined to a collection of mini-ribbons and then to a mixture of mini-stack dimers plus long ribbon segments. Complementing these models with cell culture experiments led to novel findings. Firstly, one-bolus sized WPBs are secreted faster than larger secretory granules. Secondly, microtubule depolymerization unlinks the Golgi into equal proportions of mini-stack monomers and dimers. Kinetics of binding/unbinding of mini-stack monomers underpinning the presence of stable dimers was then simulated. Assuming that stable mini-stack dimers and monomers persist within the ribbon resulted in a final model that predicts a "breathing" arrangement of the Golgi, where monomer and dimer mini-stacks within longer structures undergo continuous linking/unlinking, consistent with experimentally observed WPB size distributions. CONCLUSIONS: Hypothetical Golgi organizations were validated against a quantifiable secretory output. The best-fitting Golgi model, accounting for stable mini-stack dimers, is consistent with a highly dynamic ribbon structure, capable of rapid rearrangement. Our modeling exercise therefore predicts that at the fine-grained level the Golgi ribbon is more complex than generally thought. Future experiments will confirm whether such a ribbon organization is endothelial-specific or a general feature of vertebrate cells.

Neuronal differentiation influences progenitor arrangement in the vertebrate neuroepithelium.

Cell division, movement and differentiation contribute to pattern formation in developing tissues. This is the case in the vertebrate neural tube, in which neurons differentiate in a characteristic pattern from a highly dynamic proliferating pseudostratified epithelium. To investigate how progenitor proliferation and differentiation affect cell arrangement and growth of the neural tube, we used experimental measurements to develop a mechanical model of the apical surface of the neuroepithelium that incorporates the effect of interkinetic nuclear movement and spatially varying rates of neuronal differentiation. Simulations predict that tissue growth and the shape of lineage-related clones of cells differ with the rate of differentiation. Growth is isotropic in regions of high differentiation, but dorsoventrally biased in regions of low differentiation. This is consistent with experimental observations. The absence of directional signalling in the simulations indicates that global mechanical constraints are sufficient to explain the observed differences in anisotropy. This provides insight into how the tissue growth rate affects cell dynamics and growth anisotropy and opens up possibilities to study the coupling between mechanics, pattern formation and growth in the neural tube.

Cooperative success in epithelial public goods games.

Cancer cells obtain mutations which rely on the production of diffusible growth factors to confer a fitness benefit. These mutations can be considered cooperative, and studied as public goods games within the framework of evolutionary game theory. The population structure, benefit function and update rule all influence the evolutionary success of cooperators. We model the evolution of cooperation in epithelial cells using the Voronoi tessellation model. Unlike traditional evolutionary graph theory, this allows us to implement global updating, for which birth and death events are spatially decoupled. We compare, for a sigmoid benefit function, the conditions for cooperation to be favoured and/or beneficial for well-mixed and structured populations. We find that when population structure is combined with global updating, cooperation is more successful than if there were local updating or the population were well-mixed. Interestingly, the qualitative behaviour for the well-mixed population and the Voronoi tessellation model is remarkably similar, but the latter case requires significantly lower incentives to ensure cooperation.

Oscillations in well-mixed, deterministic feedback systems: Beyond ring oscillators.

A ring oscillator is a system in which one species regulates the next, which regulates the next and so on until the last species regulates the first. In addition, the number of the regulations which are negative, and so result in a reduction in the regulated species, is odd, making the overall feedback in the loop negative. In ring oscillators, the probability of oscillations is maximised if the degradation rates of the species are equal. When there is more than one loop in the regulatory network, the dynamics can be more complicated. Here, a systematic way of organising the characteristic equation of ODE models of regulatory networks is provided. This facilitates the identification of Hopf bifurcations. It is shown that the probability of oscillations in non-ring systems is maximised for unequal degradation rates. For example, when there is a ring and a second ring employing a subset of the genes in the first ring, then the probability of oscillations is maximised when the species in the sub-ring degrade more slowly than those outside, for a negative feedback subring. When the sub-ring forms a positive feedback loop, the optimal degradation rates are larger for the species in the sub-ring, provided the positive feedback is not too strong. By contrast, optimal degradation rates are smaller for the species in the sub-ring, when the positive feedback is very strong. Adding a positive feedback loop to a repressilator increases the probability of oscillations, provided the positive feedback is not too strong, whereas adding a negative feedback loop decreases the probability of oscillations. The work is illustrated with numerical simulations of example systems: an autoregulatory gene model in which transcription is downregulated by the protein dimer and three-species and four-species gene regulatory network examples.

Minimum Action Path Theory Reveals the Details of Stochastic Transitions Out of Oscillatory States.

Cell state determination is the outcome of intrinsically stochastic biochemical reactions. Transitions between such states are studied as noise-driven escape problems in the chemical species space. Escape can occur via multiple possible multidimensional paths, with probabilities depending nonlocally on the noise. Here we characterize the escape from an oscillatory biochemical state by minimizing the Freidlin-Wentzell action, deriving from it the stochastic spiral exit path from the limit cycle. We also use the minimized action to infer the escape time probability density function.

The CaVß Subunit Protects the I-II Loop of the Voltage-gated Calcium Channel CaV2.2 from Proteasomal Degradation but Not Oligoubiquitination.

CaVß subunits interact with the voltage-gated calcium channel CaV2.2 on a site in the intracellular loop between domains I and II (the I-II loop). This interaction influences the biophysical properties of the channel and leads to an increase in its trafficking to the plasma membrane. We have shown previously that a mutant CaV2.2 channel that is unable to bind CaVß subunits (CaV2.2 W391A) was rapidly degraded (Waithe, D., Ferron, L., Page, K. M., Chaggar, K., and Dolphin, A. C. (2011) J. Biol. Chem. 286, 9598-9611). Here we show that, in the absence of CaVß subunits, a construct consisting of the I-II loop of CaV2.2 was directly ubiquitinated and degraded by the proteasome system. Ubiquitination could be prevented by mutation of all 12 lysine residues in the I-II loop to arginines. Including a palmitoylation motif at the N terminus of CaV2.2 I-II loop was insufficient to target it to the plasma membrane in the absence of CaVß subunits even when proteasomal degradation was inhibited with MG132 or ubiquitination was prevented by the lysine-to-arginine mutations. In the presence of CaVß subunit, the palmitoylated CaV2.2 I-II loop was protected from degradation, although oligoubiquitination could still occur, and was efficiently trafficked to the plasma membrane. We propose that targeting to the plasma membrane requires a conformational change in the I-II loop that is induced by binding of the CaVß subunit.

A theoretical framework for the regulation of Shh morphogen-controlled gene expression.

How morphogen gradients govern the pattern of gene expression in developing tissues is not well understood. Here, we describe a statistical thermodynamic model of gene regulation that combines the activity of a morphogen with the transcriptional network it controls. Using Sonic hedgehog (Shh) patterning of the ventral neural tube as an example, we show that the framework can be used together with the principled parameter selection technique of approximate Bayesian computation to obtain a dynamical model that accurately predicts tissue patterning. The analysis indicates that, for each target gene regulated by Gli, which is the transcriptional effector of Shh signalling, there is a neutral point in the gradient, either side of which altering the Gli binding affinity has opposite effects on gene expression. This explains recent counterintuitive experimental observations. The approach is broadly applicable and provides a unifying framework to explain the temporospatial pattern of morphogen-regulated gene expression.

Intrinsic Noise Profoundly Alters the Dynamics and Steady State of Morphogen-Controlled Bistable Genetic Switches.

During tissue development, patterns of gene expression determine the spatial arrangement of cell types. In many cases, gradients of secreted signalling molecules-morphogens-guide this process by controlling downstream transcriptional networks. A mechanism commonly used in these networks to convert the continuous information provided by the gradient into discrete transitions between adjacent cell types is the genetic toggle switch, composed of cross-repressing transcriptional determinants. Previous analyses have emphasised the steady state output of these mechanisms. Here, we explore the dynamics of the toggle switch and use exact numerical simulations of the kinetic reactions, the corresponding Chemical Langevin Equation, and Minimum Action Path theory to establish a framework for studying the effect of gene expression noise on patterning time and boundary position. This provides insight into the time scale, gene expression trajectories and directionality of stochastic switching events between cell states. Taking gene expression noise into account predicts that the final boundary position of a morphogen-induced toggle switch, although robust to changes in the details of the noise, is distinct from that of the deterministic system. Moreover, the dramatic increase in patterning time close to the boundary predicted from the deterministic case is substantially reduced. The resulting stochastic switching introduces differences in patterning time along the morphogen gradient that result in a patterning wave propagating away from the morphogen source with a velocity determined by the intrinsic noise. The wave sharpens and slows as it advances and may never reach steady state in a biologically relevant time. This could explain experimentally observed dynamics of pattern formation. Together the analysis reveals the importance of dynamical transients for understanding morphogen-driven transcriptional networks and indicates that gene expression noise can qualitatively alter developmental patterning.

A CaV2.1 N-terminal fragment relieves the dominant-negative inhibition by an Episodic ataxia 2 mutant.

Episodic ataxia 2 (EA2) is an autosomal dominant disorder caused by mutations in the gene CACNA1A that encodes the pore-forming CaV2.1 calcium channel subunit. The majority of EA2 mutations reported so far are nonsense or deletion/insertion mutations predicted to form truncated proteins. Heterologous expression of wild-type CaV2.1, together with truncated constructs that mimic EA2 mutants, significantly suppressed wild-type calcium channel function, indicating that the truncated protein produces a dominant-negative effect (Jouvenceau et al., 2001; Page et al., 2004). A similar finding has been shown for CaV2.2 (Raghib et al., 2001). We show here that a highly conserved sequence in the cytoplasmic N-terminus is involved in this process, for both CaV2.1 and CaV2.2 channels. Additionally, we were able to interfere with the suppressive effect of an EA2 construct by mutating key N-terminal residues within it. We postulate that the N-terminus of the truncated channel plays an essential part in its interaction with the full-length CaV2.1, which prevents the correct folding of the wild-type channel. In agreement with this, we were able to disrupt the interaction between EA2 and the full length channel by co-expressing a free N-terminal peptide.

Therapeutic activity of an interleukin-4/interleukin-13 dual antagonist on oxazolone-induced colitis in mice.

Interleukin-4 (IL-4) and IL-13 are critical drivers of immune activation and inflammation in ulcerative colitis, asthma and other diseases. Because these cytokines may have redundant function, dual targeting holds promise for achieving greater efficacy. We have recently described a bifunctional therapeutic targeting IL-4 and IL-13 developed on a novel protein scaffold, generated by combining specific binding domains in an optimal configuration using appropriate linker regions. In the current study, the bifunctional IL-4/IL-13 antagonist was evaluated in the murine oxazolone-induced colitis model, which produces disease with features of ulcerative colitis. The bifunctional IL-4/IL-13 antagonist reduced body weight loss throughout the 7-day course of the model, and ameliorated the increased colon weight and decreased colon length that accompany disease. Colon tissue gene expression was modulated in accordance with the treatment effect. Concentrations of serum amyloid P were elevated in proportion to disease severity, making it an effective biomarker. Serum concentrations of the bifunctional IL-4/IL-13 antagonist were inversely proportional to disease severity, colon tissue expression of pro-inflammatory genes, and serum amyloid P concentration. Taken together, these results define a panel of biomarkers signifying engagement of the IL-4/IL-13 pathway, confirm the T helper type 2 nature of disease in this model, and demonstrate the effectiveness of dual cytokine blockade.

The Interplay Between Splicing of Two Exon Combinations Differentially Affects Membrane Targeting and Function of Human CaV2.2.

N-type calcium channels (CaV2.2) are predominantly localized in presynaptic terminals, and are particularly important for pain transmission in the spinal cord. Furthermore, they have multiple isoforms, conferred by alternatively spliced or cassette exons, which are differentially expressed. Here, we have examined alternatively spliced exon47 variants that encode a long or short C-terminus in human CaV2.2. In the Ensembl database, all short exon47-containing transcripts were associated with the absence of exon18a, therefore, we also examined the effect of inclusion or absence of exon18a, combinatorially with the exon47 splice variants. We found that long exon47, only in the additional presence of exon18a, results in CaV2.2 currents that have a 3.6-fold greater maximum conductance than the other three combinations. In contrast, cell-surface expression of CaV2.2 in both tsA-201 cells and hippocampal neurons is increased ∼4-fold by long exon47, relative to short exon47, in either the presence or the absence of exon18a. This surprising discrepancy between trafficking and function indicates that cell-surface expression is enhanced by long exon47, independently of exon18a. However, in the presence of long exon47, exon18a mediates an additional permissive effect on CaV2.2 gating. We also investigated the single-nucleotide polymorphism in exon47 that has been linked to schizophrenia and Parkinson's disease, which we found is only non-synonymous in the short exon47 C-terminal isoform, resulting in two minor alleles. This study highlights the importance of investigating the combinatorial effects of exon inclusion, rather than each in isolation, in order to increase our understanding of calcium channel function.

Regulation of long-range BMP gradients and embryonic polarity by propagation of local calcium-firing activity.

Many amniote vertebrate species including humans can form identical twins from a single embryo, but this only occurs rarely. It has been suggested that the primitive-streak-forming embryonic region emits signals that inhibit streak formation elsewhere but the signals involved, how they are transmitted and how they act has not been elucidated. Here we show that short tracks of calcium firing activity propagate through extraembryonic tissue via gap junctions and prevent ectopic primitive streak formation in chick embryos. Cross-regulation of calcium activity and an inhibitor of primitive streak formation (Bone Morphogenetic Protein, BMP) via NF-κB and NFAT establishes a long-range BMP gradient spanning the embryo. This mechanism explains how embryos of widely different sizes can maintain positional information that determines embryo polarity. We provide evidence for similar mechanisms in two different human embryo models and in Drosophila, suggesting an ancient evolutionary origin.

Coupling dynamics of 2D Notch-Delta signalling.

Understanding pattern formation driven by cell-cell interactions has been a significant theme in cellular biology for many years. In particular, due to their implications within many biological contexts, lateral-inhibition mechanisms present in the Notch-Delta signalling pathway led to an extensive discussion between biologists and mathematicians. Deterministic and stochastic models have been developed as a consequence of this discussion, some of which address long-range signalling by considering cell protrusions reaching non-neighbouring cells. The dynamics of such signalling systems reveal intricate properties of the coupling terms involved in these models. In this work, we investigate the advantages and drawbacks of a single-parameter long-range signalling model across diverse scenarios. By employing linear and multi-scale analyses, we discover that pattern selection is not only partially explained but also depends on nonlinear effects that extend beyond the scope of these analytical techniques.

The importance of cache domains in α2δ proteins and the basis for their gabapentinoid selectivity.

In this hybrid review, we have first collected and reviewed available information on the structure and function of the enigmatic cache domains in α2δ proteins. These are organized into two double cache (dCache_1) domains, and they are present in all α2δ proteins. We have also included new data on the key function of these domains with respect to amino acid and gabapentinoid binding to the universal amino acid-binding pocket, which is present in α2δ-1 and α2δ-2. We have now identified the reason why α2δ-3 and α2δ-4 do not bind gabapentinoid drugs or amino acids with bulky side chains. In relation to this, we have determined that the bulky amino acids Tryptophan and Phenylalanine prevent gabapentin from inhibiting cell surface trafficking of α2δ-1. Together, these novel data shed further light on the importance of the cache domains in α2δ proteins.