Search

Inherited retinopathies are devastating diseases that in most cases lack treatment options. Disease-modifying therapies that mitigate pathophysiology regardless of the underlying genetic lesion are desirable due to the diversity of mutations found in such diseases. We tested a systems pharmacology-based strategy that suppresses intracellular cAMP and Ca2+ activity via G protein-coupled receptor (GPCR) modulation using tamsulosin, metoprolol, and bromocriptine coadministration. The treatment improves cone photoreceptor function and slows degeneration in Pde6βrd10 and RhoP23H/WT retinitis pigmentosa mice. Cone degeneration is modestly mitigated after a 7-month-long drug infusion in PDE6A-/- dogs. The treatment also improves rod pathway function in an Rpe65-/- mouse model of Leber congenital amaurosis but does not protect from cone degeneration. RNA-sequencing analyses indicate improved metabolic function in drug-treated Rpe65-/- and rd10 mice. Our data show that catecholaminergic GPCR drug combinations that modify second messenger levels via multiple receptor actions provide a potential disease-modifying therapy against retinal degeneration., (© 2024. The Author(s).)

MicroRNAs (miRs) are short, evolutionarily conserved noncoding RNAs that canonically downregulate expression of target genes. The miR family composed of miR-204 and miR-211 is among the most highly expressed miRs 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 KO mouse lines crossed to the RPE65-ERT2-Cre driver mouse line to perform an RPE-specific conditional KO 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 noncoding 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., Competing Interests: Conflict of interest KP is a consultant for Polgenix Inc and serves on the Scientific Advisory Board at Hyperion Eye Ltd. All other authors have declared that they have no conflict of interest with the contents of this article., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Mutations in the adiponectin receptor 1 gene (AdipoR1) lead to retinitis pigmentosa and are associated with age-related macular degeneration. 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 fatty acids, and elevated ceramides in the retina. The AdipoR1 deficiency impairs peroxisome proliferator-activated receptor α signaling, which is crucial for FA metabolism, particularly affecting proteins associated with FA transport and oxidation in the retina and retinal pigmented epithelium. 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 proteins, which are crucial for forming and maintaining photoreceptor outer segments that are PUFA-enriched ciliary structures. Diminution in Bardet-Biedl syndrome-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 age-related macular degeneration., Competing Interests: Conflict of interest K. P. is a consultant for Polgenix Inc and serves on the Scientific Advisory Board at Hyperion Eye Ltd. D. S.-K. is a scientific consultant for Visgenx, Inc. All other authors declare that they have no conflict of interests with the contents of this article., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

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., Competing Interests: Competing interests statement:K.P. is a consultant for Polgenix, Inc. and serves on the Scientific Advisory Board at Hyperion Eye, Ltd. All other authors have declared that no conflict of interest exists.

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., (© 2023. Springer Nature Limited.)

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.

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 T H 1, T H 17 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.

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.

Haematopoiesis in the bone marrow (BM) maintains blood and immune cell production throughout postnatal life. Haematopoiesis first emerges in human BM at 11-12 weeks after conception 1,2 , yet almost nothing is known about how fetal BM (FBM) evolves to meet the highly specialized needs of the fetus and newborn. Here we detail the development of FBM, including stroma, using multi-omic assessment of mRNA and multiplexed protein epitope expression. We find that the full blood and immune cell repertoire is established in FBM in a short time window of 6-7 weeks early in the second trimester. FBM promotes rapid and extensive diversification of myeloid cells, with granulocytes, eosinophils and dendritic cell subsets emerging for the first time. The substantial expansion of B lymphocytes in FBM contrasts with fetal liver at the same gestational age. Haematopoietic progenitors from fetal liver, FBM and cord blood exhibit transcriptional and functional differences that contribute to tissue-specific identity and cellular diversification. Endothelial cell types form distinct vascular structures that we show are regionally compartmentalized within FBM. Finally, we reveal selective disruption of B lymphocyte, erythroid and myeloid development owing to a cell-intrinsic differentiation bias as well as extrinsic regulation through an altered microenvironment in Down syndrome (trisomy 21)., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

