A regulação endócrina dos balanços de cálcio

RESUMO: As concentrações circulantes de cálcio são notavelmente constantes a despeito das variações diárias na absorção intestinal e na eliminação renal deste elemento. A regulação da calcémia é um sistema complexo que compreende vários factores controladores (a calcémia, a fosforémia, as concentrações circulantes de paratormona (PTH) e calcitriol além de muitos outros factores como hormonas esteróides em geral, outros iões como o magnésio e outros factores hormonais) e vários órgãos alvo (glândulas paratiroideias, osso, rim e intestino). As respostas dos órgãos alvo também são muito variadas. No caso mais simples, a cristalização de sais de cálcio corresponde a uma mudança de fase em que participam moléculas orgânicas que a iniciam, aceleram ou inibem. Em geral a combinação de um factor controlador com o respectivo receptor de membrana (para polipeptídeos ou iões) ou intracelular (hormonas esteróides) é apenas o primeiro passo de uma cadeia bioquímica que introduz uma enorme amplificação na resposta. A esta variedade de mecanismos de resposta correspondem grandes diferenças nos tempos de resposta que podem ser de minutos a semanas. É hoje possível “observar” (medir) com apreciável rigor nos líquidos biológicos (sangue, urina, fezes, etc.) os factores mais importantes do sistema de regulação da calcémia (cálcio, fósforo, paratormona e calcitriol) assim como administrar estes factores em experiências agudas. Esta possibilidade reflecte – se na literatura neste campo que tem vindo a crescer. O advento das técnicas da biologia molecular tem permitido a caracterização molecular de algumas das disfunções da homeostase do cálcio e é de esperar um diagnóstico fisiopatológico cada vez mais rigoroso dessas disfunções. Com o avanço dos conhecimentos nesta área que não cessa de aumentar temos cada vez maiores capacidades para fazer diagnósticos e é cada vez mais difícil interpretar com rigor os correspondentes quadros metabólicos. A análise ou síntese de sistemas complexos é a actividade mais nobre dos engenheiros que lhes permite desenhar pontes, diques, barcos, aviões ou automóveis. Com o aparecimento de computadores de médio ou grande porte foi – lhes possível utilizar descrições matemáticas não só para desenhar sistemas como ainda para interpretar eventuais falhas na sua operação. Essas descrições matemáticas consistem numa sequência de operações realizadas num computador segundo um “programa informático” que receberam a designação genérica de modelos, por analogia com as famosas leis (equações) da física que foram deduzidas a partir de um certo número de postulados e que permitem representar matematicamente processos físicos. As famosas leis de Newton são talvez os exemplos mais famosos de “modelos” de sistemas físicos. A introdução de modelos matemáticos em biologia e particularmente em medicina só se deu recentemente.MÉTODOS No trabalho que aqui se apresenta construiu - se um modelo simplificado da homeostase do cálcio destinado ao cálculo de variáveis observáveis (concentrações de cálcio, fósforo, PTH e calcitriol) de modo a poderem comparar-se valores calculados com valores observados. A escolha dos componentes do modelo foi determinada pela nossa experiência clínica e pela informação fisiopatológica e clínica publicada. Houve a preocupação de construir o modelo de forma modular de modo a ser possível a sua expansão sem grandes transformações na descrição matemática (e informática) já existente. Na sua fase actual o modelo não pode ser usado como instrumento de diagnóstico. É antes uma ferramenta destinada a esclarecer “em princípio” mecanismos fisiopatológicos. Usou – se o modelo para simular um certo número de observações publicadas e para exemplificar a sua eventual aplicação clínica na simulação de situações hipotéticas e na análise de possíveis mecanismos fisiopatológicos responsáveis por situações de hipo ou hipercalcémias. Simultaneamente fez – se uma análise dos dados acumulados relativos a doentes vistos no Serviço de Endocrinologia do Instituto Português de Oncologia de Francisco Gentil – Centro Regional Oncológico de Lisboa, S.A. CONCLUSÕES Numa população de 894 doentes com patologias variadas do Instituto Português de Oncologia de Lisboa os valores da calcémia tiveram uma distribuição normal unimodal com uma média de 9.56 mg/dl, e um erro padrão de 0.41 mg/dl. Estas observações sugerem que a calcémia está sujeita a regulação. A partir dos resultados publicados em que o metabolismo do cálcio foi perturbado por infusões de cálcio, calcitriol ou PTH, de estudos bioquímicos e fisiológicos sobre os mecanismos de acção de factores controladores da calcémia e do estudo do comportamento de órgãos alvo (paratiroideias, intestino, osso e rim) foi possível construir um modelo matemático de parâmetros concentrados do sistema de regulação da calcémia. As expressões analíticas usadas foram baseadas na cinética enzimática de modo a que os seus parâmetros tivessem um significado físico ou fisiológico simples. O modelo revelou apreciável robustez e flexibilidade. É estável quando não perturbado e transita entre estados estacionários quando perturbado. Na sua forma actual gera simulações que reproduzem satisfatoriamente um número apreciável de dados experimentais colhidos em doentes. Isto não significa que possa ser usado como instrumento de diagnóstico aplicável a doentes individuais. O desenho do modelo comporta a adição posterior de novas relações quando surgirem situações para as quais se revele insuficiente. A utilização exaustiva do modelo permitiu explicitar aspectos do metabolismo do cálcio que ou não estão contidas na sua formulação actual – o aparecimento de hipertrofia ou de adenomas das paratiroideias e as alterações na estrutura óssea , a participação de outros factores controladores – magnésio, ou estão insuficientemente descritas – alterações do metabolismo do fósforo nos hipoparatiroidismos. A análise dos dados relativos aos doentes do Serviço de Endocrinologia do IPO permitiu o início da caracterização dos tipos de patologia que representam e de possíveis mecanismos fisiopatológicos subjacentes. Estas observações são o ponto de partida para análises futuras. São exemplos das relações encontradas: a distribuição dos doentes por dois grandes grupos conforme a calcémia é determinada pelas concentrações circulantes de PTH ou estas são determinadas pela calcémia; a distribuição sazonal das concentrações de Vit. D25. no sangue; a correlação negativa entre estas e as concentrações de PTH no sangue. Também foi possível extrair a cinética do controlo da PTH sobre a síntese de calcitriol. O estudo dos níveis circulantes de PTH no pós-operatório imediato de doentes paratiroidectomizados permitiu determinar as suas taxas de degradação metabólica. O modelo permitiu simular as relações Ca/PTH no sangue, Ca/Fracção excretada da carga tubular, Ca/P no sangue para valores normais ou altos de Ca. Foram feitas simulações de situações fisiopatológicas (em “doentes virtuais”): infusões crónicas de cálcio, PTH e calcitriol; alterações no comportamento de receptores. Estas simulações correspondem a experiências que não podem ser realizadas em humanos. São exemplos da utilização do modelo na exploração de possíveis mecanismos fisiopatológicos através da observação de resultados quantitativos inacessíveis à intuição. O modelo foi útil em duas fases do trabalho: Primeiro, durante a sua síntese implicou uma escolha criticamente selectiva de informação, sua análise quantitativa e processamento, uma explicitação rigorosa (analítica) das relações funcionais entre os controladores e as variáveis e da sua integração numa estrutura global; Segundo, a simulação de situações experimentais ou clínicas (dados do Serviço de Endocrinologia do IPO) em doentes obrigou a explicitar raciocínios fisiopatológicos habitualmente formulados em bases puramente intuitivas. Esta prática revelou comportamentos óbvios após as simulações – acção reduzida das infusões PTH (simulação de hiperparatiroidismos primários) enquanto não há inibição total da respectiva secreção, necessidade de aumento da massa secretora da paratiroideia nas insuficiências renais avançadas, etc. A síntese e utilização do modelo não implicaram uma preparação matemática avançada e foram possíveis mercê da disponibilidade de “software” interactivo especificamente desenhado para a simulação de sistemas dinâmicos em que os programas se escrevem em inglês usando a simbologia simples da álgebra elementar. A função nobre de modelos desta natureza é semelhante à dos modelos usados pelos físicos desde o século XVII: permitir explicações de carácter geral funcionando como uma ferramenta intelectual para manipulação de conceitos e para a realização de “experiências pensadas” (“thought experiments”) respeitando certos princípios físicos (princípios de conservação) que estabelecem as fronteiras da realidade. -------ABSTRACT: Calcium blood levels are remarkably constant despite great variations in calcium daily intake, intestinal absorption and renal excretion. The regulation of the calcium concentration in the blood is achieved by a complex system that includes several controller factors (mainly the serum levels of calcium, phosphorus, parathyroid hormone (PTH) and calcitriol but also of steroid hormones, ions such as magnesium and other hormonal factors) and several target organs (parathyroid glands, bone, kidney and intestine). The functional response to the controlling factors obeys a variety of kinetics. The precipitation of calcium salts is a simple phase transition in which organic molecules may provide nucleation centres or inhibit the process. The combination of a controller factor with its receptor located in the cell membrane (for peptides or ions) or in the nucleus (for steroid hormones) is only the first step of a biochemical chain that introduces a huge amplification in the response. To this great variability of response we have to add the times of response that vary from minutes to weeks. It is possible to “observe” (measure) with great accuracy in biological fluids (blood, urine, faeces, etc.) the most important factors intervening in the calcium regulation (calcium, phosphorus, PTH and calcitriol). The response of the system to acute infusions of the controlling factors has also been studied. Using molecular biology techniques it has been possible to characterize some calcium homeostasis dysfunctions and better physiopathological diagnosis are expected. With the increasingly new knowledge in this area we have better capacity to diagnose but it is harder to explain correctly the underlying metabolic mechanisms. The analysis or synthesis of complex systems is the noble activity of engineers that enables them to draw bridges, dams, boats, airplanes or cars. With the availability of medium-large frame computers it was possible to use mathematical descriptions not only to draw systems but also to explain flaws in its operations. These mathematical descriptions are generally known as models by analogy with the laws (equations) of physics that allow the mathematical description of physical processes. In practice it is not possible to find general solutions for the mathematical descriptions of complex systems but (numeric) computations for specific situations can be obtained with digital computers. The introduction of mathematical models in biology and particularly in medicine is a recent event. METHODS In this thesis a simplified model of calcium homeostasis was built that enables the computation of observable variables (concentrations of calcium, phosphorus, PTH and calcitriol) and allows the comparison between the simulation values and observed values. The choice of the model’s components was made according to our clinical experience and to the published clinical and physiopathological data. The model has a modular design that allows future expansions with minor alterations in its structure. In its present form the model cannot be used for diagnosis. It is a tool designed to enlighten physiopathological processes. To exemplify its possible clinical application in the simulation of hypothetical situations and in the analysis of possible mechanisms responsible for hypo or hypercalcemias the model was used to simulate a certain number of published observations. An analysis of clinical and laboratory data from the Endocrinology Department of the Portuguese Cancer Institute (I.P.O.F.G.-C.R.O.L.,S.A.) is also presented. CONCLUSIONS In a population of 188 patients without an identifiable disease of the calcium metabolism at the Portuguese Cancer Institute the calcemia levels had a unimodal distribution with an average of 9.56 mg/dL and a S.E.M of 0.41 mg/dL. This observation confirms that serum calcium is regulated. Using published data; in which calcium metabolism was disrupted by calcium, PTH or calcitriol infusions; from biochemical and physiological studies of the action of controller factors on the calcemia; in which the response of target organs (parathyroid glands, intestine, bone, kidney) was studied it was possible to build a mathematical model of concentrated parameters of the calcium homeostasis. Analytical expressions used were based on enzymatic kinetics. The model is flexible and robust. It is stable when not disturbed and changes between steady states when disturbed. In its present form it provides simulations that reproduce closely a number of experimental clinical data. This does not mean that it can be used as a diagnostic tool for individual patients. The exhaustive utilisation of the model revealed the need of future expansions to include aspects of the calcium metabolism not included in its present form –hypertrophy or adenomas of the parathyroid glands, bone structure changes, participation of other controller factors such as magnesium – or insufficiently described – phosphate metabolism in hypoparathyroidism. The analysis of the data collected from the I.P.O.’s Endocrinology Department allowed the initial characterization of the different pathologies represented and of their possible physiopathological mechanisms. These observations are a starting point for future analysis. As examples of the relations found were: the distribution of patients in two groups according to the dependency of calcium by PTH levels or PTH levels by calcium concentration; the seasonal distribution of the serum concentrations of D25; its negative correlation with PTH concentration. It was also possible to extract the kinetics of the control of the synthesis of calcitriol by PTH. The analysis of immediate post-surgical levels of PTH in parathyroidectomized patients allowed the determination of its metabolic clearance. The model also allowed the simulation of the relations between Ca/PTH in blood, serum Ca/Fraction of tubular load excreted and Ca/P in blood for normal and high values of calcium. Simulations were made of pathological situations (in “virtual patients”): chronic infusions of calcium, PTH and calcitriol; changes in the characteristics of receptors. These simulations are not possible in real persons. They are an example of the use of this model in exploring possible mechanisms of disease through the observation of quantitative results not accessible to simple intuition. This model was useful in two phases: Firstly, its construction required a careful choice of data, its quantitative analysis and processing, an analytical description of the relations between controller factors and variables and their integration in a global structure. Secondly, the simulation of experimental or clinical (I.P.O.’s Endocrinology Department) data implied testing physiopathological explanations that previously were based on intuition. The construction and utilisation of the model didn’t demand an advanced mathematical preparation since user-friendly interactive software was used. This software was specifically designed for the simulation of dynamic systems. The programs are written in English using elementary algebra symbols. The essential function of this type of models is identical to that of those used by physicists since the XVII century which describe quantitatively natural processes and are an intellectual tool for the manipulation of concepts and the performance of “thought experiments” based in certain physical principles (conservation principles) that are the frontiers of reality.------------------RESUMÉE: Les concentrations circulantes de calcium sont constantes même pendant des variations de l’absorption intestinale et de l’élimination rénale de cet élément. La régulation de la calcémie est un système complexe qui comprend plusieurs éléments contrôleurs (la calcémie, la phosphorémie, les concentrations circulantes de l’hormone parathyroïdienne (PTH) e du calcitriol et d’autres comme les hormones stéroïdes ou des ions comme le magnésium) et plusieurs organes (glandes parathyroïdiennes, l’os, le rein et l’intestin). Les réponses de ces organes sont variées. Dans le cas plus simple, la cristallisation des sels de calcium correspond à un changement de phase dans lequel y participent des molécules organiques que la débutent, l’accélèrent ou l’inhibent. Généralement la combinaison d’un élément contrôleur avec leur récepteur de membrane (pour les peptides ou les ions) ou intracellulaire (pour les hormones stéroïdes) n’est que le premier pas d’une chaîne biochimique qu’introduit une grande amplification de la réponse. A cette variété de réponses correspondent des grandes différences des temps de réponses qu’y vont des minuits a semaines. Il est possible « observer » (mesurer) dans les fluides biologiques (sang, urine, fèces, etc.) les éléments plus importants du système de régulation de la calcémie (calcium, phosphate, PTH et le calcitriol) et les administrer en expérimentes aigus. Cette possibilité est visible dans la littérature publiée dans ce domaine qui est en croissance permanente. L’avenir des techniques de biologie moléculaire a permis caractériser des nombreuses dysfonctions de la régulation de la calcémie et on attend un diagnostique physiopathologique de ces dysfonctions chaque fois plus rigoureuses. Les connaissances dans ce domaine s’agrandissent et on a de plus de capacités pour faire des diagnostiques et il est chaque fois plus difficile les interpréter. L’analyse ou synthèse de systèmes complexes est l’activité plus noble des ingénieurs qui les permit dessiner des ponts, bateaux, avions ou automobiles. Avec des ordinateurs de médium ou grand port il les est possible utiliser descriptions mathématiques pour dessiner les systèmes et interpréter des éventuelles fautes d’opération. Ces descriptions mathématiques sont une séquence d’opérations réalisées dans un ordinateur selon « un programme informatique » qui ont reçu la désignation générique de modèles, pour analogie avec les équations de la physique qui ont été déduits d’un nombre de postulées et qu’ont permit représenter des processus physiques en équations mathématiques. Les fameuses équations de Newton sont peut-être les exemples plus connus des systèmes physiques. L’introduction des modèles mathématiques en biologie et en particulier en médecine est un évènement récent. Dans ce travaille, on a construit un modèle simplifié de l’homéostasie du calcium pour calculer les variables observables (concentrations de calcium, phosphate, PTH et calcitriol) pour les comparer. Les choix des components a été déterminés par notre expérience clinique et par l’information physiopathologique et clinique publiée. Le modèle a été construit de façon modulaire ce que permit leur postérieur expansion sans des grandes altérations dans la description mathématique et informatique déjà existante. Dans cette forme le modèle ne peut être utilisé comme un instrument de diagnostique. Il est un outil pour éclairer la physiopathologie. Le modèle a été utilisé pour simuler un certain nombre d’observations publiées et pour exemplifier leur possible utilisation clinique dans la simulation des hypothèses et de la physiopathologie des situations d’hypo ou hypercalcémie. On a fait une analyse des éléments des procès cliniques des malades observées dans le Service d’Endocrinologie de l’IPOFG-CROL, SA. Dans une population de 894 malades avec des différentes pathologies les valeurs de calcémie on une distribution uni modale avec une Médie de 9.56 mg/dL et une erreur standard de 0.41 mg/dL. Ces observations suggèrent que la calcémie soit sujette de régulation. En utilisant des résultats de travaux publiés dans lesquels le métabolisme du calcium a été changé par des infusions de calcium, calcitriol ou PTH, des études biochimiques et physiologiques sur des mécanismes d’action des éléments contrôleurs de la calcémie et de l’étude du comportement des organes cible (parathyroïdes, intestin, rein, os), il a été possible de construire un modèle mathématique de paramètres concentrés du système de régulation de la calcémie. Les expressions analytiques utilisées ont été basées sur la cinétique enzymatique de façon à que les paramètres aient eu une signification physique ou biologique. Le modèle est stable quand il n’est pas perturbé et transit entre états stationnaires quand il est sujet a des perturbations. A ce moment il fait des simulations qui reproduisent de façon satisfaisant un nombre d’observations expérimentales. La construction du modèle permit l’addiction de nouvelles relations dans les cas ou il est insuffisant. L’utilisation exhaustive du modèle a permit expliciter des aspects du métabolisme du calcium qui y ne sont pas compris – l’hyperplasie ou la formation des adénomes des parathyroïdes, les altérations de la structure des os, la participation d’outres éléments régulateurs (magnésium), ou sont insuffisamment décrites – les altérations du métabolisme des phosphates dans l’hypoparathyroidism. L’analyse de l’information des malades du Service d’Endocrinologie a permit caractériser les pathologies représentées et leurs possibles mécanismes physiopathologiques. Ces observations sont le point de départ pour les analyses futures. Sont des exemples des relations trouvées: la distribution des malades par deux groupes: ceux dans lequel la calcémie est déterminée par la PTH ou ceux dans lesquels la PTH est déterminée par la calcémie; la distribution sazonale de la concentration de la vitamine D; la corrélation négative entre la vitamine D et la PTH. On a eu la possibilité de déduire la cinétique de control de la PTH sur la synthèse du calcitriol. L’étude des niveaux circulants de PTH sur des sujets parathyroidectomisées a permit déduire leur taux de dégradation métabolique. Le modèle a permit simuler les relations Ca/PTH dans le sang, Ca/fraction éliminée par le rein, Ca/P dans le sang pour des valeurs normales ou hautes de calcium. On a fait des simulations de situations physiopathologiques (dans “malades virtuelles”): Infusions chroniques de calcium, PTH ou calcitriol; altérations des récepteurs. Ces simulations ne peuvent pas être réalisées dans les humains. Sont des exemples d’utilisation du modèle dans l’exploration des possibles mécanismes de la physiopathologie en observant des résultats quantitatifs inaccessibles à l’intuition. Le modèle a été utile pendant deux étapes des travaux: La première, dans sa construction on a choisi l’information disponible, son analyse quantitative, l’explicitation rigoureuse (analytique) des relations fonctionnelles entre les contrôleurs et les variables et sa intégration dans une structure globale. La deuxième, la simulation de situations expérimentales ou cliniques (du Service d’Endocrinologie) a obligé d’expliciter des raisonnements physiopathologiques généralement formulés utilisant l’intuition. Cette pratique a montré des comportements – action réduite des infusions de PTH (jusqu’à l’inhibition totale de leur respective sécrétion), nécessité d’augmenter la masse sécréteuse de la parathyroïde dans les insuffisants rénales, etc. La synthèse et utilisation du modèle n’ont pas besoin d’une formation avancée en mathématique et sont possibles grâce à un programme interactif qui a été conçu pour la simulation des systèmes dynamiques dans lesquels le programme se construit en anglais en utilisant la symbolique élémentaire de l’algèbre. La fonction noble de ces modèles est semblable à celles des physiques du XVII siècle: Permettre établir explications générales en fonctionnant comme un outil intellectuel pour manipuler des concepts et pour la réalisation d’expérimentes pensées en respectant certains principes de la physique (principe de la conservation) qu’établissent les frontières de la réalité.

