Sequence analysis of NAT2 gene in Brazilians: Identification of undescribed single nucleotide polymorphisms and molecular modeling of the N-acetyltransferase 2 protein structure

N-Acetyltransferase 2 (NAT2) metabolizes a variety of xenobiotics that includes many drugs, chemicals and carcinogens. This enzyme is genetically variable in human populations and polymorphisms in the NAT2 gene have been associated with drug toxicity and efficacy as well as cancer susceptibility. Here, we have focused on the identification of NAT2 variants in Brazilian individuals from two different regions, Rio de Janeiro and Goiás, by direct sequencing, and on the characterization of new haplotypes after cloning and re-sequencing. Upon analysis of DNA samples from 404 individuals, six new SNPs (c.29T>C, c.152G>T, c.203G>A, c.228C>T, c.458C>T and c.600A>G) and seven new NAT2 alleles were identified with different frequencies in Rio de Janeiro and Goiás. All new SNPs were found as singletons (observed only once in 808 genes) and were confirmed by three independent technical replicates. Molecular modeling and structural analysis suggested that p.Gly51Val variant may have an important effect on substrate recognition by NAT2. We also observed that amino acid change p.Cys68Tyr would affect acetylating activity due to the resulting geometric restrictions and incompatibility of the functional group in the Tyr side chain with the admitted chemical mechanism for catalysis by NATs. Moreover, other variants, such like p.Thr153Ile, p.Thr193Met, p.Pro228Leu and p.Val280Met, may lead to the presence of hydrophobic residues on NAT2 surface involved in protein aggregation and/or targeted degradation. Finally, the new alleles NAT2*6H and NAT2*5N, which showed the highest frequency in the Brazilian populations considered in this study, may code for a slow activity. Functional studies are needed to clarify the mechanisms by which new SNPs interfere with acetylation.     

Functional genomics of C190T single nucleotide polymorphism in human N-acetyltransferase 2

N-acetyltransferase 2 (NAT2) catalyzes N-acetylation and O-acetylation of many drugs and environmental carcinogens. Genetic polymorphisms in the NAT2 gene have been associated with differential susceptibility to cancers and drug toxicity from these compounds. Single nucleotide polymorphisms (SNPs) have been identified in the human NAT2 coding region. A new allele, NAT2*19, possessing the C190T (R64W) exchange, was recently identified. In order to understand the effect of this new SNP, recombinant NAT2*4 (reference) and NAT2*19 were expressed in yeast (Schizosaccharomyces pombe). The C190T (R64W) SNP in NAT2*19 caused substantial reduction in the NAT2 protein level and stability, but did not cause significant reduction in transformation efficiency or mRNA level. The enzymatic activities for N-acetylation of two arylamine carcinogens (2-aminofluorene, 4-aminobiphenyl), and a sulfonamide drug (sulfamethazine) were over 100-fold lower for NAT2 19 compared to reference NAT2 4. Kinetic studies showed a reduction in Vmax but no significant change in substrate Km. In addition, the SNP caused significant reduction in the O-acetylation of the N-hydroxy-2-amino-1-methyl-6-phenylimidazo [4,5-b] pyridine. These results show that NAT2*19 possessing the C190T (R64W) SNP encodes a slow acetylator phenotype for both N- and O-acetylation, due to a reduction in the amount and stability of the NAT2 19 allozyme.     

The human BARX2 gene: genomic structure, chromosomal localization, and single nucleotide polymorphisms

The BARX genes 1 and 2 are Bar class homeobox genes expressed in craniofacial structures during development. In this report, we present the genomic structure, chromosomal localization, and polymorphic markers in BARX2. The gene has four exons, ranging in size from 85 to 1099 bp. BARX2 is localized on human chromosome 11q25, as determined by radiation hybrid mapping. In the mouse, Barx2 is coexpressed with Pitx2 in several tissues. Based on the coexpression, BARX2 was assumed to be a candidate gene for those cases of Rieger syndrome that cannot be associated with mutations of PITX2. Mutations in PITX2 cause some cases of Rieger syndrome, an autosomal dominant disorder affecting eyes, teeth, and umbilicus. DNA from Rieger patients was subjected to single-strand conformation polymorphism screening of the BARX2 coding region. Three single nucleotide polymorphisms were found in a normal population, although no etiologic mutations were detectable in over 100 cases of Rieger syndrome or in individuals with related ocular disorders.     

人类基因组SNPs的研究现状及应用前景

