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     

Dopamine-derived quinones affect the structure of the redox sensor DJ-1 through modifications at Cys-106 and Cys-53

The physiological role of DJ-1, a protein involved in familial Parkinson disease is still controversial. One of the hypotheses proposed indicates a sensor role for oxidative stress, through oxidation of a conserved cysteine residue (Cys-106). The association of DJ-1 mutations with Parkinson disease suggests a loss of function, specific to dopaminergic neurons. Under oxidative conditions, highly reactive dopamine quinones (DAQs) can be produced, which can modify cysteine residues. In cellular models, DJ-1 was found covalently modified by dopamine. We analyzed the structural modifications induced on human DJ-1 by DAQs in vitro. We described the structural perturbations induced by DAQ adduct formation on each of the three cysteine residues of DJ-1 using specific mutants. Cys-53 is the most reactive residue and forms a covalent dimer also in SH-SY5Y DJ-1-transfected cells, but modification of Cys-106 induces the most severe structural perturbations; Cys-46 is not reactive. The relevance of these covalent modifications to the several functions ascribed to DJ-1 is discussed in the context of the cell response to a dopamine-derived oxidative insult.     

The C-terminal domain is the primary determinant of histone H1 binding to chromatin in vivo

We have used a combination of kinetic measurements and targeted mutations to show that the C-terminal domain is required for high-affinity binding of histone H1 to chromatin, and phosphorylations can disrupt binding by affecting the secondary structure of the C terminus. By measuring the fluorescence recovery after photo-bleaching profiles of green fluorescent protein-histone H1 proteins in living cells, we find that the deletion of the N terminus only modestly reduces binding affinity. Deletion of the C terminus, however, almost completely eliminates histone H1.1 binding. Specific mutations of the C-terminal domain identified Thr-152 and Ser-183 as novel regulatory switches that control the binding of histone H1.1 in vivo. It is remarkable that the single amino acid substitution of Thr-152 with glutamic acid was almost as effective as the truncation of the C terminus to amino acid 151 in destabilizing histone H1.1 binding in vivo. We found that modifications to the C terminus can affect histone H1 binding dramatically but have little or no influence on the charge distribution or the overall net charge of this domain. A comparison of individual point mutations and deletion mutants, when reviewed collectively, cannot be reconciled with simple charge-dependent mechanisms of C-terminal domain function of linker histones.     

Circular permutation and deletion studies of myoglobin indicate that the correct position of its N-terminus is required for native stability and solubility but not for native-like heme binding and folding

We studied the effect of deleted and circularly permuted mutations in sperm whale myoglobin and present here results on three classes of mutants: (i) a deletion mutant, Mb(1)(-)(99), in which the C-terminal helices, G and H, were removed; (ii) two circular permutations, Mb-B_GHA, in which helix B is N-terminal and helix A is C-terminal, and Mb-C_GHAB, in which helix C is N-terminal and helices A and B are C-terminal; and (iii) a deleted circular permutation, Mb-HAB_F, in which helix H is N-terminal, helix F is C-terminal, and helix G is deleted. The conformational characteristics of the apo and holo forms of these mutants were determined at neutral pH, by spectroscopic and hydrodynamic methods. The apo form of the deleted and permuted mutants exhibited a stronger tendency to aggregate and had lower ellipticity than the wild type. The mutants retained the ability to bind heme, but only the circularly permuted holoproteins had native-like heme binding and folding. These results agree with the theory that myoglobin has a central core that is able to bind heme, but also indicate that the presence of N- and C-terminal helices is necessary for native-like heme pocket formation. Because the holopermuteins were less stable than the wild-type protein and aggregated, we propose that the native position of the N-terminus is important for the precise structural architecture of myoglobin.     

Deletion mutagenesis of Tn10 Tet repressor — localization of regions important for dimerization and inducibility in vivo

The gene for the Tn 10 Tet repressor (TetR) was subjected to deletion mutagenesis. Screening for a transdominant operator-binding negative phenotype yielded 10 mutants with internal deletions. Three deletions extend from residue D5 to residues L41, W75, or Q76, respectively, and two contain deletions of the α-helix-turn-α-helix DNA-binding motif. Five deletions range from residue K84 to residues between R87 and K98. Since residues from the N-terminus up to position 98 are not necessary for dimerization, this must take place in the C-terminal half of the protein. Ability to dimerize was probed by introducing ochre non-sense codons (oc) at residues G138, H151, E159, l174, or K202. Koc202 shows wild-type in vivo operator-binding and inducibility by tetracycline indicating that the six C-terminal residues of TetR are not important for activity. Mutants with longer C-terminal truncations are inactive and not transdominant. They show reduced steady-state protein levels and are probably impaired in folding and degraded in vivo . Two mutants (Δ151–166, Δ164–166) with deletions in a region variable in primary structure and length among Tet repressers from different resistance determinants bind tet operator efficiently, but are not inducibie by tetracycline. This result indicates that these residues are not important for dimer formation in the operator-binding form.     

