Ramu Elango

A Computational Protein Phenotype Prediction Approach to Analyze the Deleterious Mutations of Human MED12 Gene

Genetic mutations in MED12, a subunit of Mediator complex are seen in a broad spectrum of human diseases. However, the underlying basis of how these pathogenic mutations elicit protein phenotype changes in terms of 3D structure, stability and protein binding sites remains unknown. Therefore, we aimed to investigate the structural and functional impacts of MED12 mutations, using computational methods as an alternate to traditional in vivo and in vitro approaches. The MED12 gene mutations details and their corresponding clinical associations were collected from different databases and by text‐mining. Initially, diverse computational approaches were applied to categorize the different classes of mutations based on their deleterious impact to MED12. Then, protein structures for wild and mutant types built by integrative modeling were analyzed for structural divergence, solvent accessibility, stability, and functional interaction deformities. Finally, this study was able to identify that genetic mutations mapped to exon‐2 region, highly conserved LCEWAV and Catenin domains induce biochemically severe amino acid changes which alters the protein phenotype as well as the stability of MED12‐CYCC interactions. To better understand the deleterious nature of FS‐IDs and Indels, this study asserts the utility of computational screening based on their propensity towards non‐sense mediated decay. Current study findings may help to narrow down the number of MED12 mutations to be screened for mediator complex dysfunction associated genetic diseases. This study supports computational methods as a primary filter to verify the plausible impact of pathogenic mutations based on the perspective of evolution, expression and phenotype of proteins. J. Cell. Biochem. 117: 2023–2035, 2016.

In-silico analysis of inflammatory bowel disease (IBD) GWAS loci to novel connections.

Genome-wide association studies (GWASs) for many complex diseases, including inflammatory bowel disease (IBD), produced hundreds of disease-associated loci—the majority of which are noncoding. The number of GWAS loci is increasing very rapidly, but the process of translating single nucleotide polymorphisms (SNPs) from these loci to genomic medicine is lagging. In this study, we investigated 4,734 variants from 152 IBD associated GWAS loci (IBD associated 152 lead noncoding SNPs identified from pooled GWAS results + 4,582 variants in strong linkage-disequilibrium (LD) (r2 ≥0.8) for EUR population of 1K Genomes Project) using four publicly available bioinformatics tools, e.g. dbPSHP, CADD, GWAVA, and RegulomeDB, to annotate and prioritize putative regulatory variants. Of the 152 lead noncoding SNPs, around 11% are under strong negative selection (GERP++ RS ≥2); and ~30% are under balancing selection (Tajima’s D score >2) in CEU population (1K Genomes Project)—though these regions are positively selected (GERP++ RS <0) in mammalian evolution. The analysis of 4,734 variants using three integrative annotation tools produced 929 putative functional SNPs, of which 18 SNPs (from 15 GWAS loci) are in concordance with all three classifiers. These prioritized noncoding SNPs may contribute to IBD pathogenesis by dysregulating the expression of nearby genes. This study showed the usefulness of integrative annotation for prioritizing fewer functional variants from a large number of GWAS markers.

iPSC modeling of severe aplastic anemia reveals impaired differentiation and telomere shortening in blood progenitors.

Aplastic Anemia (AA) is a bone marrow failure (BMF) disorder, resulting in bone marrow hypocellularity and peripheral pancytopenia. Severe aplastic anemia (SAA) is a subset of AA defined by a more severe phenotype. Although the immunological nature of SAA pathogenesis is widely accepted, there is an increasing recognition of the role of dysfunctional hematopoietic stem cells in the disease phenotype. While pediatric SAA can be attributable to genetic causes, evidence is evolving on previously unrecognized genetic etiologies in a proportion of adults with SAA. Thus, there is an urgent need to better understand the pathophysiology of SAA, which will help to inform the course of disease progression and treatment options. We have derived induced pluripotent stem cell (iPSC) from three unaffected controls and three SAA patients and have shown that this in vitro model mimics two key features of the disease: (1) the failure to maintain telomere length during the reprogramming process and hematopoietic differentiation resulting in SAA-iPSC and iPSC-derived-hematopoietic progenitors with shorter telomeres than controls; (2) the impaired ability of SAA-iPSC-derived hematopoietic progenitors to give rise to erythroid and myeloid cells. While apoptosis and DNA damage response to replicative stress is similar between the control and SAA-iPSC-derived-hematopoietic progenitors, the latter show impaired proliferation which was not restored by eltrombopag, a drug which has been shown to restore hematopoiesis in SAA patients. Together, our data highlight the utility of patient specific iPSC in providing a disease model for SAA and predicting patient responses to various treatment modalities.

Comprehensive Computational Analysis of GWAS Loci Identifies CCR2 as a Candidate gene for Celiac Disease Pathogenesis

