Archaea
Volume 2017 (2017), Article ID 5793620, 6 pages
https://doi.org/10.1155/2017/5793620

Research Article

Expression, Purification, and Characterization of (R)-Sulfolactate Dehydrogenase (ComC) from the Rumen Methanogen Methanobrevibacter millerae SM9

Yanli Zhang, Linley R. Schofield, Carrie Sang, Debjit Dey, and Ron S. Ronimus

AgResearch Limited, Grasslands Research Centre, Tennent Drive, Private Bag 11008, Palmerston North 4442, New Zealand

Correspondence should be addressed to Linley R. Schofield; linley.schofield@agresearch.co.nz

Received 10 May 2017; Accepted 19 September 2017; Published 6 November 2017

Academic Editor: Nils-Kåre Birkeland

Copyright © 2017 Yanli Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

(R)-Sulfolactate dehydrogenase (EC 1.1.1.337), termed ComC, is a member of an NADH/NADPH-dependent oxidoreductase family of enzymes that catalyze the interconversion of 2-hydroxyacids into their corresponding 2-oxoacids. The ComC reaction is reversible and in the biosynthetic direction causes the conversion of (R)-sulfolactate to sulfopyruvate in the production of coenzyme M (2-mercaptoethanesulfonic acid). Coenzyme M is an essential cofactor required for the production of methane by the methyl-coenzyme M reductase complex. ComC catalyzes the third step in the first established biosynthetic pathway of coenzyme M and is also involved in methanopterin biosynthesis. In this study, ComC from Methanobrevibacter millerae SM9 was cloned and expressed in Escherichia coli and biochemically characterized. Sulfopyruvate was the preferred substrate using the reduction reaction, with 31% activity seen for oxaloacetate and 0.2% seen for α-ketoglutarate. Optimal activity was observed at pH 6.5. The apparent KM for coenzyme (NADH) was 55.1 μM, and for sulfopyruvate, it was 196 μM (for sulfopyruvate the Vmax was 93.9 μmol min⁻¹ mg⁻¹ and kcat was 62.8 s⁻¹). The critical role of ComC in two separate cofactor pathways makes this enzyme a potential means of developing methanogen-specific inhibitors for controlling ruminant methane emissions which are increasingly being recognized as contributing to climate change.

1. Introduction

Coenzyme M is an essential cofactor for the final reaction in the methanogenesis pathway and in the production of methane catalyzed by methyl-coenzyme M reductase (MCR) [1]. It is the smallest known cofactor, acts during catalysis as a nucleophile [1, 2], and can reach millimolar concentrations within the cell [3, 4]. In Methanobrevibacter ruminantium, the uptake of CoM is an energy-dependent reaction [3, 5].

Recently, it has been found that the biosynthesis of coenzyme M occurs through two pathways that differ in the steps leading to the production of L-sulfopyruvate, the product of the third step of the canonical pathway catalyzed by ComC [2]. In the originally characterized pathway found in the orders Methanococcales, Methanobacteriales, and Methanopyrales, four enzymes have been characterized thus far [2, 6]. These are phosphosulfolactate synthase (ComA; EC 4.4.1.19), 2-phosphosulfolactate phosphatase (ComB; EC 3.1.3.71) [7], sulfolactate dehydrogenase (ComC; EC 1.1.1.337 (formerly EC 1.1.1.272)) [8], and sulfopyruvate decarboxylase (ComDE; EC 4.1.1.79) [9]. The substrates for ComA are phosphoenolpyruvate and sulfite, and the enzyme produces (R)-phosphosulfolactate, which is then dephosphorylated by ComB to produce (R)-sulfolactate. ComC oxidizes (R)-sulfolactate to sulfopyruvate, which is then decarboxylated by ComDE to produce sulfoacetaldehyde. ComC is related by sequence to lactate/malate dehydrogenases, N-methyl-L-amino acid dehydrogenases, 2,3-diketo-L-gulonate reductases, ureidoglycolate dehydrogenases, and an uncharacterized clade of thermophilic archaeal proteins [8, 10–12]. ComC is also likely to participate in the biosynthesis of methanopterin through production of (S)-hydroxyglutaric acid, which is a component of methanopterin [11].