Oxidized DJ-1 interacts with the mitochondrial protein BCL-XL

Parkinson disease (PD)- and cancer-associated protein, DJ-1, mediates cellular protection via many signaling pathways. Deletions or mutations in the DJ-1 gene are directly linked to autosomal recessive early-onset PD. DJ-1 has potential roles in mitochondria. Here, we show that DJ-1 increases its mitochondrial distribution in response to ultraviolet B (UVB) irradiation and binds to Bcl-X(L). The interactions between DJ-1 and Bcl-X(L) are oxidation-dependent. DJ-1(C106A), a mutant form of DJ-1 that is unable to be oxidized, binds Bcl-X(L) much less than DJ-1 does. Moreover, DJ-1 stabilizes Bcl-X(L) protein level by inhibiting its ubiquitination and degradation through ubiquitin proteasome system (UPS) in response to UVB irradiation. Furthermore, under UVB irradiation, knockdown of DJ-1 leads to increases of Bcl-X(L) ubiquitination and degradation upon UVB irradiation, thereby increasing mitochondrial Bax, caspase-3 activation and PARP cleavage. These data suggest that DJ-1 protects cells against UVB-induced cell death dependent on its oxidation and its association with mitochondrial Bcl-X(L).     

Differential effects of Parkinson's disease-associated mutations on stability and folding of DJ-1

Mutations in the PARK7/DJ-1 gene cause autosomal-recessive Parkinson's disease. In some patients the gene is deleted. The molecular basis of disease in patients with point mutations is less obvious. We have investigated the molecular properties of [L166P]DJ-1 and the novel variant [E64D]DJ-1. When transfected into non-neuronal and neuronal cell lines, steady-state expression levels of [L166P]DJ-1 were dramatically lower than wild-type [WT]DJ-1 and [E64D]DJ-1. Cycloheximide and pulse-chase experiments revealed that the decreased expression levels of [L166P]DJ-1 were because of accelerated protein turnover. Proteasomal degradation was not the major pathway of DJ-1 breakdown because treatment with the proteasome inhibitor MG-132 caused only minimal accumulation of DJ-1, even of the very unstable [L166P]DJ-1 mutant. Because of the structural resemblance of DJ-1 with bacterial cysteine proteases, we considered an autoproteolytic mechanism. However, neither pharmacological inhibition nor site-directed mutagenesis of the putative active site residue Cys-106 stabilized DJ-1. To gain further insight into the structural defects of DJ-1 mutants, human [WT]DJ-1 and both mutants were expressed in Escherichia coli. As in eukaryotic cells, expression levels of [L166P]DJ-1 were dramatically reduced compared with [WT]DJ-1 and [E64D]DJ-1. Circular dichroism spectrometry revealed that the solution structures of [WT]DJ-1 and [E64D]DJ-1 are rich in beta-strand and alpha-helix conformation. Alpha-helices were more susceptible to thermal denaturation than the beta-sheet, and [WT]DJ-1 was more flexible in this regard than [E64D]DJ-1. Thus, structural defects of [E64D]DJ-1 only become apparent upon denaturing conditions, whereas the L166P mutation causes a drastic defect that leads to excessive degradation.     

Introducing Wilson disease mutations into the zinc-transporting P-type ATPase of Escherichia coli. The mutation P634L in the 'hinge' motif (GDGXNDXP) perturbs the formation of the E2P state.

