Daiki
Yamane
a,
Satsuki
Onitsuka
a,
Suyong
Re
b,
Hikaru
Isogai
a,
Rui
Hamada
a,
Tadanari
Hiramoto
a,
Eiji
Kawanishi
a,
Kenji
Mizuguchi
bc,
Naoya
Shindo
*a and
Akio
Ojida
*a
aGraduate School of Pharmaceutical Sciences, Kyushu University, 3-1-1 Maidashi, Higashi-ku, Fukuoka 812-8582, Japan. E-mail: shindo708@phar.kyushu-u.ac.jp; ojida@phar.kyushu-u.ac.jp
bArtificial Intelligence Center for Health and Biomedical Research, National Institute of Biomedical Innovation, Health and Nutrition, 7-6-8 Saito-Asagi, Ibaraki, Osaka 567-0085, Japan
cInstitute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
First published on 15th February 2022
The coronavirus disease 2019 (COVID-19) pandemic has necessitated the development of antiviral agents against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The main protease (Mpro) is a promising target for COVID-19 treatment. Here, we report an irreversible SARS-CoV-2 Mpro inhibitor possessing chlorofluoroacetamide (CFA) as a warhead for the covalent modification of Mpro. Ugi multicomponent reaction using chlorofluoroacetic acid enabled the rapid synthesis of dipeptidic CFA derivatives that identified 18 as a potent inhibitor of SARS-CoV-2 Mpro. Among the four stereoisomers, (R,R)-18 exhibited a markedly higher inhibitory activity against Mpro than the other isomers. Reaction kinetics and computational docking studies suggest that the R configuration of the CFA warhead is crucial for the rapid covalent inhibition of Mpro. Our findings highlight the prominent influence of the CFA chirality on the covalent modification of proteinous cysteines and provide the basis for improving the potency and selectivity of CFA-based covalent inhibitors.
Mpro recognizes glutamine and a hydrophobic amino acid residue at the S1 and S2 pockets in the active site, respectively, and cleaves amide bonds primarily within the Leu-Gln↓(Ser, Ala, Gly) sequence (↓: cleavage site).8 This substrate specificity is conserved across other coronaviruses, including SARS-CoV-1 and middle east respiratory syndrome coronavirus (MERS-CoV).9,10 Therefore, current efforts to develop SARS-CoV-2 Mpro inhibitors are extensively guided by the previous molecular designs, especially those targeting SARS-CoV-1 Mpro.5,6 Most of the reported SARS-CoV-2 Mpro inhibitors are peptidomimetics based on the Leu-Gln sequence, combined with an electrophilic warhead, such as α-ketoamide,11 aldehyde,12–15 ketone,16–19 vinyl sulfone,8 and nitrile,20,21 for the covalent capture of the sulfhydryl group of the catalytic Cys145. Nirmatrelvir (PF-07321332), developed by Pfizer, is one of the most advanced compounds in this category and has been recently approved by regulatory ageincies as a combination with ritonavir (Fig. 1A).21 Other than peptidomimetics, several Ugi multicomponent reaction (MCR)-generated dipeptidic compounds have been reported as inhibitors against SARS-CoV-1 and SARS-CoV-2 Mpro (Fig. 1B).22–24
Fig. 1 Chemical structures of SARS-CoV-2 Mpro inhibitors. (A) Structure of nirmatrelvir (PF-07321332) and its binding mode in the Mpro active site (the conformation was retrieved from the cocrystal structure, PDB code 7RFW). (B) The structures of dipeptidic Mpro inhibitors generated by Ugi multicomponent reaction and the binding mode of Jun8-76-3A in the Mpro active site (the conformation was retrieved from the cocrystal structure, PDB code 7KX5). |
Covalent inhibition of proteins has been a powerful strategy for achieving potent and sustained pharmacological efficacy.25 To date, a number of acrylamide-based covalent inhibitors targeting a noncatalytic cysteine of proteins have been clinically approved or under clinical trials for cancer treatment.26,27 In development of covalent inhibitors, high target selectivity should be ensured to minimize the risk of toxicity due to off-target labeling.28 This point is also crucial for designing covalent inhibitors targeting SARS-CoV-2 Mpro, as medications for infectious diseases generally require a broader safety margin than cancer treatment. We have recently introduced α-chlorofluoroacetamide (CFA) as a weakly reactive, cysteine-directed warhead for covalent inhibition of proteins including tyrosine kinases.29–31 We found that the CFA derivatives exhibited higher selectivity in covalent modification toward the targeted tyrosine kinases than the structurally related acrylamide-based inhibitors. We also demonstrated that the CFA–thiol adduct is hydrolyzed under neutral aqueous conditions to reversibly generate an unmodified cysteine, which could contribute to high target selectivity by eliminating the solvent-exposed off-target labeling. These desirable features led us to envision CFA as a valuable tool for designing selective covalent SARS-CoV-2 Mpro inhibitors.
This observation encouraged us to subsequently conduct an SAR study to improve the potency of CFA dipeptides (Fig. 3A). We first evaluated the effect of the para-substituent of the benzene ring (R1) on inhibitory activity. As this moiety occupies the S2 pocket that primarily accommodates a leucine residue, we screened compounds 6–9 with a hydrophobic para-substituent. In the assay, less polar diastereomer A was inactive (IC50 > 20 μM), whereas more polar diastereomer B showed moderate activity in all the cases. Among them, 9B possessing a pentafluorosulfanyl (SF5) group, a thermally and chemically stable bioisostere for a tert-butyl group,33 exhibited a slightly higher Mpro inhibitory activity (IC50 = 0.90 μM) than 3B. Next, we investigated the effect of the heteroaromatic ring (R2) on the activity. The replacement of the pyridine ring of 3B with a diazine increased the potency of 11–12, whereas the 3-pyridazinyl derivative 10 was almost inactive (IC50 > 20 μM as a diastereomixture). 5-Pyrimidine 12 was the only compound with separable diastereomers, and 12B exhibited the highest activity in the series (IC50 = 0.25 μM). The effect of the substituent R3 was also investigated. The introduction of a tert-butyl group, as observed in the structure of ML188, reduced the potency (14B, IC50 = 2.70 μM). Replacement of the 3-fluorophenyl group of 3 to a 3-pyridyl group slightly increased the activity (16, IC50 = 0.73 μM as a diastereomixture). The introduction of a 7-isoquinoline ring resulted in a slight loss of activity (17B, IC50 = 1.43 μM). Combining the SAR results, we synthesized compound 18 and 19 as structurally optimized CFA inhibitors (Fig. 3B). The diastereomers were separable in both cases and the active isomers 18B and 19B exhibited potent Mpro inhibitory activities (IC50 = 0.21 and 0.14 μM, respectively). Although 19B was marginally more active than 18B, we chose 18 for further investigations because the chromatograhic separation of the diastereomixture of 18 was much easier than that of 19. Fortunately, we were able to determine the relative configuration of 18B by X-ray crystallography (Fig. S1†). The analyzed structure revealed that 18B has a syn-configuration, where the 5-pyrimidinyl group and the CFA fluorine atom are oriented toward the same direction in the drawn structure (Fig. 3C). Consequently, 18A was determined to be an anti-isomer.
Fig. 3 Structure–activity relationship of chlorofluoroacetamide (CFA) derivatives. (A) Screening of the CFA derivatives with the different R1, R2, and R3 substituents. ND, not determined. A and B indicate the less and more polar diastereomers separated by silica gel column chromatography. (B) The structures and inhibitory activities of the CFA dipeptides 18 and 19. (C) The relative configuration of 18B determined by X-ray crystallography (see Fig. S1†). |
To elucidate the relationship between the CFA chirality and the inhibitory activity, we next attempted to obtain the four stereoisomers of 18 in enantiomerically pure form using optically active chlorofluoroacetic acid.34 Wakselman previously reported the preparation of optically active chlorofluoroacetyl chloride based on the chromatographic separation of the diastereomers of CFA 20 derived from (S)-(−)-1-phenethylamine.35 However, the isolation of diastereomers of 20 proved to be laborious in our attempts due to poor separation by silica gel chromatography. Thus, we established a new method employing (R)-(−)-2-phenylglycinol as a chiral resolving agent (Scheme 1). We found that the diastereomers of CFA 21, readily accessible from commercially available (R)-(−)-2-phenylglycinol and ethyl chlorofluoroacetate, were easily separated by silica gel chromatography and afforded pure 21A and 21B on a multigram scale. X-ray crystallography determined the absolute configuration at the CFA chiral center to be R in 21B (Fig. S2†). Isolated 21A and 21B were subjected to hydrolysis in refluxing 3 M H2SO4/1,4-dioxane to afford (S)- and (R)-chlorofluoroacetic acid, respectively. Finally, Ugi MCR using enantiopure chlorofluoroacetic acids and subsequent chromatographic separation of the diastereomers yielded a set of four stereoisomers of 18.
Scheme 1 The optical resolution of chlorofluoroacetic acid using (R)-(−)-2-phenylglycinol as a chiral resolving agent. See ESI† for the detail. |
Having the four stereoisomers of 18 in hand, we evaluated their inhibitory activities against SARS-CoV-2 Mpro (Fig. 4). We found that (R,R)-18 (the former R indicates the absolute configuration at the Ugi MCR-generated stereocenter, and the latter at CFA) exclusively exhibited strong activity (IC50 = 0.056 μM) among the four stereoisomers (Fig. 4B and D). Notably, this activity is much higher than that of (R,S)-18 possessing an (S)-CFA warhead (IC50 = 10.62 μM), indicating that the chirality at CFA greatly influences the inhibitory activity. The IC50 values of the other isomers (S,R)-18 and (S,S)-18 were 2.24 and 14.04 μM, respectively. To further investigate these results, we performed a kinetic analysis of Mpro inhibition by (R,R)-18 and its diastereomers (S,R)- and (R,S)-18.36 All three isomers inhibited the hydrolysis of the fluorogenic substrate by Mpro in a time-dependent manner, indicating an irreversible mode of action (Fig. 4C and S3–S5†). The kinetic parameters for irreversible inactivation of Mpro are summarized in Fig. 4D. As expected, the most potent isomer (R,R)-18 inactivated Mpro with much higher efficiency (kinact/KI = 4167 M−1 s−1) than the diastereomers (S,R)- and (R,S)-18 (kinact/KI = 66.0 and 20.0 M−1 s−1, respectively). (R,R)-18 exhibited both a higher kinact value and a lower KI value than those of (R,S)-18, suggesting that the CFA chirality largely affected the reaction rate of covalent modification as well as the binding affinity with Mpro. In comparison between (R,R)-18 and (S,R)-18, the binding constant KI was significantly different (1.34 and 64.7 μM, respectively), whereas the inactivation rates were similar (kinact = 0.0056 and 0.0043 s−1, respectively). This marked difference in KI value was consistent with previous studies of ML188 and Jun8-76-3A, where the R configuration at the α-pyridine chiral center was crucial for the inhibitory activities.22,23 We also determined the kinetic parameters of acrylamide 1 (Fig. S6†). Despite higher thiol reactivity of acrylamide than CFA,29 the kinact value of 1 was significantly smaller than (R)-CFA. This result implies the mismatched disposition of the acrylamide warhead in the 1–Mpro binding complex for the Michael addition of the thiol group of Cys145.
Fig. 4 Inhibition properties of the four stereoisomers of 18 against SARS-CoV-2 Mpro. (A) The absolute configurations of the stereoisomers of 18 (R = (3-fluorophenyl)ethyl). The R and S configurations at the CFA unit are highlighted in green and magenta, respectively. (B) The inhibitory activity of the stereoisomers against SARS-CoV-2 Mpro. Each plot represents the mean of triplicate experiments ± standard deviation. (C) Time-course plot of the Mpro-catalyzed hydrolysis of the fluorogenic substrate in the presence of various concentrations of (R,R)-18. [P] represents the concentration of the hydrolyzed product of the fluorogenic peptide substrate determined by the fluorescence intensity. Each plot represents the mean of triplicate experiments ± standard error of the mean. (D) Summary of IC50 values and kinetic parameters of the chiral stereoisomers of 18. ND, not determined. (E) and (F) The computationally predicted, most stable conformations of (R,R)- and (R,S)-18 in the Mpro active site. See ESI† for a detailed computational method. |
To gain further insight into the effect of CFA chirality, we performed docking simulations of the stereoisomers of 18 against Mpro. Based on the X-ray cocrystal structure of (R)-1 with SARS-CoV-2 Mpro (PDB code 7NW2),24 we examined the possible binding poses and affinities of each of the four isomers, (R,R)-, (R,S)-, (S,R)-, and (S,S)-18 in the same pocket of Mpro. The predicted binding affinities for the stereoisomers fell into the range of 2 kcal mol−1 difference, but there found no clear correlation with the experimental result. To find out the cause of the higher reactivity of (R,R)-18, we assessed the occurrences of “reactive conformations” in the predicted poses for the four isomers, where we defined the pose as reactive for the SN2 reaction (backside attack of Cys145 to the CFA carbon) if the C⋯S distance <4.5 Å and Cl–C⋯S angle >150°. We found a large number of reactive conformations (46 out of 183 poses) exclusively in the low binding energy region for (R,R)-18 among the four stereoisomers (Fig. S7†). The most stable pose of (R,R)-18 was SN2-reactive, where the Cl–C⋯S angle is nearly 180° (Fig. 4E). In the complex, the CFA fluorine and the carbonyl oxygen were close proximity to the backbone NH of Gly143, which is thought to constitute the oxyanion hole in the S1′ pocket.37 This result suggests the potential contribution of NH⋯F hydrogen bond to facilitate the reaction. The bound structure of (R,R)-18 was well aligned to the X-ray crystal structure of a dichloroacetamide inhibitor Jun9-62-2R covalently bound to Mpro, which was recently reported by Wang et al. (Fig. S8†).38 On the other hand, the most stable pose of (R,S)-18 was SN2-non-reactive due to the unpreferable orientation of the C–Cl bond for the nucleophilic attack by Cys145. The rotation of the CFA unit of (R,S)-18 could give an SN2-reactive pose, but the resulting rotamer lacks any stabilizing interaction with the Mpro active site. Taken together, we infer that (R,R)-18 is more likely to become SN2-reactive as the NH⋯F hydrogen bond forces the C–Cl bond to a favorable orientation in the binding complex. The NH⋯F hydrogen bond in the (R,R)-18–Mpro complex could also contribute to increasing the electrophilic reactivity of the CFA warhead.
The inhibitory activity of (R,R)-18 against other cysteine proteases was evaluated (Fig. 5A). (R,R)-18 potently inhibited SARS-CoV-1 Mpro with an IC50 value of 0.094 μM. This result was not surprising because the structures of Mpro are highly conserved between SARS-CoV-1 and SARS-CoV-2 (96% identity of the amino acid sequence). (R,R)-18 did not inhibit SARS-CoV-2 papain-like protease (PLpro) (<1% inhibition at 20 μM). Human cathepsin B and L were not inhibited by 20 μM of (R,R)-18, whereas the noncovalent inhibitor Jun8-76-3A weakly inhibited cathepsin B (55% inhibition at 20 μM). To assess the covalent binding properties of the CFA inhibitor, we designed a clickable activity-based probe (R,R)-22 (Fig. 5B). Introduction of the ethynyl group was well tolerated, and (R,R)-22 exhibited potent inhibitory activity against SARS-CoV-2 Mpro (IC50 = 0.11 μM) comparable to that of (R,R)-18 (IC50 = 0.056 μM). In the gel-based in vitro assay, the recombinant SARS-CoV-2 Mpro was treated with (R,R)-22 (37 °C, 1 h) and conjugated to rhodamine-azide by copper-catalyzed azide–alkyne cycloaddition (CuAAC). The in-gel fluorescence analysis showed that Mpro could be detected at concentrations as low as 10 nM by 2.5 μM of (R,R)-22, and 50 nM of Mpro was labeled by (R,R)-22 in a concentration-dependent manner (0.1–10 μM). Incubation of (R,R)-22 with A431 cell lysate spiked with Mpro (50 nM) showed the main band at ∼34 kDa. This band disappeared after preincubation with the covalent inhibitors PF-00835231 (ref. 16) and (R,R)-18 (5 μM, 30 min), but not with the reversible inhibitor Jun8-76-3A (Fig. 5C and S9†). These results demonstrate the ability of the CFA-type inhibitor to covalently modify Mpro in the complex human proteome. We also synthesized acrylamide probe 23 as an alkynylated analog of 1 and evaluated its proteome-wide activity under live cell conditions (Fig. 5D and S10†). Treatment of A431 cells with (R,R)-22 or 23 (1–10 μM, 37 °C, 1 h) and the subsequent in-gel fluorescence analysis revealed that acrylamide 23 exhibited a higher level of off-target labeling compared to 22, suggesting a low indiscriminate reactivity of the dipeptidic CFA inhibitor in live cells.
Footnote |
† Electronic supplementary information (ESI) available: Supplemental figures, experimental methods, and synthetic procedures. CCDC 2124399 and 2124400. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/d1sc06596c |
This journal is © The Royal Society of Chemistry 2022 |