Quantifying Reproducibility in Computational Biology: The Case of the Tuberculosis Drugome

How easy is it to reproduce the results found in a typical computational biology paper? Either through experience or intuition the reader will already know that the answer is with difficulty or not at all. In this paper we attempt to quantify this difficulty by reproducing a previously published paper for different classes of users (ranging from users with little expertise to domain experts) and suggest ways in which the situation might be improved. Quantification is achieved by estimating the time required to reproduce each of the steps in the method described in the original paper and make them part of an explicit workflow that reproduces the original results. Reproducing the method took several months of effort, and required using new versions and new software that posed challenges to reconstructing and validating the results. The quantification leads to “reproducibility maps” that reveal that novice researchers would only be able to reproduce a few of the steps in the method, and that only expert researchers with advance knowledge of the domain would be able to reproduce the method in its entirety. The workflow itself is published as an online resource together with supporting software and data. The paper concludes with a brief discussion of the complexities of requiring reproducibility in terms of cost versus benefit, and a desiderata with our observations and guidelines for improving reproducibility. This has implications not only in reproducing the work of others from published papers, but reproducing work from one’s own laboratory.

Integration of quantum key distribution in metropolitan area networks

The deployment of Quantum Key Distribution forces the development of QKD-links to be operated in current and next-generation photonic metro-access networks. These highly heterogeneous architectures determine the conditions QKD-links need to be optimized for.

Selected topics in the Hydrodynamics of floating offshore wind turbines: heave damping and second order surge motions

