Identifying quantitative operation principles in metabolic pathways: a systematic method for searching feasible enzyme activity patterns leading to cellular adaptive responses

Background: Optimization methods allow designing changes in a system so that specific goals are attained. These techniques are fundamental for metabolic engineering. However, they are not directly applicable for investigating the evolution of metabolic adaptation to environmental changes. Although biological systems have evolved by natural selection and result in well-adapted systems, we can hardly expect that actual metabolic processes are at the theoretical optimum that could result from an optimization analysis. More likely, natural systems are to be found in a feasible region compatible with global physiological requirements. Results: We first present a new method for globally optimizing nonlinear models of metabolic pathways that are based on the Generalized Mass Action (GMA) representation. The optimization task is posed as a nonconvex nonlinear programming (NLP) problem that is solved by an outer- approximation algorithm. This method relies on solving iteratively reduced NLP slave subproblems and mixed-integer linear programming (MILP) master problems that provide valid upper and lower bounds, respectively, on the global solution to the original NLP. The capabilities of this method are illustrated through its application to the anaerobic fermentation pathway in Saccharomyces cerevisiae. We next introduce a method to identify the feasibility parametric regions that allow a system to meet a set of physiological constraints that can be represented in mathematical terms through algebraic equations. This technique is based on applying the outer-approximation based algorithm iteratively over a reduced search space in order to identify regions that contain feasible solutions to the problem and discard others in which no feasible solution exists. As an example, we characterize the feasible enzyme activity changes that are compatible with an appropriate adaptive response of yeast Saccharomyces cerevisiae to heat shock Conclusion: Our results show the utility of the suggested approach for investigating the evolution of adaptive responses to environmental changes. The proposed method can be used in other important applications such as the evaluation of parameter changes that are compatible with health and disease states.

Weierstrass integrability of differential equations

The integrability problem consists in finding the class of functions a first integral of a given planar polynomial differential system must belong to. We recall the characterization of systems which admit an elementary or Liouvillian first integral. We define {\it Weierstrass integrability} and we determine which Weierstrass integrable systems are Liouvillian integrable. Inside this new class of integrable systems there are non--Liouvillian integrable systems.

Mutations in NDUFB11, Encoding a Complex I Component of the Mitochondrial Respiratory Chain, Cause Microphthalmia with Linear Skin Defects Syndrome.

Microphthalmia with linear skin defects (MLS) syndrome is an X-linked male-lethal disorder also known as MIDAS (microphthalmia, dermal aplasia, and sclerocornea). Additional clinical features include neurological and cardiac abnormalities. MLS syndrome is genetically heterogeneous given that heterozygous mutations in HCCS or COX7B have been identified in MLS-affected females. Both genes encode proteins involved in the structure and function of complexes III and IV, which form the terminal segment of the mitochondrial respiratory chain (MRC). However, not all individuals with MLS syndrome carry a mutation in either HCCS or COX7B. The majority of MLS-affected females have severe skewing of X chromosome inactivation, suggesting that mutations in HCCS, COX7B, and other as-yet-unidentified X-linked gene(s) cause selective loss of cells in which the mutated X chromosome is active. By applying whole-exome sequencing and filtering for X-chromosomal variants, we identified a de novo nonsense mutation in NDUFB11 (Xp11.23) in one female individual and a heterozygous 1-bp deletion in a second individual, her asymptomatic mother, and an affected aborted fetus of the subject's mother. NDUFB11 encodes one of 30 poorly characterized supernumerary subunits of NADH:ubiquinone oxidoreductase, known as complex I (cI), the first and largest enzyme of the MRC. By shRNA-mediated NDUFB11 knockdown in HeLa cells, we demonstrate that NDUFB11 is essential for cI assembly and activity as well as cell growth and survival. These results demonstrate that X-linked genetic defects leading to the complete inactivation of complex I, III, or IV underlie MLS syndrome. Our data reveal an unexpected role of cI dysfunction in a developmental phenotype, further underscoring the existence of a group of mitochondrial diseases associated with neurocutaneous manifestations.

