Identification and in silico analysis of functional SNPs of the BRCA1 gene

Single-nucleotide polymorphisms (SNPs) play a major role in the understanding of the genetic basis of many complex human diseases. Also, the genetics of human phenotype variation could be understood by knowing the functions of these SNPs. It is still a major challenge to identify the functional SNPs in a disease-related gene. In this work, we have analyzed the genetic variation that can alter the expression and the function of the BRCA1 gene using computational methods. Of the total 477 SNPs, 65 were found to be nonsynonymous (ns) SNPs. Among the 14 SNPs in the untranslated region, 4 were found in the 5' and 10 were found in the 3' untranslated region (UTR). It was found that 16.9% of the nsSNPs were damaging, by both the SIFT and the PolyPhen servers. The UTR Resource tool suggested that 2 of 4 SNPs in the 5' UTR and 3 of 10 SNPs in the 3' UTR might change the protein expression levels. We identified major mutations from proline to serine at positions 1776 and 1812 of the native protein of the BRCA1 gene. From a comparison of the stabilizing residues of the native and mutant proteins, we propose that an nsSNP (rs1800751) could be an important candidate for the breast cancer caused by the BRCA1 gene.     

Computational detection of deleterious SNPs and their effect on sequence and structural level of the <Emphasis Type="Italic">VHL</Emphasis> gene

In this work we have analyzed the genetic variation that can alter the expression and the function of the VHL gene using computational methods. Of 110 single nucleotide polymorphisms (SNPs), 33 were found to be nonsynonymous (nsSNPs) and 23 SNPs were found in untranslated regions. Of the 33 nsSNPs investigated, 36.3% were found to be deleterious by both SIFT and PolyPhen servers. An untranslated region (UTR) resource tool suggested that two SNPs in the 5' UTR region and six SNPs in the 3' UTR region might change the protein expression levels. It was found by both SIFT and PolyPhen servers that a mutation from histidine to arginine at position 115 of the native protein of the VHL gene was most deleterious. A structural analysis of this mutated protein and the native protein was performed and had a root mean square deviation (RMSD) of 2.78 A. Based on this work, we propose that the nsSNP with a SNPid of rs5030812 is an important candidate for the cause of von Hippel-Lindau syndrome via the VHL gene.     

Applications of computational algorithm tools to identify functional SNPs in cytokine genes

Understanding the functions of single nucleotide polymorphisms (SNPs) can greatly help to understand the genetics of the human phenotype variation and especially the genetic basis of human complex diseases. However, how to identify functional SNPs from a pool containing both functional and neutral SNPs is challenging. In this study, we analyzed the genetic variations that can alter the expression and function of a group of cytokine proteins using computational tools. As a result, we extracted 4552 SNPs from 45 cytokine proteins from SNPper database. Of particular interest, 828 SNPs were in the 5'UTR region, 961 SNPs were in the 3' UTR region, and 85 SNPs were non-synonymous SNPs (nsSNPs), which cause amino acid change. Evolutionary conservation analysis using the SIFT tool suggested that 8 nsSNPs may disrupt the protein function. Protein structure analysis using the PolyPhen tool suggested that 5 nsSNPs might alter protein structure. Binding motif analysis using the UTResource tool suggested that 27 SNPs in 5' or 3'UTR might change protein expression levels. Our study demonstrates the presence of naturally occurring genetic variations in the cytokine proteins that may affect their expressions and functions with possible roles in complex human disease, such as immune diseases.     

Computational and Structural Investigation of Deleterious Functional SNPs in Breast Cancer BRCA2 Gene

