Conditional deletion of miR-204 and miR-211 in murine retinal pigment epithelium results in retinal degeneration

MicroRNAs (miRs) are short, evolutionarily conserved non-coding RNAs that canonically downregulate expression of target genes. The miR family composed of miR-204 and miR-211 is among the most highly expressed in the retinal pigment epithelium (RPE) in both mouse and human, and also retains high sequence identity. To assess the role of this miR family in the developed mouse eye, we generated two floxed conditional knockout mouse lines crossed to the RPE65-ERT2-Cre driver mouse line to perform an RPE-specific conditional knockout of this miR family in adult mice. After Cre-mediated deletion, we observed retinal structural changes by optical coherence tomography; dysfunction and loss of photoreceptors by retinal imaging; and retinal inflammation marked by subretinal infiltration of immune cells by imaging and immunostaining. Single-cell RNA sequencing of diseased RPE and retinas showed potential miR-regulated target genes, as well as changes in non-coding RNAs in the RPE, rod photoreceptors, and Müller glia. This work thus highlights the role of miR-204 and miR-211 in maintaining RPE function and how the loss of miRs in the RPE exerts effects on the neural retina, leading to inflammation and retinal degeneration.

Restoring retinal polyunsaturated fatty acid balance and retina function by targeting ceramide in AdipoR1 deficient mice

Mutations in the adiponectin receptor 1 gene (AdipoR1) lead to retinitis pigmentosa and are associated with age-related macular degeneration (AMD). This study explores the effects of AdipoR1 gene deficiency in mice, revealing a striking decline in ω3 polyunsaturated fatty acids (PUFA), an increase in ω6 FAs, and elevated ceramides in the retina. The AdipoR1 deficiency impairs peroxisome proliferator-activated receptor α (PPARα) signaling, which is crucial for FA metabolism, particularly affecting proteins associated with FA transport and oxidation in the retina and retinal pigmented epithelium (RPE). Our lipidomic and proteomic analyses indicate changes that could affect membrane composition and viscosity through altered ω3 PUFA transport and synthesis, suggesting a potential influence of AdipoR1 on these properties. Furthermore, we noted a reduction in the Bardet-Biedl syndrome (BBS) proteins, which are crucial for forming and maintaining photoreceptor outer segments that are PUFA-enriched ciliary structures. Diminution in BBS-proteins content combined with our electron microscopic observations raises the possibility that AdipoR1 deficiency might impair ciliary function. Treatment with inhibitors of ceramide synthesis led to substantial elevation of ω3 LC-PUFAs, alleviating photoreceptor degeneration and improving retinal function. These results serve as the proof of concept for a ceramide-targeted strategy to treat retinopathies linked to PUFA deficiency, including AMD.

Distinct mouse models of Stargardt disease display differences in pharmacological targeting of ceramides and inflammatory responses

Mutations in many visual cycle enzymes in photoreceptors and retinal pigment epithelium (RPE) cells can lead to the chronic accumulation of toxic retinoid byproducts, which poison photoreceptors and the underlying RPE if left unchecked. Without a functional ATP-binding cassette, sub-family A, member 4 (ABCA4), there is an elevation of all-trans-retinal and prolonged buildup of all-trans-retinal adducts, resulting in a retinal degenerative disease known as Stargardt-1 disease. Even in this monogenic disorder, there is significant heterogeneity in the time to onset of symptoms among patients. Using a combination of molecular techniques, we studied Abca4 knockout (simulating human noncoding disease variants) and Abca4 knock-in mice (simulating human misfolded, catalytically inactive protein variants), which serve as models for Stargardt-1 disease. We compared the two strains to ascertain whether they exhibit differential responses to agents that affect cytokine signaling and/or ceramide metabolism, as alterations in either of these pathways can exacerbate retinal degenerative phenotypes. We found different degrees of responsiveness to maraviroc, a known immunomodulatory CCR5 antagonist, and to the ceramide-lowering agent AdipoRon, an agonist of the ADIPOR1 and ADIPOR2 receptors. The two strains also display different degrees of transcriptional deviation from matched WT controls. Our phenotypic comparison of the two distinct Abca4 mutant-mouse models sheds light on potential therapeutic avenues previously unexplored in the treatment of Stargardt disease and provides a surrogate assay for assessing the effectiveness for genome editing.

