Structure-Functional Examination of Cysteine Synthase A (CysK) from Limosilactobacillus reuteri LR1

TL;DR

This study explores the structure and function of CysK from Limosilactobacillus reuteri, revealing its biochemical properties and potential for biotechnology.

Contribution

The study provides new insights into the thermal stability and structural features of LreCysK, highlighting its biotechnological potential.

Findings

LreCysK exhibits similar kinetic parameters to CysKs from E. coli and Arabidopsis thaliana.

The enzyme's melting point increases to nearly 90°C with PLP, showing high stability.

LreCysK is a dimer with a unique active site featuring residues like Ile50.

Abstract

This study presents a comprehensive analysis of cysteine synthase A (CysK) from Limosilactobacillus reuteri LR1 (LreCysK), an enzyme involved in the biosynthesis of L-cysteine. This protein supports crucial cellular functions such as sulfur metabolism, antioxidant defense, detoxification, and protein synthesis. Previously, the gene encoding LreCysK was cloned, and the enzyme with His-tag on the N-terminus was obtained in active and soluble form. Here, kinetic parameters of the enzyme were determined by the previously developed high-pressure liquid chromatography (HPLC) and ninhydrin methods. It was found that LreCysK has similar KMOAS and kcat as CysKs from Escherichia coli and from the model plant Arabidopsis thaliana. The thermal stability of LreCysK was studied using differential scanning calorimetry. It was revealed that the melting point of the enzyme increases to almost 90°C when Pyridoxal-5 phosphate (PLP) is added, indicating that the stability of the enzyme complex with PLP is relatively high. Structural studies revealed that LreCysK is a dimer, and its active site is similar to those of other enzymes, but exhibits some features characteristic of lactobacilli CysKs (GISA), as well as unique residues, such as Ile50. Also, the potential biotechnological applications of LreCysK are discussed. These findings enhance our understanding of LreCysK’s biochemical versatility and its potential applications in biotechnology and medicine.

Figures

Full text

1. Introduction

Cysteine synthase A (CysK, EC 2.5.1.47) is a pyridoxal phosphate (PLP)-dependent enzyme that belongs to the transferase class and catalyzes the formation of L-cysteine from sulfide and O-acetyl-L-serine (OAS). This reaction represents a critical step in the sulfur assimilation pathway. This has a vital role in maintaining sulfur homeostasis and supporting cellular functions such as antioxidant defense, detoxification, and protein synthesis [1]. The significance of cysteine biosynthesis extends beyond individual metabolic pathways; it is a cornerstone in the synthesis of glutathione, an essential antioxidant that protects cells from oxidative stress and maintains redox balance. In the context of bacterial pathogens, cysteine biosynthesis can confer a survival advantage, particularly under conditions of oxidative stress or nutrient limitation. Therefore, CysK inhibitors are being actively studied as potential antimicrobial agents [2,3]. Additionally, CysK plays a role in contact-dependent growth inhibition, L-isoleucine production, and biofilm formation, which may suggest different possible applications of this enzyme [3,4].

L. reuteri is the probiotic strain, and its properties are well studied and include the control of insulin sensitivity, anti-inflammatory effects, synthesis of vitamins, as well as low molecular weight compounds with antibacterial effect [5,6]. The study of CysK from L. reuteri (LreCysK) is of particular interest due to its potential role in conferring resistance to environmental stressors and enhancing the probiotic effects of the bacterium, such as modulating the host immune response and inhibiting pathogenic microorganisms [2,7]. Understanding the biochemical properties of LreCysK, including its kinetics, thermostability, and structure, could provide valuable insights into its role in L. reuteri’s probiotic functions, especially in its interaction with pathogenic bacteria and potential applications in biotechnology and medicine. Due to its crucial physiological role in cysteine synthesis in various pathogenic bacteria, specific inhibitors of CysK can exhibit an antibacterial effect [3,4]. L. reuteri is a component of normal gut microbiota. Consequently, structural studies of CysK from this bacterium may help in developing specific inhibitors of this enzyme from pathogenic bacteria and in preventing the inhibition of normal microbiota growth.

