Enzymatic Kinetic Properties of the Lactate Dehydrogenase Isoenzyme C4 of the Plateau Pika (Ochotona curzoniae)

Testis-specific lactate dehydrogenase (LDH-C4) is one of the lactate dehydrogenase (LDH) isozymes that catalyze the terminal reaction of pyruvate to lactate in the glycolytic pathway. LDH-C4 in mammals was previously thought to be expressed only in spermatozoa and testis and not in other tissues. Plateau pika (Ochotona curzoniae) belongs to the genus Ochotona of the Ochotonidea family. It is a hypoxia-tolerant species living in remote mountain areas at altitudes of 3000–5000 m above sea level on the Qinghai-Tibet Plateau. Surprisingly, Ldh-c is expressed not only in its testis and sperm, but also in somatic tissues of plateau pika. To shed light on the function of LDH-C4 in somatic cells, Ldh-a, Ldh-b, and Ldh-c of plateau pika were subcloned into bacterial expression vectors. The pure enzymes of Lactate Dehydrogenase A4 (LDH-A4), Lactate Dehydrogenase B4 (LDH-B4), and LDH-C4 were prepared by a series of expression and purification processes, and the three enzymes were identified by the method of sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and native polyacrylamide gel electrophoresis (PAGE). The enzymatic kinetics properties of these enzymes were studied by Lineweaver-Burk double-reciprocal plots. The results showed the Michaelis constant (Km) of LDH-C4 for pyruvate and lactate was 0.052 and 4.934 mmol/L, respectively, with an approximate 90 times higher affinity of LDH-C4 for pyruvate than for lactate. At relatively high concentrations of lactate, the inhibition constant (Ki) of the LDH isoenzymes varied: LDH-A4 (Ki = 26.900 mmol/L), LDH-B4 (Ki = 23.800 mmol/L), and LDH-C4 (Ki = 65.500 mmol/L). These data suggest that inhibition of lactate by LDH-A4 and LDH-B4 were stronger than LDH-C4. In light of the enzymatic kinetics properties, we suggest that the plateau pika can reduce reliance on oxygen supply and enhance its adaptation to the hypoxic environments due to increased anaerobic glycolysis by LDH-C4.


Introduction
The Qinghai-Tibet Plateau, with an average elevation of over 3000 m above sea level, is the highest and also largest plateau on earth and is characterized by extremely harsh environmental conditions such as hypoxia, strong ultraviolet radiation, and cold temperature which expectedly have profound effects on human and animal survival. Over long-term adaption and evolution, both plateau-native humans and animals have developed a series of strategies to adapt to the plateau environment. Also, inhabitants have developed some unique features to avoid disadvantageous environmental aspects and effectively obtain material or energy from their survival environment as to ensure their normal growth and breeding. and Ldh-c, suggesting that the expression plasmids pCold-SUMO-Ldh-a, pET-30a-Ldh-b, and pET-30a-Ldh-c encoding the full-length pika LDHA, LDHB and LDHC, respectively, were successfully constructed.
The SDS-PAGE results of cell-free extracts (shown in Figure 2A-C) displayed the protein in the supernatant and precipitate fractions with dark bands around 40-50 kD. The recombinant LDH isoenzymes were expressed with His-tag (LDHA with His-tag and SUMO-tag) in the C-terminus, which was favorable to purification. These results demonstrated similar molecular weights to those of LDHA, LDHB, and LDHC, as we previously reported [17], indicating that expression plasmids pCold-SUMO-Ldh-a, pET-30a-Ldh-b and pET-30a-Ldh-c were successfully prepared, and expressed the proteins of interest correctly.

