A series of nondeuterated and deuterated dipeptidyl aldehyde and masked aldehyde inhibitors that incorporate in their structure a conformationally constrained cyclohexane moiety was synthesized and found to potently inhibit severe acute respiratory syndrome coronavirus-2 3CL protease in biochemical and cell-based assays. Several of the inhibitors were also found to be nanomolar inhibitors of Middle East respiratory syndrome coronavirus 3CL protease. The corresponding latent aldehyde bisulfite adducts were found to be equipotent to the precursor aldehydes. High-resolution cocrystal structures confirmed the mechanism of action and illuminated the structural determinants involved in binding. The spatial disposition of the compounds disclosed herein provides an effective means of accessing new chemical space and optimizing pharmacological activity. The cellular permeability of the identified inhibitors and lack of cytotoxicity warrant their advancement as potential therapeutics for COVID-19.
A series of nondeuterated and deuterateddipeptidyl aldehyde and masked aldehyde inhibitors that incorporate in their structure a conformationally constrained cyclohexane moiety was synthesized and found to potently inhibit severe acute respiratory syndrome coronavirus-2 3CL protease in biochemical and cell-based assays. Several of the inhibitors were also found to be nanomolar inhibitors of Middle East respiratory syndrome coronavirus 3CL protease. The corresponding latent aldehyde bisulfite adducts were found to be equipotent to the precursor aldehydes. High-resolution cocrystal structures confirmed the mechanism of action and illuminated the structural determinants involved in binding. The spatial disposition of the compounds disclosed herein provides an effective means of accessing new chemical space and optimizing pharmacological activity. The cellular permeability of the identified inhibitors and lack of cytotoxicity warrant their advancement as potential therapeutics for COVID-19.
Coronaviruses are enveloped, positive-sense, single-stranded RNA viruses that belong to the
family Coronaviridae.[1] Among humancoronaviruses,
several strains (229E, NL63, OC43, and KHU1) are the cause of mild upper respiratory
infections; however, a few coronaviruses have emerged from animals that cause severe
respiratory disease, including severe acute respiratory syndrome coronavirus (SARS-CoV),
Middle East respiratory syndrome coronavirus (MERS-CoV), and SARS-CoV-2.[2]
Of particular concern is SARS-CoV-2, the highly pathogenic causative agent of COVID-19 which
is associated with high infectivity and is a significant threat to public health
worldwide.[3,4] The
problem is further compounded by the current lack of effective vaccines or small molecule
therapeutics for the treatment of SARS-CoV-2 infections, underscoring the urgent and dire
need for the development of prophylactic and therapeutic countermeasures to combat
infections by pathogenic coronaviruses.[5−7]The SARS-CoV-2 genome is large(∼30 kb) and similar to the genomes of SARS-CoV and
MERS-CoV (∼80 and ∼50% sequence identity, respectively). It contains two open
reading frames (ORF1a and ORF1b) and encodes multiple structural and nonstructural
proteins.[1] Translation of the genomic mRNA of ORF1a yields a
polyprotein (pp1a), while a second polyprotein (pp1ab) is the product of a ribosomal
frameshift that joins ORF1a together with ORF1b. The two polyproteins are processed by a
3C-like protease (3CLpro, also referred to as main protease, Mpro) (11 cleavage
sites) and a papain-like cysteine protease (PLpro), resulting in 16 mature nonstructural
proteins which are involved in the replication–transcription complex. The two
proteases are essential for viral replication, making them attractive targets for
therapeutic intervention.[8−15]SARS-CoV-2 3CLpro is a homodimer with a catalytic Cys–His dyad (Cys145–His41)
and an extended binding cleft. Substrate specificity profiling studies[12,13] have shown that the protease displays
a strong preference for a −Y–Z–Leu–Gln–X sequence, where
X is a small amino acid, Y is a hydrophobic amino acid, and Z is solvent-exposed and fairly
diverse (V/T/K), corresponding to the amino acid residues
−P4–P3–P2–P1–P1′-of
a substrate or inhibitor.[16] Cleavage is at the
P1–P1′ scissile bond. The 3D structure of SARS-CoV-2
3CLpro is similar to that of SARS-CoV 3CLpro; however, the S2 subsite of
SARS-CoV-2 3CLpro displays considerable plasticity and can accommodate natural and unnatural
amino acids with smaller side chains.[12] Similarly, the active-site
topography of MERS-CoV 3CLpro closely resembles that of SARS-CoV-2 3CLpro.[17] High-resolution crystal structures with bound inhibitors have been
determined, enabling the use of structure-guided approaches in the design of inhibitors. In
continuing our foray in this area,[17−19] we report
herein the results of preliminary studies related to the inhibition of SARS-CoV-2 protease
by a series of inhibitors (I) (Scheme ) that
incorporate in their structure a conformationally constrained cyclohexane moiety envisaged
to exploit new chemical space and to optimally engage in favorable binding interactions with
the active site of the protease. Furthermore, several deuterated variants of the inhibitor
were also synthesized to potentially improve the PK properties and ancillary parameters of
the inhibitors.[20−22]
Scheme 1
Synthesis of Inhibitors 2(a–o) and
3(a–o)
Results and Discussion
Inhibitor Design Rationale
The design of inhibitor (I) (Scheme ) included
the use of a P1 glutamine surrogate residue and a P2 Leu residue as
recognition elements congruent with the substrate specificity of the
protease,[12,13] as
well as an aldehyde warhead or a latent aldehyde bisulfite adduct. The design of the
inhibitor (I) was further abetted by insights gained from examining the available X-ray
crystal structures of the protease with inhibitors[12,13,17] and the results of
recent studies with cyclohexyl-derived inhibitors with demonstrated efficacy in a mouse
model of MERS-CoV-2 infection and potent inhibition against SARS-CoV-2 3CL
protease.[17]
Chemistry
The synthesis of inhibitors 2(a–o) and 3(a–o) was
readily accomplished by activating the precursor primary or secondary alcohol inputs
(Scheme ) with
N,N′-disuccinimidyl carbonate[23] and coupling the mixed carbonate with the readily accessible
Leu–Gln surrogate amino alcohol to yield alcohol product 1 which was
oxidized with Dess-Martin periodinane to generate the corresponding aldehyde (Scheme , Z = CHO 2). The aldehydes were
subsequently converted to the corresponding bisulfite adducts (Scheme
, Z = CH(OH)SO3Na 3).[24]
Scheme 2
Alcohol Inputs (a–o)
Biochemical Studies
The inhibitory activity of compounds 2–3 (a–o) against
SARS-CoV-2 3CL protease[17,25] and their activity in a cell-based system were determined as described
in the Experimental Section. The IC50 values,
EC50 values for two representative inhibitors (2a/3a), and the
CC50 values in Huh-7, CRFK, or CCL1 cells[17,25] are summarized in Table , and they are the average of at least two determinations. The
inhibitory activity of compounds 2a/3a, 2f/3f, and
2k/3k against MERS-CoV 3CL protease was also determined as described
previously,[17,18,25] and the IC50 values are listed in Table .
Table 1
IC50 and CC50 Values of SARS-CoV-2 3CLpro Inhibitors
2–3 (a–o)
compound
IC50 (μM)
CC50 (μM)
compound
IC50 (μM)
CC50 (μM)
2a
0.18 ± 0.03a
>100
3h
0.37 ± 0.03
>100
3a
0.17 ± 0.02a
>100
2i
0.27 ± 0.02
21 ± 1
2b
0.29 ± 0.05
>100
3i
0.30 ± 0.02
21 ± 3
3b
0.29 ± 0.01
>100
2j
0.31 ± 0.03
20 ± 3
2c
0.31 ± 0.06
>100
3j
0.31 ± 0.06
23 ± 1
3c
0.26 ± 0.05
>100
2k
0.18 ± 0.04
>100
2d
0.28 ± 0.08
>100
3k
0.15 ± 0.04
>100
3d
0.30 ± 0.03
>100
2l
0.29 ± 0.04
>100
2e
0.26 ± 0.05
>100
3l
0.27 ± 0.04
>100
3e
0.28 ± 0.03
>100
2m
0.27 ± 0.02
>100
2f
0.14 ± 0.02
>100
3m
0.25 ± 0.04
>100
3f
0.10 ± 0.01
>100
2n
0.82 ± 0.19
>100
2g
1.90 ± 0.14
>100
3n
1.03 ± 0.24
>100
3g
1.71 ± 0.16
>100
2o
0.74 ± 0.26
>100
2h
0.39 ± 0.02
>100
3o
0.78 ± 0.18
>100
The EC50 values for inhibitors 2a and 3a
against SARS-CoV-2 in Vero E6 cells were 0.035 ± 0.001 and 0.032 ± 0.001
μM, respectively.
Table 2
IC50 Values of MERS-CoV 3CLpro Inhibitors 2a/3a,
2f/3f, and 2k/3k
compound
IC50 (μM)
2a
0.052 ± 0.001
3a
0.049 ± 0.002
2f
0.063 ± 0.003
3f
0.058 ± 0.002
2k
0.055 ± 0.002
3k
0.053 ± 0.003
The EC50 values for inhibitors 2a and 3a
against SARS-CoV-2 in Vero E6 cells were 0.035 ± 0.001 and 0.032 ± 0.001
μM, respectively.It is clearly evident from the results shown in Table that the synthesized compounds display high potency in biochemical assays,
with most IC50 values in the sub-micromolar range. Furthermore, the inhibitors
were found to be devoid of cytotoxicity and the safety index (SI), defined as the
CC50/IC50 ratio, ranged between ∼78 and 1110. The potency
of deuterated variants 2b/3b decreased ∼1.6-fold (aldehydes) and
∼1.7-fold (bisulfite adducts) as compared to the respective nondeuterated compounds
2a/3a and remained essentially the same in the case of nondeuterated2n/3n and deuterated 2o/3o inhibitors, respectively. A change
in geometry from a cyclohexene (2e/3e) to a cyclohexane (2f/3f)
resulted in a two- to threefold increase in potency. The approximately fivefold decrease
in potency of compounds 2n/3n compared to 2k/3k presumably
reflects the inimical effect on potency of the 3° hydroxyl group. Importantly, the
EC50 values of two representative inhibitors (2a/3a) against
SARS-CoV-2 in Vero E6 cells were found to be ∼4.6-fold lower (EC50 0.035
and 0.032 μM, respectively) than the corresponding IC50 values, and the
selectivity indices of compounds 2a/3a were very high (2857 and 3125,
respectively). The significance of these findings was further augmented by the notable
inhibition of MERS-CoV 3CL protease by a select number of inhibitors (Table , compounds 2a/3a, 2f/3f, and
2k/3k), demonstrating the broad spectrum of antiviral activity displayed by
this series of compounds.Emergence of viral resistance to antiviral drugs is a major concern. We previously
reported that GC376 has a high barrier to resistance to feline infectious peritonitis
virus (FIPV) in cell culture and naturally infected animals with long-term
treatment.[18] We also examined several compounds similar to the series
in this report for emergence of viral resistance by serial passaging FIPV in the presence
of each compound in cell culture. The EC50 values of the compounds did not
increase at up to 10 passage number, and the 3CLpro of viruses passaged with each compound
has the same sequence as mock-passaged viruses. These results suggest that this series of
compounds have a high barrier to resistance.
X-ray Crystallography Studies
In order to elucidate the mechanism of action of the inhibitors and identify the
structural determinants associated with the binding of inhibitors to the active site of
SARS-CoV-2 3CL protease, high-resolution cocrystal structures were determined for
inhibitors 2a, 3b, 2f, 2k,
3c, 3d and 3e. The structure of SARS-CoV-2 3CLpro
in complex with compound 2a contained a prominent difference in electron
density consistent with the inhibitor covalently bound to Cys145 in each subunit (Figure A,B). The electron density was consistent
with the inhibitor aldehydecarbon covalently bound to the Sγ atom of the catalytic
Cys145 residue and the formation of a tetrahedral hemithioacetal, confirming the mechanism
of action. Both the R- and S-enantiomers were observed at the newly formed stereocenter,
and each enantiomer was modeled with 0.5 occupancy and was observed for all structures
described here. Interestingly, 2a adopts two conformations in which the
bicyclic ring is projected away from the S4 subsite in subunit A and is
positioned in the S4 subsite in subunit B. The inhibitor engages in multiple
favorable binding interactions with the enzyme, including direct hydrogen bond
interactions with His163 and Glu166 (γ-lactam C=O and N–H,
respectively), His41, Phe140, Ser144, His164, and Gln189 (Figure C,D). The isobutyl side chain of Leu is ensconced in the
hydrophobic S2 pocket, and the γ-lactam ring of the P1 Gln
surrogate is nestled in the S1 subsite forming hydrogen bonds with His163 and
Glu166. In addition, the lipophilic bicyclic ring in subunit A is directed toward the
surface, whereas in subunit B, it is anchored in the vicinity of the hydrophobic
S4 pocket that is lined by Ala191, Leu197, and Pro168 (Figure E,F). It should be noted that the 11 sites in the pp1a
and pp1ab polyproteins cleaved by the protease are all characterized by the presence of a
P1 Gln residue, which is conserved in all known coronavirus 3CLpro cleavage
sites. Interestingly, the deuterated analogue 3b adopts the same binding mode
and superimposes nearly identical to 2a, as shown in Figure S1. The root mean square deviation (RMSD) between the Cα atoms
of 2a and 3b was 0.27 Å for 594 residues aligned.
Figure 1
Binding mode of 2a (gray) with SARS-CoV-2 3CLpro associated with subunit
A (A,C) and subunit B (B,D). Fo–Fc polder omit map (green mesh) contoured at
3σ (A,B). Hydrogen bond interactions (dashed lines) (C,D). Surface
representations showing the orientation of 2a in subunit A (E) and
subunit B (F) near the S4 subsite of SARS-CoV-2 3CLpro with neighboring
residues colored yellow (nonpolar), cyan (polar), and white (weakly polar).
Binding mode of 2a (gray) with SARS-CoV-2 3CLpro associated with subunit
A (A,C) and subunit B (B,D). Fo–Fc polder omit map (green mesh) contoured at
3σ (A,B). Hydrogen bond interactions (dashed lines) (C,D). Surface
representations showing the orientation of 2a in subunit A (E) and
subunit B (F) near the S4 subsite of SARS-CoV-2 3CLpro with neighboring
residues colored yellow (nonpolar), cyan (polar), and white (weakly polar).Likewise, the structure of 3c shows similar binding mode properties to those
observed for 2a (Figure A,C,E). The
inhibitor engages in multiple favorable binding interactions with the enzyme, including
direct hydrogen bond interactions with His163 and Glu166 (γ-lactam C=O and
N–H, respectively), His41, Ser144, His164, and Gln189. The bicyclic ring is
oriented within the hydrophobic S4 pocket, in both subunits. Likewise, the
structure of SARS-CoV-2 3CL protease with deuterated inhibitor 3d adopts a
very similar binding mode (Figure C,D,F). The
main difference is that the electron density for 3c is the most consistent
with an axial conformation of the carbon atom attached to the bicyclic ring, whereas
3d appears to adopt an equatorial orientation (Figure S2). The structures are very similar overall, and the superposition
yielded an RMSD between the Cα atoms of 0.25 Å for 596 residues aligned.
Figure 2
Binding mode of 3c (A,B,E) and its deuterated analogue 3d
(C,D,F) with SARS-CoV-2 3CLpro. Fo–Fc polder omit map (green mesh) contoured at
3σ (A,B). Hydrogen bond interactions (dashed lines) (C,D). surface
representation showing the orientation of 3c (E) and 3d (F)
near the S4 subsite of SARS-CoV-2 3CLpro with neighboring residues colored
yellow (nonpolar), cyan (polar), and white (weakly polar).
Binding mode of 3c (A,B,E) and its deuterated analogue 3d
(C,D,F) with SARS-CoV-2 3CLpro. Fo–Fc polder omit map (green mesh) contoured at
3σ (A,B). Hydrogen bond interactions (dashed lines) (C,D). surface
representation showing the orientation of 3c (E) and 3d (F)
near the S4 subsite of SARS-CoV-2 3CLpro with neighboring residues colored
yellow (nonpolar), cyan (polar), and white (weakly polar).Similarly, inhibitors 2f, 2k, and 3e in complex
with SARS-CoV-2 were found to adopt similar binding modes in the active site of the
protease, as shown in Figures and S3. Collectively, the bicyclic rings of the inhibitors span a relatively
small region in the active site and cover a space of approximately 6.3 Å (Figure ). As such, these cocrystal structures
provide valuable insights for further structure-guided multiparameter optimization.
Figure 3
Binding mode of 2f (A,B), 2k (C,D), and 3e
(E,F) with SARS-CoV-2 3CLpro. Fo–Fc polder omit map (green mesh) contoured at
3σ (A,C,E). Hydrogen bond interactions (dashed lines) (B,D,F).
Figure 4
Superposition of all seven inhibitor-bound structures 2a (red),
3b (blue), 2f (cyan), 2k (yellow),
3c (coral), 3d (magenta), and 3e (green). The
bicyclic rings cover a space of approximately 6.3 Å in the S4
subsite.
Binding mode of 2f (A,B), 2k (C,D), and 3e
(E,F) with SARS-CoV-2 3CLpro. Fo–Fc polder omit map (green mesh) contoured at
3σ (A,C,E). Hydrogen bond interactions (dashed lines) (B,D,F).Superposition of all seven inhibitor-bound structures 2a (red),
3b (blue), 2f (cyan), 2k (yellow),
3c (coral), 3d (magenta), and 3e (green). The
bicyclic rings cover a space of approximately 6.3 Å in the S4
subsite.
Conclusions
Given the major clinical importance associated with the SARS-CoV-2 pandemic and the current
paucity of effective countermeasures, the results of the studies described herein can serve
as a launching pad for conducting further preclinical studies. Most of the compounds
exhibited high potency in biochemical assays and, for two of the compounds tested, in
cellular assays. Furthermore, members of this series were also found to potently inhibit
MERS-CoV 3CL protease, suggesting that the compounds can be developed into broad-spectrum
antivirals. Since there are no known human proteases that have a primary substrate
specificity P1 residue, that is, Gln, these inhibitors could also display high
selectivity and diminished off-target effects. Furthermore, the utilization of an aldehyde
warhead, or a latent aldehyde functionality that can rapidly generate the aldehyde
in vivo, in the design of transition state inhibitors is advantageous for
several reasons, including rapid engagement with the target, leading to the reversible
formation of a covalent adduct. The high reactivity of aldehydes is generally viewed as a
toxicity alert; however, the safety indices for most of the compounds reported herein were
found to be high. Indeed, a number of pharmaceuticals that incorporate in their structure an
aldehyde functionality are currently in clinical use and, furthermore, toxicity arising from
the presence of the aldehyde is context-specific,[26] as is presumably the
case here. Finally, the present study also sought to exploit the kinetic isotope effect
associated with the H/D bioisosteric replacement[27,28] in order to dampen oxidative metabolism at the
−CH2O– metabolic soft spot in the inhibitors,[29] as well as to reduce toxicity. Thus, the availability of equipotent deuterated analogues
that display improved PK characteristics enhances further the significance of the results
reported herein. Evaluation of a select number of inhibitors in a mouse model of SARS-CoV-2infection is in progress, and the results will be reported in due course. In conclusion, a
series of potent transition state inhibitors of SARS-CoV-2 3CL protease that incorporate in
their structures a conformationally constrained cyclohexyl moiety is reported.
Experimental Section
General
Reagents and dry solvents were purchased from various chemical suppliers (Sigma-Aldrich,
Acros Organics, Chem-Impex, TCI America, Oakwood chemical, APExBIO, Cambridge Isotopes,
Alpha Aesar, Fisher, and Advanced Chemblocks) and were used as obtained. Silica gel
(230–450 mesh) used for flash chromatography was purchased from Sorbent
Technologies (Atlanta, GA). Thin layer chromatography was performed using Analtech silica
gel plates. Visualization was accomplished using UV light and/or iodine. Nuclear magnetic
resonance (NMR) spectra were recorded in CDCl3 or dimethyl sulfoxide
(DMSO)-d6 using a Varian XL-400 spectrometer.
High-resolution mass spectrometry (HRMS) was performed at the Wichita State University
mass spectrometry laboratory using an Orbitrap Velos Pro mass spectrometer (ThermoFisher,
Waltham, MA) equipped with an electrospray ion source. The purity of all final compounds
was >95% as evidenced by NMR analysis.
Synthesis of Compounds
Preparation of Compounds 1(a–o)
General Procedure
To a solution of alcohol (1 equiv) (Scheme )
in anhydrous acetonitrile (10 mL/g alcohol) was added DSC (1.2 equiv) and TEA (3.0
equiv), and the reaction mixture was stirred for 4 h at room temperature. The solvent
was removed in vacuo, and the residue was dissolved in ethyl acetate
(40 mL/g alcohol). The organic phase was washed with saturated aqueous NaHCO3
(2 × 20 mL/g alcohol), followed by brine (20 mL/g alcohol). The organic layers were
combined and dried over anhydrous Na2SO4, filtered, and
concentrated in vacuo to yield the mixed carbonate which was used in
the next step without further purification.To a solution of Leu–Gln surrogate amino alcohol (1.0 equiv) in dry methylene
chloride (10 mL/g amino alcohol) was added TEA (1.5 equiv), and the reaction mixture was
stirred for 20 min at room temperature (solution 1). In a separate flask, the mixed
carbonate was dissolved in dry methylene chloride (10 mL/g carbonate) (solution 2).
Solution 1 was added to solution 2, and the reaction mixture was stirred for 3 h at room
temperature. Methylene chloride was added to the organic phase (40 mL/g carbonate) and
then washed with saturated aqueous NaHCO3 (2 × 20 mL/g alcohol),
followed by brine (20 mL/g alcohol). The organic phase was dried over anhydrous
Na2SO4, filtered, and concentrated in vacuo.
The resultant crude product was purified by flash chromatography (hexane/ethyl acetate)
to yield dipeptidyl alcohol 1 as a white solid.
To a solution of dipeptidyl alcohol 1 (1 equiv) in anhydrous dichloromethane
(300 mL/g dipeptidyl alcohol) kept at 0–5 °C under a N2 atmosphere
was added the DMP reagent (3.0 equiv), and the reaction mixture was stirred for 3 h at
15–20 °C. The organic phase was washed with 10% aq
Na2S2O3 (2 × 100 mL/g dipeptidyl alcohol),
followed by saturated aqueous NaHCO3 (2 × 100 mL/g dipeptidyl alcohol),
distilled water (2 × 100 mL/g dipeptidyl alcohol), and brine (100 mL/g dipeptidyl
alcohol). The organic phase was dried over anhydrous Na2SO4,
filtered, and concentrated in vacuo. The resulting crude product was
purified by flash chromatography (hexane/ethyl acetate) to yield aldehyde 2
as a white solid.
To a solution of dipeptidyl aldehyde 2 (1 equiv) in ethyl acetate (10 mL/g
dipeptidyl aldehyde) was added absolute ethanol (5 mL/g dipeptidyl aldehyde) with
stirring, followed by a solution of sodium bisulfite (1 equiv) in water (1 mL/g dipeptidyl
aldehyde). The reaction mixture was stirred for 3 h at 50 °C. The reaction mixture
was allowed to cool to room temperature and then vacuum-filtered. The solid was thoroughly
washed with absolute ethanol, and the filtrate was dried over anhydrous sodium sulfate,
filtered, and concentrated to yield a white solid. The white solid was stirred with dry
ethyl ether (3 × 10 mL/g dipeptidyl aldehyde), followed by careful removal of the
solvent using a pipette and dried using a vacuum pump for 2 h to yield dipeptidyl
bisulfite adduct 3 as a white solid.
Cloning and Expression of the 3CL Protease of SARS-CoV-2 and FRET Enzyme
Assays
The codon-optimized cDNA of full length of 3CLpro of SARS-CoV-2 (GenBank number
MN908947.3) fused with sequences encoding six histidines at the N-terminal was
synthesized by Integrated DNA (Coralville, IA). The synthesized gene was subcloned into
the pET-28a(+) vector. The expression and purification of SARS-CoV-2 3CLpro were
conducted following a standard procedure described previously.[17,18,25] Briefly, a
stock solution of an inhibitor was prepared in DMSO and diluted in assay buffer composed
of 20 mM HEPES buffer, pH 8, containing NaCl (200 mM), EDTA (0.4 mM), glycerol (60%),
and 6 mM dithiothreitol. The SARS-CoV-2 protease was mixed with serial dilutions of the
inhibitor or with DMSO in 25 μL of assay buffer and incubated at 37 °C for 1
h, followed by the addition of 25 μL of assay buffer containing the substrate
(FAM-SAVLQ/SG-QXL520, AnaSpec, Fremont, CA). The substrate was derived from the cleavage
sites on the viral polyproteins of SARS-CoV. Fluorescence readings were obtained using
an excitation wavelength of 480 nm and an emission wavelength of 520 nm on a
fluorescence microplate reader (FLx800; Biotec, Winoosk, VT) for 1 h following the
addition of the substrate. Relative fluorescence units were determined by subtracting
background values (substrate-containing well without protease) from the raw fluorescence
values, as described previously.[25] The dose-dependent FRET inhibition
curves were fitted with a variable slope using GraphPad Prism software (GraphPad, La
Jolla, CA) in order to determine the IC50 values of the compounds. The
expression and purification of the 3CLpro of MERS-CoV and the FRET enzyme assays were
performed as described previously.[17,18,25]
Cell-Based Assay for Antiviral Activity
Compounds 2a and 3a were investigated for their antiviral
activity against the replication of SARS-CoV-2. Briefly, confluent Vero E6 cells were
inoculated with SARS-CoV-2 at 50–100 plaque-forming units/well, and medium
containing various concentrations of each compound and agar was applied to the cells.
After 48–72 h, plaques in each well were counted. The 50% effective concentration
(EC50) values were determined using GraphPad Prism software using a variable
slope (GraphPad, La Jolla, CA).
Confluent cells grown in 96-well plates were incubated with various concentrations
(1–100 μM) of each compound for 72 h. Cell cytotoxicity was measured using a
CytoTox 96 nonradioactive cytotoxicity assay kit (Promega, Madison, WI), and the
CC50 values were calculated using a variable slope using GraphPad Prism
software. The in vitro SI was calculated by dividing the CC50
by the IC50.
X-ray Crystallographic Studies
Crystallization and Data Collection
Purified SARS-2 3CL protease (SARS-2 3CLpro) in 100 mM NaCl and 20 mM Tris buffer, pH
8.0, was concentrated to 9.6 mg/mL (0.28 mM) for crystallization screening. All
crystallization experiments were set up using an NT8 drop-setting robot (Formulatrix Inc.)
and UVXPO MRC (Molecular Dimensions) sitting drop vapor diffusion plates at 18 °C.
Protein (100 nL) and 100 nL of crystallization solution were dispensed and equilibrated
against 50 μL of the latter. Stock solutions of the inhibitors (100 mM) were
prepared in DMSO, and the complexes were obtained by mixing 1 μL of the ligand (2
mM) with 49 μL (0.28 mM) of SARS-2 3CLpro and incubating on ice for 1 h. Crystals
were obtained in 1–2 days from the following conditions. 2a and
3b: Berkeley screen (Rigaku Reagents) condition C5 (20% (w/v) PEG 4000, 100
mM Tris pH 8.0), 2f: Index HT screen (Hampton Research) condition H6 (20%
(w/v) PEG 3350, 200 mM sodium formate), 2k: Proplex HT screen (Molecular
Dimensions) condition D7 (15% (w/v) PEG 6000, 100 mM sodium citrate pH 5.5),
3c and 3d: the Berkeley screen (Rigaku Reagents) condition D9
(20% (w/v) PEG 3350, 100 mM Bis-Tris pH 6.5, 100 mM ammonium phosphate dibasic, 5% (v/v)
2-propanol), and 3e: Index HT screen (Hampton Research) condition C5 (15%
(w/v) PEG 3350, 100 mM succinic acid pH 7.0). Samples were transferred to cryoprotectant
solutions, prior to plunging in liquid nitrogen, composed of 80% crystallization solution
and 20% (v/v) PEG 200 except for 3c and 3d for which 20% (v/v)
ethylene glycol was used as the cryoprotectant. X-ray diffraction data were collected at
the Advanced Photon Source IMCA-CAT beamline 17-ID except for the data for the complex
with 3c which were collected at the National Synchrotron Light Source II
(NSLS-II) AMX beamline 17-ID-1.
Structure Solution and Refinement
Intensities were integrated using XDS[30,31] via Autoproc,[32] and the Laue class
analysis and data scaling were performed with Aimless.[33] Structure
solution was conducted by molecular replacement with Phaser[34] using a
previously determined structure of SARS-2 3CLpro (PDB 6XMK) as the search model. Structure refinement and manual model
building were conducted with Phenix[35] and Coot,[36]
respectively. Disordered side chains were truncated to the point for which electron
density could be observed. Structure validation was conducted with MolProbity,[37] and structure analysis/figure preparation were carried out using the
CCP4mg package.[38] Crystallographic data are provided in Table S1.[39−43]
Authors: Robert L Hoffman; Robert S Kania; Mary A Brothers; Jay F Davies; Rose A Ferre; Ketan S Gajiwala; Mingying He; Robert J Hogan; Kirk Kozminski; Lilian Y Li; Jonathan W Lockner; Jihong Lou; Michelle T Marra; Lennert J Mitchell; Brion W Murray; James A Nieman; Stephen Noell; Simon P Planken; Thomas Rowe; Kevin Ryan; George J Smith; James E Solowiej; Claire M Steppan; Barbara Taggart Journal: J Med Chem Date: 2020-10-15 Impact factor: 7.446
Authors: Clemens Vonrhein; Claus Flensburg; Peter Keller; Andrew Sharff; Oliver Smart; Wlodek Paciorek; Thomas Womack; Gérard Bricogne Journal: Acta Crystallogr D Biol Crystallogr Date: 2011-03-18
Authors: Vincent B Chen; W Bryan Arendall; Jeffrey J Headd; Daniel A Keedy; Robert M Immormino; Gary J Kapral; Laura W Murray; Jane S Richardson; David C Richardson Journal: Acta Crystallogr D Biol Crystallogr Date: 2009-12-21
Authors: Niels C Pedersen; Yunjeong Kim; Hongwei Liu; Anushka C Galasiti Kankanamalage; Chrissy Eckstrand; William C Groutas; Michael Bannasch; Juliana M Meadows; Kyeong-Ok Chang Journal: J Feline Med Surg Date: 2017-09-13 Impact factor: 2.015
Authors: Bing Bai; Elena Arutyunova; Muhammad Bashir Khan; Jimmy Lu; Michael A Joyce; Holly A Saffran; Justin A Shields; Appan Srinivas Kandadai; Alexandr Belovodskiy; Mostofa Hena; Wayne Vuong; Tess Lamer; Howard S Young; John C Vederas; D Lorne Tyrrell; M Joanne Lemieux; James A Nieman Journal: RSC Med Chem Date: 2021-08-20
Authors: Tika R Malla; Lennart Brewitz; Dorian-Gabriel Muntean; Hiba Aslam; C David Owen; Eidarus Salah; Anthony Tumber; Petra Lukacik; Claire Strain-Damerell; Halina Mikolajek; Martin A Walsh; Christopher J Schofield Journal: J Med Chem Date: 2022-05-12 Impact factor: 8.039
Authors: Chamandi S Dampalla; Athri D Rathnayake; Anushka C Galasiti Kankanamalage; Yunjeong Kim; Krishani Dinali Perera; Harry Nhat Nguyen; Matthew J Miller; Trent K Madden; Hunter R Picard; Hayden A Thurman; Maithri M Kashipathy; Lijun Liu; Kevin P Battaile; Scott Lovell; Kyeong-Ok Chang; William C Groutas Journal: J Med Chem Date: 2022-05-31 Impact factor: 8.039