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.

Energy consumption and carbon dioxide emissions in rail and road freight transport in Spain: a case study of car carriers and bulk petrochemicals

This article provides a new methodology for estimating fuel consumption and emissions by enabling a correct comparison between freight transportation modes. The approach is developed and integrated as a part of an intelligent transportation system dealing with goods movement. A key issue is related to energy consumption ratios and consequent CO2 emissions. Energy consumption ratios are often used based on transport demand. However, including other ratios based on transport supply can be useful. Furthermore, it is important to indicate which factors are associated with variations in energy consumption and emissions; especially of interest are parameters that have a higher incidence and order of magnitude, in order to fairly compare and understand the difference between transport modes and sub-modes. The study finds that the use of an energy consumption equation can improve the quality of the estimates. The study proposes that coefficients that define the energy consumption equation should be tested to determine market niches and sources of improvement in energy consumption according to the category of vehicles, fuel types used, and classes of products transported.

Lacunarity of the Spatial Distributions of Soil Types in Europe.

Lacunarity as a means of quantifying textural properties of spatial distributions suggests a classification into three main classes of the most abundant soils that cover 92% of Europe. Soils with a well-defined self-similar structure of the linear class are related to widespread spatial patterns that are nondominant but ubiquitous at continental scale. Fractal techniques have been increasingly and successfully applied to identify and describe spatial patterns in natural sciences. However, objects with the same fractal dimension can show very different optical properties because of their spatial arrangement. This work focuses primary attention on the geometrical structure of the geographical patterns of soils in Europe. We made use of the European Soil Database to estimate lacunarity indexes of the most abundant soils that cover 92% of the surface of Europe and investigated textural properties of their spatial distribution. We observed three main classes corresponding to three different patterns that displayed the graphs of lacunarity functions, that is, linear, convex, and mixed. They correspond respectively to homogeneous or self-similar, heterogeneous or clustered and those in which behavior can change at different ranges of scales. Finally, we discuss the pedological implications of that classification.

YUCCA8 and YUCCA9 overexpression reveals a link between auxin signaling and lignification through the induction of ethylene biosynthesis.

Auxin is associated with the regulation of virtually every aspect of plant growth and development. Many previous genetic and biochemical studies revealed that, among the proposed routes for the production of auxin, the so-called indole-3-pyruvic acid (IPA) pathway is the main source for indole-3-acetic acid (IAA) in plants. The IPA pathway involves the action of 2 classes of enzymes, tryptophan-pyruvate aminotransferases (TRYPTOPHAN AMINOTRANSFERASE OF ARABIDOPSIS 1(TAA1)/TRYPTOPHAN AMINOTRANSFERASE RELATED (TAR)) and flavin monooxygenases (YUCCA). Both enzyme classes appear to be encoded by small gene families in Arabidopsis consisting of 5 and 11 members, respectively. We recently showed that it is possible to induce transcript accumulation of 2 YUCCA genes, YUC8 and YUC9, by methyl jasmonate treatment. Both gene products were demonstrated to contribute to auxin biosynthesis in planta.1 Here we report that the overexpression of YUC8 as well as YUC9 led to strong lignification of plant aerial tissues. Furthermore, new evidence indicates that this abnormally strong secondary growth is linked to increased levels of ethylene production.

Sight distance studies on roads: influence of digital elevation models and roadside elements

Sight distance plays an important role in road traffic safety. Two types of Digital Elevation Models (DEMs) are utilized for the estimation of available sight distance in roads: Digital Terrain Models (DTMs) and Digital Surface Models (DSMs). DTMs, which represent the bare ground surface, are commonly used to determine available sight distance at the design stage. Additionally, the use of DSMs provides further information about elements by the roadsides such as trees, buildings, walls or even traffic signals which may reduce available sight distance. This document analyses the influence of three classes of DEMs in available sight distance estimation. For this purpose, diverse roads within the Region of Madrid (Spain) have been studied using software based on geographic information systems. The study evidences the influence of using each DEM in the outcome as well as the pros and cons of using each model.

Is the polarity of content producers strongly influenced by the results of the event?

