* Biology Department, Franklin and Marshall College; and Hopkins Marine Station of Stanford University, Pacific Grove, California
Correspondence: E-mail: peter.fields{at}fandm.edu.
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Key Words: activation energy cold adaptation notothenioids protein stability site-directed mutagenesis substrate affinity
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Over evolutionary time it is likely that functional attributes of the enzyme molecule will respond to selective pressure, in order to compensate for alterations in environmental temperature (Hochachka and Somero 2002). For example, a series of studies (Johnston and Whalesby 1977; Yancey and Somero 1978; Graves and Somero 1982; Yancey and Seibenaller 1987; Dahlhoff and Somero 1993; Holland, McFall-Ngai, and Somero 1997; Fields and Somero 1997, 1998; Johns and Somero 2004) has shown that kinetic properties of muscle-type lactate dehydrogenase (A4-LDH; E.C. 1.1.1.27; lactate : NAD+ oxidoreductase) from ectothermic vertebrates are strongly correlated to environmental temperature. When the apparent Michaelis constants for the substrate pyruvate of A4-LDH orthologs from a variety of species are measured at one temperature, values increase as habitat temperature decreases (Fields 2001; Hochachka and Somero 2002). Because
is inversely related to substrate binding affinity, these results indicate that, when measured at the same temperature, A4-LDH orthologs from cold-adapted species have lower substrate binding affinities than those from warm-adapted ones. However, when measured within each species' habitat temperature range,
values are comparable, indicating that the conservation of substrate binding affinity is important in temperature adaptation (Somero 1995). Similarly, when the catalytic rate constant, or turnover number (kcat), is measured at a common temperature, orthologs from more cold-adapted species again have higher values; that is, orthologs from more cold-adapted species catalyze the conversion of pyruvate to lactate faster than those from more warm-adapted species (Fields 2001; Hochachka and Somero 2002). These differences observed in catalytic rate appear adaptive, because they are found repeatedly in orthologs of species occurring at different temperatures. Modifications of kcat likely are due to changes in activation enthalpy, or, equivalently, in activation energy barriers to catalysis (Feller and Gerday 1997; Lonhienne, Gerday, and Feller 2000; Fields et al. 2002; Feller 2003), such that cold-adapted forms have reduced activation energies and can proceed through catalysis with a lower input of thermal energy from the medium.
The studies described above provide strong evidence of compensation to environmental temperature in A4-LDH orthologs of different ectothermic vertebrates, and this offers an opportunity to study the link between functional adaptation and molecular evolution. Fortunately, A4-LDH is an ideal system in which to study the structural basis of enzymic adaptation because its catalytic mechanism and three-dimensional structure have been described in detail; in addition, A4-LDH is the predominant isoform found in white (glycolytic) muscle of teleosts, and because this muscle is involved in burst swimming during predator avoidance or prey capture, the enzyme should be under strong selective pressure to maintain optimal catalytic properties. A4-LDH active site residues involved in substrate and cofactor binding are invariant from bacteria to vertebrates, and tertiary and quaternary structures are highly conserved across taxa (Auerbach et al. 1998). Researchers also have described the catalytic mechanism of A4-LDH in detail (Holbrook et al. 1975; Holbrook and Gutfreund 1973), including both the process of hydride transfer during reduction of pyruvate to lactate (Deng et al. 1994) and the conformational rearrangements of specific structures during ligand binding and release (Dunn et al. 1991; Gerstein and Chothia 1991).
In this paper we examine how amino acid substitutions associated with cold adaptation affect and Arrhenius activation energy (Ea) in the A4-LDH ortholog of an Antarctic notothenioid fish, Chaenocephalus aceratus. An earlier study (Fields and Somero 1998) showed that A4-LDH orthologs of Antarctic notothenioids (habitat temperature 1.86 ± 2°C) have the highest
values and the highest kcat values of any vertebrate A4-LDH measured. Here we use site-directed mutagenesis to modify residues in C. aceratus A4-LDH that we determined were likely to play a role in cold adaptation of enzyme function. We tested eight nonconservative single mutants, as well as two double mutants, to determine how
and Ea changed in response to modifications in structure. The results reported here provide evidence that (1) enzymes may adapt to cold temperature through reductions in the energy barriers to catalysis, which may or may not involve concomitant reductions in substrate affinity; (2) adaptation to large changes in temperature may be achieved by one to a few amino acid substitutions (in the 332 residues composing the LDH-A monomer) outside the active site; and (3) these changes likely exert their effects by decreasing the energy input needed to alter the conformation of secondary structures that must move during catalysis, i.e., by decreasing activation energy.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Site-Directed Mutagenesis
Residues in the C. aceratus LDH-A monomer targeted for mutagenesis were selected by comparing a consensus of LDH-A amino acid sequences from nine species of Antarctic notothenioids with a consensus sequence developed from six species of warm-adapted teleosts (fig. 1; sequence data from Fields and Somero 1998). The notothenioid consensus was based on LDH-A sequences from Notothenia coriiceps, Lepidonotothen nudifrons, Gobionotothen gibberifrons, Dissostichus mawsoni, Chionodraco rastrospinosus, Chaenocephalus aceratus, Champsocephalus gunnari, Parachaenichthys charcoti, and Harpagifer antarcticus; the C. aceratus LDH-A sequence differs from the Antarctic notothenioid consensus at one position (Leu316Val). The non-notothenioid consensus was based on LDH-A orthologs we sequenced (the eelpout Austrolycos depressiceps, the lampfish Lampanyctus ritteri [unpublished data], and the goby Gillichthys mirabilis [Fields and Somero 1997]) and on sequences obtained from the literature (the killifish Fundulus heteroclitus [Quattro et al. 1995], and the barracudas Sphyraena idiastes and S. lucasana [Holland, McFall-Ngai, and Somero 1997]). The sequence of G. mirabilis LDH-A, which we used as a benchmark warm-adapted ortholog in this study (see below), differs from the non-notothenioid consensus at four positions (Met12Ser, Pro16Ala, Pro310Ala, and Glu311Asp).
|
From the 10 nonconservative substitutions found, we designed eight single mutants (all mutants are given in the cold-adapted [notothenioid] warm-adapted [non-notothenioid] direction): 75His, Ala187Gly, Gly213Thr, Gln217Lys, Gly224Ser, Glu233Met, Ala310Pro, and Gln317Val. In addition, we produced two double mutants (Glu233Met/Ala310Pro and Glu233Met/Gln317Val) by combining mutations that showed effects on A4-LDH kinetics in pilot studies, to determine whether such combinations might reveal synergistic effects among the potentially cold-adaptive substitutions.
Complementary DNA from white epaxial skeletal muscle of C. aceratus was created through reverse transcription of mRNA using Superscript II reverse transcriptase (Promega, Madison, Wis.) and polythymidine primers. The sequences of the 5' and 3' untranslated regions (UTRs) of ldh-a cDNA were determined using RACE PCR (Generacer; Invitrogen, Carlsbad, Calif.) and gene-specific primers originally derived from the sequence of barracuda ldh-a (Holland et al. 1997; for the 5' UTR: 5'-CACTCAGCTTCCAGGCCACATAG-3'; for the 3' UTR: 5'-TACACCTCCTGGGCCATTGGCATGT-3'). After sequences for the UTRs were obtained, primers designed from within the 5' UTR (5'-TCTAACTGAACAACTCAAGATGTCC-3') and 3' UTR (5'-CCCTCCTCCTGTACCTCTTACTG-3') of C. aceratus ldh-a were used to amplify the entire coding sequence. The PCR protocol used for the amplification consisted of 94°C for 2 min, followed by 35 cycles of 94°C for 30 s, 56°C for 45 s, and 72°C for 90 s; a final extension at 72°C for 10 min was followed by rapid cooling to 4°C. Amplified ldh-a cDNA was purified (Qiaquick PCR purification, Qiagen Inc., Valencia, Calif.) and cloned into pAlter-Ex1 (Altered Sites II in vitro mutagenesis system, Promega). The ldh-a insert was mutated following the manufacturer's protocol, using mismatch oligonucleotides prepared commercially (Invitrogen).
Recombinant plasmids were created using a commercial kit (Altered Sites II in vitro Mutagenesis System, Promega) and were used to transform repair-minus ES1301 Escherichia coli (mutS). After growth to log phase, plasmids were purified (Qiaprep miniprep kit, Qiagen) and used to transform E. coli JM109 cells. For each mutant, appropriate antibiotic resistance and automated nucleotide sequencing (BigDye Terminator Cycle Sequencing, Applied Biosystems, Foster City, Calif.) of the ldh-a insert confirmed mutagenesis.
A4-LDH Expression
The pAlter-Ex1 plasmid contains a tac promoter and a ribosome binding site upstream of the insert. Transformed JM109 cells were grown to log phase (OD600 = 0.4000.600) in 800 ml LB broth, and expression of recombinant A4-LDH was induced with IPTG. The culture was incubated for another 612 h, and the cells were collected by centrifugation at 4,000xg for 20 min. Cells were lysed with CelLytic B II (Sigma, St. Louis, M.O.) following the manufacturer's instructions. After centrifugation at 12,000xg and 4°C for 1 h, LDH activity was found in the supernatant. Supernatant was dialyzed (10,000 MWCO) overnight against 50 mM potassium phosphate pH 6.8, and LDH assays were performed without further purification. Activity-stained native polyacrylamide gel electrophoresis (Brewer 1970) indicated that no other isoforms of LDH were present in the supernatants (data not shown).
Measurement of KmPYR and Arrhenius Activation Energy
Kinetics assays were performed on a Shimadzu model UV1601 or a Pharmacia Ultrospec 2100pro UV-Visible spectrophotometer, each with a temperature-controlled cuvette holder (± 0.2°C). Substrate saturation curves were derived for each A4-LDH mutant, as well as for recombinant wild type C. aceratus A4-LDH (rWT), at 0°C, 5°C, 10°C, 15°C, and 20°C. Measurements were also made for A4-LDH purified (Yancey and Somero 1978) from white muscle of C. aceratus (WT) and the warm-temperate goby Gillichthys mirabilis. Three replicates of six pyruvate concentrations were used at each temperature; the cofactor NADH was saturating (150 µM) for all assays. Replicate measurements were made using a single enzyme preparation for each of the samples except recombinant wild type, for which triplicate measurements were made on four independently produced samples, totaling 12 replicates.
For each assay, conversion of NADH to NAD+ was monitored at 340 nm and converted to reaction rate. These data were used to determine three and maximal velocity (Vmax) values at each temperature, using Wilman4 Km-fitting software (Brooks and Suelter 1986) and nonlinear regression analysis (Wilkinson 1961). Activation energy was determined from Vmax data according to the Arrhenius equation, where the slope of the plot (ln [Vmax] vs. 1/K) can be related to Ea by the equation Ea = slope x R, where R is the universal gas constant (8.31441 J mol1 K1).
Statistics
Each value reported represents the mean ± SE for three replicates, except rWT, which comprises 12 replicates. To determine whether orthologs had significantly different
values, an analysis of variance was used (anova1; MATLAB; Mathworks, Natick, Mass.), followed by a Tukey-Kramer multiple comparisons test (multcompare; MATLAB) with
= 0.05. These tests were performed on the
values of each ortholog measured at 10°C. To determine whether Ea values were different among orthologs, an analysis of covariance was performed on the regressions of the Arrhenius plots (ln [Vmax] vs. 1/K). Because slope in this type of plot is proportional to Ea, statistically significant differences in slope indicate orthologs that possess different Ea values. Analysis of covariance was performed with the aoctool module of MATLAB software, and Tukey-Kramer multiple comparisons tests (
= 0.05) were performed with the multcompare module.
Molecular Modeling
To better understand how the mutations found in A4-LDH of the Antarctic notothenioids affect the catalytic process, we created a molecular model using SwissModel software (Guex and Pietsch 1997; Schwede et al. 2003). Dogfish (Abad-Zapatero et al. 1987) and pig (Dunn et al. 1991) A4-LDH structures (PDB accession codes 6LDH and 9LDT, respectively) were used as templates on which the C. aceratus LDH-A primary structure was threaded. Because of the high homology, both in primary and tertiary structure, among vertebrate A4-LDHs, we consider the resultant homology model of C. aceratus A4-LDH to be an accurate representation of the structure of the LDH-A monomer. This model was visualized with Rasmol (Sayle and Milner-White 1995) and DeepView software (Guex and Pietsch 1997) to locate the position of the mutations with respect to the active site and to mobile structures in the molecule involved in subunit interactions and catalysis.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
|
In contrast, two single mutants, shown in figure 2, did effect significant reductions in when compared to WT. Mutant Gln317Val reduced
values across the range 0°C20°C to a level midway between the C. aceratus ortholog and the warm-adapted G. mirabilis form. At 10°C the
of Gln317Val was significantly lower than that of WT, but significantly greater than that of the G. mirabilis ortholog (Fig. 3). The
of the most effective single mutant, Glu233Met, was 0.181 mM pyruvate at 10°C, a value not significantly different from that of the A4-LDH of G. mirabilis (0.151 mM pyruvate), whose habitat temperature can be more than 35°C warmer than that of C. aceratus (Barlow 1961; Fields and Somero 1997). These results indicate that a single mutation is sufficient to alter the substrate affinity of the Antarctic C. aceratus A4-LDH ortholog to a level similar to that of the warm-adapted G. mirabilis, and they suggest that Glu233Met alone, or a combination of Glu233Met and Gln317Val, are responsible for the adaptive increase in
found in notothenioid A4-LDHs.
Two double mutants also were created to determine whether there were synergistic interactions between substitutions in C. aceratus A4-LDH. The mutant Glu233Met/Ala310Pro has values not significantly different from the single mutant Glu233Met (figs. 2 and 3). This suggests that the mutation Ala310Pro exerts no effect on
whether singly or in conjunction with Glu233Met. The double mutant Gln317Val/Glu233Met produced
values between those of the single mutants Gln317Val and Glu233Met but not significantly different from either one (fig. 3).
Effects of Mutations on Arrhenius Activation Energy of Chaenocephalus aceratus A4-LDH
The Arrhenius activation energy (Ea) of each of the mutants was calculated using Vmax values collected from 0°C20°C. Figure 4 gives Ea values ±95% confidence intervals for each enzyme and shows that mutants Gly213Thr and Gly224Ser have Ea values not significantly different from WT or rWT; the other three single mutants (Glu233Met, Ala310Pro, and Gln317Val) and the two double mutants (Glu233Met/Ala310Pro and Glu233Met/Gln317Val), however, have Ea values significantly higher than that of C. aceratus A4-LDH and not significantly different from the warm-adapted G. mirabilis ortholog.
|
|
In sum, these results indicate that the mutant Glu233Met strongly affects both Ea and increasing the former and decreasing the latter; Gln317Val strongly affects Ea and moderately affects
and Ala310Pro strongly affects Ea without affecting
The changesdecreases in substrate affinity and activation energy barriersare in the direction predicted if the substitutions serve to adapt C. aceratus A4-LDH to function at low temperatures. Double mutants of these substitutions have no synergistic effects, and the other mutants tested have no significant effect on Ea or
(except Gly224Ser, which increased
).
Molecular Modeling
A three-dimensional model of one monomer of C. aceratus A4-LDH is shown in figure 6. The mutations that affect the kinetic properties of the enzyme are associated with one of two -helices (
1G-2G and
H) that form the margins of the cleft by which ligands can enter the active site. As discussed below, these helices must move during catalysis (Gerstein and Chothia 1991). Each of the effective mutants is associated with one of these helices, but the mutants that increase Ea without affecting substrate affinity (Gly224Ser and Ala310Pro) are further from the active site on these helices than the substitutions (Glu233Met and Gln317Val) that increase both kinetic parameters.
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The reasoning behind this supposition is that for many enzymes, including A4-LDH (Dunn et al. 1991), the rate-limiting step to catalysis is a conformational shift during binding of ligands or release of product. According to this view, at colder temperatures, where less thermal energy is available in the environment to drive catalytically necessary conformational changes, catalytic rates are maintained at optimal levelsi.e., temperature compensation is producedby decreasing the energy barriers to shifts between conformational states. We further have argued (Fields and Somero 1998; Fields 2001) that results of this increase in conformational entropy in an ortholog adapted to cold compared to a warm-adapted ortholog include: (1) a decrease in Ea and an increase in kcat, as lower energy barriers allow the cold-adapted form to proceed through the catalytic cycle more quickly; and (2) an increase in Km, as increased conformational flexibility causes the cold-adapted ortholog to spend more time in conformational microstates that are not optimal for substrate binding. This model has received support from a recent study by D'Amico, Gerday, and Feller (2001), who found a decrease in Km and kcat (i.e., loss of cold adaptation) in mutants of a psychrophilic -amylase that had been designed to increase thermostability through increased molecular rigidity.
In our earlier study (Fields and Somero 1998) we showed that Antarctic notothenioid A4-LDHs indeed displayed increased and kcat with respect to a group of more warm-adapted teleost orthologs. The Antarctic notothenioids, which spend their entire life cycle at extremely cold temperatures, have A4-LDH orthologs with the highest
and kcat values yet found. Here we have shown that a single amino acid substitution, Glu233Met (notothenioid
non-notothenioid), is sufficient to decrease
and increase Ea of a notothenioid A4-LDH to a level comparable to those of orthologs from much more warm-adapted teleosts.
There is, however, an alternate explanation for the differences we have found in the apparent and Ea of our recombinant A4-LDHs with respect to WT: it is possible that differences in the intracellular milieu between teleost white muscle and E. coli, such as changes in osmolality or osmolyte composition, could lead to modifications in the final conformation of the homotetramer. These subtle structural changes might affect the kinetics of the expressed enzymes, irrespective of the presence of amino acid substitutions. There is growing evidence that A4-LDH kinetics indeed can be affected by the folding environment in vivo, either in individuals of the same species acclimated to different temperatures (e.g., cod, see Zakhartsev et al. 2003; loach, see Ozernyuk, Klyachko, and Polosukhina 1994), or in confamilials whose LDH-A primary structures are identical but whose kinetics are different (Gillichthys spp., see Fields and Somero 1997; Fields et al. 2002). However, in the present study we have shown that WT C. aceratus A4-LDH (that is, enzyme extracted directly from white muscle) has identical
and Ea values to rWT A4-LDH expressed in E. coli (figs. 25
). We take this as evidence that differences in kinetics between WT and our A4-LDH mutants indeed are due to the amino acid substitutions per se, and not to differences in cytosolic composition between E. coli and C. aceratus.
Antarctic notothenioid fishes are an especially appropriate model for testing cold adaptation, because of their rapid and recent invasion of the extremely cold polar marine habitat. Before 25 MYA Antarctica was a part of the Gondwanan supercontinent, attached to South America by what is today the Antarctic peninsula. Before this time, the coastal waters around Antarctica were relatively temperate (
5°C20°C). After 25 MYA, however, Antarctica separated from South America, opening the Drake Passage and allowing formation of the circumpolar current that isolated the region from more temperate waters. From the middle Miocene to the present, water temperatures around the Antarctic continent dropped to 0°C or below (Eastman 1993; Clarke and Johnston 1996), creating a novel environment where selection for cold adaptation was strong. These geological events indicate that the speciation of the Antarctic notothenioids likely occurred around 25 MYA or later, and molecular phylogenies of the Antarctic notothenioid families suggest an even more recent radiation. Comparisons of 12S and 16S mitochondrial rRNA genes indicate that Antarctic notothenioids diverged relatively recently (
1015 MYA) in response to the rapid changes in the thermal environment around Antarctica (Bargelloni et al. 1994; Clarke and Johnston 1996). As a result of this recent speciation, the relatively few amino acid substitutions found in the Antarctic notothenioid LDH-A sequences relative to the non-notothenioid consensus (22/332, or 7%) are likely the products of positive selection for function in the cold polar environment rather than of neutral variation.
If modifications in molecular flexibility indeed are necessary during temperature adaptation, then amino acid substitutions with effects on enzyme kinetics, such as Km and Ea, should be localized to areas of the enzyme that must move during catalysis. Based on this conceptual framework, we (Fields and Somero 1998) hypothesized that temperature-adaptive changes in notothenioid A4-LDH structure were likely to reside near one or more of three specific structures in the molecule that had been described as "major movers" during catalysis (Gerstein and Chothia 1991): the "catalytic loop" region N-terminal to and including part of helix D, the region surrounding helix
1G-2G, and the region encompassing helix
H at the C-terminus of the molecule (see figs. 1 and 6). These regions form the channel through which NADH and pyruvate must pass to access the active site (fig. 6). After ligands are bound, these three regions close together, creating a catalytic vacuole separated from the surrounding medium, in which the hydride transfer takes place. These structures must open again to allow release of lactate and NAD+ and to return the enzyme to the apo form (Abad-Zapatero et al. 1987; Gerstein and Chothia 1991; Dunn et al. 1991).
Examination of figure 1 reveals that the "catalytic loop" region, although it moves the greatest amount during the catalytic cycle (15Å; Gerstein and Chothia 1991), shows no differences between the Antarctic notothenioid and the non-notothenioid teleost LDH-A consensus sequences. Therefore, changes in the flexibility of this region cannot be responsible for cold adaptation in the Antarctic notothenioid orthologs. In contrast, there are a number of substitutions between the two consensus sequences near both the
1G-2G and the
H helices. The ßJ-
1G loop region is directly N-terminal to the
1G-2G helix, a relatively disordered structure (Abad-Zapatero et al. 1987) that interacts with the solvent and may play a role in controlling the mobility of the catalytically important
1G-2G helix (Fields and Somero 1998). The 14 residues in the ßJ-
1G loop immediately N-terminal to the
1G-2G helix contain five differences between the two consensus sequences (36%, compared to 7% in the entire molecule). Of these, two (Gly213Thr and Gly 224Ser) involved replacement of polar residues in the warm-adapted forms with glycyl residues in the notothenioid orthologs. The substitution of glycyl residues in the cold-adaptive forms is suggestive, because addition of these residues increases the flexibility of the peptide backbone (Matthews, Nicholson, and Becktel 1987), and increased flexibility in the ßJ-
1G loop region might confer greater mobility on the presumably rigid
1G-2G helix, allowing it more readily to undergo catalytically necessary conformational shifts. The results of site-directed mutagenesis replacing Gly213 or Gly224 with Thr or Ser, respectively, indicate that Ea of the mutant A4-LDHs indeed increases, suggesting that the loss of flexibility in the ßJ-
1G loop region caused by either of these substitutions is sufficient to increase the activation energy barriers to catalysis. Importantly, however, neither of these mutants decreased
and Gly224Ser instead increased it, suggesting that these mutants either have no effect on substrate binding affinitypresumably because they are far from the active site in a relatively loosely constrained loopor may subtly disrupt the geometry of the binding site through effects on active site residues (e.g., Thr248) at the C-terminal end of helix
1G-2G. A third mutant from the same region, Gln217Lys, was examined in a preliminary study and showed no effect on
at 10°C. Although this mutant, which causes a loss of charge in the cold-adapted direction, might affect mobility of the region through long-range electrostatic effects or interactions with the solvent, preliminary tests suggested this is unlikely.
Interestingly, two other studies on teleost A4-LDHs implicate the ßJ-1G loop as a region important in temperature adaptation. Johns and Somero (2004) examined temperature adaptation in tropical and temperate damselfishes (Chromis spp.) and found that the relatively cold-adapted C. punctipinnis possesses an A4-LDH with lower
and kcat relative to the tropical orthologs and that these changes in kinetics can be attributed to one amino acid difference, Thr219Ala (temperate
tropical). The same substitution occurs in the notothenioids (Thr220Ala in fig. 1), but mutational analyses have not yet been performed on this residue. Holland, McFall-Ngai, and Somero (1997) examined A4-LDHs in temperate barracudas (Sphyraena spp.) and found that adaptive changes in kinetics could be explained by another single substitution, Asp8Gly, in the temperate
subtropical (S. argentea
S. lucasana) direction. This position appears far removed from the ßJ-
1G loop, but in the quaternary structure of the A4-LDH homotetramer the
A helix, in which residue 8 resides, interacts with the ßJ-
1G region of a neighboring monomer. Thus, both of these studies, along with the work presented herein, suggest strongly that the ßJ-
1G loop region is an important site for temperature adaptation in A4-LDH, likely through effects on the mobility of the neighboring
1G-2G helix.
The second mobile helix of A4-LDH, H, also had a substitution at its N-terminal end in notothenioid orthologs, with potential effects on mobility. Ala310Pro is a suggestive substitution, because prolyl residues uniquely constrain rotational freedom through the creation of a pyrrolidine ring between the
carbon and the amide nitrogen of the peptide backbone. The loss of this constraint in notothenioid A4-LDHs, with respect to the non-notothenioid consensus, suggested a mechanism by which helix
H would achieve greater structural mobility and hence reduce the activation energy necessary for catalysis. Indeed, results indicate that replacement of the notothenioid alanyl by the non-notothenioid prolyl at residue 310 increases Ea significantly (figs. 4 and 5). However, as with the addition of glycyl residues in the ßJ-
1G region, this increase in mobility has no significant effect on
of the mutant. Again, this separation of effects on activation energy from binding affinity may be due to the position of Ala310Pro, far enough from the active site to avoid modifying active site geometry but directly associated with a structure, helix
H, that must move during catalysis.
It must be noted, however, that the presence of a prolyl residue at position 310 is not universal among vertebrate A4-LDHs, and a number of orthologs that would be predicted to be "warm-adapted" have an alternate residue at this site. Examples of species with such residues include, among others, humans (310Ser; accession number 1I10_A), the green iguana (Iguana iguana; 310Ser; accession number AAN05098), and goby fishes such as the species used as a warm-adapted benchmark in this paper (G. mirabilis). It is clear, therefore, that despite the effects on of substituting proline for alanine at this position in C. aceratus A4-LDH, the loss of 310Pro alone cannot be considered indicative of cold adaptation in A4-LDH orthologs.
In contrast to Ala310Pro, two other mutants, Gln317Val and Glu233Met, each affected both and Ea of C. aceratus A4-LDH. Gln317Val had a moderate effect on both parameters, whereas Glu233Met increased both to a level comparable to those of an ortholog from the warm-temperate goby G. mirabilis (figs. 35
). Unlike the mutants described above, which were located at the ends of the mobile helices
1G-2G and
H farthest from the active site, Glu233Met and Gln317Val are positioned within each of these helices, respectively. They therefore occur closer to the active site, and they may be in a position to impact the affinity of the enzyme for substrate. Homology modeling (fig. 6) indicates that the side chains of these two residues project into the medium, rather than into the catalytic vacuole. In addition, in no other LDH studied do residues in homologous positions interact with substrate or cofactor. Thus, it is unlikely that either affects
through interactions with pyruvate or NADH, and their effects must be indirect.
A number of mechanisms potentially might explain the effects of Glu233Met and Gln317Val on the and Ea of C. aceratus A4-LDH, including (1) altered interactions of the mobile
-helixes with the aqueous medium; (2) modified interactions between these residues, positioned on either edge of the cleft leading to the active site, and the ligands pyruvate and NADH; or (3) modified long-range electrostatic interactions that lead to changes in the mobility of these structures. The relative importance of these mechanisms is unclear, and further research must be performed to determine why structural changes in these specific locations lead to the changes noted in the
and Ea of A4-LDH.
Conclusion
The data presented in this paper support the hypothesis that amino acid substitutions associated with the mobile helixes 1G-2G and
H are responsible for the cold-adaptive changes in kinetics found in notothenioid A4-LDHs. In addition, these results show that adaptive substitutions relatively distant from the active site, but associated with mobile structures (i.e., Gly213Thr, Gly224Ser, and Ala310Pro), affect activation energy barriers to catalysis, with cold-adapted orthologs having lower Ea values, but do not affect substrate affinity, as measured by
In contrast, substitutions associated with mobile structures and positioned relatively close to the active site (i.e., Glu233Met and Gly 317Val) affect both
and Ea. Further research will be necessary to determine whether similar amino acid substitutions, associated with catalytically important mobile secondary structures, can be found in orthologs of other cytosolic enzymes adapted to extreme cold.
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Abad-Zapatero, C., J. P. Griffith, J. L. Sussman, and M. G. Rossmann. 1987. Refined crystal-structure of dogfish M4 apo-lactate dehydrogenase. J. Mol. Biol. 198:445467.[ISI][Medline]
Auerbach, G., R. Ostendorp, L. Prade, I. Korndorfer, T. Dams, R. Huber, and R. Jaenicke. 1998. Lactate dehydrogenase from the hyperthermophilic bacterium Thermotoga maritima: the crystal structure at 2.1 angstrom resolution reveals strategies for intrinsic protein stabilization. Structure 6:769781.[ISI][Medline]
Bargelloni, L., P. A. Ritchie, T. Patarnello, B. Battaglia, D. M. Lambert, and A. Meyer. 1994. Molecular evolution at subzero temperaturesmitochondrial and nuclear phylogenies of fishes from Antarctica (suborder Notothenioidei), and the evolution of antifreeze glycopeptides. Mol. Biol. Evol. 11:854863.[Abstract]
Barlow, G. W. 1961. Gobies of the genus Gillichthys with comments on the sensory canals as a taxonomic tool. Copeia 1961:423437.
Brewer, G. J., and C. F. Sing. 1970. An introduction to isozyme techniques. Academic Press, New York.
Brooks, S. P. J., and C. H. Suelter. 1986. Estimating enzyme kinetic-parameters: a computer-program for linear-regression and nonparametric analysis. Int. J. Bio. Med. Comput. 19:8999.[ISI]
Clarke, A., and I. A. Johnston. 1996. Evolution and adaptive radiation of Antarctic fishes. Trends Ecol. Evol. 11:212218.[CrossRef][ISI]
Collins, T., M.-A. Meuwis, C. Gerday, and G. Feler. 2003. Activity, stability and flexibility in glycosidases adapted to extreme thermal environments. J. Mol. Biol. 328:419428.[CrossRef][ISI][Medline]
Dahlhoff, E. P., and G. N. Somero. 1993. Kinetic and structural adaptations of cytoplasmic malate dehydrogenases of eastern pacific abalone (genus Haliotis) from different thermal habitats: Biochemical correlates of biogeographical patterning. J. Exp. Biol. 185:137150.
D'Amico, S., P. Claverie, T. Collins, D. Georlette, E. Gratia, A. Hoyoux, M.-A. Meuwis, G. Feller, and C. Gerday. 2002. Molecular basis of cold adaptation. Phil. Trans. R. Soc. Lond. B 357:917925.[CrossRef][ISI][Medline]
D'Amico, S., C. Gerday, and G. Feller. 2001. Structural determinants of cold adaptation and stability in a large protein. J. Biol. Chem. 276:2579125796.
Deng, H., J. Zheng, A. Clarke, J. J. Holbrook, R. Callender, and J. W. Bergner, II. 1994. Source of catalysis in the lactate-dehydrogenase system ground-state interactions in the enzyme-substrate complex. Biochemistry 33:22972305.[ISI][Medline]
Dunn, C. R., H. M. Wilks, D. J. Halsall, T. Atkinson, A. R. Clarke, H. Muirhead, and J. J. Holbrook. 1991. Design and synthesis of new enzymes based on the lactate-dehydrogenase framework. Phil. Trans. R. Soc. Lond. B 322:177184.
Eastman, J. T. 1993. Antarctic fish biology: evolution in a unique environment. Academic Press, New York.
Feller, G. 2003. Molecular adaptations to cold in psychrophilic enzymes. Cell. Mol. Life Sci. 60:648662.[CrossRef][ISI][Medline]
Feller, G., and C. Gerday. 1997. Psychrophilic enzymes: molecular basis of cold adaptation. Cell. Mol. Life Sci. 53:830841.[CrossRef][ISI][Medline]
Fields, P. A. 2001. Protein function at thermal extremes: balancing stability and flexibility. Comp. Biochem. Physiol. A 129:417431.[ISI]
Fields, P. A., Y.-S. Kim, J. F. Carpenter, and G. N. Somero. 2002. Temperature adaptation in Gillichthys (Teleost: Gobiidae) A4-lactate dehydrogenase: identical primary structures produce subtly different conformations. J. Exp. Biol. 205:12931303.
Fields, P. A., and G. N. Somero. 1997. Amino acid sequence differences cannot fully explain interspecific variation in thermal sensitivities of gobiid fish A4-lactate dehydrogenases (A4-LDHs). J. Exp. Biol. 200:18391850.
. 1998. Hot spots in cold adaptation: localized increases in conformational flexibility in lactate dehydrogenase A4 orthologs of Antarctic notothenioid fishes. Proc. Natl. Acad. Sci. USA 95:1147611481.
Georlette, D., B. Damien, V. Blaise, E. Depiereux, V. N. Uversky, C. Gerday, and G. Feller. 2003. Structural and functional adaptations to extreme temperatures in psychrophilic, mesophilic and thermophilic DNA ligases. J. Biol. Chem. 278:3701537023.
Gerike, U., M. J. Danson, and D. W. Hough. 2001. Cold-active citrate synthase: mutagenesis of active-site residues. Protein Eng. 14:655661.[CrossRef][ISI][Medline]
Gerstein, M., and C. Chothia. 1991. Analysis of protein loop closure: 2 types of hinges produce one motion in lactate dehydrogenase. J. Mol. Biol. 220:133149.[ISI][Medline]
Graves, J. E., and G. N. Somero. 1982. Electrophoretic and functional enzymic evolution in 4 species of eastern Pacific barracudas from different thermal environments. Evolution 36:97106.[ISI]
Guex, N., and M. C. Pietsch. 1997. SWISS-MODEL and Swiss-PDB Viewer: an environment for comparative protein modeling. Electrophoresis 18:27142723.[ISI][Medline]
Hochachka, P. W., and G. N. Somero. 2002. Biochemical adaptation: mechanism and process in physiological evolution. Oxford University Press, Oxford.
Holbrook, J. J., and H. Gutfreund. 1973. Approaches to study of enzyme mechanisms: lactate dehydrogenase. FEBS Lett. 31:157169.[CrossRef][ISI][Medline]
Holbrook, J. J., A. Liljas, S. J. Steindel, and M. G. Rossmann. 1975. Lactate dehydrogenase. Pp. 191292 in P. D. Boyer, ed. The enzymes, Vol. XI. 3rd edition. Academic Press, New York.
Holland L. Z., M. McFall-Ngai, and G. N. Somero. 1997. Evolution of lactate dehydrogenase-A homologs of barracuda fishes (genus Sphyraena) from different thermal environments: differences in kinetic properties and thermal stability are due to amino acid substitutions outside the active site. Biochemistry 36:32073215.[CrossRef][ISI][Medline]
Jaenicke, R. 1991. Protein stability and molecular adaptation to extreme conditions. Eur. J. Biochem. 202:715728.[Abstract]
Johns, G. C., and G. N. Somero. 2004. Evolutionary convergence in adaptation of proteins to temperature: A4-lactate dehydrogenases of Pacific damselfishes (Chromis spp.). Mol. Biol. Evol. 21:314320.
Johnston, I. A., and N. J. Walesby. 1977. Molecular mechanisms of temperature adaptation in fish myofibrillar adenosine triphosphatases. J. Comp. Physiol. 19:195206.
Kim, S.-Y., K. Y. Hwang, S.-H. Kim, H.-C. Sung, Y. S. Han, and Y. Cho. 1999. Structural basis for cold adaptation: sequence, biochemical properties, and crystal structure of malate dehydrogenase from a psychrophile Aquaspirillium arcticum. J. Biol. Chem. 274:1176111767.
Lonhienne, T., C. Gerday, and G. Feller. 2000. Psychrophilic enzymes: revisiting the thermodynamic parameters of activation may explain localized flexibility. Biochim. Biophys. Acta 1543:110.[ISI][Medline]
Matthews, B. W., H. Nicholson, and W. J. Becktel. 1987. Enhanced protein thermostability from site-directed mutations that decrease the entropy of unfolding. Proc. Natl. Acad. Sci. USA 84:66636667.[Abstract]
Ozernyuk, O. D., O. S. Klyachko, and E. S. Polosukhina. 1994. Acclimation temperatiure affects the functional and structural properties of lactate dehydrogenase from fish (Misgurnus fossilis) skeletal muscles. Comp. Biochem. Physiol. 107B:141145.[CrossRef]
Quattro, J. M., D. D. Pollock, M. Powell, H. A. Woods, and D. A. Powers. 1995. Evolutionary relations among vertebrate muscle-type lactate dehydrogenases. Mol. Mar. Biol. Biotechnol. 4:224231.[ISI][Medline]
Sayle, R. A., and E. J. Milner-White. 1995. RasMol: biomolecular graphics for all. Trends Biochem. Sci. 20:374376.[CrossRef][ISI][Medline]
Schwede, T., J. Kopp, N. Guex, and M. C. Peitsch. 2003. SWISS-MODEL: an automated protein homology-modeling server. Nucleic Acids Res. 31:33813385.
Somero, G. N. 1995. Proteins and temperature. Ann. Rev. Physiol. 57:4368.[CrossRef][ISI][Medline]
Svingor, A., J. Kardos, I. Hajdu, A. Nemeth, and P. Zavodszky. 2001. A better enzyme to cope with cold: comparative flexibility studies on psychrotrophic, mesophilic, and thermophilic IPMDHs. J. Biol. Chem. 276:2812128125.
Tsigos, I., K. Mavromatis, M. Tzanodaskalaki, C. Pozidis, M. Kokkinidis, and V. Bouriotis. 2001. Engineering the properties of a cold active enzyme through rational redesign of the active site. Eur. J. Biochem. 268:50745080.
Wilkinson, G. N. 1961. Statistical estimations in enzyme kinetics. Biochem. J. 80:323332.
Xu, Y., G. Feller, C. Gerday, and N. Glansdorff. 2003. Metabolic enzymes from psychrophilic bacteria: challenge of adaptation to low temperatures in ornithine carbamoyltransferase form Moritella abyssi. J. Bacteriol. 185:21612168.
Yancey, P. H., and J. F. Siebenaller. 1987. Coenzyme binding ability of homologues of M4-lactate dehydrogenase in temperature adaptation. Biochim. Biophys. Acta 924:483491.[ISI][Medline]
Yancey, P. H., and G. N. Somero. 1978. Urea-requiring lactate-dehydrogenases of marine elasmobranch fishes. J. Comp. Physiol. 125:135141.[ISI]
Zakhartsev, M., T. Johansen, H. O. Portner, and R. Blust. 2003. Effects of temperature acclination on lactate dehydrogenase of cod (Gadus morhua): genetic, kinetic and thermodynamic aspects. J. Exp. Biol. 207:95112.[CrossRef][ISI]
Zar, J. H. 1984. Biostatistical analysis. Prentice-Hall, Inc., Englewood Cliffs, New Jersey.
Zavodszky, P., J. Kardos, A. Svingor, and G. A. Petsko. 1998. Adjustment of conformational flexibility is a key event in the thermal adaptation of proteins. Proc. Natl. Acad. Sci. USA 95:74067411.