Peptaibolin analogues by incorporation of α,α-dialkylglycines: synthesis and study of their membrane permeating ability

Analogues of Peptaibolin, a peptaibol with antibiotic activity, incorporating α,α-dialkylglycines (Deg, Dpg, and Ac6c) at selected positions were synthesised by MW-SPPS and fully characterized. A control analogue incorporating L-alanine was also prepared. The native peptide and the analogues were studied by fluorescence spectroscopy for their membrane permeating activity. Small unilamellar vesicles (SUVs) of egg phosphatidylcholine/ cholesterol (70:30) containing an encapsulated fluorescence probe (6-carboxyfluorescein) were used as membrane models. The assays of carboxyfluorescein release from SUVs upon peptide addition showed that Peptaibolin-Dpg and Peptaibolin-Ac6c are the most active peptides. These results indicate that the structure of the α,α-dialkylglycines is crucial for the membrane permeating ability of these Peptaibolin analogues.

The ghost of past species occurrence: improving species distribution models for presence-only data

1. Model-based approaches have been used increasingly in conservation biology over recent years. Species presence data used for predictive species distribution modelling are abundant in natural history collections, whereas reliable absence data are sparse, most notably for vagrant species such as butterflies and snakes. As predictive methods such as generalized linear models (GLM) require absence data, various strategies have been proposed to select pseudo-absence data. However, only a few studies exist that compare different approaches to generating these pseudo-absence data. 2. Natural history collection data are usually available for long periods of time (decades or even centuries), thus allowing historical considerations. However, this historical dimension has rarely been assessed in studies of species distribution, although there is great potential for understanding current patterns, i.e. the past is the key to the present. 3. We used GLM to model the distributions of three 'target' butterfly species, Melitaea didyma, Coenonympha tullia and Maculinea teleius, in Switzerland. We developed and compared four strategies for defining pools of pseudo-absence data and applied them to natural history collection data from the last 10, 30 and 100 years. Pools included: (i) sites without target species records; (ii) sites where butterfly species other than the target species were present; (iii) sites without butterfly species but with habitat characteristics similar to those required by the target species; and (iv) a combination of the second and third strategies. Models were evaluated and compared by the total deviance explained, the maximized Kappa and the area under the curve (AUC). 4. Among the four strategies, model performance was best for strategy 3. Contrary to expectations, strategy 2 resulted in even lower model performance compared with models with pseudo-absence data simulated totally at random (strategy 1). 5. Independent of the strategy model, performance was enhanced when sites with historical species presence data were not considered as pseudo-absence data. Therefore, the combination of strategy 3 with species records from the last 100 years achieved the highest model performance. 6. Synthesis and applications. The protection of suitable habitat for species survival or reintroduction in rapidly changing landscapes is a high priority among conservationists. Model-based approaches offer planning authorities the possibility of delimiting priority areas for species detection or habitat protection. The performance of these models can be enhanced by fitting them with pseudo-absence data relying on large archives of natural history collection species presence data rather than using randomly sampled pseudo-absence data.

Solid-phase synthesis and NMR studies of modified oligonucleotides forming triplex helix and of oligonucleopeptides mimicking the topoisomerase I-DNA covalent complex

Report for the scientific sojourn carried out at the Institut de Biologia Molecular de Barcelona of the CSIC –state agency – from april until september 2007. Topoisomerase I is an essential nuclear enzyme that modulates the topological status of DNA, facilitating DNA helix unwinding during replication and transcription. We have prepared the oligonucleotide-peptide conjugate Ac-NLeu-Asn-Tyr(p-3’TTCAGAAGC5’)-LeuC-CONH-(CH2)6-OH as model compound for NMR studies of the Topoisomerase I- DNA complex. Special attention was made on the synthetic aspects for the preparation of this challenging compound especially solid supports and protecting groups. The desired peptide was obtained although we did not achieve the amount of the conjugate needed for NMR studies. Most probably the low yield is due to the intrinsic sensitive to hydrolysis of the phosphate bond between oligonucleotide and tyrosine. We have started the synthesis and the structural characterization of oligonucleotides carrying intercalating compounds. At the present state we have obtained model duplex and quadruplex sequences modified with acridine and NMR studies are underway. In addition to this project we have successfully resolved the structure of a fusion peptide derived from hepatitis C virus envelope synthesized by the group of Dr. Haro and we have synthesized and started the characterization of a modified G-quadruplex.

Spatial and temporal dynamics of jasmonate synthesis and accumulation in Arabidopsis in response to wounding.