La Energía eléctrica producida mediante tecnología eólica flotante es uno de los recursos más prometedores para reducir la dependencia de energía proveniente de combustibles fósiles. Esta tecnología es de especial interés en países como España, donde la plataforma continental es estrecha y existen pocas áreas para el desarrollo de estructuras fijas. Entre los diferentes conceptos flotantes, esta tesis se ha ocupado de la tipología semisumergible. Estas plataformas pueden experimentar movimientos resonantes en largada y arfada. En largada, dado que el periodo de resonancia es largo estos puede ser inducidos por efectos de segundo orden de deriva lenta que pueden tener una influencia muy significativa en las cargas en los fondeos. En arfada las fuerzas de primer orden pueden inducir grandes movimientos y por tanto la correcta determinación del amortiguamiento es esencial para la analizar la operatividad de la plataforma. Esta tesis ha investigado estos dos efectos, para ello se ha usado como caso base el diseño de una plataforma desarrollada en el proyecto Europeo Hiprwind. La plataforma se compone de 3 columnas cilíndricas unidas mediante montantes estructurales horizontales y diagonales, Los cilindros proporcionan flotabilidad y momentos adrizante. A la base de cada columna se le ha añadido un gran “Heave Plate” o placa de cierre. El diseño es similar a otros diseños previos (Windfloat). Se ha fabricado un modelo a escala de una de las columnas para el estudio detallado del amortiguamiento mediante oscilaciones forzadas. Las dimensiones del modelo (1m diámetro en la placa de cierre) lo hacen, de los conocidos por el candidato, el mayor para el que se han publicado datos. El diseño del cilindro se ha realizado de tal manera que permite la fijación de placas de cierre planas o con refuerzo, ambos modelos se han fabricado y analizado. El modelo con refuerzos es una reproducción exacta del diseño a escala real incluyendo detalles distintivos del mismo, siendo el más importante la placa vertical perimetral. Los ensayos de oscilaciones forzadas se han realizado para un rango de frecuencias, tanto para el disco plano como el reforzado. Se han medido las fuerzas durante los ensayos y se han calculado los coeficientes de amortiguamiento y de masa añadida. Estos coeficientes son necesarios para el cálculo del fondeo mediante simulaciones en el dominio del tiempo. Los coeficientes calculados se han comparado con la literatura existente, con cálculos potenciales y por ultimo con cálculos CFD. Para disponer de información relevante para el diseño estructural de la plataforma se han medido y analizado experimentalmente las presiones en la parte superior e inferior de cada placa de cierre. Para la correcta estimación numérica de las fuerzas de deriva lenta en la plataforma se ha realizado una campaña experimental que incluye ensayos con modelo cautivo de la plataforma completa en olas bicromaticas. Pese a que estos experimentos no reproducen un escenario de oleaje realista, los mismos permiten una verificación del modelo numérico mediante la comparación de fuerzas medidas en el modelo físico y el numérico. Como resultados de esta tesis podemos enumerar las siguientes conclusiones. 1. El amortiguamiento y la masa añadida muestran una pequeña dependencia con la frecuencia pero una gran dependencia con la amplitud del movimiento. siendo coherente con investigaciones existentes. 2. Las medidas con la placa de cierre reforzada con cierre vertical en el borde, muestra un amortiguamiento significativamente menor comparada con la placa plana. Esto implica que para ensayos de canal es necesario incluir estos detalles en el modelo. 3. La masa añadida no muestra grandes variaciones comparando placa plana y placa con refuerzos. 4. Un coeficiente de amortiguamiento del 6% del crítico se puede considerar conservador para el cálculo en el dominio de la frecuencia. Este amortiguamiento es equivalente a un coeficiente de “drag” de 4 en elementos de Morison cuadráticos en las placas de cierre usadas en simulaciones en el dominio del tiempo. 5. Se han encontrado discrepancias en algunos valores de masa añadida y amortiguamiento de la placa plana al comparar con datos publicados. Se han propuesto algunas explicaciones basadas en las diferencias en la relación de espesores, en la distancia a la superficie libre y también relacionadas con efectos de escala. 6. La presión en la placa con refuerzos son similares a las de la placa plana, excepto en la zona del borde donde la placa con refuerzo vertical induce una gran diferencias de presiones entre la cara superior e inferior. 7. La máxima diferencia de presión escala coherentemente con la fuerza equivalente a la aceleración de la masa añadida distribuida sobre la placa. 8. Las masas añadidas calculadas con el código potencial (WADAM) no son suficientemente precisas, Este software no contempla el modelado de placas de pequeño espesor con dipolos, la poca precisión de los resultados aumenta la importancia de este tipo de elementos al realizar simulaciones con códigos potenciales para este tipo de plataformas que incluyen elementos de poco espesor. 9. Respecto al código CFD (Ansys CFX) la precisión de los cálculos es razonable para la placa plana, esta precisión disminuye para la placa con refuerzo vertical en el borde, como era de esperar dado la mayor complejidad del flujo. 10. Respecto al segundo orden, los resultados, en general, muestran que, aunque la tendencia en las fuerzas de segundo orden se captura bien con los códigos numéricos, se observan algunas reducciones en comparación con los datos experimentales. Las diferencias entre simulaciones y datos experimentales son mayores al usar la aproximación de Newman, que usa únicamente resultados de primer orden para el cálculo de las fuerzas de deriva media. 11. Es importante remarcar que las tendencias observadas en los resultados con modelo fijo cambiarn cuando el modelo este libre, el impacto que los errores en las estimaciones de fuerzas segundo orden tienen en el sistema de fondeo dependen de las condiciones ambientales que imponen las cargas ultimas en dichas líneas. En cualquier caso los resultados que se han obtenido en esta investigación confirman que es necesaria y deseable una detallada investigación de los métodos usados en la estimación de las fuerzas no lineales en las turbinas flotantes para que pueda servir de guía en futuros diseños de estos sistemas. Finalmente, el candidato espera que esta investigación pueda beneficiar a la industria eólica offshore en mejorar el diseño hidrodinámico del concepto semisumergible. ABSTRACT Electrical power obtained from floating offshore wind turbines is one of the promising resources which can reduce the fossil fuel energy consumption and cover worldwide energy demands. The concept is the most competitive in countries, such as Spain, where the continental shelf is narrow and does not provide space for fixed structures. Among the different floating structures concepts, this thesis has dealt with the semisubmersible one. Platforms of this kind may experience resonant motions both in surge and heave directions. In surge, since the platform natural period is long, such resonance can be excited with second order slow drift forces and may have substantial influence on mooring loads. In heave, first order forces can induce significant motion, whose damping is a crucial factor for the platform downtime. These two topics have been investigated in this thesis. To this aim, a design developed during HiPRWind EU project, has been selected as reference case study. The platform is composed of three cylindrical legs, linked together by a set of structural braces. The cylinders provide buoyancy and restoring forces and moments. Large circular heave plates have been attached to their bases. The design is similar to other documented in literature (e.g. Windfloat), which implies outcomes could have a general value. A large scale model of one of the legs has been built in order to study heave damping through forced oscillations. The final dimensions of the specimen (one meter diameter discs) make it, to the candidate’s knowledge, the largest for which data has been published. The model design allows for the fitting of either a plain solid heave plate or a flapped reinforced one; both have been built. The latter is a model scale reproduction of the prototype heave plate and includes some distinctive features, the most important being the inclusion of a vertical flap on its perimeter. The forced oscillation tests have been conducted for a range of frequencies and amplitudes, with both the solid plain model and the vertical flap one. Forces have been measured, from which added mass and damping coefficients have been obtained. These are necessary to accurately compute time-domain simulations of mooring design. The coefficients have been compared with literature, and potential flow and CFD predictions. In order to provide information for the structural design of the platform, pressure measurements on the top and bottom side of the heave discs have been recorded and pressure differences analyzed. In addition, in order to conduct a detailed investigation on the numerical estimations of the slow-drift forces of the HiPRWind platform, an experimental campaign involving captive (fixed) model tests of a model of the whole platform in bichromatic waves has been carried out. Although not reproducing the more realistic scenario, these tests allowed a preliminary verification of the numerical model based directly on the forces measured on the structure. The following outcomes can be enumerated: 1. Damping and added mass coefficients show, on one hand, a small dependence with frequency and, on the other hand, a large dependence with the motion amplitude, which is coherent with previously published research. 2. Measurements with the prototype plate, equipped with the vertical flap, show that damping drops significantly when comparing this to the plain one. This implies that, for tank tests of the whole floater and turbine, the prototype plate, equipped with the flap, should be incorporated to the model. 3. Added mass values do not suffer large alterations when comparing the plain plate and the one equipped with a vertical flap. 4. A conservative damping coefficient equal to 6% of the critical damping can be considered adequate for the prototype heave plate for frequency domain analysis. A corresponding drag coefficient equal to 4.0 can be used in time domain simulations to define Morison elements. 5. When comparing to published data, some discrepancies in added mass and damping coefficients for the solid plain plate have been found. Explanations have been suggested, focusing mainly on differences in thickness ratio and distance to the free surface, and eventual scale effects. 6. Pressures on the plate equipped with the vertical flap are similar in magnitude to those of the plain plate, even though substantial differences are present close to the edge, where the flap induces a larger pressure difference in the reinforced case. 7. The maximum pressure difference scales coherently with the force equivalent to the acceleration of the added mass, distributed over the disc surface. 8. Added mass coefficient values predicted with the potential solver (WADAM) are not accurate enough. The used solver does not contemplate modeling thin plates with doublets. The relatively low accuracy of the results highlights the importance of these elements when performing potential flow simulations of offshore platforms which include thin plates. 9. For the full CFD solver (Ansys CFX), the accuracy of the computations is found reasonable for the plain plate. Such accuracy diminishes for the disc equipped with a vertical flap, an expected result considering the greater complexity of the flow. 10. In regards to second order effects, in general, the results showed that, although the main trend in the behavior of the second-order forces is well captured by the numerical predictions, some under prediction of the experimental values is visible. The gap between experimental and numerical results is more pronounced when Newman’s approximation is considered, making use exclusively of the mean drift forces calculated in the first-order solution. 11. It should be observed that the trends observed in the fixed model test may change when the body is free to float, and the impact that eventual errors in the estimation of the second-order forces may have on the mooring system depends on the characteristics of the sea conditions that will ultimately impose the maximum loads on the mooring lines. Nevertheless, the preliminary results obtained in this research do confirm that a more detailed investigation of the methods adopted for the estimation of the nonlinear wave forces on the FOWT would be welcome and may provide some further guidance for the design of such systems. As a final remark, the candidate hopes this research can benefit the offshore wind industry in improving the hydrodynamic design of the semi-submersible concept.

