Linkage disequilibrium grouping of single nucleotide polymorphisms (SNPs) reflecting haplotype phylogeny for efficient selection of tag SNPs

Single nucleotide polymorphisms (SNPs) have been proposed to be grouped into haplotype blocks harboring a limited number of haplotypes. Within each block, the portion of haplotypes is expected to be tagged by a selected subset of SNPs; however, none of the proposed selection algorithms have been definitive. To address this issue, we developed a tag SNP selection algorithm based on grouping of SNPs by the linkage disequilibrium (LD) coefficient r(2) and examined five genes in three ethnic populations--the Japanese, African Americans, and Caucasians. Additionally, we investigated ethnic diversity by characterizing 979 SNPs distributed throughout the genome. Our algorithm could spare 60% of SNPs required for genotyping and limit the imprecision in allele-frequency estimation of nontag SNPs to 2% on average. We discovered the presence of a mosaic pattern of LD plots within a conventionally inferred haplotype block. This emerged because multiple groups of SNPs with strong intragroup LD were mingled in their physical positions. The pattern of LD plots showed some similarity, but the details of tag SNPs were not entirely concordant among three populations. Consequently, our algorithm utilizing LD grouping allows selection of a more faithful set of tag SNPs than do previous algorithms utilizing haplotype blocks.     

Biophysical defects in voltage.gated sodium channels associated with long QT and Brugada syndromes

The activity of voltage-gated sodium channels contributes to onset and duration of the cardiac action potential through an intricate balance with the activity of other ion channels. Activation of sodium channels leads to membrane depolarization and Phase 0 of the cardiac action potential. Sodium channel fast inactivation contributes to Phase 1, the initial repolarization. Slow inactivation and closed state fast inactivation determine channel availability and, thus, overall membrane excitability. Defects in any of these biophysical states or transitions between them, imparted by (over 170 reported thus far, including both Long QT3 and Brugada syndromes) mutations in the (over 2000) amino acids that compose the sodium channel protein, can lead to channel dysfunction that manifests as an abnormal cardiac action potential and electrocardiogram. A causal relationship between several such abnormalities and the panoply of sodium channel mutations have led to a greater understanding of the molecular underpinnings of cardiac arrhythmias as well as a deeper appreciation for the intricacies of sodium channel function. Here, we review the literature regarding these causal relationships from a perspective of the biophysical properties of sodium channels.     

In silico analysis of Single Nucleotide Polymorphisms (SNPs) in human BRAF gene

BRAF gene mutations are frequently seen in both inherited and somatic diseases. However, the harmful mutations for BRAF gene have not been predicted in silico. Owing to the importance of BRAF gene in cell division, differentiation and secretion processes, the functional analysis was carried out to explore the possible association between genetic mutations and phenotypic variations. Genomic analysis of BRAF was initiated with SIFT followed by PolyPhen and SNPs&GO servers to retrieve the 85 deleterious non-synonymous SNPs (nsSNPs) from dbSNP. A total of 5 mutations i.e. c.406T>G (S136A), c.1446G>T (R462I), c.1556 A>G (K499E), c.1860 T>A (V600E) and c.2352 C>T (P764L) that are found to exert benign effects on the BRAF protein structure and function were chosen for further analysis. Protein structural analysis with these amino acid variants was performed by using I-Mutant, FOLD-X, HOPE, NetSurfP, Swiss PDB viewer, Chimera and NOMAD-Ref servers to check their solvent accessibility, molecular dynamics and energy minimization calculations. Our in silico analysis suggested that S136A and P764L variants of BRAF could directly or indirectly destabilize the amino acid interactions and hydrogen bond networks thus explain the functional deviations of protein to some extent. Screening for BRAF, S136A and P764Lvariants may be useful for disease molecular diagnosis and also to design the molecular inhibitors of BRAF pathways.     

Phenotype guided characterization and molecular analysis of Indian patients with long QT syndromes

