Structuring research methods and data with the research object model: genomics workflows as a case study

Background: One of the main challenges for biomedical research lies in the computer-assisted integrative study of large and increasingly complex combinations of data in order to understand molecular mechanisms. The preservation of the materials and methods of such computational experiments with clear annotations is essential for understanding an experiment, and this is increasingly recognized in the bioinformatics community. Our assumption is that offering means of digital, structured aggregation and annotation of the objects of an experiment will provide necessary meta-data for a scientist to understand and recreate the results of an experiment. To support this we explored a model for the semantic description of a workflow-centric Research Object (RO), where an RO is defined as a resource that aggregates other resources, e.g., datasets, software, spreadsheets, text, etc. We applied this model to a case study where we analysed human metabolite variation by workflows. Results: We present the application of the workflow-centric RO model for our bioinformatics case study. Three workflows were produced following recently defined Best Practices for workflow design. By modelling the experiment as an RO, we were able to automatically query the experiment and answer questions such as “which particular data was input to a particular workflow to test a particular hypothesis?”, and “which particular conclusions were drawn from a particular workflow?”. Conclusions: Applying a workflow-centric RO model to aggregate and annotate the resources used in a bioinformatics experiment, allowed us to retrieve the conclusions of the experiment in the context of the driving hypothesis, the executed workflows and their input data. The RO model is an extendable reference model that can be used by other systems as well.

Combinatorial interaction between AtbZIP44 and other bZIP transcription factors in the regulation of Arabidopsis thaliana seed maturation and germination

El gen AtbZIP44 se expresa en el endospermo micropilar y en la radícula de semillas en germinación. El mRNA del gen AtbZIP9 se localiza también en la radícula. Se realizan experimentos de interacción proteína-proteína entre los correspondientes factores transcripcionales.

Auscultación de presas de hormigón con láser escáner terrestre y modelización de las deformaciones del paramento