In this work, we have analyzed the genetic variation that can alter the expression and the function in BRCA2 gene using computational methods. Out of the total 534 SNPs, 101 were found to be non synonymous (nsSNPs). Among the 7 SNPs in the untranslated region, 3 SNPs were found in 5′ and 4 SNPs were found in 3′ un-translated regions (UTR). Of the nsSNPs 20.7% were found to be damaging by both SIFT and PolyPhen server among the 101 nsSNPs investigated. UTR resource tool suggested that 2 SNPs in the 5′ UTR region and 4 SNPs in the 3′ UTR regions might change the protein expression levels. The mutation from asparagine to isoleucine at the position 3124 of the native protein of BRCA2 gene was most deleterious by both SIFT and PolyPhen servers. A structural analysis of this mutated protein and the native protein was made which had an RMSD value of 0.301 nm. Based on this work, we proposed that this most deleterious nsSNP with an SNPid rs28897759 is an important candidate for the cause of breast cancer by BRCA2 gene.     

Folding of bovine pancreatic trypsin inhibitor (BPTI) variants in which almost half the residues are alanine

Recent studies indicate that a fraction of the information contained in an amino acid sequence may be sufficient for specifying a native protein structure. An earlier alanine-scanning experiment conducted on bovine pancreatic trypsin inhibitor (BPTI; 58 residues) suggested that if cumulative mutations have additive effects on protein stability, a native protein structure could be built from BPTI sequences that contained many alanine residues distributed throughout the protein. To test this hypothesis, we designed and produced six BPTI mutants containing from 21 to 29 alanine residues. We found that the melting temperature of mutants containing up to 27 alanine residues (48 % of the total number of residues) could be predicted quite well by the sum of the change in melting temperature for the single mutations. Additionally, these same mutants folded into a native-like structure, as judged by their cooperative thermal denaturation curves and heteronuclear multiple quantum correlation (HMQC) NMR spectra. A BPTI mutant containing 22 alanine residues was further shown by 2D and 3D-NMR to fold into a structure very similar to that of native BPTI, and to be a functional trypsin inhibitor. These results provide insight into the extent to which native protein structure and function can be achieved with a highly simplified amino acid sequence.     

In silico analysis of single nucleotide polymorphism (SNPs) in human β-globin gene

Single amino acid substitutions in the globin chain are the most common forms of genetic variations that produce hemoglobinopathies--the most widespread inherited disorders worldwide. Several hemoglobinopathies result from homozygosity or compound heterozygosity to beta-globin (HBB) gene mutations, such as that producing sickle cell hemoglobin (HbS), HbC, HbD and HbE. Several of these mutations are deleterious and result in moderate to severe hemolytic anemia, with associated complications, requiring lifelong care and management. Even though many hemoglobinopathies result from single amino acid changes producing similar structural abnormalities, there are functional differences in the generated variants. Using in silico methods, we examined the genetic variations that can alter the expression and function of the HBB gene. Using a sequence homology-based Sorting Intolerant from Tolerant (SIFT) server we have searched for the SNPs, which showed that 200 (80%) non-synonymous polymorphism were found to be deleterious. The structure-based method via PolyPhen server indicated that 135 (40%) non-synonymous polymorphism may modify protein function and structure. The Pupa Suite software showed that the SNPs will have a phenotypic consequence on the structure and function of the altered protein. Structure analysis was performed on the key mutations that occur in the native protein coded by the HBB gene that causes hemoglobinopathies such as: HbC (E→K), HbD (E→Q), HbE (E→K) and HbS (E→V). Atomic Non-Local Environment Assessment (ANOLEA), Yet Another Scientific Artificial Reality Application (YASARA), CHARMM-GUI webserver for macromolecular dynamics and mechanics, and Normal Mode Analysis, Deformation and Refinement (NOMAD-Ref) of Gromacs server were used to perform molecular dynamics simulations and energy minimization calculations on β-Chain residue of the HBB gene before and after mutation. Furthermore, in the native and altered protein models, amino acid residues were determined and secondary structures were observed for solvent accessibility to confirm the protein stability. The functional study in this investigation may be a good model for additional future studies.     

内含肽介导的氯离子通道蛋白CFTR的反式剪接

