RNA "COURIERs": Enabling synthetic cell-to-cell communication in human cells.

The development of molecular couriers to selectively package, export, and recover RNA molecules within human cells is a significant challenge. In this issue of Cell, Horns et al.1 introduce cellular RNA exporters, termed COURIERs, that package, secrete, and protect RNA cargo and establish the foundation for sophisticated cell-to-cell RNA communication.

What can we learn when fitting a simple telegraph model to a complex gene expression model?

In experiments, the distributions of mRNA or protein numbers in single cells are often fitted to the random telegraph model which includes synthesis and decay of mRNA or protein, and switching of the gene between active and inactive states. While commonly used, this model does not describe how fluctuations are influenced by crucial biological mechanisms such as feedback regulation, non-exponential gene inactivation durations, and multiple gene activation pathways. Here we investigate the dynamical properties of four relatively complex gene expression models by fitting their steady-state mRNA or protein number distributions to the simple telegraph model. We show that despite the underlying complex biological mechanisms, the telegraph model with three effective parameters can accurately capture the steady-state gene product distributions, as well as the conditional distributions in the active gene state, of the complex models. Some effective parameters are reliable and can reflect realistic dynamic behaviors of the complex models, while others may deviate significantly from their real values in the complex models. The effective parameters can also be applied to characterize the capability for a complex model to exhibit multimodality. Using additional information such as single-cell data at multiple time points, we provide an effective method of distinguishing the complex models from the telegraph model. Furthermore, using measurements under varying experimental conditions, we show that fitting the mRNA or protein number distributions to the telegraph model may even reveal the underlying gene regulation mechanisms of the complex models. The effectiveness of these methods is confirmed by analysis of single-cell data for E. coli and mammalian cells. All these results are robust with respect to cooperative transcriptional regulation and extrinsic noise. In particular, we find that faster relaxation speed to the steady state results in more precise parameter inference under large extrinsic noise.

Impact of variability in cell cycle periodicity on cell population dynamics.

The cell cycle consists of a series of orchestrated events controlled by molecular sensing and feedback networks that ultimately drive the duplication of total DNA and the subsequent division of a single parent cell into two daughter cells. The ability to block the cell cycle and synchronize cells within the same phase has helped understand factors that control cell cycle progression and the properties of each individual phase. Intriguingly, when cells are released from a synchronized state, they do not maintain synchronized cell division and rapidly become asynchronous. The rate and factors that control cellular desynchronization remain largely unknown. In this study, using a combination of experiments and simulations, we investigate the desynchronization properties in cervical cancer cells (HeLa) starting from the G1/S boundary following double-thymidine block. Propidium iodide (PI) DNA staining was used to perform flow cytometry cell cycle analysis at regular 8 hour intervals, and a custom auto-similarity function to assess the desynchronization and quantify the convergence to an asynchronous state. In parallel, we developed a single-cell phenomenological model the returns the DNA amount across the cell cycle stages and fitted the parameters using experimental data. Simulations of population of cells reveal that the cell cycle desynchronization rate is primarily sensitive to the variability of cell cycle duration within a population. To validate the model prediction, we introduced lipopolysaccharide (LPS) to increase cell cycle noise. Indeed, we observed an increase in cell cycle variability under LPS stimulation in HeLa cells, accompanied with an enhanced rate of cell cycle desynchronization. Our results show that the desynchronization rate of artificially synchronized in-phase cell populations can be used a proxy of the degree of variance in cell cycle periodicity, an underexplored axis in cell cycle research.

Transcriptomics and solid tumors: The next frontier in precision cancer medicine.

