Computational exploration of polymorphisms in 5-hydoxytryptamine 5-HT₁A and 5-HT₂A receptors associated with psychiatric disease

The huge polymorphic data have been prioritized towards a specific disease based on sequence and structure homology tools to a large extent. In this study, we have explored the potential non-synonymous Single Nucleotide Polymorphism (nsSNP) in serotonin (5-HT) receptors involved in psychotic syndromes and their response pathway. The most damaging point mutations were screened from 12 classes of serotonin receptors comprising 7743 variants. In 5HT(1A) receptor, two alleles were found to be highly deleterious located at ligand binding extracellular-2 and one at intracellular loop-3 domains. Similarly, we found two alleles predicted to be highly damaging in 5HT(2A) residing at N and C-Terminal domains. The above alleles were further confirmed based on their flexibility and stability difference using the molecular dynamic simulation analysis. Integrating these results appeared promising for being able to filter out potential non-synonymous Single Nucleotide Polymorphisms for neuropsychiatric disorders.     

Interactions of the Melanocortin-4 Receptor with the Peptide Agonist NDP-MSH

Melanocortin-4 receptor (MC4R) has an important regulatory role in energy homeostasis and food intake. Peptide agonists of the MC4R are characterized by the conserved sequence His₆-Phe₇-Arg₈-Trp₉, which is crucial for their interaction with the receptor. This investigation utilized the covalent attachment approach to identify receptor residues in close proximity to the bound ligand [Nle₄,d-Phe₇]melanocyte-stimulating hormone (NDP-MSH), thereby differentiating between residues directly involved in ligand binding and those mutations that compromise ligand binding by inducing conformational changes in the receptor. Also, recent X-ray structures of G-protein-coupled receptors were utilized to refine a model of human MC4R in the active state (R^?), which was used to generate a better understanding of the binding mode of the ligand NDP-MSH at the atomic level.The mutation of residues in the human MC4R—such as Leu106 of extracellular loop 1, and Asp122, Ile125, and Asp126 of transmembrane (TM) helix 3, His264 (TM6), and Met292 (TM7)—to Cys residues produced definitive indications of proximity to the side chains of residues in the core region of the peptide ligand. Of particular interest was the contact between d-Phe₇ on the ligand and Ile125 of TM3 on the MC4R. Additionally, Met292 (TM7) equivalent to Lys(7.45) (Ballesteros numbering scheme) involved in covalently attaching retinal in rhodopsin is shown to be in close proximity to Trp₉.For the first time, the interactions between the terminal regions of NDP-MSH and the receptor are described. The amino-terminus appears to be adjacent to a series of hydrophilic residues with novel interactions at Cys196 (TM5) and Asp189 (extracellular loop 2). These interactions are reminiscent of sequential ligand binding exhibited by the β₂-adrenergic receptor, with the former interaction being equivalent to the known interaction involving Ser204 of the β₂-adrenergic receptor.     

Single nucleotide polymorphisms in genes that are common targets of luteotropin and luteolysin in primate corpus luteum: Computational exploration

Luteal insufficiency affects fertility and hence study of mechanisms that regulate corpus luteum (CL) function is of prime importance to overcome infertility problems. Exploration of human genome sequence has helped to study the frequency of single nucleotide polymorphisms (SNPs). Clinical benefits of screening SNPs in infertility are being recognized well in recent times. Examining SNPs in genes associated with maintenance and regression of CL may help to understand unexplained luteal insufficiency and related infertility. Publicly available microarray gene expression databases reveal the global gene expression patterns in primate CL during the different functional state. We intend to explore computationally the deleterious SNPs of human genes reported to be common targets of luteolysin and luteotropin in primate CL. Different computational algorithms were used to dissect out the functional significance of SNPs in the luteinizing hormone sensitive genes. The results raise the possibility that screening for SNPs might be integrated to evaluate luteal insufficiency associated with human female infertility for future studies.     