研究利用内含肽(intein)的蛋白质反式剪接功能在大肠杆菌中对囊性纤维化跨膜传导调节因子(cystic fibrosis transmembrane regulator, CFTR)的反式剪接作用.CFTR基因突变导致一种常染色体隐性遗传疾病囊性纤维化(cystic fibrosis, CF).将CFTR的cDNA于剪接反应所需的保守性氨基酸残基Ser-660前断裂为N端和C端,分别与split mini Ssp DnaB 内含肽的106个氨基酸残基的N端和48个氨基酸残基的C端编码序列融合,构建到原核表达载体pBV220 诱导表达后SDS-PAGE可见预期大小剪接形成的CFTR蛋白条带,Western印迹用CFTR特异性抗体进一步证明为剪接所产生的CFTR蛋白,表明内含肽可有效催化CFTR的反式剪接.     

Divergent CFTR orthologs respond differently to the channel inhibitors CFTRinh-172, glibenclamide, and GlyH-101

Comparison of diverse orthologs is a powerful tool to study the structure and function of channel proteins. We investigated the response of human, killifish, pig, and shark cystic fibrosis transmembrane conductance regulator (CFTR) to specific inhibitors of the channel: CFTR(inh)-172, glibenclamide, and GlyH-101. In three systems, including organ perfusion of the shark rectal gland, primary cultures of shark rectal gland tubules, and expression studies of each ortholog in cRNA microinjected Xenopus laevis oocytes, we observed fundamental differences in the sensitivity to inhibition by these channel blockers. In organ perfusion studies, shark CFTR was insensitive to inhibition by CFTR(inh)-172. This insensitivity was also seen in short-circuit current experiments with cultured rectal gland tubular epithelial cells (maximum inhibition 4 ± 1.3%). In oocyte expression studies, shark CFTR was again insensitive to CFTR(inh)-172 (maximum inhibition 10.3 ± 2.5% at 25 μM), pig CFTR was insensitive to glibenclamide (maximum inhibition 18.4 ± 4.4% at 250 μM), and all orthologs were sensitive to GlyH-101. The amino acid residues considered responsible by previous site-directed mutagenesis for binding of the three inhibitors are conserved in the four CFTR isoforms studied. These experiments demonstrate a profound difference in the sensitivity of different orthologs of CFTR proteins to inhibition by CFTR blockers that cannot be explained by mutagenesis of single amino acids. We believe that the potency of the inhibitors CFTR(inh)-172, glibenclamide, and GlyH-101 on the CFTR chloride channel protein is likely dictated by the local environment and the three-dimensional structure of additional residues that form the vestibules, the chloride pore, and regulatory regions of the channel.     

Structural basis of typhoid: Salmonella typhi type IVb pilin (PilS) and cystic fibrosis transmembrane conductance regulator interaction

The type IVb pilus of the enteropathogenic bacteria Salmonella typhi is a major adhesion factor during the entry of this pathogen into gastrointestinal epithelial cells. Its target of adhesion is a stretch of 10 residues from the first extracellular domain of cystic fibrosis transmembrane conductance regulator (CFTR). The crystal structure of the N‐terminal 25 amino acid deleted S. typhi native PilS protein (ΔPilS), which makes the pilus, was determined at 1.9 Å resolution by the multiwavelength anomalous dispersion method. Also, the structure of the complex of ΔPilS and a target CFTR peptide, determined at 1.8 Å, confirms that residues 113–117 (NKEER) of CFTR are involved in binding with the pilin protein and gives us insight on the amino acids that are essential for binding. Furthermore, we have also explored the role of a conserved disulfide bridge in pilus formation. The subunit structure and assembly architecture are crucial for understanding pilus functions and designing suitable therapeutics against typhoid. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

乳腺癌易感蛋白BRCA1的BRCT1结构域染色质伸展活性的定位