Transcriptomics, which encompasses assessments of alternative splicing and alternative polyadenylation, identification of fusion transcripts, explorations of noncoding RNAs, transcript annotation, and discovery of novel transcripts, is a valuable tool for understanding cancer mechanisms and identifying biomarkers. Recent advances in high-throughput technologies have enabled large-scale gene expression profiling. Importantly, RNA expression profiling of tumor tissue has been successfully used to determine clinically actionable molecular alterations. The WINTHER precision medicine clinical trial was the first prospective trial in diverse solid malignancies that assessed both genomics and transcriptomics to match treatments to specific molecular alterations. The use of transcriptome analysis in WINTHER and other trials increased the number of targetable -omic changes compared to genomic profiling alone. Other applications of transcriptomics involve the evaluation of tumor and circulating noncoding RNAs as predictive and prognostic biomarkers, the improvement of risk stratification by the use of prognostic and predictive multigene assays, the identification of fusion transcripts that drive tumors, and an improved understanding of the impact of DNA changes as some genomic alterations are silenced at the RNA level. Finally, RNA sequencing and gene expression analysis have been incorporated into clinical trials to identify markers predicting response to immunotherapy. Many issues regarding the complexity of the analysis, its reproducibility and variability, and the interpretation of the results still need to be addressed. The integration of transcriptomics with genomics, proteomics, epigenetics, and tumor immune profiling will improve biomarker discovery and our understanding of disease mechanisms and, thereby, accelerate the implementation of precision oncology.

Plasmonic LAMP: Improving the Detection Specificity and Sensitivity for SARS-CoV-2 by Plasmonic Sensing of Isothermally Amplified Nucleic Acids.

The ability to detect pathogens specifically and sensitively is critical to combat infectious diseases outbreaks and pandemics. Colorimetric assays involving loop-mediated isothermal amplification (LAMP) provide simple readouts yet suffer from the intrinsic non-template amplification. Herein, a highly specific and sensitive assay relying on plasmonic sensing of LAMP amplicons via DNA hybridization, termed as plasmonic LAMP, is developed for the severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) RNA detection. This work has two important advances. First, gold and silver (Au-Ag) alloy nanoshells are developed as plasmonic sensors that have 4-times stronger extinction in the visible wavelengths and give a 20-times lower detection limit for oligonucleotides over Au counterparts. Second, the integrated method allows cutting the complex LAMP amplicons into short repeats that are amendable for hybridization with oligonucleotide-functionalized Au-Ag nanoshells. In the SARS-CoV-2 RNA detection, plasmonic LAMP takes ≈75 min assay time, achieves a detection limit of 10 copies per reaction, and eliminates the contamination from non-template amplification. It also shows better detection specificity and sensitivity over commercially available LAMP kits due to the additional sequence identification. This work opens a new route for LAMP amplicon detection and provides a method for virus testing at its early representation.

Uncoupling gene expression noise along the central dogma using genome engineered human cell lines.

Eukaryotic protein synthesis is an inherently stochastic process. This stochasticity stems not only from variations in cell content between cells but also from thermodynamic fluctuations in a single cell. Ultimately, these inherently stochastic processes manifest as noise in gene expression, where even genetically identical cells in the same environment exhibit variation in their protein abundances. In order to elucidate the underlying sources that contribute to gene expression noise, we quantify the contribution of each step within the process of protein synthesis along the central dogma. We uncouple gene expression at the transcriptional, translational, and post-translational level using custom engineered circuits stably integrated in human cells using CRISPR. We provide a generalized framework to approximate intrinsic and extrinsic noise in a population of cells expressing an unbalanced two-reporter system. Our decomposition shows that the majority of intrinsic fluctuations stem from transcription and that coupling the two genes along the central dogma forces the fluctuations to propagate and accumulate along the same path, resulting in increased observed global correlation between the products.

Coevolutionary Couplings Unravel PAM-Proximal Constraints of CRISPR-SpCas9.

The clustered regularly interspaced short palindromic repeats (CRISPR) system, an immune system analog found in prokaryotes, allows a single-guide RNA to direct a CRISPR-associated protein (Cas) with combined helicase and nuclease activity to DNA. The presence of a specific protospacer adjacent motif (PAM) next to the DNA target site plays a crucial role in determining both efficacy and specificity of gene editing. Herein, we introduce a coevolutionary framework to computationally unveil nonobvious molecular interactions in CRISPR systems and experimentally probe their functional role. Specifically, we use direct coupling analysis, a statistical inference framework used to infer direct coevolutionary couplings, in the context of protein/nucleic acid interactions. Applied to Streptococcus pyogenes Cas9, a Hamiltonian metric obtained from coevolutionary relationships reveals, to our knowledge, novel PAM-proximal nucleotide preferences at the seventh position of S. pyogenes Cas9 PAM (5'-NGRNNNT-3'), which was experimentally confirmed by in vitro and functional assays in human cells. We show that coevolved and conserved interactions point to specific clues toward rationally engineering new generations of Cas9 systems and may eventually help decipher the diversity of this family of proteins.