This paper presents an approach to compare two types of data, subjective data (Polarity of Pan American Games 2011 event by country) and objective data (the number of medals won by each participating country), based on the Pearson corre- lation. When dealing with events described by people, knowledge acquisition is difficult because their structure is heterogeneous and subjective. A first step towards knowing the polarity of the information provided by people consists in automatically classifying the posts into clusters according to their polarity. The authors carried out a set of experiments using a corpus that consists of 5600 posts extracted from 168 Internet resources related to a specific event: the 2011 Pan American games. The approach is based on four components: a crawler, a filter, a synthesizer and a polarity analyzer. The PanAmerican approach automatically classifies the polarity of the event into clusters with the following results: 588 positive, 336 neutral, and 76 negative. Our work found out that the polarity of the content produced was strongly influenced by the results of the event with a correlation of .74. Thus, it is possible to conclude that the polarity of content is strongly affected by the results of the event. Finally, the accuracy of the PanAmerican approach is: .87, .90, and .80 according to the precision of the three classes of polarity evaluated.

Factor-specific modulation of CREB-binding protein acetyltransferase activity

CREB-binding proteins (CBP) and p300 are essential transcriptional coactivators for a large number of regulated DNA-binding transcription factors, including CREB, nuclear receptors, and STATs. CBP and p300 function in part by mediating the assembly of multiprotein complexes that contain additional cofactors such as p300/CBP interacting protein (p/CIP), a member of the p160/SRC family of coactivators, and the p300/CBP associated factor p/CAF. In addition to serving as molecular scaffolds, CBP and p300 each possess intrinsic acetyltransferase activities that are required for their function as coactivators. Here we report that the adenovirus E1A protein inhibits the acetyltransferase activity of CBP on binding to the C/H3 domain, whereas binding of CREB, or a CREB/E1A fusion protein to the KIX domain, fails to inhibit CBP acetyltransferase activity. Surprisingly, p/CIP can either inhibit or stimulate CBP acetyltransferase activity depending on the specific substrate evaluated and the functional domains present in the p/CIP protein. While the CBP interaction domain of p/CIP inhibits acetylation of histones H3, H4, or high mobility group by CBP, it enhances acetylation of other substrates, such as Pit-1. These observations suggest that the acetyltransferase activities of CBP/p300 and p/CAF can be differentially modulated by factors binding to distinct regions of CBP/p300. Because these interactions are likely to result in differential effects on the coactivator functions of CBP/p300 for different classes of transcription factors, regulation of CBP/p300 acetyltransferase activity may represent a mechanism for integration of diverse signaling pathways.

In vivo function of mutated spliced leader RNAs in Caenorhabditis elegans

The role of spliced leader RNA (SL RNA) in trans-splicing in Caenorhabditis elegans has been studied through a combination of in vitro mutagenesis and in vivo complementation of rrs-1 mutant nematodes, which lack endogenous SL1 RNA. Three classes of mutant SL1 RNAs have been found—those that rescue the lethal phenotype at low concentration of transforming DNA, those that rescue at high but not low concentration, and those that do not rescue at all. These studies showed that some mutations in the otherwise highly conserved 22-nt spliced leader are tolerated for splicing and post-splicing events. A longer spliced leader also can be tolerated but only when present in high copy number. Changes in the first 16 nucleotides result in the appearance of no SL RNA, consistent with the in vitro studies by others showing that the SL1 RNA promoter partly resides within the spliced leader sequence.

Rapid evolution of sex-related genes in Chlamydomonas

Biological speciation ultimately results in prezygotic isolation—the inability of incipient species to mate with one another–but little is understood about the selection pressures and genetic changes that generate this outcome. The genus Chlamydomonas comprises numerous species of unicellular green algae, including numerous geographic isolates of the species C. reinhardtii. This diverse collection has allowed us to analyze the evolution of two sex-related genes: the mid gene of C. reinhardtii, which determines whether a gamete is mating-type plus or minus, and the fus1 gene, which dictates a cell surface glycoprotein utilized by C. reinhardtii plus gametes to recognize minus gametes. Low stringency Southern analyses failed to detect any fus1 homologs in other Chlamydomonas species and detected only one mid homolog, documenting that both genes have diverged extensively during the evolution of the lineage. The one mid homolog was found in C. incerta, the species in culture that is most closely related to C. reinhardtii. Its mid gene carries numerous nonsynonymous and synonymous codon changes compared with the C. reinhardtii mid gene. In contrast, very high sequence conservation of both the mid and fus1 sequences is found in natural isolates of C. reinhardtii, indicating that the genes are not free to drift within a species but do diverge dramatically between species. Striking divergence of sex determination and mate recognition genes also has been encountered in a number of other eukaryotic phyla, suggesting that unique, and as yet unidentified, selection pressures act on these classes of genes during the speciation process.

