CYR61 and TAZ Upregulation and Focal Epithelial to Mesenchymal Transition May Be Early Predictors of Barrett’s Esophagus Malignant Progression

Barrett's esophagus is the major risk factor for esophageal adenocarcinoma. It has a low but non-neglectable risk, high surveillance costs and no reliable risk stratification markers. We sought to identify early biomarkers, predictive of Barrett's malignant progression, using a meta-analysis approach on gene expression data. This in silico strategy was followed by experimental validation in a cohort of patients with extended follow up from the Instituto Português de Oncologia de Lisboa de Francisco Gentil EPE (Portugal). Bioinformatics and systems biology approaches singled out two candidate predictive markers for Barrett's progression, CYR61 and TAZ. Although previously implicated in other malignancies and in epithelial-to-mesenchymal transition phenotypes, our experimental validation shows for the first time that CYR61 and TAZ have the potential to be predictive biomarkers for cancer progression. Experimental validation by reverse transcriptase quantitative PCR and immunohistochemistry confirmed the up-regulation of both genes in Barrett's samples associated with high-grade dysplasia/adenocarcinoma. In our cohort CYR61 and TAZ up-regulation ranged from one to ten years prior to progression to adenocarcinoma in Barrett's esophagus index samples. Finally, we found that CYR61 and TAZ over-expression is correlated with early focal signs of epithelial to mesenchymal transition. Our results highlight both CYR61 and TAZ genes as potential predictive biomarkers for stratification of the risk for development of adenocarcinoma and suggest a potential mechanistic route for Barrett's esophagus neoplastic progression.

Développement, validation et nouvelles applications d’un modèle d’analyse des modes normaux basé sur la séquence et la structure de protéines

Les protéines existent sous différents états fonctionnels régulés de façon précise par leur environnement afin de maintenir l‘homéostasie de la cellule et de l‘organisme vivant. La prévalence de ces états protéiques est dictée par leur énergie libre de Gibbs alors que la vitesse de transition entre ces états biologiquement pertinents est déterminée par le paysage d‘énergie libre. Ces paramètres sont particulièrement intéressants dans un contexte thérapeutique et biotechnologique, où leur perturbation par la modulation de la séquence protéique par des mutations affecte leur fonction. Bien que des nouvelles approches expérimentales permettent d‘étudier l‘effet de mutations en haut débit pour une protéine, ces méthodes sont laborieuses et ne couvrent qu‘une fraction de l‘ensemble des structures primaires d‘intérêt. L‘utilisation de modèles bio-informatiques permet de tester et générer in silico différentes hypothèses afin d‘orienter les approches expérimentales. Cependant, ces méthodes basées sur la structure se concentrent principalement sur la prédiction de l‘enthalpie d‘un état, alors que plusieurs évidences expérimentales ont démontré l‘importance de la contribution de l‘entropie. De plus, ces approches ignorent l‘importance de l‘espace conformationnel protéique dicté par le paysage énergétique cruciale à son fonctionnement. Une analyse des modes normaux peut être effectuée afin d‘explorer cet espace par l‘approximation que la protéine est dans une conformation d‘équilibre où chaque acide aminé est représenté par une masse régie par un potentiel harmonique. Les approches actuelles ignorent l‘identité des résidus et ne peuvent prédire l‘effet de mutations sur les propriétés dynamiques. Nous avons développé un nouveau modèle appelé ENCoM qui pallie à cette lacune en intégrant de l‘information physique et spécifique sur les contacts entre les atomes des chaînes latérales. Cet ajout permet une meilleure description de changements conformationnels d‘enzymes, la prédiction de l‘effet d‘une mutation allostérique dans la protéine DHFR et également la prédiction de l‘effet de mutations sur la stabilité protéique par une valeur entropique. Comparativement à des approches spécifiquement développées pour cette application, ENCoM est plus constant et prédit mieux l‘effet de mutations stabilisantes. Notre approche a également été en mesure de capturer la pression évolutive qui confère aux protéines d‘organismes thermophiles une thermorésistance accrue.

Bioinformatics analysis of the Pedobacter sp. NL19 genome