Regulating the Uptake of Viral Nanoparticles in Macrophage and Cancer Cells via a pH Switch.

Controlling the uptake of nanomaterials into phagocytes is a challenging problem. We describe an approach to inhibit the cellular uptake by macrophages and HeLa cells of nanoparticles derived from bacteriophage Qß by conjugating negatively charged terminal hexanoic acid moieties onto its surface. Additionally, we show hydrazone linkers can be installed between the surface of Qß and the terminal hexanoic acid moieties, resulting in a pH-responsive conjugate that, in acidic conditions, can release the terminal hexanoic acid moiety and allow for the uptake of the Qß nanoparticle. The installation of the "pH switch" did not change the structure-function properties of the hexanoic acid moiety and the uptake of the Qß conjugates by macrophages.

Nitroxyl Modified Tobacco Mosaic Virus as a Metal-Free High-Relaxivity MRI and EPR Active Superoxide Sensor.

Superoxide overproduction is known to occur in multiple disease states requiring critical care; yet, noninvasive detection of superoxide in deep tissue remains a challenge. Herein, we report a metal-free magnetic resonance imaging (MRI) and electron paramagnetic resonance (EPR) active contrast agent prepared by "click conjugating" paramagnetic organic radical contrast agents (ORCAs) to the surface of tobacco mosaic virus (TMV). While ORCAs are known to be reduced in vivo to an MRI/EPR silent state, their oxidation is facilitated specifically by reactive oxygen species-in particular, superoxide-and are largely unaffected by peroxides and molecular oxygen. Unfortunately, single molecule ORCAs typically offer weak MRI contrast. In contrast, our data confirm that the macromolecular ORCA-TMV conjugates show marked enhancement for T1 contrast at low field (<3.0 T) and T2 contrast at high field (9.4 T). Additionally, we demonstrated that the unique topology of TMV allows for a "quenchless fluorescent" bimodal probe for concurrent fluorescence and MRI/EPR imaging, which was made possible by exploiting the unique inner and outer surface of the TMV nanoparticle. Finally, we show TMV-ORCAs do not respond to normal cellular respiration, minimizing the likelihood for background, yet still respond to enzymatically produced superoxide in complicated biological fluids like serum.

Guide RNA engineering for versatile Cas9 functionality.

The Clustered Regularly Interspaced Short Palindromic Repeats system allows a single guide RNA (sgRNA) to direct a protein with combined helicase and nuclease activity to the DNA. Streptococcus pyogenes Cas9 (SpCas9), a CRISPR-associated protein, has revolutionized our ability to probe and edit the human genome in vitro and in vivo Arguably, the true modularity of the Cas9 platform is conferred through the ease of sgRNA programmability as well as the degree of modifications the sgRNA can tolerate without compromising its association with SpCas9 and function. In this review, we focus on the properties and recent engineering advances of the sgRNA component in Cas9-mediated genome targeting.

Discriminating direct and indirect connectivities in biological networks.

Reverse engineering of biological pathways involves an iterative process between experiments, data processing, and theoretical analysis. Despite concurrent advances in quality and quantity of data as well as computing resources and algorithms, difficulties in deciphering direct and indirect network connections are prevalent. Here, we adopt the notions of abstraction, emulation, benchmarking, and validation in the context of discovering features specific to this family of connectivities. After subjecting benchmark synthetic circuits to perturbations, we inferred the network connections using a combination of nonparametric single-cell data resampling and modular response analysis. Intriguingly, we discovered that recovered weights of specific network edges undergo divergent shifts under differential perturbations, and that the particular behavior is markedly different between topologies. Our results point to a conceptual advance for reverse engineering beyond weight inference. Investigating topological changes under differential perturbations may address the longstanding problem of discriminating direct and indirect connectivities in biological networks.

CRISPR-based self-cleaving mechanism for controllable gene delivery in human cells.