When the Methanocaldococcus jannaschii genome sequence became available, two malate dehydrogenase genes (MJ1425 and MJ0490) were annotated, and the two enzymes were expressed and characterized by Graupner et al. [11]. Although both enzymes possessed malate dehydrogenase activity (converting oxaloacetate to (S)-malate) [13], MJ1425 was identified as likely to be the biologically relevant ComC due to its kinetic parameters (higher Vmax and Vmax/KM for oxidation of (R)-sulfolactate by NAD⁺ to sulfopyruvate) and was classified as a sulfolactate dehydrogenase (EC 1.1.1.337) [10, 11]. More recently, the M. jannaschii ComC (MJ1425) has been shown to also catalyze the NAD-dependent oxidation of 2-hydroxy-4-mercaptobutyric acid to 4-mercapto-2-oxobutyric acid, a precursor of the natural product 3-mercaptopropionic acid [14]. 3-mercaptopropionic acid is proposed to act as an alternative coenzyme M. Using a structure-based amino acid alignment, Irimia et al. [10] identified three archaeal homologs to the M. jannaschii ComC (MJ1425) that should act as true sulfolactate dehydrogenases. One, assigned as a malate dehydrogenase from Methanothermus fervidus, has also been characterized (Mfer0389; MdhI or MdhIII or MDH; [8, 11, 15]). The other two are MTH1205 from Methanothermobacter thermautotrophicus and MK0392 from Methanopyrus kandleri. Additionally, MdhI from Methanothermobacter marburgensis (former species name of Methanobacterium thermoautotrophicum strain Marburg; MTBMAc15830) has been characterized and corresponds to MTH1205 [16]. In the alternative pathway, which is predicted to be present in the orders Methanosarcinales and Methanomicrobiales, the steps that lead to sulfopyruvate production are carried out by only two enzymes, the first being cysteate synthase that converts L-phosphoserine to L-cysteate (related to threonine synthase) [2, 6, 17]. An aspartate amino acid transferase then carries out the conversion of cysteate to sulfopyruvate [2], whereby the two coenzyme M pathways converge, with both utilizing ComDE. The final steps in the biosynthesis of coenzyme M, proposed to involve reductive thiolation, have not yet to our knowledge been positively confirmed in any species [1, 7, 18].

Structural analysis of the ComC from the hyperthermophilic methanogen Methanocaldococcus jannaschii (MJ1425) has revealed that the NADH binding pocket is not a typical Rossmann-fold type as found in the malate/lactate dehydrogenase family and represents a new class of NADH-based dehydrogenase [10]. In addition, analyses of potential substrates suggest that electrostatic interactions are important for substrate recognition [10]. In the structure which forms a tight dimer and with both monomers contributing to coenzyme binding, the coenzyme is in an extended conformation with the nicotinamide ring placed to enable transfer of the pro-(S) hydrogen [10]. Only (S)-isomers of substrates are indicated to be converted in the reduction reaction [10].

In this paper, we describe the purification and characterization of ComC from Methanobrevibacter millerae SM9, a representative of the dominant species of rumen methanogens [19–21]. Methane emissions from ruminants are responsible for approximately a quarter of total man-made methane emissions, or 12–17% of total global methane emissions, [22] and are increasingly being recognized for their contribution to climate change. As coenzyme M is essential for methanogenesis, its biosynthesis potentially represents a specific methanogen target for inhibitor development to control methane emissions from ruminants. Sulfopyruvate reduction by ComC is NADH-dependent and therefore readily amenable to development as an absorbance-based screening assay for rapidly screening compound libraries.

2. Materials and Methods

2.1. Materials