One of the most challenging problems in modern medicine is the emergence of pathogenic bacteria with multiple drug resistance. That is why many researchers try to find new sources of antibacterial substances. Recent studies have shown that L. reuteri LR1 can inhibit the growth of Klebsiella pneumoniae, a pathogenic bacterium known for its multidrug resistance [8]. Exoproteome analysis demonstrated that L. reuteri responds to K. pneumoniae by secreting new enzymes of various groups, including peptidases, nucleosidases, and metabolic enzymes. One of the proteins found was metabolic enzyme cysteine synthase A (CysK), which may contribute to its antagonistic effect [8]. Additionally, recent research showed the possible presence of such activity [9]. This discovery suggests a potential mechanism through which L. reuteri could exert a competitive advantage over pathogenic bacteria by modulating its metabolic pathways. Since K. pneumoniae is a pathogenic organism that causes sepsis, pneumonia, inflammation of the urinary system, and problems with the liver and kidneys, new approaches to combat it are needed, especially since resistance to classic antibiotics is increasing every year [10].

The literature contains data on the structural organization of CysK from several bacteria (6Z4N, 6Z4N, 2Q3C, 5XOQ, 4I1Y, 2EGU) and protozoa (4JBN, 3SPX, 8B9Y). It is known that for most structures, CysK is a dimer consisting of two domains, located closer to the N- and C-termini of each subunit. There is extremely little data on the structure of CysK from lactobacilli in the scientific literature. PLP binds to the amino acid residue of lysine by a covalent bond similar to other PLP-dependent enzymes to form a Schiff base [11]. Previously, the gene encoding cysteine synthase A from L. reuteri was cloned, and we obtained the enzyme in an active and soluble form [12]. Additionally, a method for determining activity using Hydrophilic Interaction Liquid Chromatography (HILIC) chromatography was developed [12]. In this paper, we present data on the biochemical and structural characteristics of this enzyme.

2. Results and Discussion

2.1. Obtaining the LreCysK Enzyme with His–Tag

CysK with a His–tag is expressed in E. coli cells in a soluble form, with a yield of approximately 40 mg per liter of nutrient medium. The degree of purity of the enzyme was confirmed using SDS–PAGE (Supplementary Figure S1). Tandem MALDI-TOF/TOF mass spectroscopy confirmed the sequence of LreCysK. The percentage after MALDI analysis, using the theoretical sequence, was 92% (Supplementary Figure S2). The enzyme was obtained in a highly purified and soluble state with a degree of purity of at least 95%. The theoretical mass of the LreCysK monomer with His–tag is 33.12 kDa, which is also consistent with the SDS–PAGE.

2.2. Determination of the Oligomeric Composition of the Enzyme

The oligomeric composition of CysK was determined using analytical gel filtration. The calibration curve with LreCysK plotted on it is presented in Figure 1.

The theoretical monomer mass of LreCysK with His-tag at N-terminus, calculated using the amino acid sequence, is 33.12 kDa. Parameter K_av_ was calculated using the formula K_av_ = (V_e_ − V_O_)/(V_C_ − V_O_), where V_O_ is the column void volume, V_e_ is the elution volume, and V_C_ is the geometric column volume. For all proteins used for the calibration curve, the parameter K_av_ was calculated (Table S1, Supplementary Figure S3). The M_r_ for CysK was calculated from the dependence of K_av_ on lgM_r_.

The molecular mass of LreCysK, as determined by gel filtration, was 54 kDa, corresponding to a dimer. This correlates with the structures of different CysK enzymes, which are typically dimers.

2.3. Kinetic Parameters of LreCysK

The widely used method for studying CysK kinetics is the spectrophotometric method, which utilizes the ninhydrin acidic reagent (NAR). In this method, cysteine binds with ninhydrin to form a pink product, measured at 560 nm. This method has several disadvantages, including low stability of the colored reaction product (color saturation decreases over time), limited specificity for cysteine (ninhydrin can interact with other amino acids), and a relatively high detection limit (approximately 0.5–1 mM). In this study, we employed the ninhydrin method and developed a novel approach using HILIC for the separation of cysteine from other components of the mixture and its subsequent determination [12]. Using HPLC generally provides better results for this task due to its ability to separate the target analyte from other substances. HILIC is a relatively new method for measuring enzymatic activity, providing better and faster separation of polar substances. To separate a mixture of complex composition: two substrates (OAS and sulfide), a potential product (cysteine), and components of a buffer solution (HEPES), the following composition of the mobile phase was chosen: 12 mM phosphate buffer, pH 6.5 (aqueous) and acetonitrile (22:78 vol.%) according to [12]. It is essential to be able to separate cysteine and HEPES, as well as to separate them from the two substrates.