Controllable gene delivery via vector-based systems remains a formidable challenge in mammalian synthetic biology and a desirable asset in gene therapy applications. Here, we introduce a methodology to control the copies and residence time of a gene product delivered in host human cells but also selectively disrupt fragments of the delivery vehicle. A crucial element of the proposed system is the CRISPR protein Cas9. Upon delivery, Cas9 guided by a custom RNA sequence cleaves the delivery vector at strategically placed targets thereby inactivating a co-expressed gene of interest. Importantly, using experiments in human embryonic kidney cells, we show that specific parameters of the system can be adjusted to fine-tune the delivery properties. We envision future applications in complex synthetic biology architectures, gene therapy and trace-free delivery.

MiR-192-Mediated Positive Feedback Loop Controls the Robustness of Stress-Induced p53 Oscillations in Breast Cancer Cells.

The p53 tumor suppressor protein plays a critical role in cellular stress and cancer prevention. A number of post-transcriptional regulators, termed microRNAs, are closely connected with the p53-mediated cellular networks. While the molecular interactions among p53 and microRNAs have emerged, a systems-level understanding of the regulatory mechanism and the role of microRNAs-forming feedback loops with the p53 core remains elusive. Here we have identified from literature that there exist three classes of microRNA-mediated feedback loops revolving around p53, all with the nature of positive feedback coincidentally. To explore the relationship between the cellular performance of p53 with the microRNA feedback pathways, we developed a mathematical model of the core p53-MDM2 module coupled with three microRNA-mediated positive feedback loops involving miR-192, miR-34a, and miR-29a. Simulations and bifurcation analysis in relationship to extrinsic noise reproduce the oscillatory behavior of p53 under DNA damage in single cells, and notably show that specific microRNA abrogation can disrupt the wild-type cellular phenotype when the ubiquitous cell-to-cell variability is taken into account. To assess these in silico results we conducted microRNA-perturbation experiments in MCF7 breast cancer cells. Time-lapse microscopy of cell-population behavior in response to DNA double-strand breaks, together with image classification of single-cell phenotypes across a population, confirmed that the cellular p53 oscillations are compromised after miR-192 perturbations, matching well with the model predictions. Our study via modeling in combination with quantitative experiments provides new evidence on the role of microRNA-mediated positive feedback loops in conferring robustness to the system performance of stress-induced response of p53.

Synthetic mammalian transgene negative autoregulation.

Biological networks contain overrepresented small-scale topologies, typically called motifs. A frequently appearing motif is the transcriptional negative-feedback loop, where a gene product represses its own transcription. Here, using synthetic circuits stably integrated in human kidney cells, we study the effect of negative-feedback regulation on cell-wide (extrinsic) and gene-specific (intrinsic) sources of uncertainty. We develop a theoretical approach to extract the two noise components from experiments and show that negative feedback results in significant total noise reduction by reducing extrinsic noise while marginally increasing intrinsic noise. We compare the results to simple negative regulation, where a constitutively transcribed transcription factor represses a reporter protein. We observe that the control architecture also reduces the extrinsic noise but results in substantially higher intrinsic fluctuations. We conclude that negative feedback is the most efficient way to mitigate the effects of extrinsic fluctuations by a sole regulatory wiring.

NeurostimML: a machine learning model for predicting neurostimulation-induced tissue damage.

Objective. The safe delivery of electrical current to neural tissue depends on many factors, yet previous methods for predicting tissue damage rely on only a few stimulation parameters. Here, we report the development of a machine learning approach that could lead to a more reliable method for predicting electrical stimulation-induced tissue damage by incorporating additional stimulation parameters.Approach. A literature search was conducted to build an initial database of tissue response information after electrical stimulation, categorized as either damaging or non-damaging. Subsequently, we used ordinal encoding and random forest for feature selection, and investigated four machine learning models for classification: Logistic Regression, K-nearest Neighbor, Random Forest, and Multilayer Perceptron. Finally, we compared the results of these models against the accuracy of the Shannon equation.Main Results. We compiled a database with 387 unique stimulation parameter combinations collected from 58 independent studies conducted over a period of 47 years, with 195 (51%) categorized as non-damaging and 190 (49%) categorized as damaging. The features selected for building our model with a Random Forest algorithm were: waveform shape, geometric surface area, pulse width, frequency, pulse amplitude, charge per phase, charge density, current density, duty cycle, daily stimulation duration, daily number of pulses delivered, and daily accumulated charge. The Shannon equation yielded an accuracy of 63.9% using akvalue of 1.79. In contrast, the Random Forest algorithm was able to robustly predict whether a set of stimulation parameters was classified as damaging or non-damaging with an accuracy of 88.3%.Significance. This novel Random Forest model can facilitate more informed decision making in the selection of neuromodulation parameters for both research studies and clinical practice. This study represents the first approach to use machine learning in the prediction of stimulation-induced neural tissue damage, and lays the groundwork for neurostimulation driven by machine learning models.