Quantitative modeling of stochastic systems in molecular biology by using stochastic Petri nets

An integrated understanding of molecular and developmental biology must consider the large number of molecular species involved and the low concentrations of many species in vivo. Quantitative stochastic models of molecular interaction networks can be expressed as stochastic Petri nets (SPNs), a mathematical formalism developed in computer science. Existing software can be used to define molecular interaction networks as SPNs and solve such models for the probability distributions of molecular species. This approach allows biologists to focus on the content of models and their interpretation, rather than their implementation. The standardized format of SPNs also facilitates the replication, extension, and transfer of models between researchers. A simple chemical system is presented to demonstrate the link between stochastic models of molecular interactions and SPNs. The approach is illustrated with examples of models of genetic and biochemical phenomena where the UltraSAN package is used to present results from numerical analysis and the outcome of simulations.

Adjustment of conformational flexibility is a key event in the thermal adaptation of proteins

3-Isopropylmalate dehydrogenase (IPMDH, E.C. 1.1.1.85) from the thermophilic bacterium Thermus thermophilus HB8 is homologous to IPMDH from the mesophilic Escherichia coli, but has an approximately 17°C higher melting temperature. Its temperature optimum is 22–25°C higher than that of the E. coli enzyme; however, it is hardly active at room temperature. The increased conformational rigidity required to stabilize the thermophilic enzyme against heat denaturation might explain its different temperature-activity profile. Hydrogen/deuterium exchange studies were performed on this thermophilic-mesophilic enzyme pair to compare their conformational flexibilities. It was found that Th. thermophilus IPMDH is significantly more rigid at room temperature than E. coli IPMDH, whereas the enzymes have nearly identical flexibilities under their respective optimal working conditions, suggesting that evolutionary adaptation tends to maintain a “corresponding state” regarding conformational flexibility. These observations confirm that conformational fluctuations necessary for catalytic function are restricted at room temperature in the thermophilic enzyme, suggesting a close relationship between conformational flexibility and enzyme function.