Materials including sulfopyruvate, oxaloacetate, α-ketoglutarate, NADPH, NADH, 3-(N-morpholino)propanesulfonic acid (Mops), 1,3-bis(tris(hydroxymethyl)methylamino)propane or Bis-Tris propane (BTP), tris(2-carboxyethyl)phosphine (TCEP), and dithiothreitol (DTT) were purchased from Sigma-Aldrich (USA). Other common chemicals were obtained from ThermoFisher Scientific (NZ).

2.2. Bacteria, Plasmids, and Strains

The ComC gene of Methanobrevibacter millerae SM9 [21] was amplified using forward primer 5-CACCATGAAGATAATGAAGGATAACGAAA and reverse primer 5-TCATTAATCTTCAAGATAAGAATCTATATC with the reverse primer containing two stop codons. The PCR reaction utilized high-fidelity Hercules II DNA polymerase (1.0 μL; Stratagene, USA) in a 50 μL reaction with 0.2 μM of each primer, 0.3 μM dNTP, approximately 20 ng M. millerae SM9 DNA, and 1 × buffer. The PCR cycling parameters had an initial denaturation of 95°C for 2 min, followed by 35 cycles of 94°C for 30 s, 56°C for 30 s, and 68°C for 40 s. The PCR product was purified using agarose DNA electrophoresis and a Wizard SV Gel and PCR kit (Promega, USA). It was then inserted into pET151D using TOPO cloning in chemically competent Escherichia coli strain TOP 10F according to the manufacturer’s instructions (Invitrogen, USA). Colonies were screened by colony PCR using pET151D T7 forward primer and the ComC reverse primer using 2.5 U Taq polymerase (Roche, NZ), and then the recombinant plasmids were isolated using alkaline lysis and purification with a Wizard SV Gel and PCR Clean-up kit (Promega, USA). The plasmid used for expression was sequenced to verify that the gene was in frame and that the sequence was identical to the reference sequence and then transformed into E. coli Rosetta 2 cells (Novagen, USA).

2.3. Protein Expression and Purification

We followed the methods of Schofield et al. [23] to express ComC in E. coli and purify the protein using nickel affinity chromatography, except for the following conditions. The lysis buffer was slightly different (50 mM Tris pH 7.5 containing 1 mM DTT, 300 mM NaCl, 10 mM imidazole, 1% (v/v) Triton X-100, 20% (v/v) glycerol, 2 mM CaCl₂, and 2 mM MgCl₂), and lysis was performed on ice. Cell debris was removed by centrifugation (17,400, 20 min, 4°C), and the supernatant was filtered (0.8 and 0.22 μm). Buffer was exchanged to 20 mM Mops pH 7.0 containing 2 mM TCEP. Glycerol (10% v/v) was added to the purified protein; it was snap frozen in liquid nitrogen and stored at −85°C until further use.

2.4. Assays of ComC Activity

Spectrophotometric measurements and calculation of initial velocity were performed using a Cary 100 UV-vis spectrophotometer (Agilent Technologies, USA) with a thermostatted cuvette holder, using 1 cm path length quartz cuvettes. The consumption of NADH (366 nm, ε 3070 M⁻¹ cm⁻¹) during the reduction of sulfopyruvate by ComC was monitored. Activity was measured at 37°C. M. millerae SM9 was isolated from the rumen of a sheep and has an optimal growth temperature of 38°C [21]. One unit of activity (U) is defined as the conversion of one μmol of NADH to NAD⁺ per minute under standard assay conditions.

Standard assay conditions are 0.08–0.25 μM ComC, 300 μM NADH, 500 μM sulfopyruvate, 400 mM KCl, and 50 mM BTP pH 6.5. The standard assay involved incubation of the above solution without sulfopyruvate at 37°C for 4 min and then initiation of the reaction by the addition of substrate (sulfopyruvate). Kinetics for NADH required variable concentration of the coenzyme (10–500 μM) and 500 μM of sulfopyruvate substrate. Sulfopyruvate kinetics with variable concentration (50–600 μM) were carried out using 400 μM NADH. A concentration of 0.08–0.25 μM ComC in the assay was chosen for accuracy and so as to obtain absorbance changes of about 0.1 to 0.2 per minute. The total volume of all assays was 200 μL. Assays were carried out in triplicate. Kinetic parameters were determined by fitting the data to the Michaelis-Menten equation using GraFit [24].