STAT6-targeting antisense oligonucleotides against solitary fibrous tumor.

Solitary fibrous tumor (SFT) is a rare, non-hereditary soft tissue sarcoma thought to originate from fibroblastic mesenchymal stem cells. The etiology of SFT is thought to be due to an environmental intrachromosomal gene fusion between NGFI-A-binding protein 2 (NAB2) and signal transducer and activator protein 6 (STAT6) genes on chromosome 12, wherein the activation domain of STAT6 is fused with the DNA-binding domain of NAB2 resulting in the oncogenesis of SFT. All NAB2-STAT6 fusion variations discovered in SFTs contain the C-terminal of STAT6 transcript, and thus can serve as target site for antisense oligonucleotides (ASOs)-based therapies. Indeed, our in vitro studies show the STAT6 3' untranslated region (UTR)-targeting ASO (ASO 993523) was able to reduce expression of NAB2-STAT6 fusion transcripts in multiple SFT cell models with high efficiency (half-maximal inhibitory concentration: 116-300 nM). Encouragingly, in vivo treatment of SFT patient-derived xenograft mouse models with ASO 993523 resulted in acceptable tolerability profiles, reduced expression of NAB2-STAT6 fusion transcripts in xenograft tissues (21.9%), and, importantly, reduced tumor growth (32.4% decrease in tumor volume compared with the untreated control). Taken together, our study established ASO 993523 as a potential agent for the treatment of SFTs.

Clinical gene therapy development for the central nervous system: Candidates and challenges for AAVs.

Many diseases affecting the central nervous system (CNS) are deadly but less understood, leading to impaired mental and motor capabilities and poor patient prospects. Gene therapy is a promising therapeutic modality for correcting many genetic disorders, expanding in breadth and scope with further advances. This review summarizes the candidate CNS disorders for gene therapy, mechanisms of gene therapy, and recent clinical advances and limitations of gene therapy in CNS disorders. We highlight that improving delivery across CNS barriers, safety, monitoring techniques, and multiplexing therapies are predominant factors in advancing long-term outcomes from gene therapy.

NeurostimML: A machine learning model for predicting neurostimulation-induced tissue damage.

Objective: The safe delivery of electrical current to neural tissue depends on many factors, yet previous methods for predicting tissue damage rely on only a few stimulation parameters. Here, we report the development of a machine learning approach that could lead to a more reliable method for predicting electrical stimulation-induced tissue damage by incorporating additional stimulation parameters. Approach: A literature search was conducted to build an initial database of tissue response information after electrical stimulation, categorized as either damaging or non-damaging. Subsequently, we used ordinal encoding and random forest for feature selection, and investigated four machine learning models for classification: Logistic Regression, K-nearest Neighbor, Random Forest, and Multilayer Perceptron. Finally, we compared the results of these models against the accuracy of the Shannon equation. Main Results: We compiled a database with 387 unique stimulation parameter combinations collected from 58 independent studies conducted over a period of 47 years, with 195 (51%) categorized as non-damaging and 190 (49%) categorized as damaging. The features selected for building our model with a Random Forest algorithm were: waveform shape, geometric surface area, pulse width, frequency, pulse amplitude, charge per phase, charge density, current density, duty cycle, daily stimulation duration, daily number of pulses delivered, and daily accumulated charge. The Shannon equation yielded an accuracy of 63.9% using a k value of 1.79. In contrast, the Random Forest algorithm was able to robustly predict whether a set of stimulation parameters was classified as damaging or non-damaging with an accuracy of 88.3%. Significance: This novel Random Forest model can facilitate more informed decision making in the selection of neuromodulation parameters for both research studies and clinical practice. This study represents the first approach to use machine learning in the prediction of stimulation-induced neural tissue damage, and lays the groundwork for neurostimulation driven by machine learning models.