Purification of ribonucleotide reductase subunits Y1, Y2, Y3, and Y4 from yeast: Y4 plays a key role in diiron cluster assembly

Ribonucleotide reductases (RNRs) catalyze the conversion of nucleotides to deoxynucleotides. Class I RNRs are composed of two types of subunits: RNR1 contains the active site for reduction and the binding sites for the nucleotide allosteric effectors. RNR2 contains the diiron-tyrosyl radical (Y⋅) cofactor essential for the reduction process. Studies in yeast have recently identified four RNR subunits: Y1 and Y3, Y2 and Y4. These proteins have been expressed in Saccharomyces cerevisiae and in Escherichia coli and purified to ≈90% homogeneity. The specific activity of Y1 isolated from yeast and E. coli is 0.03 μmol⋅min−1⋅mg−1 and of (His)6-Y2 [(His)6-Y2-K387N] from yeast is 0.037 μmol⋅min−1⋅mg−1 (0.125 μmol⋅min−1⋅mg−1). Y2, Y3, and Y4 isolated from E. coli have no measurable activity. Efforts to generate Y⋅ in Y2 or Y4 using Fe2+, O2, and reductant have been unsuccessful. However, preliminary studies show that incubation of Y4 and Fe2+ with inactive E. coli Y2 followed by addition of O2 generates Y2 with a specific activity of 0.069 μmol⋅min−1⋅mg−1 and a Y⋅. A similar experiment with (His)6-Y2-K387N, Y4, O2, and Fe2+ results in an increase in its specific activity to 0.30 μmol⋅min−1⋅mg−1. Studies with antibodies to Y4 and Y2 reveal that they can form a complex in vivo. Y4 appears to play an important role in diiron-Y⋅ assembly of Y2.

Identification of a gene encoding an acyl CoA:diacylglycerol acyltransferase, a key enzyme in triacylglycerol synthesis