spinDrop: a droplet microfluidic platform to maximise single-cell sequencing information content

Droplet microfluidic methods have massively increased the throughput of single-cell sequencing campaigns. The benefit of scale-up is, however, accompanied by increased background noise when processing challenging samples and the overall RNA capture efficiency is lower. These drawbacks stem from the lack of strategies to enrich for high-quality material or specific cell types at the moment of cell encapsulation and the absence of implementable multi-step enzymatic processes that increase capture. Here we alleviate both bottlenecks using fluorescence-activated droplet sorting to enrich for droplets that contain single viable cells, intact nuclei, fixed cells or target cell types and use reagent addition to droplets by picoinjection to perform multi-step lysis and reverse transcription. Our methodology increases gene detection rates fivefold, while reducing background noise by up to half. We harness these properties to deliver a high-quality molecular atlas of mouse brain development, despite starting with highly damaged input material, and provide an atlas of nascent RNA transcription during mouse organogenesis. Our method is broadly applicable to other droplet-based workflows to deliver sensitive and accurate single-cell profiling at a reduced cost.

Stress resilience-enhancing drugs preserve tissue structure and function in degenerating retina via phosphodiesterase inhibition

Chronic, progressive retinal diseases, such as age-related macular degeneration (AMD), diabetic retinopathy, and retinitis pigmentosa, arise from genetic and environmental perturbations of cellular and tissue homeostasis. These disruptions accumulate with repeated exposures to stress over time, leading to progressive visual impairment and, in many cases, legal blindness. Despite decades of research, therapeutic options for the millions of patients suffering from these disorders remain severely limited, especially for treating earlier stages of pathogenesis when the opportunity to preserve the retinal structure and visual function is greatest. To address this urgent, unmet medical need, we employed a systems pharmacology platform for therapeutic development. Through integrative single-cell transcriptomics, proteomics, and phosphoproteomics, we identified universal molecular mechanisms across distinct models of age-related and inherited retinal degenerations, characterized by impaired physiological resilience to stress. Here, we report that selective, targeted pharmacological inhibition of cyclic nucleotide phosphodiesterases (PDEs), which serve as critical regulatory nodes that modulate intracellular second messenger signaling pathways, stabilized the transcriptome, proteome, and phosphoproteome through downstream activation of protective mechanisms coupled with synergistic inhibition of degenerative processes. This therapeutic intervention enhanced resilience to acute and chronic forms of stress in the degenerating retina, thus preserving tissue structure and function across various models of age-related and inherited retinal disease. Taken together, these findings exemplify a systems pharmacology approach to drug discovery and development, revealing a new class of therapeutics with potential clinical utility in the treatment or prevention of the most common causes of blindness.

Tim-3 adapter protein Bat3 acts as an endogenous regulator of tolerogenic dendritic cell function

Dendritic cells (DCs) sense environmental cues and adopt either an immune-stimulatory or regulatory phenotype, thereby fine-tuning immune responses. Identifying endogenous regulators that determine DC function can thus inform the development of therapeutic strategies for modulating the immune response in different disease contexts. Tim-3 plays an important role in regulating immune responses by inhibiting the activation status and the T cell priming ability of DC in the setting of cancer. Bat3 is an adaptor protein that binds to the tail of Tim-3; therefore, we studied its role in regulating the functional status of DCs. In murine models of autoimmunity (experimental autoimmune encephalomyelitis) and cancer (MC38-OVA-implanted tumor), lack of Bat3 expression in DCs alters the T cell compartment-it decreases TH1, TH17 and cytotoxic effector cells, increases regulatory T cells, and exhausted CD8+ tumor-infiltrating lymphocytes, resulting in the attenuation of autoimmunity and acceleration of tumor growth. We found that Bat3 expression levels were differentially regulated by activating versus inhibitory stimuli in DCs, indicating a role for Bat3 in the functional calibration of DC phenotypes. Mechanistically, loss of Bat3 in DCs led to hyperactive unfolded protein response and redirected acetyl-coenzyme A to increase cell intrinsic steroidogenesis. The enhanced steroidogenesis in Bat3-deficient DC suppressed T cell response in a paracrine manner. Our findings identified Bat3 as an endogenous regulator of DC function, which has implications for DC-based immunotherapies.

Inhibition of ceramide accumulation in AdipoR1-/- mice increases photoreceptor survival and improves vision