2.5. Molecular Mass Determination

We followed the methods of Schofield et al. [23] to determine the native molecular mass of ComC using gel filtration chromatography. However, the filtered sample of ComC (400 μl) was at a concentration of 1 mg·ml⁻¹, and the elution buffer was 50 mM Mops pH 7.0 containing 2 mM TCEP and 0.5 M KCl.

2.6. General Methods

We followed the general methods of Schofield et al. [23] for electrophoresis and the determination of protein concentration and pH values of buffers. However, spectrophotometric measurements were performed using a Cary 100 UV-vis spectrophotometer (Agilent Technologies, USA).

3. Results and Discussion

In this study, the M. millerae SM9 ComC was expressed in E. coli, purified using nickel affinity chromatography, and characterized. The apparent molecular mass of ComC was 75 kDa, determined by gel filtration chromatography. As the predicted molecular mass of His-tagged ComC is 40104 Da (376 amino acids), this indicates that ComC is dimeric in solution. Both the M. jannaschii ComC (MJ1425) [11] and the M. fervidus enzyme (Mfer0389) [15] were also shown to be dimeric.

The effect of pH on the activity of sulfopyruvate reduction by ComC, in the presence of 400 mM KCl, was investigated. Two different buffers (pH 4.0–6.5 citrate buffer and pH 6.5–9.5 BTP buffer) were used. Optimal activity was observed at pH 6.5 (Figure 1). For comparison, the pH used for characterizing the M. jannaschii ComC in the reduction reaction was 8.0 [11], the M. fervidus enzyme was assayed at pH 7.4 or 8.0 [8, 11, 15], the M. marburgensis enzyme (MTBMAc15830) was assayed at pH 7.6 [16] and the Chromobacter salexigens ComC was characterized at pH 9.0 [18]. No pH optimization was reported in any of these studies.

Figure 1: Effect of pH on the specific activity of ComC in the reduction reaction. Two buffers were used; ○ = 50 mM citrate buffer and ● = 50 mM BTP. Other assay conditions were standard; 37°C and 400 mM KCl. Assays were carried out in at least triplicate. One unit of activity (U) is defined as the conversion of one μmol of NADH to NAD⁺ per minute under standard assay conditions.

The activity of sulfopyruvate reduction by ComC at pH 6.5 was affected by the KCl concentration with the enzyme showing low activity in the absence of added KCl. The specific activity increased approximately 7-fold with the increasing KCl concentration in the assay and reached a maximum at 400 mM KCl, after which a progressive decrease in activity was observed (Figure 2). No data on the effect of ionic strength or KCl on ComC has yet been reported to our knowledge. However, the M. jannaschii ComC and the M. fervidus enzyme were characterized in 100 mM potassium phosphate [11, 15]. High K⁺ and phosphate concentrations, corresponding to the intracellular medium of M. fervidus cells, enhanced the thermostability of the M. fervidus enzyme [15]. The intracellular concentration of potassium in M. millerae SM9 is unknown, but it is not likely to significantly exceed the ionic strength typical for rumen fluid, which is similar to that of blood (0.15 M) [25].

Figure 2: Effect of KCl concentration on the specific activity of ComC in the reduction reaction. Other assay conditions were standard; 37°C and pH 6.5. Assays were carried out in triplicate. One unit of activity (U) is defined as the conversion of one μmol of NADH to NAD⁺ per minute under standard assay conditions.

Sulfopyruvate, oxaloacetate, and α-ketoglutarate were tested as ComC substrates using the reduction reaction in the presence of NADH. Sulfopyruvate was the preferred substrate, so was used in the standard ComC assay. Oxaloacetate showed 31% of the specific activity seen for sulfopyruvate, and α-ketoglutarate showed approximately 0.2%. A similar substrate preference was seen for the M. jannaschii ComC which showed a 3-fold lower KM for sulfopyruvate compared to oxaloacetate (40 μM and 130 μM, resp.) and a 30-fold higher Vmax/KM [11]. While for α-ketoglutarate, the M. jannaschii ComC had a 50-fold lower KM (40 μM and 1900 μM, resp.) and a 400-fold higher Vmax/KM [11]. Furthermore, the M. fervidus enzyme had a slightly lower KM for sulfopyruvate compared to oxaloacetate and 7-fold higher Vmax/KM [11].

