When the optimal is not the best: parameter estimation in complex biological models

Background

The vast computational resources that became available during the past decade enabled the development and simulation of increasingly complex mathematical models of cancer growth. These models typically involve many free parameters whose determination is a substantial obstacle to model development. Direct measurement of biochemical parameters in vivo is often difficult and sometimes impracticable, while fitting them under data-poor conditions may result in biologically implausible values.

Results

We discuss different methodological approaches to estimate parameters in complex biological models. We make use of the high computational power of the Blue Gene technology to perform an extensive study of the parameter space in a model of avascular tumor growth. We explicitly show that the landscape of the cost function used to optimize the model to the data has a very rugged surface in parameter space. This cost function has many local minima with unrealistic solutions, including the global minimum corresponding to the best fit.

Conclusions

The case studied in this paper shows one example in which model parameters that optimally fit the data are not necessarily the best ones from a biological point of view. To avoid force-fitting a model to a dataset, we propose that the best model parameters should be found by choosing, among suboptimal parameters, those that match criteria other than the ones used to fit the model. We also conclude that the model, data and optimization approach form a new complex system and point to the need of a theory that addresses this problem more generally.     

MicroRNA-210 Regulates Mitochondrial Free Radical Response to Hypoxia and Krebs Cycle in Cancer Cells by Targeting Iron Sulfur Cluster Protein ISCU

Background

Hypoxia in cancers results in the upregulation of hypoxia inducible factor 1 (HIF-1) and a microRNA, hsa-miR-210 (miR-210) which is associated with a poor prognosis.

Methods and Findings

In human cancer cell lines and tumours, we found that miR-210 targets the mitochondrial iron sulfur scaffold protein ISCU, required for assembly of iron-sulfur clusters, cofactors for key enzymes involved in the Krebs cycle, electron transport, and iron metabolism. Down regulation of ISCU was the major cause of induction of reactive oxygen species (ROS) in hypoxia. ISCU suppression reduced mitochondrial complex 1 activity and aconitase activity, caused a shift to glycolysis in normoxia and enhanced cell survival. Cancers with low ISCU had a worse prognosis.

Conclusions

Induction of these major hallmarks of cancer show that a single microRNA, miR-210, mediates a new mechanism of adaptation to hypoxia, by regulating mitochondrial function via iron-sulfur cluster metabolism and free radical generation.     

Genome Rearrangements Detected by SNP Microarrays in Individuals with Intellectual Disability Referred with Possible Williams Syndrome

Background

Intellectual disability (ID) affects 2–3% of the population and may occur with or without multiple congenital anomalies (MCA) or other medical conditions. Established genetic syndromes and visible chromosome abnormalities account for a substantial percentage of ID diagnoses, although for ∼50% the molecular etiology is unknown. Individuals with features suggestive of various syndromes but lacking their associated genetic anomalies pose a formidable clinical challenge. With the advent of microarray techniques, submicroscopic genome alterations not associated with known syndromes are emerging as a significant cause of ID and MCA.

Methodology/Principal Findings

High-density SNP microarrays were used to determine genome wide copy number in 42 individuals: 7 with confirmed alterations in the WS region but atypical clinical phenotypes, 31 with ID and/or MCA, and 4 controls. One individual from the first group had the most telomeric gene in the WS critical region deleted along with 2 Mb of flanking sequence. A second person had the classic WS deletion and a rearrangement on chromosome 5p within the Cri du Chat syndrome (OMIM:123450) region. Six individuals from the ID/MCA group had large rearrangements (3 deletions, 3 duplications), one of whom had a large inversion associated with a deletion that was not detected by the SNP arrays.

Conclusions/Significance

Combining SNP microarray analyses and qPCR allowed us to clone and sequence 21 deletion breakpoints in individuals with atypical deletions in the WS region and/or ID or MCA. Comparison of these breakpoints to databases of genomic variation revealed that 52% occurred in regions harboring structural variants in the general population. For two probands the genomic alterations were flanked by segmental duplications, which frequently mediate recurrent genome rearrangements; these may represent new genomic disorders. While SNP arrays and related technologies can identify potentially pathogenic deletions and duplications, obtaining sequence information from the breakpoints frequently provides additional information.     

