Practical limits for reverse engineering of dynamical systems: a statistical analysis of sensitivity and parameter inferability in systems biology models

The size and complexity of cellular systems make building predictive models an extremely difficult task. In principle dynamical time-course data can be used to elucidate the structure of the underlying molecular mechanisms, but a central and recurring problem is that many and very different models can be fitted to experimental data, especially when the latter are limited and subject to noise. Even given a model, estimating its parameters remains challenging in real-world systems. Here we present a comprehensive analysis of 180 systems biology models, which allows us to classify the parameters with respect to their contribution to the overall dynamical behaviour of the different systems. Our results reveal candidate elements of control in biochemical pathways that differentially contribute to dynamics. We introduce sensitivity profiles that concisely characterize parameter sensitivity and demonstrate how this can be connected to variability in data. Systematically linking data and model sloppiness allows us to extract features of dynamical systems that determine how well parameters can be estimated from time-course measurements, and associates the extent of data required for parameter inference with the model structure, and also with the global dynamical state of the system. The comprehensive analysis of so many systems biology models reaffirms the inability to estimate precisely most model or kinetic parameters as a generic feature of dynamical systems, and provides safe guidelines for performing better inferences and model predictions in the context of reverse engineering of mathematical models for biological systems.     

Demography and dispersal ability of a threatened saproxylic beetle: a mark-recapture study of the Rosalia Longicorn (Rosalia alpina)

The Rosalia longicorn or Alpine longhorn (Coleoptera: Cerambycidae) is an endangered and strictly protected icon of European saproxylic biodiversity. Despite its popularity, lack of information on its demography and mobility may compromise adoption of suitable conservation strategies. The beetle experienced marked retreat from NW part of its range; its single population survives N of the Alps and W of the Carpathians. The population inhabits several small patches of old beech forest on hill-tops of the Ralska Upland, Czech Republic. We performed mark-recapture study of the population and assessed its distribution pattern. Our results demonstrate the high mobility of the beetle, including dispersal between hills (up to 1.6 km). The system is thus interconnected; it contained ～2000 adult beetles in 2008. Estimated population densities were high, ranging between 42 and 84 adult beetles/hectare a year. The population survives at a former military-training ground despite long-term isolation and low cover of mature beech forest (～1%). Its survival could be attributed to lack of forestry activities between the 1950s and 1990s, slow succession preventing canopy closure and undergrowth expansion, and probably also to the distribution of habitat patches on conspicuous hill-tops. In order to increase chances of the population for long term survival, we propose to stop clear-cuts of old beech forests, increase semi-open beech woodlands in areas currently covered by conifer plantations and active habitat management at inhabited sites and their wider environs.     

A comparison of the therapeutic and reactivating efficacy of newly developed bispyridinium compounds (K206, K269) with currently available oximes against tabun in rats and mice

The potency of newly developed bispyridinium compounds (K206, K269) in reactivating tabun-inhibited acetylcholinesterase and eliminating tabun-induced lethal toxic effects was compared with commonly used oximes (obidoxime, trimedoxime, the oxime HI-6) using in vivo methods. Studies which determined percentage of reactivation of tabun-inhibited blood and tissue AChE in poisoned rats showed that the reactivating efficacy of both newly developed oximes is comparable with obidoxime and trimedoxime in blood but lower than the reactivating potency of trimedoxime and obidoxime in the diaphragm and brain. Nevertheless, the differences in reactivating efficacy of obidoxime, trimedoxime and K206 was not significant while the potency of K269 to reactivate tabun-inhibited acetylcholinesterase was significantly lower. Both newly developed oximes were also found to be relatively efficacious in elimination of the lethal toxic effects in tabun-poisoned mice. Their therapeutic efficacy corresponds to the therapeutic potency of obidoxime. The oxime HI-6, relatively efficacious against soman, did not seem to be an adequately effective oxime in reactivation of tabun-inhibited AChE and to counteract lethal effects of tabun. Both newly developed oximes (K206, K269) are significantly more efficacious in reactivating tabun-inhibited AChE in rats and to eliminate lethal toxic effects of tabun in mice than the oxime HI-6 but their reactivating and therapeutic potency does not prevail over the effectiveness of currently available obidoxime and trimedoxime and, therefore, they are not suitable for their replacement of commonly used oximes for the treatment of acute tabun poisoning.     

Discovery of pyrazolo[1,5-a]pyrimidine-based CHK1 inhibitors: a template-based approach--part 1