乳腺癌易感基因BRCA1（Breast cancer susceptibility gene 1）在乳腺癌的发生、发展中起重要作用。BRCA1 C末端含有2个BRCT结构域（BRCT1和BRCT2）,许多乳腺癌易感突变发生在BRCA1的BRCT结构域中。利用染色质结构检测技术表明,BRCT结构域具有染色质伸展活性。本文利用缺失突变技术构建了6种BRCT1结构域（1642-1736 aa）缺失突变体并将BRCT1结构域中与染色质伸展相关的重要区域定位到1691-1721之间的氨基酸残基;用丙氨酸扫描技术构建了10种BRCT1结构域丙氨酸扫描突变体并将重要氨基酸残基序列定位到1707-1711之间的IAGGK。利用定位的重要区域进行Blast分析,结果找到一新型同源蛋白质。BRCT1结构域的定位有助于预测BRCT1结构域突变后发生乳腺癌的风险,也为进一步研究BRCT1结构域的功能机制提供了有用的材料。     

In-Silico Analyses of Nonsynonymous Variants in the BRCA1 Gene

BReast CAncer gene 1 (BRCA1)—a tumor suppressor gene plays an important role in the DNA repair mechanism. Several BRCA1 variants perturb its structure and function, including synonymous and nonsynonymous single nucleotide polymorphisms (SNPs). In the present study, we performed in-silico analyses of nonsynonymous SNPs (nsSNPs) of the BRCA1 gene. In total, 122 nsSNPs were retrieved from the NCBI SNP database and in-silico analyses were performed using computational prediction tools: SIFT, PROVEAN, Mutation Taster, PolyPhen-2, MutPred, and ConSurf. Of these tools, SIFT, PROVEAN, and Mutation Taster predicted 61 out of 122 nsSNPs as “damaging”, based on structural homology analysis. PolyPhen-2 classified 22 nsSNPs as “probably damaging”. These nsSNPs were further analyzed by MutPred to predict basic molecular mechanisms of amino acid alteration. ConSurf analysis predicted eleven conserved amino acid residues with structural and functional consequences. We identified five amino acid residues in the RING finger domain (L22, C39, H41, C44, and C47) and two in the BRCT domain (P1771 and I1707) with the potential to deter the BRCA1 protein function. This study provides insights into the effect of nsSNPs and amino acid substitutions in BRCA1.

     

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.     

Phosphorylation-induced conformational changes of cystic fibrosis transmembrane conductance regulator monitored by attenuated total reflection-Fourier transform IR spectroscopy and fluorescence spectroscopy

Cystic fibrosis transmembrane conductance regulator (CFTR) is a member of the ABC protein superfamily. Phosphorylation of a regulatory domain of this protein is a prerequisite for activity. We analyzed the effect of protein kinase A (PKA) phosphorylation on the structure of purified and reconstituted CFTR protein. 1H/2H exchange monitored by attenuated total reflection Fourier transform IR spectroscopy demonstrates that CFTR is highly accessible to aqueous medium. Phosphorylation of the regulatory (R) domain by PKA further increases this accessibility. More specifically, fluorescence quenching of cytosolic tryptophan residues revealed that the accessibility of the cytoplasmic part of the protein is modified by phosphorylation. Moreover, the combination of polarized IR spectroscopy with 1H/2H exchange suggested an increase of the accessibility of the transmembrane domains of CFTR. This suggests that CFTR phosphorylation can induce a large conformational change that could correspond either to a displacement of the R domain or to long range conformational changes transmitted from the phosphorylation sites to the nucleotide binding domains and the transmembrane segments. Such structural changes may provide better access for the solutes to the nucleotide binding domains and the ion binding site.     

Identification of deleterious SNPs and their effects on BCL11A,the master regulator of fetal hemoglobin expression

The B-cell lymphoma/leukemia 11A protein (encoded by BCL11A gene) is a key regulator of fetal-to-adult hemoglobin switching as seen in post-natal life. Although genetic polymorphisms like SNPs in BCL11A gene are expected to affect fetal hemoglobin (HbF) expression levels, yet their implications are poorly studied. This study utilizes a computational approach to identify the deleterious SNPs which may affect the structure and function of BCL11A protein. The study also generated a 3D structure of native and mutants. The analysis identified two SNPs in BCL11A as highly deleterious: N391K and C414S which are expected to affect structure and stability of the protein. According to conservation analysis, both residues N391 and C414 were identified as highly conserved. Additionally, post-translational modification sites were predicted at both sites. Ligand binding sites were also predicted in N391 and C414. Therefore, N391K and C414S in BCL11A can considered as important candidates to mediate HbF variation.     