Adiponectin receptor 1 (ADIPOR1) is a lipid and glucose metabolism regulator that possesses intrinsic ceramidase activity. Mutations of the ADIPOR1 gene have been associated with nonsyndromic and syndromic retinitis pigmentosa. Here, we show that the absence of AdipoR1 in mice leads to progressive photoreceptor degeneration, significant reduction of electroretinogram amplitudes, decreased retinoid content in the retina, and reduced cone opsin expression. Single-cell RNA-Seq results indicate that ADIPOR1 encoded the most abundantly expressed ceramidase in mice and one of the 2 most highly expressed ceramidases in the human retina, next to acid ceramidase ASAH1. We discovered an accumulation of ceramides in the AdipoR1^–/– retina, likely due to insufficient ceramidase activity for healthy retina function, resulting in photoreceptor death. Combined treatment with desipramine/L-cycloserine (DC) lowered ceramide levels and exerted a protective effect on photoreceptors in AdipoR1^–/– mice. Moreover, we observed improvement in cone-mediated retinal function in the DC-treated animals. Lastly, we found that prolonged DC treatment corrected the electrical responses of the primary visual cortex to visual stimuli, approaching near-normal levels for some parameters. These results highlight the importance of ADIPOR1 ceramidase in the retina and show that pharmacological inhibition of ceramide generation can provide a therapeutic strategy for ADIPOR1-related retinopathy.

Cycling cancer persister cells arise from lineages with distinct programs

Non-genetic mechanisms have recently emerged as important drivers of cancer therapy failure¹, where some cancer cells can enter a reversible drug-tolerant persister state in response to treatment². Although most cancer persisters remain arrested in the presence of the drug, a rare subset can re-enter the cell cycle under constitutive drug treatment. Little is known about the non-genetic mechanisms that enable cancer persisters to maintain proliferative capacity in the presence of drugs. To study this rare, transiently resistant, proliferative persister population, we developed Watermelon, a high-complexity expressed barcode lentiviral library for simultaneous tracing of each cell's clonal origin and proliferative and transcriptional states. Here we show that cycling and non-cycling persisters arise from different cell lineages with distinct transcriptional and metabolic programs. Upregulation of antioxidant gene programs and a metabolic shift to fatty acid oxidation are associated with persister proliferative capacity across multiple cancer types. Impeding oxidative stress or metabolic reprogramming alters the fraction of cycling persisters. In human tumours, programs associated with cycling persisters are induced in minimal residual disease in response to multiple targeted therapies. The Watermelon system enabled the identification of rare persister lineages that are preferentially poised to proliferate under drug pressure, thus exposing new vulnerabilities that can be targeted to delay or even prevent disease recurrence.

TIM-3 restrains anti-tumour immunity by regulating inflammasome activation

T cell immunoglobulin and mucin-containing molecule 3 (TIM-3), first identified as a molecule expressed on interferon-γ producing T cells1, is emerging as an important immune-checkpoint molecule, with therapeutic blockade of TIM-3 being investigated in multiple human malignancies. Expression of TIM-3 on CD8+ T cells in the tumour microenvironment is considered a cardinal sign of T cell dysfunction; however, TIM-3 is also expressed on several other types of immune cell, confounding interpretation of results following blockade using anti-TIM-3 monoclonal antibodies. Here, using conditional knockouts of TIM-3 together with single-cell RNA sequencing, we demonstrate the singular importance of TIM-3 on dendritic cells (DCs), whereby loss of TIM-3 on DCs-but not on CD4+ or CD8+ T cells-promotes strong anti-tumour immunity. Loss of TIM-3 prevented DCs from expressing a regulatory program and facilitated the maintenance of CD8+ effector and stem-like T cells. Conditional deletion of TIM-3 in DCs led to increased accumulation of reactive oxygen species resulting in NLRP3 inflammasome activation. Inhibition of inflammasome activation, or downstream effector cytokines interleukin-1β (IL-1β) and IL-18, completely abrogated the protective anti-tumour immunity observed with TIM-3 deletion in DCs. Together, our findings reveal an important role for TIM-3 in regulating DC function and underscore the potential of TIM-3 blockade in promoting anti-tumour immunity by regulating inflammasome activation.

Inherited myeloproliferative neoplasm risk affects haematopoietic stem cells