基因组DNA是生物体各种生理、病理性状的物质基础,人类DNA序列变异约90%表现为单核苷酸多态性(singlenucleotidepolymorphisms,SNPs),这是一种常见的遗传变异类型,在人类基因组中广泛存在,被认为是人类疾病易感性和药物反应的决定性因素。本文主要介绍了SNPs的分类及特点、人类基因组SNPs的研究现状、SNPs在实践中的应用,以及SNPs在遗传作图、医药、遗传易感性、个体化医疗等方面的研究前景,并探讨了当前SNPs研究中存在的问题。     

An in silico analysis of deleterious single nucleotide polymorphisms and molecular dynamics simulation of disease linked mutations in genes responsible for neurodegenerative disorder

Abstract

Mutation in two genes deglycase gene (DJ-1) and retromer complex component gene (VPS35) are linked with neurodegenerative disorder such as Parkinson's disease, Huntington's disease, and Alzheimer's disease. DJ-1 gene located at 1p36 chromosomal position and involved in PD pathogenesis through many pathways including mitochondrial dysfunction and oxidative injury. VPS35 gene located at 16q13-q21 chromosomal position and the two pathways, the Wnt signaling pathway, and retromer-mediated DMT1 missorting are proposed for basis of VPS35 related PD. The study focuses on identifying most deleterious SNPs through computational analysis. Result obtained from various bioinformatics tools shows that D149A is most deleterious in DJ-1 and A54W, R365H, and V717M are most deleterious in VPS35. To understand the functionality of protein comparative modeling of DJ-1 and VPS35 native and mutants was done by MODELLER. The generated structures are validated by two web servers–ProSa and RAMPAGE. Molecular dynamic simulation (MDS) analysis done for the most validated structures to know the functional and structural nature of native and mutants protein of DJ-1 and VPS35. Native structure of DJ-1 and VPS35 show more flexibility through MDS analysis. DJ-1 D149A mutant structures become more compact which shows the structural perturbation and loss of DJ-1 protein function which in turn are probable cause for PD. A54W, R365H, and V717M mutant protein of VPS35 also shows compactness which cause structure perturbation and absence of retromer function which likely to be linked to PD pathogenesis. This in silico study may provide a new insight for fundamental molecular mechanism involved in Parkinson’s disease.

Communicated by Ramaswamy H. Sarma     

Cost-effective procedures for genotyping of human FCN2 gene single nucleotide polymorphisms

L-ficolin (ficolin-2) is a complement-activating pattern-recognition lectin taking part in the innate immune response. Both its serum concentration and sugar binding capacity are influenced by single nucleotide polymorphisms (SNP) of the corresponding FCN2 gene. Cost-effective and simple procedures, based on polymerase chain reaction (PCR) or PCR-restriction fragment length polymorphism for an investigation of four FCN2 SNPs are proposed: ?64 A?>?C (rs7865453), ?4 A?>?G (rs17514136; both located in the promoter region), +6359 C?>?T (rs17549193), +6424 G?>?T (rs7851696; both in exon 8). Variant alleles of ?64 and +6424 (in strong linkage disequlibrium) are known to be associated with low L-ficolin level or activity. In contrast, variant alleles at positions ?4 and +6359 (also in strong linkage disequlibrium) correspond to higher values. Since several L-ficolin clinical associations have been reported, FCN2 genotyping seems to be a valuable tool for disease association studies.     

Fluorescence detection of single nucleotide polymorphisms using a universal molecular beacon

We present a simple and novel assay—employing a universal molecular beacon (MB) in the presence of Hg²⁺—for the detection of single nucleotide polymorphisms (SNPs) based on Hg²⁺–DNA complexes inducing a conformational change in the MB. The MB (T₇-MB) contains a 19-mer loop and a stem of a pair of seven thymidine (T) bases, a carboxyfluorescein (FAM) unit at the 5′-end, and a 4-([4-(dimethylamino)phenyl]azo)benzoic acid (DABCYL) unit at the 3′-end. Upon formation of Hg²⁺–T₇-MB complexes through T–Hg²⁺–T bonding, the conformation of T₇-MB changes from a random coil to a folded structure, leading to a decreased distance between the FAM and DABCYL units and, hence, increased efficiency of fluorescence resonance energy transfer (FRET) between the FAM and DABCYL units, resulting in decreased fluorescence intensity of the MB. In the presence of complementary DNA, double-stranded DNA complexes form (instead of the Hg²⁺–T₇-MB complexes), with FRET between the FAM and DABCYL units occurring to a lesser extent than in the folded structure. Under the optimal conditions (20 nM T₇-MB, 20 mM NaCl, 1.0 μM Hg²⁺, 5.0 mM phosphate buffer solution, pH 7.4), the linear plot of the fluorescence intensity against the concentration of perfectly matched DNA was linear over the range 2–30 nM (R² = 0.991), with a limit of detection of 0.5 nM at a signal-to-noise ratio of 3. This new probe provides higher selectivity toward DNA than that exhibited by conventional MBs.     