Integrability of a linear center perturbed by a fourth degree homogeneous polynomial

In this work we study the integrability of a two-dimensional autonomous system in the plane with linear part of center type and non-linear part given by homogeneous polynomials of fourth degree. We give sufficient conditions for integrability in polar coordinates. Finally we establish a conjecture about the independence of the two classes of parameters which appear in the system; if this conjecture is true the integrable cases found will be the only possible ones.

Integrability of a linear center perturbed by a fifth degree homogeneous polynomial

In this work we study the integrability of two-dimensional autonomous system in the plane with linear part of center type and non-linear part given by homogeneous polynomials of fifth degree. We give a simple characterisation for the integrable cases in polar coordinates. Finally we formulate a conjecture about the independence of the two classes of parameters which appear on the system; if this conjecture is true the integrable cases found will be the only possible ones.

Steady-state global optimization of metabolic non-linear dynamic models through recasting into power-law canonical models

Background: Design of newly engineered microbial strains for biotechnological purposes would greatly benefit from the development of realistic mathematical models for the processes to be optimized. Such models can then be analyzed and, with the development and application of appropriate optimization techniques, one could identify the modifications that need to be made to the organism in order to achieve the desired biotechnological goal. As appropriate models to perform such an analysis are necessarily non-linear and typically non-convex, finding their global optimum is a challenging task. Canonical modeling techniques, such as Generalized Mass Action (GMA) models based on the power-law formalism, offer a possible solution to this problem because they have a mathematical structure that enables the development of specific algorithms for global optimization. Results: Based on the GMA canonical representation, we have developed in previous works a highly efficient optimization algorithm and a set of related strategies for understanding the evolution of adaptive responses in cellular metabolism. Here, we explore the possibility of recasting kinetic non-linear models into an equivalent GMA model, so that global optimization on the recast GMA model can be performed. With this technique, optimization is greatly facilitated and the results are transposable to the original non-linear problem. This procedure is straightforward for a particular class of non-linear models known as Saturable and Cooperative (SC) models that extend the power-law formalism to deal with saturation and cooperativity. Conclusions: Our results show that recasting non-linear kinetic models into GMA models is indeed an appropriate strategy that helps overcoming some of the numerical difficulties that arise during the global optimization task.

Definition of linear color models in the RGB vector color space to detect red peaches in orchard images taken under natural illumination

This work proposes the detection of red peaches in orchard images based on the definition of different linear color models in the RGB vector color space. The classification and segmentation of the pixels of the image is then performed by comparing the color distance from each pixel to the different previously defined linear color models. The methodology proposed has been tested with images obtained in a real orchard under natural light. The peach variety in the orchard was the paraguayo (Prunus persica var. platycarpa) peach with red skin. The segmentation results showed that the area of the red peaches in the images was detected with an average error of 11.6%; 19.7% in the case of bright illumination; 8.2% in the case of low illumination; 8.6% for occlusion up to 33%; 12.2% in the case of occlusion between 34 and 66%; and 23% for occlusion above 66%. Finally, a methodology was proposed to estimate the diameter of the fruits based on an ellipsoidal fitting. A first diameter was obtained by using all the contour pixels and a second diameter was obtained by rejecting some pixels of the contour. This approach enables a rough estimate of the fruit occlusion percentage range by comparing the two diameter estimates.

Analysis of vineyard differential management zones and relation to vine development, grape maturity and quality

