Mathematical modeling of septic shock based on clinical data.

BACKGROUND: Mathematical models of diseases may provide a unified approach for establishing effective treatment strategies based on fundamental pathophysiology. However, models that are useful for clinical practice must overcome the massive complexity of human physiology and the diversity of patients' environmental conditions. With the aim of modeling a complex disease, we choose sepsis, which is highly complex, life-threatening systemic disease with high mortality. In particular, we focused on septic shock, a subset of sepsis in which underlying circulatory and cellular/metabolic abnormalities are profound enough to substantially increase mortality. Our model includes cardiovascular, immune, nervous system models and a pharmacological model as submodels and integrates them to create a sepsis model based on pathological facts. RESULTS: Model validation was done in two steps. First, we established a model for a standard patient in order to confirm the validity of our approach in general aspects. For this, we checked the correspondence between the severity of infection defined in terms of pathogen growth rate and the ease of recovery defined in terms of the intensity of treatment required for recovery. The simulations for a standard patient showed good correspondence. We then applied the same simulations to a patient with heart failure as an underlying disease. The model showed that spontaneous recovery would not occur without treatment, even for a very mild infection. This is consistent with clinical experience. We next validated the model using clinical data of three sepsis patients. The model parameters were tuned for these patients based on the model for the standard patient used in the first part of the validation. In these cases, the simulations agreed well with clinical data. In fact, only a handful parameters need to be tuned for the simulations to match with the data. CONCLUSIONS: We have constructed a model of septic shock and have shown that it can reproduce well the time courses of treatment and disease progression. Tuning of model parameters for each patient could be easily done. This study demonstrates the feasibility of disease models, suggesting the possibility of clinical use in the prediction of disease progression, decisions on the timing of drug dosages, and the estimation of time of infection.

A mathematical model of type 1 diabetes involving leptin effects on glucose metabolism.

Leptin, a hormone released from fat cells in adipose tissues, was recently found to be capable of normalizing glucose metabolism in animals. Clinical data on patients with lipodystrophy indicates that leptin may have a positive effect on glucose metabolism in individuals with diabetes. There are growing expectations that leptin can improve the current insulin treatment for patients with type 1 diabetes. We investigated this possibility through in silico experiments based on a mathematical model of diabetes, which is currently the only mode of research that eliminates human risk. A model of the brain-centered glucoregulatory system, in which leptin plays a central role, was constructed and integrated within a conventional model of insulin/glucose dynamics. The model has been validated using experimental data from animal studies. The in silico combination experiments showed excellent therapeutic performance over insulin monotherapy.

Direct uptake and degradation of DNA by lysosomes.

Lysosomes contain various hydrolases that can degrade proteins, lipids, nucleic acids and carbohydrates. We recently discovered "RNautophagy," an autophagic pathway in which RNA is directly taken up by lysosomes and degraded. A lysosomal membrane protein, LAMP2C, a splice variant of LAMP2, binds to RNA and acts as a receptor for this pathway. In the present study, we show that DNA is also directly taken up by lysosomes and degraded. Like RNautophagy, this autophagic pathway, which we term "DNautophagy," is dependent on ATP. The cytosolic sequence of LAMP2C also directly interacts with DNA, and LAMP2C functions as a receptor for DNautophagy, in addition to RNautophagy. Similarly to RNA, DNA binds to the cytosolic sequences of fly and nematode LAMP orthologs. Together with the findings of our previous study, our present findings suggest that RNautophagy and DNautophagy are evolutionarily conserved systems in Metazoa.

Ubiquitin C-terminal hydrolase L1 (UCH-L1) acts as a novel potentiator of cyclin-dependent kinases to enhance cell proliferation independently of its hydrolase activity.

Dysregulation of cell proliferation and the cell cycle are associated with various diseases, such as cancer. Cyclin-dependent kinases (CDKs) play central roles in cell proliferation and the cell cycle. Ubiquitin C-terminal hydrolase L1 (UCH-L1) is expressed in a restricted range of tissues, including the brain and numerous types of cancer. However, the molecular functions of UCH-L1 remain elusive. In this study, we found that UCH-L1 physically interacts with CDK1, CDK4, and CDK5, enhancing their kinase activity. Using several mutants of UCH-L1, we showed that this enhancement is dependent upon interaction levels between UCH-L1 and CDKs but is independent of the known ubiquitin-related functions of UCH-L1. Gain- and loss-of-function studies revealed that UCH-L1 enhances proliferation of multiple cell types, including human cancer cells. Inhibition of the interaction between UCH-L1 and cell cycle-associated CDK resulted in the abolishment of UCH-L1-induced enhancement of cell proliferation. RNA interference of UCH-L1 reduced the growth of human xenograft tumors in mice. We concluded that UCH-L1 is a novel regulator of the kinase activities of CDKs. We believe that our findings from this study will significantly contribute to our understanding of cell cycle-associated diseases.