A new metabolite profiling approach combined with an ultrarapid sample preparation procedure was used to study the temporal and spatial dynamics of the wound-induced accumulation of jasmonic acid (JA) and its oxygenated derivatives in Arabidopsis thaliana. In addition to well known jasmonates, including hydroxyjasmonates (HOJAs), jasmonoyl-isoleucine (JA-Ile), and its 12-hydroxy derivative (12-HOJA-Ile), a new wound-induced dicarboxyjasmonate, 12-carboxyjasmonoyl-l-isoleucine (12-HOOCJA-Ile) was discovered. HOJAs and 12-HOOCJA-Ile were enriched in the midveins of wounded leaves, strongly differentiating them from the other jasmonate metabolites studied. The polarity of these oxylipins at physiological pH correlated with their appearance in midveins. When the time points of accumulation of different jasmonates were determined, JA levels were found to increase within 2-5 min of wounding. Remarkably, these changes occurred throughout the plant and were not restricted to wounded leaves. The speed of the stimulus leading to JA accumulation in leaves distal to a wound is at least 3 cm/min. The data give new insights into the spatial and temporal accumulation of jasmonates and have implications in the understanding of long-distance wound signaling in plants.

Synthesis, binding affinity and SAR of new benzolactam derivatives as dopamine D3 receptor ligands

A series of new benzolactam derivatives was synthesized and the derivatives were evaluated for theiraffinities at the dopamine D1, D2, and D3 receptors. Some of these compounds showed high D2 and/orD3 affinity and selectivity over the D1 receptor. The SAR study of these compounds revealed structuralcharacteristics that decisively influenced their D2 and D3 affinities. Structural models of the complexesbetween some of the most representative compounds of this series and the D2 and D3 receptors wereobtained with the aim of rationalizing the observed experimental results. Moreover, selected compoundsshowed moderate binding affinity on 5-HT2A which could contribute to reducing the occurrence of extrapyramidalside effects as potential antipsychotics.

Lipid metabolism in myelinating glial cells: lessons from human inherited disorders and mouse models.

The integrity of central and peripheral nervous system myelin is affected in numerous lipid metabolism disorders. This vulnerability was so far mostly attributed to the extraordinarily high level of lipid synthesis that is required for the formation of myelin, and to the relative autonomy in lipid synthesis of myelinating glial cells because of blood barriers shielding the nervous system from circulating lipids. Recent insights from analysis of inherited lipid disorders, especially those with prevailing lipid depletion and from mouse models with glia-specific disruption of lipid metabolism, shed new light on this issue. The particular lipid composition of myelin, the transport of lipid-associated myelin proteins, and the necessity for timely assembly of the myelin sheath all contribute to the observed vulnerability of myelin to perturbed lipid metabolism. Furthermore, the uptake of external lipids may also play a role in the formation of myelin membranes. In addition to an improved understanding of basic myelin biology, these data provide a foundation for future therapeutic interventions aiming at preserving glial cell integrity in metabolic disorders.

Do habitat suitability models reliably predict the recovery areas of threatened species?

1. Identifying those areas suitable for recolonization by threatened species is essential to support efficient conservation policies. Habitat suitability models (HSM) predict species' potential distributions, but the quality of their predictions should be carefully assessed when the species-environment equilibrium assumption is violated.2. We studied the Eurasian otter Lutra lutra, whose numbers are recovering in southern Italy. To produce widely applicable results, we chose standard HSM procedures and looked for the models' capacities in predicting the suitability of a recolonization area. We used two fieldwork datasets: presence-only data, used in the Ecological Niche Factor Analyses (ENFA), and presence-absence data, used in a Generalized Linear Model (GLM). In addition to cross-validation, we independently evaluated the models with data from a recolonization event, providing presences on a previously unoccupied river.3. Three of the models successfully predicted the suitability of the recolonization area, but the GLM built with data before the recolonization disagreed with these predictions, missing the recolonized river's suitability and badly describing the otter's niche. Our results highlighted three points of relevance to modelling practices: (1) absences may prevent the models from correctly identifying areas suitable for a species spread; (2) the selection of variables may lead to randomness in the predictions; and (3) the Area Under Curve (AUC), a commonly used validation index, was not well suited to the evaluation of model quality, whereas the Boyce Index (CBI), based on presence data only, better highlighted the models' fit to the recolonization observations.4. For species with unstable spatial distributions, presence-only models may work better than presence-absence methods in making reliable predictions of suitable areas for expansion. An iterative modelling process, using new occurrences from each step of the species spread, may also help in progressively reducing errors.5. Synthesis and applications. Conservation plans depend on reliable models of the species' suitable habitats. In non-equilibrium situations, such as the case for threatened or invasive species, models could be affected negatively by the inclusion of absence data when predicting the areas of potential expansion. Presence-only methods will here provide a better basis for productive conservation management practices.