Aplicação de um modelo de programação linear inteira mista (PLIM) na expansão da rede de distribuição de energia elétrica em Angola

Em Angola, apenas cerca de 30% da população tem acesso à energia elétrica, nível que decresce para valores inferiores a 10% em zonas rurais mais remotas. Este problema é agravado pelo facto de, na maioria dos casos, as infraestruturas existentes se encontrarem danificadas ou não acompanharem o desenvolvimento da região. Em particular na capital angolana, Luanda que, sendo a menor província de Angola, é a que regista atualmente a maior densidade populacional. Com uma população de cerca de 5 milhões de habitantes, não só há frequentemente problemas relacionados com a falha do fornecimento de energia elétrica como há ainda uma percentagem considerável de municípios onde a rede elétrica ainda nem sequer chegou. O governo de Angola, no seu esforço de crescimento e aproveitamento das suas enormes potencialidades, definiu o setor energético como um dos fatores críticos para o desenvolvimento sustentável do país, tendo assumido que este é um dos eixos prioritários até 2016. Existem objetivos claros quanto à reabilitação e expansão das infraestruturas do setor elétrico, aumentando a capacidade instalada do país e criando uma rede nacional adequada, com o intuito não só de melhorar a qualidade e fiabilidade da rede já existente como de a aumentar. Este trabalho de dissertação consistiu no levantamento de dados reais relativamente à rede de distribuição de energia elétrica de Luanda, na análise e planeamento do que é mais premente fazer relativamente à sua expansão, na escolha dos locais onde é viável localizar novas subestações, na modelação adequada do problema real e na proposta de uma solução ótima para a expansão da rede existente. Depois de analisados diferentes modelos matemáticos aplicados ao problema de expansão de redes de distribuição de energia elétrica encontrados na literatura, optou-se por um modelo de programação linear inteira mista (PLIM) que se mostrou adequado. Desenvolvido o modelo do problema, o mesmo foi resolvido por recurso a software de otimização Analytic Solver e CPLEX. Como forma de validação dos resultados obtidos, foi implementada a solução de rede no simulador PowerWorld 8.0 OPF, software este que permite a simulação da operação do sistema de trânsito de potências.