M. millerae SM9 ComC showed less than 1% specific activity in the reduction reaction with NADPH compared to NADH. This bias for NADH was also found for the M. marburgensis enzyme which had negligible specific activity using NADPH compared to NADH in the presence of oxaloacetate (<0.4%) [16]. NADH was also preferred by the M. jannaschii ComC which showed a 5-fold lower KM for sulfopyruvate in the presence of NADH compared to NADPH (40 μM and 210 μM, resp.) and a 60-fold higher Vmax/KM [11]. Additionally, the M. fervidus enzyme had a 3-fold lower KM for sulfopyruvate in the presence of NADH compared to NADPH [11] and an almost 3-fold lower KM for NADH compared to NADPH in the presence of oxaloacetate [15].

The kinetic parameters of the M. millerae SM9 ComC (Table 1) were obtained by plots of substrate concentration versus specific activity and their respective double-reciprocal plots (Figure 3). The rates followed typical Michaelis-Menten kinetics for both sulfopyruvate (KM 196 μM) and NADH (KM 55.1 μM) with the sulfopyruvate value being 5-fold that found for M. jannaschii ComC (40 μM) [11]. The M. millerae SM9 ComC KM for sulfopyruvate is also higher than that for the M. fervidus enzyme (70 μM) [11]. The Vmax of the M. millerae SM9 ComC for sulfopyruvate was 93.9 U·mg⁻¹, and for NADH, it was 73.1 U·mg⁻¹. The Vmax value is lower than that for sulfopyruvate of the M. jannaschii ComC (370 U·mg⁻¹) [11] and similar to that for sulfopyruvate of the M. fervidus enzyme (120 U·mg⁻¹) [11]. The M. jannaschii ComC and the M. fervidus enzymes show substrate inhibition at very low sulfopyruvate concentrations (100 μM, i.e., 2.5-fold or 1.4-fold the sulfopyruvate KM values of these enzymes, resp.) [11]. The results in Figure 3 may indicate a similar effect for the M. millerae SM9 ComC, with possible substrate inhibition at 600 μM, that is, 3-fold the sulfopyruvate KM value.

Table 1: Kinetic parameters for ComC, for sulfopyruvate, and for NADH. Standard assay conditions in the reduction reaction were used; 400 mM KCl, pH 6.5, and 37°C. Assays were carried out in triplicate. One unit of activity (U) is defined as the conversion of one μmol of NADH to NAD⁺ per minute under standard assay conditions.

KM (μM)
Vmax (U·mg⁻¹)
kcat (s⁻¹)

Sulfopyruvate
196 ± 71
93.9 ± 14.5
62.8 ± 9.7

NADH
55.1 ± 7.6
73.1 ± 3.2
48.9 ± 2.2

Figure 3: Michaelis-Menten plots for ComC. Michaelis-Menten plots were used to calculate kinetic parameters. The respective Lineweaver-Burk plots are inset. Standard assay conditions in the reduction reaction were used; 400 mM KCl, pH 6.5, and 37°C. Assays were carried out in triplicate. One unit of activity (U) is defined as the conversion of one μmol of NADH to NAD⁺ per minute under standard assay conditions.