Molecular cloning and functional analysis of SUT-1, a sulfate transporter from human high endothelial venules

High endothelial venules (HEV) are specialized postcapillary venules found in lymphoid organs and chronically inflamed tissues that support high levels of lymphocyte extravasation from the blood. One of the major characteristics of HEV endothelial cells (HEVEC) is their capacity to incorporate large amounts of sulfate into sialomucin-type counter-receptors for the lymphocyte homing receptor L-selectin. Here, we show that HEVEC express two functional classes of sulfate transporters defined by their differential sensitivity to the anion-exchanger inhibitor 4,4′-diisothiocyanostilbene-2,2′-disulfonic acid (DIDS), and we report the molecular characterization of a DIDS-resistant sulfate transporter from human HEVEC, designated SUT-1. SUT-1 belongs to the family of Na+-coupled anion transporters and exhibits 40–50% amino acid identity with the rat renal Na+/sulfate cotransporter, NaSi-1, as well as with the human and rat Na+/dicarboxylate cotransporters, NaDC-1/SDCT1 and NaDC-3/SDCT2. Functional expression studies in cRNA-injected Xenopus laevis oocytes showed that SUT-1 mediates high levels of Na+-dependent sulfate transport, which is resistant to DIDS inhibition. The SUT-1 gene mapped to human chromosome 7q33. Northern blotting analysis revealed that SUT-1 exhibits a highly restricted tissue distribution, with abundant expression in placenta. Reverse transcription–PCR analysis indicated that SUT-1 and the diastrophic dysplasia sulfate transporter (DTD), one of the two known human DIDS-sensitive sulfate transporters, are coexpressed in HEVEC. SUT-1 and DTD could correspond, respectively, to the DIDS-resistant and DIDS-sensitive components of sulfate uptake in HEVEC. Together, these results demonstrate that SUT-1 is a distinct human Na+-coupled sulfate transporter, likely to play a major role in sulfate incorporation in HEV.

Direct interaction of Gβγ with a C-terminal Gβγ-binding domain of the Ca2+ channel α1 subunit is responsible for channel inhibition by G protein-coupled receptors

Several classes of voltage-gated Ca2+ channels (VGCCs) are inhibited by G proteins activated by receptors for neurotransmitters and neuromodulatory peptides. Evidence has accumulated to indicate that for non-L-type Ca2+ channels the executing arm of the activated G protein is its βγ dimer (Gβγ). We report below the existence of two Gβγ-binding sites on the A-, B-, and E-type α1 subunits that form non-L-type Ca2+ channels. One, reported previously, is in loop 1 connecting transmembrane domains I and II. The second is located approximately in the middle of the ca. 600-aa-long C-terminal tails. Both Gβγ-binding regions also bind the Ca2+ channel β subunit (CCβ), which, when overexpressed, interferes with inhibition by activated G proteins. Replacement in α1E of loop 1 with that of the G protein-insensitive and Gβγ-binding-negative loop 1 of α1C did not abolish inhibition by G proteins, but the exchange of the α1E C terminus with that of α1C did. This and properties of α1E C-terminal truncations indicated that the Gβγ-binding site mediating the inhibition of Ca2+ channel activity is the one in the C terminus. Binding of Gβγ to this site was inhibited by an α1-binding domain of CCβ, thus providing an explanation for the functional antagonism existing between CCβ and G protein inhibition. The data do not support proposals that Gβγ inhibits α1 function by interacting with the site located in the loop I–II linker. These results define the molecular mechanism by which presynaptic G protein-coupled receptors inhibit neurotransmission.

Identification of two distinct human SMC protein complexes involved in mitotic chromosome dynamics