The kinetic parameters are presented in Table 1. This table also shows the kinetic parameters of other CysK enzymes. Parameters for the enzymes obtained in this work are in bold. Supplementary Figure S4 shows the dependence of the enzymatic reaction rate on the substrate concentration used to calculate the kinetic parameters for the LreCysK using HILIC and ninhydrin methods.

Results demonstrate that both methods yield similar K_M_ values for OAS; however, the k_cat_ measured by the ninhydrin method is 2.5 times higher than that measured by our new method. This may be due to the nonselective reaction between NAR and other reaction components, except cysteine. When comparing the kinetic parameters of the LreCysK enzyme obtained in this work with those from other sources, it is notable that the Michaelis constant for OAS is typically 10 to 100 times higher than that for sulfide. The only exception is the enzyme from L. casei FAM18110. Interestingly, although both enzymes are derived from lactobacteria, their characteristics differ. LreCysK is closer to CysK from E. coli, S. typhimorium, and A. thaliana, which allows us to assume their common evolutionary origin. AthCysK is indeed evolutionarily close to LreCysK (Supplementary Figure S5), but EcoCysK is an evolutionarily distant enzyme. The similarity in the kinetic parameters of EcoCysK and LreCysK can be attributed to their shared habitat—the human gastrointestinal tract.

2.4. Temperature Stability of LreCysK

Differential scanning calorimetry (DSC) was used to study the temperature stability of LreCysK. During the experiment, the temperature of the sample cell and the comparison cell was linearly increased, and the change in heat capacity was monitored. Since protein denaturation is a phase transition, the experimental curve shows an increase in the change in heat capacity of the sample with a peak at the temperature at which the maximum rate of denaturation occurs, and which is used as a characteristic of the protein’s thermal stability (the so-called phase transition temperature Tm). DSC data for both variants of LreCysK are presented in Figure 2 and Table 2.

DSC analysis revealed two well-distinguishable calorimetric domains with melting points of 72.3 °C and 83.5 °C in the apo form of CysK, as well as two poorly distinguishable calorimetric domains in the holo form (in the presence of PLP) with melting points of 80.1 and 88.9 °C, respectively. A cardinal shift in melting point and redistribution of denaturation enthalpy between domains indicated that the holo form of CysK is more stable than the apo form. This effect is often observed in PLP-dependent enzymes and has been demonstrated for serine dehydratase, ornithine δ-aminotransferase, and tyrosine aminotransferase [22,23]. The observed thermal stability of LreCysK is likely due to a combination of factors, including a highly hydrophobic core, extensive hydrogen bonding, and ionic interactions. The presence of proline residues in strategic locations could also contribute to the enzyme’s rigidity and resistance to thermal denaturation. Further studies involving mutagenesis could help identify specific residues or regions responsible for this thermostability, providing insights for engineering more stable enzyme variants. The high thermal stability observed is a desirable trait for industrial applications, where enzymes are often exposed to high temperatures.

2.5. Structural Analysis of LreCysK