Discovery of a novel type of autophagy targeting RNA.

Regulated degradation of cellular components by lysosomes is essential to maintain biological homeostasis. In mammals, three forms of autophagy, macroautophagy, microautophagy and chaperone-mediated autophagy (CMA), have been identified. Here, we showed a novel type of autophagy, in which RNA is taken up directly into lysosomes for degradation. This pathway, which we term "RNautophagy," is ATP-dependent, and unlike CMA, is independent of HSPA8/Hsc70. LAMP2C, a lysosomal membrane protein, serves as a receptor for this pathway. The cytosolic tail of LAMP2C specifically binds to almost all total RNA derived from mouse brain. The cytosolic sequence of LAMP2C and its affinity for RNA are evolutionarily conserved from nematodes to humans. Our findings shed light on the mechanisms underlying RNA homeostasis in higher eukaryotes.

Parkinson's disease-associated mutations in α-synuclein and UCH-L1 inhibit the unconventional secretion of UCH-L1.

Ubiquitin carboxy-terminal hydrolase L1 (UCH-L1) is an intracellular protein abundantly expressed in neurons, and a mutation in UCH-L1 has been identified in familial Parkinson's disease. UCH-L1 has been detected in human cerebrospinal fluid, raising the possibility that UCH-L1 is secreted from neurons. In the present study, we showed that a portion of UCH-L1 is secreted from cultured cells. The secretion of D30K UCH-L1, which lacks ubiquitin binding activity, was decreased compared with that of wild-type UCH-L1, while the secretion of C90S UCH-L1, which lacks hydrolase activity, was not. Treatment with Brefeldin A, an inhibitor of vesicle transport from the endoplasmic reticulum to the Golgi, did not block the secretion of UCH-L1, indicating that UCH-L1 is secreted by an unconventional pathway. The UCH-L1 sequence from Leu-32 to Leu-39 is similar to the unconventional secretory signal sequence of engrailed 2, and substitution of the leucines within this region (L32S/L32A/L34S/L34A/L39S/L39A) reduced the secretion of UCH-L1. We found that the Parkinson's disease-associated mutation I93M in UCH-L1 decreased the secretion of I93M UCH-L1. In addition, Parkinson's disease-linked α-synuclein mutants reduced the secretion of endogenous UCH-L1. Our results indicate that the hydrolase activity is not necessary for the unconventional secretion of UCH-L1, and suggest that the ubiquitin binding activity and the sequence between Leu-32 and Leu-39 are involved in the secretion. Moreover, the secretion of UCH-L1 could be involved in the pathology of Parkinson's disease.

Formulas for intrinsic noise evaluation in oscillatory genetic networks.

The linear noise approximation is a useful method for stochastic noise evaluations in genetic regulatory networks, where the covariance equation described as a Lyapunov equation plays a central role. We discuss the linear noise approximation method for evaluations of an intrinsic noise in autonomously oscillatory genetic networks; in such oscillatory networks, the covariance equation becomes a periodic differential equation that provides generally an unbounded covariance matrix, so that the standard method of noise evaluation based on the covariance matrix cannot be adopted directly. In this paper, we develop a new method of noise evaluation in oscillatory genetic networks; first, we investigate structural properties, e.g., orbital stability and periodicity, of the solutions to the covariance equation given as a periodic Lyapunov differential equation by using the Floquet-Lyapunov theory, and propose a global measure for evaluating stochastic amplitude fluctuations on the periodic trajectory; we also derive an evaluation formula for the period fluctuation. Finally, we apply our method to a model of circadian oscillations based on negative auto-regulation of gene expression, and show validity of our method by comparing the evaluation results with stochastic simulations.

Identification of a novel chemical potentiator and inhibitors of UCH-L1 by in silico drug screening.