ZntA, a bacterial zinc-transporting P-type ATPase, is homologous to two human ATPases mutated in Menkes and Wilson diseases. To explore the roles of the bacterial ATPase residues homologous to those involved in the human diseases, we have introduced several point mutations into ZntA. The mutants P401L, D628A and P634L correspond to the Wilson disease mutations P992L, D1267A and P1273L, respectively. The mutations D628A and P634L are located in the C-terminal part of the phosphorylation domain in the so-called hinge motif conserved in all P-type ATPases. P401L resides near the N-terminal portion of the phosphorylation domain whereas the mutations H475Q and P476L affect the heavy metal ATPase-specific HP motif in the nucleotide binding domain. All mutants show reduced ATPase activity corresponding 0-37% of the wild-type activity. The mutants P401L, H475Q and P476L are poorly phosphorylated by both ATP and P(i). Their dephosphorylation rates are slow. The D628A mutant is inactive and cannot be phosphorylated at all. In contrast, the mutant P634L six residues apart in the same domain shows normal phosphorylation by ATP. However, phosphorylation by P(i) is almost absent. In the absence of added ADP the P634L mutant dephosphorylates much more slowly than the wild-type, whereas in the presence of ADP the dephosphorylation rate is faster than that of the wild-type. We conclude that the mutation P634L affects the conversion between the states E1P and E2P so that the mutant favors the E1 or E1P state.     

Mutational analysis of the cytoplasmic tail of jaagsiekte sheep retrovirus envelope protein

Jaagsiekte sheep retrovirus (JSRV) is the etiologic agent of a transmissible lung cancer in sheep, ovine pulmonary adenocarcinoma. JSRV is unique in that the envelope protein functions as an oncogene, since it can morphologically transform fibroblast and epithelial cells in culture and can induce lung tumors in mice. Previous studies indicated that the transmembrane (TM) protein is essential for transformation, and particular attention has focused on a YXXM motif in the cytoplasmic tail. In this study, we carried out systematic mutagenesis of the cytoplasmic tail of JSRV Env. Alanine scanning mutagenesis revealed four classes of mutants: mutants in which transformation was abrogated, those in which transformation was not affected, those with reduced transformation, and those with increased transformation (supertransformers). In general, the alanine mutations did not affect Env protein production or its localization to the plasma membrane. Three functional domains of the cytoplasmic tail were identified: an amphipathic helix at the N-terminal (juxtamembrane) side, a nonessential C-terminal region, and an internal region (including the YXXM motif) where mutations resulted in abrogation, decreases, or increases in transformation. Alanine mutations in the amphipathic helix in both the hydrophobic and hydrophilic faces generally abolished transformation. The mutation R591A showed partial transformation that was consistent with loss of signaling through the Akt-mTOR pathway and signaling predominantly through the Ras-Raf-MEK1/2-extracellular signal-regulated kinase 1/2 pathway. The supertransforming mutants generally showed increased signaling through Akt and reduced activation of p38 MAPK that is inhibitory for transformation. These mutants provide further insight into the role of the TM cytoplasmic tail in JSRV transformation.     

A yeast mitochondrial membrane methyltransferase-like protein can compensate for oxa1 mutations

Members of the Oxa1p/Alb3/YidC family mediate the insertion of various organelle or bacterial hydrophobic proteins into membranes. They present at least five transmembrane segments (TM) linked by hydrophilic domains located on both sides of the membrane. To examine how Oxa1p structure relates to its function, we have introduced point mutations and large deletions into various domains of the yeast mitochondrial protein. These mutants allowed us to show the importance of the first TM domain as well as a synergistic interaction between the first loop and the C-terminal tail, which both protrude into the matrix. These mutants also led to the isolation of a high copy suppressor, OMS1, which encodes a member of the methyltransferase family. Overexpression of OMS1 seems to increase the steady-state level of both the mutant and wild-type Oxa1p. We show that Oms1p is a mitochondrial inner membrane protein inserted independently of Oxa1p. Oms1p presents one TM and a N-in C-out topology with the C-terminal domain carrying the methyltransferase-like domain. A conserved motif within this domain is essential for the suppression of oxa1 mutations. We discuss the possible role of Oms1p on Oxa1p intermembrane space domain.     

Aberrant regulation of methylesterase activity in cheD chemotaxis mutants of Escherichia coli