The synthesis and hit-to-lead SAR development of a pyrazolo[1,5-a]pyrimidine hit 4 is described leading to a series of potent, selective CHK1 inhibitors such as compound 17r. In the Letter, the further utility of the pyrazolo[1,5-a]pyrimidine template for the development of potent, selective kinase inhibitors is detailed.     

Synthesis and characterization of a cell-permeable bimodal contrast agent targeting β-galactosidase

Noninvasive monitoring of intracellular targets such as enzymes, receptors, or mRNA by means of magnetic resonance imaging (MRI) is increasingly gaining relevance in various research areas. A vital prerequisite for their visualization is the development of cell-permeable imaging probes, which can specifically interact with the target that characterizes the cellular or molecular process of interest. Here, we describe a dual-labeled probe, Gd-DOTA-k(FR)-Gal-CPP, designed to report the presence of intracellular β-galactosidase (β-gal) enzyme by MRI. This conjugate consists of a galactose based core serving as cleavable spacer, incorporated between the cell-penetrating peptide D-Tat(49-57) and reporter moieties (Gd-DOTA, fluorescein (FR)). We employed a facile building block approach to obtain our bimodal probe, Gd-DOTA-k(FR)-Gal-CPP. This strategy involved the preparation of the building blocks and their subsequent assembly using Fmoc-mediated solid phase synthesis, followed by the complexation of ligand 14 with GdCl(3). Gd-DOTA-k(FR)-Gal-CPP showed a considerably higher relaxivity enhancement (16.8±0.6 mM(-1)s(-1), 123 MHz, ～21°C) relative to the commercial Gd-DOTA (4.0±0.12 mM(-1)s(-1), 123MHz, ～21 °C). The activation of Gd-DOTA-k(FR)-Gal-CPP was based on a cellular retention strategy that required enzymatic cleavage of the delivery vector from galactose moiety following the cell internalization to achieve a prolonged accumulation of the reporter components (Gd-DOTA/FR) in the β-gal expressing cells. Cellular uptake of Gd-DOTA-k(FR)-Gal-CPP in β-gal expressing C6/LacZ and enzyme deficient parental C6 rat glioma cells was confirmed by fluorescence spectroscopy, MR imaging and ICP-AES measurements. All methods showed higher accumulation of measured reporters in C6/LacZ cells compared to enzyme deficient parental C6 cells. Fluorescence microscopy of cells labeled with Gd-DOTA-k(FR)-Gal-CPP indicated a predominantly vesicular localization of the green fluorescent conjugate around cell nuclei. This cellular distribution was most likely responsible for the observed non-specific background signal in the enzyme deficient C6 cells. Even though the specific accumulation of our bimodal probe has to be further improved, it could be already used for cell imaging by MRI and optical modalities.     

The open-access dataset for insilico cardiotoxicity prediction system

Drug cardiotoxicity is one of the main reasons of fatal drug related problem events and the subsequent withdrawals. Therefore, its early assessment is a crucial element of the drug development process. For the drug driven hERG inhibition assessment, which is assumed to be the main reason for toxicity, in vitro techniques are used. Gold standards are based on the Patch Clamp method with the use of various cell models but due to its low throughput, insilico models have become more appreciated. To develop a reliable empirical QSAR model, wide dataset containing a variety of cases has to be available. In this article, a freely available for download, set of data is described. It is based on literature peer-reviewed reports and contains hERG inhibition information expressed as IC50 value for 263 molecules described in 642 records. All studies were done with the use of three cell models (XO, CHO, HEK) and other elements describe the electrophysiological settings of the in vitro study. The above mentioned set was used for the successful development of the predictive models.     

The changes in pectin metabolism in flax infected with Fusarium

Fusarium culmorum and Fusarium oxysporum are the most common fungal pathogens of flax (Linum usitatissimum L.), thus leading to the greatest losses in crop yield. A subtractive cDNA library was constructed from flax seedlings exposed for two days to F. oxysporum. This revealed a set of genes that are potentially involved in the flax defense responses. Two of those genes directly participate in cell wall sugar polymer metabolism: UDP-d-glucuronate 4-epimerase (GAE; EC 5.1.3.6) and formate dehydrogenase (FDH; EC 1.2.1.2). GAE delivers the main substrate for pectin biosynthesis, and decreases were detected in its mRNA level after Fusarium infection. FDH participates in the metabolism of formic acid, and the expression level of its gene increased after Fusarium infection. However, metabolite profiling analysis disclosed that the pectin content in the infected plants remained unchanged, but that there were reductions in both the levels of the soluble sugars that serve as pectin precursors, and in the level of formic acid. Since formic acid is the product of pectin demethylesterification, the level of mRNAs coding for pectin methylesterase (EC 3.1.1.11) in the infected flax was measured, revealing a decrease in its expression upon plant infection. Transgenic flax plants overexpressing fungal polygalacturonase (EC 3.2.1.15) and rhamnogalacturonase (EC 3.2.1.-) showed a decrease in the pectin content and an elevated level of formic acid, but the level of expression of the FDH gene remained unchanged. It is suspected that the expression of the formate dehydrogenase gene is directly controlled by the pathogen in the early stage of infection, and additionally by pectin degradation in the later stages.     

