Automated search of natively folded protein fragments for high-throughput structure determination in structural genomics

Structural genomic projects envision almost routine protein structure determinations, which are currently imaginable only for small proteins with molecular weights below 25,000 Da. For larger proteins, structural insight can be obtained by breaking them into small segments of amino acid sequences that can fold into native structures, even when isolated from the rest of the protein. Such segments are autonomously folding units (AFU) and have sizes suitable for fast structural analyses. Here, we propose to expand an intuitive procedure often employed for identifying biologically important domains to an automatic method for detecting putative folded protein fragments. The procedure is based on the recognition that large proteins can be regarded as a combination of independent domains conserved among diverse organisms. We thus have developed a program that reorganizes the output of BLAST searches and detects regions with a large number of similar sequences. To automate the detection process, it is reduced to a simple geometrical problem of recognizing rectangular shaped elevations in a graph that plots the number of similar sequences at each residue of a query sequence. We used our program to quantitatively corroborate the premise that segments with conserved sequences correspond to domains that fold into native structures. We applied our program to a test data set composed of 99 amino acid sequences containing 150 segments with structures listed in the Protein Data Bank, and thus known to fold into native structures. Overall, the fragments identified by our program have an almost 50% probability of forming a native structure, and comparable results are observed with sequences containing domain linkers classified in SCOP. Furthermore, we verified that our program identifies AFU in libraries from various organisms, and we found a significant number of AFU candidates for structural analysis, covering an estimated 5 to 20% of the genomic databases. Altogether, these results argue that methods based on sequence similarity can be useful for dissecting large proteins into small autonomously folding domains, and such methods may provide an efficient support to structural genomics projects.     

A large deletion causes apparent homozygosity for the D1152H mutation in the cystic fibrosis transmembrane regulator (CFTR) gene

We report the case of a patient with an apparent homozygosity for the D1152H mutation located in exon 18 of the cystic fibrosis transmembrane conductance regulator (CFTR) gene. The parents had no personal history of cystic fibrosis (CF) and referred to our laboratory after the diagnosis of fetal bowel hyperechogenicity. The proband presented with meconium ileus and normal sweat chloride test. Sequencing of the CFTR exon 18 together with quantitative genomic assays, such as real-time PCR and the multiplex ligation probe amplification (MLPA) techniques, were performed and revealed that the father was heterozygous for the D1152H mutation and the mother carried a large deletion of the CFTR gene encompassing the genomic sequence including the same mutation. The child inherited D1152H from his father and the large deletion of the CFTR gene from his mother. We suggest that D1152H likely acts as a mild mutation with a dominant effect on the severe deletion of exon 18, considering that after 3 years of clinical examinations the child shows no classical signs and symptoms of CF. Not testing for large deletions in subjects with apparent homozygosity for a mutated CFTR allele could lead to the misidentification of CFTR mutation carrier status.     

Intra-specific hybridization: generator of genetic diversification and heterosis in Andrographis paniculata Nees. A bridge from extinction to survival

Andrographis paniculata (AP) has been stated as a low-diverse, endangered and red-listed plant species. Self-pollinated mating system, being an introduced species and experiencing a bottleneck as well as over exploitation cause such a consequence. Inter and intra-specific hybridizations have been suggested as essential techniques for generating genetic diversity. To test the effect of intra-specific hybridization on diversification and heterosis of AP, seven accessions were outcrossed manually in all 21 possible combinations. Three types of markers including morphological, phytochemical and RAPD markers were employed to evaluate the mentioned hypothesis. The results revealed that hybridization acted as a powerful engine for diversification of AP as it caused heterotic expression of the studied traits, simultaneously. Initially, it seems that additive and non-additive gene effects both can be considered as the genetic basis of heterosis in AP for the investigated traits. Agronomic and morphological traits were differentiated from each other, while positive heterosis was recorded mainly for agronomic traits but not for the morphological traits. Intra-specific hybridization increased the genetic diversity in AP population. Nevertheless, a part of this variation could also be attributed to the negative heterosis. The current exploration demonstrated the first ever conducted manual intra-specific hybridization among AP accessions in a mass scale. However, the 17 RAPD primers produced a monomorph pattern, but perhaps increasing the number of markers can feature a new genetic profile in this plant.     

Association between Glu298Asp/677C-T single nucleotide polymorphism in the eNOS/MTHRF gene and blood stasis syndrome of ischemic stroke

Blood stasis syndrome of ischemic stroke (BSS-IS) is a common clinical phenotype that may be affected by certain mutagenic environmental factors or chemotherapeutic drugs; however, the role of susceptibility genes remains unclear. Previous studies have shown that ischemic stroke (IS) was closely associated with the Glu298Asp polymorphism in the eNOS gene and the 677C-T (Ala→Val) polymorphism in methylenetetrahydrofolate reductase (MTHRF) gene. Therefore, these two single nucleotide polymorphisms (SNPs) were selected to detect their associations with BSS-IS in this study. A SNP chip was employed to screen the SNP variation between both groups, and the results were verified using denaturing high-performance liquid chromatography (DHPLC) and restriction fragment length polymorphism (RFLP). The results confirmed that the TT genotype of Glu298Asp in the eNOS gene may be one of the risk factors associated with BSS-IS, while the genotype of 677C-T (Ala→Val) in the MTHRF gene may not be relevant to BSS-IS.     

Differential dependence on DNA ligase of type II restriction enzymes: a practical way toward ligase-free DNA automaton