Mutational spectrum in the recent human genome inferred by single nucleotide polymorphisms

So far, there is no genome-wide estimation of the mutational spectrum in humans. In this study, we systematically examined the directionality of the point mutations and maintenance of GC content in the human genome using approximately 1.8 million high-quality human single nucleotide polymorphisms and their ancestral sequences in chimpanzees. The frequency of C-->T (G-->A) changes was the highest among all mutation types and the frequency of each type of transition was approximately fourfold that of each type of transversion. In intergenic regions, when the GC content increased, the frequency of changes from G or C increased. In exons, the frequency of G:C-->A:T was the highest among the genomic categories and contributed mainly by the frequent mutations at the CpG sites. In contrast, mutations at the CpG sites, or CpG-->TpG/CpA mutations, occurred less frequently in the CpG islands relative to intergenic regions with similar GC content. Our results suggest that the GC content is overall not in equilibrium in the human genome, with a trend toward shifting the human genome to be AT rich and shifting the GC content of a region to approach the genome average. Our results, which differ from previous estimates based on limited loci or on the rodent lineage, provide the first representative and reliable mutational spectrum in the recent human genome and categorized genomic regions.     

Discovery of single nucleotide polymorphisms and mutations by pyrosequencing

Comparative genomics, analyzing variation among individual genomes, is an area of intense investigation. DNA sequencing is usually employed to look for polymorphisms and mutations. Pyrosequencing, a real-time DNA sequencing method, is emerging as a popular platform for comparative genomics. Here we review the use of this technology for mutation scanning, polymorphism discovery and chemical haplotyping. We describe the methodology and accuracy of this technique and discuss how to reduce the cost for large-scale analysis.     

Detection of single nucleotide polymorphisms

Single nucleotide polymorphism (SNP) detection technologies are used to scan for new polymorphisms and to determine the allele(s) of a known polymorphism in target sequences. SNP detection technologies have evolved from labor intensive, time consuming, and expensive processes to some of the most highly automated, efficient, and relatively inexpensive methods. Driven by the Human Genome Project, these technologies are now maturing and robust strategies are found in both SNP discovery and genotyping areas. The nearly completed human genome sequence provides the reference against which all other sequencing data can be compared. Global SNP discovery is therefore only limited by the amount of funding available for the activity. Local, target, SNP discovery relies mostly on direct DNA sequencing or on denaturing high performance liquid chromatography (dHPLC). The number of SNP genotyping methods has exploded in recent years and many robust methods are currently available. The demand for SNP genotyping is great, however, and no one method is able to meet the needs of all studies using SNPs. Despite the considerable gains over the last decade, new approaches must be developed to lower the cost and increase the speed of SNP detection.     

Structural analysis of human CCR2b and primate CCR2b by molecular modeling and molecular dynamics simulation

CCR2b, a chemokine receptor for MCP-1, -2, -3, -4, plays an important role in a variety of diseases involving infection, inflammation, and/or injury, as well as being a coreceptor for HIV-1 infection. Two models of human CCR2b (hCCR2b) were generated by homology modeling and 1 ns restrained molecular dynamics (MD) simulation. In one only C113-C190 forms a disulfide bond (SS model); in another the potential C32-C277 disulfide bond was formed (2SS model). Analysis of the structures and averaged displacements of Calpha atoms of the N-terminal residues shows that the main differences between the SS and 2SS models lie in a region D25YDYGAPCHKFD36; in the extracellular part of the 2SS model the accessible surfaces of N12, F23, Y26, Y28 and F35 are obviously raised and a more stable H-bond net is formed. The potential energy of the 2SS-water assembly finally fluctuated around -43,020 kJ x mol(-1), which is about 302 kJ x mol(-1) lower than that of the SS-water assembly. All these results suggest that the 2SS model is more favorable. The CCR2b genes of 17 primates were sequenced and four CCR2b models for primates Ateles paniscus (A. pan), Hylobates leucogyneus(H. leu), Papio cynocephalus (P. cyn) and Trachypithecus francoist ( T. fra) were generated based on the 2SS model. A comparison of hCCR2b with primate CCR2b also supports the importance of the region D25YDYGAPCHKFD36. Electrostatic potential maps of human and primate CCR2b all display the dipolar characteristics of CCR2b with the negative pole located in the extracellular part and a strong positive pole in the cytoplasmic part. Based on the CCR2b model, we suggest that the main functional residues fall in the D25YDYGAPCHKFD36 region, and the negative electrostatic feature is a non-specific, but necessary, factor for ligands or gp120/CD4 binding.     

