Disease Mutations in the Human Mitochondrial DNA Polymerase Thumb Subdomain Impart Severe Defects in Mitochondrial DNA Replication

Forty-five different point mutations in POLG, the gene encoding the catalytic subunit of the human mitochondrial DNA polymerase (pol γ), cause the early onset mitochondrial DNA depletion disorder, Alpers syndrome. Sequence analysis of the C-terminal polymerase region of pol γ revealed a cluster of four Alpers mutations at highly conserved residues in the thumb subdomain (G848S, c.2542g→a; T851A, c.2551a→g; R852C, c.2554c→t; R853Q, c.2558g→a) and two Alpers mutations at less conserved positions in the adjacent palm subdomain (Q879H, c.2637g→t and T885S, c.2653a→t). Biochemical characterization of purified, recombinant forms of pol γ revealed that Alpers mutations in the thumb subdomain reduced polymerase activity more than 99% relative to the wild-type enzyme, whereas the palm subdomain mutations retained 50–70% wild-type polymerase activity. All six mutant enzymes retained physical and functional interaction with the pol γ accessory subunit (p55), and none of the six mutants exhibited defects in misinsertion fidelity in vitro. However, differential DNA binding by these mutants suggests a possible orientation of the DNA with respect to the polymerase during catalysis. To our knowledge this study represents the first structure-function analysis of the thumb subdomain in pol γ and examines the consequences of mitochondrial disease mutations in this region.As the only DNA polymerase found in animal cell mitochondria, DNA polymerase γ (pol γ)³ bears sole responsibility for DNA synthesis in all replication and repair transactions involving mitochondrial DNA (1, 2). Mammalian cell pol γ is a heterotrimeric complex composed of one catalytic subunit of 140 kDa (p140) and two 55-kDa accessory subunits (p55) that form a dimer (3). The catalytic subunit contains an N-terminal exonuclease domain connected by a linker region to a C-terminal polymerase domain. Whereas the exonuclease domain contains essential motifs I, II, and III for its activity, the polymerase domain comprising the thumb, palm, and finger subdomains contains motifs A, B, and C that are crucial for polymerase activity. The catalytic subunit is a family A DNA polymerase that includes bacterial pol I and T7 DNA polymerase and possesses DNA polymerase, 3′ → 5′ exonuclease, and 5′-deoxyribose phosphate lyase activities (for review, see Refs. 1 and 2). The 55-kDa accessory subunit (p55) confers processive DNA synthesis and tight binding of the pol γ complex to DNA (4, 5).Depletion of mtDNA as well as the accumulation of deletions and point mutations in mtDNA have been observed in several mitochondrial disorders (for review, see Ref. 6). mtDNA depletion syndromes are caused by defects in nuclear genes responsible for replication and maintenance of the mitochondrial genome (7). Mutation of POLG, the gene encoding the catalytic subunit of pol γ, is frequently involved in disorders linked to mutagenesis of mtDNA (8, 9). Presently, more than 150 point mutations in POLG are linked with a wide variety of mitochondrial diseases, including the autosomal dominant (ad) and recessive forms of progressive external ophthalmoplegia (PEO), Alpers syndrome, parkinsonism, ataxia-neuropathy syndromes, and male infertility (tools.niehs.nih.gov/polg) (9).Alpers syndrome, a hepatocerebral mtDNA depletion disorder, and myocerebrohepatopathy are rare heritable autosomal recessive diseases primarily affecting young children (10–12). These diseases generally manifest during the first few weeks to years of life, and symptoms gradually develop in a stepwise manner eventually leading to death. Alpers syndrome is characterized by refractory seizures, psychomotor regression, and hepatic failure (11, 12). Mutation of POLG was first linked to Alpers syndrome in 2004 (13), and to date 45 different point mutations in POLG (18 localized to the polymerase domain) are associated with Alpers syndrome (9, 14, 15). However, only two Alpers mutations (A467T and W748S, both in the linker region) have been biochemically characterized (16, 17).During the initial cloning and sequencing of the human, Drosophila, and chicken pol γ genes, we noted a highly conserved region N-terminal to motif A in the polymerase domain that was specific to pol γ (18). This region corresponds to part of the thumb subdomain that tracks DNA into the active site of both Escherichia coli pol I and T7 DNA polymerase (19–21). A high concentration of disease mutations, many associated with Alpers syndrome, is found in the thumb subdomain.Here we investigated six mitochondrial disease mutations clustered in the N-terminal portion of the polymerase domain of the enzyme (Fig. 1A). Four mutations (G848S, c.2542g→a; T851A, c.2551a→g; R852C, c.2554c→t; R853Q, c.2558g→a) reside in the thumb subdomain and two (Q879H, c.2637g→t and T885S, c.2653a→t) are located in the palm subdomain. These mutations are associated with Alpers, PEO, mitochondrial myopathy, encephalopathy, lactic acidosis, and stroke-like episodes (MELAS), ataxia-neuropathy syndrome, Leigh syndrome, and myocerebrohepatopathy (

Nickel-based Enzyme Systems

Of the eight known nickel enzymes, all but glyoxylase I catalyze the use and/or production of gases central to the global carbon, nitrogen, and oxygen cycles. Nickel appears to have been selected for its plasticity in coordination and redox chemistry and is able to cycle through three redox states (1+, 2+, 3+) and to catalyze reactions spanning ∼1.5 V. This minireview focuses on the catalytic mechanisms of nickel enzymes, with an emphasis on the role(s) of the metal center. The metal centers vary from mononuclear to complex metal clusters and catalyze simple hydrolytic to multistep redox reactions.Seven of the eight known nickel enzymes (1). CODH² interconverts CO and CO₂; ACS utilizes CO; the nickel ARD produces CO; hydrogenase generates/utilizes hydrogen gas; MCR generates methane; urease produces ammonia; and SOD generates O₂.

TABLE 1

Nickel-containing enzymes

Chimeric Nitrogenase-like Enzymes of (Bacterio)chlorophyll Biosynthesis

Nitrogenase-like light-independent protochlorophyllide oxidoreductase (DPOR) is involved in chlorophyll biosynthesis. Bacteriochlorophyll formation additionally requires the structurally related chlorophyllide oxidoreductase (COR). During catalysis, homodimeric subunit BchL₂ or ChlL₂ of DPOR transfers electrons to the corresponding heterotetrameric catalytic subunit, (BchNB)₂ or (ChlNB)₂. Analogously, subunit BchX₂ of the COR enzymes delivers electrons to subunit (BchYZ)₂. Various chimeric DPOR enzymes formed between recombinant subunits (BchNB)₂ and BchL₂ from Chlorobaculum tepidum or (ChlNB)₂ and ChlL₂ from Prochlorococcus marinus and Thermosynechococcus elongatus were found to be enzymatically active, indicating a conserved docking surface for the interaction of both DPOR protein subunits. Biotin label transfer experiments revealed the interaction of P. marinus ChlL₂ with both subunits, ChlN and ChlB, of the (ChlNB)₂ tetramer. Based on these findings and on structural information from the homologous nitrogenase system, a site-directed mutagenesis approach yielded 10 DPOR mutants for the characterization of amino acid residues involved in protein-protein interaction. Surface-exposed residues Tyr¹²⁷ of subunit ChlL, Leu⁷⁰ and Val¹⁰⁷ of subunit ChlN, and Gly⁶⁶ of subunit ChlB were found essential for P. marinus DPOR activity. Next, the BchL₂ or ChlL₂ part of DPOR was exchanged with electron-transferring BchX₂ subunits of COR and NifH₂ of nitrogenase. Active chimeric DPOR was generated via a combination of BchX₂ from C. tepidum or Roseobacter denitrificans with (BchNB)₂ from C. tepidum. No DPOR activity was observed for the chimeric enzyme consisting of NifH₂ from Azotobacter vinelandii in combination with (BchNB)₂ from C. tepidum or (ChlNB)₂ from P. marinus and T. elongatus, respectively.Chlorophyll and bacteriochlorophyll biosynthesis, as well as nitrogen fixation, are essential biochemical processes developed early in the evolution of life (1). During biological fixation of nitrogen, nitrogenase catalyzes the reduction of atmospheric dinitrogen to ammonia (2). Enzyme systems homologous to nitrogenase play a crucial role in the formation of the chlorin and bacteriochlorin ring system of chlorophylls (Chl)2 and bacteriochlorophylls (Bchl) (3, 4) (Fig. 1a). For the synthesis of both Chl and Bchl, the stereospecific reduction of the C-17-C-18 double bond of ring D of protochlorophyllide (Pchlide) catalyzed by the nitrogenase-like enzyme light-independent (dark-operative) protochlorophyllide oxidoreductase (DPOR) results in the formation of chlorophyllide (Chlide) (Fig. 1a, left) (5, 6). DPOR enzymes consist of three protein subunits which are designated BchN, BchB and BchL in Bchl-synthesizing organisms and ChlN, ChlB and ChlL in Chl-synthesizing organisms. A second reduction step at ring B (C-7-C-8) unique to the synthesis of Bchl converts the chlorin Chlide into a bacteriochlorin ring structure to form bacteriochlorophyllide (Bchlide) (Fig. 1a, right, Bchlide). This reaction is catalyzed by another nitrogenase-like enzyme, termed chlorophyllide oxidoreductase (COR) (7). COR enzymes are composed of subunits BchY, BchZ, and BchX.Open in a separate windowFIGURE 1.Comparison of the three subunit enzymes DPOR, COR, and nitrogenase. a, during Chl and Bchl biosynthesis, ring D is stereospecifically reduced by the nitrogenase-like enzyme DPOR (subunit composition BchL₂/(BchNB)₂ or ChlL₂/(ChlNB)₂) leading to the chlorin Chlide. Subunits N, B, and L are named ChlN, ChlB, and ChlL in Chl-synthesizing organisms and BchN, BchB, and BchL in Bchl-synthesizing organisms. The synthesis of Bchl additionally requires the stereospecific B ring reduction by a second nitrogenase-like enzyme called COR, with the subunit composition BchX₂/(BchYZ)₂. COR catalyzes the formation of the bacteriochlorin Bchlide. Subunits Y, Z, and X of the COR enzyme are named BchY, BchZ, and BchX. b, the homologous nitrogenase complex has the subunit composition NifH₂/(NifD/NifK)₂. Rings A–E and the carbon atoms are designated according to IUPAC nomenclature (41). R is either a vinyl or an ethyl moiety. The position marked by an asterisk indicates either a vinyl or a hydroxyethyl moiety (42).All subunits share significant amino acid sequence homology to the corresponding subunits of nitrogenase, which are designated NifD, NifK, and NifH, respectively (1) (compare Fig. 1, a and b). Whereas subunits BchL or ChlL, BchX and NifH exhibit a sequence identity at the amino acid level of ∼33%, subunits BchN or ChlN, BchY, NifD, and BchB or ChlB, BchZ, and NifK, respectively, show lower sequence identities of ∼15% (1). For all enzymes a common oligomeric protein architecture has been proposed consisting of the heterotetrameric complexes (BchNB)₂ or (ChlNB)₂, (BchYZ)₂, and (NifD/NifK)₂, which are completed by a homodimeric protein subunit BchL₂ or ChlL₂, BchX₂, and NifH₂, respectively (compare Fig. 1, a and b) (3, 7, 8).Nitrogenase is a well characterized protein complex that catalyzes the reduction of nitrogen to ammonia in a reaction that requires at least 16 molecules of MgATP (2, 9, 10). During nitrogenase catalysis, subunit NifH₂ (Fe protein) associates with and dissociates from the (NifD/NifK)₂ complex (MoFe protein). Binding, hydrolysis of MgATP and structural rearrangements are coupled to sequential intersubunit electron transfer. For this purpose, NifH₂ contains an ATP-binding motif and an intersubunit [4Fe-4S] cluster coordinated by two cysteine residues from each NifH monomer (1, 11). Electrons from this [4Fe-4S] cluster are transferred via a [8Fe-7S] cluster (P-cluster) onto the [1Mo-7Fe-9S-X-homocitrate] cluster (MoFe cofactor). Both of the latter clusters are located on (NifD/NifK)₂, where dinitrogen is reduced to ammonia (10). Three-dimensional structures of NifH₂ in complex with (NifD/NifK)₂ revealed a detailed picture of the dynamic interaction of both subcomplexes (8, 12).Based on biochemical and bioinformatic approaches, it has been proposed that the initial steps of DPOR reaction strongly resemble nitrogenase catalysis. Key amino acid residues essential for DPOR function have been identified by mutagenesis of the enzyme from Chlorobaculum tepidum (formerly denoted as Chlorobium tepidum) (3). The catalytic mechanism of DPOR includes the electron transfer from a “plant-type” [2Fe-2S] ferredoxin onto the dimeric DPOR subunit, BchL₂, carrying an intersubunit [4Fe-4S] redox center coordinated by Cys⁹⁷ and Cys¹³¹ in C. tepidum. Analogous to nitrogenase, Lys¹⁰ in the phosphate-binding loop (P-loop) and Leu¹²⁶ in the switch II region of DPOR were found essential for DPOR catalysis. Moreover, it was shown that the BchL₂ protein from C. tepidum does not form a stable complex with the catalytic (BchNB)₂ subcomplex. Therefore, a transient interaction responsible for the electron transfer onto protein subunit (BchNB)₂ has been proposed (3).The subsequent [Fe-S] cluster-dependent catalysis and the specific substrate recognition at the active site located on subunit (BchNB)₂ are unrelated to nitrogenase. The (BchNB)₂ subcomplex was shown to carry a second [4Fe-4S] cluster, which was proposed to be ligated by Cys²¹, Cys⁴⁶, and Cys¹⁰³ of the BchN subunit and Cys⁹⁴ of subunit BchB (C. tepidum numbering) (3). No evidence for any type of additional cofactor was obtained from biochemical and EPR spectroscopic analyses (5, 13). Thus, despite the same common oligomeric architecture, the catalytic subunits (BchNB)₂ and (ChlNB)₂ clearly differ from the corresponding nitrogenase complex, as no molybdenum-containing cofactor or P-cluster equivalent is employed (5, 14). From these results it was concluded that electrons from the [4Fe-4S] cluster of (BchNB)₂ or (ChlNB)₂ are transferred directly onto the Pchlide substrate at the active site of DPOR.The second nitrogenase-like enzyme, COR, catalyzes the reduction of ring B of Chlide during the biosynthesis of Bchl (7). Therefore, an accurate discrimination of the ring systems of the individual substrates is required. COR subunits share an overall amino acid sequence identity of 15–22% for BchY and BchZ and 31–35% for subunit BchX when compared with the corresponding DPOR subunits (supplemental Figures S2–S4). In amino acid sequence alignments of BchX proteins with the closely related BchL or ChlL subunits of DPOR, both cysteinyl ligands responsible for [4Fe-4S] cluster formation and residues for ATP binding are conserved (1). Furthermore, all cysteinyl residues characterized as ligands for a catalytic [4Fe-4S] cluster in (BchNB)₂ or (ChlNB)₂ are conserved in the sequences of subunits BchY and BchZ of COR (7). These findings correspond to a recent EPR study in which a characteristic signal for a [4Fe-4S] cluster was obtained for the COR subunit BchX₂ as well as for subunit (BchYZ)₂ (15). These results indicate that the catalytic mechanism of COR strongly resembles DPOR catalysis. In vitro assays for nitrogenase as well as for DPOR and COR make use of the artificial electron donor dithionite in the presence of high concentrations of ATP (7, 16, 17).

TABLE 1

Amino acid sequence identities of the individual subunits of DPOR, COR, and nitrogenaseAmino acid sequences of the individual subunits of DPOR, COR, and nitrogenase employed in the present study (compareFig. 3A) were aligned by using the ClustalW method in MegAlign (DNASTAR), and sequence identities were calculated.

Correlation of Fragile Histidine Triad (Fhit) Protein Structural Features with Effector Interactions and Biological Functions

We have previously shown that Fhit tumor suppressor protein interacts with Hsp60 chaperone machinery and ferredoxin reductase (Fdxr) protein. Fhit-effector interactions are associated with a Fhit-dependent increase in Fdxr stability, followed by generation of reactive oxygen species and apoptosis induction under conditions of oxidative stress. To define Fhit structural features that affect interactions, downstream signaling, and biological outcomes, we used cancer cells expressing Fhit mutants with amino acid substitutions that alter enzymatic activity, enzyme substrate binding, or phosphorylation at tyrosine 114. Gastric cancer cell clones stably expressing mutants that do not bind substrate or cannot be phosphorylated showed decreased binding to Hsp60 and Fdxr and reduced mitochondrial localization. Expression of Fhit or mutants that bind interactor proteins results in oxidative damage and accumulation of cells in G₂/M or sub-G₁ fractions after peroxide treatment; noninteracting mutants are defective in these biological effects. Gastric cancer clones expressing noncomplexing Fhit mutants show reduction of Fhit tumor suppressor activity, confirming that substrate binding, interaction with heat shock proteins, mitochondrial localization, and interaction with Fdxr are important for Fhit tumor suppressor function.Fhit protein is a powerful tumor suppressor that is frequently lost or reduced in cancer cells because of rearrangement of the exquisitely DNA damage-sensitive fragile FHIT gene. Restoration of Fhit expression suppresses tumorigenicity of cancer cells of various types, and the ability to induce apoptosis in cancer cells in vitro is reduced by specific Fhit mutations (1, 2).Through studies of signal pathways affected by Fhit expression, by searches for Fhit protein effectors, and by in vitro analyses of Fhit activity, we and others have defined Fhit enzymatic activity in vitro (3), apoptotic activity in cells and tumors (4–6), and most recently identification of a Fhit protein complex that affects Fhit stability, mitochondrial localization, and interaction with ferredoxin reductase (Fdxr)⁵ (7). The complex includes Hsp60 and Hsp10 that mediate Fhit stability and may affect import into mitochondria, where Fhit interacts with Fdxr, which is responsible for transferring electrons from NADPH to cytochrome P450 via ferredoxin. Virally mediated Fhit restoration in Fhit-deficient cancer cells increases production of intracellular reactive oxygen species (ROS), followed by increased apoptosis of cancer cells under oxidative stress conditions; conversely, Fhit-negative cells escape apoptosis, likely carrying oxidative DNA damage that contributes to accumulation of mutations.The Fhit protein sequence, showing high homology to the histidine triad (HIT) family of proteins, suggested that the protein product would hydrolyze diadenosine tetraphosphate or diadenosine triphosphate (Ap₃A) (8), and in vitro studies showed that Ap₃A was cleaved into ADP and AMP by Fhit. The catalytic histidine triad within Fhit was essential for catalytic activity (3), and a Fhit mutant that substituted Asn for His at the central histidine (H96N mutant) was catalytically inactive, although it bound substrate well (3). Early tumor suppression studies showed that cancer cells stably transfected with wild type (WT) or H96N mutant Fhit were suppressed for tumor growth in nude mice. This suggested the hypothesis that the Fhit-substrate complex sends the tumor suppression signal (9, 10). To test this hypothesis, a series of FHIT alleles was designed to reduce substrate-binding and/or hydrolytic rates and was characterized by quantitative cell-death assays on cancer cells virally infected with each allele. The allele series covered defects as great as 100,000-fold in k_cat and increases as large as 30-fold in K_m. Mutants with 2–7-fold increases in K_m had significantly reduced apoptotic indices and the mutant with a 30-fold increase in K_m retained little apoptotic function. Thus, the proapoptotic function of Fhit, which is likely associated with tumor suppressor function, is limited by substrate binding and is unrelated to substrate hydrolysis (11).Fhit, a homodimeric protein of 147 amino acids, is a target of tyrosine phosphorylation by the Src family protein kinases, which can phosphorylate Tyr-114 of Fhit in vitro and in vivo (12). After co-expression of Fhit with the Elk tyrosine kinase in Escherichia coli to generate phosphorylated forms of Fhit, unphosphorylated, mono-, and diphosphorylated Fhit were purified, and enzyme kinetics studies showed that monophosphorylated Fhit exhibited monophasic kinetics with K_m and k_cat values ∼2- and ∼7-fold lower, respectively, than for unphosphorylated Fhit. Diphosphorylated Fhit exhibited biphasic kinetics; one site had K_m and k_cat values ∼2- and ∼140-fold lower, respectively, than for unphosphorylated Fhit; the second site had a K_m ∼60-fold higher and a k_cat ∼6-fold lower than for unphosphorylated Fhit (13). Thus, it was possible that the alterations in K_m and k_cat values for phosphorylated forms of Fhit might favor formation and lifetime of the Fhit-Ap₃A complex and enhance tumor suppressor activity (see Fhit forms

---

Kinetic parameters

---

% Sub-G₁

---

Direct binding

---

Subcellular location

---

Co-IP in vivo

---

8-OHdG

---

Apoptosis

---

Tumor suppressor

---

K_m (mm)k_cat (s^–1)A549MKN74Hsp60FdxrHsp60Fdxr Fhit WT 1.6 +/– 0.19 2.7 +/– 0.95 43 24 Yes Yes Cyt & mito Yes Yes Yes Yes Yes Catalyt mutants    H96D Up 2-fold Down >2 × 10⁴ 29 NT NT NT Cyt & mito Yes Yes NT Yes NT    H96N Up 2-fold Down >5 × 10⁵ 31 14.4 NT NT Cyt & mito Yes Yes Yes Yes Yes Loop mutants    Y114A Up 23-fold Down 2-fold 3.7 NT NT NT Cyt +/– +/– +/– No No    Y114D NT NT 2.9 6 NT NT Cyt +/– +/– – No –/+    Y114E NT NT NT NT NT NT Cyt & mito –/+ –/+ – No NT    Y114F Up 5-fold Up 1.1-fold 11.5 3 NT NT Cyt & mito –/+ –/+ – No No    Y114W Up 5-fold Up 1.4-fold NT NT NT NT Cyt & mito –/+ – – NT NT    del113–117 Up 10-fold Down 38-fold 5 NT NT NT NT NT NT – No NT Other mutants    L25W Up 7-fold Down 4-fold 15 NT NT NT Cyt – – – NT –/+    I10W,L25W Up 32-fold Down 6-fold 11 NT NT NT NT NT NT NT NT NT    F5W Up 3.3 fold NT NT 5 NT NT NT NT NT +/– No NT Purified pFhit    pFhit Down 0.4-fold Down 7-fold NA NA –/+ Yes NA NA NA NA NA NA    ppFhit Down 0.4-fold Down > 100-fold NA NA –/+ Yes NA NA NA NA NA NA Up 60-fold Down 6-fold Open in a separate windowTo explore the in vivo importance of the Tyr-114 phosphorylation site and define Fhit-mediated signaling events, Semba et al. (14) compared the differential biological effects of Ad-FHIT WT and Ad-FHIT Tyr-114 mutant expression in human lung cancer cells. Caspase-dependent apoptosis was effectively induced only by WT Fhit protein. However, the biological significance of phosphorylation at Tyr-114 has been difficult to study because the endogenous phosphorylated forms have very short half-lives; activation of epidermal growth facto receptor family members induces Fhit phosphorylation by Src and proteasome degradation of phosphorylated Fhit (15).Although there are possible connections among the various pathways known to be altered in Fhit-deficient cells, apoptosis, DNA damage-response checkpoint activation, ROS production, and related biological effects of Fhit loss or overexpression, details of the pathway(s) leading from Fhit overexpression to cell death and tumor suppression have not been delineated. Now that a Fhit signaling complex has been identified, we set out to examine which structural features of Fhit protein might participate in individual steps of the pathway leading from Fhit overexpression through complex formation, subcellular localization, interaction with mitochondrial Fdxr, DNA damage induction, cell cycle changes, apoptosis, and ultimately tumor suppression. The underlying hypotheses were as follows: substrate-binding mutants would behave similarly to WT; nonsubstrate-binding mutants would be defective in some step of the pathway, perhaps complexing with heat shock proteins or Fdxr or perhaps induction of DNA damage; and Tyr-114 mutants, which also affect formation or stability of the enzyme-substrate complex, would also be defective in executing some step of the Fhit overexpression pathway to cell death. One goal was to identify specific mutants that exhibited deficiency in specific steps of the pathway, so that such mutants could be used to dissect each step in more detail. Using in vitro Fhit and Fhit-effector protein interactions, we aimed to determine the following: 1) which proteins of the complex interact directly with Fhit, and 2) the biological role of these interactions in vivo. Using cancer cells expressing exogenous WT and mutant Fhit proteins, we were able to examine the structural features of Fhit that affect the direct interaction with its effectors, participate in ROS production, and are necessary for tumor suppression activity.     