BackgroundLong QT syndromes (LQTS) are characterized by prolonged QTc interval on electrocardiogram (ECG) and manifest with syncope, seizures or sudden cardiac death. Long QT 1–3 constitute about 75% of all inherited LQTS. We classified a cohort of Indian patients for the common LQTS based on T wave morphology and triggering factors to prioritize the gene to be tested. We sought to identify the causative mutations and mutation spectrum, perform genotype-phenotype correlation and screen family members.MethodsThirty patients who fulfilled the criteria were enrolled. The most probable candidate gene among KCNQ1, KCNH2 and SCN5A were sequenced.ResultsOf the 30 patients, 22 were classified at LQT1, two as LQT2 and six as LQT3. Mutations in KCNQ1 were identified in 17 (77%) of 22 LQT1 patients, KCNH2 mutation in one of two LQT2 and SCN5A mutations in two of six LQT3 patients. We correlated the presence of the specific ECG morphology in all mutation positive cases. Eight mutations in KCNQ1 and one in SCN5A were novel and predicted to be pathogenic by in-silico analysis. Of all parents with heterozygous mutations, 24 (92%) of 26 were asymptomatic. Ten available siblings of nine probands were screened and three were homozygous and symptomatic, five heterozygous and asymptomatic.ConclusionsThis study in a cohort of Asian Indian patients highlights the mutation spectrum of common Long QT syndromes. The clinical utility for prevention of unexplained sudden cardiac deaths is an important sequel to identification of the mutation in at-risk family members.     

In silico structural and functional analysis of the human cytomegalovirus (HHV5) genome

The open reading frames of human cytomegalovirus (human herpesvirus-5, HHV5) encode some 213 unique proteins with mostly unknown functions. Using the threading program, ProCeryon, we calculated possible matches between the amino acid sequences of these proteins and the Protein Data Bank library of three-dimensional structures. Thirty-six proteins were fully identified in terms of their structure and, often, function; 65 proteins were recognized as members of narrow structural/functional families (e.g. DNA-binding factors, cytokines, enzymes, signaling particles, cell surface receptors etc.); and 87 proteins were assigned to broad structural classes (e.g. all-beta, 3-layer-alphabetaalpha, multidomain, etc.). Genes encoding proteins with similar folds, or containing identical structural traits (extreme sequence length, runs of unstructured (Pro and/or Gly-rich) residues, transmembrane segments, etc.) often formed tandem clusters throughout the genome. In the course of this work, benchmarks on about 20 known folds were used to optimize adjustable parameters of threading calculations, i.e. gap penalty weights used in sequence/structure alignments; new scores obtained as simple combinations of existing scoring functions; and number of threading runs conducive to meaningful results. An introduction of summed, per-residue-normalized scores has been essential for discovery of subdomains (EGF-like, SH2, SH3) in longer protein sequences, such as the eight "open sandwich" cytokine domains, 60-70 amino acids long and having the 3beta1alpha fold with one or two disulfide bridges, present in otherwise unrelated proteins.     

In silico screening of deleterious single nucleotide polymorphisms (SNPs) and molecular dynamics simulation of disease associated mutations in gene responsible for oculocutaneous albinism type 6 (OCA 6) disorder

Solute carrier family 24 member 5 (SLC24A5) is a gene that is associated with oculocutaneous albinism type 6 (OCA6) disorder and is involved in skin and hair pigmentation. It is involved in the maturation of melanosomes and melanin synthesis. SLC24A5 gene is located in the chromosomal position of 15q21.1. The present study involves the use of computational techniques in order to obtain a detailed picture of the most probable mutations that are associated with SLC24A5. From the observed result it was found that the mutation S145F is most deleterious and disease associated is predicted using several bioinformatics tools. The 3-D structures of native and mutant (S145F) were modeled in order to understand protein functionality using ab initio Robetta server. The modeled structure validation was done with ERRAT, Verify-3D, Procheck and RAMPAGE Ramachandran plot analysis. The most validated structure undergoes molecular dynamics simulations (MDS) study to understand the structural and functional behaviour of the native and mutant proteins. The MDS result showed the more flexibility in the native SLC24A5 structure. Due to mutation in the SLC24A5 protein structure it became more rigid and might disturb the conformational changes and glycosylation function of protein structure and might play role in inducing the OCA6. This study provides a significant insight into the underlying molecular mechanism involved in albinism associated with OCA6. It further helps scientists to develop a drug therapy against OCA 6 disease.