LDH-A4, LDH-B4 and LDH-C4 Purification and Activity Measurement
The native polyacrylamide gel electropheresis (native PAGE) results shown in Figure 3P1 indicate that LDH-A4 was observed in 100, 250 and 500 mmol/L imidazole buffer elutions with LDH activity, and LDH-B4 and LDH-C4 were observed in 250 and 500 mmol/L imidazole buffers elutions with LDH activity. LDH-A4, LDH-B4 and LDH-C4 bound so completely to Ni-NTA that no active protein was assayed in the remainder of the purified supernatant. As shown in Figure 3P2, the purified LDH-A4, LDH-B4, and LDH-C4, were represented with only one band as observed by SDS-PAGE, and the molecular weights were consistent with LDHA, LDHB, and LDHC, respectively. Figure 3P3 presents the LDH activity of the purified LDH-A4, LDH-B4, and LDH-C4 as determined by native PAGE. These results indicate that LDH-A4, LDH-B4, and LDH-C4 of plateau pika were successfully purified. As shown in Figure 4, the respective enzyme activity (n = 5) of LDH-A4, LDH- Ldh-c, suggesting that the expression plasmids pCold-SUMO-Ldh-a, pET-30a-Ldh-b, and pET-30a-Ldh-c encoding the full-length pika LDHA, LDHB and LDHC, respectively, were successfully constructed.
The SDS-PAGE results of cell-free extracts (shown in Figure 2A-C) displayed the protein in the supernatant and precipitate fractions with dark bands around 40-50 kD. The recombinant LDH isoenzymes were expressed with His-tag (LDHA with His-tag and SUMO-tag) in the C-terminus, which was favorable to purification. These results demonstrated similar molecular weights to those of LDHA, LDHB, and LDHC, as we previously reported [17], indicating that expression plasmids pCold-SUMO-Ldh-a, pET-30a-Ldh-b and pET-30a-Ldh-c were successfully prepared, and expressed the proteins of interest correctly.  and Ldh-c, suggesting that the expression plasmids pCold-SUMO-Ldh-a, pET-30a-Ldh-b, and pET-30a-Ldh-c encoding the full-length pika LDHA, LDHB and LDHC, respectively, were successfully constructed.
The SDS-PAGE results of cell-free extracts (shown in Figure 2A-C) displayed the protein in the supernatant and precipitate fractions with dark bands around 40-50 kD. The recombinant LDH isoenzymes were expressed with His-tag (LDHA with His-tag and SUMO-tag) in the C-terminus, which was favorable to purification. These results demonstrated similar molecular weights to those of LDHA, LDHB, and LDHC, as we previously reported [17], indicating that expression plasmids pCold-SUMO-Ldh-a, pET-30a-Ldh-b and pET-30a-Ldh-c were successfully prepared, and expressed the proteins of interest correctly.

LDH-A4, LDH-B4 and LDH-C4 Purification and Activity Measurement
The native polyacrylamide gel electropheresis (native PAGE) results shown in Figure 3P1 indicate that LDH-A4 was observed in 100, 250 and 500 mmol/L imidazole buffer elutions with LDH activity, and LDH-B4 and LDH-C4 were observed in 250 and 500 mmol/L imidazole buffers elutions with LDH activity. LDH-A4, LDH-B4 and LDH-C4 bound so completely to Ni-NTA that no active protein was assayed in the remainder of the purified supernatant. As shown in Figure 3P2, the purified LDH-A4, LDH-B4, and LDH-C4, were represented with only one band as observed by SDS-PAGE, and the molecular weights were consistent with LDHA, LDHB, and LDHC, respectively. Figure 3P3 presents the LDH activity of the purified LDH-A4, LDH-B4, and LDH-C4 as determined by native PAGE. These results indicate that LDH-A4, LDH-B4, and LDH-C4 of plateau pika were successfully purified. As shown in Figure 4, the respective enzyme activity (n = 5) of LDH-A4, LDH-B4, and LDH-C4 was 2241 ± 9.3 U/g protein, 3918 ± 37.9 U/g protein, and 10,741 ± 80.9 U/g protein,

LDH-A 4 , LDH-B 4 and LDH-C 4 Purification and Activity Measurement
The native polyacrylamide gel electropheresis (native PAGE) results shown in Figure 3P 1 indicate that LDH-A 4 was observed in 100, 250 and 500 mmol/L imidazole buffer elutions with LDH activity, and LDH-B 4 and LDH-C 4 were observed in 250 and 500 mmol/L imidazole buffers elutions with LDH activity. LDH-A 4 , LDH-B 4 and LDH-C 4 bound so completely to Ni-NTA that no active protein was assayed in the remainder of the purified supernatant. As shown in Figure 3P 2 , the purified LDH-A 4 , LDH-B 4 , and LDH-C 4 , were represented with only one band as observed by SDS-PAGE, and the molecular weights were consistent with LDHA, LDHB, and LDHC, respectively. Figure 3P 3 presents the LDH activity of the purified LDH-A 4 , LDH-B 4 , and LDH-C 4 as determined by native PAGE.