DNA computing study is a new paradigm in computer science and biological computing fields. As one of DNA computing approaches, DNA automaton is composed of the hardware, input DNA molecule and state transition molecules. By now restriction enzymes are key hardware for DNA computing automaton. It has been found that DNA computing efficiency may be independent on DNA ligases when type IIS restriction enzymes like FokI are used as hardware. In this study, we compared FokI with four other distinct enzymes HgaI, BsmFI, BbsI, and BseMII, and found their differential independence on T4 DNA ligase when performing automaton reactions. Since DNA automaton is a potential powerful tool to tackle gene relationship in genomic network scale, the feasible ligase-free DNA automaton may set an initial base to develop functional DNA automata for various DNA technology development and implications in genetics study in the near future.     

Characterization of the first angiotensin-converting like enzyme in bacteria: Ancestor ACE is already active

Angiotensin-converting enzyme (ACE) is a metallopeptidase that converts angiotensin I into angiotensin II. ACE is crucial in the control of cardiovascular and renal homeostasis and fertility in mammals. In vertebrates, both transmembrane and soluble ACE, containing one or two active sites, have been characterized. So far, only soluble, single domain ACEs from invertebrates have been cloned, and these have been implicated in reproduction in insects. Furthermore, an ACE-related carboxypeptidase was recently characterized in Leishmania, a unicellular eukaryote, suggesting the existence of ACE in more distant organisms. Interestingly, in silico databank analysis revealed that bacterial DNA sequences could encode putative ACE-like proteins, strikingly similar to vertebrates' enzymes. To gain more insight into the bacterial enzymes, we cloned the putative ACE from the phytopathogenic bacterium, Xanthomonas axonopodis pv. citri, named XcACE. The 2 kb open reading frame encodes a 672-amino-acid soluble protein containing a single active site. In vitro expression and biochemical characterization revealed that XcACE is a functional 72 kDa dipeptidyl-carboxypeptidase. As in mammals, this metalloprotease hydrolyses angiotensin I into angiotensin II. XcACE is sensitive to ACE inhibitors and chloride ions concentration. Variations in the active site residues, highlighted by structural modelling, can account for the different substrate selectivity and inhibition profile compared to human ACE. XcACE characterization demonstrates that ACE is an ancestral enzyme, provoking questions about its appearance and structure/activity specialisation during the course of evolution.     

Multiple evolutionary mechanisms drive papillomavirus diversification

The circular, double-stranded 8-kb DNA genome of papillomaviruses (PVes) consists mainly of 4 large genes, E1, E2, L2, and L1. Approximately 150 papillomavirus genomes have been sequenced to date. We analyzed a representative sample of 53 PVes genomes using maximum likelihood, Bayesian inference, maximum parsimony, and distance-based methods both on nucleotide (nt) and on amino acid (aa) alignments. When the 4 genes were analyzed separately, aa-inferred phylogenies contradicted each other less than nt-inferred trees (judged by partition homogeneity tests). In particular, gene combinations including the L2 gene generated significant incongruence (P < 0.001). Combined analyses of the remaining genes E1-E2-L1 produced a well-supported phylogeny including supertaxon beta + gamma + pi + xi-PVes (infecting Artiodactyla, Carnivora, Primates, and Rodentia) and supertaxon kappa + lambda + mu + nu + sigma-PVes (infecting Carnivora, Lagomorpha, Primates, and Rodentia). Based on the tree topology, host-linked evolution appears plausible at shallow, rather than deeper, taxonomic levels. Diversification within PVes may also involve adaptive radiation establishing different niches (within a single-host species) and recombination events (within single-host cells). Heterogeneous groups of closely related PVes infecting, for example, humans and domestic animals such as hamster, dog, and cattle suggest multiple infections across species borders. Additional evolutionary phenomena such as strong codon usage preferences, and computational biases including reconstruction artifacts and insufficient taxon sampling, may contribute to the incomplete resolution of deep phylogenetic nodes. The molecular data globally supports a complex evolutionary scenario for PVes, which is driven by multiple mechanisms but not exclusively by coevolution with corresponding hosts.     

MOSBIE: a tool for comparison and analysis of rule-based biochemical models

Background

Mechanistic models that describe the dynamical behaviors of biochemical systems are common in computational systems biology, especially in the realm of cellular signaling. The development of families of such models, either by a single research group or by different groups working within the same area, presents significant challenges that range from identifying structural similarities and differences between models to understanding how these differences affect system dynamics.

Results

We present the development and features of an interactive model exploration system, MOSBIE, which provides utilities for identifying similarities and differences between models within a family. Models are clustered using a custom similarity metric, and a visual interface is provided that allows a researcher to interactively compare the structures of pairs of models as well as view simulation results.

Conclusions

We illustrate the usefulness of MOSBIE via two case studies in the cell signaling domain. We also present feedback provided by domain experts and discuss the benefits, as well as the limitations, of the approach.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-316) contains supplementary material, which is available to authorized users.     

MSDA,a proteomics software suite for in‐depth Mass Spectrometry Data Analysis using grid computing

One of the major bottlenecks in the proteomics field today resides in the computational interpretation of the massive data generated by the latest generation of high‐throughput MS instruments. MS/MS datasets are constantly increasing in size and complexity and it becomes challenging to comprehensively process such huge datasets and afterwards deduce most relevant biological information. The Mass Spectrometry Data Analysis (MSDA, https://msda.unistra.fr ) online software suite provides a series of modules for in‐depth MS/MS data analysis. It includes a custom databases generation toolbox, modules for filtering and extracting high‐quality spectra, for running high‐performance database and de novo searches, and for extracting modified peptides spectra and functional annotations. Additionally, MSDA enables running the most computationally intensive steps, namely database and de novo searches, on a computer grid thus providing a net time gain of up to 99% for data processing.