Se está produciendo en la geodesia un cambio de paradigma en la concepción de los modelos digitales del terreno, pasando de diseñar el modelo con el menor número de puntos posibles a hacerlo con cientos de miles o millones de puntos. Este cambio ha sido consecuencia de la introducción de nuevas tecnologías como el escáner láser, la interferometría radar y el tratamiento de imágenes. La rápida aceptación de estas nuevas tecnologías se debe principalmente a la gran velocidad en la toma de datos, a la accesibilidad por no precisar de prisma y al alto grado de detalle de los modelos. Los métodos topográficos clásicos se basan en medidas discretas de puntos que considerados en su conjunto forman un modelo; su precisión se deriva de la precisión en la toma singular de estos puntos. La tecnología láser escáner terrestre (TLS) supone una aproximación diferente para la generación del modelo del objeto observado. Las nubes de puntos, producto del escaneo con TLS, pasan a ser tratadas en su conjunto mediante análisis de áreas, de forma que ahora el modelo final no es el resultado de una agregación de puntos sino la de la mejor superficie que se adapta a las nubes de puntos. Al comparar precisiones en la captura de puntos singulares realizados con métodos taquimétricos y equipos TLS la inferioridad de estos últimos es clara; sin embargo es en el tratamiento de las nubes de puntos, con los métodos de análisis basados en áreas, se han obtenido precisiones aceptables y se ha podido considerar plenamente la incorporación de esta tecnología en estudios de deformaciones y movimientos de estructuras. Entre las aplicaciones del TLS destacan las de registro del patrimonio, registro de las fases en la construcción de plantas industriales y estructuras, atestados de accidentes y monitorización de movimientos del terreno y deformaciones de estructuras. En la auscultación de presas, comparado con la monitorización de puntos concretos dentro, en coronación o en el paramento de la presa, disponer de un modelo continuo del paramento aguas abajo de la presa abre la posibilidad de introducir los métodos de análisis de deformaciones de superficies y la creación de modelos de comportamiento que mejoren la comprensión y previsión de sus movimientos. No obstante, la aplicación de la tecnología TLS en la auscultación de presas debe considerarse como un método complementario a los existentes. Mientras que los péndulos y la reciente técnica basada en el sistema de posicionamiento global diferencial (DGPS) dan una información continua de los movimientos de determinados puntos de la presa, el TLS permite ver la evolución estacional y detectar posibles zonas problemáticas en todo el paramento. En este trabajo se analizan las características de la tecnología TLS y los parámetros que intervienen en la precisión final de los escaneos. Se constata la necesidad de utilizar equipos basados en la medida directa del tiempo de vuelo, también llamados pulsados, para distancias entre 100 m y 300 m Se estudia la aplicación del TLS a la modelización de estructuras y paramentos verticales. Se analizan los factores que influyen en la precisión final, como el registro de nubes, tipo de dianas y el efecto conjunto del ángulo y la distancia de escaneo. Finalmente, se hace una comparación de los movimientos dados por los péndulos directos de una presa con los obtenidos del análisis de las nubes de puntos correspondientes a varias campañas de escaneos de la misma presa. Se propone y valida el empleo de gráficos patrón para relacionar las variables precisión o exactitud con los factores distancia y ángulo de escaneo en el diseño de trabajos de campo. Se expone su aplicación en la preparación del trabajo de campo para la realización de una campaña de escaneos dirigida al control de movimientos de una presa y se realizan recomendaciones para la aplicación de la técnica TLS a grandes estructuras. Se ha elaborado el gráfico patrón de un equipo TLS concreto de alcance medio. Para ello se hicieron dos ensayos de campo en condiciones reales de trabajo, realizando escaneos en todo el rango de distancias y ángulos de escaneo del equipo. Se analizan dos métodos para obtener la precisión en la modelización de paramentos y la detección de movimientos de estos: el método del “plano de mejor ajuste” y el método de la “deformación simulada”. Por último, se presentan los resultados de la comparación de los movimientos estacionales de una presa arco-gravedad entre los registrados con los péndulos directos y los obtenidos a partir de los escaneos realizados con un TLS. Los resultados muestran diferencias de milímetros, siendo el mejor de ellos del orden de un milímetro. Se explica la metodología utilizada y se hacen consideraciones respecto a la densidad de puntos de las nubes y al tamaño de las mallas de triángulos. A shift of paradigm in the conception of the survey digital models is taking place in geodesy, moving from designing a model with the fewer possible number of points to models of hundreds of thousand or million points. This change has happened because of the introduction of new technologies like the laser scanner, the interferometry radar and the processing of images. The fast acceptance of these new technologies has been due mainly to the great speed getting the data, to the accessibility as reflectorless technique, and to the high degree of detail of the models. Classic survey methods are based on discreet measures of points that, considered them as a whole, form a model; the precision of the model is then derived from the precision measuring the single points. The terrestrial laser scanner (TLS) technology supposes a different approach to the model generation of the observed object. Point cloud, the result of a TLS scan, must be treated as a whole, by means of area-based analysis; so, the final model is not an aggregation of points but the one resulting from the best surface that fits with the point cloud. Comparing precisions between the one resulting from the capture of singular points made with tachometric measurement methods and with TLS equipment, the inferiority of this last one is clear; but it is in the treatment of the point clouds, using area-based analysis methods, when acceptable precisions have been obtained and it has been possible to consider the incorporation of this technology for monitoring structures deformations. Among TLS applications it have to be emphasized those of registry of the cultural heritage, stages registry during construction of industrial plants and structures, police statement of accidents and monitorization of land movements and structures deformations. Compared with the classical dam monitoring, approach based on the registry of a set of points, the fact having a continuous model of the downstream face allows the possibility of introducing deformation analysis methods and behavior models that would improve the understanding and forecast of dam movements. However, the application of TLS technology for dam monitoring must be considered like a complementary method with the existing ones. Pendulums and recently the differential global positioning system (DGPS) give a continuous information of the movements of certain points of the dam, whereas TLS allows following its seasonal evolution and to detect damaged zones of the dam. A review of the TLS technology characteristics and the factors affecting the final precision of the scanning data is done. It is stated the need of selecting TLS based on the direct time of flight method, also called pulsed, for scanning distances between 100m and 300m. Modelling of structures and vertical walls is studied. Factors that influence in the final precision, like the registry of point clouds, target types, and the combined effect of scanning distance and angle of incidence are analyzed. Finally, a comparison among the movements given by the direct pendulums of a dam and the ones obtained from the analysis of point clouds is done. A new approach to obtain a complete map-type plot of the precisions of TLS equipment based on the direct measurement of time of flight method at midrange distances is presented. Test were developed in field-like conditions, similar to dam monitoring and other civil engineering works. Taking advantage of graphic semiological techniques, a “distance - angle of incidence” map based was designed and evaluated for field-like conditions. A map-type plot was designed combining isolines with sized and grey scale points, proportional to the precision values they represent. Precisions under different field conditions were compared with specifications. For this purpose, point clouds were evaluated under two approaches: the standar "plane-of-best-fit" and the proposed "simulated deformation”, that showed improved performance. These results lead to a discussion and recommendations about optimal TLS operation in civil engineering works. Finally, results of the comparison of seasonal movements of an arc-gravity dam between the registered by the direct pendulums ant the obtained from the TLS scans, are shown. The results show differences of millimeters, being the best around one millimeter. The used methodology is explained and considerations with respect to the point cloud density and to the size of triangular meshes are done.