The crystal structure of the apoenzyme (LreCysKapo) was solved at a resolution of 2.2 Å. The asymmetric unit contained two subunits related by 2-fold non-crystallographic symmetry (NCS). (RMSD between Cα atoms between subunits is 0.3 Å) organized in a dimer whose buried area comprised 16% of the total surface area of each subunit (Figure 3). According to gel filtration, LreCysK is also a dimer in solution. The subunit of LreCysKapo consists of large (residues 21–43 and 161–313) and small (residues 1–20 and 44–152) domains, similar to those of other CysK proteins described in [24,25]. We assume that these domains are the N- and C-terminus domains, which often occur in the structure of cysteine synthases [26]. The full structure is also typical for cysteine synthases: every domain consists of α/β-folds with the central β-sheet surrounded by α-spirals [27]. The closest structural homologs of the LreCysKapo are cysteine synthase from Mycobacterium ulcerans (MulCysKapo, PDB ID 4I1Y, RMSD between Cα atoms of subunits is 2.03 Å^2^, Table 3) and O-acetylserine sulfhydrylase from Geobacillus kaustophilus HTA426 (GkaCysKapo, PDB ID 2EGU, RMSD is 2.33 Å). Structure superposition (Figure 3B, C) and RMSD showed severe differences between LreCysKapo and MulCysKapo, with GkaCysKapo, but the RMSD of Cα atoms between subunits of MulCysKapo and GkaCysKapo is 1.3 Å (95% of the aligned residues). To elucidate the reasons for this difference, subunits were superimposed by domains. RMSD between large and small domains of LreCysKapo and MulCysKapo are 0.78 Å (97% of aligned residues) and 1.17 Å (98% of aligned residues), respectively (Figure 3D), indicating that the large and small domains of LreCysKapo are similar to those from MulCysKapo. The high RMSD between subunits was caused by differences in the relative orientation of the large and small domains in LreCysKapo compared to the reference structures.

Moreover, α-helix (residues 162–172) (Figure 4) in LreCysKapo has a different conformation than in MulCysKapo (residues 152–166) and GkaCysKapo (residues 157–166) (Figure 4). In LreCysKapo, residues of this α-helix (E170 and D174) form polar contacts with N-terminus His–tag residues (H6 and H7) from the adjacent subunit (Figure 4A). The addition of His–tag to the N-terminus of the enzyme could change the conformation of this α-helix. We assume that these changes don’t influence the catalytic activity since this α-helix is not part of the active site and presumably only stabilizes the conformation. Residues 214–240 and 213–235 in subunits A and B, respectively, have no electron density. Moreover, this region is absent in MulCysKapo (residues 183–204) and GkaCysKapo (residues 207–229). This region has no electron density in all apo form structures due to the high mobility. In StyCysK holo form structure (PDB ID 6Z4N, StyCysKholo), this region (residues 206–240, black in Figure 4B) is clearly visible and shields the active site from the solvent and fixes the PLP in the active site.

Table 3 shows the comparison of the LreCysKapo structure with the known crystal structures of CysK from other organisms. All mentioned structures have one molecule per asymmetric unit. We can indicate that, except for GkaCysK and AthCysK, practically all the mentioned structures are dimers.

Since the crystal structure of our LreCysKapo protein contained two unresolved loops (positions 153–160 and 214–240), we modelled these loops. Our further structural investigations of LreCysKapo were made on this model structure. In Figure 5, the active site of LreCysKapo, featuring highly conserved motifs, is presented based on this model structure. Analysis of amino acid sequences of well-studied CysKs was already done by us and published in [28]. In this paper, we focused on certain sequences of enzymes, closely related to LreCysK, based on their biochemical characteristics (Table 2) and phylogenetic analysis (Figure S5). The alignment of the similar CysKs based on these factors is presented in Figure 6.

The conserved lysine residue (Lys51 in LreCysKapo) plays a pivotal role in forming a Schiff base with the PLP cofactor, a highly conserved feature among PLP-dependent enzymes. This interaction is crucial for catalytic activity, as it facilitates the nucleophilic attack on the OAS substrate. The pyrimidine ring of PLP fits the small hydrophobic cavity of the active site and is usually stabilized by a Val residue, although in the case of LreCysK, Ile50 plays this role. Glycine-rich region ^185^GTGGT^189^ creates a loop near the phosphate group of the PLP, making hydrogen bonds with the oxygen atoms. The deprotonated 3′-OH group of the PLP creates a hydrogen bond with the amide nitrogen of the Asn81, and the nitrogen of the pyrimidine ring of PLP creates a hydrogen bond with the Ser273, thus being the last point of the coenzyme stabilization. The catalytic site in the case of LreCysK consists of amide atoms of a highly conserved region ^77^TSGNT^82^, which correlates with sequences of practically all known CysK structures. The other conserved region crucial for catalysis is ^229^GISA^232^ (which is usually GIGA, as can be seen from other structures). We assumed that GISA is a structural characteristic of lactobacilli CysKs as both LreCysK and LcaCysK possess it. This assumption is supported by research on amino acid sequences from different lactobacilli. Most of them have the GISA, not the GIGA (Supplementary Figure S6). During catalysis, when the intermediate product is formed, these two loop regions close the active site, limiting solvent access to the highly hydrophobic region of the active site and thereby protecting the intermediate product.

