A knowledge-driven protocol for prediction of proteins of interest with an emphasis on biosynthetic pathways.

This protocol describes a stepwise process to identify proteins of interest from a query proteome derived from NGS data. We implemented this protocol on Moringa oleifera transcriptome to identify proteins involved in secondary metabolite and vitamin biosynthesis and ion transport. This knowledge-driven protocol identifies proteins using an integrated approach involving sensitive sequence search and evolutionary relationships. We make use of functionally important residues (FIR) specific for the query protein family identified through its homologous sequences and literature. We screen protein hits based on the clustering with true homologues through phylogenetic tree reconstruction complemented with the FIR mapping. The protocol was validated for the protein hits through qRT-PCR and transcriptome quantification. Our protocol demonstrated a higher specificity as compared to other methods, particularly in distinguishing cross-family hits. This protocol was effective in transcriptome data analysis of M. oleifera as described in Pasha et al.•Knowledge-driven protocol to identify secondary metabolite synthesizing protein in a highly specific manner.•Use of functionally important residues for screening of true hits.•Beneficial for metabolite pathway reconstruction in any (species, metagenomics) NGS data.

Dataset for the combined transcriptome assembly of M. oleifera and functional annotation.

In this paper, we present the data acquired during transcriptome analysis of the plant Moringa oleifera [1] from five different tissues (root, stem, leaf, flower and seed) by RNA sequencing. A total of 271 million reads were assembled with an N50 of 2094 bp. The combined transcriptome was assessed for transcript abundance across five tissues. The protein coding genes identified from the transcripts were annotated and used for orthology analysis. Further, enzymes involved in the biosynthesis of select medicinally important secondary metabolites, vitamins and ion transporters were identified and their expression levels across tissues were examined. The data generated by RNA sequencing has been deposited to NCBI public repository under the accession number PRJNA394193 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA394193).

The transcriptome enables the identification of candidate genes behind medicinal value of Drumstick tree (Moringa oleifera).

Moringa oleifera is a plant well-known for its nutrition value, drought resistance and medicinal properties. cDNA libraries from five different tissues (leaf, root, stem, seed and flower) of M. oleifera cultivar Bhagya were generated and sequenced. We developed a bioinformatics pipeline to assemble transcriptome, along with the previously published M. oleifera genome, to predict 17,148 gene models. Few candidate genes related to biosynthesis of secondary metabolites, vitamins and ion transporters were identified. Expressions were further confirmed by real-time quantitative PCR experiments for few promising leads. Quantitative estimation of metabolites, as well as elemental analysis, was also carried out to support our observations. Enzymes in the biosynthesis of vitamins and metabolites like quercetin and kaempferol are highly expressed in leaves, flowers and seeds. The expression of iron transporters and calcium storage proteins were observed in root and leaves. In general, leaves retain the highest amount of small molecules of interest.

Molecular mechanisms and structural features of cardiomyopathy-causing troponin T mutants in the tropomyosin overlap region.

Point mutations in genes encoding sarcomeric proteins are the leading cause of inherited primary cardiomyopathies. Among them are mutations in the TNNT2 gene that encodes cardiac troponin T (TnT). These mutations are clustered in the tropomyosin (Tm) binding region of TnT, TNT1 (residues 80-180). To understand the mechanistic changes caused by pathogenic mutations in the TNT1 region, six hypertrophic cardiomyopathy (HCM) and two dilated cardiomyopathy (DCM) mutants were studied by biochemical approaches. Binding assays in the absence and presence of actin revealed changes in the affinity of some, but not all, TnT mutants for Tm relative to WT TnT. HCM mutants were hypersensitive and DCM mutants were hyposensitive to Ca2+ in regulated actomyosin ATPase activities. To gain better insight into the disease mechanism, we modeled the structure of TNT1 and its interactions with Tm. The stability predictions made by the model correlated well with the affinity changes observed in vitro of TnT mutants for Tm. The changes in Ca2+ sensitivity showed a strong correlation with the changes in binding affinity. We suggest the primary reason by which these TNNT2 mutations between residues 92 and 144 cause cardiomyopathy is by changing the affinity of TnT for Tm within the TNT1 region.