Ubiquitin-C-terminal hydrolase L1 (UCH-L1) is a de-ubiquitinating enzyme expressed in the brain and reproductive tissues as well as certain cancers. The hydrolase activity of UCH-L1 has been implicated in Alzheimer's disease and cancer invasion; therefore, it may represent a therapeutic target for these diseases. The present study was undertaken to identify novel chemical modulators for the hydrolase activity of UCH-L1. To identify chemicals that bind to the active site of UCH-L1, we carried out in silico structure-based drug screening using human UCH-L1 crystal structure data (PDB ID: 2ETL) and virtual compound libraries containing 26,891 and 304,205 compounds. Among the compounds with the highest binding scores, we identified one that potentiates the hydrolase activity of UCH-L1, and six that inhibit the activity in enzymatic assays. These compounds may be useful for research on UCH-L1 function, and could lead to candidate therapeutics for UCH-L1-associated diseases.

Aberrant molecular properties shared by familial Parkinson's disease-associated mutant UCH-L1 and carbonyl-modified UCH-L1.

Parkinson's disease (PD) is a neurodegenerative disorder characterized by loss of dopaminergic neurons. The I93M mutation in ubiquitin C-terminal hydrolase L1 (UCH-L1) is associated with familial PD, and we have previously shown that the I93M UCH-L1-transgenic mice exhibit dopaminergic cell loss. Over 90% of neurodegenerative diseases, including PD, occur sporadically. However, the molecular mechanisms underlying sporadic PD as well as PD associated with I93M UCH-L1 are largely unknown. UCH-L1 is abundant (1-5% of total soluble protein) in the brain and is a major target of oxidative/carbonyl damage associated with sporadic PD. As well, abnormal microtubule dynamics and tubulin polymerization are associated with several neurodegenerative diseases including frontotemporal dementia and parkinsonism linked to chromosome 17. Here we show that familial PD-associated mutant UCH-L1 and carbonyl-modified UCH-L1 display shared aberrant properties: compared with wild-type UCH-L1, they exhibit increased insolubility and elevated interactions with multiple proteins, which are characteristics of several neurodegenerative diseases-linked mutants. Circular dichroism analyses suggest similar structural changes in both UCH-L1 variants. We further report that one of the proteins interacting with UCH-L1 is tubulin, and that aberrant interaction of mutant or carbonyl-modified UCH-L1 with tubulin modulates tubulin polymerization. These findings may underlie the toxic gain of function by mutant UCH-L1 in familial PD. Our results also suggest that the carbonyl modification of UCH-L1 and subsequent abnormal interactions of carbonyl-modified UCH-L1 with multiple proteins, including tubulin, constitute one of the causes of sporadic PD.

Auto-regulation of the circadian slave oscillator component AtGRP7 and regulation of its targets is impaired by a single RNA recognition motif point mutation.

The clock-regulated RNA-binding protein AtGRP7 (Arabidopsis thaliana glycine-rich RNA-binding protein) influences circadian oscillations of its transcript by negative feedback at the post-transcriptional level. Here we show that site-specific mutation of one conserved arginine to glutamine within the RNA recognition motif impairs binding of recombinant AtGRP7 to its pre-mRNA in vitro. This correlates with the loss of the negative auto-regulation in vivo: in transgenic plants constitutively overexpressing AtGRP7 (AtGRP7-ox), a shift occurs to an alternatively spliced AtGRP7 transcript that decays rapidly, and thus does not accumulate to high levels. In contrast, constitutive ectopic overexpression of the AtGRP7-RQ mutant does not lead to alternative splicing of the endogenous AtGRP7 transcript and concomitant damping of the oscillations. This highlights the importance of AtGRP7 binding to its own transcript for the negative auto-regulatory circuit. Moreover, regulation of AtGRP7 downstream targets also depends on its RNA-binding activity, as AtGRP8 and other targets identified by transcript profiling of wild-type and AtGRP7-ox plants using fluorescent differential display are negatively affected by AtGRP7 but not by AtGRP7-RQ. In mutants impaired in the nonsense-mediated decay (NMD) components UPF1 or UPF3, levels of the alternatively spliced AtGRP7 and AtGRP8 transcripts that contain premature termination codons are strongly elevated, implicating UPF1 and UPF3 in the decay of these clock-regulated transcripts.

A temperature-compensated model for circadian rhythms that can be entrained by temperature cycles.