Mutational analysis of the Saccharomyces cerevisiae ATP-binding cassette transporter protein Ycf1p

Ycf1p is a member of the ATP-binding cassette transporter family of membrane proteins. Strong sequence similarity has been observed between Ycf1p, the cystic fibrosis transmembrane conductance regulator (CFTR) and multidrug resistance protein (MRP). In this work, we have examined the functional significance of several of the conserved amino acid residues and the genetic requirements for Ycf1p subcellular localization. Biochemical fractionation experiments have established that Ycf1p, expressed at single-copy gene levels, co-fractionates with the vacuolar membrane and that this co-fractionation is independent of vps15 , vps34 or end3 gene function. Several cystic fibrosis-associated alleles of the CFTR were introduced into Ycf1p and found to elicit defects analogous to those seen in the CFTR. An amino-terminal extension shared between Ycf1p and MRP, but absent from CFTR, was found to be required for Ycf1p function, but not its subcellular localization. Mutant forms of Ycf1p were also identified that exhibited enhanced biological function relative to the wild-type protein. These studies indicate that Ycf1p will provide a simple, genetically tractable model system for the study of the trafficking and function of ATP-binding cassette transporter proteins, such as the CFTR and MRP.     

Construction of a synthetic gene for an R-plasmid-encoded dihydrofolate reductase and studies on the role of the N-terminus in the protein

R67 dihydrofolate reductase (DHFR) is a novel protein that provides clinical resistance to the antibacterial drug trimethoprim. The crystal structure of a dimeric form of R67 DHFR indicates the first 16 amino acids are disordered [Matthews et al. (1986) Biochemistry 25, 4194-4204]. To investigate whether these amino acids are necessary for protein function, the first 16 N-terminal residues have been cleaved off by chymotrypsin. The truncated protein is fully active with kcat = 1.3 s-1, Km(NADPH) = 3.0 microM, and Km(dihydrofolate) = 5.8 microM. This result suggests the functional core of the protein resides in the beta-barrel structure defined by residues 27-78. To study this protein further, synthetic genes coding for full-length and truncated R67 DHFRs were constructed. Surprisingly, the gene coding for truncated R67 DHFR does not produce protein in vivo or confer trimethoprim resistance upon Escherichia coli. Therefore, the relative stabilities of native and truncated R67 DHFR were investigated by equilibrium unfolding studies. Unfolding of dimeric native R67 DHFR is protein concentration dependent and can be described by a two-state model involving native dimer and unfolded monomer. Using absorbance, fluorescence, and circular dichroism techniques, an average delta GH2O of 13.9 kcal mol-1 is found for native R67 DHFR. In contrast, an average delta GH2O of 11.3 kcal mol-1 is observed for truncated R67 DHFR. These results indicate native R67 DHFR is 2.6 kcal mol-1 more stable than truncated protein. This stability difference may be part of the reason why protein from the truncated gene is not found in vivo in E. coli.     

Structure of the His44 --> Ala single point mutant of the distal finger motif of HIV-1 nucleocapsid protein: a combined NMR, molecular dynamics simulation, and fluorescence study