Myeloproliferative neoplasms (MPNs) are blood cancers that are characterized by the excessive production of mature myeloid cells and arise from the acquisition of somatic driver mutations in haematopoietic stem cells (HSCs). Epidemiological studies indicate a substantial heritable component of MPNs that is among the highest known for cancers1. However, only a limited number of genetic risk loci have been identified, and the underlying biological mechanisms that lead to the acquisition of MPNs remain unclear. Here, by conducting a large-scale genome-wide association study (3,797 cases and 1,152,977 controls), we identify 17 MPN risk loci (P < 5.0 × 10-8), 7 of which have not been previously reported. We find that there is a shared genetic architecture between MPN risk and several haematopoietic traits from distinct lineages; that there is an enrichment for MPN risk variants within accessible chromatin of HSCs; and that increased MPN risk is associated with longer telomere length in leukocytes and other clonal haematopoietic states-collectively suggesting that MPN risk is associated with the function and self-renewal of HSCs. We use gene mapping to identify modulators of HSC biology linked to MPN risk, and show through targeted variant-to-function assays that CHEK2 and GFI1B have roles in altering the function of HSCs to confer disease risk. Overall, our results reveal a previously unappreciated mechanism for inherited MPN risk through the modulation of HSC function.

Cumulus provides cloud-based data analysis for large-scale single-cell and single-nucleus RNA-seq

Massively parallel single-cell and single-nucleus RNA sequencing has opened the way to systematic tissue atlases in health and disease, but as the scale of data generation is growing, so is the need for computational pipelines for scaled analysis. Here we developed Cumulus-a cloud-based framework for analyzing large-scale single-cell and single-nucleus RNA sequencing datasets. Cumulus combines the power of cloud computing with improvements in algorithm and implementation to achieve high scalability, low cost, user-friendliness and integrated support for a comprehensive set of features. We benchmark Cumulus on the Human Cell Atlas Census of Immune Cells dataset of bone marrow cells and show that it substantially improves efficiency over conventional frameworks, while maintaining or improving the quality of results, enabling large-scale studies.

Integrated scRNA-Seq Identifies Human Postnatal Thymus Seeding Progenitors and Regulatory Dynamics of Differentiating Immature Thymocytes

During postnatal life, thymopoiesis depends on the continuous colonization of the thymus by bone-marrow-derived hematopoietic progenitors that migrate through the bloodstream. The current understanding of the nature of thymic immigrants is largely based on data from pre-clinical models. Here, we employed single-cell RNA sequencing (scRNA-seq) to examine the immature postnatal thymocyte population in humans. Integration of bone marrow and peripheral blood precursor datasets identified two putative thymus seeding progenitors that varied in expression of CD7; CD10; and the homing receptors CCR7, CCR9, and ITGB7. Whereas both precursors supported T cell development, only one contributed to intrathymic dendritic cell (DC) differentiation, predominantly of plasmacytoid dendritic cells. Trajectory inference delineated the transcriptional dynamics underlying early human T lineage development, enabling prediction of transcription factor (TF) modules that drive stage-specific steps of human T cell development. This comprehensive dataset defines the expression signature of immature human thymocytes and provides a resource for the further study of human thymopoiesis.

Transcriptional Atlas of Intestinal Immune Cells Reveals that Neuropeptide α-CGRP Modulates Group 2 Innate Lymphoid Cell Responses

Signaling abnormalities in immune responses in the small intestine can trigger chronic type 2 inflammation involving interaction of multiple immune cell types. To systematically characterize this response, we analyzed 58,067 immune cells from the mouse small intestine by single-cell RNA sequencing (scRNA-seq) at steady state and after induction of a type 2 inflammatory reaction to ovalbumin (OVA). Computational analysis revealed broad shifts in both cell-type composition and cell programs in response to the inflammation, especially in group 2 innate lymphoid cells (ILC2s). Inflammation induced the expression of exon 5 of Calca, which encodes the alpha-calcitonin gene-related peptide (α-CGRP), in intestinal KLRG1⁺ ILC2s. α-CGRP antagonized KLRG1⁺ ILC2s proliferation but promoted IL-5 expression. Genetic perturbation of α-CGRP increased the proportion of intestinal KLRG1⁺ ILC2s. Our work highlights a model where α-CGRP-mediated neuronal signaling is critical for suppressing ILC2 expansion and maintaining homeostasis of the type 2 immune machinery.