Neuronal ceroid lipofuscinosis type CLN2: A new rationale for the construction of phenotypic subgroups based on a survey of 25 cases in South America

Tripeptidyl-peptidase 1 (TPP1) null or residual activity occurs in neuronal ceroid lipofuscinosis (NCL) with underlying TPP1/CLN2 mutations. A survey of 25 South American CLN2 affected individuals enabled the differentiation of two phenotypes: classical late-infantile and variant juvenile, each in approximately 50% of patients, with residual TPP1 activity occurring in approximately 32%. Each individual was assigned to one of three subgroups: (I) n = 11, null TPP1 activity in leukocytes; (II) n = 8, residual TPP1 activity of 0.60–15.85 nmol/h/mg (nr 110–476); (III) n = 6, activity not measured in leukocytes. Curvilinear bodies (CB) appeared in almost all studied CLN2 subjects; the only exceptions occurred in cases of subgroup II: two individuals had combined CBs/fingerprints (FPs), and one case had pure FPs. There were 15 mutations (4 first published in this paper, 3 previously observed in South America by our group, and 8 previously observed by others). In subgroup I, mutations were either missense or nonsense; in subgroups II and III, mutations prevailed at the non-conserved intronic site, c.887 − 10A>G (intron 7), and to a lesser extent at c.89 + 5G>C (intron 2), in heterozygous combinations. Grouping phenotypically and genetically known individuals on the basis of TPP1 activity supported the concept that residual enzyme activity underlies a protracted disease course. The prevalence of intronic mutations at non-conserved sites in subgroup II individuals indicates that some alternative splicing might allow some residual TPP1 activity.     

Exploration of deleterious single nucleotide polymorphisms in late-onset Alzheimer disease susceptibility genes

Non-synonymous single nucleotide polymorphisms (nsSNPs) are considered as biomarkers to disease susceptibility. In the present study, nsSNPs in CLU, PICALM and BIN1 genes were screened for their functional impact on concerned proteins and their plausible role in Alzheimer disease (AD) susceptibility. Initially, SNPs were retrieved from dbSNP database, followed by identification of potentially deleterious nsSNPs and prediction of their effect on proteins by PolyPhen and SIFT. Protein stability and the probability of mutation occurrence were predicted using I-Mutant and PANTHER respectively. SNPs3D and FASTSNP were used for the functional analysis of nsSNPs. The functional impact on the 3D structure of proteins was evaluated by SWISSPDB viewer and NOMAD-Ref server. On analysis, 3 nsSNPs with IDs rs12800974 (T158P) of PICALM and rs11554585 (R397C) and rs11554585 (N106D) of BIN1 were predicted to be functionally significant with higher scores of I-Mutant, SIFT, PolyPhen, PANTHER, FASTSNP and SNPs3D. The mutant models of these nsSNPs also showed very high energies and RMSD values compared to their native structures. Current study proposes that the three nsSNPs identified in this study constitute a unique resource of potential genetic factors for AD susceptibility.     

Systems biology approach for mutational and site-specific structural investigation of DNA repair genes for xeroderma pigmentosum

Xeroderma pigmentosum (XP) is a rare genetic skin disorder caused due to the extreme sensitivity for ultraviolet (UV) radiations. On its exposure, DNA acquires damages leading to skin and often neurological abnormalities. The DNA repair implicated in fixing UV-induced damages is NER and mutations in genes involved in NER and TLS form the basis of XP. The analyses of such mutations are vital for understanding XP and involved cancer genetics to facilitate the identification of crucial biomarkers and anticancer therapeutics. We detected the deleterious nsSNPs and examined them at structure-level by altering the structure, estimating secondary structure, solvent accessibility and performing site specific analysis. Crucial phosphorylation sites were also identified for their role in the disorder. These mutational and structural analyses offer valuable insight to the fundamental association of genetic mutations with phenotypic variations in XP and will assist experimental biologists to evaluate the mutations and their impact on genome.     

