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.

Comparative proteome analysis of psychrophilic versus mesophilic bacterial species: Insights into the molecular basis of cold adaptation of proteins.

BACKGROUND: Cold adapted or psychrophilic organisms grow at low temperatures, where most of other organisms cannot grow. This adaptation requires a vast array of sequence, structural and physiological adjustments. To understand the molecular basis of cold adaptation of proteins, we analyzed proteomes of psychrophilic and mesophilic bacterial species and compared the differences in amino acid composition and substitution patterns to investigate their likely association with growth temperatures. RESULTS: In psychrophilic bacteria, serine, aspartic acid, threonine and alanine are overrepresented in the coil regions of secondary structures, whilst glutamic acid and leucine are underrepresented in the helical regions. Compared to mesophiles, psychrophiles comprise a significantly higher proportion of amino acids that contribute to higher protein flexibility in the coil regions of proteins, such as those with tiny/small or neutral side chains. Amino acids with aliphatic, basic, aromatic and hydrophilic side chains are underrepresented in the helical regions of proteins of psychrophiles. The patterns of amino acid substitutions between the orthologous proteins of psychrophiles versus mesophiles are significantly different for several amino acids when compared to their substitutions in orthologous proteins of within the mesophiles or psychrophiles. CONCLUSION: Current results provide quantitative substitution preferences (or avoidance) of amino acids that lead to the adaptation of proteins to cold temperatures. These finding would help future efforts in selecting mutations for rational design of proteins with enhanced psychrophilic properties.

Modeling and analysis of molecularinteraction between Smurf1-WW2 domain and various isoforms of LIM mineralization protein.

LIM Mineralization Protein-1 (LMP-1) has been cloned and shown to be osteoinductive. Our efforts to understand the mode of action of LMP-1 led to the determination that LMP-1 interacts with Smad Ubiquitin Regulatory Factor-1 (Smurf1). Smurf1 targets osteogenic Smads, Smad1/5, for ubiquitin-mediated proteasomal degradation. Smurf1 interaction with LMP-1 or Smads is based on the presence of unique WW-domain interacting motif in these target molecules. By performing site-directed mutagenesis and binding studies in vitro on purified recombinant proteins, we identified a specific motif within the osteogenic region of several LMP isoforms that is necessary for Smurf1 interaction. Similarly, we have identified that the WW2 domain of Smurf1 is necessary for target protein interaction. Here, we present a homology-based modeling of the Smurf1 WW2 domain and its interacting motif of LMP-1. We performed computational docking of the interacting domains in Smurf1 and LMPs to identify the key amino acid residues involved in their binding regions. In support of the computational predictions, we also present biochemical evidence supporting the hypothesis that the physical interaction of Smurf1 and osteoinductive forms of LMP may prevent Smurf1 from targeting osteogenic Smads by ubiquitin-mediated proteasomal degradation.

Molecular interaction between Smurf1 WW2 domain and PPXY motifs of Smad1, Smad5, and Smad6--modeling and analysis.

The ubiquitin-proteasome proteolytic pathway is essential for various important biological processes including cell cycle progression, gene transcription, and signal transduction. One of the important regulatory mechanisms by which the bone-inducing activity of the bone morphogenetic protein (BMP) signaling is modulated involves ubiquitin-mediated proteasomal degradation. The BMP induced receptor signal is transmitted intracellularly by phosphorylation of Smad proteins by the activated receptor I. The phosphorylated Smads 1, 5, and 8 (R-Smads) oligomerize with the co-Smad (Smad4). The complex, thus, formed translocates to the nucleus and interacts with other cofactors to regulate the expression of downstream target genes. R-Smads contain PPXY motif in the linker region that interacts with Smad ubiquitin regulatory factor 1 (Smurf1), an E3 ubiquitin ligase that catalyzes ubiquitination of target proteins for proteasomal degradation. Smurf1 contains a HECT domain, a C2 domain, and 2 WW domains (WW1, WW2). The PPXY motif in target proteins and its interaction with Smurf1 may form the basis for regulation of steady-state levels of Smads in controlling BMP-responsiveness of cells. Here, we present a homology-based model of the Smurf1 WW2 domain and the target octa-peptides containing PPXY motif of Smurf1-interacting Smads. We carried out docking of Smurf1 WW2 domain with the PPXY motifs of Smad1, Smad5, and Smad6 and identified the key amino acid residues involved in interaction. Furthermore, we present experimental evidence that WW2 domain of Smurf1 does indeed interact with the Smad proteins and that the deletion of WW2 domain of Smurf1 results in loss of its binding to Smads using the purified recombinant proteins. Finally, we also present data confirming that the deletion of WW2 domain in Smurf1 abolishes its ubiquitination activity on Smad1 in an in vitro ubiquitination assay. It shows that the interaction between the WW domain and Smad PPXY motif is a key step in Smurf1-mediated ubiquitination of its natural targets such as Smad1, Smad5, and Smad6. This work facilitates further strategies to unravel the biological function of such interactions and help in designing effective mimetic compounds that either mimic or disrupt the specific interaction.