The last decades of the 20th century defined the genetic engineering advent, climaxing in the development of techniques, such as PCR and Sanger sequencing. This, permitted the appearance of new techniques to sequencing whole genomes, identified as next-generation sequencing. One of the many applications of these techniques is the in silico search for new secondary metabolites, synthesized by microorganisms exhibiting antimicrobial properties. The peptide antibiotics compounds can be classified in two classes, according to their biosynthesis, in ribosomal or nonribosomal peptides. Lanthipeptides are the most studied ribosomal peptides and are characterized by the presence of lanthionine and methylanthionine that result from posttranslational modifications. Lanthipeptides are divided in four classes, depending on their biosynthetic machinery. In class I, a LanB enzyme dehydrate serine and threonine residues in the C-terminus precursor peptide. Then, these residues undergo a cyclization step performed by a LanC enzyme, forming the lanthionine rings. The cleavage and the transport of the peptide is achieved by the LanP and LanT enzymes, respectively. Although, in class II only one enzyme, LanM, is responsible for the dehydration and cyclization steps and also only one enzyme performs the cleavage and transport, LanT. Pedobacter sp. NL19 is a Gram-negative bacterium, isolated from sludge of an abandon uranium mine, in Viseu (Portugal). Antibacterial activity in vitro was detected against several Gram-positive and Gram-negative bacteria. Sequencing and in silico analysis of NL19 genome revealed the presence of 21 biosynthetic clusters for secondary metabolites, including nonribosomal and ribosomal peptides biosynthetic clusters. Four lanthipeptides clusters were predicted, comprising the precursor peptides, the modifying enzymes (LanB and LanC), and also a bifunctional LanT. This result revealed the hybrid nature of the clusters, comprising characteristics from two distinct classes, which are poorly described in literature. The phylogenetic analysis of their enzymes showed that they clustered within the bacteroidetes clade. Furthermore, hybrid gene clusters were also found in other species of this phylum, revealing that it is a common characteristic in this group. Finally, the analysis of NL19 colonies by MALDI-TOF MS allowed the identification of a 3180 Da mass that corresponds to the predicted mass of a lanthipeptide encoded in one of the clusters. However, this result is not fully conclusive and further experiments are needed to understand the full potential of the compounds encoded in this type of clusters. In conclusion, it was determined that NL19 strain has the potential to produce diverse secondary metabolites, including lanthipeptides that were not functionally characterized so far.

The Repellent DEET Potentiates Carbamate Effects via Insect Muscarinic Receptor Interactions: An Alternative Strategy to Control Insect Vector-Borne Diseases.

International audience

Complementation of the embryo-lethal T-DNA insertion mutant of AUXIN-BINDING-PROTEIN 1 (ABP1) with abp1 point mutated versions reveals crosstalk of ABP1 and phytochromes

The function of the extracytoplasmic AUXIN-BINDING-PROTEIN1 (ABP1) is largely enigmatic. We complemented a homozygous T-DNA insertion null mutant of ABP1 in Arabidopsis thaliana Wassilewskia with three mutated and one wild-type (wt) ABP1 cDNA, all tagged C-terminally with a strepII-FLAG tag upstream the KDEL signal. Based on in silico modelling, the abp1 mutants were predicted to have altered geometries of the auxin binding pocket and calculated auxin binding energies lower than the wt. Phenotypes linked to auxin transport were compromised in these three complemented abp1 mutants. Red light effects, such as elongation of hypocotyls in constant red (R) and far-red (FR) light, in white light supplemented by FR light simulating shade, and inhibition of gravitropism by R or FR, were all compromised in the complemented lines. Using auxin-or light-induced expression of marker genes, we showed that auxininduced expression was delayed already after 10 min, and light-induced expression within 60 min, even though TIR1/AFB or phyB are thought to act as receptors relevant for gene expression regulation. The expression of marker genes in seedlings responding to both auxin and shade showed that for both stimuli regulation of marker gene expression was altered after 10-20 min in the wild type and phyB mutant. The rapidity of expression responses provides a framework for the mechanics of functional interaction of ABP1 and phyB to trigger interwoven signalling pathways.

Modelling Immunological Memory

Accurate immunological models offer the possibility of performing highthroughput experiments in silico that can predict, or at least suggest, in vivo phenomena. In this chapter, we compare various models of immunological memory. We first validate an experimental immunological simulator, developed by the authors, by simulating several theories of immunological memory with known results. We then use the same system to evaluate the predicted effects of a theory of immunological memory. The resulting model has not been explored before in artificial immune systems research, and we compare the simulated in silico output with in vivo measurements. Although the theory appears valid, we suggest that there are a common set of reasons why immunological memory models are a useful support tool; not conclusive in themselves.

Multi-body simulation of a canine hind limb: model development, experimental validation and calculation of ground reaction forces