Communicated by Ramaswamy H. Sarma     

In silico studies on bacterial xylanase enzyme: Structural and functional insight

Xylans are the second most abundant form of hemicelluloses and are the second most abundant polysaccharide in nature after cellulose. To degrade xylan, microbes produce mainly xylanase enzyme. Wide range of microorganisms like fungi, bacteria, yeast, marine algae etc. are capable of producing xylanase. Main source of xylanase is fungi but industrial production of bacterial xylanase is low cost, easy downstream process and high production rate. To understand primary, secondary and tertiary structure of xylanase, in silico composition of amino acids, basic physiological characteristics; viz., pI, molecular weight, instability index, GRAVY, molar extinction coefficient, secondary structure, presence of functional domain and motifs, phylogenetic tree, salt bridge compositions are determined. In silico study of xylanase focused on 36 different bacterial sources are performed by retrieving FASTA and PDB sequences using RCSB PDB. FASTA and PDB files are proceed further in ExPASy-ProtParam, RAMPAGE, QMEAN, MEME, PSIPRED, InterProScan, MOTIF scan, ERRAT, Peptide cutter, ESBRI and MEGA 7. The instability index range (16.90–38.78) clearly indicates that the protein is highly stable. α-helix mean value (27.11%) infers the protein is dominated by α-helix region. The aliphatic index (39.80–90.68) gives information that the protein is highly thermostable, prevalence by alanine amino acid in aliphatic side chain. No transmembrane domain was found in the protein which confirms the enzyme is extracellular in nature. Ancestor chart analysis confirmed that it is a part of carbohydrate metabolic process and more specifically a member of glycoside hydrolase super family.     

In vitro and In silico studies on inhibitory effects of curcumin on multi drug resistance associated protein (MRP1) in retinoblastoma cells

Multi Drug Resistance (MDR) is one of the major causes of chemotherapy failure in human malignancies. Curcumin, the active constituent of Curcuma longa is a proven anticancer agent potentially modulating the expression and function of these MDR proteins. In this study, we attempted to test curcumin for its potential to inhibit the expression and function of multidrug resistance associated protein 1 (MRP1) in retinoblastoma (RB) cell lines through western blot, RT-PCR and functional assays. In silico analysis were also performed to understand the molecular interactions conferred by curucmin on MRP1 in RB cells. Western blot and RTPCR analysis did not show any correlation of MRP1 expression with increase in concentration of curcumin. However, inhibitory effect of curcumin on MRP1 function was observed as a decrease in the efflux of fluorescent substrate. Moreover, Curcumin did not affect 8-azido-ATP-biotin binding to MRP1 and it also showed inhibition of ATP-hydrolysis stimulated by quercetin, which is indicative of curcumin's interaction with the substrate binding site of MRP1. Furthermore, homology modelling and docking simulation studies of MRP1 also provided deeper insights into the molecular interactions, thereby inferring the potential binding mode of curcumin into the substrate binding site of MRP1.     

In vivo gene transfer of Kv1.5 normalizes action potential duration and shortens QT interval in mice with long QT phenotype

Mutations in cardiac voltage-gated K+ channels cause long QT syndrome (LQTS) and sudden death. We created a transgenic mouse with a long QT phenotype (Kv1DN) by overexpression of a truncated K+ channel in the heart and investigated whether the dominant negative effect of the transgene would be overcome by the direct injection of adenoviral vectors expressing wild-type Kv1.5 (AV-Kv1.5) into the myocardium. End points at 3-10 days included electrophysiology in isolated cardiomyocytes, surface ECG, programmed stimulation of the right ventricle, and in vivo optical mapping of action potentials and repolarization gradients in Langendorff-perfused hearts. Overexpression of Kv1.5 reconstituted a 4-aminopyridine-sensitive outward K+ current, shortened the action potential duration, eliminated early afterdepolarizations, shortened the QT interval, decreased dispersion of repolarization, and increased the heart rate. Each of these changes is consistent with a physiologically significant primary effect of adenoviral expression of Kv1.5 on ventricular repolarization of Kv1DN mice.     