Insight in protein S deficiency from mouse models

Protein S (ProS) is an important negative regulator of blood coagulation. Its physiological importance is evident in purpura fulminans and other life-threatening thrombotic disorders typical of ProS deficient patients. Our previous characterization of ProS deficiency in mouse models has shown similarities with the human phenotypes: heterozygous ProS-deficient mice (Pros+/-) had increased thrombotic risk whereas homozygous deficiency in ProS (Pros-/-) was incompatible with life (Blood 2009; 114:2307-2314). In tissues, ProS exerts cellular functions by binding to and activating tyrosine kinase receptors of the Tyro3 family (TAM) on the cell surface.To extend the analysis of coagulation defects beyond the Pros-/- phenotype and add new insights into the sites of synthesis ProS and its action, we generated mice with inactivated ProS in hepatocytes (Proslox/loxAlbCre+) as well as in endothelial and hematopoietic cells (Proslox/loxTie2Cre+). Both models resulted in significant reduction of circulating ProS levels and in a remarkable increased thrombotic risk in vivo. In a model of tissue factor (TF)-induced venous thromboembolism (VTE), only 17% of Proslox/loxAlbCre+ mice (n=12) and only 13% of Proslox/loxTie2Cre+ mice (n=14) survived, compared with 86% of Proslox/lox mice (n=14; P<0.001).To mimic a severe acquired ProS deficiency, ProS gene was inactivated at the adult stage using the polyI:C-inducible Mx1-Cre system (Proslox/loxMx1Cre+). Ten days after polyI:C treatment, Proslox/loxMx1Cre+ mice developed disseminated intravascular coagulation with extensive lung and liver thrombosis.It is worth noting that no skin lesions compatible with purpura fulminans were observed in any of the above-described models of partial ProS deficiency. In order to shed light on the pathogenesis of purpura fulminans, we exposed the different ProS-deficient mice to warfarin (0.2 mg/day). We observed that Pros+/-, Proslox/loxAlbCre+ and Proslox/loxTie2Cre+ mice developed retiform purpura (characterized by erythematous and necrotic lesions of the genital region and extremities) and died after 3 to 5 days after the first warfarin administration.In human, ProS is also synthesized by megakaryocytes and hence stored at high concentrations in circulating platelets (pProS). The role of pProS has been investigated by generating megakaryocyte ProS-deficient model using the PF4 promoter as Cre driver (Proslox/loxPf4Cre+). In the TF-induced VTE model, Proslox/loxPf4Cre+ (n=15) mice showed a significant increased risk of thrombosis compared to Proslox/lox controls (n=14; survival rate 47% and 86%, respectively; P<0.05). Furthermore, preliminary results suggest survival to be associated with higher circulating ProS levels. In order to evaluate the potential role of pProS in thrombus formation, we investigated the thrombotic response to intravenous injection of collagen-epinephrine in vivo and platelet function in vitro. Both in vivo and in vitro experiments showed similar results between Proslox/loxPf4Cre+ and Proslox/lox, indicating that platelet reactivity was not influenced by the absence of pProS. These data suggest that pProS is delivered at the site of thrombosis to inhibit thrombin generation.We further investigated the ability of ProS to function as a ligand of TAM receptors, by using homozygous and heterozygous deficient mice for both the TAM ligands ProS and Gas6. Gas6-/-Pros-/- mice died in utero and showed comparable dramatic bleeding and thrombotic phenotype as described for Pros-/- embryos.In conclusion, like complete ProS deficiency, double deficiency in ProS and Gas6 was lethal, whereas partial ProS deficiency was not. Mice partially deficient in ProS displayed a prothrombotic phenotype, including those with only deficiency in pProS. Purpura fulminans did not occur spontaneously in mice with partial Pros deficiency but developed upon warfarin administration.Thus, the use of different mice models of ProS deficiency can be instrumental in the study of its highly variable thrombotic phenotype and in the investigation of additional roles of ProS in inflammation and autoimmunity through TAM signaling.

Synthesis of triheptanoin and formulation as a solid diet for rodents