Amino acid alignment of LreCysK with several other CysK enzymes showed a high number of aligned regions with over 90% of identity (Table 4, Figure 6). For alignment, we selected the enzyme from another Lactobacillus species (L. casei) and two other enzymes that share similar kinetic properties with LreCysK (Table 2). Across the alignment, several highly conserved regions are observed, particularly in key functional motifs related to substrate binding and catalysis as discussed earlier. These regions are represented by loops that have identical or similar residues across all organisms, indicating their importance in maintaining CysKs enzymatic function. Curiously, the two lactobacterial enzymes LreCysK and LcaCysK (the latter being from L. casei FAM18110) have a difference in one catalytically important region. This region is ^347^GISA^350^ (numbered as in alignment), and it closes the catalytic cavity during the reaction as discussed earlier. Other CysK enzymes have GIGA in this region instead. Apparently, the substitution of glycine for serine in the GISA loop, characteristic of lactobacilli, does not significantly affect the Michaelis constant for o-acetylserine. Since the K_M_^OAS^ in LreCysK is similar to the corresponding values for CysK from A. thaliana and E. coli, which have GIGA in this region. This allows us to assume that LcaCysK has some other important amino acid residues that affect the binding of o-acetylserine.

3. Materials and Methods

3.1. Obtaining LreCysK and Protein Expression

The gene encoding CysK enzyme was cloned previously from the L. reuteri LR1 bacteria (the strain was kindly provided by the All–Russian Dairy Research Institute (VNIMI)) with the addition of a fragment encoding a six-histidine residue sequence (His–tag) to the N-terminus of the enzyme for acceleration and simplifying the enzyme purification process [12].

Sequencing results showed that the genetic construct contains only the CysK gene with the target insert (His–tag at the N-terminus of the enzyme). The size of the CysK gene with a His–tag was 939 bp, the gene encodes a protein with a length of 313 amino acid residues and a size of 33.12 kDa, the accession number of the CysK enzyme from L. reuteri without His–tag is MBU5982312.1.

The plasmid, containing the gene encoding LreCysK, was used to transform E. coli BL21(DE3) CodonPlus/pLysS cells resistant to chloramphenicol. Expression of the enzyme was performed according to the algorithm from previous work [29]. The resulting cell suspensions were stored at –20 °C until purification.

3.2. Enzyme Purification

E. coli cells containing the CysK enzyme were lysed on an ultrasonic disintegrator. The cell debris was settled on an Eppendorf 5804 R centrifuge at 4 °C and 5000 rpm, and the supernatant was transferred into clean plastic tubes with a volume of 50 mL. Purification was carried out on an AKTA Start chromatography system using a HisTrap FF 1 mL column (Cytiva, Marlborough, MA, USA). The column was pre-equilibrated with buffer A (0.05 M Tris–HCl, 0.5 M NaCl, 0.02 M imidazole, pH 7.5). The supernatant after cell disruption in buffer A was applied to the column at a flow of 0.5 mL min*^−^*^1^. After the disappearance of the absorption peak at 280 nm of impurity proteins that did not bind to the column, the target enzyme was eluted with a linear gradient of increasing concentration of buffer B (0.05 M Tris–HCl, 0.5 M NaCl, 0.5 M imidazole, pH 7.5). The resulting enzyme solution was desalted using size-exclusion chromatography on a Sephadex G25 (Cytiva, Marlborough, MA, USA) column into a solution containing 0.05 M Tris-HCl, pH 7.5. The purity of the resulting enzyme solution was confirmed by SDS–PAGE as described in [30].

Confirmation of the amino acid sequence of the CysK was performed using MALDI-TOF/TOF mass spectrometry as described previously [29].

3.3. Determination of Oligomeric Composition of the Enzyme by Size-Exclusion Chromatography

