Title: An in vitro platform for the enzymatic characterization of the rhomboid protease RHBDL4
Rhomboid proteases are ubiquitous intramembrane serine proteases that can cleave transmembrane substrates within lipid bilayers. They exhibit many and diverse functions, such as but not limited to, growth factor signaling, immune and inflammatory response, protein quality control, and parasitic invasion. Human rhomboid protease RHBDL4 has been demonstrated to play a critical role in removing misfolded proteins from the endoplasmic reticulum and is implicated in severe diseases such as various cancers and Alzheimer's disease. Therefore, RHBDL4 is expected to constitute an important therapeutic target for such devastating diseases. Despite its critical role in many biological processes, the enzymatic properties of RHBDL4 remain largely unknown. To enable a comprehensive characterization of RHBDL4's kinetics, catalytic parameters, substrate specificity, and binding modality, we expressed and purified recombinant RHBDL4 and employed it in a Förster resonance energy transfer-based cleavage assay. Until now, kinetic studies have been limited mostly to bacterial rhomboid proteases. Our in vitro platform offers a new method for studying RHBDL4's enzymatic function and substrate preferences. Furthermore, we developed and tested potential inhibitors using our assay and successfully identified peptidyl α-ketoamide inhibitors of RHBDL4 that are highly effective against recombinant RHBDL4. We utilize ensemble docking and molecular dynamics simulations to explore the binding modality of substrate-derived peptides bound to RHBDL4. Our analysis focused on key interactions and dynamic movements within RHBDL4's active site that contributed to binding stability, offering valuable insights for optimizing the nonprime side of RHBDL4 ketoamide inhibitors. In summary, our study offers fundamental insights into RHBDL4's catalytic activities and substrate preferences, laying the foundation for downstream applications such as drug inhibitor screenings and structure-function studies, which will enable the identification of lead drug compounds for RHBDL4.  more » « less
Award ID(s):
2047391
PAR ID:
10578292
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Neal, Sonya E
Date Published:
Journal Name:
Journal of Biological Chemistry
Volume:
301
Issue:
3
ISSN:
0021-9258
Page Range / eLocation ID:
108275
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The main protease (Mpro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of Mpro, a cysteine protease, have been determined, facilitating structure-based drug design. Mpro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41-Cys145, Mpro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nu-cleophile Cys145 have been debated in previous studies of SARS-CoV Mpro, but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 Mpro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of Mpro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α-ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored Nδ (HD) and Nϵ (HE) protonation of His41 and His164, respectively, the α-ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 Mpro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts. 
    more » « less
  2. The main protease (M pro ) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of M pro , a cysteine protease, have been determined, facilitating structure-based drug design. M pro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41–Cys145, M pro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nucleophile Cys145 have been debated in previous studies of SARS-CoV M pro , but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 M pro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of M pro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α-ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored N δ (HD) and N ϵ (HE) protonation of His41 and His164, respectively, the α-ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 M pro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts. 
    more » « less
  3. null (Ed.)
    Aromatase (CYP19) catalyzes the last biosynthetic step of estrogens in mammals and is a primary drug target for hormone-related breast cancer. However, treatment with aromatase inhibitors is often associated with adverse effects and drug resistance. In this study, we used virtual screening targeting a predicted cytochrome P450 reductase binding site on aromatase to discover four novel non-steroidal aromatase inhibitors. The inhibitors have potencies comparable to the noncompetitive tamoxifen metabolite, endoxifen. Our two most potent inhibitors, AR11 and AR13, exhibit both mixed-type and competitive-type inhibition. The cytochrome P450 reductase-CYP19 coupling interface likely acts as a transient binding site. Our modeling shows that our inhibitors bind better at different sites near the catalytic site. Our results predict the location of multiple ligand binding sites on aromatase. The combination of modeling and experimental results supports the important role of the reductase binding interface as a low affinity, promiscuous ligand binding site. Our new inhibitors may be useful as alternative chemical scaffolds that may show different adverse effects profiles than current clinically used aromatase inhibitors. 
    more » « less
  4. Chemical inhibitors bind to enzymes, thereby inhibiting their catalytic activity. While many enzymes catalyze reactions with a single substrate, others, like DNA polymerase, can act on multiple related substrates. Substrate-selective inhibitors (SSIs) target these multisubstrate enzymes to modulate their specificity. Although SSIs hold promise as therapeutics, our theoretical understanding of how different inhibitors influence enzyme specificity remains limited. In this study, we examine enzyme selectivity within kinetic networks corresponding to known inhibition mechanisms. We demonstrate that competitive and uncompetitive inhibitors do not affect substrate specificity, regardless of rate constants. In contrast, noncompetitive and mixed inhibition can alter specificity and can lead to nonmonotonic responses to the inhibitor. We show that mixed and noncompetitive inhibitors achieve substrate-selective inhibition by altering the effective free-energy barriers of product formation pathways that are enabled by the inhibitor's presence. We then apply this framework to the Sirtuin-family deacylase SIRT2, showing that the suicide inhibitor thiomyristoyl lysine (TM) cannot influence substrate specificity unless there is a direct substrate exchange reaction or biochemical constraints are relaxed. These findings provide insights into engineering systems where cofactor binding modulates metabolic flux ratios. 
    more » « less
  5. Abstract 1,4-Dioxane is a widespread groundwater contaminant frequently co-occurring with chlorinated solvents. The group-6 propane monooxygenase (PRM) from Mycobacterium dioxanotrophicus PH-06 degrades dioxane efficiently, yet the molecular determinants underlying its broad substrate spectrum and inhibition behavior remain unresolved. Here, we combined AlphaFold2-based structure prediction, molecular docking, and alchemical free energy calculations to systematically characterize ligand binding in PH-06 PRM. The catalytic center (residues 87–265), coordinated by two histidines and four glutamates, exhibited calculated binding free energies for five substrates (C2–C4 alkanes and cyclic ethers) that significantly correlated with experimentally determined degradation rates (p < 0.05), supporting the ability of ligand binding thermodynamics to capture major substrate-recognition trends in PH-06 PRM. In addition, we identified a remote binding site (residues 329–446) approximately 30 Å from the catalytic center. Thermodynamic analysis indicates that this remote site can attract methane and ethane, potentially limiting their access to the catalytic center. Comparative binding analyses of three common chlorinated solvent inhibitors (i.e., trichloroethylene [TCE], 1,1-dichloroethylene [DCE], and 1,1,1-trichloroethane [TCA]) revealed distinct site preferences consistent with their experimentally observed inhibition modes. Noncompetitive inhibitors (DCE and TCA) displayed stronger affinity for the remote site, whereas the competitive inhibitor TCE exhibited comparable affinity for both binding regions. Sequence conservation, contact probability analysis, and Gene Ontology enrichment further suggest the structural relevance of the remote site and its possible role in ligand recognition and inhibitor selectivity. These findings demonstrate that the AlphaFold-guided free energy calculations can capture experimentally observed substrate and inhibitor trends in PH-06 PRM, while identifying a remote ligand-binding pocket that warrants further mechanistic investigation. 
    more » « less