The myosin mesa and the basis of hypercontractility caused by hypertrophic cardiomyopathy mutations.

Hypertrophic cardiomyopathy (HCM) is primarily caused by mutations in ß-cardiac myosin and myosin-binding protein-C (MyBP-C). Changes in the contractile parameters of myosin measured so far do not explain the clinical hypercontractility caused by such mutations. We propose that hypercontractility is due to an increase in the number of myosin heads (S1) that are accessible for force production. In support of this hypothesis, we demonstrate myosin tail (S2)-dependent functional regulation of actin-activated human ß-cardiac myosin ATPase. In addition, we show that both S2 and MyBP-C bind to S1 and that phosphorylation of either S1 or MyBP-C weakens these interactions. Importantly, the S1-S2 interaction is also weakened by four myosin HCM-causing mutations but not by two other mutations. To explain these experimental results, we propose a working structural model involving multiple interactions, including those with myosin's own S2 and MyBP-C, that hold myosin in a sequestered state.

Multidimensional structure-function relationships in human ß-cardiac myosin from population-scale genetic variation.

Myosin motors are the fundamental force-generating elements of muscle contraction. Variation in the human ß-cardiac myosin heavy chain gene (MYH7) can lead to hypertrophic cardiomyopathy (HCM), a heritable disease characterized by cardiac hypertrophy, heart failure, and sudden cardiac death. How specific myosin variants alter motor function or clinical expression of disease remains incompletely understood. Here, we combine structural models of myosin from multiple stages of its chemomechanical cycle, exome sequencing data from two population cohorts of 60,706 and 42,930 individuals, and genetic and phenotypic data from 2,913 patients with HCM to identify regions of disease enrichment within ß-cardiac myosin. We first developed computational models of the human ß-cardiac myosin protein before and after the myosin power stroke. Then, using a spatial scan statistic modified to analyze genetic variation in protein 3D space, we found significant enrichment of disease-associated variants in the converter, a kinetic domain that transduces force from the catalytic domain to the lever arm to accomplish the power stroke. Focusing our analysis on surface-exposed residues, we identified a larger region significantly enriched for disease-associated variants that contains both the converter domain and residues on a single flat surface on the myosin head described as the myosin mesa. Notably, patients with HCM with variants in the enriched regions have earlier disease onset than patients who have HCM with variants elsewhere. Our study provides a model for integrating protein structure, large-scale genetic sequencing, and detailed phenotypic data to reveal insight into time-shifted protein structures and genetic disease.

Genome sequencing of herb Tulsi (Ocimum tenuiflorum) unravels key genes behind its strong medicinal properties.

BACKGROUND: Krishna Tulsi, a member of Lamiaceae family, is a herb well known for its spiritual, religious and medicinal importance in India. The common name of this plant is 'Tulsi' (or 'Tulasi' or 'Thulasi') and is considered sacred by Hindus. We present the draft genome of Ocimum tenuiflurum L (subtype Krishna Tulsi) in this report. The paired-end and mate-pair sequence libraries were generated for the whole genome sequenced with the Illumina Hiseq 1000, resulting in an assembled genome of 374 Mb, with a genome coverage of 61 % (612 Mb estimated genome size). We have also studied transcriptomes (RNA-Seq) of two subtypes of O. tenuiflorum, Krishna and Rama Tulsi and report the relative expression of genes in both the varieties. RESULTS: The pathways leading to the production of medicinally-important specialized metabolites have been studied in detail, in relation to similar pathways in Arabidopsis thaliana and other plants. Expression levels of anthocyanin biosynthesis-related genes in leaf samples of Krishna Tulsi were observed to be relatively high, explaining the purple colouration of Krishna Tulsi leaves. The expression of six important genes identified from genome data were validated by performing q-RT-PCR in different tissues of five different species, which shows the high extent of urosolic acid-producing genes in young leaves of the Rama subtype. In addition, the presence of eugenol and ursolic acid, implied as potential drugs in the cure of many diseases including cancer was confirmed using mass spectrometry. CONCLUSIONS: The availability of the whole genome of O.tenuiflorum and our sequence analysis suggests that small amino acid changes at the functional sites of genes involved in metabolite synthesis pathways confer special medicinal properties to this herb.