Bases para o estudo da genética de populações dos Hymenoptera em geral e dos Apinae sociais em particular

This paper deals with problems on population genetics in Hymenoptera and particularly in social Apidae. 1) The studies on populations of Hymenoptera were made according to the two basic types of reproduction: endogamy and panmixia. The populations of social Apinae have a mixed method of reproduction with higher percentage of panmixia and a lower of endogamy. This is shown by the following a) males can enter any hive in swarming time; b) males of Meliponini are expelled from hives which does not need them, and thus, are forced to look for some other place; c) Meliponini males were seen powdering themselves with pollen, thus becoming more acceptable in any other hive. The panmixia is not complete owing to the fact that the density of the breeding population as very low, even in the more frequent species as low as about 2 females and 160 males per reproductive area. We adopted as selection values (or survival indices) the expressions according to Brieger (1948,1950) which may be summarised as follows; a population: p2AA + ²pq Aa + q2aa became after selection: x p2AA + 2pq Aa + z q²aa. For alge-braics facilities Brieger divided the three selective values by y giving thus: x/y p2 AA + y/y 2 pq Aa + z/y q²aa. He called x/y of RA and z/y of Ra, that are survival or selective index, calculated in relation to the heterozygote. In our case all index were calculated in relation to the heterozygote, including the ones for haploid males; thus we have: RA surveval index of genotype AA Ra surveval index of genotype aa R'A surveval index of genotype A R'a surveval index of genotype a 1 surveval index of genotype Aa The index R'A ande R'a were equalized to RA and Ra, respectively, for facilities in the conclusions. 2) Panmitic populations of Hymenoptera, barring mutations, migrations and selection, should follow the Hardy-Weinberg law, thus all gens will be present in the population in the inicial frequency (see Graphifc 1). 3) Heterotic genes: If mutation for heterotic gene ( 1 > RA > Ra) occurs, an equilibrium will be reached in a population when: P = R A + Ra - 2R²a _____________ (9) 2(R A + Ra - R²A - R²a q = R A + Ra - 2R²A _____________ (10) 2(R A + Ra - R²A - R²a A heterotic gene in an hymenopteran population may be maintained without the aid of new mutation only if the survival index of the most viable mutant (RA) does not exced the limiting value given by the formula: R A = 1 + √1+Ra _________ 4 If RA has a value higher thah the one permitted by the formula, then only the more viable gene will remain present in the population (see Graphic 10). The only direct proof for heterotic genes in Hymenoptera was given by Mackensen and Roberts, who obtained offspring from Apis mellefera L. queens fertilized by their own sons. Such inbreeding resulted in a rapid loss of vigor the colony; inbred lines intercrossed gave a high hybrid vigor. Other fats correlated with the "heterosis" problem are; a) In a colony M. quadrifasciata Lep., which suffered severely from heat, the percentage of deths omong males was greater .than among females; b) Casteel and Phillips had shown that in their samples (Apis melifera L). the males had 7 times more abnormalities tian the workers (see Quadros IV to VIII); c) just after emerging the males have great variation, but the older ones show a variation equal to that of workers; d) The tongue lenght of males of Apis mellifera L., of Bombus rubicundus Smith (Quadro X), of Melipona marginata Lep. (Quadro XI), and of Melipona quadrifasciata Lep. Quadro IX, show greater variationthan that of workers of the respective species. If such variation were only caused by subviables genes a rapid increasse of homozigoty for the most viable alleles should be expected; then, these .wild populations, supposed to be in equilibrium, could .not show such variability among males. Thus we conclude that heterotic genes have a grat importance in these cases. 4) By means of mathematical models, we came to the conclusion tht isolating genes (Ra ^ Ra > 1), even in the case of mutations with more adaptability, have only the opor-tunity of survival when the population number is very low (thus the frequency of the gene in the breeding population will be large just after its appearence). A pair of such alleles can only remain present in a population when in border regions of two races or subspecies. For more details see Graphics 5 to 8. 5) Sex-limited genes affecting only females, are of great importance toHymenoptera, being subject to the same limits and formulas as diploid panmitic populations (see formulas 12 and 13). The following examples of these genes were given: a) caste-determining genes in the genus Melipona; b) genes permiting an easy response of females to differences in feeding in almost all social Hymenoptera; c) two genes, found in wild populations, one in Trigona (Plebéia) mosquito F. SMITH (quadro XII) and other in Melipona marginata marginata LEP. (Quadro XIII, colonies 76 and 56) showing sex-limited effects. Sex-limited genes affecting only males do not contribute to the plasticity or genie reserve in hymenopteran populations (see formula 14). 6) The factor time (life span) in Hymenoptera has a particular importance for heterotic genes. Supposing one year to be the time unit and a pair of heterotic genes with respective survival indice equal to RA = 0, 90 and Ra = 0,70 to be present; then if the life time of a population is either one or two years, only the more viable gene will remain present (see formula 11). If the species has a life time of three years, then both alleles will be maintained. Thus we conclude that in specis with long lif-time, the heterotic genes have more importance, and should be found more easily. 7) The colonies of social Hymenoptera behave as units in competition, thus in the studies of populations one must determine the survival index, of these units which may be subdivided in indice for egg-laying, for adaptive value of the queen, for working capacity of workers, etc. 8) A study of endogamic hymenopteran populations, reproduced by sister x brother mating (fig. 2), lead us to the following conclusions: a) without selection, a population, heterozygous for one pair of alleles, will consist after some generations (theoretically after an infinite number of generation) of females AA fecundated with males A and females aa fecundated with males a (see Quadro I). b) Even in endogamic population there is the theoretical possibility of the presence of heterotic genes, at equilibrium without the aid of new mutations (see Graphics 11 and 12), but the following! conditions must be satisfied: I - surveval index of both homozygotes (RA e Ra) should be below 0,75 (see Graphic 13); II - The most viable allele must riot exced the less viable one by more than is permited by the following formula (Pimentel Gomes 1950) (see Gra-fic 14) : 4 R5A + 8 Ra R4A - 4 Ra R³A (Ra - 1) R²A - - R²a (4 R²a + 4 Ra - 1) R A + 2 R³a < o Considering these two conditions, the existance of heterotic genes in endogamic populations of Hymenoptera \>ecames very improbable though not - impossible. 9) Genie mutation offects more hymenopteran than diploid populations. Thus we have for lethal genes in diploid populations: u = q2, and in Hymenoptera: u = s, being u the mutation ratio and s the frequency of the mutant in the male population. 10) Three factors, important to competition among species of Meliponini were analysed: flying capacity of workers, food gathering capacity of workers, egg-laying of the queen. In this connection we refer to the variability of the tongue lenght observed in colonies from several localites, to the method of transporting the pollen in the stomach, from some pots (Melliponi-ni storage alveolus) to others (e. g. in cases of pillage), and to the observation that the species with the most populous hives are almost always the most frequent ones also. 11) Several defensive ways used for Meliponini to avoid predation are cited, but special references are made upon the camouflage of both hive (fig. 5) and hive entrance (fig. 4) and on the mimetism (see list in page ). Also under the same heading we described the method of Lestrimelitta for pillage. 12) As mechanisms important for promoting genetic plasticity of hymenopteran species we cited: a) cytological variations and b) genie reserve. As to the former, duplications and numerical variations of chromosomes were studied. Diprion simile ATC was cited as example for polyploidy. Apis mellife-ra L. (n •= 16) also sugests polyploid origen since: a) The genus Melipona, which belongs to a" related tribe, presents in all species so far studied n = 9 chromosomes and b) there occurs formation of dyads in the firt spermatocyte division. It is su-gested that the origin of the sex-chromosome of Apis mellifera It. may be related to the possible origin of diplo-tetraploidy in this species. With regards to the genie reserve, several possible types of mutants were discussed. They were classified according to their survival indices; the heterotic and neutral mutants must be considered as more important for the genie reserve. 13) The mean radius from a mother to a daghter colony was estimated as 100 meters. Since the Meliponini hives swarm only once a year we may take 100 meters a year as the average dispersion of female Meliponini in ocordance to data obtained from Trigona (tetragonisca) jaty F. SMITH and Melipona marginata LEP., while other species may give different values. For males the flying distance was roughly estimated to be 10 times that for females. A review of the bibliography on Meliponini swarm was made (pg. 43 to 47) and new facts added. The population desity (breeding population) corresponds in may species of Meliponini to one male and one female per 10.000 square meters. Apparently the males are more frequent than the females, because there are sometimes many thousands, of males in a swarm; but for the genie frequency the individuals which have descendants are the ones computed. In the case of Apini and Meliponini, only one queen per hive and the males represented by. the spermatozoos in its spermateca are computed. In Meliponini only one male mate with the queen, while queens of Apis mellijera L. are fecundated by an average of about 1, 5 males. (Roberts, 1944). From the date cited, one clearly sees that, on the whole, populations of wild social bees (Meliponini) are so small that the Sewall Wright effect may become of great importance. In fact applying the Wright's formula: f = ( 1/aN♂ + 1/aN♀) (1 - 1/aN♂ + 1/aN♀) which measures the fixation and loss of genes per generation, we see that the fixation or loss of genes is of about 7% in the more frequent species, and rarer species about 11%. The variation in size, tergite color, background color, etc, of Melipona marginata Lep. is atributed to this genetic drift. A detail, important to the survival of Meliponini species, is the Constance of their breeding population. This Constance is due to the social organization, i. e., to the care given to the reproductive individuals (the queen with its sperm pack), to the way of swarming, to the food storage intended to control variations of feeding supply, etc. 14) Some species of the Meliponini are adapted to various ecological conditions and inhabit large geographical areas (e. g. T. (Tetragonisca jaty F. SMITH), and Trigona (Nanno-trigona testaceicornis LEP.) while others are limited to narrow regions with special ecological conditions (e. g. M. fuscata me-lanoventer SCHWARZ). Other species still, within the same geographical region, profit different ecological conditions, as do M. marginata LEP. and M. quadrifasciata LEP. The geographical distribution of Melipona quadrifasciata LEP. is different according to the subspecies: a) subsp anthidio-des LEP. (represented in Fig. 7 by black squares) inhabits a region fron the North of the S. Paulo State to Northeastern Brazil, ,b) subspecies quadrifasciata LEP., (marked in Fig. 7 with black triangles) accurs from the South of S. Paulo State to the middle of the State of Rio Grande do Sul (South Brazil). In the margined region between these two areas of distribution, hi-brid colonies were found (Fig. 7, white circles); they are shown with more details in fig. 8, while the zone of hybridization is roughly indicated in fig. 9 (gray zone). The subspecies quadrifasciata LEP., has 4 complete yellow bands on the abdominal tergites while anthidioides LEP. has interrupted ones. This character is determined by one or two genes and gives different adaptative properties to the subspecies. Figs. 10 shows certains meteorological isoclines which have aproximately the same configuration as the limits of the hybrid zone, suggesting different climatic adaptabilities for both genotypes. The exis-tance of a border zone between the areas of both subspecies, where were found a high frequency of hybrids, is explained as follows: being each subspecies adapted to a special climatic zone, we may suppose a poor adaptation of either one in the border region, which is also a region of intermediate climatic conditions. Thus, the hybrids, having a combination of the parent qualities, will be best adapted to the transition zone. Thus, the hybrids will become heterotic and an equilibrium will be reached with all genotypes present in the population in the border region.