The oligomeric composition of the enzyme was determined using analytical size-exclusion chromatography (gel filtration) on an AKTA Start chromatograph on a HiLoad 16/600 Superdex 200 pg column (Cytiva, Marlborough, MA, USA). To construct a calibration curve, a set of molecular weight standards for gel filtration and analytical electrophoresis, Gel Filtration Calibration Kit HMW (Cytiva, Marlborough, MA, USA) was used. Column equilibration and calibration were performed according to the manufacturer’s protocol. The test sample was added to the column in an amount of 1 mg, and its retention time on the column was observed. Next, using a previously constructed calibration curve, the molecular mass of the sample was determined and compared with the theoretical mass of the monomer, from which a conclusion was drawn about the oligomeric composition of the enzyme.

3.4. Carrying out an Enzymatic Reaction

The enzymatic reaction was carried out using 99% pure O-acetyl-L-serine (OAS), sodium sulfide nonahydrate from Sigma–Aldrich (St. Louis, MO, USA), PLP monohydrate from Loba-Chemie Austranal-Praparate (Fischamend, Austria), and water purified by a MilliQ unit. The reaction was carried out as follows: the desired stock concentrations of PLP, OAS and sulfide were prepared by weighing and dissolving reagents in 10 mM HEPES buffer pH 7.5, then both substrates were added at the required concentrations (diluted with HEPES buffer when needed) to plastic test tubes (total volume 1.5 mL) to the 1000 µL volume, PLP was added to the final concentration 5 µM. For each concentration of substrate, three replicates were made. Next, 10 µL of purified CysK solution at a concentration of approximately 290 µg mL*^−^^1^ (for HILIC) and 77 µg mL^−^^1^ (for ninhydrin method) was added to each sample (the final concentration in the solution is 2,9 (for HILIC) and 0.77 (for ninhydrin method) µg mL^−^*^1^) and stirred. At certain time intervals, the enzymatic reaction was stopped by adding 10 µL of concentrated HCl at different times of the reaction for conducting analysis with HILIC, as previously developed [9], or the ninhydrin method. The prepared samples were then analyzed. For measuring the kinetics of the enzyme, one substrate was added at a fixed, saturating concentration to the mixture, while the concentration of the other was varied.

3.5. Conducting the Analysis Using Nynhidrin Reaction

For analysis, 300 µL of modified Ninhydrin Acid Reagent (NAR, comprising 250 mg of ninhydrin, 11 mL of concentrated acetic acid, and 4 mL of concentrated HCl) was added to the reaction mixture, resulting in a final volume of 1.02 mL. The mixtures were heated to 95 °C for 10 minutes. After this, mixtures were cooled on ice for 1–2 min, and 650 µL of ethanol was added to them. The presence of cysteine in mixtures was gained by absorption measurement at a wavelength of 560 nm. Cysteine concentrations were calculated using a calibration curve that was linear in the range of 0.1 mM to 0.7 mM.

3.6. Conducting the Analysis Using HILIC

For analysis, 99% pure cysteine from Sigma–Aldrich (St. Louis, MO, USA), water purified using a MilliQ unit, and HPLC gradient grade acetonitrile (Panreac, Darmstadt, Germany) were used. The experiments were conducted using a Dionex UltiMate 3000 Liquid Chromatography system (Dionex, a part of Thermo Scientific, Sunnyvale, CA, USA), equipped with a gradient pump, autosampler, column thermostat, and diode array detector. Chromatograms were recorded using a personal computer and the Chromeleon 7.0 software package (Thermo Fisher Scientific, Waltham, MA, USA). A 100 × 3.0 mm i.d. stainless steel column was packed with a hydrophilic stationary phase based on 3-aminopropylsilica with amide groups in a functional layer, which was prepared via multicomponent Ugi reaction according to [12,31,32] by means of 2-acetylpyrrole, 2-morpholinoethyl isocyanide, and glycolic acid. This column exhibited high selectivity towards amino acids [31], and a column with a similar functional layer [32] demonstrated high stability during enzymatic activity determination. The separation and analysis were carried out as described in [12]. A set of “blank” injections of a substrate mixture solution in HEPES and HEPES buffer solution alone was also conducted to check for corresponding impurities. Kinetic parameters of the reaction were determined using the recently described HILIC method, and calculations were carried out in Origin 2018.