Background: Among other causes the long-term result of hip prostheses in dogs is determined by aseptic loosening. A prevention of prosthesis complications can be achieved by an optimization of the tribological system which finally results in improved implant duration. In this context a computerized model for the calculation of hip joint loadings during different motions would be of benefit. In a first step in the development of such an inverse dynamic multi-body simulation (MBS-) model we here present the setup of a canine hind limb model applicable for the calculation of ground reaction forces. Methods: The anatomical geometries of the MBS-model have been established using computer tomography- (CT-) and magnetic resonance imaging- (MRI-) data. The CT-data were collected from the pelvis, femora, tibiae and pads of a mixed-breed adult dog. Geometric information about 22 muscles of the pelvic extremity of 4 mixed-breed adult dogs was determined using MRI. Kinematic and kinetic data obtained by motion analysis of a clinically healthy dog during a gait cycle (1 m/s) on an instrumented treadmill were used to drive the model in the multi-body simulation. Results and Discussion: As a result the vertical ground reaction forces (z-direction) calculated by the MBS-system show a maximum deviation of 1.75%BW for the left and 4.65%BW for the right hind limb from the treadmill measurements. The calculated peak ground reaction forces in z- and y-direction were found to be comparable to the treadmill measurements, whereas the curve characteristics of the forces in y-direction were not in complete alignment. Conclusion: In conclusion, it could be demonstrated that the developed MBS-model is suitable for simulating ground reaction forces of dogs during walking. In forthcoming investigations the model will be developed further for the calculation of forces and moments acting on the hip joint during different movements, which can be of help in context with the in silico development and testing of hip prostheses.

"Mutações no gene OCT4 em células eritroleucêmicas humanas: implicações sobre o fenótipo de resistência a múltiplas drogas (MDR)”

O fator de transcrição OCT4 é um importante marcador de células tronco e tem sido relacionado com o conceito de células tronco tumorais (CTTs). Recentemente, ele tem sido também relacionado ao fenótipo de resistência a múltiplas drogas (MDR). O objetivo deste trabalho foi testar se o OCT4 está ligado ao fenótipo MDR em células eritroleucêmicas derivadas da linhagem LMC-K562. Para isso, foi realizada a análise de expressão de genes associados à superfamília de transportadores ABC (ATP Binding Cassette) e, também, de fatores de transcrição relacionados a células tronco. Os primeiros resultados apontaram uma relação direta entre OCT4 e transportadores ABC na linhagem MDR derivada de K562 (Lucena). O sequenciamento de promotores ABC não revelou qualquer mutação que pudesse explicar a expressão diferenciada do OCT4 na linhagem MDR. Posteriormente, o sequenciamento da região contendo o domínio homeobox do gene OCT4 de ambas as linhagens celulares evidenciou, pela primeira vez, que este fator de transcrição é alvo de mutações que podem estar relacionados com o fenótipo MDR. As mutações encontradas implicam substituições de vários aminoácidos em ambas as linhagens celulares. K562 teve sete substituições (três exclusivas), enquanto Lucena teve 13 (nove exclusivas). Além disso, um busca in silico por motivos de fosforilação dentro da sequência de aminoácidos comparada, revelou que o OCT4 humano normal possui sete motivos potenciais de fosforilação. Entretanto, K562 perdeu um motivo e Lucena dois deles. Moléculas com diferentes padrões de fosforilação podem ter sua função modificada. Estes resultados indicam o OCT4 com uma alvo potencial para o tratamento do câncer, especialmente aqueles resistentes à quimioterapia.

Improving the use of G-CSF during chemotherapy using physiological mathematical modelling : a quantitative systems pharmacology approach