Enzyme Kinetics Properties of LDH-A4, LDH-B4, and LDH-C4
The Km values for LDH-A4, LDH-B4 and LDH-C4 were determined from Lineweaver-Burk double-reciprocal plots. The activities of the LDH isoenzymes in forward reaction (pyruvate + NADH + H + → lactate + NAD + ) increased with increasing concentration of the substrate pyruvate, as shown in Figure 5. The kinetic parameters of the forward reaction for LDH-A4, LDH-B4, and LDH-C4, calculated from the Lineweaver-Burk double-reciprocal plots of reciprocal pyruvate concentration vs. reciprocal velocity, had Km values of 0.260, 0.172, and 0.052 mmol/L, respectively. The activities of LDH isoenzymes in reverse reaction (lactate + NAD + → pyruvate + NADH + H + ) also increased with The enzyme activity of LDH-A 4 , LDH-B 4 , and LDH-C 4 was 2241˘9.3 U/g protein, 3918˘37.9 U/g protein, 10741˘80.9 U/g protein, each sample size was 5, ** p < 0.01. The enzyme activity of LDH-C 4 was significantly higher that of LDH-A 4 and LDH-B 4 (p < 0.01).

Enzyme Kinetics Properties of LDH-A 4 , LDH-B 4 , and LDH-C 4
The Km values for LDH-A 4 , LDH-B 4 and LDH-C 4 were determined from Lineweaver-Burk double-reciprocal plots. The activities of the LDH isoenzymes in forward reaction (pyruvate + NADH + H + Ñ lactate + NAD + ) increased with increasing concentration of the substrate pyruvate, as shown in Figure 5. The kinetic parameters of the forward reaction for LDH-A 4 , LDH-B 4 , and LDH-C 4 , calculated from the Lineweaver-Burk double-reciprocal plots of reciprocal pyruvate concentration vs. reciprocal velocity, had Km values of 0.260, 0.172, and 0.052 mmol/L, respectively. The activities of LDH isoenzymes in reverse reaction (lactate + NAD + Ñ pyruvate + NADH + H + ) also increased with increasing concentration of the substrate lactate, as shown in Figure 6. The kinetic parameters of the reverse reaction for LDH-A 4 , LDH-B 4 , and LDH-C 4 , calculated from the Lineweaver-Burk double-reciprocal plots of reciprocal lactate concentration vs.   The double reciprocal plots of initial velocities at different pyruvate concentrations and the inhibitory effect of lactate at different concentrations on pika LDH isozymes are shown in Figure 7. In all cases, the competitive inhibition was observed. Initial velocities were determined in the presence and absence of lactate of different concentrations. The Ki values calculated from the Lineweaver-Burk double-reciprocal plots of reciprocal pyruvate concentration vs. reciprocal velocity. Reciprocals of initial velocity with substrate concentrations below the above-mentioned concentrations gave linear plots with these substrates. The results showed that the Ki of lactate for LDH-C 4 (Ki = 65.500 mmol/L) was almost three times higher than that for LDH-A 4 (Ki = 26.900 mmol/L) and LDH-B 4 (Ki = 23.800 mmol/L). The inhibitory effects of lactate at different concentrations on the enzyme reaction velocity of LDH isozymes are shown in Figure 8. Under the experimental conditions employed in the present study, lactate at higher concentrations exhibited a strong inhibition of LDH-A 4 and LDH-B 4 activity, but only slightly inhibited LDH-C 4 activity. The results indicated that LDH-C 4 was less inhibited by lactate, which was likely beneficial in catalyzing the conversion of pyruvate to lactate even at high concentrations of lactate.   The double reciprocal plots of initial velocities at different pyruvate concentrations and the inhibitory effect of lactate at different concentrations on pika LDH isozymes are shown in Figure 7. In all cases, the competitive inhibition was observed. Initial velocities were determined in the presence and absence of lactate of different concentrations. The Ki values calculated from the Lineweaver-Burk double-reciprocal plots of reciprocal pyruvate concentration vs. reciprocal velocity. Reciprocals of initial velocity with substrate concentrations below the above-mentioned concentrations gave linear plots with these substrates. The results showed that the Ki of lactate for LDH-C4 (Ki = 65.500 mmol/L) was almost three times higher than that for LDH-A4 (Ki = 26.900 mmol/L) and LDH-B4 (Ki = 23.800 mmol/L). The inhibitory effects of lactate at different concentrations on the enzyme reaction velocity of LDH isozymes are shown in Figure 8. Under the experimental conditions employed in the present study, lactate at higher concentrations exhibited a strong inhibition of LDH-A4 and LDH-B4 activity, but only slightly inhibited LDH-C4 activity. The results indicated that LDH-C4 was less inhibited by lactate, which was likely beneficial in catalyzing the conversion of pyruvate to lactate even at high concentrations of lactate.

