Intramembrane rhomboid proteases are of particular interest because of their function to hydrolyze a peptide bond of a substrate buried in the membrane. Crystal structures of the bacterial rhomboid protease GlpG have revealed a catalytic dyad (Ser201-His254) and oxyanion hole (His150/Asn154/the backbone amide of Ser201) surrounded by the protein matrix and contacting a narrow water channel. Although multiple crystal structures have been solved, the catalytic mechanism of GlpG is not completely understood. Because it is a serine protease, hydrogen bonding interactions between the active site residues are thought to play a critical role in the catalytic cycle. Here, we dissect the interaction energies among the active site residues His254, Ser201, and Asn154 of Escherichia coli GlpG, which form a hydrogen bonding network. We combine double mutant cycle analysis with stability measurements using steric trapping. In mild detergent, the active site residues are weakly coupled with interaction energies (ΔΔGInter) of ‒1.4 kcal/mol between His254 and Ser201 and ‒0.2 kcal/mol between Ser201 and Asn154. Further, by analyzing the propagation of single mutations of the active site residues, we find that these residues are important not only for function but also for the folding cooperativity of GlpG. The weak interaction between Ser and His in the catalytic dyad may partly explain the unusually slow proteolysis by GlpG compared with other canonical serine proteases. Our result suggests that the weak hydrogen bonds in the active site are sufficient to carry out the proteolytic function of rhomboid proteases.
Rhomboid proteases are a unique class of membrane-integrated enzymes that mediate site-specific proteolysis of the membrane-embedded region of integral membrane proteins where water is scarce (Freeman, 2014). Rhomboids are implicated in a variety of regulatory processes by releasing membrane-bound effector proteins including growth factors, transcription factors, or enzymes, which activate them; in Drosophila melanogaster, Rhomboid-1 regulates early embryonic development by the cleavage of the membrane-bound epidermal growth factor Spitz (Wasserman et al., 2000; Lee et al., 2001). In yeast, rhomboid Rdb1p regulates mitochondrial remodeling and fusion by the cleavage of GTPase Mgm1p (Sesaki et al., 2003). The mammalian homologue, PARL, serves as an antiapoptotic protein by releasing Opa1 protein into the intermembrane space of the mitochondria to maintain a proper level of cytochrome c (Cipolat et al., 2006; Shi and McQuibban, 2017). In Toxoplasma gondii, TgROM rhomboids mediate host–cell invasion processes by cleaving adhesins in internal micronemes, which enables their trafficking to the posterior of parasites (Brossier et al., 2005). In the Gram-negative bacterium Providencia stuartii, the cleavage of twin-arginine translocase A (TatA), a component of the Tat protein secretion pathway, by AarA rhomboid is required for the production of secreted signaling molecules for quorum sensing (Stevenson et al., 2007).
Mutational and inhibitory studies of Rhomboid-1 first suggested that rhomboids are serine proteases (Urban et al., 2001). Rhomboids possess a unique serine–histidine catalytic dyad (Lemberg and Freeman, 2007) rather than the canonical Ser–His–Asp triad found in other serine proteases (Hedstrom, 2002). Crystallographic studies of the rhomboid protease GlpG from Escherichia coli and Haemophilus influenzae have revealed that the active site is buried in the protein matrix ∼10 Å below the membrane surface (Wang et al., 2006; Wu et al., 2006; Ben-Shem et al., 2007; Lemieux et al., 2007). These studies also have shown a narrow water-filled cavity contacting the catalytic dyad (Fig. 1). It has been suggested that this cavity serves as a water retention site, which provides water molecules necessary for catalysis (Zhou et al., 2012). Although the molecular details of how rhomboid proteases carry out proteolysis have not been confirmed, a mechanism has been proposed on the basis of inhibitory, crystallographic, and enzyme kinetic studies along with comparison to canonical serine proteases (Hedstrom, 2002; Ha, 2007; Vinothkumar et al., 2010, 2013; Xue and Ha, 2012, 2013; Brooks and Lemieux, 2013; Dickey et al., 2013; Zoll et al., 2014; Cho et al., 2016; Tichá et al., 2017). First, the hydrogen bond between catalytic His254 and Ser201 (numbering based on E. coli GlpG; dHis,Nε2⋅⋅⋅Ser,Oγ = 2.6 Å) activates the hydroxyl group of the serine for a nucleophilic attack on the carbonyl carbon of the substrate peptide bond to create the first anionic tetrahedral intermediate (Lemberg et al., 2005; Wang et al., 2006). Emerging early as an important residue that may stabilize the intermediate was the conserved asparagine (Asn154) that presumably forms a weak hydrogen bond with the backbone amide group of Ser201 (dAsn,Oδ1⋅⋅⋅Ser,N = 3.3 Å). Crystal structures with peptide inhibitors confirmed that the tetrahedral intermediate is stabilized by the interaction with the unique oxyanion triad composed of the conserved Asn154, His 150, and the backbone of the catalytic Ser201 (Xue and Ha, 2012, 2013; Cho et al., 2016). Next, the intermediate is collapsed, resulting in the formal cleavage of the peptide bond and the formation of an acyl enzyme (Vinothkumar et al., 2010, 2013; Xue and Ha, 2012; Brooks and Lemieux, 2013). His254 may then activate a water molecule to initiate the formation of the second tetrahedral intermediate and active site regeneration (Ha, 2007; Xue and Ha, 2012; Brooks and Lemieux, 2013; Vinothkumar et al., 2013). His254 is separated too far from the crystallographic water molecules near the active site (4–5 Å) to mediate the formation of the second tetrahedral intermediate, such that a large conformational change may occur for an optimal rearrangement of the active site (Vinothkumar et al., 2013).
In canonical serine proteases, the catalytic triad (Ser–His–Asp) forms a tight hydrogen bond network, which coordinates a charge relay necessary for catalysis (Hedstrom, 2002). For example, in chymotrypsin, the strong hydrogen bond between His57 and Asp102 is known to facilitate the nucleophilic attack of Ser195 on the substrate peptide bond and stabilizes the doubly protonated form of His57 (Frey et al., 1994). However, it has been argued whether strong hydrogen bonds are necessary between Ser195 and His57 as well as between His57 and Asp102 for catalysis, or whether weak hydrogen bonds are sufficient (Warshel and Papazyan, 1996; Ash et al., 1997; Stratton et al., 2001; Fuhrmann et al., 2006; Ishida, 2006; Tamada et al., 2009; Petrillo et al., 2012; Agback and Agback, 2018). Therefore, measuring the hydrogen bond strengths in the active site network has been a focus of numerous studies (Frey et al., 1994; Markley and Westler, 1996; Ash et al., 1997; Cleland et al., 1998; Lin et al., 1998; Lau and Bruice, 1999; Frey, 2004; Zheng et al., 2006). Rhomboid proteases lack aspartate, which implies that His254 alone should be sufficient to carry out the activation of Ser201 as well as the subsequent charge relay as a general base. Therefore, measuring the strengths of the active site hydrogen bonds is an important task toward understanding the mechanism of proteolysis mediated by rhomboid proteases.
In this study, we determined the pairwise interaction energies between the active site residues (His254, Ser201, and Asn154) of E. coli GlpG using double mutant thermodynamic cycles combined with stability measurement directly under mild n-dodecyl-β-D-maltopyranoside (DDM) micellar conditions without using chemical denaturants. Interestingly, we show that the catalytic residues in a rhomboid protease are engaged by weak hydrogen bonding interactions compared with those in the canonical serine proteases. The bacterial rhomboid GlpG may represent a unique example among serine proteases in which the strong hydrogen bond network is not required for catalytic proteolysis.
Materials and methods
Expression and purification of GlpG
The transmembrane (TM) domain of GlpG (residues 87–276) encoded by pET15b vector was expressed in E. coli BL21(DE3)-RP strain (Agilent Technologies) with an N-terminal His6-tag (Guo et al., 2016). Cells were grown at 37°C in LB media containing 100 mg/l ampicillin until OD600nm = 0.6 was reached. Then the culture was cooled down on ice for 20 min, and 0.5 mM isopropyl β-thiogalactopyranoside (GoldBio) was added to induce protein expression. The cells were further cultured at 15°C for 16 h. After cell lysis, GlpG was purified from the total membrane fraction obtained by ultracentrifugation (50,000 g for 2 h; Type 45 Ti rotor; Beckman Coulter) using Ni2+-NTA affinity chromatography (Qiagen) after solubilization with 1% DDM (Anatrace). Site-directed mutagenesis was performed using the QuikChange Site-Directed Mutagenesis kit (Agilent Technologies).
Biotin labeling of GlpG
GlpG was labeled with the thiol-reactive biotin derivative possessing pyrene fluorophore, N-(5-(2-iodoacetamido)-6-oxo-6-(2-(+)-Biotin hydrazinyl)hexyl)-4-(pyren-1-yl)butanamide (BtnPyr-IA; Btn: biotin, Pyr: pyrene, IA: iodoacetmaide; Fig. 2 a; Guo et al., 2016). Purified double-cysteine cysteine variant (P95C/G172C or G172C/V267C; hereafter, P95C/G172C is denoted as 95/172N, in which N indicates that two cysteine residues are introduced in the N-terminal half of GlpG, and G172C/V267C is denoted as 172/267C, in which C indicates that two cysteine residues are introduced in the C-terminal half) in 0.2% DDM, 50 mM TrisHCl, and 200 mM NaCl, pH 8.0, buffer were diluted to ∼50 µM and incubated with a 10-fold molar excess Tris(2-carboxyethyl)phosphine–HCl (TCEP; Pierce) for 1 h at room temperature. Then, 40 times molar excess of BtnPyr-IA dissolved in DMSO (∼20 mg/ml) was added to the mixture with gentle vortexing. The labeling reaction was incubated in the dark at room temperature overnight. Excess free label was removed by binding the protein to Ni2+-NTA affinity resin and washing the bound protein with 0.2% DDM, 50 mM TrisHCl, and 200 mM NaCl, pH 8.0, buffer. Labeled GlpG was passed through a desalting column (EconoPac 10DG; Bio-Rad), which was preequilibrated with 0.2% DDM, 50 mM TrisHCl, and 200 mM NaCl, pH 8.0, buffer to remove imidazole. Typically, the labeling efficiency of BtnPyr-IA ranged from 1.5–2.2 as estimated from the concentration of BtnPyr determined by pyrene absorbance (ε346nm = 43,000 M⋅cm−1) and the concentration of GlpG determined by Detergent Compatible protein assay (Bio-Rad; Guo et al., 2016).
Preparation of monovalent streptavidin (mSA)
WT mSA and its variants mSA-S27A and mSA-S45A, in which the mutations were made on the active subunit of tetrameric mSA, were prepared as described previously (Howarth et al., 2006; Hong et al., 2013). Each variant additionally contained a single-cysteine mutation S83C in the active subunit, to which the thiol-reactive dabcyl quencher (DABCYL Plus C2 maleimide; Anaspec) was conjugated for a binding assay between mSA and GlpG-BtnPyr2 using FRET (Guo et al., 2016).
Expression and purification of GlpG substrate SN-LYTM2
For the functional assay of GlpG (Fig. S1), we used the second TM domain of the lactose permease of E. coli (Akiyama and Maegawa, 2007) fused to staphylococcal nuclease (SN-LYTM2) as a model substrate (Fig. S1 a; Lemmon et al., 1992; Hong et al., 2010). SN-LYTM2 containing the SN domain and C-terminal His6-tag was encoded in pET30a vector (Guo et al., 2016). In the LYTM2 region, a single cysteine was engineered at the five residues upstream from the scissile bond (i.e., P5 position) for labeling with thiol-reactive, environment-sensitive fluorophore iodoacetyl-7-nitrobenz-2-oxa-1,3-diazol (IA-NBD amide; Setareh Biotech). The construct was expressed in BL21(DE3)-RP E. coli strain. The protein was expressed, purified, and labeled as described previously (Hong et al., 2013). The cleavage reaction of LYTM2 by GlpG (1 µM) was initiated by addition of 10 times molar excess of NBD-labeled SN-LYTM2 to purified GlpG variants. Time-dependent changes of NBD fluorescence (Fig. S1 b) were monitored in a 96-well plate using a SpectraMax M5e plate reader (Molecular Devices) with excitation and emission wavelengths of 485 nm and 535 nm, respectively. The decrease in fluorescence intensity, which indicates the transfer of the environment-sensitive NBD fluorophore from the hydrophobic micellar phase to the aqueous phase upon cleavage, was normalized to a control sample with NBD-SN-LYTM2 only. The initial slope of fluorescence change versus time represents the substrate cleavage rate.
Construction of binding isotherm to determine thermodynamic stability of GlpG using steric trapping
Thermodynamic stability (ΔGoU) of GlpG in DDM micelles (20 mM) was determined by measuring the attenuated second binding of mSA labeled with dabcyl quencher (mSADAB) to GlpG doubly labeled with BtnPyr (95/172N-BtnPyr2 or 172/267C-BtnPyr2) at room temperature (Guo et al., 2016). mSADAB binding was monitored by quenching of pyrene fluorescence from BtnPyr labels by FRET. 1 µM of 95/172N-BtnPyr2 or 172/267C-BtnPyr2 was titrated with mSADAB variant with a reduced biotin binding affinity, mSADAB-S45A (Kd,biotin = 9.0 ± 4.3 nM) or mSADAB-S27A (Kd,biotin = 1.4 ± 0.9 nM) in 20 mM DDM, 0.25 mM Tris(2-carboxyethyl)phosphine, 20 mM sodium phosphate, and 200 mM NaCl (pH 7.5; Guo et al., 2016). The use of mSA variants was necessary to achieve the reversibility of the second mSA binding to obtain ΔGoU. The titrated samples were transferred to a 96-well UV-compatible microplate, sealed with a polyolefin tape, and incubated for 5 d (for 95/172N-BtnPyr2) or 2 d (for 172/267C-BtnPyr2) at room temperature. Quenching of pyrene-monomer fluorescence at 390 nm was monitored with an excitation wavelength of 345 nm on a SpectraMax M5e plate reader. Data were averaged from three readings. Nonspecific FRET between pyrene and dabcyl was negligible (Guo et al., 2016).
Fitting of binding isotherm to determine thermodynamic stability of GlpG
F and U denote the folded and unfolded state, respectively. KU is the equilibrium constant for unfolding of GlpG, and Kd,biotin is the dissociation constant for unhindered biotin binding of mSA. The fitting equation for the second mSA binding phase in the binding isotherm was
where FL is measured fluorescence intensity, and FL0 and FL∞ are the fluorescence intensities at [mSA] = 0 and at the saturated binding level, respectively. [mSA] is the total mSA concentration. After obtaining the fitted KU, the thermodynamic stability was calculated using the equation ΔGoU = –RT lnKU.
Cooperativity profiling of the active site residues
The following is the method to identify cooperative and localized side-chain interactions that contribute to the protein stability (Guo et al., 2016). To apply this method to GlpG, we first made a single mutation to perturb a specific side-chain interaction in the background of double biotin variants of GlpG, 95/172N-BtnPyr2 and 172/267C-BtnPyr2. Next, using steric trapping, the stability changes induced by the mutation are measured with two different biotin pairs that are located in the N and C subdomains, respectively. The differential effect of the same mutation on the stability of each subdomain (ΔΔΔGoU) is quantified (ΔΔΔGoU) using Eq. 4:
ΔΔGoU,95/172N-BtnPyr2(WT-Mut) and ΔΔGoU,172/267C-BtnPyr2(WT-Mut) designate the stability changes caused by the same mutation in the backgrounds of 95/172N-BtnPyr2 and 172/267C-BtnPyr2, respectively. If the mutation causes a similar degree of destabilization for both double-biotin variants with a difference smaller than thermal fluctuation energy (∣ΔΔΔGoU∣ ≤ RT = 0.6 kcal/mol; R: gas constant; T = 298K), the mutated site engages in a “cooperative” interaction. Among the cases where∣ΔΔΔGoU∣ RT, if the mutation preferentially destabilizes the subdomain containing it, the perturbed interactions are “localized” within that subdomain. If the mutation of a residue, which makes its side-chain contacts only with the subdomain containing it, preferentially destabilizes the other subdomain, we define that perturbation as “over-propagated.”
Online supplemental material
Figure S1 displays an activity assay of GlpG WT and active site mutants. Figure S2 displays a comparison of the GlpG structures determined in detergent and bicelles. Table S1 displays stability changes induced by single Ala mutations on the active site residues.
Rationale of double mutant cycle analysis
A double mutant cycle involves WT protein, two single mutants, and the corresponding double mutant. If the change in thermodynamic stability (ΔGoU) upon the double mutation differs from the sum of the changes caused by the single mutations, the two residues in WT are coupled, and the magnitude of the difference (interaction energy: ΔΔGInter) is related to the strength of interaction between them (Horovitz, 1996), such that
where X and Y denote a residue pair of interest in WT, and A designates Ala. The thermodynamic scheme using a double mutant is advantageous for quantifying the strength of a specific inter-residue interaction in the context of the native WT structure: a single mutation disrupts not only the interaction between a specific residue pair of interest but also the interaction between the mutated residue and its environment. Also, the single mutation may induce global or local structural relaxation. These energetic contributions other than the specific inter-residue interaction can be subtracted out by measuring the stability change induced by the same mutation in the absence of the partner residue (Eq. 5; Fersht et al., 1992). Double mutant cycle analysis has been widely used to measure the strengths of intramolecular and intermolecular side-chain interactions for both globular and membrane proteins (Serrano et al., 1990, 1991; Doura and Fleming, 2004; Hong et al., 2006, 2007; Harel et al., 2007; Joh et al., 2008; Sokolovski et al., 2017). To minimize the possibility of creating new interactions after mutation, we replaced each active site residue with alanine (Horovitz, 1996).
Mild destabilization by single alanine mutations in the active site residues
To calculate ΔΔGInter using double mutant cycle analysis, we measured ΔGoU’s of WT, single-Ala, and double-Ala mutants of GlpG using a steric trapping, which couples transient unfolding of a doubly biotinylated protein to double binding of bulky mSA (52 kD; see Fig. 2 a for a more detailed description of the principle; Hong et al., 2010; Guo et al., 2016). Compared with conventional stability measurements using chemical denaturants, this method is advantageous because the protein stability can be directly measured in a native solvent and lipid environment. Previously, we have identified optimal sites of thiol-specific biotinylation on GlpG for steric trapping, P95C/G172C (95/172N: N indicates the N-terminal subdomain where a biotin pair is located; Fig. 2 b). The unfolded state trapped with this biotin pair is globally denatured and reversibly refolds to the native state upon addition of excess free biotin to dissociate bound mSA (Guo et al., 2016). For measuring ΔGoU of GlpG, a binding isotherm between doubly biotinylated GlpG and mSA is obtained by using the thiol-reactive biotin derivative with a pyrene fluorophore (BtnPyr) and mSADAB (Guo et al., 2016; Fig. 2 a). When an mSADAB variant with a reduced biotin binding affinity is used, the binding isotherm monitored by quenching of pyrene fluorescence displays two-phase mSA binding after the binding equilibrium has been reached: the first mSA tightly binds to either biotin label with an intrinsic binding affinity, and the second mSA binds with a weaker affinity because of the coupling to GlpG unfolding. ΔGoU of GlpG is determined by fitting the second binding phase to Eq. 3 (see “Materials and methods”).
The binding isotherms using weaker biotin-binding mSA variants (mSADAB-S27A or mSADAB-S45A; Guo et al., 2016) are shown in Fig. 2 b. The single Ala mutations at the active site residues Ser201, His254, and Asn154 completely abolished GlpG activity, as previously reported (Clemmer et al., 2006; Baker and Urban, 2012; Dickey et al., 2013; Fig. S1 b). In the crystal structures determined in detergent, these residues are completely buried in the protein matrix (Fig. S2 a). ΔG°U of the double-biotin variant without additional Ala mutation (i.e., WT) was 5.6 ± 0.1 kcal/mol. Single mutants S201A and H254A were mildly destabilized with ΔΔG°U’s of 1.1 ± 0.1 and 0.7 ± 0.1 kcal/mol, respectively. Ala mutation at Asn154 induced larger destabilization with ΔΔG°U = 1.5 ± 0.1 kcal/mol. Overall, mutations in the active site did not induce substantial destabilization relative to other previously characterized Ala mutations in the buried region of GlpG, for which ΔΔG°U can be as large as ∼4 kcal/mol (Guo et al., 2016). Overall, the mild destabilization by the single active site mutations obtained by steric trapping agrees with previous studies using SDS-induced or irreversible thermal denaturation (Table S1; Baker and Urban, 2012; Paslawski et al., 2015).
Next, the stabilities of double-Ala mutants were measured. Interestingly, the double mutation on the catalytic dyad (S201A/H254A) yielded a smaller decrease in the stability (ΔΔG°U = 0.4 ± 0.2 kcal/mol) than individual single mutations (ΔΔG°U = 0.7–1.1 kcal/mol), indicating that the double mutation caused a certain extent of structural relaxation. The double mutations N154A/H254A and N154A/S201A induced larger destabilization (ΔΔG°U = 1.8‒2.4 kcal/mol) than individual single mutations, implying an additive effect of the single mutations.
Weak interaction energies between the active site residues of GlpG
Next, we determined the interaction strengths (ΔΔGInter) between the active site residue pairs using double mutant cycles (Fig. 3; Horovitz, 1996). From this analysis, Ser201 and His254, which form the catalytic dyad and are engaged in a close hydrogen bond (dHis,Nε2⋅⋅⋅Ser,Oγ = 2.6 Å), favorably interacted (ΔΔGInter = ‒1.4 ± 0.2 kcal/mol), whereas the interaction between Ser201 and Asn154, which form a more distant side-chain–backbone hydrogen bond (dAsnOδ1⋅⋅⋅SerN = 3.3 Å), was not significant (ΔΔGInter = ‒0.2 ± 0.2 kcal/mol). ΔΔGInter between His254 and Asn154, which are apparently not engaged in any interaction, was not significant either (‒0.4 ± 0.2 kcal/mol).
Hydrogen bonds can be categorized according to the strength of their interaction: weak or conventional (2‒12 kcal/mol), strong or low-barrier (12‒24 kcal/mol), and very strong or single-well (24 kcal/mol; Frey et al., 1994). Although favorable, the measured hydrogen bond strength of the His-Ser catalytic dyad of GlpG (ΔΔGInter ≈ ‒1.4 kcal/mol) is regarded as weak. This interaction is not only substantially weaker than the His57-Asp102 interaction in the active site of chymotrypsin or α-lytic protease (‒5 to ‒7 kcal/mol before the formation of the first tetrahedral intermediate and ‒7 to ‒10 kcal/mol in the first tetrahedral intermediate; Ash et al., 1997; Frey, 2004), which has been suggested to form a low-barrier hydrogen bond, but also weaker than the Ser195-His57 interaction in chymotrypsinogen, classified as moderately strong at an acidic pH (~–13 kcal/mol; Frey et al., 1994; Markley and Westler, 1996).
Active site residues are involved in cooperative interactions
Finally, we analyzed the contribution of each active site residue to the folding cooperativity of GlpG using the steric trapping-based cooperativity profiling (Fig. 4; Guo et al., 2016). This method is based on the principle that steric trapping captures the transient unfolding of the tertiary interactions in the region to which a specific biotin pair is conjugated. Thus, the local stability of a protein can be measured, and how the local sequence perturbation caused by mutation is propagated throughout the protein structure can be quantified. Briefly, the effect of a specific mutation on the stability (ΔΔGoU) is measured with two biotin pairs located in different regions. If the difference in the measured stability changes (ΔΔΔGoU) is smaller than thermal fluctuation energy (i.e., ∣ΔΔΔGoU∣ ≤ RT = 0.6 kcal/mol), it indicates that the side-chain perturbation by the mutation is propagated evenly throughout the protein and the mutated side chain is engaged in cooperative interactions. If the mutation preferentially destabilizes the subdomain that includes the mutation site (∣ΔΔΔGoU∣ RT), the mutated side chain is engaged in localized interactions. If the mutation preferentially destabilizes the subdomain that does not include the mutation site with ∣ΔΔΔGoU∣ RT, the mutated side chain is engaged in over-propagated interactions.
To apply this method to the active site residues of GlpG, the stability changes upon each single alanine mutation were measured at the biotin pairs 95/172N-BtnPyr2 (Fig. 2 b) and 172/267C-BtnPyr2 (C indicates the C-terminal subdomain in which the biotin pair is located), respectively (Fig. 4 a; Guo et al., 2016). We have shown that the global stability of GlpG is not affected by the presence of these biotin pairs (Guo et al., 2016). S201A mutation at the subdomain interface similarly destabilized the N and C subdomains (ΔΔGoU = 1.1 ± 0.1 kcal/mol and 0.6 ± 0.2 kcal/mol, respectively), yielding ∣ΔΔΔGoU∣ = 0.5 ± 0.2 RT. Thus, Ser201 is classified as cooperatively engaged (Fig. 4 b). Interestingly, H254A mutation in the C subdomain induced destabilization of the N subdomain (ΔΔGoU = 0.7 ± 0.1 kcal/mol), whereas the same mutation stabilized the C subdomain containing the mutated site (ΔΔGoU = −0.8 ± 0.2 kcal/mol), yielding ∣ΔΔΔGoU∣ = 1.5 ± 0.2 kcal/mol (2RT). Thus, we assign His254 interactions as highly overpropagated. We reason that the stabilization of the C subdomain by the mutation H254A is due to the global structural relaxation induced by the large changes in the side-chain volume and polarity (Fig. 2 b and Fig. 4 b). The mutation N154A preferentially destabilized the N subdomain, where the mutation resides. The resulting ∣ΔΔΔGoU∣of 0.7 ± 0.2 kcal/mol was slightly larger than RT. Thus, we assign Asn154 interactions as moderately localized. Our analysis indicates that the absolutely conserved catalytic dyad Ser201–His254 is not only critical for function but also highly communicative with its environment to maintain the folding cooperativity of GlpG.
Here we have shown that the hydrogen bond network in the active site of GlpG is maintained by weak side-chain interactions for its assembly and function. Notably, the strengths of the buried side-chain hydrogen bonds obtained in this study fall into the range of those measured in globular and membrane proteins (0–2.0 kcal/mol; Fleming and Engelman, 2001; Gratkowski et al., 2001; Takano et al., 2003; Hong et al., 2006, 2010; Stanley et al., 2006; Joh et al., 2008; Bowie, 2011; Baker and Urban, 2012). Whereas the previous efforts for membrane proteins have been mainly concerned with “structural” hydrogen bonds that are not directly involved in function, we present the first example of measuring the strengths of “functional” hydrogen bonds conserved in the rhomboid protease family. Our result provides important chemical insights into the initial step of the proteolysis mechanism by rhomboid proteases, i.e., the activation of Ser201 by His254 and the stabilization of the anionic tetrahedral intermediate by Asn154 forming a part of the oxyanion hole. The weak hydrogen bond between His254 and Ser201 must be sufficient to activate Ser201 for the nucleophilic attack on the peptide bond. In addition, the negligible interaction between Asn154 and Ser201 implies that Asn154 will easily gain the flexibility to be adapted to a conformation that can stabilize the oxyanion intermediate.
The kinetic study of GlpG function in the lipid bilayers has shown that the proteolytic activity of GlpG can be characterized as a low substrate-binding affinity (KM ∼135 µM or 0.001 mol fraction, [substrate]/[lipid]), a slow catalytic rate (kcat = 0.0063 s−1), and a low efficiency (kcat/KM ∼47 M−1s−1) compared with those of other well-studied serine proteases such as chymotrypsin, trypsin, elastase, and α-lytic protease (Dickey et al., 2013). These other proteases, whose catalytic mechanisms are known to be driven by tightly coupled catalytic triads, possess the kinetic parameters of KM ∼101–104 μM, kcat ∼3–30 s−1, and kcat/KM ∼103–107 M−1s−1 (Brothers and Kostić, 1990; Tsu and Craik, 1996; Coombs et al., 1999). The slow catalytic rate (kcat) by GlpG has been mainly attributed to the slow opening of the gating helix TM5, which controls the access of the substrate to the active site in the membrane (Dickey et al., 2013). Although it has been suggested that the hydrolysis reaction is not a rate-determining step responsible for the low kcat, the weak hydrogen bond between His254 and Ser201 demonstrated in this study may partly explain the vastly low catalytic ability of GlpG (104–105 fold lower kcat than other robust serine proteases).
Here, the side-chain hydrogen bonds were measured in the detergent micellar phase. The possibility remains that the strengths of the same interaction can be different in the lipid bilayers, probably because of the effect of the lateral packing pressure of lipid molecules (Cantor, 1997) as well as the possible difference in the water dynamics that can compete with the active site hydrogen bonds. The thermodynamic stability of helical bundle membrane proteins has been measured in the lipid bilayer environments for the membrane transporter LeuT using urea denaturation (in liposomes), the proton pump bacteriorhodopsin using steric trapping (in DMPC/CHAPSO bicelles), and GlpG using single-molecule force spectroscopy (in DMPC/CHAPSO bicelles; Chang and Bowie, 2014; Min et al., 2015; Sanders et al., 2018). However, the contribution of specific side-chain interactions to the stability has not been quantified in the bilayers. Although the reversible folding of GlpG has not yet been achieved in the bilayer using steric trapping, we expect that the strengths of the hydrogen bonds determined in detergent may not be much different from those in the lipid bilayers for the following reasons. (a) From the crystal structures of GlpG determined in detergent and bilayers (i.e., bicelles), the three active site residues are largely buried in the protein matrix, implying that the inter-residue interactions would be maintained to a similar extent regardless of the surrounding lipid environment (Fig. S2 a; Wang et al., 2006; Vinothkumar, 2011; Cho et al., 2016). (b) The structures of GlpG in detergent and bicelles are virtually identical (Fig. S2, b and c; Wang et al., 2006; Vinothkumar, 2011; Cho et al., 2016). Although the structure of WT in bicelles is not available, the structural comparison of the inactive mutant S201T determined in detergent to that in bicelles indicates that the geometry of the active site residues and the network of the water-mediated hydrogen bonds are similar in the two environments (Fig. S2 b; Vinothkumar, 2011). (c) Although it has been reported that the lipid bilayers enhance the activity of GlpG by ∼twofold relative to detergent, the detergent environment well support the proteolytic activity in a similar timescale (tens of minutes) for known rhomboid substrates (Strisovsky et al., 2009; Baker and Urban, 2012; Moin and Urban, 2012; Xue and Ha, 2013). (d) For the single membrane–spanning TM helices, the strengths of the intermolecular hydrogen bonds determined in detergents are similar to those in the lipid bilayers (1.0 ± 0.5 kcal/mol; Bowie, 2011). For example, the hydrogen bond mediated by glutamate residue in the TM helix of a fibroblast growth factor receptor contributes to the dimer stability by ∼0.7 kcal/mol per monomer in the natural cell membranes (Li et al., 2006). The side chain–backbone intermolecular hydrogen bond mediated by a threonine residue stabilizes the dimer of the glycophorin A TM helix by 1 kcal/mol per monomer (Hong et al., 2010). However, to clarify this issue, it is still necessary to determine the thermodynamic stability of GlpG and the strengths of the hydrogen bonds in a lipid bilayer. Steric trapping is a promising tool to measure the membrane protein stability without disrupting the bilayers. Currently, developing such a protocol for achieving the reversible folding of GlpG in the bilayer is under progress.
We thank Ruiqiong Guo and Miyeon Kim for providing the reagent BtnPyr-IA and critically reading the manuscript.
This work was supported by National Institutes of Health grant R01GM118685, a Michigan State University start-up fund, and the Hunt for a Cure foundation grant RG070103 to H. Hong.
The authors declare no competing financial interests.
Author contributions: K. Gaffney and H. Hong conceived the project and designed experiments. K. Gaffney carried out the experiments and data analysis. K. Gaffney and H. Hong prepared the figures and wrote the manuscript.
José D. Faraldo-Gómez served as editor.
This work is part of the special collection entitled “Molecular Physiology of the Cell Membrane: An Integrative Perspective from Experiment and Computation.”