Dynamic balance between vesicle transport and microtubule growth enables neurite outgrowth

From transcriptomics to digital twins of organ function

Cell level functions underlie tissue and organ physiology. Gene expression patterns offer extensive views of the pathways and processes within and between cells. Single cell transcriptomics provides detailed information on gene expression within cells, cell types, subtypes and their relative proportions in organs. Functional pathways can be scalably connected to physiological functions at the cell and organ levels. Integrating experimentally obtained gene expression patterns with prior knowledge of pathway interactions enables identification of networks underlying whole cell functions such as growth, contractility, and secretion. These pathways can be computationally modeled using differential equations to simulate cell and organ physiological dynamics regulated by gene expression changes. Such computational systems can be thought of as parts of digital twins of organs. Digital twins, at the core, need computational models that represent in detail and simulate how dynamics of pathways and networks give rise to whole cell level physiological functions. Integration of transcriptomic responses and numerical simulations could simulate and predict whole cell functional outputs from transcriptomic data. We developed a computational pipeline that integrates gene expression timelines and systems of coupled differential equations to generate cell-type selective dynamical models. We tested our integrative algorithm on the eicosanoid biosynthesis network in macrophages. Converting transcriptomic changes to a dynamical model allowed us to predict dynamics of prostaglandin and thromboxane synthesis and secretion by macrophages that matched published lipidomics data obtained in the same experiments. Integration of cell-level system biology simulations with genomic and clinical data using a knowledge graph framework will allow us to create explicit predictive models that mechanistically link genomic determinants to organ function. Such integration requires a multi-domain ontological framework to connect genomic determinants to gene expression and cell pathways and functions to organ level phenotypes in healthy and diseased states. These integrated scalable models of tissues and organs as accurate digital twins predict health and disease states for precision medicine.

A spatially specified systems pharmacology therapy for axonal recovery after injury

There are no known drugs or drug combinations that promote substantial central nervous system axonal regeneration after injury. We used systems pharmacology approaches to model pathways underlying axonal growth and identify a four-drug combination that regulates multiple subcellular processes in the cell body and axons using the optic nerve crush model in rats. We intravitreally injected agonists HU-210 (cannabinoid receptor-1) and IL-6 (interleukin 6 receptor) to stimulate retinal ganglion cells for axonal growth. We applied, in gel foam at the site of nerve injury, Taxol to stabilize growing microtubules, and activated protein C to clear the debris field since computational models predicted that this drug combination regulating two subcellular processes at the growth cone produces synergistic growth. Physiologically, drug treatment restored or preserved pattern electroretinograms and some of the animals had detectable visual evoked potentials in the brain and behavioral optokinetic responses. Morphology experiments show that the four-drug combination protects axons or promotes axonal regrowth to the optic chiasm and beyond. We conclude that spatially targeted drug treatment is therapeutically relevant and can restore limited functional recovery.

Synaptojanin 1 Modulates Functional Recovery After Incomplete Spinal Cord Injury in Male Apolipoprotein E Epsilon 4 Mice

Apolipoprotein E epsilon 4 (ApoE4) is the second most common variant of ApoE, being present in ∼14% of the population. Clinical reports identify ApoE4 as a genetic risk factor for poor outcomes after traumatic spinal cord injury (SCI) and spinal cord diseases such as cervical myelopathy. To date, there is no intervention to promote recovery of function after SCI/spinal cord diseases that is specifically targeted at ApoE4-associated impairment. Studies in the human and mouse brain link ApoE4 to elevated levels of synaptojanin 1 (synj1), a lipid phosphatase that degrades phosphoinositol 4,5-bisphosphate (PIP2) into inositol 4-monophosphate. Synj1 regulates rearrangements of the cytoskeleton as well as endocytosis and trafficking of synaptic vesicles. We report here that, as compared to ApoE3 mice, levels of synj1 messenger RNA and protein were elevated in spinal cords of healthy ApoE4 mice associated with lower PIP2 levels. Using a moderate-severity model of contusion SCI in mice, we found that genetic reduction of synj1 improved locomotor function recovery at 14 days after SCI in ApoE4 mice without altering spared white matter. Genetic reduction of synj1 did not alter locomotor recovery of ApoE3 mice after SCI. Bulk RNA sequencing revealed that at 14 days after SCI in ApoE4 mice, genetic reduction of synj1 upregulated genes involved in glutaminergic synaptic transmission just above and below the lesion. Overall, our findings provide evidence for a link between synj1 to poor outcomes after SCI in ApoE4 mice, up to 14 days post-injury, through mechanisms that may involve the function of excitatory glutaminergic neurons.