Reduction of Tumor Growth with RNA-Targeting Treatment of the NAB2-STAT6 Fusion Transcript in Solitary Fibrous Tumor Models.

Solitary fibrous tumor (SFT) is a rare soft-tissue sarcoma. This nonhereditary cancer is the result of an environmental intrachromosomal gene fusion between NAB2 and STAT6 on chromosome 12, which fuses the activation domain of STAT6 with the repression domain of NAB2. Currently there is not an approved chemotherapy regimen for SFTs. The best response on available pharmaceuticals is a partial response or stable disease for several months. The purpose of this study is to investigate the potential of RNA-based therapies for the treatment of SFTs. Specifically, in vitro SFT cell models were engineered to harbor the characteristic NAB2-STAT6 fusion using the CRISPR/SpCas9 system. Cell migration as well as multiple cancer-related signaling pathways were increased in the engineered cells as compared to the fusion-absent parent cells. The SFT cell models were then used for evaluating the targeting efficacies of NAB2-STAT6 fusion-specific antisense oligonucleotides (ASOs) and CRISPR/CasRx systems. Our results showed that fusion specific ASO treatments caused a 58% reduction in expression of fusion transcripts and a 22% reduction in cell proliferation after 72 h in vitro. Similarly, the AAV2-mediated CRISPR/CasRx system led to a 59% reduction in fusion transcript expressions in vitro, and a 55% reduction in xenograft growth after 29 days ex vivo.

Non-viral nitric oxide-based gene therapy improves perfusion and liposomal doxorubicin sonopermeation in neuroblastoma models.

Neuroblastoma (NB) is a pediatric malignancy that accounts for 15% of cancer-related childhood mortality. High-risk NB requires an aggressive chemoradiotherapy regimen that causes significant off-target toxicity. Despite this invasive treatment, many patients either relapse or do not respond adequately. Recent studies suggest that improving tumor perfusion can enhance drug accumulation and distribution within the tumor tissue, potentially augmenting treatment effects without inflicting systemic toxicity. Accordingly, methods that transiently increase tumor perfusion prior to treatment may help combat this disease. Here, we show the use of gene therapy to confer inducible nitric oxide synthase (iNOS) expression solely in the tumor space, using focused ultrasound targeting. NOS catalyzes the reaction that generates nitric oxide (NO), a potent endogenous vasodilator. This study reports the development of a targeted non-viral image-guided platform to deliver iNOS-expressing plasmid DNA (pDNA) to vascular endothelial cells encasing tumor blood vessels. Following transfection, longitudinal quantitative contrast-enhanced ultrasound (qCEUS) imaging revealed an increase in tumor perfusion over 72 h, attributed to elevated intratumoral iNOS expression. Methods: To construct a gene delivery vector, cationic ultrasound-responsive agents (known as "microbubbles") were employed to carry pDNA in circulation and transfect tumor vascular endothelial cells in vivo using focused ultrasound (FUS) energy. This was followed by liposomal doxorubicin (L-DOX) treatment. The post-transfection tumor response was monitored longitudinally using qCEUS imaging to determine relative changes in blood volumes and perfusion rates. After therapy, ex vivo analysis of tumors was performed to examine the bioeffects associated with iNOS expression. Results: By combining FUS therapy with cationic ultrasound contrast agents (UCAs), we achieved selective intratumoral transfection of pDNA encoding the iNOS enzyme. While transitory, the degree of expression was sufficient to induce significant increases in tumoral perfusion, to appreciably enhance the chemotherapeutic payload and to extend survival time in an orthotopic xenograft model. Conclusion: We have demonstrated the ability of a novel targeted non-viral gene therapy strategy to enhance tumor perfusion and improve L-DOX delivery to NB xenografts. While our results demonstrate that transiently increasing tumor perfusion improves liposome-encapsulated chemotherapeutic uptake and distribution, we expect that our iNOS gene delivery paradigm can also significantly improve radio and immunotherapies by increasing the delivery of radiosensitizers and immunomodulators, potentially improving upon current NB treatment without concomitant adverse effects. Our findings further suggest that qCEUS imaging can effectively monitor changes in tumor perfusion in vivo, allowing the identification of an ideal time-point to administer therapy.