Triheptanoin-enriched diets have been successfully used in the experimental treatment of various metabolic disorders. Maximal therapeutic effect is achieved in the context of a ketogenic diet where triheptanoin oil provides 3040% of the daily caloric intake. However, pre-clinical studies using triheptanoin-rich diets are hindered by the difficulty of administering to laboratory animals as a solid foodstuff. In the present study, we successfully synthesized triheptanoin to the highest standards of purity from glycerol and heptanoic acid, using sulfonated charcoal as a catalyst. Triheptanoin oil was then formulated as a solid, stable and palatable preparation using a ketogenic base and a combination of four commercially available formulation agents: hydrophilic fumed silica, hydrophobic fumed silica, microcrystalline cellulose, and talc. Diet compliance and safety was tested on C57Bl/6 mice over a 15-week period, comparing overall status and body weight change. Practical applications: This work provides a complete description of (i) an efficient and cost-effective synthesis of triheptanoin and (ii) its formulation as a solid, stable, and palatable ketogenic diet (triheptanoin-rich; 39% of the caloric intake) for rodents. Triheptanoin-rich diets will be helpful on pre-clinical experiments testing the therapeutic efficacy of triheptanoin in different rodent models of human diseases. In addition, using the same solidification procedure, other oils could be incorporated into rodent ketogenic diet to study their dosage and long-term effects on mammal health and development. This approach could be extremely valuable as ketogenic diet is widely used clinically for epilepsy treatment.

Synthesis of triheptanoin and formulation as a solid diet for rodents

Triheptanoin-enriched diets have been successfully used in the experimental treatment of various metabolic disorders. Maximal therapeutic effect is achieved in the context of a ketogenic diet where triheptanoin oil provides 3040% of the daily caloric intake. However, pre-clinical studies using triheptanoin-rich diets are hindered by the difficulty of administering to laboratory animals as a solid foodstuff. In the present study, we successfully synthesized triheptanoin to the highest standards of purity from glycerol and heptanoic acid, using sulfonated charcoal as a catalyst. Triheptanoin oil was then formulated as a solid, stable and palatable preparation using a ketogenic base and a combination of four commercially available formulation agents: hydrophilic fumed silica, hydrophobic fumed silica, microcrystalline cellulose, and talc. Diet compliance and safety was tested on C57Bl/6 mice over a 15-week period, comparing overall status and body weight change. Practical applications: This work provides a complete description of (i) an efficient and cost-effective synthesis of triheptanoin and (ii) its formulation as a solid, stable, and palatable ketogenic diet (triheptanoin-rich; 39% of the caloric intake) for rodents. Triheptanoin-rich diets will be helpful on pre-clinical experiments testing the therapeutic efficacy of triheptanoin in different rodent models of human diseases. In addition, using the same solidification procedure, other oils could be incorporated into rodent ketogenic diet to study their dosage and long-term effects on mammal health and development. This approach could be extremely valuable as ketogenic diet is widely used clinically for epilepsy treatment.

Model studies on the synthesis of madangamine alkaloids. Assembly of the macrocyclic rings

Using simplified model derivatives, the assembly of the macrocyclic rings of madangamines, including the 13- and 14-membered D rings of madangamines C-E, the all-cis-triunsaturated 15-membered D ring of madangamine A, and the (Z,Z)-unsaturated 11-membered E ring common to madangamines A-E, has been studied.

Synthesis of triheptanoin and formulation as a solid diet for rodents

Triheptanoin-enriched diets have been successfully used in the experimental treatment of various metabolic disorders. Maximal therapeutic effect is achieved in the context of a ketogenic diet where triheptanoin oil provides 3040% of the daily caloric intake. However, pre-clinical studies using triheptanoin-rich diets are hindered by the difficulty of administering to laboratory animals as a solid foodstuff. In the present study, we successfully synthesized triheptanoin to the highest standards of purity from glycerol and heptanoic acid, using sulfonated charcoal as a catalyst. Triheptanoin oil was then formulated as a solid, stable and palatable preparation using a ketogenic base and a combination of four commercially available formulation agents: hydrophilic fumed silica, hydrophobic fumed silica, microcrystalline cellulose, and talc. Diet compliance and safety was tested on C57Bl/6 mice over a 15-week period, comparing overall status and body weight change. Practical applications: This work provides a complete description of (i) an efficient and cost-effective synthesis of triheptanoin and (ii) its formulation as a solid, stable, and palatable ketogenic diet (triheptanoin-rich; 39% of the caloric intake) for rodents. Triheptanoin-rich diets will be helpful on pre-clinical experiments testing the therapeutic efficacy of triheptanoin in different rodent models of human diseases. In addition, using the same solidification procedure, other oils could be incorporated into rodent ketogenic diet to study their dosage and long-term effects on mammal health and development. This approach could be extremely valuable as ketogenic diet is widely used clinically for epilepsy treatment.