Celiac disease (CD) is a gluten intolerance disorder with known genetic contribution. The recent fine mapping and genome‐wide association studies (GWAS) have identified up to 57 non‐HLA CD susceptibility SNPs, majority of which are non‐coding variants lacking any functional annotation. Therefore, we adopted multidimensional computational approach for uncovering the plausible mechanisms through which these GWAS SNPs are connected to CD pathogenesis. At initial phase, we identified that 25 (43.85%) out of 57 CD‐SNPs lies in evolutionarily constrained genetic element regions. In follow‐up phases, through computational (CADD, GWAVA, and FATHMM algorithms) deleterious intensity measurements, we have discovered that 42 (3.94%) out of 1065 variants (57 CD‐lead and 1008‐linked SNPs; r2 ≥ 0.8) are differentially deleterious in nature to CD. Further functional scrutinization of these CD variants by public domain eQTL mapping, gene expression, knockout mouse model, and pathway analyses revealed that deleterious SNPs of CCR2 gene influences its expression levels and may also elicit a cascade of T‐cell‐mediated immunological events leading to intestinal gluten intolerance in genetically susceptible individuals. This study demonstrates the utility of integrated in silico analysis of annotations, gene expression, and pathways in prioritizing the potential complex disease variants from large‐scale open source genomic data. J. Cell. Biochem. 118: 2193–2207, 2017.

Replication of GWAS Coding SNPs Implicates MMEL1 as a Potential Susceptibility Locus among Saudi Arabian Celiac Disease Patients

Celiac disease (CD), a gluten intolerance disorder, was implicated to have 57 genetic susceptibility loci for Europeans but not for culturally and geographically distinct ethnic populations like Saudi Arabian CD patients. Therefore, we genotyped Saudi CD patients and healthy controls for three polymorphisms, that is, Phe196Ser in IRAK1, Trp262Arg in SH2B3, and Met518Thr in MMEL1 genes. Single locus analysis identified that carriers of the 518 Thr/Thr (MMEL1) genotype conferred a 1.6-fold increased disease risk compared to the noncarriers (OR = 2.6; 95% CI: 1.22–5.54; P < 0.01). This significance persisted even under allelic (OR = 1.55; 95% CI: 1.05–2.28; P = 0.02) and additive (OR = 0.35; 95% CI: 0.17–0.71; P = 0.03) genetic models. However, frequencies for Trp262Arg (SH2B3) and Phe196Ser (IRAK1) polymorphisms were not significantly different between patients and controls. The overall best MDR model included Met518Thr and Trp262Arg polymorphisms, with a maximal testing accuracy of 64.1% and a maximal cross-validation consistency of 10 out of 10 (P = 0.0156). Allelic distribution of the 518 Thr/Thr polymorphism in MMEL1 primarily suggests its independent and synergistic contribution towards CD susceptibility among Saudi patients. Lack of significant association of IRAK and SH2B3 gene polymorphisms in Saudi patients but their association in European groups suggests the genetic heterogeneity of CD.

Replication of GWAS loci revealed the moderate effect of TNRC6B locus on susceptibility of Saudi women to develop uterine leiomyomas

Uterine leiomyomas (UL) are smooth muscular nodes, whose growth is dependant up on the complex interplay of hormones with genes and uterine physiology. Global statistics indicate the role of ethnic and racial background as contributory factors for UL development. Owing to the lack of data, this study aimed to examine the association between genetic polymorphisms and susceptibility of Arab women of developing UL.

Induced pluripotent stem cell modelling of HLHS underlines the contribution of dysfunctional NOTCH signalling to impaired cardiogenesis

Abstract Hypoplastic left heart syndrome (HLHS) is among the most severe forms of congenital heart disease. Although the consensus view is that reduced flow through the left heart during development is a key factor in the development of the condition, the molecular mechanisms leading to hypoplasia of left heart structures are unknown. We have generated induced pluripotent stem cells (iPSC) from five HLHS patients and two unaffected controls, differentiated these to cardiomyocytes and identified reproducible in vitro cellular and functional correlates of the HLHS phenotype. Our data indicate that HLHS-iPSC have a reduced ability to give rise to mesodermal, cardiac progenitors and mature cardiomyocytes and an enhanced ability to differentiate to smooth muscle cells. HLHS-iPSC-derived cardiomyocytes are characterised by a lower beating rate, disorganised sarcomeres and sarcoplasmic reticulum and a blunted response to isoprenaline. Whole exome sequencing of HLHS fibroblasts identified deleterious variants in NOTCH receptors and other genes involved in the NOTCH signalling pathway. Our data indicate that the expression of NOTCH receptors was significantly downregulated in HLHS-iPSC-derived cardiomyocytes alongside NOTCH target genes confirming downregulation of NOTCH signalling activity. Activation of NOTCH signalling via addition of Jagged peptide ligand during the differentiation of HLHS-iPSC restored their cardiomyocyte differentiation capacity and beating rate and suppressed the smooth muscle cell formation. Together, our data provide firm evidence for involvement of NOTCH signalling in HLHS pathogenesis, reveal novel genetic insights important for HLHS pathology and shed new insights into the role of this pathway during human cardiac development.

Protein phenotype diagnosis of autosomal dominant calmodulin mutations causing irregular heart rhythms