Decoding human fetal liver haematopoiesis

Definitive haematopoiesis in the fetal liver supports self-renewal and differentiation of haematopoietic stem cells and multipotent progenitors (HSC/MPPs) but remains poorly defined in humans. Here, using single-cell transcriptome profiling of approximately 140,000 liver and 74,000 skin, kidney and yolk sac cells, we identify the repertoire of human blood and immune cells during development. We infer differentiation trajectories from HSC/MPPs and evaluate the influence of the tissue microenvironment on blood and immune cell development. We reveal physiological erythropoiesis in fetal skin and the presence of mast cells, natural killer and innate lymphoid cell precursors in the yolk sac. We demonstrate a shift in the haemopoietic composition of fetal liver during gestation away from being predominantly erythroid, accompanied by a parallel change in differentiation potential of HSC/MPPs, which we functionally validate. Our integrated map of fetal liver haematopoiesis provides a blueprint for the study of paediatric blood and immune disorders, and a reference for harnessing the therapeutic potential of HSC/MPPs.

A Cellular Taxonomy of the Bone Marrow Stroma in Homeostasis and Leukemia

Stroma is a poorly defined non-parenchymal component of virtually every organ with key roles in organ development, homeostasis, and repair. Studies of the bone marrow stroma have defined individual populations in the stem cell niche regulating hematopoietic regeneration and capable of initiating leukemia. Here, we use single-cell RNA sequencing (scRNA-seq) to define a cellular taxonomy of the mouse bone marrow stroma and its perturbation by malignancy. We identified seventeen stromal subsets expressing distinct hematopoietic regulatory genes spanning new fibroblastic and osteoblastic subpopulations including distinct osteoblast differentiation trajectories. Emerging acute myeloid leukemia impaired mesenchymal osteogenic differentiation and reduced regulatory molecules necessary for normal hematopoiesis. These data suggest that tissue stroma responds to malignant cells by disadvantaging normal parenchymal cells. Our taxonomy of the stromal compartment provides a comprehensive bone marrow cell census and experimental support for cancer cell crosstalk with specific stromal elements to impair normal tissue function and thereby enable emergent cancer.

Optimal-Transport Analysis of Single-Cell Gene Expression Identifies Developmental Trajectories in Reprogramming

Understanding the molecular programs that guide differentiation during development is a major challenge. Here, we introduce Waddington-OT, an approach for studying developmental time courses to infer ancestor-descendant fates and model the regulatory programs that underlie them. We apply the method to reconstruct the landscape of reprogramming from 315,000 single-cell RNA sequencing (scRNA-seq) profiles, collected at half-day intervals across 18 days. The results reveal a wider range of developmental programs than previously characterized. Cells gradually adopt either a terminal stromal state or a mesenchymal-to-epithelial transition state. The latter gives rise to populations related to pluripotent, extra-embryonic, and neural cells, with each harboring multiple finer subpopulations. The analysis predicts transcription factors and paracrine signals that affect fates and experiments validate that the TF Obox6 and the cytokine GDF9 enhance reprogramming efficiency. Our approach sheds light on the process and outcome of reprogramming and provides a framework applicable to diverse temporal processes in biology.

Hill kinetics as a noise filter: the role of transcription factor autoregulation in gene cascades

An intuition based on deterministic models of chemical kinetics is that population heterogeneity of transcription factor levels in cells is transmitted unchanged downstream to the target genes. We use a stochastic model of a two-gene cascade with a self-regulating upstream gene to show that, counter to the intuition, there is no simple mapping (bimodal to bimodal, unimodal to unimodal) between the shapes of the distributions of transcription factor numbers and target protein numbers in cells. Due to the presence of the two regulations, the system contains two nonlinear transfer functions, defined by the Hill kinetics of transcription factor binding. The transfer function of the regulator can "interfere" with the transfer function of the target, converting the bimodal input into a unimodal output or vice versa. We show that this effect can be predicted by a geometric construction. As an example application of the method, we present a case study of a system of several downstream genes of different sensitivities, controlled by a common transcription factor which also regulates its own transcription. We show that a single regulator can induce qualitatively different patterns (binary or graded) of responses to a signal in different downstream genes, depending on whether the sensitivity regions of the transfer functions of the upstream and downstream genes overlap or not. Alternatively, the same model can be interpreted as describing a single downstream gene that has different sensitivities in different cell lines due to mutations. Our model shows, therefore, a possible kinetic mechanism by which different genes can interpret the same biological signal in a different manner.