Methods and Techniques for the Generation and Eficient Exploitation of RDB2RDF Mappings

RDB to RDF Mapping Language (R2RML) es una recomendación del W3C que permite especificar reglas para transformar bases de datos relacionales a RDF. Estos datos en RDF se pueden materializar y almacenar en un sistema gestor de tripletas RDF (normalmente conocidos con el nombre triple store), en el cual se pueden evaluar consultas SPARQL. Sin embargo, hay casos en los cuales la materialización no es adecuada o posible, por ejemplo, cuando la base de datos se actualiza frecuentemente. En estos casos, lo mejor es considerar los datos en RDF como datos virtuales, de tal manera que las consultas SPARQL anteriormente mencionadas se traduzcan a consultas SQL que se pueden evaluar sobre los sistemas gestores de bases de datos relacionales (SGBD) originales. Para esta traducción se tienen en cuenta los mapeos R2RML. La primera parte de esta tesis se centra en la traducción de consultas. Se propone una formalización de la traducción de SPARQL a SQL utilizando mapeos R2RML. Además se proponen varias técnicas de optimización para generar consultas SQL que son más eficientes cuando son evaluadas en sistemas gestores de bases de datos relacionales. Este enfoque se evalúa mediante un benchmark sintético y varios casos reales. Otra recomendación relacionada con R2RML es la conocida como Direct Mapping (DM), que establece reglas fijas para la transformación de datos relacionales a RDF. A pesar de que ambas recomendaciones se publicaron al mismo tiempo, en septiembre de 2012, todavía no se ha realizado un estudio formal sobre la relación entre ellas. Por tanto, la segunda parte de esta tesis se centra en el estudio de la relación entre R2RML y DM. Se divide este estudio en dos partes: de R2RML a DM, y de DM a R2RML. En el primer caso, se estudia un fragmento de R2RML que tiene la misma expresividad que DM. En el segundo caso, se representan las reglas de DM como mapeos R2RML, y también se añade la semántica implícita (relaciones de subclase, 1-N y M-N) que se puede encontrar codificada en la base de datos. Esta tesis muestra que es posible usar R2RML en casos reales, sin necesidad de realizar materializaciones de los datos, puesto que las consultas SQL generadas son suficientemente eficientes cuando son evaluadas en el sistema gestor de base de datos relacional. Asimismo, esta tesis profundiza en el entendimiento de la relación existente entre las dos recomendaciones del W3C, algo que no había sido estudiado con anterioridad. ABSTRACT. RDB to RDF Mapping Language (R2RML) is a W3C recommendation that allows specifying rules for transforming relational databases into RDF. This RDF data can be materialized and stored in a triple store, so that SPARQL queries can be evaluated by the triple store. However, there are several cases where materialization is not adequate or possible, for example, if the underlying relational database is updated frequently. In those cases, RDF data is better kept virtual, and hence SPARQL queries over it have to be translated into SQL queries to the underlying relational database system considering that the translation process has to take into account the specified R2RML mappings. The first part of this thesis focuses on query translation. We discuss the formalization of the translation from SPARQL to SQL queries that takes into account R2RML mappings. Furthermore, we propose several optimization techniques so that the translation procedure generates SQL queries that can be evaluated more efficiently over the underlying databases. We evaluate our approach using a synthetic benchmark and several real cases, and show positive results that we obtained. Direct Mapping (DM) is another W3C recommendation for the generation of RDF data from relational databases. While R2RML allows users to specify their own transformation rules, DM establishes fixed transformation rules. Although both recommendations were published at the same time, September 2012, there has not been any study regarding the relationship between them. The second part of this thesis focuses on the study of the relationship between R2RML and DM. We divide this study into two directions: from R2RML to DM, and from DM to R2RML. From R2RML to DM, we study a fragment of R2RML having the same expressive power than DM. From DM to R2RML, we represent DM transformation rules as R2RML mappings, and also add the implicit semantics encoded in databases, such as subclass, 1-N and N-N relationships. This thesis shows that by formalizing and optimizing R2RML-based SPARQL to SQL query translation, it is possible to use R2RML engines in real cases as the resulting SQL is efficient enough to be evaluated by the underlying relational databases. In addition to that, this thesis facilitates the understanding of bidirectional relationship between the two W3C recommendations, something that had not been studied before.