Adaptive Divergence in Experimental Populations of Pseudomonas fluorescens. IV. Genetic Constraints Guide Evolutionary Trajectories in a Parallel Adaptive Radiation

The capacity for phenotypic evolution is dependent upon complex webs of functional interactions that connect genotype and phenotype. Wrinkly spreader (WS) genotypes arise repeatedly during the course of a model Pseudomonas adaptive radiation. Previous work showed that the evolution of WS variation was explained in part by spontaneous mutations in wspF, a component of the Wsp-signaling module, but also drew attention to the existence of unknown mutational causes. Here, we identify two new mutational pathways (Aws and Mws) that allow realization of the WS phenotype: in common with the Wsp module these pathways contain a di-guanylate cyclase-encoding gene subject to negative regulation. Together, mutations in the Wsp, Aws, and Mws regulatory modules account for the spectrum of WS phenotype-generating mutations found among a collection of 26 spontaneously arising WS genotypes obtained from independent adaptive radiations. Despite a large number of potential mutational pathways, the repeated discovery of mutations in a small number of loci (parallel evolution) prompted the construction of an ancestral genotype devoid of known (Wsp, Aws, and Mws) regulatory modules to see whether the types derived from this genotype could converge upon the WS phenotype via a novel route. Such types—with equivalent fitness effects—did emerge, although they took significantly longer to do so. Together our data provide an explanation for why WS evolution follows a limited number of mutational pathways and show how genetic architecture can bias the molecular variation presented to selection.UNDERSTANDING—and importantly, predicting—phenotypic evolution requires knowledge of the factors that affect the translation of mutation into phenotypic variation—the raw material of adaptive evolution. While much is known about mutation rate (e.g., Drake et al. 1998; Hudson et al. 2002), knowledge of the processes affecting the translation of DNA sequence variation into phenotypic variation is minimal.Advances in knowledge on at least two fronts suggest that progress in understanding the rules governing the generation of phenotypic variation is possible (Stern and Orgogozo 2009). The first stems from increased awareness of the genetic architecture underlying specific adaptive phenotypes and recognition of the fact that the capacity for evolutionary change is likely to be constrained by this architecture (Schlichting and Murren 2004; Hansen 2006). The second is the growing number of reports of parallel evolution (e.g., Pigeon et al. 1997; ffrench-Constant et al. 1998; Allender et al. 2003; Colosimo et al. 2004; Zhong et al. 2004; Boughman et al. 2005; Shindo et al. 2005; Kronforst et al. 2006; Woods et al. 2006; Zhang 2006; Bantinaki et al. 2007; McGregor et al. 2007; Ostrowski et al. 2008)—that is, the independent evolution of similar or identical features in two or more lineages—which suggests the possibility that evolution may follow a limited number of pathways (Schluter 1996). Indeed, giving substance to this idea are studies that show that mutations underlying parallel phenotypic evolution are nonrandomly distributed and typically clustered in homologous genes (Stern and Orgogozo 2008).While the nonrandom distribution of mutations during parallel genetic evolution may reflect constraints due to genetic architecture, some have argued that the primary cause is strong selection (e.g., Wichman et al. 1999; Woods et al. 2006). A means of disentangling the roles of population processes (selection) from genetic architecture is necessary for progress (Maynard Smith et al. 1985; Brakefield 2006); also necessary is insight into precisely how genetic architecture might bias the production of mutations presented to selection.Despite their relative simplicity, microbial populations offer opportunities to advance knowledge. The wrinkly spreader (WS) morphotype is one of many different niche specialist genotypes that emerge when experimental populations of Pseudomonas fluorescens are propagated in spatially structured microcosms (Rainey and Travisano 1998). Previous studies defined, via gene inactivation, the essential phenotypic and genetic traits that define a single WS genotype known as LSWS (Spiers et al. 2002, 2003) (Figure 1). LSWS differs from the ancestral SM genotype by a single nonsynonymous nucleotide change in wspF. Functionally (see Figure 2), WspF is a methyl esterase and negative regulator of the WspR di-guanylate cyclase (DGC) (Goymer et al. 2006) that is responsible for the biosynthesis of c-di-GMP (Malone et al. 2007), the allosteric activator of cellulose synthesis enzymes (Ross et al. 1987). The net effect of the wspF mutation is to promote physiological changes that lead to the formation of a microbial mat at the air–liquid interface of static broth microcosms (Rainey and Rainey 2003).Open in a separate windowFigure 1.—Outline of experimental strategy for elucidation of WS-generating mutations and their subsequent identity and distribution among a collection of independently evolved, spontaneously arising WS genotypes. The strategy involves, first, the genetic analysis of a specific WS genotype (e.g., LSWS) to identify the causal mutation, and second, a survey of DNA sequence variation at specific loci known to harbor causal mutations among a collection of spontaneously arising WS genotypes. For example, suppressor analysis of LSWS using a transposon to inactivate genes necessary for expression of the wrinkly morphology delivered a large number of candidate genes (top left) (Spiers et al. 2002). Genetic and functional analysis of these candidate genes (e.g., Goymer et al. 2006) led eventually to the identity of the spontaneous mutation (in wspF) responsible for the evolution of LSWS from the ancestral SM genotype (Bantinaki et al. 2007). Subsequent analysis of the wspF sequence among 26 independent WS genotypes (bottom) showed that 50% harbored spontaneous mutations (of different kinds; see Open in a separate windowFigure 2.—Network diagram of DGC-encoding pathways underpinning the evolution of the WS phenotype and their regulation. Overproduction of c-di-GMP results in overproduction of cellulose and other adhesive factors that determine the WS phenotype. The ancestral SBW25 genome contains 39 putative DGCs, each in principle capable of synthesizing the production of c-di-GMP, and yet WS genotypes arise most commonly as a consequence of mutations in just three DGC-containing pathways: Wsp, Aws, and Mws. In each instance, the causal mutations are most commonly in the negative regulatory component: wspF, awsX, and the phosphodiesterase domain of mwsR (see text).To determine whether spontaneous mutations in wspF are a common cause of the WS phenotype, the nucleotide sequence of this gene was obtained from a collection of 26 spontaneously arising WS genotypes (WS_A-Z) taken from 26 independent adaptive radiations, each founded by the same ancestral SM genotype (Figure 1): 13 contained mutations in wspF (Bantinaki et al. 2007). The existence of additional mutational pathways to WS provided the initial motivation for this study.

TABLE 1

Mutational causes of WS

Mode of Action of cGMP-dependent Protein Kinase-specific Inhibitors Probed by Photoaffinity Cross-linking Mass Spectrometry

The inhibitor peptide DT-2 (YGRKKRRQRRRPPLRKKKKKH) is the most potent and selective inhibitor of the cGMP-dependent protein kinase (PKG) known today. DT-2 is a construct of a PKG tight binding sequence (W45, LRKKKKKH, K_I = 0.8 μm) and a membrane translocating sequence (DT-6, YGRKKRRQRRRPP, K_I = 1.1 μm), that combined strongly inhibits PKG catalyzed phosphorylation (K_I = 12.5 nm) with ∼1000-fold selectivity toward PKG over protein kinase A, the closest relative of PKG. However, the molecular mechanism behind this inhibition is not entirely understood. Using a combination of photoaffinity labeling, stable isotope labeling, and mass spectrometry, we have located the binding sites of PKG-specific substrate and inhibitor peptides. Covalent linkage of a PKG-specific substrate analogue was localized in the catalytic core on residues 356–372, also known as the glycine-rich loop, essential for ATP binding. By analogy, the individual inhibitor peptides W45 and DT-6 were also found to cross-link near the glycine-rich loop, suggesting these are both substrate competitive inhibitors. A bifunctional photoreactive analogue of DT-2 was found to generate dimers of PKG. This cross-linking induced covalent PKG dimerization was not observed for an N-terminal deletion mutant of PKG, which lacks the dimerization domain. In addition, non-covalent mass spectrometry was used to determine binding stoichiometry and binding order of the inhibitor peptides. Dimeric PKG binds two W45 and DT-6 peptides, whereas only one DT-2 molecule was observed to bind to the dimeric PKG. Taken together, these findings imply that (i) the two individual components making up DT-2 are both targeted against the substrate-binding site and (ii) binding of a single DT-2 molecule inactivates both PKG monomers simultaneously, which is an indication that (iii) in cGMP-activated PKG the catalytic centers of both subunits may be in each other''s proximity.Among the superfamily of protein kinases the two cyclic nucleotide-regulated protein kinases, cAMP-dependent protein kinase and cGMP-dependent protein kinase, form a closely related subfamily of serine/threonine protein kinases (1–4). Both proteins share several structural elements, such as the N-terminal dimerization domain, an autoinhibition site, two in-tandem cyclic nucleotide-binding sites, and a highly conserved catalytic core (Fig. 1, A and B). Despite these similarities, these two enzymes display differences, which account for their unique properties. Whereas PKA² is nearly ubiquitous, PKG is primarily found in the lung, cerebellum, and smooth muscles (5, 6). From a structural point of view these cyclic nucleotide-dependent protein kinases differ as well. The holoenzyme of PKA is a tetramer composed of two regulatory and two catalytic subunits. The catalytic subunits are non-covalently attached to the regulatory subunit dimer. Upon interaction with cAMP, the catalytic subunits dissociate from the holoenzyme and are free to catalyze heterophosphorylation (Fig. 1C). The mammalian type I PKGs are homodimeric cytosolic proteins containing two identical polypeptides of ∼76 kDa. Alternative mRNA splicing produces type Iα and type Iβ PKG, which are identical proteins apart from their first ∼100 N-terminal residues (7). Each PKG subunit is composed of a regulatory and a catalytic domain on a single polypeptide chain. Consequently, when cGMP activates PKG, the catalytic and regulatory components remain physically attached (Fig. 1D). Within the catalytic domain PKA and PKG share a strong primary sequence homology (8). Not surprisingly, these enzymes also exhibit overlapping substrate specificities, a feature that often interferes with efforts to elucidate their distinct biological pathways. Peptide substrates with a primary amino acid sequence motif RRX(S/T)X are in general recognized by both PKA and PKG (9). Besides this strong overlapping substrate specificity, several studies report on subtle differences in determinants that discriminate for PKA and PKG substrate specificity (10–16). To specifically discriminate between PKG and PKA activity in biological assays a highly specific PKG peptide inhibitor was developed (17). This peptide, YGRKKRRQRRRPPLRKKKKKH (DT-2), is the most potent and selective PKG inhibitor known today. Recently, the validity of DT-2 as a superior inhibitor of PKG in terms of potency, selectivity, and membrane permeability has been demonstrated (18–24). The inhibitor is a construct of a substrate competitive sequence, LRKKKKKH (W45), derived from a library screen that selected for tight PKG binding sequences, with a significant specificity toward PKG over PKA, and a membrane translocating signal peptide, YGRKKRRQRRRPP (DT-6). DT-2 strongly inhibits PKG-catalyzed phosphorylation (K_i = 12.5 nm), however, the molecular nature of DT-2 inhibition is not entirely understood (25). Because high resolution structural data are not available for PKG, one of our goals is to elucidate binding sites for PKG-specific substrates and inhibitors in more detail using a combination of mass spectrometric techniques and photoaffinity labeling. To further delineate the nature of inhibition we have developed photoaffinity analogues of DT-2 and related inhibitory peptides, as well as a high affinity peptide substrate. The method of photoaffinity labeling enables the direct probing of target proteins through a covalent bond, which is photochemically introduced between a ligand and its specific receptor (26). In combination with modern mass spectrometric techniques this is a powerful approach for the characterization of peptide-protein interactions (27). Substrate and inhibitor peptides containing photoactivatable analogues of phenylalanine, 4-benzoyl-l-phenylalanine (Phe(Bz)) or 4′-(3-(trifluoromethyl)-3H-diazirin-3-yl)-l-phenylalanine (Phe(Tmd)) were synthesized and used to locate their substrate/inhibitor-binding sites on PKG. These measurements indicate that the substrate peptide resides near the glycine-rich loop within the catalytic domain and that the inhibitor peptides are directed similarly toward this substrate-binding site, thereby acting as competitive inhibitors. In addition, nanoflow electrospray ionization time of flight mass spectrometry (ESI-TOF-MS) was performed to study the interaction between DT-2 and PKG in more detail. ESI-MS has proven to be a useful tool to analyze the non-covalent interaction of proteins with ligands, oligonucleotides, peptides, or other proteins (28–31). Using this technique, important information on conformational changes (32–35), measurement of relative dissociation constants (36, 37), and sequential binding order and cooperativity (38, 39) can be obtained. ESI-MS confirms that PKG is primarily a homodimer and is able to bind four cGMP molecules. Binding of DT-2 was strongly enhanced in the presence of cGMP. Surprising is the observation that only one DT-2 molecule binds to dimeric PKG. The information derived from these measurements allows for molecular modeling and structural refinements of the next generation of PKG-selective inhibitors.Open in a separate windowFIGURE 1.Linear arrangement of the functional domains of the regulatory and catalytic subunit of PKA (A) and PKG (B) type I and schematic representation of the current working models of the activation process of PKA (C) and PKG (D) type 1. Binding of cAMP to the PKA induces a conformational change that results in the dissociation of the catalytic subunits. Binding of cGMP to PKG also induces a conformational change, which exposes the catalytic domains, but both catalytic domains remain near each other via the N-terminal dimerization domain. (Images adapted from Scholten et al. (4).)

TABLE 1

Inhibition contants (K_I) of PKA- or PKG-specific peptide inhibitors and the PKA/PKG specificity index