Evidence of colorectal cancer risk associated variant Lys25Ser in the proximity of human bone morphogenetic protein 2

Colorectal cancer (CRC) is the third most prevalent cancer and fourth leading cause of cancer-related deaths globally. It has been shown that the nsSNP variants play an important role in diseases, however it remained unclear how these variants are associated with the disease. Recently, several CRC risk associated SNPs have been discovered, however rs961253 (Lys25Arg at 20p12.3) located in the proximity of bone morphogenetic protein 2 (Bmp2) and fermitin family homolog 1 Fermt1 genes have been reported to be highly associated with the CRC risk. Here we provide evidence for the first time in silico biological functional and structural implications of non-synonymous (nsSNPs) CRC disease-associated variant Lys25Arg via molecular dynamic (MD) simulation. Protein structural analysis was performed with a particular variant allele (A/C, Lys25Arg) and compared with the predicted native protein structure. Our results showed that this nsSNP will cause changes in the protein structure and as a result is associated with the disease. In addition to the native and mutant 3D structures of CRC associated risk allele protein domain (CRAPD), they were also analyzed using solvent accessibility models for further protein stability confirmation. Taken together, this study confirmed that this variant has functional effect and structural impact on the CRAPD and may play an important role in CRC disease progression; hence it could be a reasonable approach for studying the effect of other deleterious variants in future studies.     

Structure of the Third Intracellular Loop of the Vasopressin V2 Receptor and Conformational Changes upon Binding to gC1qR

The V2 vasopressin receptor is a G-protein-coupled receptor that regulates the renal antidiuretic response. Its third intracellular loop is involved in the coupling not only with the GαS protein but also with gC1qR, a potential chaperone of G-protein-coupled receptors. In this report, we describe the NMR solution structure of the V2 i3 loop under a cyclized form (i3_cyc) and characterize its interaction with gC1qR. i3_cyc formed a left-twisted α-helical hairpin structure. The building of a model of the entire V2 receptor including the i3_cyc NMR structure clarified the side-chain orientation of charged residues, in agreement with literature mutagenesis reports. In the model, the i3 loop formed a rigid helical column, protruding deep inside the cytoplasm, as does the i3 loop in the recently elucidated structure of squid rhodopsin. However, its higher packing angle resulted in a different structural motif at the intracellular interface, which may be important for the specific recognition of GαS. Moreover, we could estimate the apparent K_d of the i3_cyc/gC1qR complex by anisotropy fluorescence. Using a shorter and more soluble version of i3_cyc, which encompassed the putative site of gC1qR binding, we showed by NMR saturation transfer difference spectroscopy that the binding surface corresponded to the central arginine cluster. Binding to gC1qR induced the folding of the otherwise disordered short peptide into a spiral-like path formed by a succession of I and IV turns. Our simulations suggested that this folding would rigidify the arginine cluster in the entire i3 loop and would alter the conformation of the cytosolic extensions of TM V and TM VI helices. In agreement with this conformational rearrangement, we observed that binding of gC1qR to the full-length receptor modifies the intrinsic tryptophan fluorescence binding curves of V2 to an antagonist.     

Congenital hyperinsulinism: Clinical and molecular analysis of a large Italian cohort

Congenital hyperinsulinism (CHI) is a genetic disorder characterized by profound hypoglycemia related to an inappropriate insulin secretion. It is a heterogeneous disease classified into two major subgroups: “channelopathies” due to defects in ATP-sensitive potassium channel, encoded by ABCC8 and KCNJ11 genes, and “metabolopathies” caused by mutation of several genes (GLUD1, GCK, HADH, SLC16A1, HNF4A and HNF1A) and involved in different metabolic pathways. To elucidate the genetic etiology of CHI in the Italian population, we conducted an extensive sequencing analysis of the CHI-related genes in a large cohort of 36 patients: Twenty-nine suffering from classic hyperinsulinism (HI) and seven from hyperinsulinism–hyperammonemia (HI/HA). Seventeen mutations have been found in fifteen HI patients and five mutations in five HI/HA patients. Our data confirm the major role of ATP-sensitive potassium channel in the pathogenesis of Italian cases (~ 70%) while the remaining percentage should be attributed to other. A better knowledge of molecular basis of CHI would lead to improve strategies for genetic screening and prenatal diagnosis. Moreover, genetic analysis might also help to distinguish the two histopathological forms of CHI, which would lead to a clear improvement in the treatment and in genetic counseling.     