TLR4 single nucleotide polymorphisms (SNPs) associated with Salmonella shedding in pigs

Toll-like receptor 4 (TLR4) is a key factor in the innate immune recognition of lipopolysaccharide (LPS) from Gram-negative bacteria. Previous studies from our group identified differences in the expression profile of TLR4 and genes affected by the TLR4 signaling pathway among pigs that shed varying levels of Salmonella, a Gram-negative bacterium. Therefore, genetic variation in this gene may be involved with the host’s immune response to bacterial infections. The current study screened for single nucleotide polymorphisms (SNPs) in the TLR4 gene and tested their association with Salmonella fecal shedding. Pigs (n?=?117) were intranasally challenged at 7 weeks of age with 1?×?10⁹ CFU of S. Typhimurium χ4232 and were classified as low or persistent Salmonella shedders based on the levels of Salmonella being excreted in fecal material. Salmonella fecal shedding was determined by quantitative bacteriology on days 2, 7, 14, and 20/21 post exposure, and the cumulative levels of Salmonella were calculated to identify the low (n?=?20) and persistent (n?=?20) Salmonella shedder pigs. From those 40 animals, the TLR4 region was sequenced, and 18 single nucleotide polymorphisms (SNPs) in TLR4 were identified. Twelve SNPs have been previously described and six are novel SNPs of which five are in the 5′ untranslated region and one is in intron 2. Single marker association test identified 13 SNPs associated with the qualitative trait of Salmonella fecal shedding, and seven of those SNPs were also associated with a quantitative measurement of fecal shedding (P?<?0.05). Using a stepwise regression process, a haplotype composed of SNPs rs80787918 and rs80907449 (P?≤?4.0?×?10^?3) spanning a region of 4.9 Kb was identified, thereby providing additional information of the influence of those SNPs on Salmonella fecal shedding in pigs.     

Divergent enhancer haplotype of ebony on inversion In(3R)Payne associated with pigmentation variation in a tropical population of Drosophila melanogaster

The pattern and intensity of pigmentation have direct impact on individual fitness through various ecological factors. In a Drosophila melanogaster population from southern Japan, thoracic trident pigmentation intensity of most of the strains could be classified into Dark or Light‐type. The expression level variation of the ebony gene correlated well with this phenotype and the allelic differences in expression indicated that the variation is partly due to cis‐regulatory changes. In the ～13 kb gene region, we identified 17 nucleotide sites and 2 indels that were in complete association with the thoracic trident pigmentation intensity. Interestingly, 11 out of 19 sites located within ～0.5 kb of the core epidermis enhancer. These sites had no obvious association with the abdominal pigmentation intensity in the previously analysed African populations from Uganda and Kenya, which suggested that multiple potential mutational pathways in the cis‐regulatory control region of a single gene could lead to similar phenotypic variation within this species. We also found that the Light‐type enhancer haplotype is strongly linked to a cosmopolitan inversion, In(3R)Payne, which is predominant in warmer climatic regions in both hemispheres. The sequence pattern suggested that the strong linkage may be due to selective forces related to thermal adaptation. The inferred selection for lighter pigmentation in the Japanese population is in the opposite direction of the previously reported case of selection for darker individuals in African populations. Nevertheless, both adaptive changes involved cis‐regulatory changes of ebony, which shows that this gene is likely to be a common target of natural selection.     

Large-scale single-nucleotide polymorphism (SNP) and haplotype analyses, using dense SNP Maps, of 199 drug-related genes in 752 subjects: the analysis of the association between uncommon SNPs within haplotype blocks and the haplotypes constructed with haplotype-tagging SNPs