Genome inventory and analysis of nuclear hormone receptors in Tetraodon nigroviridis.

Nuclear hormone receptors (NRs) form a large superfamily of ligand-activated transcription factors, which regulate genes underlying a wide range of (patho) physiological phenomena. Availability of the full genome sequence of Tetraodon nigroviridis facilitated a genome wide analysis of the NRs in fish genome. Seventy one NRs were found in Tetraodon and were compared with mammalian and fish NR family members. In general, there is a higher representation of NRs in fish genomes compared to mammalian ones. They showed high diversity across classes as observed by phylogenetic analysis. Nucleotide substitution rates show strong negative selection among fish NRs except for pregnane x receptor (PxR), estrogen receptor (ER) and liver x receptor (LxR). This may be attributed to crucial role played by them in metabolism and detoxification of xenobiotic and endobiotic compounds and might have resulted in slight positive selection. Chromosomal mapping and pairwise comparisons of NR distribution in Tetraodon and humans led to the identification of nine syntenic NR regions, of which three are common among fully sequenced vertebrate genomes. Gene structure analysis shows strong conservation of exon structures among orthologoues. Whereas paralogous members show different splicing patterns with intron gain or loss and addition or substitution of exons played a major role in evolution of NR superfamily.

Genome wide survey of G protein-coupled receptors in Tetraodon nigroviridis.

BACKGROUND: The G-protein-coupled receptors (GPCRs) constitute one of the largest and most ancient superfamilies of membrane proteins. They play a central role in physiological processes affecting almost all aspects of the life cycle of an organism. Availability of the complete sets of putative members of a family from diverse species provides the basis for cross genome comparative studies. RESULTS: We have defined the repertoire of GPCR superfamily of Tetraodon complement with the availability of complete sequence of the freshwater puffer fish Tetraodon nigroviridis. Almost all 466 Tetraodon GPCRs (Tnig-GPCRs) identified had a clear human homologue. 189 putative human and Tetraodon GPCR orthologous pairs could be identified. Tetraodon GPCRs are classified into five GRAFS families, by phylogenetic analysis, concurrent with human GPCR classification. CONCLUSION: Direct comparison of GPCRs in Tetraodon and human genomes displays a high level of orthology and supports large-scale gene duplications in Tetraodon. Examples of lineage specific gene expansions were also observed in opsin and odorant receptors. The human and Tetraodon GPCR sequences are analogous in terms of GPCR subfamilies but display disproportionate numbers of receptors at the subfamily level. The teleost genome with its expanded set of GPCRs provides additional and interesting comparators to study both evolution and function of these receptors.

Cross genome phylogenetic analysis of human and Drosophila G protein-coupled receptors: application to functional annotation of orphan receptors.

BACKGROUND: The cell-membrane G-protein coupled receptors (GPCRs) are one of the largest known superfamilies and are the main focus of intense pharmaceutical research due to their key role in cell physiology and disease. A large number of putative GPCRs are 'orphans' with no identified natural ligands. The first step in understanding the function of orphan GPCRs is to identify their ligands. Phylogenetic clustering methods were used to elucidate the chemical nature of receptor ligands, which led to the identification of natural ligands for many orphan receptors. We have clustered human and Drosophila receptors with known ligands and orphans through cross genome phylogenetic analysis and hypothesized higher relationship of co-clustered members that would ease ligand identification, as related receptors share ligands with similar structure or class. RESULTS: Cross-genome phylogenetic analyses were performed to identify eight major groups of GPCRs dividing them into 32 clusters of 371 human and 113 Drosophila proteins (excluding olfactory, taste and gustatory receptors) and reveal unexpected levels of evolutionary conservation across human and Drosophila GPCRs. We also observe that members of human chemokine receptors, involved in immune response, and most of nucleotide-lipid receptors (except opsins) do not have counterparts in Drosophila. Similarly, a group of Drosophila GPCRs (methuselah receptors), associated in aging, is not present in humans. CONCLUSION: Our analysis suggests ligand class association to 52 unknown Drosophila receptors and 95 unknown human GPCRs. A higher level of phylogenetic organization was revealed in which clusters with common domain architecture or cellular localization or ligand structure or chemistry or a shared function are evident across human and Drosophila genomes. Such analyses will prove valuable for identifying the natural ligands of Drosophila and human orphan receptors that can lead to a better understanding of physiological and pathological roles of these receptors.