La diminution des doses administrées ou même la cessation complète d'un traitement chimiothérapeutique est souvent la conséquence de la réduction du nombre de neutrophiles, qui sont les globules blancs les plus fréquents dans le sang. Cette réduction dans le nombre absolu des neutrophiles, aussi connue sous le nom de myélosuppression, est précipitée par les effets létaux non spécifiques des médicaments anti-cancéreux, qui, parallèlement à leur effet thérapeutique, produisent aussi des effets toxiques sur les cellules saines. Dans le but d'atténuer cet impact myélosuppresseur, on administre aux patients un facteur de stimulation des colonies de granulocytes recombinant humain (rhG-CSF), une forme exogène du G-CSF, l'hormone responsable de la stimulation de la production des neutrophiles et de leurs libération dans la circulation sanguine. Bien que les bienfaits d'un traitement prophylactique avec le G-CSF pendant la chimiothérapie soient bien établis, les protocoles d'administration demeurent mal définis et sont fréquemment déterminés ad libitum par les cliniciens. Avec l'optique d'améliorer le dosage thérapeutique et rationaliser l'utilisation du rhG-CSF pendant le traitement chimiothérapeutique, nous avons développé un modèle physiologique du processus de granulopoïèse, qui incorpore les connaissances actuelles de pointe relatives à la production des neutrophiles des cellules souches hématopoïétiques dans la moelle osseuse. À ce modèle physiologique, nous avons intégré des modèles pharmacocinétiques/pharmacodynamiques (PK/PD) de deux médicaments: le PM00104 (Zalypsis®), un médicament anti-cancéreux, et le rhG-CSF (filgrastim). En se servant des principes fondamentaux sous-jacents à la physiologie, nous avons estimé les paramètres de manière exhaustive sans devoir recourir à l'ajustement des données, ce qui nous a permis de prédire des données cliniques provenant de 172 patients soumis au protocol CHOP14 (6 cycles de chimiothérapie avec une période de 14 jours où l'administration du rhG-CSF se fait du jour 4 au jour 13 post-chimiothérapie). En utilisant ce modèle physio-PK/PD, nous avons démontré que le nombre d'administrations du rhG-CSF pourrait être réduit de dix (pratique actuelle) à quatre ou même trois administrations, à condition de retarder le début du traitement prophylactique par le rhG-CSF. Dans un souci d'applicabilité clinique de notre approche de modélisation, nous avons investigué l'impact de la variabilité PK présente dans une population de patients, sur les prédictions du modèle, en intégrant des modèles PK de population (Pop-PK) des deux médicaments. En considérant des cohortes de 500 patients in silico pour chacun des cinq scénarios de variabilité plausibles et en utilisant trois marqueurs cliniques, soient le temps au nadir des neutrophiles, la valeur du nadir, ainsi que l'aire sous la courbe concentration-effet, nous avons établi qu'il n'y avait aucune différence significative dans les prédictions du modèle entre le patient-type et la population. Ceci démontre la robustesse de l'approche que nous avons développée et qui s'apparente à une approche de pharmacologie quantitative des systèmes (QSP). Motivés par l'utilisation du rhG-CSF dans le traitement d'autres maladies, comme des pathologies périodiques telles que la neutropénie cyclique, nous avons ensuite soumis l'étude du modèle au contexte des maladies dynamiques. En mettant en évidence la non validité du paradigme de la rétroaction des cytokines pour l'administration exogène des mimétiques du G-CSF, nous avons développé un modèle physiologique PK/PD novateur comprenant les concentrations libres et liées du G-CSF. Ce nouveau modèle PK a aussi nécessité des changements dans le modèle PD puisqu’il nous a permis de retracer les concentrations du G-CSF lié aux neutrophiles. Nous avons démontré que l'hypothèse sous-jacente de l'équilibre entre la concentration libre et liée, selon la loi d'action de masse, n'est plus valide pour le G-CSF aux concentrations endogènes et mènerait en fait à la surestimation de la clairance rénale du médicament. En procédant ainsi, nous avons réussi à reproduire des données cliniques obtenues dans diverses conditions (l'administration exogène du G-CSF, l'administration du PM00104, CHOP14). Nous avons aussi fourni une explication logique des mécanismes responsables de la réponse physiologique aux deux médicaments. Finalement, afin de mettre en exergue l’approche intégrative en pharmacologie adoptée dans cette thèse, nous avons démontré sa valeur inestimable pour la mise en lumière et la reconstruction des systèmes vivants complexes, en faisant le parallèle avec d’autres disciplines scientifiques telles que la paléontologie et la forensique, où une approche semblable a largement fait ses preuves. Nous avons aussi discuté du potentiel de la pharmacologie quantitative des systèmes appliquées au développement du médicament et à la médecine translationnelle, en se servant du modèle physio-PK/PD que nous avons mis au point.

Control HSMC para regulación de glucosa sanguínea

Este artículo propone una nueva estrategia de control basada en medidas continuas de glucosa y un controlador por modo deslizante que se habitúa (HSMC). El HSMC es desarrollado, combinando la ley de control por modo deslizante y los principios de control por habituación. El HSMC aplicado a la regulación de glucosa sanguínea en la unidad de cuidados intensivos, incluye tanto entrada de glucosa, como de infusión de insulina intravasculares a fin de proveer el suministro de nutrición y mejorar el rechazo a la perturbación. El estudio basado en simulaciones (in silico), usando un modelo fisiológico de la dinámica glucosa-insulina, muestra que la estrategia de control propuesta funciona apropiadamente. Finalmente, se compara el desempeño del controlador propuesto con respecto a un controlador PID estándar.