Methanobrevibacter is indicated to be the dominant genus of methanogens in the rumen of sheep and cows [19, 20]. Coenzyme M is required for methanogenesis to occur, and therefore targeting enzymes for coenzyme M synthesis represents a valid approach for finding novel small molecule compounds for controlling ruminant methane emissions [20]. Analysis of the genome of coenzyme M-requiring Methanobrevibacter ruminantium strain DSM 1093 has revealed, interestingly, that this rumen methanogen does not contain an intact coenzyme M biosynthesis pathway [26]. It contains genes for ComB and ComC but no recognizable genes for ComA or ComDE [26]. Thus, M. ruminantium is likely to use coenzyme M derived from other methanogens within the rumen. Lovley et al. [27] have described the isolation of coenzyme-requiring rumen methanogens and those that do not, with the methanogens that can synthesize coenzyme M apparently having faster growth rates.

4. Conclusions

We have reported here the cloning, expression, and biochemical characterization of ComC from a representative rumen methanogen. The activity of M. millerae SM9 ComC was assayed in the reverse reaction of coenzyme M biosynthesis using sulfopyruvate. The only other ComC to have been extensively characterized is from the thermophilic methanogen M. jannaschii [8, 11]. Due to its role in two separate cofactor biosynthesis pathways and its dependence on NADH, the M. millerae SM9 ComC represents a potential means of screening compounds for their ability to inhibit methanogens and in so doing, help mitigate methane emissions from ruminants.

Conflicts of Interest

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank the Pastoral Greenhouse Gas Research Consortium for funding and Mark Aspin for support. They thank Bill Kelly and Graeme Attwood for supplying the gene sequence for the M. millerae SM9 ComC gene prior to publishing the genome.

References

D. E. Graham and R. H. White, “Elucidation of methanogenic coenzyme biosyntheses: from spectroscopy to genomics,” Natural Product Reports, vol. 19, no. 2, pp. 133–147, 2002.

D. E. Graham, “2-oxoacid metabolism in methanogenic CoM and CoB biosynthesis,” Methods in Enzymology, vol. 494, pp. 301–326, 2011.

W. E. Balch and R. S. Wolfe, “Transport of coenzyme M (2-mercaptoethanesulfonic acid) in Methanobacterium ruminantium,” Journal of Bacteriology, vol. 137, no. 1, pp. 264–273, 1979.

C. E. Hand and J. F. Honek, “Biological chemistry of naturally occurring thiols of microbial and marine origin,” Journal of Natural Products, vol. 68, no. 2, pp. 293–308, 2005.

N. Santoro and J. Konisky, “Characterization of bromoethanesulfonate-resistant mutants of Methanococcus voltae: evidence of a coenzyme M transport system,” Journal of Bacteriology, vol. 169, no. 2, pp. 660–665, 1987.

Y. Liu, L. L. Beer, and W. B. Whitman, “Sulfur metabolism in archaea reveals novel processes,” Environmental Microbiology, vol. 14, no. 10, pp. 2632–2644, 2012.

D. E. Graham, M. Graupner, H. Xu, and R. H. White, “Identification of coenzyme M biosynthetic 2-phosphosulfolactate phosphatase: a member of a new class of Mg²⁺-dependent acid phosphatases,” European Journal of Biochemistry, vol. 268, no. 19, pp. 5176–5188, 2001.

M. Graupner and R. H. White, “The first examples of (S)-2-hydroxyacid dehydrogenases catalyzing the transfer of the pro-4S hydrogen of NADH are found in the archaea,” Biochimica et Biophysica Acta (BBA) - Protein Structure and Molecular Enzymology, vol. 1548, no. 1, pp. 169–173, 2001.

M. Graupner, H. Xu, and R. H. White, “Identification of the gene encoding sulfopyruvate decarboxylase, an enzyme involved in biosynthesis of coenzyme M,” Journal of Bacteriology, vol. 182, no. 17, pp. 4862–4867, 2000.

A. Irimia, D. Madern, G. Zaccaï, and F. M. Vellieux, “Methanoarchaeal sulfolactate dehydrogenase: prototype of a new family of NADH-dependent enzymes,” The EMBO Journal, vol. 23, no. 6, pp. 1234–1244, 2004.

M. Graupner, H. Xu, and R. H. White, “Identification of an archaeal 2-hydroxy acid dehydrogenase catalyzing reactions involved in coenzyme biosynthesis in methanoarchaea,” Journal of Bacteriology, vol. 182, no. 13, pp. 3688–3692, 2000.