3.7. Kinetic Parameters Determination for LreCysK

The dependence of the enzymatic reaction rate on the substrate concentration was obtained as described above. This dependence was analyzed using Origin Pro 8.5. The Michaelis constant KM was determined using nonlinear regression in the range of substrate concentrations of 0.5 KM–5 KM. From the same dependence, the maximum speed of the enzymatic reaction, Vmax, was determined. The concentration of the purified enzyme was assayed by measurement of solution absorption at 280 nm (A^280^) and calculated using the equation C_protein_ = A^280^/(ε_protein_ × l), where ε_protein_ was 11,920 M^−1^·cm^−1^. The A^280^ value was corrected by absorption of PLP at the same wavelength. Using Vmax and Cprotein, the rate constant of the enzymatic reaction, k_cat_, was calculated from the equation Vmax = k_cat_ × Cprotein.

3.8. Study of Enzyme Thermostability

To study the thermal stability of the LreCysK enzyme, differential scanning calorimetry (DSC) was used. DSC was carried out using a MicroCal VP-Capillary differential scanning microcalorimeter (Malvern Instruments, Northampton, MA, USA) equipped with 140 µL tantalum capillary cells. The experimental protocol involved a heating rate of 1 °C·min^−1^, with an overpressure maintained at 3 bars to prevent boiling inside the cells. Enzyme samples were prepared at concentrations ranging from 1 to 2 mg·mL^-1^ in 0.1 M Na-phosphate buffer (pH 7.0). For samples containing PLP, the cofactor was added at a concentration of 1.0 mM. The instrumental baseline was established by recording buffer-to-buffer scans, which were subsequently subtracted from the protein data to correct for background signals. The resulting calorimetric peaks corresponding to individual domains were deconvoluted using the Zubov-Markov method [33], implemented through an original procedure developed for the MATLAB R2025b computing environment (MathWorks, Natick, MA, USA).

3.9. Multiple Sequence Alignment (MSA)

For MSA, Clustal-Omega was selected using the program Clustalo (version 1.2.4) https://www.ebi.ac.uk/jdispatcher/msa/clustalo (accessed on 7 September 2025). Visualization of MSA was performed using the Jalview 2.11.5.0 program. For visualization, Clustal colors were chosen in Jalview 2.11.5.0. [software] https://www.jalview.org/.

3.10. Crystallization and X-Ray Diffraction Data Collection