Non-genetic mechanisms have recently emerged as important drivers of cancer therapy failure 1 , where some cancer cells can enter a reversible drug-tolerant persister state in response to treatment 2 . 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., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

T cell immunoglobulin and mucin-containing molecule 3 (TIM-3), first identified as a molecule expressed on interferon-γ producing T cells 1 , 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.

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 cancers 1 . 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.

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.

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., Competing Interests: Declaration of Interests A.R. is an Scientific Advistory Board member of Thermo Fisher Scientific, Neogene Therapeutics, Asimov, and Syros Pharmaceuticals; an equity holder of Immunitas; and a founder of and equity holder in Celsius Therapeutics. The remaining authors declare no competing interests., (Copyright © 2020 Elsevier Inc. All rights reserved.)

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., (Copyright © 2019 Elsevier Inc. All rights reserved.)

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.

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., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Over 90% of genetic variants associated with complex human traits map to non-coding regions, but little is understood about how they modulate gene regulation in health and disease. One possible mechanism is that genetic variants affect the activity of one or more cis-regulatory elements leading to gene expression variation in specific cell types. To identify such cases, we analyzed ATAC-seq and RNA-seq profiles from stimulated primary CD4 + T cells in up to 105 healthy donors. We found that regions of accessible chromatin (ATAC-peaks) are co-accessible at kilobase and megabase resolution, consistent with the three-dimensional chromatin organization measured by in situ Hi-C in T cells. Fifteen percent of genetic variants located within ATAC-peaks affected the accessibility of the corresponding peak (local-ATAC-QTLs). Local-ATAC-QTLs have the largest effects on co-accessible peaks, are associated with gene expression and are enriched for autoimmune disease variants. Our results provide insights into how natural genetic variants modulate cis-regulatory elements, in isolation or in concert, to influence gene expression.

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.

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.

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 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.

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., (Copyright © 2015 Elsevier B.V. All rights reserved.)

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.

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.

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 Mw = 12.91 kDa) and 1550 s (e.g. for PepA6 of Mw = 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.

Taylor Dispersion Analysis (TDA) has been performed for analytes moving at high flow rates in long, coiled capillaries. A thin injection zone of the analyte is stretched by the flow and final distribution of concentration of the analyte at the end of the capillary has the gaussian shape. The high flow rates in coiled capillary generate vortices. They convectively mix the analyte across the capillary. This mixing reduces the width of the gaussian distribution several times in comparison to the width obtained in a straight capillary in standard TDA. We have determined an empirical, scaling equation for the width as a function of the flow rate, molecular diffusion coefficient of the analyte, viscosity of the carrier phase, internal radius of the cylindrical capillary, and external radius of the coiled capillary. This equation can be used for different sizes of capillaries in a wide range of parameters without an additional calibration procedure. Our experimental results of flow in the coiled capillary could not be explained by current models based on approximate solutions of the Navier-Stokes equation. We applied the technique to determine the diffusion coefficients of the following analytes: salts, drugs, single amino acids, peptides (from dipeptides to hexapeptides), and proteins.

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.

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.

Bimodality of gene expression, as a mechanism contributing to phenotypic diversity, enhances the survival of cells in a fluctuating environment. To date, the bimodal response of a gene regulatory system has been attributed to the cooperativity of transcription factor binding or to feedback loops. It has remained unclear whether noncooperative binding of transcription factors can give rise to bimodality in an open-loop system. We study a theoretical model of gene expression in a two-step cascade (a deterministically monostable system) in which the regulatory gene produces transcription factors that have a nonlinear effect on the activity of the target gene. We show that a unimodal distribution of transcription factors over the cell population can generate a bimodal steady-state output without cooperative transcription factor binding. We introduce a simple method of geometric construction that allows one to predict the onset of bimodality. The construction only involves the parameters of bursting of the regulatory gene and the dose-response curve of the target gene. Using this method, we show that the gene expression may switch between unimodal and bimodal as the concentration of inducers or corepressors is varied. These findings may explain the experimentally observed bimodal response of cascades consisting of a fluorescent protein reporter controlled by the tetracycline repressor. The geometric construction provides a useful tool for designing experiments and for interpretation of their results. Our findings may have important implications for understanding the strategies adopted by cell populations to survive in changing environments.

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.

We measured the viscosity of poly(ethylene glycol) (PEG 6000, 12,000, 20,000) in water using capillary electrophoresis and fluorescence correlation spectroscopy with nanoscopic probes of different diameters (from 1.7 to 114 nm). For a probe of diameter smaller than the radius of gyration of PEG (e.g. rhodamine B or lyzozyme) the measured nanoviscosity was orders of magnitude smaller than the macroviscosity. For sizes equal to (or larger than) the polymer radius of gyration, macroscopic value of viscosity was measured. A mathematical relation for macro and nanoviscosity was found as a function of PEG radius of gyration, R(g), correlation length in semi-dilute solution, xi, and probe size, R. For R < R(g), the nanoviscosity (normalized by water viscosity) is given by exp(b(R/xi)a), and for R > R(g), both nano and macroviscosity follow the same curve, exp(b(R/xi)a), where a and b are two constants close to unity. This mathematical relation was shown to equally well describe rhodamine (of size 1.7 nm) in PEG 20,000 and the macroviscosity of PEG 8,000,000, whose radius of gyration exceeds 200 nm. Additionally, for the smallest probes (rhodamine B and lysozyme) we have verified, using capillary electrophoresis and fluorescence correlation spectroscopy, that the Stokes-Einstein (SE) relation holds, providing that we use a size-dependent viscosity in the formula. The SE relation is correct even in PEG solutions of very high viscosity (three orders of magnitude larger than that of water).

Understanding a biological module involves recognition of its structure and the dynamics of its principal components. In this report we present an analysis of the dynamics of the repression module within the regulation of the trp operon in Escherichia coli. We combine biochemical data for reaction rate constants for the trp repressor binding to trp operator and in vivo data of a number of tryptophan repressors (TrpRs) that bind to the operator. The model of repression presented in this report greatly differs from previous mathematical models. One, two or three TrpRs can bind to the operator and repress the transcription. Moreover, reaction rates for detachment of TrpRs from the operator strongly depend on tryptophan (Trp) concentration, since Trp can also bind to the repressor-operator complex and stabilize it. From the mathematical modeling and analysis of reaction rates and equilibrium constants emerges a high-quality, accurate and effective module of trp repression. This genetic switch responds accurately to fast consumption of Trp from the interior of a cell. It switches with minimal dispersion when the concentration of Trp drops below a thousand molecules per cell.

Little is known about the role of protists in the saturated subsurface. Porous media microcosms, containing bacteria and protists, were used to determine whether flagellates from an organically contaminated aquifer could substantively affect the number of free-living bacteria (FLB). When flagellates were present, the 3-40% maximum breakthrough of fluorescently labelled FLB injected into the microcosms was much lower than the 60-130% observed for killed controls. Grazing and clearance rates (3-27 FLB flag-1 h-1 and 12-23 nl flag-1 h-1, respectively) calculated from the data were in the range reported for flagellates in other aqueous environments. The data provide evidence that flagellate bacterivory is an important control on groundwater FLB populations.

The effects of particulate matter from reconstituted carbenicillin disodium on cultures of human macrophages was studied. Particulate matter from commercially available carbenicillin disodium was collected on 0.1 micrometer membrane filters. Monocytes/macrophages obtained from human peripheral blood were cultivated on the particle bearing surface of the filter in culture medium. The cells were incubated in tissue cultures dishes at a concentration of 2 x 10(6) cells per dish. Experimental and control cultures were observed over six days for changes in morphology, cell membrane integrity, phagocytic capacity, and viability. Exposure of the cells to the particles was not cytotoxic. Cell membrane integrity remained intact as shown by no increase of LDH in the extracellular fluid. The phagocytic capacity of the cell system was not impaired as assessed by observing the intracellular levels of glucose-6-phosphate dehydrogenase. Morphological changes as monitored by phase contrast and bright field microscopy were consistent with increased cellular activity. Particulate matter from carbenicillin disodium may not be toxic to human macrophages in vitro for up to six days. At the concentrations used in this study, particulate matter did not lower the phagocytic capacity of human macrophages in vitro.