Membrane structure of the human immunodeficiency virus gp41 fusion domain by molecular dynamics simulation

The structures of the 16-residue fusion domain (or fusion peptide, FP) of the human immunodeficiency virus gp41 fusion protein, two of its mutants, and a shortened peptide (5-16) were studied by molecular dynamics simulation in an explicit palmitoyloleoylphosphoethanolamine bilayer. The simulations showed that the active wild-type FP inserts into the bilayer approximately 44 degrees +/- 6 degrees with respect to the bilayer normal, whereas the inactive V2E and L9R mutants and the inactive 5 to 16 fragment lie on the bilayer surface. This is the first demonstration by explicit molecular dynamics of the oblique insertion of the fusion domain into lipid bilayers, and provides correlation between the mode of insertion and the fusogenic activity of these peptides. The membrane structure of the wild-type FP is remarkably similar to that of the influenza HA(2) FP as determined by nuclear magnetic resonance and electron spin resistance power saturation. The secondary structures of the wild-type FP and the two inactive mutants are quite similar, indicating that the secondary structure of this fusion domain plays little or no role in affecting the fusogenic activity of the fusion peptide. The insertion of the wild-type FP increases the thickness of the interfacial area of the bilayer by disrupting the hydrocarbon chains and extending the interfacial area toward the head group region, an effect that was not observed in the inactive FPs.     

Assessment of population structure by single nucleotide polymorphisms (SNPs) in goat breeds

Single nucleotide polymorphisms (SNPs) may be used in biodiversity studies and commercial tasks like traceability, paternity testing and selection for suitable genotypes. Twenty-seven SNPs were characterized and genotyped on 250 individuals belonging to eight Italian goat breeds. Multilocus genotype data were used to infer population structure and assign individuals to populations. To estimate the number of groups (K) to test in population structure analysis we used likelihood values and variance of the bootstrap samples, deriving optimal K from a drop in the likelihood and a rise in the variance plots against K.     

Identification and characterization of functional single nucleotide polymorphisms (SNPs) in Axin 1 gene: a molecular dynamics approach

Wnt signaling pathway has been reported to play crucial role in intestinal crypt formation and deregulation of this pathway is responsible for colorectal cancer initiation and progression. Axin 1, a scaffold protein, play pivotal role in the regulation of Wnt/β-catenin signaling pathway and has been found to be mutated in several cancers; primarily in colon cancer. Considering its crucial role, a structural and functional analysis of missense mutations in Axin 1 gene was performed in this study. Initially, one hundred non-synonymous single nucleotide polymorphisms in the coding regions of Axin 1 gene were selected for in silico analysis. Six variants (G820S, G856S, E830K, L811V, L847V, and R767C) were predicted to be deleterious by combinatorial prediction. Further investigation of structural attributes confirmed two highly deleterious single nucleotide polymorphisms (G820S and G856S). Molecular dynamics simulation demonstrated variation in different structural attributes between native and two highly deleterious Axin 1 mutant models. Finally, docking analysis showed variation in binding affinity of mutant Axin 1 proteins with two destruction complex members, GSK3β and adenomatous polyposis. The results collectively showed the deleterious effect of the above predicted single nucleotide polymorphisms on the Axin 1 protein structure and could prove to be an adjunct in the disease genotype-phenotype correlation studies.     

单核苷酸多态概述

单核苷酸多态ＳＮＰ是遍布于基因组中的一种ＤＮＡ序列变化类型,人类基因组中平均约每一千碱基中有一个。单核苷酸多态是一种双等位型多态,群体中出现的频率大于１％或２％者视为多态,低于１％或２％的则视为突变。由于其具有高信息量、高密度又便于自动化操作的特点,单核苷酸多态在遗传性疾病基因的克隆和药物的设计与开发方面具有广阔的应用前景。本文对单核苷酸的概念、特点、应用前景,及其研究应用的一些问题作一综述。