Initial crystallization screening was done by the “sitting drop” vapor diffusion method on a robotic system (Oryx4, East Garston, Berkshire, UK) with 96-well VDX plates (Hampton Research, Aliso Viejo, CA, USA) and commercial crystallization screens from Hampton Research (Aliso Viejo, CA, USA) and Molecular Dimensions Inc. (Holland, OH, USA). The N-terminal His-tagged LreCysK apo form was concentrated to 15 mg·mL^-1^ in 0.05 M Tris-HCl buffer, pH 7,5. The protein and crystallization solution were mixed in the ratios 2:1, 1:2 (using a 0.3 µL drop volume), and 1:1 (using a 0.2 µL drop volume). The reservoir contained 50 µL of the precipitant solution. The initial crystallization hit was observed under the following conditions: 0.15 M Lithium sulfate monohydrate, 0.1 M Citric acid, pH 3.5, 18% PEG6000 at all ratios at 288 K. During the optimization (done using “hanging drop” vapor-diffusion method in 24-well VDX plates (Hampton Research) the drop and the precipitant solution volumes were increased by 10 times.

Crystals of the LreCysK were briefly soaked in a mother liquor containing 20% ethylene glycol as a cryoprotectant, and subsequently flash-frozen in liquid nitrogen. The dataset was collected at 100 K at Rigaku OD XtaLAB Synergy-S (IOC RAS, Moscow, Russia). The datasets for the enzyme were indexed, integrated, and scaled using the CrysAlisPro software CrysAlisPro 1.0.43 (Oxford Diffraction/Agilent Technologies UK Ltd., Yarnton, UK). Space group was suggested by Pointless [34] as P2_1_ (Table 4).

The structure of LreCysKapo was obtained using the MOLREP program [35] by molecular replacement with the atomic coordinates of CysK from Planctomyces limnophilus (PDB ID: 5XOQ) as the starting model. The REFMAC5 program of the CCP4 suite (Harwell Science and Innovation Campus, Didcot, UK) [36] was used for the refinement with isotropic B-factors. The model rebuilding was carried out using the COOT interactive graphics program [37].

In the final model, an asymmetric unit contained two copies of the protein (555 visible residues), 350 water molecules, two PEG molecules, 6 ethylene glycol molecules, and 4 SO_4_^2−^ ions from the crystallization solution. 1–4 N-terminal, 152–161 and 213–241 residues in subunit A and 1–5 N-terminal, 155–160 and 212–236 residues in subunit B have no electron density.

3.11. Structure Analysis

The visual inspection was performed with the PyMOL Molecular Graphics System, Version 4.6 (Schrödinger, New York, NY, USA) and the COOT program. Structures were compared and superposed with the PDBeFOLD 2.58 (EMBL-EBI, Hinxton, UK) program [38]. The contact analysis was done using PDBePISA 1.48 (EMBL-EBI, Hinxton, UK) [39].

3.12. Structure Modelling

The model of LreCysK, which contains all missing residues in the crystal structure, was generated using the following procedure. Missing loops were created manually with the Discovery Studio 3.0 software suite (Accelrys, San Diego, CA, USA) using the insertion tool of residues according to the protein sequence received by MALDI-TOF/TOF for correct visualization of the active site, because we can’t see GISA, the active site’s one of the main loops, in our 3D structure. The inserted loop was refined with the “Loop refinement” protocol as implemented in Discovery Studio. Missing residues according to the protein sequence received by MALDI-TOF/TOF were manually added using the insertion tool implemented in the Discovery Studio suite (Accelrys, San Diego, CA, USA). Loop refinement was prepared by a multi-step CHARMm protocol implemented in Discovery Studio. The calculations are based on the ab initio loop prediction algorithm LOOPER. A set of low-energy conformations is generated for the specified loop region and scored based on their CHARMm energy [40].

The structures were visualized using the PyMOL Molecular Graphics System, Version 2.0, Schrödinger, LLC.

4. Conclusions

This study provides comprehensive insights into the kinetic properties, thermal stability, and structural characteristics of cysteine synthase A from L. reuteri LR1. The kinetic parameters were determined using two methods: the previously developed HILIC method [9] and the ninhydrin method, which is commonly employed for most cysteine synthases. A comparative analysis of the data obtained showed that the results of measurements of the catalytic parameters for o-acetylserine are similar. A new method for determining the activity of cysteine synthase can be used for similar enzymes. Structural analysis of LreCysK and its detailed comparison with structures from other organisms revealed that traditionally isolated regions of the amino acid sequence may not be sufficient for a comprehensive structure-functional analysis. Therefore, based on the obtained results, we assume that CysK enzymes have previously uncharacterized amino acid regions that are involved in the enzyme activity. LreCysK’s stability and function in varying environmental conditions make it a promising candidate for medical and industrial applications.

LreCysK’s unique properties, such as high thermal stability and catalytic activity, make it a promising candidate for biotechnological applications. In the industry, it can be used as a producer of cysteine [41]. We think that in the food industry, the enzyme could be incorporated into probiotic formulations to enhance the nutritional and health benefits of fermented products. Its ability to function efficiently in the human gut suggests that it could help to study the mechanism of cysteine and other sulfur-containing compounds production in detail. These compounds are essential for maintaining redox balance and cellular health [26]. This study may help develop new, more efficient methods for producing these compounds. Additionally, the enzyme can serve as a target for antibiotics and as an antimicrobial agent itself [8,42].

Future studies should investigate the engineering of LreCysK to enhance its activity or stability, thereby fully realizing its biotechnological potential. Also, this research on LreCysK could focus on several key areas, including the development of small-molecule inhibitors or activators that could modulate its activity. We assume that these compounds could be useful in studying the enzyme’s role in cellular metabolism or in developing new therapeutic strategies targeting sulfur metabolism in pathogenic bacteria. Additionally, exploring the role of LreCysK in the gut microbiome, particularly its interactions with other microbial species and its impact on host health, could provide valuable insights into the broader implications of cysteine metabolism in the gut.