Synergistic interaction between leptin and cholecystokinin to reduce short-term food intake in lean mice

Leptin is a circulating protein involved in the long-term regulation of food intake and body weight. Cholecystokinin (CCK) is released postprandially and elicits satiety signals. We investigated the interaction between leptin and CCK-8 in the short-term regulation of food intake induced by 24-hr fasting in lean mice. Leptin, injected intraperitoneally (i.p.) at low doses (4–120 μg/kg), which did not influence feeding behavior for the first 3 hr postinjection, decreased food intake dose dependently by 47–83% during the first hour when coinjected with a subthreshold dose of CCK. Such an interaction was not observed between leptin and bombesin. The food-reducing effect of leptin injected with CCK was not associated with alterations in gastric emptying or locomotor behavior. Leptin–CCK action was blocked by systemic capsaicin at a dose inducing functional ablation of sensory afferent fibers and by devazepide, a CCK-A receptor antagonist but not by the CCK-B receptor antagonist, L-365,260. The decrease in food intake which occurs 5 hr after i.p. injection of leptin alone was also blunted by devazepide. Coinjection of leptin and CCK enhanced the number of Fos-positive cells in the hypothalamic paraventricular nucleus by 60%, whereas leptin or CCK alone did not modify Fos expression. These results indicate the existence of a functional synergistic interaction between leptin and CCK leading to early suppression of food intake which involves CCK-A receptors and capsaicin-sensitive afferent fibers.

Structural details of an interaction between cardiolipin and an integral membrane protein

Anionic lipids play a variety of key roles in biomembrane function, including providing the immediate environment for the integral membrane proteins that catalyze photosynthetic and respiratory energy transduction. Little is known about the molecular basis of these lipid–protein interactions. In this study, x-ray crystallography has been used to examine the structural details of an interaction between cardiolipin and the photoreaction center, a key light-driven electron transfer protein complex found in the cytoplasmic membrane of photosynthetic bacteria. X-ray diffraction data collected over the resolution range 30.0–2.1 Å show that binding of the lipid to the protein involves a combination of ionic interactions between the protein and the lipid headgroup and van der Waals interactions between the lipid tails and the electroneutral intramembrane surface of the protein. In the headgroup region, ionic interactions involve polar groups of a number of residues, the protein backbone, and bound water molecules. The lipid tails sit along largely hydrophobic grooves in the irregular surface of the protein. In addition to providing new information on the immediate lipid environment of a key integral membrane protein, this study provides the first, to our knowledge, high-resolution x-ray crystal structure for cardiolipin. The possible significance of this interaction between an integral membrane protein and cardiolipin is considered.

Interaction between RGS7 and polycystin

Regulators of G protein signaling (RGS) proteins accelerate the intrinsic GTPase activity of certain Gα subunits and thereby modulate a number of G protein-dependent signaling cascades. Currently, little is known about the regulation of RGS proteins themselves. We identified a short-lived RGS protein, RGS7, that is rapidly degraded through the proteasome pathway. The degradation of RGS7 is inhibited by interaction with a C-terminal domain of polycystin, the protein encoded by PKD1, a gene involved in autosomal-dominant polycystic kidney disease. Furthermore, membranous expression of C-terminal polycystin relocalized RGS7. Our results indicate that rapid degradation and interaction with integral membrane proteins are potential means of regulating RGS proteins.

Energetics of the interaction between water and the helical peptide group and its role in determining helix propensities

The alanine helix provides a model system for studying the energetics of interaction between water and the helical peptide group, a possible major factor in the energetics of protein folding. Helix formation is enthalpy-driven (−1.0 kcal/mol per residue). Experimental transfer data (vapor phase to aqueous) for amides give the enthalpy of interaction with water of the amide group as ≈−11.5 kcal/mol. The enthalpy of the helical peptide hydrogen bond, computed for the gas phase by quantum mechanics, is −4.9 kcal/mol. These numbers give an enthalpy deficit for helix formation of −7.6 kcal/mol. To study this problem, we calculate the electrostatic solvation free energy (ESF) of the peptide groups in the helical and β-strand conformations, by using the delphi program and parse parameter set. Experimental data show that the ESF values of amides are almost entirely enthalpic. Two key results are: in the β-strand conformation, the ESF value of an interior alanine peptide group is −7.9 kcal/mol, substantially less than that of N-methylacetamide (−12.2 kcal/mol), and the helical peptide group is solvated with an ESF of −2.5 kcal/mol. These results reduce the enthalpy deficit to −1.5 kcal/mol, and desolvation of peptide groups through partial burial in the random coil may account for the remainder. Mutant peptides in the helical conformation show ESF differences among nonpolar amino acids that are comparable to observed helix propensity differences, but the ESF differences in the random coil conformation still must be subtracted.