Statistical modeling of bowing control applied to violin sound synthesis

Excitation-continuous music instrument control patterns are often not explicitly represented in current sound synthesis techniques when applied to automatic performance. Both physical model-based and sample-based synthesis paradigmswould benefit from a flexible and accurate instrument control model, enabling the improvement of naturalness and realism. Wepresent a framework for modeling bowing control parameters inviolin performance. Nearly non-intrusive sensing techniques allow for accurate acquisition of relevant timbre-related bowing control parameter signals.We model the temporal contour of bow velocity, bow pressing force, and bow-bridge distance as sequences of short Bézier cubic curve segments. Considering different articulations, dynamics, and performance contexts, a number of note classes are defined. Contours of bowing parameters in a performance database are analyzed at note-level by following a predefined grammar that dictates characteristics of curve segment sequences for each of the classes in consideration. As a result, contour analysis of bowing parameters of each note yields an optimal representation vector that is sufficient for reconstructing original contours with significant fidelity. From the resulting representation vectors, we construct a statistical model based on Gaussian mixtures suitable for both the analysis and synthesis of bowing parameter contours. By using the estimated models, synthetic contours can be generated through a bow planning algorithm able to reproduce possible constraints caused by the finite length of the bow. Rendered contours are successfully used in two preliminary synthesis frameworks: digital waveguide-based bowed stringphysical modeling and sample-based spectral-domain synthesis.

Study of the incorporation and the anti-tumour efficacy of radio-iododeoxyuridine according to the cellular cycle and in presence of synthesis inhibitor of thymidine