Enforcement of developmental lineage specificity by transcription factor Oct1

Embryonic stem cells co-express Oct4 and Oct1, a related protein with similar DNA-binding specificity. To study the role of Oct1 in ESC pluripotency and transcriptional control, we constructed germline and inducible-conditional Oct1-deficient ESC lines. ESCs lacking Oct1 show normal appearance, self-renewal and growth but manifest defects upon differentiation. They fail to form beating cardiomyocytes, generate neurons poorly, form small, poorly differentiated teratomas, and cannot generate chimeric mice. Upon RA-mediated differentiation, Oct1-deficient cells induce lineage-appropriate developmentally poised genes poorly while lineage-inappropriate genes, including extra-embryonic genes, are aberrantly expressed. In ESCs, Oct1 co-occupies a specific set of targets with Oct4, but does not occupy differentially expressed developmental targets. Instead, Oct1 occupies these targets as cells differentiate and Oct4 declines. These results identify a dynamic interplay between Oct1 and Oct4, in particular during the critical window immediately after loss of pluripotency when cells make the earliest developmental fate decisions.

DOI:http://dx.doi.org/10.7554/eLife.20937.001

Humans and most other animals are composed of hundreds of different types of cell, including nerve cells, muscle cells and blood cells. Despite performing many different roles, these cells all develop from a single fertilized egg, which divides to make a particular group of cells that when studied in the laboratory are called embryonic stem cells (or ESCs for short).

The ability of a cell to become a different cell type is defined as “potency”. ESCs are unique because they can specialize into any type of cell present in the adult organism, and they are therefore called “pluripotent”. However, as the embryo develops, its ESCs gradually lose their potency, and become more and more specialized. The activity of a great number of genes must be regulated during the transition from pluripotent to specialized cells, and some of the mechanisms involved in this transition are still unclear.

ESCs are known to need a gene-regulating protein called Oct4 to remain pluripotent and Shen, Kang, Shakya et al. now show that a similar protein named Oct1 is essential for their transition to becoming more specialized. When the gene for Oct1 was deleted from mouse ECSs, they behaved largely like “normal” ESCs, but could not properly mature into certain cell types such as heart and nerve cells. Molecular analyses revealed that Oct4 and Oct1 compete to regulate the activity of many common genes with opposing outcomes: Oct4 keeps ESCs pluripotent while Oct1 leads them to specialize. The Oct4 protein is abundant in ESCs and prevails over Oct1, but as the cells mature, the levels of Oct4 drop, and Oct1 takes over in the regulation of their common target genes.

Going forward, a better understanding of how ESCs become specialized will help basic research in the laboratory and allow scientists to tackle new questions about how the human body develops and how our organs work. In the longer-term, these findings might also have applications in the field of regenerative medicine, which aims to repair or replace a person’s cells, tissues or organs to improve their health.

DOI:http://dx.doi.org/10.7554/eLife.20937.002

The Hinge Region Strengthens the Nonspecific Interaction between Lac-Repressor and DNA: A Computer Simulation Study

LacI is commonly used as a model to study the protein-DNA interaction and gene regulation. The headpiece of the lac-repressor (LacI) protein is an ideal system for investigation of nonspecific binding of the whole LacI protein to DNA. The hinge region of the headpiece has been known to play a key role in the specific binding of LacI to DNA, whereas its role in nonspecific binding process has not been elucidated. Here, we report the results of explicit solvent molecular dynamics simulation and continuum electrostatic calculations suggesting that the hinge region strengthens the nonspecific interaction, accounting for up to 50% of the micro-dissociation free energy of LacI from DNA. Consequently, the rate of microscopic dissociation of LacI from DNA is reduced by 2~3 orders of magnitude in the absence of the hinge region. We find the hinge region makes an important contribution to the electrostatic energy, the salt dependence of electrostatic energy, and the number of salt ions excluded from binding of the LacI-DNA complex.

The Hinge Region Strengthens the Nonspecific Interaction between Lac-Repressor and DNA: A Computer Simulation Study