Mechanistic heterogeneity in contractile properties of α-tropomyosin (TPM1) mutants associated with inherited cardiomyopathies.

The most frequent known causes of primary cardiomyopathies are mutations in the genes encoding sarcomeric proteins. Among those are 30 single-residue mutations in TPM1, the gene encoding α-tropomyosin. We examined seven mutant tropomyosins, E62Q, D84N, I172T, L185R, S215L, D230N, and M281T, that were chosen based on their clinical severity and locations along the molecule. The goal of our study was to determine how the biochemical characteristics of each of these mutant proteins are altered, which in turn could provide a structural rationale for treatment of the cardiomyopathies they produce. Measurements of Ca(2+) sensitivity of human ß-cardiac myosin ATPase activity are consistent with the hypothesis that hypertrophic cardiomyopathies are hypersensitive to Ca(2+) activation, and dilated cardiomyopathies are hyposensitive. We also report correlations between ATPase activity at maximum Ca(2+) concentrations and conformational changes in TnC measured using a fluorescent probe, which provide evidence that different substitutions perturb the structure of the regulatory complex in different ways. Moreover, we observed changes in protein stability and protein-protein interactions in these mutants. Our results suggest multiple mechanistic pathways to hypertrophic and dilated cardiomyopathies. Finally, we examined a computationally designed mutant, E181K, that is hypersensitive, confirming predictions derived from in silico structural analysis.

Myosinome: a database of myosins from select eukaryotic genomes to facilitate analysis of sequence-structure-function relationships.

Myosins are one of the largest protein superfamilies with 24 classes. They have conserved structural features and catalytic domains yet show huge variation at different domains resulting in a variety of functions. Myosins are molecules driving various kinds of cellular processes and motility until the level of organisms. These are ATPases that utilize the chemical energy released by ATP hydrolysis to bring about conformational changes leading to a motor function. Myosins are important as they are involved in almost all cellular activities ranging from cell division to transcriptional regulation. They are crucial due to their involvement in many congenital diseases symptomatized by muscular malfunctions, cardiac diseases, deafness, neural and immunological dysfunction, and so on, many of which lead to death at an early age. We present Myosinome, a database of selected myosin classes (myosin II, V, and VI) from five model organisms. This knowledge base provides the sequences, phylogenetic clustering, domain architectures of myosins and molecular models, structural analyses, and relevant literature of their coiled-coil domains. In the current version of Myosinome, information about 71 myosin sequences belonging to three myosin classes (myosin II, V, and VI) in five model organisms (Homo Sapiens, Mus musculus, D. melanogaster, C. elegans and S. cereviseae) identified using bioinformatics surveys are presented, and several of them are yet to be functionally characterized. As these proteins are involved in congenital diseases, such a database would be useful in short-listing candidates for gene therapy and drug development. The database can be accessed from http://caps.ncbs.res.in/myosinome.

Structural attributes for the recognition of weak and anomalous regions in coiled-coils of myosins and other motor proteins.

BACKGROUND: Coiled-coils are found in different proteins like transcription factors, myosin tail domain, tropomyosin, leucine zippers and kinesins. Analysis of various structures containing coiled-coils has revealed the importance of electrostatic and hydrophobic interactions. In such domains, regions of different strength of interactions need to be identified since they could be biologically relevant. FINDINGS: We have updated our coiled-coil validation webserver, now called COILCHECK+, where new features were added to efficiently identify the strength of interaction at the interface region and measure the density of charged residues and hydrophobic residues. We have examined charged residues and hydrophobic ladders, using a new algorithm called CHAHO, which is incorporated within COILCHECK + server. CHAHO permits the identification of spatial charged residue patches and the continuity of hydrophobic ladder which stabilizes and destabilizes the coiled-coil structure. CONCLUSIONS: The availability of such computational tools should be useful to understand the importance of spatial clustering of charged residues and the continuity of hydrophobic residues at the interface region of coiled-coil dimers. COILCHECK + is a structure based tool to validate coiled-coil stability; it can be accessed at http://caps.ncbs.res.in/coilcheckplus.