Detection of EPO doping and blood doping: the haematological module of the Athlete Biological Passport.

The increase of the body's capacity to transport oxygen is a prime target for doping athletes in all endurance sports. For this pupose, blood transfusions or erythropoiesis stimulating agents (ESA), such as erythropoietin, NESP, and CERA are used. As direct detection of such manipulations is difficult, biomarkers that are connected to the haematopoietic system (haemoglobin concentration, reticulocytes) are monitored over time (Athlete Biological Passport (ABP)) and analyzed using mathematical models to identify patterns suspicious of doping. With this information, athletes can either be sanctioned directly based on their profile or targeted with conventional doping tests. Key issues for the appropriate use of the ABP are correct targeting and use of all available information (e.g. whereabouts, cross sectional population data) in a forensic manner. Future developments of the passport include the correction of all concentration-based variables for shifts in plasma volume, which might considerably increase sensitivity. New passport markers from the genomic, proteomic, and metabolomic level might add further information, but need to be validated before integration into the passport procedure. A first assessment of blood data of federations that have implemented the passport show encouraging signs of a decreased blood-doping prevalence in their athletes, which adds scientific credibility to this innovative concept in the fight against ESA- and blood doping. Copyright © 2012 John Wiley & Sons, Ltd.

Study and implementation of some quantitative trading models