Leigh syndrome associated with mitochondrial complex I deficiency due to novel mutations In NDUFV1 and NDUFS2

Leigh syndrome (LS) is a progressive neurodegenerative disease caused by either mitochondrial or nuclear DNA mutations resulting in dysfunctional mitochondrial energy metabolism. Mutations in genes encoding for subunits of the respiratory chain or assembly factors of respiratory chain complexes are often documented in LS cases. Nicotinamide adenine dinucleotide (NADH):ubiquinone oxidoreductase (complex I) enzyme deficiencies account for a significant proportion of mitochondrial disorders, including LS. In an attempt to expand the repertoire of known mutations accounting for LS, we describe the clinical, radiological, biochemical and molecular data of six patients with LS found to have novel mutations in two complex I subunits (NDUFV1 and NDUFS2). Two siblings were homozygous for the previously undescribed R386C mutation in NDUFV1, one patient was a compound heterozygote for the R386C mutation in NDUFV1 and a frameshift mutation in the same gene, one patient was a compound heterozygote for the R88G and R199P mutations in NDUFV1, and two siblings were compound heterozygotes for an undescribed E104A mutation in NDUFS2. After the novel mutations were identified, we employed prediction models using protein conservation analysis (SIFT, PolyPhen and UCSC genome browser) to determine pathogenicity. The R386C, R88G, R199P, and E104A mutations were found to be likely pathogenic, and thus presumably account for the LS phenotype. This case series broadens our understanding of the etiology of LS by identifying new molecular defects that can result in complex I deficiency and may assist in targeted diagnostics and/or prenatal diagnosis of LS in the future.     

Involvement of the second extracellular loop and transmembrane residues of CCR5 in inhibitor binding and HIV-1 fusion: insights into the mechanism of allosteric inhibition

C-C chemokine receptor 5 (CCR5), a member of G-protein-coupled receptors, serves as a coreceptor for human immunodeficiency virus type 1 (HIV-1). In the present study, we examined the interactions between CCR5 and novel CCR5 inhibitors containing the spirodiketopiperazine scaffolds AK530 and AK317, both of which were lodged in the hydrophobic cavity located between the upper transmembrane domain and the second extracellular loop (ECL2) of CCR5. Although substantial differences existed between the two inhibitors—AK530 had 10-fold-greater CCR5-binding affinity (K_d = 1.4 nM) than AK317 (16.7 nM)—their antiviral potencies were virtually identical (IC₅₀ = 2.1 nM and 1.5 nM, respectively). Molecular dynamics simulations for unbound CCR5 showed hydrogen bond interactions among transmembrane residues Y108, E283, and Y251, which were crucial for HIV-1-gp120/sCD4 complex binding and HIV-1 fusion. Indeed, AK530 and AK317, when bound to CCR5, disrupted these interhelix hydrogen bond interactions, a salient molecular mechanism enabling allosteric inhibition. Mutagenesis and structural analysis showed that ECL2 consists of a part of the hydrophobic cavity for both inhibitors, although AK317 is more tightly engaged with ECL2 than AK530, explaining their similar anti-HIV-1 potencies despite the difference in K_d values. We also found that amino acid residues in the β-hairpin structural motif of ECL2 are critical for HIV-1-elicited fusion and binding of the spirodiketopiperazine-based inhibitors to CCR5. The direct ECL2-engaging property of the inhibitors likely produces an ECL2 conformation, which HIV-1 gp120 cannot bind to, but also prohibits HIV-1 from utilizing the “inhibitor-bound” CCR5 for cellular entry—a mechanism of HIV-1's resistance to CCR5 inhibitors. The data should not only help delineate the dynamics of CCR5 following inhibitor binding but also aid in designing CCR5 inhibitors that are more potent against HIV-1 and prevent or delay the emergence of resistant HIV-1 variants.     