To optimize the strategies for population-based pharmacogenetic studies, we extensively analyzed single-nucleotide polymorphisms (SNPs) and haplotypes in 199 drug-related genes, through use of 4,190 SNPs in 752 control subjects. Drug-related genes, like other genes, have a haplotype-block structure, and a few haplotype-tagging SNPs (htSNPs) could represent most of the major haplotypes constructed with common SNPs in a block. Because our data included 860 uncommon (frequency <0.1) SNPs with frequencies that were accurately estimated, we analyzed the relationship between haplotypes and uncommon SNPs within the blocks (549 SNPs). We inferred haplotype frequencies through use of the data from all htSNPs and one of the uncommon SNPs within a block and calculated four joint probabilities for the haplotypes. We show that, irrespective of the minor-allele frequency of an uncommon SNP, the majority (mean +/- SD frequency 0.943+/-0.117) of the minor alleles were assigned to a single haplotype tagged by htSNPs if the uncommon SNP was within the block. These results support the hypothesis that recombinations occur only infrequently within blocks. The proportion of a single haplotype tagged by htSNPs to which the minor alleles of an uncommon SNP were assigned was positively correlated with the minor-allele frequency when the frequency was <0.03 (P<.000001; n=233 [Spearman's rank correlation coefficient]). The results of simulation studies suggested that haplotype analysis using htSNPs may be useful in the detection of uncommon SNPs associated with phenotypes if the frequencies of the SNPs are higher in affected than in control populations, the SNPs are within the blocks, and the frequencies of the SNPs are >0.03.     

In silico analysis of non-synonymous missense SNPs (nsSNPs) in CPE,GNAS genes and experimental validation in type II diabetes mellitus through Next Generation Sequencing

Non-synonymous missense SNPs (nsSNPs) in CPE and GNAS genes were investigated computationally. In silico identified nsSNPs were experimentally validated in type II diabetes mellitus (T2DM) in Pakistani Pathan population using next generation sequencing (NGS).Sixty two high-risk nsSNPs in CPE and 44 in GNAS were identified. Only 12 in GNAS were clinically significant. Thirty six high-risk nsSNPs in CPE and 08 clinically significant nsSNPs in GNAS lies in the most conserved regions. I-mutant predicted that nsSNPs decrease the proteins stability and ModPred predicted 20 and 12 post-translational modification sites in CPE and GNAS proteins respectively. Ramachandran plot showed 88.7% residues are in the most favored region of protein models. By experimentation, none of the nsSNPs were found to be associated with T2DM.In conclusion, this study differentiates the deleterious nsSNPs from the neutral ones. Although nsSNPs are not associated with T2DM, they can be targeted in other CPE and GNAS genes related disorders.     

Preparation, functional characterization, and NMR studies of human KCNE1, a voltage-gated potassium channel accessory subunit associated with deafness and long QT syndrome

KCNE1, also known as minK, is a member of the KCNE family of membrane proteins that modulate the function of KCNQ1 and certain other voltage-gated potassium channels (KV). Mutations in human KCNE1 cause congenital deafness and congenital long QT syndrome, an inherited predisposition to potentially life-threatening cardiac arrhythmias. Although its modulation of KCNQ1 function has been extensively characterized, many questions remain regarding KCNE1's structure and location within the channel complex. In this study, KCNE1 was overexpressed in Escherichia coli and purified. Micellar solutions of the protein were then microinjected into Xenopus oocytes expressing KCNQ1 channels, followed by electrophysiological recordings aimed at testing whether recombinant KCNE1 can co-assemble with the channel. Nativelike modulation of channel properties was observed following injection of KCNE1 in lyso-myristoylphosphatidylglycerol (LMPG) micelles, indicating that KCNE1 is not irreversibly misfolded and that LMPG is able to act as a vehicle for delivering membrane proteins into the membranes of viable cells. 1H-15N TROSY NMR experiments indicated that LMPG micelles are well-suited for structural studies of KCNE1, leading to assignment of its backbone resonances and to relaxation studies. The chemical shift data confirmed that KCNE1's secondary structure includes several alpha-helices and demonstrated that its distal C-terminus is disordered. Surprisingly, for KCNE1 in LMPG micelles, there appears to be a break in alpha-helicity at sites 59-61, near the middle of the transmembrane segment, a feature that is accompanied by increased local backbone mobility. Given that this segment overlaps with sites 57-59, which are known to play a critical role in modulating KCNQ1 channel activation kinetics, this unusual structural feature likely has considerable functional relevance.     

Vibrational (FT-IR and FT-Raman) spectral investigations of 7-aminoflavone with density functional theoretical simulations

FT-Raman and FT-IR spectra of the 7-aminoflavone have been recorded and analysed. The detailed interpretation of the vibrational spectra has been carried out with the aid of normal coordinate analysis following the scaled quantum mechanical force field methodology. The various intramolecular interactions that are responsible for stabilisation of the molecule were revealed by natural bond orbital analysis. The obtained vibrational wavenumbers and optimised geometric parameters were observed to be in good agreement with the experimental data. The carbonyl stretching vibrations have been lowered due to conjugation and hydrogen bonding in the molecules.     

In silico analyses of penicillin binding proteins in Burkholderia pseudomallei uncovers SNPs with utility for phylogeography,species differentiation,and sequence typing

Burkholderia pseudomallei causes melioidosis. Sequence typing this pathogen can reveal geographical origin and uncover epidemiological associations. Here, we describe B. pseudomallei genes encoding putative penicillin binding proteins (PBPs) and investigate their utility for determining phylogeography and differentiating closely related species. We performed in silico analysis to characterize 10 PBP homologs in B. pseudomallei 1026b. As PBP active site mutations can confer β-lactam resistance in Gram-negative bacteria, PBP sequences in two resistant B. pseudomallei strains were examined for similar alterations. Sequence alignments revealed single amino acid polymorphisms (SAAPs) unique to the multidrug resistant strain Bp1651 in the transpeptidase domains of two PBPs, but not directly within the active sites. Using BLASTn analyses of complete assembled genomes in the NCBI database, we determined genes encoding PBPs were conserved among B. pseudomallei (n = 101) and Burkholderia mallei (n = 26) strains. Within these genes, single nucleotide polymorphisms (SNPs) useful for predicting geographic origin of B. pseudomallei were uncovered. SNPs unique to B. mallei were also identified. Based on 11 SNPs identified in two genes encoding predicted PBP-3s, a dual-locus sequence typing (DLST) scheme was developed. The robustness of this typing scheme was assessed using 1,523 RefSeq genomes from B. pseudomallei (n = 1,442) and B. mallei (n = 81) strains, resulting in 32 sequence types (STs). Compared to multi-locus sequence typing (MLST), the DLST scheme demonstrated less resolution to support the continental separation of Australian B. pseudomallei strains. However, several STs were unique to strains originating from a specific country or region. The phylogeography of Western Hemisphere B. pseudomallei strains was more highly resolved by DLST compared to internal transcribed spacer (ITS) typing, and all B. mallei strains formed a single ST. Conserved genes encoding PBPs in B. pseudomallei are useful for strain typing, can enhance predictions of geographic origin, and differentiate strains of closely related Burkholderia species.     

Computational and artificial neural network based study of functional SNPs of human LEPR protein associated with reproductive function

Genetic polymorphisms are mostly associated with inherited diseases, detecting and analyzing the biological significance of functional single-nucleotide polymorphisms (SNPs) using wet laboratory experiments is an arduous task hence the computational analysis of putative SNPs is essential before conducting a study on a large population. SNP in the leptin receptor (LEPR) could result in the retention of intracellular signalling due to the structural and functional instability of the receptor causing abnormal reproductive function in human. In this first comprehensive computational analysis of LEPR gene mutation, we have identified and analyzed the functional consequence and structural significance of the SNPs in LEPR using recently developed several computational algorithms. Thirteen deleterious mutations such as W13C, S93G, I232R, Q307H, Y354C, E497A, Q571H, R612H, K656N, T690A, T699M V741M, and L760R were identified in the LEPR gene coding region. Backpropagation algorithm has been developed to forestall the deleterious nature of SNP and to validate the outcome of the tested computational tools. From ConSurf prediction three SNPs (Q571H, R612H, and T699M) were highly conserved on LEPR protein and the most deleterious variant R612H had one hydrogen bond abolished and severely reduced protein stability. Molecular docking suggested that the mutant (R612H) LEPR had lowest binding energy than native LEPR with the ligand molecule. Thus the energetically destructive changeover of ARG to HIS in R612H could possibly affect the LEPR protein structural stability and functional constancy due to interruption in the amino acid interactions and could result in reproductive disorders in human and increases the complication in obstetric and pregnancy outcome.