The objective of research was to analyse the potential of Normalized Difference Vegetation Index (NDVI) maps from satellite images, yield maps and grapevine fertility and load variables to delineate zones with different wine grape properties for selective harvesting. Two vineyard blocks located in NE Spain (Cabernet Sauvignon and Syrah) were analysed. The NDVI was computed from a Quickbird-2 multi-spectral image at veraison (July 2005). Yield data was acquired by means of a yield monitor during September 2005. Other variables, such as the number of buds, number of shoots, number of wine grape clusters and weight of 100 berries were sampled in a 10 rows × 5 vines pattern and used as input variables, in combination with the NDVI, to define the clusters as alternative to yield maps. Two days prior to the harvesting, grape samples were taken. The analysed variables were probable alcoholic degree, pH of the juice, total acidity, total phenolics, colour, anthocyanins and tannins. The input variables, alone or in combination, were clustered (2 and 3 Clusters) by using the ISODATA algorithm, and an analysis of variance and a multiple rang test were performed. The results show that the zones derived from the NDVI maps are more effective to differentiate grape maturity and quality variables than the zones derived from the yield maps. The inclusion of other grapevine fertility and load variables did not improve the results.

Signalling by neurotrophins and hepatocyte growth factor regulates axon morphogenesis by differential β-catenin phosphorylation

Tyrosine phosphorylation of ß-catenin, a component of adhesion complexes and the Wnt pathway, affects cell adhesion, migration and gene transcription. By reducing ßcatenin availability using shRNA-mediated gene silencing or expression of intracellular N-cadherin, we show that ß-catenin is required for axon growth downstream of Brain Derived Neurotrophic Factor (BDNF) and Hepatocyte Growth Factor (HGF) signalling. We demonstrate that receptor tyrosine kinases (RTK) Trk and Met interact with and phosphorylate ß-catenin. Neurotrophins (NT) stimulation of Trk receptors results in phosphorylation of ß-catenin at residue Y654 and increased axon growth and branching. Conversely, pharmacological inhibition of Trk or a Y654F mutant blocks these effects. ß-catenin phospho(P)-Y654 colocalizes with the cytoskeleton at growth cones. However, HGF that also increases axon growth and branching, induces ß-catenin phosphorylation at Y142 and a nuclear localization. Interestingly, dominant negative ΔN-TCF4 abolishes the effects of HGF in axon growth and branching, but not of NT. We conclude that NT and HGF signalling differentially phosphorylate ß-catenin, targeting ß-catenin to distinct compartments to regulate axon morphogenesis by TCF4-transcription-dependent and independent mechanisms. These results place ß-catenin downstream of growth factor/RTK signalling in axon differentiation.

THE IMPACT OF ARCHITECTURAL DESIGN ON SOFTWARE DEVELOPMENT

Tässä työssä tutkitaan ohjelmistoarkkitehtuurisuunnitteluominaisuuksien vaikutusta erään client-server –arkkitehtuuriin perustuvan mobiilipalvelusovelluksen suunnittelu- ja toteutusaikaan. Kyseinen tutkimus perustuu reaalielämän projektiin, jonka kvalitatiivinen analyysi paljasti arkkitehtuurikompponenttien välisten kytkentöjen merkittävästi vaikuttavan projektin työmäärään. Työn päätavoite oli kvantitatiivisesti tutkia yllä mainitun havainnon oikeellisuus. Tavoitteen saavuttamiseksi suunniteltiin ohjelmistoarkkitehtuurisuunnittelun mittaristo kuvaamaan kyseisen järjestelmän alijärjestelmien arkkitehtuuria ja luotiin kaksi suunniteltua mittaristoa käyttävää, työmäärää (komponentin suunnittelu-, toteutus- ja testausaikojen summa) arvioivaa mallia, joista toinen on lineaarinen ja toinen epälineaarinen. Näiden mallien kertoimet sovitettiin optimoimalla niiden arvot epälineaarista gloobaalioptimointimenetelmää, differentiaalievoluutioalgoritmia, käyttäen, niin että mallien antamat arvot vastasivat parhaiten mitattua työmäärää sekä kaikilla ominaisuuksilla eli attribuuteilla että vain osalla niistä (yksi jätettiin vuorotellen pois). Kun arkkitehtuurikompenttien väliset kytkennät jätettiin malleista pois, mitattujen ja arvoitujen työmäärien välinen ero (ilmaistuna virheenä) kasvoi eräässä tapauksessa 367 % entisestä tarkoittaen sitä, että näin muodostettu malli vastasi toteutusaikoja huonosti annetulla ainestolla. Tämä oli suurin havaitu virhe kaikkien poisjätettyjen ominaisuuksien kesken. Saadun tuloksen perusteella päätettiin, että kyseisen järjestelmän toteutusajat ovat vahvasti riippuvaisia kytkentöjen määrästä, ja näin ollen kytkentöjen määrä oli mitä todennäköisemmin kaikista tärkein työmäärään vaikuttava tekijä tutkitun järjestelmän arkkitehtuurisuunnittelussa.

