Abstract
Single nucleotide polymorphisms (SNPs) are the most common source of human genetic variation, occurring at a frequency of 1 in every 1000 nucleotide bases of the 3 billion bases human genome. SNPs can be intronic or exonic and can thus affect mRNA transcription, structure, splicing and translation, protein expression, structure, enzymatic activity and stability. SNPs are linked large number of genetic defects and also associated with many diseases. In this study we analyzed structural and dynamical effects of some of the important SNPs from Arylamine N-acetyltransferase 2 (NAT2) and the transcription factor pituitary homeobox protein 2 (PITX2). NAT2 is an important catalytic enzyme that metabolizes the carcinogenic arylamines, hydrazine drugs and chemicals. This enzyme is highly polymorphic in different human populations. Several polymorphisms of NAT2, including the single amino acid substitutions R64Q, I114T, D122N, L137F, Q145P, R197Q, and G286E, are classified as slow acetylators, whereas the wild-type NAT2 is classified as a fast acetylator. The slow acetylators are often associated with drug toxicity and efficacy as well as cancer susceptibility. The biological functions of these 7 mutations have previously been characterized, but the structural basis behind the reduced catalytic activity and reduced protein level is not clear. We performed multiple molecular dynamics simulations of these mutants as well as NAT2 to investigate the structural and dynamical effects throughout the protein structure, specifically the catalytic triad, cofactor binding site, and the substrate binding pocket. None of these mutations induced unfolding; instead, their effects were confined to the inter-domain, domain 3 and 17-residue insert region, where the flexibility was significantly reduced relative to the wild-type. Structural effects of these mutations propagate through space and cause a change in catalytic triad conformation, cofactor binding site, substrate binding pocket size/shape and electrostatic potential. Our results showed that the dynamical properties of all the mutant structures, especially in inter-domain, domain 3 and 17-residue insert region were affected in the same manner. Similarly, the electrostatic potential of all the mutants were altered and also the functionally important regions such as catalytic triad, cofactor binding site, and substrate binding pocket adopted different orientation and/or conformation relative to the wild-type that may affect the functions of the mutants. Overall, our study may provide the structural basis for reduced catalytic activity and protein level, as was experimentally observed for these polymorphisms. PITX2 is involved in genetic control of development. Mutations in PITX2, most in the homeodomain, cause the autosomal-dominant disorder Rieger syndrome. The mutants L16Q, K50E and R53P destabilize the structure and disrupt DNA-binding activity. The biological functions of these mutants have been characterized but not the structural basis behind the loss of DNA-binding activity. We performed multiple molecular dynamics simulations at 37º C to investigate the structural and dynamic effects of these mutants. Compared with the wild type (WT), the L16Q mutant induces a kink in the alpha 3 helix, which is stabilized by the hydrogen bond of Q21-R59. The disruption in backbone hydrogen bonds of V47-N51 and W48-R52 leads to a kink formation in the α3 helix of K50E. The R53P mutant alters the relative orientation of helices, which is apparently stabilized by the formation of new hydrogen bonds of T38-Q11, T38-Q12, T38-R2, N39-R2, L40-Q1, L40-R2, and T41-Q4. The hydrophobic core residues F8, L13, L40 and V45 change their positions in all mutants to break the hydrophobic core. Thus, changes in helical orientations and hydrophobic core cause rearrangement of the DNA-binding surface and disrupt DNA-binding activity in the mutants. The structural and molecular dynamics properties of 3 PITX2 homeodomain mutants differ from those of the WT, especially in formation of a kink in the recognition helix, change in the packing of helices and disruption of the hydrophobic core. Overall, our study may provide the structural basis for reduced catalytic activity and protein level, as was experimentally observed for these polymorphisms in NAT2 and the structural basis for the loss of DNA-binding activity for these polymorphisms in PITX2 homeodomain.