The structural maintenance of chromosomes (SMC) family member proteins previously were shown to play a critical role in mitotic chromosome condensation and segregation in yeast and Xenopus. Other family members were demonstrated to be required for DNA repair in yeast and mammals. Although several different SMC proteins were identified in different organisms, little is known about the SMC proteins in humans. Here, we report the identification of four human SMC proteins that form two distinct heterodimeric complexes in the cell, the human chromosome-associated protein (hCAP)-C and hCAP-E protein complex (hCAP-C/hCAP-E), and the human SMC1 (hSMC1) and hSMC3 protein complex (hSMC1/hSMC3). The hCAP-C/hCAP-E complex is the human ortholog of the Xenopus chromosome-associated protein (XCAP)-C/XCAP-E complex required for mitotic chromosome condensation. We found that a second complex, hSMC1/hSMC3, is required for metaphase progression in mitotic cells. Punctate vs. diffuse distribution patterns of the hCAP-C/hCAP-E and hSMC1/hSMC3 complexes in the interphase nucleus indicate independent behaviors of the two complexes during the cell cycle. These results suggest that two distinct classes of SMC protein complexes are involved in different aspects of mitotic chromosome organization in human cells.

Design of compounds that increase the absorption of polar molecules

Hydrophilic drugs are often poorly absorbed when administered orally. There has been considerable interest in the possibility of using absorption enhancers to promote absorption of polar molecules across membrane surfaces. The bile acids are one of the most widely investigated classes of absorption enhancers, but there is disagreement about what features of bile acid enhancers are responsible for their efficacy. We have designed a class of glycosylated bile acid derivatives to evaluate how increasing the hydrophilicity of the steroid nucleus affects the ability to transport polar molecules across membranes. Some of the glycosylated molecules are significantly more effective than taurocholate in promoting the intestinal absorption of a range of drugs, showing that hydrophobicity is not a critical parameter in transport efficacy, as previously suggested. Furthermore, the most effective glycosylated compound is also far less damaging to membranes than the best bile acid absorption promoters, presumably because it is more hydrophilic. The results reported here show that it is possible to decouple absorption-promoting activity from membrane damage, a finding that should spark interest in the design of new compounds to facilitate the delivery of polar drugs.

Photorepair mutants of Arabidopsis

UV radiation induces two major DNA damage products, the cyclobutane pyrimidine dimer (CPD) and, at a lower frequency, the pyrimidine (6–4) pyrimidinone dimer (6–4 product). Although Escherichia coli and Saccharomyces cerevisiae produce a CPD-specific photolyase that eliminates only this class of dimer, Arabidopsis thaliana, Drosophila melanogaster, Crotalus atrox, and Xenopus laevis have recently been shown to photoreactivate both CPDs and 6–4 products. We describe the isolation and characterization of two new classes of mutants of Arabidopsis, termed uvr2 and uvr3, that are defective in the photoreactivation of CPDs and 6–4 products, respectively. We demonstrate that the CPD photolyase mutation is genetically linked to a DNA sequence encoding a type II (metazoan) CPD photolyase. In addition, we are able to generate plants in which only CPDs or 6–4 products are photoreactivated in the nuclear genome by exposing these mutants to UV light and then allowing them to repair one or the other class of dimers. This provides us with a unique opportunity to study the biological consequences of each of these two major UV-induced photoproducts in an intact living system.

A common mechanism for the biosynthesis of methoxy and cyclopropyl mycolic acids in Mycobacterium tuberculosis

Mycobacterium tuberculosis produces three classes of mycolic acids that differ primarily in the presence and nature of oxygen-containing substituents in the distal portion of the meromycolate branch. The methoxymycolate series has a methoxy group adjacent to a methyl branch, in addition to a cyclopropane in the proximal position. Using the gene for the enzyme that introduces the distal cyclopropane (cma1) as a probe, we have cloned and sequenced a cluster of genes coding for four highly homologous methyl transferases (mma1–4). When introduced into Mycobacterium smegmatis, this gene cluster conferred the ability to synthesize methoxymycolates. By determining the structure of the mycolic acids produced following expression of each of these genes individually and in combination, we have elucidated the biosynthetic steps responsible for the production of the major series of methoxymycolates. The mma4 gene product (MMAS-4) catalyzes an unusual S-adenosyl-l-methionine-dependent transformation of the distal cis-olefin into a secondary alcohol with an adjacent methyl branch. MMAS-3 O-methylates this secondary alcohol to form the corresponding methyl ether, and MMAS-2 introduces a cis-cyclopropane in the proximal position of the methoxy series. The similarity of these reactions and the enzymes that catalyze them suggests that some of the structural diversity of mycolic acids results from different chemical fates of a common cationic intermediate, which in turn results from methyl group addition to an olefinic mycolate precursor.