Quantitative or algorithmic trading is the automatization of investments decisions obeying a fixed or dynamic sets of rules to determine trading orders. It has increasingly made its way up to 70% of the trading volume of one of the biggest financial markets such as the New York Stock Exchange (NYSE). However, there is not a signi cant amount of academic literature devoted to it due to the private nature of investment banks and hedge funds. This projects aims to review the literature and discuss the models available in a subject that publications are scarce and infrequently. We review the basic and fundamental mathematical concepts needed for modeling financial markets such as: stochastic processes, stochastic integration and basic models for prices and spreads dynamics necessary for building quantitative strategies. We also contrast these models with real market data with minutely sampling frequency from the Dow Jones Industrial Average (DJIA). Quantitative strategies try to exploit two types of behavior: trend following or mean reversion. The former is grouped in the so-called technical models and the later in the so-called pairs trading. Technical models have been discarded by financial theoreticians but we show that they can be properly cast into a well defined scientific predictor if the signal generated by them pass the test of being a Markov time. That is, we can tell if the signal has occurred or not by examining the information up to the current time; or more technically, if the event is F_t-measurable. On the other hand the concept of pairs trading or market neutral strategy is fairly simple. However it can be cast in a variety of mathematical models ranging from a method based on a simple euclidean distance, in a co-integration framework or involving stochastic differential equations such as the well-known Ornstein-Uhlenbeck mean reversal ODE and its variations. A model for forecasting any economic or financial magnitude could be properly defined with scientific rigor but it could also lack of any economical value and be considered useless from a practical point of view. This is why this project could not be complete without a backtesting of the mentioned strategies. Conducting a useful and realistic backtesting is by no means a trivial exercise since the \laws" that govern financial markets are constantly evolving in time. This is the reason because we make emphasis in the calibration process of the strategies' parameters to adapt the given market conditions. We find out that the parameters from technical models are more volatile than their counterpart form market neutral strategies and calibration must be done in a high-frequency sampling manner to constantly track the currently market situation. As a whole, the goal of this project is to provide an overview of a quantitative approach to investment reviewing basic strategies and illustrating them by means of a back-testing with real financial market data. The sources of the data used in this project are Bloomberg for intraday time series and Yahoo! for daily prices. All numeric computations and graphics used and shown in this project were implemented in MATLAB^R scratch from scratch as a part of this thesis. No other mathematical or statistical software was used.