Assembly of the neutrophil respiratory burst oxidase: A direct interaction between p67PHOX and cytochrome b558

Activation of the phagocyte NADPH oxidase complex requires the assembly of the cytosolic factors p47PHOX, p67PHOX, p40PHOX, and Rac1 or Rac2, with the membrane-bound cytochrome b558. Whereas the interaction of p47PHOX with cytochrome b558 is well established, an interaction between p67PHOX and cytochrome b558 has never been investigated. We report here a direct interaction between p67PHOX and cytochrome b558. First, labeled p67PHOX recognizes a 91-kDa band in specific granules from a normal patient but not from a cytochrome b558-deficient patient. Second, p67PHOX binds to cytochrome b558 that has been bound to nitrocellulose. Third, GTP-p67PHOX bound to glutathione agarose is able to pull down cytochrome b558. Rac1-GTP or Rac1-GDP increased the binding of p67PHOX to cytochrome b558, suggesting that at least one of the oxidase-related functions of Rac1 is to promote the interaction between p67PHOX and cytochrome b558.

The Interaction between Cold and Light Controls the Expression of the Cold-Regulated Barley Gene cor14b and the Accumulation of the Corresponding Protein1

We report the expression of the barley (Hordeum vulgare L.) COR (cold-regulated) gene cor14b (formerly pt59) and the accumulation of its chloroplast-localized protein product. A polyclonal antibody raised against the cor14b-encoded protein detected two chloroplast COR proteins: COR14a and COR14b. N-terminal sequencing of COR14a and expression of cor14b in Arabidopsis plants showed that COR14a is not encoded by the cor14b sequence, but it shared homology with the wheat (Triticum aestivum L.) WCS19 COR protein. The expression of cor14b was strongly impaired in the barley albino mutant an, suggesting the involvement of a plastidial factor in the control of gene expression. Low-level accumulation of COR14b was induced by cold treatment in etiolated plants, although cor14b expression and protein accumulation were enhanced after a short light pulse. Light quality was a determining factor in regulating gene expression: red or blue but not far-red or green light pulses were able to promote COR14b accumulation in etiolated plants, suggesting that phytochrome and blue light photoreceptors may be involved in the control of cor14b gene expression. Maximum accumulation of COR14b was reached only when plants were grown and/or hardened under the standard photoperiod. The effect of light on the COR14b stability was demonstrated by using transgenic Arabidopsis. These plants constitutively expressed cor14b mRNAs regardless of temperature and light conditions; nevertheless, green plants accumulated about twice as much COR14b protein as etiolated plants.

An isoleucine/leucine residue in the carboxyltransferase domain of acetyl-CoA carboxylase is critical for interaction with aryloxyphenoxypropionate and cyclohexanedione inhibitors

cDNA fragments encoding the carboxyltransferase domain of the multidomain plastid acetyl-CoA carboxylase (ACCase) from herbicide-resistant maize and from herbicide-sensitive and herbicide-resistant Lolium rigidum were cloned and sequenced. A Leu residue was found in ACCases from herbicide-resistant plants at a position occupied by Ile in all ACCases from sensitive grasses studied so far. Leu is present at the equivalent position in herbicide-resistant ACCases from other eukaryotes. Chimeric ACCases containing a 1000-aa fragment of two ACCase isozymes found in a herbicide-resistant maize were expressed in a yeast ACC1 null mutant to test herbicide sensitivity of the enzyme in vivo and in vitro. One of the enzymes was resistant/tolerant, and one was sensitive to haloxyfop and sethoxydim, rendering the gene-replacement yeast strains resistant and sensitive to these compounds, respectively. The sensitive enzyme has an Ile residue, and the resistant one has a Leu residue at the putative herbicide-binding site. Additionally, a single Ile to Leu replacement at an equivalent position changes the wheat plastid ACCase from sensitive to resistant. The effect of the opposite substitution, Leu to Ile, makes Toxoplasma gondii apicoplast ACCase resistant to haloxyfop and clodinafop. In this case, inhibition of the carboxyltransferase activity of ACCase (second half-reaction) of a large fragment of the Toxoplasma enzyme expressed in Escherichia coli was tested. The critical amino acid residue is located close to a highly conserved motif of the carboxyltransferase domain, which is probably a part of the enzyme active site, providing the basis for the activity of fop and dim herbicides.