LacI is commonly used as a model to study the protein-DNA interaction and gene regulation. The headpiece of the lac-repressor (LacI) protein is an ideal system for investigation of nonspecific binding of the whole LacI protein to DNA. The hinge region of the headpiece has been known to play a key role in the specific binding of LacI to DNA, whereas its role in nonspecific binding process has not been elucidated. Here, we report the results of explicit solvent molecular dynamics simulation and continuum electrostatic calculations suggesting that the hinge region strengthens the nonspecific interaction, accounting for up to 50% of the micro-dissociation free energy of LacI from DNA. Consequently, the rate of microscopic dissociation of LacI from DNA is reduced by 2~3 orders of magnitude in the absence of the hinge region. We find the hinge region makes an important contribution to the electrostatic energy, the salt dependence of electrostatic energy, and the number of salt ions excluded from binding of the LacI-DNA complex.

Small Crowders Slow Down Kinesin-1 Stepping by Hindering Motor Domain Diffusion

The dimeric motor protein kinesin-1 moves processively along microtubules against forces of up to 7 pN. However, the mechanism of force generation is still debated. Here, we point to the crucial importance of diffusion of the tethered motor domain for the stepping of kinesin-1: small crowders stop the motor at a viscosity of 5 mPa·s-corresponding to a hydrodynamic load in the sub-fN (~10^{-4} pN) range-whereas large crowders have no impact even at viscosities above 100 mPa·s. This indicates that the scale-dependent, effective viscosity experienced by the tethered motor domain is a key factor determining kinesin's functionality. Our results emphasize the role of diffusion in the kinesin-1 stepping mechanism and the general importance of the viscosity scaling paradigm in nanomechanics.

Motion of nanoprobes in complex liquids within the framework of the length-scale dependent viscosity model

This paper deals with the recent phenomenological model of the motion of nanoscopic objects (colloidal particles, proteins, nanoparticles, molecules) in complex liquids. We analysed motion in polymer, micellar, colloidal and protein solutions and the cytoplasm of living cells using the length-scale dependent viscosity model. Viscosity monotonically approaches macroscopic viscosity as the size of the object increases and thus gives a single, coherent picture of motion at the nano and macro scale. The model includes interparticle interactions (solvent-solute), temperature and the internal structure of a complex liquid. The depletion layer ubiquitously occurring in complex liquids is also incorporated into the model. We also discuss the biological aspects of crowding in terms of the length-scale dependent viscosity model.

Method for the analysis of contribution of sliding and hopping to a facilitated diffusion of DNA-binding protein: Application to in vivo data

DNA-binding protein searches for its target, a specific site on DNA, by means of diffusion. The search process consists of many recurrent steps of one-dimensional diffusion (sliding) along the DNA chain and three-dimensional diffusion (hopping) after dissociation of a protein from the DNA chain. Here we propose a computational method that allows extracting the contribution of sliding and hopping to the search process in vivo from the measurements of the kinetics of the target search by the lac repressor in Escherichia coli [P. Hammar et al., Science 336, 1595 (2012)]. The method combines lattice Monte Carlo simulations with the Brownian excursion theory and includes explicitly steric constraints for hopping due to the helical structure of DNA. The simulation results including all experimental data reveal that the in vivo target search is dominated by sliding. The short-range hopping to the same base pair interrupts one-dimensional sliding while long-range hopping does not contribute significantly to the kinetics of the search of the target in vivo.

Transcriptional leakage versus noise: a simple mechanism of conversion between binary and graded response in autoregulated genes

We study the response of an autoregulated gene to a range of concentrations of signal molecules. We show that transcriptional leakage and noise due to translational bursting have the opposite effects. In a positively autoregulated gene, increasing the noise converts the response from graded to binary, while increasing the leakage converts the response from binary to graded. Our findings support the hypothesis that, being a common phenomenon, leaky expression may be a relatively easy way for evolutionary tuning of the type of gene response without changing the type of regulation from positive to negative.

Quantitative influence of macromolecular crowding on gene regulation kinetics