Epidemic models with an infected-infectious period

The introduction of an infective-infectious period on the geographic spread of epidemics is considered in two different models. The classical evolution equations arising in the literature are generalized and the existence of epidemic wave fronts is revised. The asymptotic speed is obtained and improves previous results for the Black Death plague

Avoiding transcription factor competition at promoter level increases the chances of obtaining oscillations

Background: The ultimate goal of synthetic biology is the conception and construction of genetic circuits that are reliable with respect to their designed function (e.g. oscillators, switches). This task remains still to be attained due to the inherent synergy of the biological building blocks and to an insufficient feedback between experiments and mathematical models. Nevertheless, the progress in these directions has been substantial. Results: It has been emphasized in the literature that the architecture of a genetic oscillator must include positive (activating) and negative (inhibiting) genetic interactions in order to yield robust oscillations. Our results point out that the oscillatory capacity is not only affected by the interaction polarity but by how it is implemented at promoter level. For a chosen oscillator architecture, we show by means of numerical simulations that the existence or lack of competition between activator and inhibitor at promoter level affects the probability of producing oscillations and also leaves characteristic fingerprints on the associated period/amplitude features. Conclusions: In comparison with non-competitive binding at promoters, competition drastically reduces the region of the parameters space characterized by oscillatory solutions. Moreover, while competition leads to pulse-like oscillations with long-tail distribution in period and amplitude for various parameters or noisy conditions, the non-competitive scenario shows a characteristic frequency and confined amplitude values. Our study also situates the competition mechanism in the context of existing genetic oscillators, with emphasis on the Atkinson oscillator.