Eosinophilic fasciitis: Typical abnormalities, variants and differential diagnosis of fasciae abnormalities using MR imaging.

Eosinophilic fasciitis is a rare condition. It is generally limited to the distal parts of the arms and legs. MRI is the ideal imaging modality for diagnosing and monitoring this condition. MRI findings typically evidence only fascial involvement but on a less regular basis signal abnormalities may be observed in neighboring muscle tissue and hypodermic fat. Differential diagnosis of eosinophilic fasciitis by MRI requires the exclusion of several other superficial and deep soft tissue disorders.

Differential monocular vs. binocular pupil responses from melanopsin-based photoreception in patients with anterior ischemic optic neuropathy.

We examined the effect of anterior ischemic optic neuropathy (AION) on the activity of intrinsically photosensitive retinal ganglion cells (ipRGCs) using the pupil as proxy. Eighteen patients with AION (10 unilateral, 8 bilateral) and 29 age-matched control subjects underwent chromatic pupillometry. Red and blue light stimuli increasing in 0.5 log steps were presented to each eye independently under conditions of dark and light adaptation. The recorded pupil contraction was plotted against stimulus intensity to generate scotopic and photopic response curves for assessment of synaptically-mediated ipRGC activity. Bright blue light stimuli presented monocularly and binocularly were used for melanopsin activation. The post-stimulus pupil size (PSPS) at the 6th second following stimulus offset was the marker of intrinsic ipRGC activity. Finally, questionnaires were administered to assess the influence of ipRGCs on sleep. The pupil response and PSPS to all monocularly-presented light stimuli were impaired in AION eyes, indicating ipRGC dysfunction. To binocular light stimulation, the PSPS of AION patients was similar to that of controls. There was no difference in the sleep habits of the two groups. Thus after ischemic injury to one or both optic nerves, the summated intrinsic ipRGC activity is preserved when both eyes receive adequate light exposure.

Langevin equations with multiplicative noise: Application to domain growth

Langevin Equations of Ginzburg-Landau form, with multiplicative noise, are proposed to study the effects of fluctuations in domain growth. These equations are derived from a coarse-grained methodology. The Cahn-Hiliard-Cook linear stability analysis predicts some effects in the transitory regime. We also derive numerical algorithms for the computer simulation of these equations. The numerical results corroborate the analytical predictions of the linear analysis. We also present simulation results for spinodal decomposition at large times.

Conjugation of a Ru(II) Arene Complex to Neomycin or to Guanidinoneomycin Leads to Compounds with Differential Cytotoxicities and Accumulation between Cancer and Normal Cells

