Screening, Synthesis and Biochemical Characterization of SARS-CoV-2 Protease Inhibitors

The severe acute respiratory syndrome-causing coronavirus 2 (SARS-CoV-2) papain-like protease (PLpro) and main protease (Mpro) play an important role in viral replication events and are important targets for anti-coronavirus drug discovery. In search of these protease inhibitors, we screened a library of 1300 compounds using a fluorescence thermal shift assay (FTSA) and identified 53 hits that thermally stabilized or destabilized PLpro. The hit compounds structurally belonged to two classes of small molecules: thiazole derivatives and symmetrical disulfide compounds. Compound dissociation constants (Kd) were determined using an enzymatic inhibition method. Seven aromatic disulfide compounds were identified as efficient PLpro inhibitors with Kd values in the micromolar range. Two disulfides displayed six-fold higher potency for PLpro (Kd = 0.5 µM) than for Mpro. The disulfide derivatives bound covalently to both proteases, as confirmed through mass spectrometry. The identified compounds can serve as lead compounds for further chemical optimization toward anti-COVID-19 drugs.


Introduction
SARS-CoV-2 detected in China in December 2019 spread rapidly worldwide and was accompanied by high infectivity and mortality rates [1].The COVID-19 pandemic has revealed that despite advances in the medical biotechnology and pharmaceutical industry sectors, there is still a great need to study coronavirus mutation, spread and proteins to understand their action and to develop an effective treatment to stop future epidemic outbreaks [2][3][4].
Since the coronavirus infection outbreak, several vaccines have been approved, and vaccination has emerged as one of the most effective ways to prevent the spread of the coronavirus [5].However, spike protein mutations led to the emergence of new strains of the coronavirus, against which vaccines were less effective.Another strategy for treating the acute respiratory disease caused by coronavirus is the development of antiviral drugs, which could suppress virus multiplication in the human body and thus alleviate the symptoms of the disease.Numerous drug discovery projects initiated by biopharmaceutical companies and public sector organizations have emerged to find promising antiviral lead compounds, and more than 700 agents with anti-SARS-CoV-2 activities reported in preclinical or clinical studies were reviewed in [6][7][8].All of the viral proteins involved in its life cycle were possible targets for anti-coronavirus drug development.However, papain-like protease (PL pro ) and main protease (M pro , also known as 3CL pro ) have emerged as the most important targets to inhibit viral activity.
Both PL pro and M pro are cysteine hydrolases capable of cleaving polyproteins pp1a and pp1b at multiple sites to release sixteen non-structural proteins (nsp).PL pro cleaves sites between nsp1 and 2, nsp2 and 3, and nsp3 and 4, while M pro cuts the other 11 sites between the remaining nsps, resulting in a total of 16 nsps [9].These nsps assemble and form the viral replication and transcription complex on the host cell membrane.During SARS-CoV-2 infections, such cleavage of proteins led to enhanced cytokine production and the inflammatory response observed in COVID-19 patients [10,11].
The active site of SARS-CoV-2 PL pro features a catalytic triad formed by Cys111, His272, and Asp286 that cleaves the LXGG↓XX motif in viral proteins (nsp1-nsp4).Cys111 acts as the nucleophile, His272 as a general acid-base, and Asp286 favors the His alignment, thus promoting Cys111 deprotonation [12].In addition, PL pro can act as deubiquitinase.The enzyme cleaves the C-terminal LXGG sequence of ubiquitin and ubiquitin-like interferon-stimulated gene 15 (ISG15) protein leading to the suppression of the innate immune response.Thus, targeting PL pro can inhibit viral replication and promote antiviral immunity [13].
The active site of M pro consists of a catalytic dyad comprised of a nucleophilic cysteine Cys145 and nearby histidine His41 [14].M pro cleavage motif Gln↓(Ser/Ala/Asn) is highly specific for this enzyme, although some human cysteine proteases (cathepsins B, K, L, S) can also cleave after glutamine [15].The dimerization of M pro is essential for catalytic activity as the N-finger of each protomer interacts with the Glu166 of the other protomer, forming the S1 pocket of the substrate-binding site [16].
Several of the most recent reviews provide an overview of the SARS-CoV-2 protease inhibitors that have been approved or are being investigated for infection treatment, highlighting their chemical structures and binding modes [17][18][19].The known PL pro inhibitors were mainly developed based on GRL0617.This naphthalene scaffold-containing compound was originally developed as a non-covalent inhibitor of SARS-CoV PL pro , and in 2020, was found to inhibit the SARS-CoV-2 PL pro protease.The efficiency of this compound in inhibiting PL pro activity is not high (IC 50 about 2 µM) [13,20].The surface plasmon resonance (SPR) results showed weak GRL0617 binding to PL pro (K d 10.8 µM).Several new compounds with naphthyl subunits were synthesized, showing more potent binding to PL pro (K d of 2.6 µM) [21].The restricted binding pockets at substrate binding sites (Gly-Gly recognition) were a key reason for the lack of potent PL pro inhibitors.Due to their insufficient efficacy, none of the PL pro inhibitors have entered the clinical trial stage.Recently, several covalent inhibitors were discovered with a modified methyl group of GRL0617 by an electrophile capable of reacting with Cys111 of PL pro .These covalent inhibitors seem to be the most efficient inhibitors of PL pro discovered to date, with IC 50 exceeding 0.1 µM [22].
A significantly larger number of designed compounds with inhibitory properties toward M pro than PL pro have been reported, including covalent peptidomimetic and nonpeptidomimetic and non-covalent inhibitors targeting substrate binding and allosteric sites [6, [23][24][25][26].The most successful antiviral agent, Paxlovid, is the only FDA-approved drug for the treatment of high-risk patients with a confirmed SARS-CoV-2 infection.One of Paxlovid's components is nirmatrelvir, a covalent peptidomimetic compound in which a nitrile warhead forms a covalent bond with Cys145 of M pro .Nirmatrelvir demonstrates high affinity (IC 50 = 4 nM) and selectivity for M pro versus human cysteine proteases [27].Another non-covalent and non-peptidomimetic compound, S-217622, was found to be a nanomolar inhibitor of M pro (K d 13 nM) [28] and is currently in a phase 2/3 study in Japan.The treatment of patients with mild-to-moderate COVID-19 with S-217622 (ensitrelvir) demonstrated rapid and favorable antiviral efficacy and an acceptable safety profile [29].The efficacy of ensitrelvir in a wide range of patients, including older people, with COVID-19 will be further assessed in a phase 3 multinational study.
Many compounds were in silico screening as potential dual inhibitors of SARS-CoV-2 proteases.However, the dual inhibitory effect was only proven for some compounds through in vitro experiments.One compound with reactive α-chloroketone moiety was identified as a dual-acting protease inhibitor with a micromolar affinity for both proteases (IC 50 for M pro 1.72 µM, IC 50 for PL pro 0.67 µM) [30].In the study by Meewan et al. [31], several thiuram disulfide or dithiobis-(thioformate) derivatives inhibited three proteases-M pro , PL pro , and human cathepsin L-in the sub-micromolar range.
In this work, we screened a library of 1300 compounds against SARS-CoV-2 PL pro protease and tested several of the most potent hit compounds' efficacy in inhibiting PL pro enzymatic activity.The two resultant groups, the thiazole-and disulfide-containing compounds that inhibit PL pro , were tested as M pro inhibitors to probe the possibility of dual action by these compounds.

Results
A library of the available low molecular weight (MW) compounds was screened using a fluorescent thermal shift assay (FTSA) to detect their binding to a full-length recombinant protein PL pro by monitoring its melting temperature (T m ) change upon compound addition.To accurately measure the PL pro melting temperatures, we used several dyes, including ANS and Glomelt, which increase the fluorescence upon the appearance of exposed hydrophobic protein amino acids when the protein unfolds.We estimated the lowest suitable protein concentration (5 µM) and the final dye concentration (50 µM of ANS) to be optimal for compound screening.All of the PL pro melting curves demonstrated a single melting transition.Some hit compounds exhibited thermal stabilization of the protein, thus indicating compound binding.
Compounds were selected for screening from 1300 small molecules available from our in-house library.Most of the compounds were primary or secondary sulfonamides that have been synthesized in our laboratory or in collaboration with other laboratories and designed as carbonic anhydrase inhibitors.Furthermore, many resorcinol-bearing, heterocyclic compounds recently investigated as possible Hsp90, HDAC, or MMPs inhibitors were included in this library.The structures of over 550 compounds from this library and their interactions with carbonic anhydrases and Hsp90 were provided in the Protein-Ligand Binding Database (PLBD), available at https://plbd.org(accessed on 15 July 2023).
The compounds were screened at concentrations of 200 µM and 50 µM.Compounds that stabilized PL pro and induced a melting temperature shift of >1.5 • C or had a large destabilizing effect were selected for follow-up studies.Compounds that produced poorly defined melting curves were excluded from further analysis.More than 53 positives, called "hits", were identified based on this threshold value.
We tested all of the confirmed hits in dose-response biochemical inhibition experiments.The affinity (K d ) of the hit compounds for PL pro was evaluated using Z-Arg-Leu-Arg-Gly-Gly-AMC acetate as a substrate.The fluorescence of the released product, namely C-terminal AMC dye (7-amido-4-methylcoumarin), generated by PL pro proteolysis was monitored at an excitation wavelength of 340 nm and an emission wavelength of 450 nm.Compounds that did not show inhibition at the highest used concentration of 200 µM were defined as inactive (K d ≥ 200 µM).Hits that inhibited less than 50% of the enzymatic activity at 50 µM of the compound were not further investigated in the compound dose-response assay due to weak potency.However, we determined the exact K d values for several such compounds to elucidate the structure-activity relationship (SAR) of the inhibitors and to identify the functional groups that govern enhanced affinity.Structurally similar hits were grouped into compounds containing a naphthalene ring (Figure 1A) and ligands containing a disulfide bond (Figure 1B).Structural similarities between the compounds in each group guided the SAR studies.
(Figure 1A) and ligands containing a disulfide bond (Figure 1B).Structural similarities between the compounds in each group guided the SAR studies.

Naphthalene-Based PL pro and M pro Inhibitors
Based on the modifications of 3-[naphthalene-1-yl(1,3-thiazol-2-yl)amino]propanoic acid (compound 3), the compounds were grouped into compounds bearing (Figure 2):    (Figure 1A) and ligands containing a disulfide bond (Figure 1B).Structural similarities between the compounds in each group guided the SAR studies.Nineteen compounds inhibited PL pro activity with a K d value of ≤50 µM (Table 1, Figures S19 and S20).We performed a detailed structure-activity correlation analysis for this group of compounds and found that the parent compound, 3-(naphtylamino]propanoic acid (compound 1), weakly inhibited the PL pro activity (K d = 68 µM).The attachment of the thioamide functional group to the nitrogen atom (compound 2) led to a 2-fold increase in the binding affinity (K d = 29 µM), while the introduction of a thiazole ring (compound 3) enhanced the affinity by 5-fold (K d = 14 µM) compared to the starting compound 1 (Figure 3).Examining the binding of the first group of compounds modified in the 4-position of the thiazole ring, it is evident that all of the modifications decreased the binding affinity.For example, the addition of the methyl group (compound 4) weakened the interaction by 8-fold, while carbonyl (compound 5) resulted in an inactive compound (K d > 200 µM).Only one compound, 6e, with a 4-cyanophenyl substituent at the 4-position, bound PL pro with similar affinity to the parent compound 3 (K d 20 µM and 14 µM, respectively).Examining the binding of the first group of compounds modified in the 4-position of the thiazole ring, it is evident that all of the modifications decreased the binding affinity.For example, the addition of the methyl group (compound 4) weakened the interaction by 8-fold, while carbonyl (compound 5) resulted in an inactive compound (Kd > 200 µM).Only one compound, 6e, with a 4-cyanophenyl substituent at the 4-position, bound PL pro with similar affinity to the parent compound 3 (Kd 20 µM and 14 µM, respectively).
We next determined whether our PL pro inhibitors inhibited the M pro activity using ([5-FAM]-AVLQSGFR-[Lys(Dabcyl)]-K-amide as a substrate.Interestingly, six thiazole derivatives were inhibitors of M pro , with 13a and 13c being the most potent (K d 6.5 µM and 8.6 µM, respectively) (Table 1).The K d values indicated that 13a was 3.7-fold more potent for M pro inhibition than PL pro , and 13c bound both proteases with the same affinity (K d s are 8.6 µM, Figure 4).
Overall, the SAR results emphasized the thiazole derivatives as a promising scaffold for the further development of dual PL pro and M pro inhibitors.

Disulfide Derivatives as PL pro and M pro Inhibitors
We found seven benzene disulfide compounds with varied functional groups (Figure 6) that were efficient PL pro inhibitors, with K d values between 0.43 and 2.5 µM (Table 2).The highest affinities (K d near 0.5 µM) were determined for compounds 32, 34, and 37.Moreover, two compounds-32 and 34-inhibited M pro with K d values of 2.8 and 3.5 µM, respectively.Compound 34 was eight times more potent for PL pro than M pro (Figure 7).

Disulfide Derivatives as PL and M Inhibitors
We found seven benzene disulfide compounds with varied functional group ure 6) that were efficient PL pro inhibitors, with Kd values between 0.43 and 2.5 µM 2).The highest affinities (Kd near 0.5 µM) were determined for compounds 32, 34, a Moreover, two compounds-32 and 34-inhibited M pro with Kd values of 2.8 and 3 respectively.Compound 34 was eight times more potent for PL pro than M pro (Figure Figure 6.Chemical structures of disulfide compounds discovered as PL pro inhibitors.The sy of the compounds and their characterization have previously been described in [32].Table 2.The dissociation constants K d (µM) for disulfide compounds' interaction with PL pro and M pro were determined using the enzymatic inhibition assay at 25 It was recently confirmed through native mass spectrometry that disulfiram and tetraethylthiuram disulfide interact covalently with cysteine residues of M pro and the inhibition potency highly depends on the reducing agents, such as DTT or glutathione [33].However, Wang et al. [34] synthesized unsymmetrical aromatic disulfide derivatives, inhibiting SARS-CoV M pro through the reversible binding mechanism.We performed mass spectrometry experiments to elucidate the covalent binding mechanism of disulfide compounds with the proteases in the absence of reducing agents.The covalent binding of 32 with PL pro and M pro was confirmed through MS experiments and the data are shown in Figure 8. Compound 32 (MW (32) = 248.37)incubated with PL pro showed a major peak relating to protein-compound adduct, with a mass shift equal to the compound molecular mass of half moiety of the disulfide (∆MW = 38,507.79− 38,384.66= 123.13)and a minor set of peaks relating to two or three inhibitor molecules bound to the protein.The MS data also showed the covalent binding of compound 34 to PLpro (Figure S21).The molecular weight of PL pro changes upon the addition of ligands 32 and 34 by approximately 8.6 Da (from 38,393.27 to 38,384.66), it is most likely that with this addition, modification occurs in the presence of disulfide compounds.Table 2.The dissociation constants Kd (µM) for disulfide compounds' interaction with PL pro and M pro were determined using the enzymatic inhibition assay at 25 °C and pH 7.5.1.9 >200

Compound
It was recently confirmed through native mass spectrometry that disulfiram and tetraethylthiuram disulfide interact covalently with cysteine residues of M pro and the inhibition potency highly depends on the reducing agents, such as DTT or glutathione [33].However, Wang et al. [34] synthesized unsymmetrical aromatic disulfide derivatives, inhibiting SARS-CoV M pro through the reversible binding mechanism.We performed mass spectrometry experiments to elucidate the covalent binding mechanism of disulfide compounds with the proteases in the absence of reducing agents.The covalent binding of 32 with PL pro and M pro was confirmed through MS experiments and the data are shown in Figure 8. Compound 32 (MW (32) = 248.37)incubated with PL pro showed a major peak relating to protein-compound adduct, with a mass shift equal to the compound molecular mass of half moiety of the disulfide (ΔMW = 38,507.79− 38,384.66= 123.13)and a minor set of peaks relating to two or three inhibitor molecules bound to the protein.The MS data also showed the covalent binding of compound 34 to PLpro (Figure S21).The The molecular weight of M pro after incubation with compound 32 had a major shift corresponding to two bound inhibitor molecules with their half moieties (ΔMW = 34,040.75-33,795.02= 245.7).Also, minor peaks relating to one and three inhibitor molecules bound to the protein were detected.The MS results confirmed that symmetrical disulfide compounds can modify the catalytic cysteine (C111 of PL pro and C145 of M pro ) and bind to other cysteine residues in the proteases.To confirm the covalent nature of the disulfide compound binding to PL pro , we produced mutant PL pro C111S and performed FTSA experiments with disulfide 32 and GRL0617.The binding affinities of wt PL pro and the mutant PL pro C111S proteins for GRL0617, a well-known noncovalent inhibitor, were similar (Kd 1.2 and 2.5 µM, respectively).The addition of the disulfide compound 32 strongly destabilized the native PL pro , but not the mutant PL pro C111S (Figure 9A,B).This shows that the binding of disulfide 32 The molecular weight of M pro after incubation with compound 32 had a major shift corresponding to two bound inhibitor molecules with their half moieties (∆MW = 34,040.75− 33,795.02= 245.7).Also, minor peaks relating to one and three inhibitor molecules bound to the protein were detected.The MS results confirmed that symmetrical disulfide compounds can modify the catalytic cysteine (C111 of PL pro and C145 of M pro ) and bind to other cysteine residues in the proteases.
To confirm the covalent nature of the disulfide compound binding to PL pro , we produced mutant PL pro C111S and performed FTSA experiments with disulfide 32 and GRL0617.The binding affinities of wt PL pro and the mutant PL pro C111S proteins for GRL0617, a well-known noncovalent inhibitor, were similar (K d 1.2 and 2.5 µM, respectively).The addition of the disulfide compound 32 strongly destabilized the native PL pro , but not the mutant PL pro C111S (Figure 9A,B).This shows that the binding of disulfide 32 involves catalytic cysteine and the covalent modification of cysteine would lead to protein degradation.

Discussion
In the search for compounds that would exhibit inhibitory properties against any enzymes participating in any disease, a library of compounds is often screened to discover whether any of the compounds bind with sufficient affinity.Various assays are used for screening and may be based on recombinantly prepared target proteins or cells containing the proteins.Both inhibition and binding assays may be employed.Here, we chose the fluorescence-based thermal shift assay (FTSA, also known as differential scanning fluorimetry, DSF) as it requires low protein consumption and simple PCR equipment that allows for fluorescence measurements over increased temperatures.Even weakly binding ligands affect protein stability, and based on a shift in the melting temperature (Tm) of the target protein, hit compounds can be identified.Additional molecular interactions between ligand and protein usually stabilizes the protein, thus reducing the Gibbs energy [36,37].Compounds that bind more strongly will exhibit greater ther-

Discussion
In the search for compounds that would exhibit inhibitory properties against any enzymes participating in any disease, a library of compounds is often screened to discover whether any of the compounds bind with sufficient affinity.Various assays are used for screening and may be based on recombinantly prepared target proteins or cells containing the proteins.Both inhibition and binding assays may be employed.Here, we chose the fluorescence-based thermal shift assay (FTSA, also known as differential scanning fluorimetry, DSF) as it requires low protein consumption and simple PCR equipment that allows for fluorescence measurements over increased temperatures.Even weakly binding ligands affect protein stability, and based on a shift in the melting temperature (T m ) of the target protein, hit compounds can be identified.Additional molecular interactions between ligand and protein usually stabilizes the protein, thus reducing the Gibbs energy [36,37].Compounds that bind more strongly will exhibit greater thermal stabilization of the proteins than weakly binding compounds.
However, in addition to thermal stabilization, the same target protein may exhibit thermal destabilization in the presence of some ligands, shifting the T m downward.It is thought that such ligands bind more strongly to the unfolded than to the native protein state.However, the destabilizing effect is still poorly understood [38].Compounds that destabilize a protein may lead to more rapid protein degradation.Therefore, such destabilizers are also worth investigating [39].Our studies confirmed that protein destabilizers, such as symmetrical disulfide compounds, could be potent SARS-CoV-2 PL pro and M pro inhibitors.However, their use as SARS-CoV-2 drugs could be limited due to the reactive nature of such compounds.The disulfide compounds not only modify amino acids at the active site, but can most likely bind to nucleophylic amino acids localized in other sites of the protein.
The structure-activity relationship of the naphthalene-bearing thiazole derivatives revealed that the naphthalene group is not a necessary structural element for PL pro inhibition.It was previously shown through the X-ray crystallography of GRL0617 and its derivatives in complex with PL pro that the naphthyl ring occupies the hydrophobic pocket of the enzyme and forms hydrophobic interactions with the aromatic amino acids Y264 and Y268 [40].The replacement of naphthyl by the benzylidene group had similar effects to the binding affinity.Compounds such as 13c with benzylidene moiety could inhibit PL pro with K d 8.6 µM, while the best naphthalene-group bearing compound, 22, inhibited PL pro with K d 5.2 µM.It should be emphasized that the binding affinity to both proteases was highly dependent on the nature, length, and hydrophobicity of the substituents on the thiazole ring.The identified hit compounds were quite varied in their structures; therefore, it was difficult to perform detailed SAR analysis.We believe thiazole-bearing compounds can be optimized further by modifying the carboxy-ethyl group.In our library of compounds, there were no modifications at this site, but further optimization at the 2'-position on the thiazole ring might improve the affinity of compounds for both SARS-CoV-2 proteases.

Chemistry
Commercially available solvents and reagents were used without further purification unless otherwise mentioned.Reagents and solvents were obtained from Sigma-Aldrich (St. Louis, MO, USA) and used without further purification.The reaction course and purity of the synthesized compounds were monitored through TLC using aluminium plates precoated with silica gel with F 254 nm (Merck KGaA, Darmstadt, Germany).Melting points were determined with a B-540 melting point analyzer (BüchiCorporation, New Castle, DE, USA) and were uncorrected.NMR spectra were recorded on a Brucker Avance III (400, 101 MHz) spectrometer.Chemical shifts were reported in (δ) ppm relative to tetramethyl silane (TMS) with the residual solvent as an internal reference (DMSO-d 6 , δ = 2.50 ppm for 1 H and δ = 39.5 ppm for 13 C).Data were reported as follows: chemical shift, multiplicity, coupling constant (Hz), integration, and assignment.IR spectra (ν, cm −1 ) were recorded on a Perkin-Elmer Spectrum BX FT-IR spectrometer (Perkin-Elmer Inc., Waltham, MA, USA) using KBr pellets.Elemental analyses (C, H, N) were conducted using the Elemental Analyzer CE-440 (Exeter Analytical, Inc., Chelmsford, MA, USA); their results were found to be in good agreement (±0.3%) with the calculated values.
Synthetic details for the preparation of compounds 1-38 and the characterization data of these compounds are given in the Supplementary Material (Figures S1-S18).
Site-directed mutagenesis of plasmid encoding full-length PL pro was performed with primers (Supplementary Material Table S1) using the standard PCR procedure and verified through sequencing analysis (Supplementary Material).Expression and purification conditions for the PL pro C111S mutant were the same as for the native protein.
Recombinant full-length M pro protein was purchased from Sigma-Aldrich.

Fluorescence Thermal Shift Assay (FTSA)
FTSA experiments were performed with the real-time PCR instrument "QIAGEN Rotor-Gene Q".The PL pro solution in the absence and presence of compounds was heated from 25 • C to 99 • C through increasing the temperature by 1 • C per minute.We optimized the experimental conditions for the screening of compounds: assay buffer, PL pro concentration, and fluorescence dye (1-anilino-8-naphthalenesulfonate (ANS) or Glomelt).The fluorescence excitation wavelengths for ANS and Glomelt were set at 365 or 468 nm and the detection wavelengths were 460 or 507 nm, respectively.The optimal PL pro concentration was 5 µM and ANS concentration 50 µM.In addition, 10 mM stock solutions of 1300 compounds from our in-house library were made in DMSO, then diluted in 50 mM HEPES containing 100 mM NaCl at pH 7.5, with the final compound concentration of 50 and 200 µM.Protein solutions were prepared in 5 mM HEPES (pH 7.5) containing 10 mM NaCl.Protein T m values were obtained by fitting the melting curves to a two-state model.
To evaluate the binding affinity, 2-fold serial dilutions of the ligand stock in DMSO were performed, and then the prepared samples were 50-fold diluted with the buffer.Data analysis and compound dissociation constants were determined as previously described [35].

Enzymatic Inhibition Assay
The protease enzymatic activity of SARS-CoV-2 PL pro was measured using a FRETbased enzymatic activity assay.Experiments were carried out in non-binding 96-well plates at 25 • C, containing a final reaction volume of 100 µL.Ligand solutions (200-0 µM) were prepared in 50 mM HEPES buffer containing 100 mM NaCl, 2% (v/v) DMSO, pH 7.5, and mixed with PL pro (final concentration 40 nM).After 30 min of incubation, peptidomimetic substrate Z RLRGG-AMC (250 µM) was added to initiate the reaction.Fluorescence was monitored for 30 min on a BioTek Synergy H4 Hybrid plate reader (λex = 340 nm; λem = 450 nm).Enzyme activity was equalized to the slope of the reaction progress linear curve and relative activity was calculated from the slopes' ratio of inhibited and control (without ligand) samples.Three independent experiments with each ligand (except for low potency compounds, K d > 50 µM) were used to calculate mean relative enzyme activity.
Mean relative activity values were applied for the evaluation of binding affinities using the Morrison Equation (1).Data precision was evaluated using standard deviation.
Enzyme activity (%) = Inhibition experiments of compounds 32-38 were performed with PL pro dialyzed into 20 mM HEPES buffer containing 50 mM NaCl, 5% glycerol, pH 7.5 to avoid DTT impact for disulfide compounds stability disruption.

Mass Spectrometry
Mass spectrometry experiments were performed with an electrospray ionization timeof-flight mass spectrometer (Q-TOF) to show the covalent binding of compound 32 to proteases.The 0.1 mg/mL protease solution was prepared in the absence or presence of 32 (1:5 PL pro : compound molar ratio and 1:10 M pro : compound molar ratio) and the solution was incubated for 1 h at room temperature before analysis.The final DMSO concentration was 1% (v/v).

Conclusions
A new class of inhibitors of SARS-causing coronavirus CoV-2 protease PL pro were discovered using the thermal shift assay.From a library of 1300 available compounds, 53 hits were selected for the in vitro enzyme inhibition assay of PL pro .Several thiazolebearing compounds inhibited PL pro with micromolar affinity.Further optimization at the 2 -position on the thiazole ring is required to improve the activity of the compounds as SARS-CoV-2 proteases inhibitors.Several symmetrical disulfide compounds with a dissociation constant in the sub-micromolar range (K d = 0.5 µM) were identified.Interestingly, two disulfides also demonstrated promising M pro inhibition, with K d values of 3.5 µM.The mass spectrometry data revealed the covalent binding mechanism of disulfides toward both proteases.The identified compounds with thiazole or aromatic disulfide scaffolds can serve as lead compounds for anti-coronavirus drug development.

Figure 1 .
Figure 1.The discovered PL pro hits: naphthalene-based (A) and disulfide bond-bearing compounds (B).Structure of GRL0617, one of the first compounds discovered as a selective SARS-CoV-2 PL pro inhibitor (C).

Figure 2 .
Figure 2. Chemical structures of thiazole derivatives.Compound synthesis is described in the Supplementary Material.

Figure 1 .
Figure 1.The discovered PL pro hits: naphthalene-based (A) and disulfide bond-bearing compounds (B).Structure of GRL0617, one of the first compounds discovered as a selective SARS-CoV-2 PL pro inhibitor (C).

Figure 2 .
Figure 2. Chemical structures of thiazole derivatives.Compound synthesis is described in the Supplementary Material.

Figure 2 .
Figure 2. Chemical structures of thiazole derivatives.Compound synthesis is described in the Supplementary Material.

Figure 3 .
Figure 3. Compound binding affinity for PL pro .The addition of thiazole group (1-3) increased the binding affinity for PL pro by approximately 5-fold.

Figure 3 .
Figure 3. Compound binding affinity for PL pro .The addition of thiazole group (1-3) increased the binding affinity for PL pro by approximately 5-fold.

Figure 4 .
Figure 4. Inhibition of PL pro (A,B) and M pro (C,D) activity by compound 13c.(A,C) Time-dependent product formation in the absence and presence of inhibitor.(B,D) Dose-response curve of PL pro and M pro inhibition with 13c obtained after fitting curves in panels (A,C).

Figure 4 .
Figure 4. Inhibition of PL pro (A,B) and M pro (C,D) activity by compound 13c.(A,C) Time-dependent product formation in the absence and presence of inhibitor.(B,D) Dose-response curve of PL pro and M pro inhibition with 13c obtained after fitting curves in panels (A,C).

Figure 4 .
Figure 4. Inhibition of PL pro (A,B) and M pro (C,D) activity by compound 13c.(A,C) Time-dependent product formation in the absence and presence of inhibitor.(B,D) Dose-response curve of PL pro and M pro inhibition with 13c obtained after fitting curves in panels (A,C).

Figure 6 .
Figure 6.Chemical structures of disulfide compounds discovered as PL pro inhibitors.The synthesis of the compounds and their characterization have previously been described in[32].

Figure 7 .
Figure 7. Inhibition of PL pro (A,B) and M pro (C,D) activity by compound 34.(A,C) Time-dependent product formation in the absence and presence of inhibitor.(B,D) Dose-response curve of PL pro and M pro inhibition with 34 obtained after fitting curves in panels (A,C).

Figure 7 .
Figure 7. Inhibition of PL pro (A,B) and M pro (C,D) activity by compound 34.(A,C) Time-dependent product formation in the absence and presence of inhibitor.(B,D) Dose-response curve of PL pro and M pro inhibition with 34 obtained after fitting curves in panels (A,C).

Figure 8 .
Figure 8. Mass spectra of native PL pro , M pro and protease-compound complexes.(A)-MS spectra of native PL pro corresponds to the molecular mass of the protein equal to 38,393.27.(B)-MS spectra of PL pro incubated with compound 32 (protein/ligand molar ratio 1:5).The difference compared to the molecular weight of PL pro is equal to 123.13 Da. (C)-MS spectra of native M pro corresponds to the molecular mass of the protein equal to 33,797.54.(D)-spectra of M pro incubated with 32 (protein/ligand molar ratio 1:10).The difference in molecular weight of M pro -32 complex (major peak) compared to the molecular weight of native M pro is equal to 245.73 Da.

Figure 8 .
Figure 8. Mass spectra of native PL pro , M pro and protease-compound complexes.(A)-MS spectra of native PL pro corresponds to the molecular mass of the protein equal to 38,393.27.(B)-MS spectra of PL pro incubated with compound 32 (protein/ligand molar ratio 1:5).The difference compared to the molecular weight of PL pro is equal to 123.13 Da. (C)-MS spectra of native M pro corresponds to the molecular mass of the protein equal to 33,797.54.(D)-spectra of M pro incubated with 32 (protein/ligand molar ratio 1:10).The difference in molecular weight of M pro -32 complex (major peak) compared to the molecular weight of native M pro is equal to 245.73 Da.
Int. J. Mol.Sci.2023, 24, x FOR PEER REVIEW 11 of 16involves catalytic cysteine and the covalent modification of cysteine would lead to protein degradation.

Figure 9 .
Figure 9. Compound binding measured by FTSA.Raw data of 32 binding to wild-type PL pro (A)and mutant PL pro C111S (B).The thermal destabilization effect (Tm decrease) induced by 32 is seen only for wild-type PL pro (C), indicating that the destabilization is attributed to the direct binding of the ligand to the Cys111.(D) Dependence of Tm of wild-type and mutant PL pro on concentrations of the added GRL0617 (red squares-wt PL pro , blue squares-PL pro C111S).Dosing curves were fitted according to[35].

Figure 9 .
Figure 9. Compound binding measured by FTSA.Raw data of 32 binding to wild-type PL pro (A) and mutant PL pro C111S (B).The thermal destabilization effect (T m decrease) induced by 32 is seen only for wild-type PL pro (C), indicating that the destabilization is attributed to the direct binding of the ligand to the Cys111.(D) Dependence of T m of wild-type and mutant PL pro on concentrations of the added GRL0617 (red squares-wt PL pro , blue squares-PL pro C111S).Dosing curves were fitted according to [35].

Table 1 .
The dissociation constants K d (µM) for compound interaction with PL pro and M pro were determined by enzymatic inhibition assay at 25 • C and pH 7.5.
>200 Uncertainty of measurement is approximately 1.8-fold of the K d .
is approximately 1.8-fold of the Kd.