Whole cell response to receptor stimulation involves many deep and distributed subcellular biochemical processes

Neurite outgrowth is an integrated whole cell response triggered by the cannabinoid-1 receptor. We sought to identify the many different biochemical pathways that contribute to this whole cell response. To understand underlying mechanisms, we identified subcellular processes (SCPs) composed of one or more biochemical pathways and their interactions required for this response. Differentially expressed genes and proteins were obtained from bulk transcriptomics and proteomic analysis of extracts from cells stimulated with a cannabinoid-1 receptor agonist. We used these differentially expressed genes and proteins to build networks of interacting SCPs by combining the expression data with prior pathway knowledge. From these SCP networks, we identified additional genes that when ablated, experimentally validated the SCP involvement in neurite outgrowth. Our experiments and informatics modeling allowed us to identify diverse SCPs such as those involved in pyrimidine metabolism, lipid biosynthesis, and mRNA splicing and stability, along with more predictable SCPs such as membrane vesicle transport and microtubule dynamics. We find that SCPs required for neurite outgrowth are widely distributed among many biochemical pathways required for constitutive cellular functions, several of which are termed ‘deep’, since they are distal to signaling pathways and the key SCPs directly involved in extension of the neurite. In contrast, ‘proximal’ SCPs are involved in microtubule growth and membrane vesicle transport dynamics required for neurite outgrowth. From these bioinformatics and dynamical models based on experimental data, we conclude that receptor-mediated regulation of subcellular functions for neurite outgrowth is both distributed, that is, involves many different biochemical pathways, and deep.

FKBP51 modulates hippocampal size and function in post-translational regulation of Parkin

FK506-binding protein 51 (encoded by Fkpb51, also known as Fkbp5) has been associated with stress-related mental illness. To investigate its function, we studied the morphological consequences of Fkbp51 deletion. Artificial Intelligence-assisted morphological analysis revealed that male Fkbp51 knock-out (KO) mice possess more elongated dentate gyrus (DG) but shorter hippocampal height in coronal sections when compared to WT. Primary cultured Fkbp51 KO hippocampal neurons were shown to exhibit larger dendritic outgrowth than wild-type (WT) controls and pharmacological manipulation experiments suggest that this may occur through the regulation of microtubule-associated protein. Both in vitro primary culture and in vivo labeling support a role for FKBP51 in the regulation of microtubule-associated protein expression. Furthermore, Fkbp51 KO hippocampi exhibited decreases in βIII-tubulin, MAP2, and Tau protein levels, but a greater than 2.5-fold increase in Parkin protein. Overexpression and knock-down FKBP51 demonstrated that FKBP51 negatively regulates Parkin in a dose-dependent and ubiquitin-mediated manner. These results indicate a potential novel post-translational regulatory mechanism of Parkin by FKBP51 and the significance of their interaction on disease onset. KO has more flattened hippocampus using AI-assisted measurement Both pyramidal cell layer (PCL) of CA and granular cell layer (GCL) of DG distinguishable as two layers: deep cell layer and superficial layer. Distinct MAP2 expression between deep and superficial layer between KO and WT, Higher Parkin expression in KO brain Mechanism of FKBP51 inhibition resulting in Parkin, MAP2, Tau, and Tubulin expression differences between KO and WT mice, and resulting neurite outgrowth differences.

RhoA-GTPase Modulates Neurite Outgrowth by Regulating the Expression of Spastin and p60-Katanin

RhoA-GTPase (RhoA) is widely regarded as a key molecular switch to inhibit neurite outgrowth by rigidifying the actin cytoskeleton. However, during neurite outgrowth, whether and how microtubule dynamics are regulated by RhoA remains to be elucidated. Herein, CT04 and Y27632 were used to inactivate RhoA and its downstream effector Rho-associated coiled coil-forming kinase (ROCK), while the RhoAQ63L lentiviral vector was utilized to overexpress the constitutively activated RhoA in dorsal root ganglion (DRG) neurons or neuronal differentiated PC12 cells. The current data illustrate that the RhoA signaling pathway negatively modulates neurite outgrowth and elevates the expression of Glu-tubulin (a marker for a stabilized microtubule). Meanwhile, the microtubule-severing proteins spastin and p60-katanin were downregulated by the RhoA signaling pathway. When spastin and p60-katanin were knocked down, the effects of RhoA inhibition on neurite outgrowth were significantly reversed. Taken together, this study demonstrates that the RhoA pathway-mediated inhibition of neurite outgrowth is not only related to the modulation of microfilament dynamics but is also attributable to the regulation of the expression of spastin and p60-katanin and thus influences microtubule dynamics.