Genetic association of KCNJ10 rs1130183 with seizure susceptibility and computational analysis of deleterious non-synonymous SNPs of KCNJ10 gene

Establishing genetic basis of Idiopathic generalized epilepsies (IGE) is challenging because of their complex inheritance pattern and genetic heterogeneity. Kir4.1 inwardly rectifying channel (KCNJ10) is one of the independent genes reported to be associated with seizure susceptibility. In the current study we have performed a comprehensive in silico analysis of genetic variants in KCNJ10gene at functional and structural level along with a case–control analysis for the association ofrs1130183 (R271C) polymorphism in Indian patients with IGE. Age and sex matched 108epileptic patients and normal healthy controls were examined. Genotyping of KCNJ10rs1130183variation was performed using PCR-RFLP method. The risk association was determined by using odds ratio and 95% confidence interval. Functional effects of non-synonymous SNPs (nsSNPs) in KCNJ10 gene were analyzed using SIFT PolyPhen-2, I-Mutant 2.0, PANTHER and FASTSNP. Subsequently, homology modeling of protein three dimensional (3D) structures was performed using Modeller tool (9.10v) and compared the native protein with mutant for assessment of structure and stability. SIFT, PolyPhen-2, I-Mutant 2.0 and PANTHER collectively showed rs1130183, rs1130182 and rs137853073 SNPs inKCNJ10 gene affect protein structure and function. There was a considerable variation in the Root Mean Square Deviation (RMSD) value between the native and mutant structure (1.17?). Association analysis indicate KCNJ10rs1130183 did not contribute to risk of seizure susceptibility in Indian patients with IGE (OR- 0.38; 95%CI, 0.07–2.05) and T allele frequency (0.02%) was in concordance with dbSNP reports. This study identifies potential SNPs that may contribute to seizure susceptibility and further studies with the selected SNPs in larger number of samples and their functional analysis is required for understanding the variants of KCNJ10with seizure susceptibility.     

Exome sequencing identifies RDH12 compound heterozygous mutations in a family with severe retinitis pigmentosa

Objective

Retinitis pigmentosa (RP) is the most prevalent type of inherited retinal degeneration and one of the commonest causes of genetically determined visual dysfunction worldwide. To date, approximately 35 genes have been associated with nonsyndromic autosomal recessive RP (arRP), however the small contribution of each gene to the total prevalence of arRP and the lack of a clear genotype–phenotype correlation complicate the genetic analysis in affected patients. Next generation sequencing technologies are powerful and cost-effective methods for detecting causative mutations in both sporadic and familial RP cases.

Methods

A Mexican family with 5 members affected from arRP was studied. All patients underwent a complete ophthalmologic examination. Molecular methods included genome-wide SNP homozygosity mapping, exome sequencing analysis, and Sanger-sequencing confirmation of causal mutations.

Results

No regions of shared homozygosity among affected subjects were identified. Exome sequencing in a single patient allowed the detection of two missense mutations in the RDH12 gene: a c.446T>C transition predicting a novel p.L149P substitution, and a c.295C>A transversion predicting a previously reported p.L99I replacement. Sanger sequencing confirmed that all affected subjects carried both RDH12 mutations.

Conclusions

This study adds to the molecular spectrum of RDH12-related retinopathy and offers an additional example of the power of exome sequencing in the diagnosis of recessively inherited retinal degenerations.     