A straightforward methodology for the synthesis of conjugates between a cytotoxic organometallic ruthenium(II) complex and amino- and guanidinoglycosides, as potential RNA-targeted anticancer compounds, is described. Under microwave irradiation, the imidazole ligand incorporated on the aminoglycoside moiety (neamine or neomycin) was found to replace one triphenylphosphine ligand from the ruthenium precursor [(η6-p-cym)RuCl(PPh3)2]+, allowing the assembly of the target conjugates. The guanidinylated analogue was easily prepared from the neomycin-ruthenium conjugate by reaction with N,N′-di-Boc-N″-triflylguanidine, a powerful guanidinylating reagent that was compatible with the integrity of the metal complex. All conjugates were purified by semipreparative high-performance liquid chromatography (HPLC) and characterized by electrospray ionization (ESI) and matrix-assisted laser desorption-ionization time-of-flight (MALDI-TOF) mass spectrometry (MS) and NMR spectroscopy. The cytotoxicity of the compounds was tested in MCF-7 (breast) and DU-145 (prostate) human cancer cells, as well as in the normal HEK293 (Human Embryonic Kidney) cell line, revealing a dependence on the nature of the glycoside moiety and the type of cell (cancer or healthy). Indeed, the neomycin-ruthenium conjugate (2) displayed moderate antiproliferative activity in both cancer cell lines (IC50 ≈ 80 μM), whereas the neamine conjugate (4) was inactive (IC50 ≈ 200 μM). However, the guanidinylated analogue of the neomycin-ruthenium conjugate (3) required much lower concentrations than the parent conjugate for equal effect (IC50 = 7.17 μM in DU-145 and IC50 = 11.33 μM in MCF-7). Although the same ranking in antiproliferative activity was found in the nontumorigenic cell line (3 2 > 4), IC50 values indicate that aminoglycoside-containing conjugates are about 2-fold more cytotoxic in normal cells (e.g., IC50 = 49.4 μM for 2) than in cancer cells, whereas an opposite tendency was found with the guanidinylated conjugate, since its cytotoxicity in the normal cell line (IC50 = 12.75 μM for 3) was similar or even lower than that found in MCF-7 and DU-145 cancer cell lines, respectively. Cell uptake studies performed by ICP-MS with conjugates 2 and 3 revealed that guanidinylation of the neomycin moiety had a positive effect on accumulation (about 3-fold higher in DU-145 and 4-fold higher in HEK293), which correlates well with the higher antiproliferative activity of 3. Interestingly, despite the slightly higher accumulation in the normal cell than in the cancer cell line (about 1.4-fold), guanidinoneomycin-ruthenium conjugate (3) was more cytotoxic to cancer cells (about 1.8-fold), whereas the opposite tendency applied for neomycin-ruthenium conjugate (2). Such differences in cytotoxic activity and cellular accumulation between cancer and normal cells open the way to the creation of more selective, less toxic anticancer metallodrugs by conjugating cytotoxic metal-based complexes such as ruthenium(II) arene derivatives to guanidinoglycosides.

Absoluuttisten solmukoordinaattien menetelmä kuorimaisten kappaleiden mallinnuksessa

Työn tavoitteena oli kehittää nopeasti konvergoiva kuorielementti epälineaarisesti joustavien kappaleiden analysointiin. Kuorielementti perustuu absoluuttisten solmukoordinaattien menetelmään ja se hyödyntää kaarevuuden kuvausta elastisten voimien määrityksessä. Kehitettyä elementtiä verrattiin kontinuumimekaniikalla kehitettyyn kuorielementtiin ja kaupallisen elementtimenetelmän kuorielementtiin. Yksinkertaisimman kuormitustapauksen tuloksia verrattiin teknisen taivutusteorian mukaiseen analyyttiseen ratkaisuun. Staattisten testien tulokset tässä työssä kehitetyllä kuorielementillä vastasivat hyvin kaupallisella elementtimenetelmällä saatuja tuloksia. Deformaatioiden ollessa geometrisesti lineaarisella alueella, kehitetyllä kuorielementillä saadut tulokset vastasivat paremmin sekä analyyttistä ratkaisua että kaupallisella elementtimenetelmällä saatuja tuloksia kuin aiemman kontinuumimekaniikkaan perustuvan kuorielementin tulokset. Kehitetyn kuorielementin ongelmana verrattuna kontinuumimekaniikkaan perustuvaan elementtiin on monimutkaisempi kinematiikan kuvaus. Tästä on seurauksena laskenta-ajan huomattava kasvaminen. Jatkossa kannattaisi keskittyä numeeristen ratkaisumenetelmien kehittämiseen.