Triacylglycerols are quantitatively the most important storage form of energy for eukaryotic cells. Acyl CoA:diacylglycerol acyltransferase (DGAT, EC 2.3.1.20) catalyzes the terminal and only committed step in triacylglycerol synthesis, by using diacylglycerol and fatty acyl CoA as substrates. DGAT plays a fundamental role in the metabolism of cellular diacylglycerol and is important in higher eukaryotes for physiologic processes involving triacylglycerol metabolism such as intestinal fat absorption, lipoprotein assembly, adipose tissue formation, and lactation. DGAT is an integral membrane protein that has never been purified to homogeneity, nor has its gene been cloned. We identified an expressed sequence tag clone that shared regions of similarity with acyl CoA:cholesterol acyltransferase, an enzyme that also uses fatty acyl CoA as a substrate. Expression of a mouse cDNA for this expressed sequence tag in insect cells resulted in high levels of DGAT activity in cell membranes. No other acyltransferase activity was detected when a variety of substrates, including cholesterol, were used as acyl acceptors. The gene was expressed in all tissues examined; during differentiation of NIH 3T3-L1 cells into adipocytes, its expression increased markedly in parallel with increases in DGAT activity. The identification of this cDNA encoding a DGAT will greatly facilitate studies of cellular glycerolipid metabolism and its regulation.

Nitric oxide plays a key role in adaptive control of locomotion in cat

Diverse roles in cellular functions have been ascribed to nitric oxide (NO), and its involvement in induction of long-term depression in cerebellar Purkinje cells has been demonstrated. Manipulations of NO concentration or its synthesis in cerebellar tissues therefore provide a means for investigating roles of NO in cerebellar functions at both cellular and behavioral levels. We tested adaptive control of locomotion to perturbation in cats, and found that this form of motor learning was abolished by application of either an inhibitor of NO synthase or a scavenger of NO to the cerebellar cortical locomotion area. This finding supports the view that NO in the cerebellum plays a key role in motor learning.

Reciprocal interactions between β1-integrin and epidermal growth factor receptor in three-dimensional basement membrane breast cultures: A different perspective in epithelial biology

Anchorage and growth factor independence are cardinal features of the transformed phenotype. Although it is logical that the two pathways must be coregulated in normal tissues to maintain homeostasis, this has not been demonstrated directly. We showed previously that down-modulation of β1-integrin signaling reverted the malignant behavior of a human breast tumor cell line (T4–2) derived from phenotypically normal cells (HMT-3522) and led to growth arrest in a three-dimensional (3D) basement membrane assay in which the cells formed tissue-like acini (14). Here, we show that there is a bidirectional cross-modulation of β1-integrin and epidermal growth factor receptor (EGFR) signaling via the mitogen-activated protein kinase (MAPK) pathway. The reciprocal modulation does not occur in monolayer (2D) cultures. Antibody-mediated inhibition of either of these receptors in the tumor cells, or inhibition of MAPK kinase, induced a concomitant down-regulation of both receptors, followed by growth-arrest and restoration of normal breast tissue morphogenesis. Cross-modulation and tissue morphogenesis were associated with attenuation of EGF-induced transient MAPK activation. To specifically test EGFR and β1-integrin interdependency, EGFR was overexpressed in nonmalignant cells, leading to disruption of morphogenesis and a compensatory up-regulation of β1-integrin expression, again only in 3D. Our results indicate that when breast cells are spatially organized as a result of contact with basement membrane, the signaling pathways become coupled and bidirectional. They further explain why breast cells fail to differentiate in monolayer cultures in which these events are mostly uncoupled. Moreover, in a subset of tumor cells in which these pathways are misregulated but functional, the cells could be “normalized” by manipulating either pathway.

Cloning of chlorophyllase, the key enzyme in chlorophyll degradation: Finding of a lipase motif and the induction by methyl jasmonate