The adaptation process in several cheD chemotaxis mutants, which carry defects in tsr, the serine transducer gene, was examined. cheD mutants are smooth swimming and generally nonchemotactic; the defect is dominant to the wild-type tsr gene (J. S. Parkinson, J. Bacteriol. 142:953-961, 1980). All classes of methyl-accepting chemotaxis proteins synthesized in unstimulated cheD strains are overmethylated relative to the wild type. We found that the steady-state rate of demethylation in cheD mutants was low; this may explain their overmethylated phenotype. In addition, all cheD mutants showed diminished responsiveness of methylesterase activity to attractant and repellent stimuli transduced by either the Tsr or Tar protein, and they did not adapt. These results suggest that the dominant nature of the cheD mutations is manifested as a general defect in the regulation of demethylation. Some of these altered properties of methylesterase activity in cheD mutants were exhibited in wild-type cells that were treated with saturating concentrations of serine. The mutant Tsr protein thus seems to be locked into a signaling mode that suppresses tumbling and inhibits methylesterase activity in a global fashion. We found that the Tar and mutant Tsr proteins synthesized in cheD strains were methylated and deamidated at the correct sites and that the mutations were not located in the methylated peptides. Thus, the signaling properties of the transducers may be controlled at sites distinct from the methyl-accepting sites.     

rctB mutations that increase copy number of Vibrio cholerae oriCII in Escherichia coli

RctB serves as the initiator protein for replication from oriCII, the origin of replication of Vibrio cholerae chromosome II. RctB is conserved between members of Vibrionaceae but shows no homology to known replication initiator proteins and has no recognizable sequence motifs. We used an oriCII based minichromosome to isolate copy-up mutants in Escherichia coli. Three point mutations rctB(R269H), rctB(L439H) and rctB(Y381N) and one IS10 insertion in the 3'-end of the rctB gene were obtained. We determined the maximal C-terminal deletion that still gave rise to a functional RctB protein to be 165 amino acids. All rctB mutations led to decreased RctB-RctB interaction indicating that the monomer is the active form of the initiator protein. All mutations also showed various defects in rctB autoregulation. Loss of the C-terminal part of RctB led to overinitiation by reducing binding of RctB to both rctA and inc regions that normally serve to limit initiation from oriCII. Overproduction of RctB(R269H) and RctB(L439H) led to a rapid increase in oriCII copy number. This suggests that the initiator function of the two mutant proteins is increased relative to the wild-type.     

Effect of single amino acid substitution on oxidative modifications of the Parkinson's disease-related protein, DJ-1

Mutations in the gene encoding DJ-1 have been identified in patients with familial Parkinson's disease (PD) and are thought to inactivate a neuroprotective function. Oxidation of the sulfhydryl group to a sulfinic acid on cysteine residue C106 of DJ-1 yields the "2O " form, a variant of the protein with enhanced neuroprotective function. We hypothesized that some familial mutations disrupt DJ-1 activity by interfering with conversion of the protein to the 2O form. To address this hypothesis, we developed a novel quantitative mass spectrometry approach to measure relative changes in oxidation at specific sites in mutant DJ-1 as compared with the wild-type protein. Treatment of recombinant wild-type DJ-1 with a 10-fold molar excess of H(2)O(2) resulted in a robust oxidation of C106 to the sulfinic acid, whereas this modification was not detected in a sample of the familial PD mutant M26I exposed to identical conditions. Methionine oxidized isoforms of wild-type DJ-1 were depleted, presumably as a result of misfolding and aggregation, under conditions that normally promote conversion of the protein to the 2O form. These data suggest that the M26I familial substitution and methionine oxidation characteristic of sporadic PD may disrupt DJ-1 function by disfavoring a site-specific modification required for optimal neuroprotective activity. Our findings indicate that a single amino acid substitution can markedly alter a protein's ability to undergo oxidative modification, and they imply that stimulating the conversion of DJ-1 to the 2O form may be therapeutically beneficial in familial or sporadic PD.     

Aberrant folding of pathogenic Parkin mutants: aggregation versus degradation

Loss-of-function mutations in the Parkin gene (PARK2) are responsible for the majority of autosomal recessive Parkinson disease. A growing body of evidence indicates that misfolding and aggregation of Parkin is a major mechanism of Parkin inactivation, accounting for the loss-of-function phenotype of various pathogenic Parkin mutants. Remarkably, wild-type Parkin is also prone to misfolding under certain cellular conditions, suggesting a more general role of Parkin in the pathogenesis of Parkinson disease. We now show that misfolding of Parkin can lead to two phenotypes: the formation of detergent-insoluble, aggregated Parkin, or destabilization of Parkin resulting in an accelerated proteasomal degradation. By combining two pathogenic Parkin mutations, we could demonstrate that destabilization of Parkin is dominant over the formation of detergent-insoluble Parkin aggregates. Furthermore, a comparative analysis with HHARI, an E3 ubiquitin ligase with an RBR domain highly homologous to that of Parkin, revealed that folding of Parkin is specifically dependent on the integrity of the C-terminal domain, but not on the presence of a putative PDZ-binding motif at the extreme C terminus.     