Caracterização e comparação de genes expressos em ápices meristemáticos de cana-de-açúcar cultivados em SP e RN

In this work, we used sugarcane as a model due to its importance for sugar and ethanol production. Unlike the current plant models, sugarcane presents a complex genetics and an enormous allelic variation. Here, we report the analysis of SAGE libraries produced using the shoot apical meristem from contrasted genotypes by flowering induction (non-flowering vs. early-flowering varieties) grown under São Paulo state conditions. The expression pattern was analyzed using samples from São Paulo (SP) and Rio Grande do Norte (RN) states. These results showed that cDNAs identified by SAGE libraries had differential expression only in São Paulo state samples. Furthermore, the cDNA identified CYP (Citocrome P450) was chosen for in silico and genome characterization because it was found in SAGE libraries and subtractive libraries from samples from RN. Phylogenetic trees showed the relationship for these sequences. Furthermore, the qRT-PCR for CYP showed a potential role as flowering indutor for RN samples considering different isophorms. Considering the results present here, it can be consider that CYP gene may be used as molecular marker

Polimorfismos nos genes ERBB2 e PARK2/PACRG e sua associação com a hanseníase em populações do Rio Grande do Norte e Pará

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior

Caracterização funcional de dois cDNAs homólogos à ciclofilina e à proteína reprimida por auxina (ARP) identificados em bibliotecas subtrativas a para floração em tomateiro

Flowering is a fundamental process in the life cycle for plant. This process is marked by vegetative to reproductive apical meristem conversion, due to interactions between several factors, both internal and external to plant. Therefore, eight subtractive libraries were constructed using apical meristem induced or not induced for two contrasting species: Solanum lycopersicum cv. Micro-Tom and Solanum pimpinellifolium. Several cDNAs were identified and among these, were selected two cDNAs: one homologous cDNA to cyclophilin (LeCYP1) and the other to Auxin repressed protein (ARP). It has observed that LeCYP1 and ARP genes are important in the developmental process to plants. In silico analysis, were used several databases with the exclusion criterion E-value <1.0x10-15. As a result, conservation was observed for proteins analyzed by means of multiple alignments and the presence of functional domains. Then, overexpression cassettes were constructed for the ARP cDNA in sense and antisense orientations. For this step, it was used the CaMV35S promoter. The cDNA orientation (sense or antisense) in relation to the promoter was determined by restriction enzymes and sequencing. Then, this cassette was transferred to binary vector pZP211 and these cassettes were transferred into Agrobacterium tumefaciens LBA4404. S. lycopersicum cv. Micro-Tom (MT) and MT-Rg1 plants were transformed. In addition, seedlings were subjected to hormone treatments using a synthetic auxin (- naphthalene acetic acid) and cyclosporin A (cyclophilin inhibitor) treatments and it was found that the hormone treatment there were changes in development of lateral roots pattern, probably related to decreases in auxin signaling caused by reduction of LeCYP1 in MT-dgt plants while cyclosporin A treatments, there was a slight delay in flowering in cv. MT plants. Furthermore, assay with real-time PCR (RT-qPCR) were done for expression level analysis from LeCYP1 and ARP in order to functionally characterize these sequences in tomato plants.

Análise filogenética e funcional de dois genes de reparo homólogos a AP endonuclease em cana-de-açúcar: ScARP1 e ScARP3

The genome of all organisms constantly suffers the influence of mutagenic factors from endogenous and/or exogenous origin, which may result in damage for the genome. In order to keep the genome integrity there are different DNA repair pathway to detect and correct these lesions. In relation to the plants as being sessile organisms, they are exposed to this damage frequently. The Base Excision DNA Repair (BER) is responsible to detect and repair oxidative lesions. Previous work in sugarcane identified two sequences that were homologous to Arabidopsis thaliana: ScARP1 ScARP3. These two sequences were homologous to AP endonuclease from BER pathway. Then, the aim of this work was to characterize these two sequence using different approaches: phylogenetic analysis, in silico protein organelle localization and by Nicotiana tabacum transgenic plants with overexpression cassette. The in silico data obtained showed a duplication of this sequence in sugarcane and Poaceae probably by a WGD event. Furthermore, in silico analysis showed a new localization in nuclei for ScARP1 protein. The data obtained with transgenic plants showed a change in development and morphology. Transgenic plants had slow development when compared to plants not transformed. Then, these results allowed us to understand better the potential role of this sequence in sugarcane and in plants in general. More work is important to be done in order to confirm the protein localization and protein characterization for ScARP1 and ScARP3