2-D protein maps of rat gastrocnemius and soleus muscles: a tool for muscle plasticity assessment

Functional characterization of muscle fibers relies on ATPase activity and on differential measurements of metabolic proteins, including mitochondrial and glycolytic enzymes, glucose, lactate and lactic acid transporters, calcium cycling proteins and components of the contractile machinery. The recent introduction of microarray technology has enabled detailed gene expression studies under different physiological and pathological conditions, thus generating novel hypotheses on muscle function. However, microarray approaches are limited by the incomplete genome coverage of currently available chips, and by poor correlation between mRNA concentration and protein expression level. We have used 2-DE and MS to build a reference map of proteins from rat mixed gastrocnemius and soleus muscle, and to assess qualitative and quantitative differences in protein distribution between these two functionally dissimilar muscles. More than 800 spots on each gel were detected by silver staining, of which 167 were excised, digested in-gel with trypsin and analyzed by ESI-MS/MS. One hundred and twenty eight distinct gene products were identified, including metabolic, transport and contractile proteins. Forty one spots displayed differences in relative expression level between mixed gastrocnemius and soleus samples. These data not only enable differentiation of functionally distinct slow-twitch and fast-twitch fiber types, but also provide tools for investigating muscle plasticity in response to physiological and environmental conditions such as aging or hypoxia.     

The MitoDrome database annotates and compares the OXPHOS nuclear genes of Drosophila melanogaster, Drosophila pseudoobscura and Anopheles gambiae

The oxidative phosphorylation (OXPHOS) is the primary energy-producing process of all aerobic organisms and the only cellular function under the dual control of both the mitochondrial and the nuclear genomes. Functional characterization and evolutionary study of the OXPHOS system is of great importance for the understanding of many as yet unclear aspects of nucleus-mitochondrion genomic co-evolution and co-regulation gene networks. The MitoDrome database is a web-based database which provides genomic annotations about nuclear genes of Drosophila melanogaster encoding for mitochondrial proteins. Recently, MitoDrome has included a new section annotating genomic information about OXPHOS genes in Drosophila pseudoobscura and Anopheles gambiae and their comparative analysis with their Drosophila melanogaster and human counterparts. The introduction of this new comparative annotation section into MitoDrome is expected to be a useful resource for both functional and structural genomics related to the OXPHOS system.     

Down-regulation of the strawberry Bet v 1-homologous allergen in concert with the flavonoid biosynthesis pathway in colorless strawberry mutant

Proteomic screening of strawberry (Fragaria ananassa) yielded a 58% success rate in protein identification in spite of the fact that no genomic sequence is available for this species. This was achieved by a combination of MALDI-MS/MS de novo sequencing of double-derivatized peptides and indel-tolerant searching against local protein databases built on both EST and full-length nucleotide sequences. The amino acid sequence of a strawberry allergen, homologous to the well-known major birch pollen allergen Bet v 1, was partially determined. This strawberry allergen, named Fra a 1 according to the nomenclature for allergen proteins, showed sequence identity of 54 and 77%, respectively, with corresponding allergens from birch and apple. Differential expression, as evaluated by 2-D DIGE, occurred in 10% of protein spots when red strawberries were compared to a colorless (white) strawberry mutant. White strawberries, known to be tolerated by individuals affected by allergy, were found to be virtually free from the strawberry allergen. Also several enzymes in the pathway for biosynthesis of flavonoids, to which the red color pelargonidin belongs, were down-regulated. This approach to assess differential protein expression without access to genomic sequence information can also be applied to other crop plants and phenotypic traits.     

Leaching of Mercury from Carbonated and Non-Carbonated Cement-Solidified Dredged Sediments

This study investigated traditional cement-based and non-conventional (using accelerated carbonation) solidification/stabilization to treat 2 dredged sediments contaminated with mercury from two different locations in UK. Canal and estuarine-derived sediments were mixed with blended binders and powdered activated carbon. Fresh mixtures of sediment and cement were exposed to gaseous carbon dioxide and were allowed to carbonate for fixed time periods, after which they were cured for 28 days. Following curing, samples were leach tested to evaluate the fixation of mercury in the treated products. The results obtained indicated that both conventional and accelerated carbonated treatments were capable of reducing the concentration of mercury in the eluates to acceptable limits.