Driver scheduling problem modelling

The Drivers Scheduling Problem (DSP) consists of selecting a set of duties for vehicle drivers, for example buses, trains, plane or boat drivers or pilots, for the transportation of passengers or goods. This is a complex problem because it involves several constraints related to labour and company rules and can also present different evaluation criteria and objectives. Being able to develop an adequate model for this problem that can represent the real problem as close as possible is an important research area.The main objective of this research work is to present new mathematical models to the DSP problem that represent all the complexity of the drivers scheduling problem, and also demonstrate that the solutions of these models can be easily implemented in real situations. This issue has been recognized by several authors and as important problem in Public Transportation. The most well-known and general formulation for the DSP is a Set Partition/Set Covering Model (SPP/SCP). However, to a large extend these models simplify some of the specific business aspects and issues of real problems. This makes it difficult to use these models as automatic planning systems because the schedules obtained must be modified manually to be implemented in real situations. Based on extensive passenger transportation experience in bus companies in Portugal, we propose new alternative models to formulate the DSP problem. These models are also based on Set Partitioning/Covering Models; however, they take into account the bus operator issues and the perspective opinions and environment of the user.We follow the steps of the Operations Research Methodology which consist of: Identify the Problem; Understand the System; Formulate a Mathematical Model; Verify the Model; Select the Best Alternative; Present the Results of theAnalysis and Implement and Evaluate. All the processes are done with close participation and involvement of the final users from different transportation companies. The planner s opinion and main criticisms are used to improve the proposed model in a continuous enrichment process. The final objective is to have a model that can be incorporated into an information system to be used as an automatic tool to produce driver schedules. Therefore, the criteria for evaluating the models is the capacity to generate real and useful schedules that can be implemented without many manual adjustments or modifications. We have considered the following as measures of the quality of the model: simplicity, solution quality and applicability. We tested the alternative models with a set of real data obtained from several different transportation companies and analyzed the optimal schedules obtained with respect to the applicability of the solution to the real situation. To do this, the schedules were analyzed by the planners to determine their quality and applicability. The main result of this work is the proposition of new mathematical models for the DSP that better represent the realities of the passenger transportation operators and lead to better schedules that can be implemented directly in real situations.

Exactly solvable phase oscillator models with syncrhonization dynamics

Populations of phase oscillators interacting globally through a general coupling function f(x) have been considered. We analyze the conditions required to ensure the existence of a Lyapunov functional giving close expressions for it in terms of a generating function. We have also proposed a family of exactly solvable models with singular couplings showing that it is possible to map the synchronization phenomenon into other physical problems. In particular, the stationary solutions of the least singular coupling considered, f(x) = sgn(x), have been found analytically in terms of elliptic functions. This last case is one of the few nontrivial models for synchronization dynamics which can be analytically solved.

Freeway Work Zone Lane Capacity, TPF-5(081), 2009

