Genome bioinformatic analysis of nonsynonymous SNPs

Background: Genome-wide association studies of common diseases for common, low penetrance causal variants are underway. A proportion of these will alter protein sequences, the most common of which is the non-synonymous single nucleotide polymorphism (nsSNP). It would be an advantage if the functional effects of an nsSNP on protein structure and function could be predicted, both for the final identification process of a causal variant in a disease-associated chromosome region, and in further functional analyses of the nsSNP and its disease-associated protein. Results: In the present report we have compared and contrasted structure-and sequence-based methods of prediction to over 5500 genes carrying nearly 24,000 nsSNPs, by employing an automatic comparative modelling procedure to build models for the genes. The nsSNP information came from two sources, the OMIM database which are rare (minor allele frequency, MAF, 0.05, have no known link to a disease. For over 40% of the nsSNPs, structure-based methods predicted which of these sequence changes are likely to either disrupt the structure of the protein or interfere with the function or interactions of the protein. For the remaining 60%, we generated sequence-based predictions. Conclusion: We show that, in general, the prediction tools are able distinguish disease causing mutations from those mutations which are thought to have a neutral affect. We give examples of mutations in genes that are predicted to be deleterious and may have a role in disease. Contrary to previous reports, we also show that rare mutations are consistently predicted to be deleterious as often as commonly occurring nsSNPs.RIGHTS : This article is licensed under the BioMed Central licence at http://www.biomedcentral.com/about/license which is similar to the 'Creative Commons Attribution Licence'. In brief you may : copy, distribute, and display the work; make derivative works; or make commercial use of the work - under the following conditions: the original author must be given credit; for any reuse or distribution, it must be made clear to others what the license terms of this work are

Model of the complex of Parathyroid hormone-2 receptor and Tuberoinfundibular peptide of 39 residues

<p>Abstract</p> <p>Background</p> <p>We aim to propose interactions between the parathyroid hormone-2 receptor (PTH2R) and its ligand the tuberoinfundibular peptide of 39 residues (TIP39) by constructing a homology model of their complex. The two related peptides parathyroid hormone (PTH) and parathyroid hormone related protein (PTHrP) are compared with the complex to examine their interactions.</p> <p>Findings</p> <p>In the model, the hydrophobic N-terminus of TIP39 is buried in a hydrophobic part of the central cavity between helices 3 and 7. Comparison of the peptide sequences indicates that the main discriminator between the agonistic peptides TIP39 and PTH and the inactive PTHrP is a tryptophan-phenylalanine replacement. The model indicates that the smaller phenylalanine in PTHrP does not completely occupy the binding site of the larger tryptophan residue in the other peptides. As only TIP39 causes internalisation of the receptor and the primary difference being an aspartic acid in position 7 of TIP39 that interacts with histidine 396 in the receptor, versus isoleucine/histidine residues in the related hormones, this might be a trigger interaction for the events that cause internalisation.</p> <p>Conclusions</p> <p>A model is constructed for the complex and a trigger interaction for full agonistic activation between aspartic acid 7 of TIP39 and histidine 396 in the receptor is proposed.</p

Insights into the Binding of Phenyltiocarbamide (PTC) Agonist to Its Target Human TAS2R38 Bitter Receptor

Humans' bitter taste perception is mediated by the hTAS2R subfamily of the G protein-coupled membrane receptors (GPCRs). Structural information on these receptors is currently limited. Here we identify residues involved in the binding of phenylthiocarbamide (PTC) and in receptor activation in one of the most widely studied hTAS2Rs (hTAS2R38) by means of structural bioinformatics and molecular docking. The predictions are validated by site-directed mutagenesis experiments that involve specific residues located in the putative binding site and trans-membrane (TM) helices 6 and 7 putatively involved in receptor activation. Based on our measurements, we suggest that (i) residue N103 participates actively in PTC binding, in line with previous computational studies. (ii) W99, M100 and S259 contribute to define the size and shape of the binding cavity. (iii) W99 and M100, along with F255 and V296, play a key role for receptor activation, providing insights on bitter taste receptor activation not emerging from the previously reported computational models

Vascular smooth muscle cells remodel collagen matrices by long-distance action and anisotropic interaction

While matrix remodeling plays a key role in vascular physiology and pathology, the underlying mechanisms have remained incompletely understood. We studied the remodeling of collagen matrices by individual vascular smooth muscle cells (SMCs), clusters and monolayers. In addition, we focused on the contribution of transglutaminase 2 (TG2), which plays an important role in the remodeling of small arteries. Single SMCs displaced fibers in collagen matrices at distances up to at least 300 μm in the course of 8–12 h. This process involved both ‘hauling up’ of matrix by the cells and local matrix compaction at a distance from the cells, up to 200 μm. This exceeded the distance over which cellular protrusions were active, implicating the involvement of secreted enzymes such as TG2. SMC isolated from TG2 KO mice still showed compaction, with changed dynamics and relaxation. The TG active site inhibitor L682777 blocked local compaction by wild type cells, strongly reducing the displacement of matrix towards the cells. At increasing cell density, cells cooperated to establish compaction. In a ring-shaped collagen matrix, this resulted in preferential displacement in the radial direction, perpendicular to the cellular long axis. This process was unaffected by inhibition of TG2 cross-linking. These results show that SMCs are capable of matrix remodeling by prolonged, gradual compaction along their short axis. This process could add to the 3D organization and remodeling of blood vessels based on the orientation and contraction of SMCs

HIV-1 Nef Protein Structures Associated with Brain Infection and Dementia Pathogenesis

The difference between regional rates of HIV-associated dementia (HAD) in patients infected with different subtypes of HIV suggests that genetic determinants exist within HIV that influence the ability of the virus to replicate in the central nervous system (in Uganda, Africa, subtype D HAD rate is 89%, while subtype A HAD rate is 24%). HIV-1 nef is a multifunctional protein with known toxic effects in the brain compartment. The goal of the current study was to identify if specific three-dimensional nef structures may be linked to patients who developed HAD. HIV-1 nef structures were computationally derived for consensus brain and non-brain sequences from a panel of patients infected with subtype B who died due to varied disease pathologies and consensus subtype A and subtype D sequences from Uganda. Site directed mutation analysis identified signatures in brain structures that appear to change binding potentials and could affect folding conformations of brain-associated structures. Despite the large sequence variation between HIV subtypes, structural alignments confirmed that viral structures derived from patients with HAD were more similar to subtype D structures than to structures derived from patient sequences without HAD. Furthermore, structures derived from brain sequences of patients with HAD were more similar to subtype D structures than they were to their own non-brain structures. The potential finding of a brain-specific nef structure indicates that HAD may result from genetic alterations that alter the folding or binding potential of the protein

Parallel molecular routes to cold adaptation in eight genera of New Zealand stick insects.

The acquisition of physiological strategies to tolerate novel thermal conditions allows organisms to exploit new environments. As a result, thermal tolerance is a key determinant of the global distribution of biodiversity, yet the constraints on its evolution are not well understood. Here we investigate parallel evolution of cold tolerance in New Zealand stick insects, an endemic radiation containing three montane-occurring species. Using a phylogeny constructed from 274 orthologous genes, we show that stick insects have independently colonized montane environments at least twice. We compare supercooling point and survival of internal ice formation among ten species from eight genera, and identify both freeze tolerance and freeze avoidance in separate montane lineages. Freeze tolerance is also verified in both lowland and montane populations of a single, geographically widespread, species. Transcriptome sequencing following cold shock identifies a set of structural cuticular genes that are both differentially regulated and under positive sequence selection in each species. However, while cuticular proteins in general are associated with cold shock across the phylogeny, the specific genes at play differ among species. Thus, while processes related to cuticular structure are consistently associated with adaptation for cold, this may not be the consequence of shared ancestral genetic constraints

Genome-wide association analysis reveals QTL and candidate mutations involved in white spotting in cattle

International audienceAbstractBackgroundWhite spotting of the coat is a characteristic trait of various domestic species including cattle and other mammals. It is a hallmark of Holstein–Friesian cattle, and several previous studies have detected genetic loci with major effects for white spotting in animals with Holstein–Friesian ancestry. Here, our aim was to better understand the underlying genetic and molecular mechanisms of white spotting, by conducting the largest mapping study for this trait in cattle, to date.ResultsUsing imputed whole-genome sequence data, we conducted a genome-wide association analysis in 2973 mixed-breed cows and bulls. Highly significant quantitative trait loci (QTL) were found on chromosomes 6 and 22, highlighting the well-established coat color genes KIT and MITF as likely responsible for these effects. These results are in broad agreement with previous studies, although we also report a third significant QTL on chromosome 2 that appears to be novel. This signal maps immediately adjacent to the PAX3 gene, which encodes a known transcription factor that controls MITF expression and is the causal locus for white spotting in horses. More detailed examination of these loci revealed a candidate causal mutation in PAX3 (p.Thr424Met), and another candidate mutation (rs209784468) within a conserved element in intron 2 of MITF transcripts expressed in the skin. These analyses also revealed a mechanistic ambiguity at the chromosome 6 locus, where highly dispersed association signals suggested multiple or multiallelic QTL involving KIT and/or other genes in this region.ConclusionsOur findings extend those of previous studies that reported KIT as a likely causal gene for white spotting, and report novel associations between candidate causal mutations in both the MITF and PAX3 genes. The sizes of the effects of these QTL are substantial, and could be used to select animals with darker, or conversely whiter, coats depending on the desired characteristics

Differential Modulation of Beta-Adrenergic Receptor Signaling by Trace Amine-Associated Receptor 1 Agonists

Trace amine-associated receptors (TAAR) are rhodopsin-like G-protein-coupled receptors (GPCR). TAAR are involved in modulation of neuronal, cardiac and vascular functions and they are potentially linked with neurological disorders like schizophrenia and Parkinson's disease. Subtype TAAR1, the best characterized TAAR so far, is promiscuous for a wide set of ligands and is activated by trace amines tyramine (TYR), phenylethylamine (PEA), octopamine (OA), but also by thyronamines, dopamine, and psycho-active drugs. Unfortunately, effects of trace amines on signaling of the two homologous β-adrenergic receptors 1 (ADRB1) and 2 (ADRB2) have not been clarified yet in detail. We, therefore, tested TAAR1 agonists TYR, PEA and OA regarding their effects on ADRB1/2 signaling by co-stimulation studies. Surprisingly, trace amines TYR and PEA are partial allosteric antagonists at ADRB1/2, whereas OA is a partial orthosteric ADRB2-antagonist and ADRB1-agonist. To specify molecular reasons for TAAR1 ligand promiscuity and for observed differences in signaling effects on particular aminergic receptors we compared TAAR, tyramine (TAR) octopamine (OAR), ADRB1/2 and dopamine receptors at the structural level. We found especially for TAAR1 that the remarkable ligand promiscuity is likely based on high amino acid similarity in the ligand-binding region compared with further aminergic receptors. On the other hand few TAAR specific properties in the ligand-binding site might determine differences in ligand-induced effects compared to ADRB1/2. Taken together, this study points to molecular details of TAAR1-ligand promiscuity and identified specific trace amines as allosteric or orthosteric ligands of particular β-adrenergic receptor subtypes

Chemogenomic Analysis of G-Protein Coupled Receptors and Their Ligands Deciphers Locks and Keys Governing Diverse Aspects of Signalling

Understanding the molecular mechanism of signalling in the important super-family of G-protein-coupled receptors (GPCRs) is causally related to questions of how and where these receptors can be activated or inhibited. In this context, it is of great interest to unravel the common molecular features of GPCRs as well as those related to an active or inactive state or to subtype specific G-protein coupling. In our underlying chemogenomics study, we analyse for the first time the statistical link between the properties of G-protein-coupled receptors and GPCR ligands. The technique of mutual information (MI) is able to reveal statistical inter-dependence between variations in amino acid residues on the one hand and variations in ligand molecular descriptors on the other. Although this MI analysis uses novel information that differs from the results of known site-directed mutagenesis studies or published GPCR crystal structures, the method is capable of identifying the well-known common ligand binding region of GPCRs between the upper part of the seven transmembrane helices and the second extracellular loop. The analysis shows amino acid positions that are sensitive to either stimulating (agonistic) or inhibitory (antagonistic) ligand effects or both. It appears that amino acid positions for antagonistic and agonistic effects are both concentrated around the extracellular region, but selective agonistic effects are cumulated between transmembrane helices (TMHs) 2, 3, and ECL2, while selective residues for antagonistic effects are located at the top of helices 5 and 6. Above all, the MI analysis provides detailed indications about amino acids located in the transmembrane region of these receptors that determine G-protein signalling pathway preferences