H. Muramatsu, H. Mihara, M. Goto et al., “A new family of NAD(P)H-dependent oxidoreductases distinct from conventional Rossmann-fold proteins,” Journal of Bioscience and Bioengineering, vol. 99, no. 6, pp. 541–547, 2005.

C. R. Goward and D. J. Nicholls, “Malate dehydrogenase: a model for structure, evolution, and catalysis,” Protein Science, vol. 3, no. 10, pp. 1883–1888, 1994.

K. D. Allen and R. H. White, “Occurrence and biosynthesis of 3-mercaptopropionic acid in Methanocaldococcus jannaschii,” FEMS Microbiology Letters, vol. 363, no. 19, 2016.

E. Honka, S. Fabry, T. Niermann, P. Palm, and R. Hensel, “Properties and primary structure of the L-malate dehydrogenase from the extremely thermophilic archaebacterium Methanothermus fervidus,” European Journal of Biochemistry, vol. 188, no. 3, pp. 623–632, 1990.

H. Thompson, A. Tersteegen, R. K. Thauer, and R. Hedderich, “Two malate dehydrogenases in Methanobacterium thermoautotrophicum,” Archives of Microbiology, vol. 170, no. 1, pp. 38–42, 1998.

D. E. Graham, S. M. Taylor, R. Z. Wolf, and S. C. Namboori, “Convergent evolution of coenzyme M biosynthesis in the Methanosarcinales: cysteate synthase evolved from an ancestral threonine synthase,” Biochemical Journal, vol. 424, no. 3, pp. 467–478, 2009.

K. Denger and A. M. Cook, “Racemase activity effected by two dehydrogenases in sulfolactate degradation by Chromohalobacter salexigens: purification of (S)-sulfolactate dehydrogenase,” Microbiology, vol. 156, no. 3, pp. 967–974, 2010.

G. Henderson, F. Cox, S. Kittelmann et al., “Effect of DNA extraction methods and sampling techniques on the apparent structure of cow and sheep rumen microbial communities,” PLoS One, vol. 8, no. 9, 2013.

T. Knight, R. S. Ronimus, D. Dey et al., “Chloroform decreases rumen methanogenesis and methanogen populations without altering rumen function in cattle,” Animal Feed Science and Technology, vol. 166-167, pp. 101–112, 2011.

W. J. Kelly, D. M. Pacheco, D. Li, G. T. Attwood, E. Altermann, and S. C. Leahy, “The complete genome sequence of the rumen methanogen Methanobrevibacter millerae SM9,” Standards in Genomic Sciences, vol. 11, no. 1, 2016.

K. R. Lassey, “Livestock methane emission and its perspective in the global methane cycle,” Australian Journal of Experimental Agriculture, vol. 48, no. 2, pp. 114–118, 2008.

L. R. Schofield, A. K. Beattie, C. M. Tootill, D. Dey, and R. S. Ronimus, “Biochemical characterisation of phage pseudomurein endoisopeptidases PeiW and PeiP using synthetic peptides,” Archaea, vol. 2015, Article ID 828693, 12 pages, 2015.

R. J. Leatherbarrow, GraFit Version 7, Erithacus Software Ltd., Horley, UK, 2009.

R. A. Kohn and T. F. Dunlap, “Calculation of the buffering capacity of bicarbonate in the rumen and in vitro,” Journal of Animal Science, vol. 76, no. 6, pp. 1702–1709, 1998.

S. C. Leahy, W. J. Kelly, E. Altermann et al., “The genome sequence of the rumen methanogen Methanobrevibacter ruminantium reveals new possibilities for controlling ruminant methane emissions,” PLoS One, vol. 5, no. 1, 2010.

D. R. Lovley, R. C. Greening, and J. G. Ferry, “Rapidly growing rumen methanogenic organism that synthesizes coenzyme M and has a high affinity for formate,” Applied and Environmental Microbiology, vol. 48, no. 1, pp. 81–87, 1984.