Chlorophyllase (Chlase) is the first enzyme involved in chlorophyll (Chl) degradation and catalyzes the hydrolysis of ester bond to yield chlorophyllide and phytol. In the present study, we isolated the Chlase cDNA. We synthesized degenerate oligo DNA probes based on the internal amino acid sequences of purified Chlase from Chenopodium album, screened the C. album cDNA library, and cloned a cDNA (CaCLH, C. album chlorophyll-chlorophyllido hydrolase). The deduced amino acid sequence (347 aa residues) had a lipase motif overlapping with an ATP/GTP-binding motif (P-loop). CaCLH possibly was localized in the extraplastidic part of the cell, because a putative signal sequence for endoplasmic reticulum is at the N terminus. The amino acid sequence shared 37% identity with a function-unknown gene whose mRNA is inducible by coronatine and methyl jasmonate (MeJA) in Arabidopsis thaliana (AtCLH1). We expressed the gene products of AtCLH1 and of CaCLH in Escherichia coli, and they similarly exhibited Chlase activity. Moreover, we isolated another full-length cDNA based on an Arabidopsis genomic fragment and expressed it in E. coli, demonstrating the presence of the second Arabidopsis CLH gene (AtCLH2). No typical feature of signal sequence was identified in AtCLH1, whereas AtCLH2 had a typical signal sequence for chloroplast. AtCLH1 mRNA was induced rapidly by a treatment of MeJA, which is known to promote senescence and Chl degradation in plants, and a high mRNA level was maintained up to 9 h. AtCLH2, however, did not respond to MeJA.

A targeted mutation in the murine gene encoding the high density lipoprotein (HDL) receptor scavenger receptor class B type I reveals its key role in HDL metabolism

Plasma high density lipoprotein (HDL), which protects against atherosclerosis, is thought to remove cholesterol from peripheral tissues and to deliver cholesteryl esters via a selective uptake pathway to the liver (reverse cholesterol transport) and steroidogenic tissues (e.g., adrenal gland for storage and hormone synthesis). Despite its physiologic and pathophysiologic importance, the cellular metabolism of HDL has not been well defined. The class B, type I scavenger receptor (SR-BI) has been proposed to play an important role in HDL metabolism because (i) it is a cell surface HDL receptor which mediates selective cholesterol uptake in cultured cells, (ii) its physiologically regulated expression is most abundant in the liver and steroidogenic tissues, and (iii) hepatic overexpression dramatically lowers plasma HDL. To test directly the normal role of SR-BI in HDL metabolism, we generated mice with a targeted null mutation in the SR-BI gene. In heterozygous and homozygous mutants relative to wild-type controls, plasma cholesterol concentrations were increased by ≈31% and 125%, respectively, because of the formation of large, apolipoprotein A-I (apoA-I)-containing particles, and adrenal gland cholesterol content decreased by 42% and 72%, respectively. The plasma concentration of apoA-I, the major protein in HDL, was unchanged in the mutants. This, in conjunction with the increased lipoprotein size, suggests that the increased plasma cholesterol in the mutants was due to decreased selective cholesterol uptake. These results provide strong support for the proposal that in mice the gene encoding SR-BI plays a key role in determining the levels of plasma lipoprotein cholesterol (primarily HDL) and the accumulation of cholesterol stores in the adrenal gland. If it has a similar role in controlling plasma HDL in humans, SR-BI may influence the development and progression of atherosclerosis and may be an attractive candidate for therapeutic intervention in this disease.

Signaling through fibroblast growth factor receptor 2b plays a key role in the development of the exocrine pancreas

The development of the pancreas depends on epithelial-mesenchymal interactions. Fibroblast growth factors (FGFs) and their receptors (FGFRs 1–4) have been identified as mediators of epithelial-mesenchymal interactions in different organs. We show here that FGFR-2 IIIb and its ligands FGF-1, FGF-7, and FGF-10 are expressed throughout pancreatic development. We also show that in mesenchyme-free cultures of embryonic pancreatic epithelium FGF-1, FGF-7, and FGF-10 stimulate the growth, morphogenesis, and cytodifferentiation of the exocrine cells of the pancreas. The role of FGFs signaling through FGFR-2 IIIb was further investigated by inhibiting FGFR-2 IIIb signaling in organocultures of pancreatic explants (epithelium + mesenchyme) by using either antisense FGFR-2 IIIb oligonucleotides or a soluble recombinant FGFR-2 IIIb protein. Abrogation of FGFR-2 IIIb signaling resulted in a considerable reduction in the size of the explants and in a 2-fold reduction of the development of the exocrine cells. These results demonstrate that FGFs signaling through FGFR-2 IIIb play an important role in the development of the exocrine pancreas.