Polymorphisms of stress-related genes and the risk of nonsyndromic cleft lip with or without cleft palate

BACKGROUND: Nonsyndromic cleft lip with or without cleft palate (NCL/P) is a common structural malformation with a complex and multifactorial etiology. It has been shown that maternal psychological stress in the periconceptional period can contribute to an increase in the risk of NCL/P affecting pregnancy. METHODS: Twenty‐four single nucleotide polymorphisms of 11 stress‐related genes (COMT, CRHR1, FKBP5, GABRA6, HSD11β2, MAOA, NPY, NR3C1, SERPINA6, SLC6A4, and TPH2) were investigated in 220 healthy mothers of children with facial clefts and 210 matched controls using restriction fragment‐length polymorphism and high‐resolution melting analysis. RESULTS: We found that polymorphisms in SLC6A4, TPH2, and SERPINA6 appear to be maternal factors increasing the risk of having a child with facial clefts. The closest correlations with NCL/P were found for the SLC6A4 rs2020942 and TPH2 rs10879357 gene variants (odds ratio [OR], 1.720; 95% confidence interval [CI], 1.158–2.553; p = 0.0069; p_trend = 0.0036; and OR, 1.837; 95% CI, 1.226–2.753, p = 0.0030, p_trend = 0.0057; respectively). Moreover, haplotype analysis revealed that several combinations of markers in SLC6A4, TPH2, and SERPINA6 might be significantly associated with the risk of NCL/P affected pregnancies. However, these associations were not statistically significant after correction for multiple testing. CONCLUSION: This study suggests that nucleotide variants of genes encoding components of the hypothalamus‐pituitary‐adrenal axis and serotoninergic system have a role in the etiology of NCL/P in the Polish population. SLC6A4, TPH2, and SERPINA6 might be novel candidate genes for this common congenital anomaly. Birth Defects Research (Part A), 2011. © 2011 Wiley‐Liss, Inc.     

MHC influences infection with parasites and winter survival in the root vole Microtus oeconomus

Selective pressure from parasites is thought to maintain the polymorphism of major histocompatibility complex (MHC) genes. Although a number of studies have shown a relationship between the MHC and parasitic infections, the fitness consequences of such associations are less well documented. In the present paper, we characterised the variation in exon 2 of MHC class II DRB gene in the root vole and examined the effects of that gene on parasite prevalence and winter survival. We identified 18 unique exon 2 sequences, which translated into 10 unique amino acid sequences. Phylogenetic analysis revealed the presence of three distinct clusters, and allele distributions among these individuals suggested that the clusters correspond to three different loci. Although the rate of synonymous substitutions (d_S) exceeded the rate of nonsynonymous substitutions (d_N) across sequences, implying purifying selection, d_N was significantly elevated at antigen-binding sites, suggesting that these sites could be under positive selection. Screening for parasites revealed a moderate prevalence of infection with gastrointestinal parasites (24 % infected), but a high infection rate for blood parasites (56 % infected). Infection with the blood parasite Babesia ssp. decreased survival almost twofold (25.7 vs. 13.9 %). Animals possessing the amino acid sequence AA*08 survived better than others (44.9 vs. 22 %), and they were infected with Babesia ssp. less often (13.9 vs 25.7 %). In contrast, individuals carrying allele AA*05 were infected more often (31.7 vs. 15.3 %). Heterozygosity at one of the putative loci was associated with a lower probability of infection with Babesia ssp., but at the other locus, the association was reversed. The unexpected latter result could be at least partly explained by the increased frequency of the susceptible allele AA*05 among heterozygotes. Overall, we demonstrate that infection with Babesia ssp. is a strong predictor of winter survival and that MHC genes are important predictors of infection status as well as survival in the root vole.