Parkinson's disease-associated DJ-1 mutations impair mitochondrial dynamics and cause mitochondrial dysfunction

Mitochondrial dysfunction represents a critical event during the pathogenesis of Parkinson's disease (PD) and expanding evidences demonstrate that an altered balance in mitochondrial fission/fusion is likely an important mechanism leading to mitochondrial and neuronal dysfunction/degeneration. In this study, we investigated whether DJ-1 is involved in the regulation of mitochondrial dynamics and function in neuronal cells. Confocal and electron microscopic analysis demonstrated that M17 human neuroblastoma cells over-expressing wild-type DJ-1 (WT DJ-1 cells) displayed elongated mitochondria while M17 cells over-expressing PD-associated DJ-1 mutants (R98Q, D149A and L166P) (mutant DJ-1 cells) showed significant increase of fragmented mitochondria. Similar mitochondrial fragmentation was also noted in primary hippocampal neurons over-expressing PD-associated mutant forms of DJ-1. Functional analysis revealed that over-expression of PD-associated DJ-1 mutants resulted in mitochondria dysfunction and increased neuronal vulnerability to oxidative stress (H(2) O(2)) or neurotoxin. Further immunoblot studies demonstrated that levels of dynamin-like protein (DLP1), also known as Drp1, a regulator of mitochondrial fission, was significantly decreased in WT DJ-1 cells but increased in mutant DJ-1 cells. Importantly, DLP1 knockdown in these mutant DJ-1 cells rescued the abnormal mitochondria morphology and all associated mitochondria/neuronal dysfunction. Taken together, these studies suggest that DJ-1 is involved in the regulation of mitochondrial dynamics through modulation of DLP1 expression and PD-associated DJ-1 mutations may cause PD by impairing mitochondrial dynamics and function.     

Sequential mutations in the interleukin-3 (IL3)/granulocyte-macrophage colony-stimulating factor/IL5 receptor beta-subunit genes are necessary for the complete conversion to growth autonomy mediated by a truncated beta C subunit.

An amino-terminally truncated beta C receptor (beta C-R) subunit of the interleukin-3 (IL3)/granulocyte-macrophage colony-stimulating factor/IL5 receptor complex mediates factor-independent and tumorigenic growth in two spontaneous mutants of a promyelocytic cell line. The constitutive activation of the JAK2 protein kinase in these mutants confirms that signaling occurs through the truncated receptor protein. Noteworthily, in addition to a 10-kb deletion in the beta C-R subunit gene encoding the truncated receptor, several secondary and independent mutations that result in the deletion or functional inactivation of the allelic beta C-R subunit and the closely related beta IL3-R subunit genes were observed in both mutants, suggesting that such mutations are necessary for the full oncogenic penetrance of the truncated beta C-R subunit. Reversion of these mutations by the expression of the wild-type beta C-R in the two mutants resulted in a fivefold decrease in cloning efficiency of the mutants in the absence of IL3, confirming a functional interaction between the wild-type and truncated proteins. Furthermore, expression of the truncated beta C-R subunit in factor-dependent myeloid cells did not immediately render the cells autonomous but increased the spontaneous frequency to factor-independent growth by 4 orders of magnitude. Implications for both leukemogenic progression and receptor-subunit interaction and signaling are discussed.     

Contribution of the serine 129 of histone H2A to chromatin structure

Phosphorylation of a yeast histone H2A at C-terminal serine 129 has a central role in double-strand break repair. Mimicking H2A phosphorylation by replacement of serine 129 with glutamic acid (hta1-S129E) suggested that phosphorylation destabilizes chromatin structures and thereby facilitates the access of repair proteins. Here we have tested chromatin structures in hta1-S129 mutants and in a C-terminal tail deletion strain. We show that hta1-S129E affects neither nucleosome positioning in minichromosomes and genomic loci nor supercoiling of minichromosomes. Moreover, hta1-S129E has no effect on chromatin stability measured by conventional nuclease digestion, nor does it affect DNA accessibility and repair of UV-induced DNA lesions by nucleotide excision repair and photolyase in vivo. Similarly, deletion of the C-terminal tail has no effect on nucleosome positioning and stability. These data argue against a general role for the C-terminal tail in chromatin organization and suggest that phosphorylated H2A, gamma-H2AX in higher eukaryotes, acts by recruitment of repair components rather than by destabilizing chromatin structures.