The focus of this report is a capacity analysis of two long-term urban freeway Work Zones. Work Zone #1 tapered four mainline lanes to two, using two separate tapers; Work Zone #2 tapered two mainline lanes to one. Work Zone throughput was analyzed throughout the day over multiple days and traffic operations conditions were analyzed up to a distance of five miles upstream of the Work Zone entrance. Historical data from pavement-embedded detectors were used to analyze traffic conditions. The database consisted of five-minute volume, speed and occupancy data collected from 78 detectors for a total of 50 days. Congestion during each analyzed Work Zone existed for more than fourteen hours each day; Work Zone impacts adversely affected freeway operations over distances of 3.7 to 4.2 miles. Speed and occupancy conditions further upstream were, however, not affected, or even improved due to significant trip diversion. Work Zone capacity was defined based on the maximum traffic flows observed over a one-hour period; throughput values were also compiled over longer periods of time when traffic was within 90% of the maximum observed one-hour flows, as well as over the multi-hour mid-day period. The Highway Capacity Manual freeway capacity definition based on the maximum observed 15-min period was not used, since it would have no practical application in estimating Work Zone throughput when congested conditions prevail for the majority of the hours of the day. Certain noteworthy changes took place for the duration of the analyzed Work Zones: per-lane throughput dropped; morning peak periods started earlier, evening peak periods ended later and lasted longer; mid-day volumes dropped accompanied by the highest occupancies of the day. Trip diversion was evident in lower volumes entering the analyzed freeway corridor, higher volumes using off-ramps and lower volumes using onramps upstream of the Work Zones. The majority of diverted traffic comprised smaller vehicles (vehicles up to 21 feet in length); combination truck volumes increased and their use of the median lane increased, contrary to smaller vehicles that shifted toward a heavier use of the shoulder lane.

Missouri Work Zone Capacity: Results of Field Data Analysis, TPF-5(081), 2011

This report presents the results of work zone field data analyzed on interstate highways in Missouri to determine the mean breakdown and queue-discharge flow rates as measures of capacity. Several days of traffic data collected at a work zone near Pacific, Missouri with a speed limit of 50 mph were analyzed in both the eastbound and westbound directions. As a result, a total of eleven breakdown events were identified using average speed profiles. The traffic flows prior to and after the onset of congestion were studied. Breakdown flow rates ranged between 1194 to 1404 vphpl, with an average of 1295 vphpl, and a mean queue discharge rate of 1072 vphpl was determined. Mean queue discharge, as used by the Highway Capacity Manual 2000 (HCM), in terms of pcphpl was found to be 1199, well below the HCM’s average capacity of 1600 pcphpl. This reduced capacity found at the site is attributable mainly to narrower lane width and higher percentage of heavy vehicles, around 25%, in the traffic stream. The difference found between mean breakdown flow (1295 vphpl) and queue-discharge flow (1072 vphpl) has been observed widely, and is due to reduced traffic flow once traffic breaks down and queues start to form. The Missouri DOT currently uses a spreadsheet for work zone planning applications that assumes the same values of breakdown and mean queue discharge flow rates. This study proposes that breakdown flow rates should be used to forecast the onset of congestion, whereas mean queue discharge flow rates should be used to estimate delays under congested conditions. Hence, it is recommended that the spreadsheet be refined accordingly.

Dynamic models of viral replication and latency.

PURPOSE OF REVIEW: HIV targets primary CD4(+) T cells. The virus depends on the physiological state of its target cells for efficient replication, and, in turn, viral infection perturbs the cellular state significantly. Identifying the virus-host interactions that drive these dynamic changes is important for a better understanding of viral pathogenesis and persistence. The present review focuses on experimental and computational approaches to study the dynamics of viral replication and latency. RECENT FINDINGS: It was recently shown that only a fraction of the inducible latently infected reservoirs are successfully induced upon stimulation in ex-vivo models while additional rounds of stimulation make allowance for reactivation of more latently infected cells. This highlights the potential role of treatment duration and timing as important factors for successful reactivation of latently infected cells. The dynamics of HIV productive infection and latency have been investigated using transcriptome and proteome data. The cellular activation state has shown to be a major determinant of viral reactivation success. Mathematical models of latency have been used to explore the dynamics of the latent viral reservoir decay. SUMMARY: Timing is an important component of biological interactions. Temporal analyses covering aspects of viral life cycle are essential for gathering a comprehensive picture of HIV interaction with the host cell and untangling the complexity of latency. Understanding the dynamic changes tipping the balance between success and failure of HIV particle production might be key to eradicate the viral reservoir.

Application of theoretical models to financial innovation

The present paper is aimed at providing a general strategic overview of the existing theoretical models that have applications in the field of financial innovation. Whereas most financialdevelopments have relied upon traditional economic tools, a new stream of research is defining a novel paradigm in which mathematical models from diverse scientific disciplines are being applied to conceptualize and explain economic and financial behavior. Indeed, terms such as ‘econophysics’ or ‘quantum finance’ have recently appeared to embrace efforts in this direction. As a first contact with such research, the project will present a brief description of some of the main theoretical models that have applications in finance and economics, and will try to present, if possible, potential new applications to particular areas in financial analysis, or new applicable models. As a result, emphasiswill be put on the implications of this research for the financial sector and its future dynamics.

Monte Carlo Methods in Parameter Estimation of Nonlinear Models

Yksi keskeisimmistä tehtävistä matemaattisten mallien tilastollisessa analyysissä on mallien tuntemattomien parametrien estimointi. Tässä diplomityössä ollaan kiinnostuneita tuntemattomien parametrien jakaumista ja niiden muodostamiseen sopivista numeerisista menetelmistä, etenkin tapauksissa, joissa malli on epälineaarinen parametrien suhteen. Erilaisten numeeristen menetelmien osalta pääpaino on Markovin ketju Monte Carlo -menetelmissä (MCMC). Nämä laskentaintensiiviset menetelmät ovat viime aikoina kasvattaneet suosiotaan lähinnä kasvaneen laskentatehon vuoksi. Sekä Markovin ketjujen että Monte Carlo -simuloinnin teoriaa on esitelty työssä siinä määrin, että menetelmien toimivuus saadaan perusteltua. Viime aikoina kehitetyistä menetelmistä tarkastellaan etenkin adaptiivisia MCMC menetelmiä. Työn lähestymistapa on käytännönläheinen ja erilaisia MCMC -menetelmien toteutukseen liittyviä asioita korostetaan. Työn empiirisessä osuudessa tarkastellaan viiden esimerkkimallin tuntemattomien parametrien jakaumaa käyttäen hyväksi teoriaosassa esitettyjä menetelmiä. Mallit kuvaavat kemiallisia reaktioita ja kuvataan tavallisina differentiaaliyhtälöryhminä. Mallit on kerätty kemisteiltä Lappeenrannan teknillisestä yliopistosta ja Åbo Akademista, Turusta.

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.