Résumé La iododeoxyuridine (IdUrd), une fois marqué au 123I ou au 125I, est un agent potentiel pour des thérapies par rayonnements Auger. Cependant, des limitations restreignent son incorporation dans l'ADN. Afin d'augmenter celle-ci, différents groupes ont étudié la fluorodeoxyuridine (FdUrd), qui favorise l'incorporation d'analogue de la thymidine, sans toutefois parvenir à une toxicité associé plus importante. Dans notre approche, 3 lignées cellulaires de glioblastomes humains et une lignée de cancer ovarien ont été utilisées. Nous avons observé, 16 à 24 h après un court pré-traitement à la FdUrd, un fort pourcentage de cellules s'accumulant en phase S. Plus qu'une accumulation, c'était une synchronisation des cellules, celles-ci restant capables d'incorporer la radio-IdIrd et repartant dans le cycle cellulaire. De plus, ces cellules accumulées après un pré-traitement à la FdUrd étaient plus radio-sensibles. Après le même intervalle de 16 à 24 h suivant la FdUrd, les 4 lignées cellulaires ont incorporé des taux plus élevés de radio-IdUrd que sans ce prétraitement. Une corrélation temporelle entre l'accumulation des cellules en phase S et la forte incorporation de radio-IdUrd a ainsi été révélée 16 à 24 h après pré-traitement à la FdUrd. Les expériences de traitement par rayonnements Auger sur les cellules accumulées en phase S ont montré une augmentation significative de l'efficacité thérapeutique de 125I-IdUrd comparé aux cellules non prétraitées à la FdUrd. Une première estimation a permis de déterminer que 100 désintégrations de 125I par cellules étant nécessaires afin d'atteindre l'efficacité thérapeutique. De plus, p53 semble jouer un rôle dans l'induction directe de mort cellulaire après des traitements par rayonnements Auger, comme indiqué par les mesures par FACS d'apoptose et de nécrose 24 et 48 h après le traitement. Concernant les expériences in vivo, nous avons observé une incorporation marquée de la radio-IdUrd dans l'ADN après un pré-traitement à la FdUrd dans un model de carcinomatose ovarienne péritonéale. Une augmentation encore plus importante a été observée après injection intra-tumorale dans des transplants sous-cutanés de glioblastomes sur des souris nues. Ces modèles pourraient être utilisés pour de plus amples études de diffusion de radio-IdUrd et de thérapie par rayonnement Auger. En conclusion, ce travail montre une première application réussie de la FdUrd afin d'accroître l'efficacité de la radio-IdUrd par traitements aux rayonnements Auger. La synchronisation des cellules en phase S combinée avec la forte incorporation de radio-IdUrd dans l'ADN différées après un pré-traitement à la FdUrd ont montré le gain thérapeutique attendu in vitro. De plus, des études in vivo sont tout indiquées après les observations encourageantes d'incorporation de radio-IdUrd dans les models de transplants sous-cutanés de glioblastomes et de tumeurs péritonéales ovariennes. Summary Iododeoxyuridine (IdUrd), labelled with 123I or 125I, could be a potential Auger radiation therapy agent. However, limitations restrict its DNA incorporation in proliferating cells. Therefore, fluorodeoxyuridine (FdUrd), which favours incorporation of thymidine analogues, has been studied by different groups in order to increase radio-IdUrd DNA incorporation, however therapeutic efficacy increase could not be reached. In our approach, 3 human glioblastoma cell lines with different p53 expression and one ovarian cancer line were pre-treated with various FdUrd conditions. We observed a high percentage of cells accumulating in early S phase 16 to 24 h after a short and non-toxic FdUrd pre-treatment. More than an accumulation, this was a synchronization, cells remaining able to incorporate radio-IdUrd and re-entering the cell cycle. Furthermore, the S phase accumulated cells post FdUrd pre-treatment were more radiosensitive. After the same delay of 16 to 24 h post FdUrd pre-treatment, the 4 cell lines were incorporating higher rates of radio-IdUrd compared with untreated cells. A time correlation between S phase accumulation and high radio-IdUrd incorporation was therefore revealed 16 to 24 h post FdUrd pre-treatment. Auger radiation treatment experiments performed on S phase enriched cells showed a significant increase of killing efficacy of 125I-IdUrd compared with cells not pre-treated with FdUrd. A first estimation indicates further that about 100 125I decays were required to reach killing in the targeted cells. Moreover, p53 might play a role on the direct induction of cell death pathways after Auger radiation treatments, as indicated by differential apoptosis and necrosis induction measured by FACS 24 and 48 h after treatment initiation. Concerning in vivo results, we observed a marked DNA incorporation increase of radio-IdUrd after FdUrd pre-treatment in peritoneal carcinomatosis in SCID mice. Even higher incorporation increase was observed after intra-tumoural injection of radio-IdUrd in subcutaneous glioblastoma transplants in nude mice. These tumour models might be further useful for diffusion of radio-IdUrd and Auger radiation therapy studies. In conclusion, these data show a first successful application of thymidine synthesis inhibition able to increase the efficacy of radio-IdUrd Auger radiation treatment. The S phase synchronization combined with a high percentage DNA incorporation of radio-IdUrd delayed post FdUrd pre-treatment provided the expected therapeutic gain in vitro. Further in vivo studies are indicated after the observations of encouraging radio-IdUrd uptake experiments in glioblastoma subcutaneous xenografts and in an ovarian peritoneal carcinomatosis model.

Identifying the preferred subset of enzymatic profiles in nonlinear kinetic metabolic models via multiobjective global optimization and pareto filters

Optimization models in metabolic engineering and systems biology focus typically on optimizing a unique criterion, usually the synthesis rate of a metabolite of interest or the rate of growth. Connectivity and non-linear regulatory effects, however, make it necessary to consider multiple objectives in order to identify useful strategies that balance out different metabolic issues. This is a fundamental aspect, as optimization of maximum yield in a given condition may involve unrealistic values in other key processes. Due to the difficulties associated with detailed non-linear models, analysis using stoichiometric descriptions and linear optimization methods have become rather popular in systems biology. However, despite being useful, these approaches fail in capturing the intrinsic nonlinear nature of the underlying metabolic systems and the regulatory signals involved. Targeting more complex biological systems requires the application of global optimization methods to non-linear representations. In this work we address the multi-objective global optimization of metabolic networks that are described by a special class of models based on the power-law formalism: the generalized mass action (GMA) representation. Our goal is to develop global optimization methods capable of efficiently dealing with several biological criteria simultaneously. In order to overcome the numerical difficulties of dealing with multiple criteria in the optimization, we propose a heuristic approach based on the epsilon constraint method that reduces the computational burden of generating a set of Pareto optimal alternatives, each achieving a unique combination of objectives values. To facilitate the post-optimal analysis of these solutions and narrow down their number prior to being tested in the laboratory, we explore the use of Pareto filters that identify the preferred subset of enzymatic profiles. We demonstrate the usefulness of our approach by means of a case study that optimizes the ethanol production in the fermentation of Saccharomyces cerevisiae.