Crystal structure of the GluR0 ligand-binding core from Nostoc punctiforme in complex with L-glutamate: structural dissection of the ligand interaction and subunit interface

GluR0 from Nostoc punctiforme (NpGluR0) is a bacterial homologue of the ionotropic glutamate receptor (iGluR). We have solved the crystal structure of the ligand-binding core of NpGluR0 in complex with l-glutamate at a resolution of 2.1 Å. The structure exhibits a noncanonical ligand interaction and two distinct subunit interfaces. The side-chain guanidium group of Arg80 forms a salt bridge with the γ-carboxyl group of bound l-glutamate: in GluR0 from Synechocystis (SGluR0) and other iGluRs, the equivalent residues are Asn or Thr, which cannot form a similar interaction. We suggest that the local positively charged environment and the steric constraint created by Arg80 mediate the selectivity of l-glutamate binding by preventing the binding of positively charged and hydrophobic amino acids. In addition, the NpGluR0 ligand-binding core forms a new subunit interface in which the two protomers are arranged differently than the known iGluR and SGluR0 dimer interfaces. The significance of there being two different dimer interfaces was investigated using analytical ultracentrifugation analysis.     

Computational identification and analysis of functional polymorphisms involved in the activation and detoxification genes implicated in endometriosis

Endometriosis is a complex disorder of the female reproductive system where endometrial tissue embeds and grows at extrauterine location leading to inflammation and pain. Hundreds of polymorphisms in several genes have been studied as probable risk factors of this debilitating disease. Bioinformatics tools have come a long way in augmenting the search for putative functional polymorphisms in human diseases. In this study we have explored 16 genes involved in the detoxification of xenobiotic chemicals that are implicated in endometriosis by utilising publically available programs like SIFT, Polyphen, Panther, FastSNP, SNPeffect and PhosSNP. The variations among different ethnic populations of the SNPs were studied. We then calculated the extent to which bioinformatics based predictions are concurrent with real world epidemiological, genotyping studies using a set of SNPs that have been studied in endometriosis case–control studies. Our study shows that there is a significant positive correlation (r = 0.569, p < 0.005) between the summary of the predicted scores taken from 4 different servers and the odds ratio found from epidemiological studies. This report has identified and catalogued various deleterious SNPs that could be important in endometriosis and could aid in further analysis by in vitro and in vivo methods for the better understanding of the disease pathophysiology.     

Thermostabilization of the Neurotensin Receptor NTS1

Structural studies on G-protein-coupled receptors have been hampered for many years by their instability in detergent solution and by the number of potential conformations that receptors can adopt. Recently, the structures of the β₁ and β₂ adrenergic receptors and the adenosine A_2a receptor were determined in the antagonist-bound state, a receptor conformation that is thought to be more stable than the agonist-bound state. In contrast to these receptors, the neurotensin (NT) receptor NTS1 is much less stable in detergent solution. We have therefore used a systematic mutational approach coupled with activity assays to identify receptor mutants suitable for crystallization, both alone and in complex with the peptide agonist NT. The best receptor mutant NTS1-7m contained four point mutations. It showed increased stability compared to the wild-type receptor, in the absence of ligand, after solubilization with a variety of detergents. In addition, NTS1-7m bound to NT was more stable than unliganded NTS1-7m. Of the four thermostabilizing mutations, only one residue (A86L) is predicted to be in the lipid environment. In contrast, I260A appears to be buried within the transmembrane helix bundle, F342A may form a distant part of the putative ligand-binding site, whereas F358A is likely to be in a region that is important for receptor activation. NTS1-7m binds NT with a similar affinity for the wild-type receptor. However, agonist dissociation was slower, and NTS1-7m activated G-proteins poorly. The affinity of NTS1-7m for the antagonist SR48692 was also lower than that of the wild-type receptor. Thus, we have successfully stabilized NTS1 in an agonist-binding conformation that does not efficiently couple to G-proteins.