Discussion
LDH is the terminal enzyme of glycolysis with the role of reducing pyruvate to lactate with limiting oxygen [43]. LDH-C4 catalyzes the terminal reaction in the glycolytic pathway and displays a unique structure and functional properties [43,44].
Ultimately, the biochemical properties that distinguish LDH-C4 from the other LDH isoforms may be attributed to its high glycolytic rate. Studies have shown that since the same chemical reaction is being catalyzed, in spite of different kinetic constants, the equilibrium constant (Keq) is the same for all LDH isoenzymes, as stated in the Haldane equation [45,46]. The complex Haldane kinetics mechanism may be reduced to the Michaelis-Menten mechanism more tractably under special conditions of substrate-enzyme interactions [46]. Also, the basic Michaelis-Menten description can be modified in a logistical manner to fit the reaction parameters with in vitro conditions for the enzymecatalyzed reactions [47]. In the sense it is even hyperbolic to logistic kinetics that kinetic constants may be considered as being the same, however with the amendment they would respectively correspond to various fractions from maximum velocity of outtake [45,47]. Based on this theory, the LDH isoenzyme pattern in itself may not change the equilibrium of the LDH reaction or influence the tissue lactate concentration [45]. However, the situation could be quite different during transition in energy metabolism where especially the glycolytic flux can undergo great rapid changes. During fast metabolic transitions involving significant pyruvate concentrations changes, primarily resulting from major changes glycolytic rate, the LDH isoenzyme pattern may be of great importance in the compound metabolic response to the altered energy metabolism [45].Therefore, we presumed that ATP was produced rapidly through anaerobic glycolysis by LDH-C4 in plateau pika, which stirs up the metabolic response to the changed energy metabolism, enhancing its adaptation to the hypoxic environment.
Previous studies have investigated in detail the enzymatic kinetics properties of LDH-C4 in other species [31,48,49]. The biochemical properties distinguishing LDH-C4 from the other LDH isoforms may be conducible to the high glycolytic rate. Compared with LDH-A4, pika LDH-C4 has a high Km for lactate (~2.0 mmol/L) and a low Km for pyruvate (~0.030 mmol/L) [31,[48][49][50]. This finding suggests that LDH-C4 has a 60-fold higher affinity for pyruvate than that for lactate. The finding also implies that pyruvate conversion to lactate may occur even at high concentrations of extracellular or endogenous lactate [31,[48][49][50]. This theory has been supported by an experiment that addition of excess lactate in the environment (50-fold excess in relation to pyruvate) did not affect adenosine triphosphate (ATP) production in capacitating spermatozoa [51]. In the present study, to investigate

Discussion
LDH is the terminal enzyme of glycolysis with the role of reducing pyruvate to lactate with limiting oxygen [43]. LDH-C 4 catalyzes the terminal reaction in the glycolytic pathway and displays a unique structure and functional properties [43,44].
Ultimately, the biochemical properties that distinguish LDH-C 4 from the other LDH isoforms may be attributed to its high glycolytic rate. Studies have shown that since the same chemical reaction is being catalyzed, in spite of different kinetic constants, the equilibrium constant (K eq ) is the same for all LDH isoenzymes, as stated in the Haldane equation [45,46]. The complex Haldane kinetics mechanism may be reduced to the Michaelis-Menten mechanism more tractably under special conditions of substrate-enzyme interactions [46]. Also, the basic Michaelis-Menten description can be modified in a logistical manner to fit the reaction parameters with in vitro conditions for the enzyme-catalyzed reactions [47]. In the sense it is even hyperbolic to logistic kinetics that kinetic constants may be considered as being the same, however with the amendment they would respectively correspond to various fractions from maximum velocity of outtake [45,47]. Based on this theory, the LDH isoenzyme pattern in itself may not change the equilibrium of the LDH reaction or influence the tissue lactate concentration [45]. However, the situation could be quite different during transition in energy metabolism where especially the glycolytic flux can undergo great rapid changes. During fast metabolic transitions involving significant pyruvate concentrations changes, primarily resulting from major changes glycolytic rate, the LDH isoenzyme pattern may be of great importance in the compound metabolic response to the altered energy metabolism [45].Therefore, we presumed that ATP was produced rapidly through anaerobic glycolysis by LDH-C 4 in plateau pika, which stirs up the metabolic response to the changed energy metabolism, enhancing its adaptation to the hypoxic environment.
Previous studies have investigated in detail the enzymatic kinetics properties of LDH-C 4 in other species [31,48,49]. The biochemical properties distinguishing LDH-C 4 from the other LDH isoforms may be conducible to the high glycolytic rate. Compared with LDH-A 4 , pika LDH-C 4 has a high Km for lactate (~2.0 mmol/L) and a low Km for pyruvate (~0.030 mmol/L) [31,[48][49][50]. This finding suggests that LDH-C 4 has a 60-fold higher affinity for pyruvate than that for lactate. The finding also implies that pyruvate conversion to lactate may occur even at high concentrations of extracellular or endogenous lactate [31,[48][49][50]. This theory has been supported by an experiment that addition of excess lactate in the environment (50-fold excess in relation to pyruvate) did not affect adenosine triphosphate (ATP) production in capacitating spermatozoa [51]. In the present study, to investigate the enzyme kinetics of LDH-A 4 , LDH-B 4 , and LDH-C 4 of plateau pika, the expression plasmids pCold-SUMO-Ldh-a, pET-30a-Ldh-b, and pET-30a-Ldh-c, encoding the full-length of pika LDHA, LDHB, and LDHC, were constructed and expressed. SDS-PAGE was used to determine the molecular weight and confirm the recombinant proteins. Figure 2A-C demonstrates that their molecular weights were about 50 kD, 40-45 kD and 40-45 kD for LDHA, LDHB, and LDHC, respectively. According to our previous paper [17], the open reading frames (ORF) of pika Ldh-a, Ldh-b and Ldh-c (accession numbers HQ704676, HQ704677, and HQ704678 in GenBank) were 999 bp, 1005 bp and 999 bp, respectively, encoding 332, 334 and 332 amino acids proteins, with their molecular weights of 36,557.5, 36,464.3 and 36,052.9 Da. Thus the molecular weight of the recombinant LDHA, LDHB, and LDHC in the current study (Figure 2A-C, Figure 3P 2 ) was higher than the value calculated from their amino acid sequences deduced from the encoding cDNA. The addition of the His-tag resulted in a 5-10 kD deviation in the molecular weight of fusion proteins as shown in SDS-PAGE, and also indicated in other studies [52][53][54]. In addition, fusion proteins with the SUMO-tag added 12 kD to its molecular weight [55].
The recombinant expression shuttles (pET-30a-Ldh-b, pET-30a-Ldh-c, pCold-SUMO-Ldh-a) were transformed into E. coli BL21 cells and the enzymatic proteins (LDH-A 4 , LDH-B 4 , and LDH-C 4 ) with functional activity were purified from protein fluid extracts. Our results indicate that the enzyme activity of the pika LDH-C 4 was significantly higher than that of A 4 and B 4 , that LDH-C 4 was less sensitive to lactate inhibition than A 4 and B 4 (p < 0.01), and the Km values of LDH-C 4 for pyruvate and lactate were also higher than that of A 4 and B 4 . Compared with LDH-A 4 and LDH-B 4 , LDH-C 4 had a low Km for pyruvate (~0.052 mmol/L) and a high Km for lactate (~4.934 mmol/L); and the affinity of LDH-C 4 for pyruvate is 90-fold higher than that for lactate. These properties of pika LDH-C 4 were beneficial to catalyze the conversion of pyruvate to lactate even at high concentration of lactate.
Although the conversion mediated by LDH does not generate ATP, the reaction depends on NADH as a cofactor that accompanies oxidation to NAD + ; the concentration of NAD + is a rate-limiting factor in glycolysis and necessary for continued glycolysis [56]. Previous studies demonstrated that knockout of Ldh-c or inhibition of LDH-C 4 in sperm led to rapid decline in sperm ATP levels [56,57], a decrease in progressive motility, and a failure to develop hyper-activated motility. It was revealed that all consumed 13 C-labeled pyruvate added in sperm culture medium was converted to lactate rather than oxidized in the tricarboxylic acid (TCA) cycle in metabolic tracing experiments; in the presence of exogenous pyruvate, the ATP concentration was increased by more than 50% [51]. When carbonyl cyanide m-chlorophenylhydrazone (CCCP) and NaCN was applied to inhibit the oxidative phosphorylation in mitochondria, the amount of ATP was kept at the equivalent level to that without CCCP and the vigorous motility of sperm was maintained [51,58]. These results suggest that LDH-C 4 is essential in sperm glycolysis which has an important role in providing abundant ATP for sperm motility [56,58,59], and it is related to LDH-C 4 enzymatic kinetics properties.
Plateau pika has a strong adaptability to the hypoxic plateau environment of the Qinghai-Tibet Plateau, given the sperm-specific lactate dehydrogenase (Ldh-c) gene is also expressed in their somatic cells. In light of the enzymatic kinetics properties, we propose that the pika could reduce dependence on oxygen and enhance the adaptation to the hypoxic environments due to increased anaerobic glycolysis by LDH-C 4 .

Reagents and Animal Procedures
All reagents were from Sangon (Sangon, Shanghai, China) unless otherwise noted. Plateau pikas were live-trapped from Haibei Alpine Meadow Ecosystem Research Station at an altitude of 3800 m in Qinghai Province, China. All animals were first anesthetized with sodium pentobarbital (5%), and then sacrificed by cervical dislocation immediately before dissection. Skeletal muscle, kidney, and testis were rapidly removed and frozen in liquid nitrogen for long-term storage. All procedures involved in the handling and care of animals were approved by the China Zoological Society according to the China Practice for the Care and Use of Laboratory Animals (permit number: GB 14923-2010).

Protein Expression, Purification, Analysis and Enzyme Activity Measurement
E. coli BL21 (DE3) cells were cultured in media of Lysogeny broth (LB). E. coli clones with plasmids (pCold-SUMO-Ldh-a, pET-30a-Ldh-b, and pET-30a-Ldh-c) were inoculated into a liter of medium with 100 µg/mL ampicillin or 50 µg/mL kanamycin and cultured at 37˝C. When A600 value arrived at 0.6, isopropyl-β-D-thiogalactopyranoside (IPTG) reagent was added till the final concentration of 1 mmol/L, followed by further culturing at 25˝C for 10 h to induce the expression of Ldh-a, Ldh-b, and Ldh-c. To determine LDHA, LDHB and LDHC expression in the supernatant and precipitate of E. coli cells, bacterial cells were collected by centrifugation and resuspended in 20 mmol/L phosphate buffer saline (PBS) (pH 8.0), and disrupted by ultrasonication (36ˆ5 s pulses with 5 s intervals). The supernatant was then collected; the precipitate was resuspended and washed in 20 mmol/L PBS (pH 8.0) with 8 mol/L urea. SDS-PAGE was performed to assay the collected supernatant and treated precipitate with a steady current voltage electrophoresis instrument (Beijing Liuyi Instrument Factory, Beijing, China) as follows: 6 µL samples were loaded on a 12% (w/v) separating gel with a 5% (w/v) stacking gel in Tris-glycine buffer (pH = 8.3). The electric voltage was 80 V in the stacking gel and 120 V in the separating gel. The gel was stained in with Coomassie staining solution (10% methyl alcohol, 10% acetic acid, 0.2% R-250) for 4 h, and finally stored in the destaining solution (10% methyl alcohol, 10% acetic acid).
To obtain the active enzyme protein in the supernatant of E. coli cells, bacterial cells were collected with following procedures: 4000 r/min centrifugation (15 min), resuspension in 20 mmol/L PBS (pH = 7.0), and disruption by repeated freeze-thaw cycles (in liquid nitrogen and 37˝C water bath) for 10 times. In freeze-thaw cycles, Triton X-100 (final concentration of 1%), lysozyme (final concentration of 1%) and DNA enzyme (final concentration of 0.2%) were added. The lysates were centrifuged at 15,000 r/min under a temperature of 4˝C for 10 min, and the obtained supernatant was used for purification. Ni-NTA resin (QIAGEN, Hilden, Germany) was used to purify the recombinant LDHs. The resins were washed twice with increasing concentrations of imidazole at 10, 20, 50, 100, 250 and 500 mmol/L in PBS gradient (50 mmol/L NaH 2 PO 4 , 300 mmol/L NaCl, pH = 8.0). Each fraction was collected using individual collection tubes and detected by native PAGE. Elutions with LDH activity were totally collected. Enzyme liquids were obtained from elutions after removing imidazole by ultrafiltration (Ultrafiltration Tubes, 50 kD, 15 mL, Merck Millipore, Billerica, MA, USA). On average, a final amount of 3 mg LDH isoenzymes were purified from every liter of bacterial culture following this approach. LDH isoenzymes stocks were finally stored in PBS (pH = 8.0) at 4˝C until following steps. The LDH isoenzymes were prepared simultaneously under the same conditions [60].
SDS-PAGE was used to detect the protein purification and native PAGE was used to test and verify the enzyme activity. Native PAGE experiment was also performed with a voltage electrophoresis apparatus (Beijing Liuyi Instrument Factory, Beijing, China) with the following conditions: 8% (w/v) separating gel with a 4% (w/v) stacking gel (non-denaturing polyacrylamide gel) in Tris-glycine (pH = 8.3), and 6 µL samples were loaded. LDH bands were stained at 37˝C in a specific staining reagent (a mixture of 0.5 mol/L phosphate buffer (pH = 7.5), 10 mL of 1 mg/mL nitrobenzenethiocyanate chloride (NBT), 4 mL of 5 mg/mL NAD + , 2.5 mL of 1 mol/L sodium lactate, 2.5 mL of 0.1 mol/L NaCl and 1 mL of 1 mg/mL phenazinemethosulfate (PMS)) for half an hour in the dark. The stained gels were then rinsed with ddH 2 O and stored in the storage liquid (10% glycerol and 7% acetic acid). The reference marker for native PAGE was made with a mixture of pika kidney and skeletal muscle, tissues were homogenized on ice as a 1:4 (w/v) dilution in 0.9% physiological saline, and centrifuged at 4000 r/min at 4˝C for 10 min, and finally the supernatant was collected.

LDH-A 4 , LDH-B 4 and LDH-C 4 Kinetics Properties
The activity of LDH was measured using UV spectrophotometer (Unicon 2800, Shanghai, China) by recording the change of 340 nm absorbance produced by NADH oxidation. Assays were kept at consistent temperature of 37˝C. Calculation of Km value for pyruvate and lactate was resulted from Lineweaver-Burk plots. Ki value of lactate was determined from Km and Vmax, which was obtained with or without inhibitor added to reaction buffer via various concentrations of pyruvate at a constant inhibitor concentration, followed by plotting the slope (Km/Vmax value) against inhibitor concentrations. The mixed reagent of reaction to from pyruvate lactate contained 0.15 mmol/L NADH, 50 mmol/L PBS (pH8.0) and the sodium pyruvate gradients as substrate were performed at 0.05, 0.1, 0.2, 0.4, 0.6 and 0.8 mmol/L for LDH-A 4 , LDH-B 4 , and LDH-C 4 , and the enzyme persisted at a constant inhibitor concentration. The reaction from lactate to pyruvate contained 0.5 mmol/L NAD + , 50 mmol/LPBS, (pH 8.0); sodium lactate gradients as substrate were performed at 2, 5, 10, 20, 30 and 40 mmol/L for LDH-A 4 , LDH-B 4 , and LDH-C 4 . The concentration gradients of sodium lactate as inhibitor were added as 10, 20, and 40 mmol/L, respectively; and the substrates, coenzymes, and other buffers followed the method of the reaction from pyruvate to lactate as outlined above. A ∆E 340 value of 0.06-0.07 per minute was provided in the experiment. The enzyme preparation was diluted with PBS (pH = 8.0) and assayed in a 1 cm light path. Before starting the reaction by addition of substrate, the coenzyme was incubated with the buffer used in the assay for 10 min at 37˝C. In total, 3 mL volumes of reagents were added in the process [61]. All chemicals were analytically pure from Sangon (Sangon, Shanghai, China). All measurements in the study were repeated 5 times for the average values.

Statistics
The data are presented as mean˘standard deviation (SD). Statistical analysis was performed by one-way analysis of variance (ANOVA) followed by Duncan's test using SPSS 22.0 (SPSS Inc., Chicago, IL, USA). Values of p < 0.01 was considered highly significant and p < 0.05 was considered statistically significant.

Conclusions
In light of the enzymatic kinetics properties, we suggest that the pika could reduce dependence on oxygen and enhance adaptation to the hypoxic environment due to its increasing anaerobic glycolysis to produce ATP rapidly by LDH-C 4, which is often the role of LDH-A 4 in most mammal species.