Comparison of Analysis Tools for miRNA High Throughput Sequencing Using Nerve Crush as a Model.

Recent advances in sample preparation and analysis for next generation sequencing have made it possible to profile and discover new miRNAs in a high throughput manner. In the case of neurological disease and injury, these types of experiments have been more limited. Possibly because tissues such as the brain and spinal cord are inaccessible for direct sampling in living patients, and indirect sampling of blood and cerebrospinal fluid are affected by low amounts of RNA. We used a mouse model to examine changes in miRNA expression in response to acute nerve crush. We assayed miRNA from both muscle tissue and blood plasma. We examined how the depth of coverage (the number of mapped reads) changed the number of detectable miRNAs in each sample type. We also found that samples with very low starting amounts of RNA (mouse plasma) made high depth of mature miRNA coverage more difficult to obtain. Each tissue must be assessed independently for the depth of coverage required to adequately power detection of differential expression, weighed against the cost of sequencing that sample to the adequate depth. We explored the changes in total mapped reads and differential expression results generated by three different software packages: miRDeep2, miRNAKey, and miRExpress and two different analysis packages, DESeq and EdgeR. We also examine the accuracy of using miRDeep2 to predict novel miRNAs and subsequently detect them in the samples using qRT-PCR.

Virtual screening and selection of drug-like compounds to block noggin interaction with bone morphogenetic proteins.

Noggin is a major natural extracellular antagonist to bone morphogenetic proteins (BMPs) which binds to BMPs and blocks binding of them to BMP-specific receptors and thus negatively regulates BMP-induced osteoblastic differentiation. Bone morphogenetic proteins (BMPs) signal through heteromeric protein complexes composed of type I and type II serine/threonine kinase receptors. Preventing the BMP-2/noggin interaction will preserve free BMP-2 and enhance the efficacy of BMP-2 to induce bone formation. This work is an attempt to use the current understanding of BMP-2, and its interaction with its receptors and antagonist to design an inhibitor of BMP-2/noggin interaction with the goal of lowering the dose of BMP-2 required in clinical applications. The crystal structure of the BMP-7/noggin complex, the BMP-2/BMP receptor IA ectodomain complex and the extracellular domain of BMP receptor II monomer are known. We modeled the BMP-2 based on the structure of its homologue BMP-7 and its binding complex with noggin. We also modeled a complex of BMP-2/BMPRIA/BMPRII by modeling BMPRII and replacing ActRIIB in the BMP-2/BMPRIA/ActRIIB complex. We then identified the binding region of noggin with BMP-2 and the receptors with BMP-2. From the analysis of structures of these complexes and modeling we identified the key amino acids present in the entire interacting surfaces among these proteins that play important physiological role in the regulation of cell differentiation and bone metabolism. By in silico screening we selected and ranked several compounds that have high theoretical scores to bind to noggin to block BMP-noggin interaction.

Genome inventory and analysis of nuclear hormone receptors in Tetraodon nigroviridis.

Nuclear hormone receptors (NRs) form a large superfamily of ligand-activated transcription factors, which regulate genes underlying a wide range of (patho) physiological phenomena. Availability of the full genome sequence of Tetraodon nigroviridis facilitated a genome wide analysis of the NRs in fish genome. Seventy one NRs were found in Tetraodon and were compared with mammalian and fish NR family members. In general, there is a higher representation of NRs in fish genomes compared to mammalian ones. They showed high diversity across classes as observed by phylogenetic analysis. Nucleotide substitution rates show strong negative selection among fish NRs except for pregnane x receptor (PxR), estrogen receptor (ER) and liver x receptor (LxR). This may be attributed to crucial role played by them in metabolism and detoxification of xenobiotic and endobiotic compounds and might have resulted in slight positive selection. Chromosomal mapping and pairwise comparisons of NR distribution in Tetraodon and humans led to the identification of nine syntenic NR regions, of which three are common among fully sequenced vertebrate genomes. Gene structure analysis shows strong conservation of exon structures among orthologoues. Whereas paralogous members show different splicing patterns with intron gain or loss and addition or substitution of exons played a major role in evolution of NR superfamily.