The nucleocapsid protein (NCp7) of human immunodeficiency virus type 1 (HIV-1) contains two highly conserved CCHC zinc fingers that strongly bind Zn(2+) through coordination of one His and three Cys residues. It has been suggested that NCp7 function is conformation specific since substitution of any of the zinc coordinating residues in the zinc finger motifs leads to subsequent loss of viral infectivity. To further determine the structural requirements necessary for this specific conformation, we investigated by (1)H 2D NMR and molecular dynamics simulations the structure of the distal finger motif of NCp7 in which the zinc coordinating amino acid, His 44, was substituted by a noncoordinating Ala residue. While the fold of the N-terminal part of this mutated peptide was similar to that of the native peptide, an increased lability and significant conformational changes were observed in the vicinity of the His-to-Ala mutation. Moreover, molecular dynamics simulations suggested a mechanism by which the variant peptide can bind zinc ion even though one zinc-coordinating amino acid was lacking. Using the fluorescence of the naturally occurring Trp37 residue, the binding affinity of the variant peptide to the (TG)(3) model oligonucleotide was found to be decreased by about 2 orders of magnitude with respect with the native peptide. Modeling of the DNA:NCp7 complex using structures of the variant peptide suggests that the residues forming a hydrophobic cleft in the native protein are improperly oriented for efficient DNA binding by the variant peptide.     

Development and characterization of synthetic antibodies binding to the cystic fibrosis conductance regulator

Cystic fibrosis transmembrane conductance regulator (CFTR) is a chloride channel in the apical surface of epithelial cells in the airway and gastrointestinal tract, and mutation of CFTR is the underlying cause of cystic fibrosis. However, the precise molecular details of the structure and function of CFTR in native and disease states remains elusive and cystic fibrosis researchers are hindered by a lack of high specificity, high affinity binding reagents for use in structural and biological studies. Here, we describe a panel of synthetic antigen-binding fragments (Fabs) isolated from a phage-displayed library that are specific for intracellular domains of CFTR that include the nucleotide-binding domains (NBD1 and NBD2), the R-region, and the regulatory insertion loop of NBD1. Binding assays performed under conditions that promote the native fold of the protein demonstrated that all Fabs recognized full-length CFTR. However, only the NBD1-specific Fab recognized denatured CFTR by western blot, suggesting a conformational epitope requirement for the other Fabs. Surface plasmon resonance experiments showed that the R-region Fab binds with high affinity to both the phosphorylated and unphosphorylated R-region. In addition, NMR analysis of bound versus unbound R-region revealed a distinct conformational effect upon Fab binding. We further defined residues involved with antibody recognition using an overlapping peptide array. In summary, we describe methodology complementary to previous hybridoma-based efforts to develop antibody reagents to CFTR, and introduce a synthetic antibody panel to aid structural and biological studies.     

R-domain interactions with distal regions of CFTR lead to phosphorylation and activation

Cystic fibrosis is caused by the aberrant function of the cystic fibrosis transmembrane conductance regulator (CFTR) protein. We examined whether intramolecular binding interactions involving the regulatory (R) domain contributed to CFTR regulation and function. When the R-domain (amino acids 596-836) was coexpressed with Delta1-836 CFTR (a carboxyl hemi-CFTR beginning immediately after the R-domain), strong binding between the two polypeptides was exhibited. The R-domain that co-immunoprecipitated with Delta1-836 exhibited a slower mobility on SDS-PAGE that resulted from phosphorylation of the protein. A larger CFTR polypeptide that included the R-domain (M837X) also exhibited a phosphorylation-dependent mobility shift when coexpressed with Delta1-836. Moreover, coexpression of M837X and Delta1-836 led to enhanced halide permeability in living cells. The activity, unlike in full-length CFTR, was present without forskolin activation, but still sensitive to the PKA inhibitor, Rp-8-CPT-cAMPS. This PKA inhibition of activity was found to be dependent on the carboxy region of the R-domain, amino acids 723-836. Our results indicate that the R-domain binds CFTR residues after amino acid 836 and that this binding facilitates phosphorylation and CFTR activation. We have also characterized a subdomain within CFTR (residues 723-837) that is necessary for PKA-dependent constitutive activation. Finally, these experiments demonstrate that constitutive CFTR activity can be accomplished by at least two mechanisms: (1) direct modulation of the R-domain to abrogate PKA regulation and (2) modifications that increase R-domain susceptibility to steady-state phosphorylation through PKA.