From the viewpoint that reaction rates will change with temperature, we present a general method to build circadian clock models that generate circadian oscillations with almost constant period under different constant ambient temperature, and propose an algorithm estimating the parameter condition for compensated period against the change of temperature based on the PER single-feedback loop model of Goldbeter [1995. A model for circadian oscillations in the Drosophila period protein (PER). Proc. R. Soc. London Ser. B 261, 319-324] for Drosophila. We show that the model with derived parameters can realize the temperature compensation over a wide range of temperature, and simultaneously can realize the entrainment to temperature cycles.

Improvement of biomolecule quantification precision and use of a single-element aspheric objective lens in fluorescence correlation spectroscopy.

We found a way to increase the precision with which biomolecules present at concentrations below 10(-10) M can be quantified by fluorescence correlation spectroscopy (FCS). The effectiveness of the way was demonstrated experimentally by using a single-element aspheric objective lens, which was newly developed to reduce the cost of FCS instruments. In the first part of this paper, the relative standard deviation (RSD) of FCS-based concentration measurements is estimated theoretically by an analytical approximation assuming the detection volume profiles in FCS setups to be Gaussian and by molecular simulations in which more realistic profiles are calculated from physical parameters of the measurement setups. In a limit of infinitely bright molecules and zero background emission, the analytical approximation predicts that the RSD at a concentration is minimized when the mean number of molecules in a detection volume is approximately 0.5. A detection volume of the order of 10(-13) L thus gives smaller RSD values for concentrations from 10(-11) to 10(-10) M than does one of the order of 10(-15) L, which is widely used in FCS. This prediction is supported by the molecular simulations, taking into account the finite molecule brightness and background emission. In the second part of the paper, the RSD is evaluated experimentally with an FCS setup with a detection volume of 1.1 x 10(-13) L. The newly developed objective lens, serving as the bottom of the sample cell in this setup, has a large numerical aperture (0.9) without using immersion liquid. When a calibration line was made by 30-s FCS measurements of Cy3-labeled, 112-mer single-stranded DNA solutions, the RSD roughly agreed with the simulation result and was less than 0.1 for DNA concentrations from 2 x 10(-11) to 10(-10) M.

The A and B loci in tobacco regulate a network of stress response genes, few of which are associated with nicotine biosynthesis.

Nicotine biosynthesis in Nicotiana tabacum is under genetic control by the A and B loci. Plants containing semi-dominant mutations at both the A and B loci (i.e. aabb genotype) have lower nicotine levels, reduced nicotine biosynthetic enzyme activities, and reduced mRNA levels of the corresponding biosynthetic genes. The A and B loci therefore appear to be coordinate regulators of several nicotine biosynthetic genes and define a group of co-regulated genes called the A-B regulon. To investigate the composition of genes in the A-B regulon, a fluorescent differential display (FDD) screen was used to randomly sample the transcriptomes of wild type and mutant aabb roots. This resulted in the isolation of 64 FDD clones, representing 49 unique genes or gene families. Four genes associated with nicotine biosynthesis were identified, whereas most of the other FDD clones were homologous with an assortment of stress response genes. Thirty-three genes or gene families showed reproducible aabb genotype effects, representing seven distinct mRNA expression patterns in response media treatments that increase the mRNA levels of known alkaloid biosynthetic genes. Thus, the A and B loci regulate the mRNA levels of some target genes differently than others. Eleven genes or gene families showed only treatment-specific effects, representing four mRNA accumulation patterns. These results indicate the A-B regulon is complex network of differentially expressed stress response genes, only a small subset of which are involved in nicotine biosynthesis.

The identification of candidate genes for a reverse genetic analysis of development and function in the Arabidopsis gynoecium.

The screening for mutants and their subsequent molecular analysis has permitted the identification of a number of genes of Arabidopsis involved in the development and functions of the gynoecium. However, these processes remain far from completely understood. It is clear that in many cases, genetic redundancy and other factors can limit the efficiency of classical mutant screening. We have taken the alternative approach of a reverse genetic analysis of gene function in the Arabidopsis gynoecium. A high-throughput fluorescent differential display screen performed between two Arabidopsis floral homeotic mutants has permitted the identification of a number of genes that are specifically or preferentially expressed in the gynoecium. Here, we present the results of this screen and a detailed characterization of the expression profiles of the genes identified. Our expression analysis makes novel use of several Arabidopsis floral homeotic mutants to provide floral organ-specific gene expression profiles. The results of these studies permit the efficient targeting of effort into a functional analysis of gynoecium-expressed genes.