We introduce macromolecular crowding quantitatively into the model for kinetics of gene regulation in Escherichia coli. We analyse and compute the specific-site searching time for 180 known transcription factors (TFs) regulating 1300 operons. The time is between 160 s (e.g. for SoxS M_w = 12.91 kDa) and 1550 s (e.g. for PepA₆ of M_w = 329.28 kDa). Diffusion coefficients for one-dimensional sliding are between for large proteins up to for small monomers or dimers. Three-dimensional diffusion coefficients in the cytoplasm are 2 orders of magnitude larger than 1D sliding coefficients, nevertheless the sliding enhances the binding rates of TF to specific sites by 1–2 orders of magnitude. The latter effect is due to ubiquitous non-specific binding. We compare the model to experimental data for LacI repressor and find that non-specific binding of the protein to DNA is activation- and not diffusion-limited. We show that the target location rate by LacI repressor is optimized with respect to microscopic rate constant for association to non-specific sites on DNA. We analyse the effect of oligomerization of TFs and DNA looping effects on searching kinetics. We show that optimal searching strategy depends on TF abundance.

Biologistics--diffusion coefficients for complete proteome of Escherichia coli

Motivation: Biologistics provides data for quantitative analysis of transport (diffusion) processes and their spatio-temporal correlations in cells. Mobility of proteins is one of the few parameters necessary to describe reaction rates for gene regulation. Although understanding of diffusion-limited biochemical reactions in vivo requires mobility data for the largest possible number of proteins in their native forms, currently, there is no database that would contain the complete information about the diffusion coefficients (DCs) of proteins in a given cell type.

Results: We demonstrate a method for the determination of in vivo DCs for any molecule—regardless of its molecular weight, size and structure—in any type of cell. We exemplify the method with the database of in vivo DC for all proteins (4302 records) from the proteome of K12 strain of Escherichia coli, together with examples of DC of amino acids, sugars, RNA and DNA. The database follows from the scale-dependent viscosity reference curve (sdVRC). Construction of sdVRC for prokaryotic or eukaryotic cell requires ~20 in vivo measurements using techniques such as fluorescence correlation spectroscopy (FCS), fluorescence recovery after photobleaching (FRAP), nuclear magnetic resonance (NMR) or particle tracking. The shape of the sdVRC would be different for each organism, but the mathematical form of the curve remains the same. The presented method has a high predictive power, as the measurements of DCs of several inert, properly chosen probes in a single cell type allows to determine the DCs of thousands of proteins. Additionally, obtained mobility data allow quantitative study of biochemical interactions in vivo.

Contact:rholyst@ichf.edu.pl

Supplementary information:Supplementary data are available at Bioinformatics Online.

Characterization of Caulobacter crescentus FtsZ protein using dynamic light scattering

The self-assembly of the tubulin homologue FtsZ at the mid-cell is a critical step in bacterial cell division. We introduce dynamic light scattering (DLS) spectroscopy as a new method to study the polymerization kinetics of FtsZ in solution. Analysis of the DLS data indicates that the FtsZ polymers are remarkably monodisperse in length, independent of the concentrations of GTP, GDP, and FtsZ monomers. Measurements of the diffusion coefficient of the polymers demonstrate that their length is remarkably stable until the free GTP is consumed. We estimated the mean size of the FtsZ polymers within this interval of stable length to be between 9 and 18 monomers. The rates of FtsZ polymerization and depolymerization are likely influenced by the concentration of GDP, as the repeated addition of GTP to FtsZ increased the rate of polymerization and slowed down depolymerization. Increasing the FtsZ concentration did not change the size of FtsZ polymers; however, it increased the rate of the depolymerization reaction by depleting free GTP. Using transmission electron microscopy we observed that FtsZ forms linear polymers in solutions which rapidly convert to large bundles upon contact with surfaces at time scales as short as several seconds. Finally, the best studied small molecule that binds to FtsZ, PC190723, had no stabilizing effect on Caulobacter crescentus FtsZ filaments in vitro, which complements previous studies with Escherichia coli FtsZ and confirms that this class of small molecules binds Gram-negative FtsZ weakly.

Binary and graded evolution in time in a simple model of gene induction

We solve analytically the model of gene expression induction which consists of three steps: gene activation, gene products synthesis, and product degradation. The solution is given as a time-dependent probability distribution for gene products. Following the distribution in time from the inactive state of the gene to the stationary state we observe binary or graded response depending solely on the ratio r of the gene activation rate to the rate of the gene product degradation. If r << 1 the response is binary and the continuous transition from binary to graded response occurs between r=0.1 and r=1. Therefore, if binary response is observed during relaxation to steady state, then the activation rate constant must be smaller than the degradation rate constant.