The life‐threatening group of irregular cardiac rhythmic disorders also known as Cardiac Arrhythmias (CA) are caused by mutations in highly conserved Calmodulin (CALM/CaM) genes. Herein, we present a multidimensional approach to diagnose changes in phenotypic, stability, and Ca2+ ion binding properties of CA‐causing mutations. Mutation pathogenicity was determined by diverse computational machine learning approaches. We further modeled the mutations in 3D protein structure and analyzed residue level phenotype plasticity. We have also examined the influence of torsion angles, number of H‐bonds, and free energy dynamics on the stability, near‐native simulation dynamic potential of residue fluctuations in protein structures, Ca2+ ion binding potentials, of CaM mutants. Our study recomends to use M‐CAP method for measuring the pathogenicity of CA causing CaM variants. Interestingly, most CA‐causing variants we analyzed, exists in either third (V/H‐96, S/I‐98, V‐103) or fourth (G/V‐130, V/E/H‐132, H‐134, P‐136, G‐141, and L‐142) EF‐hands located in carboxyl domains of the CaM molecule. We observed that the minor structural fluctuations caused by these variants are likely tolerable owing to the highly flexible nature of calmodulins globular domains. However, our molecular docking results supports that these variants disturb the affinity of CaM toward Ca2+ ions and corroborate previous findings from functional studies. Taken together, these computational findings can explain the molecular reasons for subtle changes in structure, flexibility, and stability aspects of mutant CaM molecule. Our comprehensive molecular scanning approach demonstrates the utility of computational methods in quick preliminary screening of CA‐ CaM mutations before undertaking time consuming and complicated functional laboratory assays.

Whole exome sequencing of a consanguineous family identifies the possible modifying effect of a globally rare AK5 allelic variant in celiac disease development among Saudi patients

Celiac disease (CD), a multi-factorial auto-inflammatory disease of the small intestine, is known to occur in both sporadic and familial forms. Together HLA and Non-HLA genes can explain up to 50% of CD’s heritability. In order to discover the missing heritability due to rare variants, we have exome sequenced a consanguineous Saudi family presenting CD in an autosomal recessive (AR) pattern. We have identified a rare homozygous insertion c.1683_1684insATT, in the conserved coding region of AK5 gene that showed classical AR model segregation in this family. Sequence validation of 200 chromosomes each of sporadic CD cases and controls, revealed that this extremely rare (EXac MAF 0.000008) mutation is highly penetrant among general Saudi populations (MAF is 0.62). Genotype and allelic distribution analysis have indicated that this AK5 (c.1683_1684insATT) mutation is negatively selected among patient groups and positively selected in the control group, in whom it may modify the risk against CD development [p<0.002]. Our observation gains additional support from computational analysis which predicted that Iso561 insertion shifts the existing H-bonds between 400th and 556th amino acid residues lying near the functional domain of adenylate kinase. This shuffling of amino acids and their H-bond interactions is likely to disturb the secondary structure orientation of the polypeptide and induces the gain-of-function in nucleoside phosphate kinase activity of AK5, which may eventually down-regulates the reactivity potential of CD4+ T-cells against gluten antigens. Our study underlines the need to have population-specific genome databases to avoid false leads and to identify true candidate causal genes for the familial form of celiac disease.

Distribution of CYP2C8 and CYP2C9 Amino Acid Substitution Alleles in South Indian Diabetes Patients: A Genotypic and Computational Protein Phenotype Study

The CYP2C8 and CYP2C9 are two major isoforms of the cytochrome P450 enzyme family, which is involved in drug response, detoxification, and disease development. This study describes the differential distribution of amino acid substitution variants of CYP2C8 (*2‐I269F & *3‐R139K) and CYP2C9 (*2‐C144R & *3‐L359A) genes in 234 type 2 diabetes mellitus (T2DM) patients and 218 healthy controls from Andhra Pradesh, South India. Single locus genotype analysis has revealed that homozygous recessive genotypes of 2C8*2‐TT (P ≤ .03), 2C9*2‐TT (P ≤ .02), and heterozygous 2C9*3‐AC (P ≤ .006) are seen to be increasingly present in the case group, indicating a significant level of their association with diabetes in Andhra population. The statistical significance of these recessive genotypes has persisted even under their corresponding allelic forms (P ≤ .01). Genotype association results were further examined by computational protein structure and stability analysis to assess the deleteriousness of the amino acid changes. The mutant CYP 2C8 and 2C9 (both *2 and *3) proteins showed structural drifts at both amino acid residue (range 0.43Å‐0.77Å), and polypeptide chain levels (range 0.68Å‐1.81Å) compared to their wild‐type counterparts. Furthermore, the free energy value differences (range –0.915 to –1.38 Kcal/mol) between mutant and native protein structures suggests the deleterious and destabilizing potential of amino acid substitution polymorphisms of CYP genes. The present study confirms the variable distribution of CYP2C8 (*2 and *3) and CYP2C9 (*2 and *3) allelic polymorphisms among South Indian diabetic populations and further warrants the serious attention of CYP gene family, as a putative locus for disease risk assessment and therapy.