Unsymmetrical aromatic disulfides as SARS-CoV-2 Mpro inhibitors: Molecular docking, molecular dynamics, and ADME scoring investigations

Graphical abstract


Introduction
The ongoing COVID-19 outbreak caused by novel coronavirus, which emerged in Wuhan by the end of December 2019, is still posing a serious problem to global health. At the time of writing this paper, the disease has caused more than 452 million affected cases and more than 6 million deaths worldwide (http://covid19.who.in). This crisis has pushed government and scientist throughout the world to intensify their efforts to develop treatments and to prevent the rapid spreading of this disease. Unfortunately, there is no available drug treatment against COVID-19 approved yet. Luckily, several vaccines candidates are now available, and WHO authorize some of them for emergency use.
SARS-CoV-2 agent is closely related to that of SARS-CoV, and this allows previously defined protein structures to be used to quickly build a drug discovery model for the novel coronavirus (Hui et al., 2020). The main protease (Mpro), which have a significant importance in the replication process of the virus, has been considered as a potential target in the development of novel SARS inhibitors (Anand et al., 2003;Holmes, 2003;Leng and Bentwich, 2003;Marra et al., 2003;Poutanen et al., 2003;Snijder et al., 2003;Thiel et al., 2003;Rota et al., 2003). Despite, MERS-CoV, SARS-CoV and SARS-CoV-2 agent belong to beta-coronaviruses category; they are in fact slightly distinguished from each other. Recent studies have shown that SARS-CoV-2 shares approximately 80% nucleotide identity and 89.10% nucleotide similarity with SARS-CoV. Therefore, the main protease of SARS-CoV, 3CLpro, has been the target of several in silico investigations to develop novel inhibitors candidates (Tahir UlQamar et al., 2020;Aouidate et al., 2020;Fouedjou et al., 2021). The 3CLpro has a high sequence identity rate between nCoV and nCoV-2; hence, their 3CLpro are likely homologous and have identical structures and functions. Moreover, SARS-CoV and SARS-CoV-2 agents infect cells in the same way and use the same protein machinery to multiply inside the host cell.
Asymmetric disulfides have become powerful resources in combinatorial chemistry and pharmaceutical research (Chtita et al.; Ameta et al.), due to their various biological activities. In fact, they exhibit anticancer activity (Branowska et al.), antibacterial activity (Kumar et al.), as well as antiviral activity by blocking the interaction between arenavirus Z protein and myeloid cell leukemia protein (Garcia et al.). These compounds are known also for their herbicidal activity (Wang et al.).
Keeping in view these features mentioned above, we propose herein to study the behavior of SARS-CoV 3CLpro inhibitors towards SARS-CoV-2 Mpro. In fact, we have performed an in silico study using docking, Dynamic, and ADME of 40 aromatic disulfides known as inhibitors of SARS-CoV 3CLpro (Wang et al., 2017;Chtita et al., 2021). The study was realized in respect of the active site of the recently crystallized SARS-CoV-2 Mpro (Protein Data Bank [PDB]: 6LU7) (Liu et al., 2020).

Material and methods
(S3 of the Supplementary Material).

Results and discussion
In this investigation, a set of 40 unsymmetrical aromatic disulfides was studied in silico to highlight their possible binding energy and interaction modes with the active site of SARS-CoV-2 Mpro, using the AutoDock Vina software. The estimated binding energy of the binding site of the 6LU7 enzyme structure is summarized in Table 1 for all the studied compounds. According to the outcomes obtained from docking screening, six compounds (6, 7, 8, 9, 10, and 11) with the strongest binding energies were chosen to depict the binding mode of the disulfide inhibitors.
Before docking molecules 6,7,8,9,10,11 into the Mpro receptor pocket, we first validated the 3D model of the Mpro protein. For this, we compared the types of residues with which the N3 inhibitor interacted in the 6LU7 complex obtained from the protein core library (https://www.rcsb.org/structure/6LU7) ( Fig. S2-A) with the residues with which the rigid structure of the N3 inhibitor interacted after it was prepared in PDBQT format ( Fig. S2-B) following the same protocol of validation reported in a study performed by Daoui et al.
From the two-dimensional visualizations shown in Fig. S2-A and S2-B, we can notice that the prepared N3 ligand interacts with the same residue that the crystallized N3 ligand interacted with. Therefore, we can consider that the model of the Mpro protein complexed with the N3 inhibitor is valid and that molecular docking on the Mpro protein can be performed in this work.
The interactions between the inhibitor and surrounding residues of SARS-CoV-2 Mpro are illustrated in the 2D and 3D schematics, which they were obtained by importing docking results into the Discovery Studio Visualizer. Fig. 1 and Table 2 show the amino acids that participate in the pattern of interactions between the ligand and enzyme with a significant partial energy contribution to total energy of interaction. Most of these interactions include Van der Waals interactions, hydrophobic contacts (Pi-Pistacked, Pi-Pi T-shaped, amide-Pistacked, alkyl Pi-Sigma, and Pi-alkyl), hydrogen bonds (hydrogen bond, carbon, and Pidonor), electrostatic, carbonyl, and specific atom-aromatic ring (Pi-cation, Pi-donor, Pi-halogen, and Pi-carbon) and provide insight into understanding molecular recognition (Du et al., 2016). Fig. 1 depict the docking interactions of the most active molecules (6, 7, 8, 9, 10, and 11) based on docking studies. Fig. 2 show Docked conformation poses of ligands 6, 7, 8, 9, 10 and 11 inside the Mpro receptor pocket and the superimposed view of these ligands.
Moreover, the results show that compound8is the most promising ligand, which bound with SARS-CoV-2 Mpro via multiple hydrogen bonding and hydrophobic contacts. The binding site is mostly located in a hydrophobic cleft lined by the following amino acids: His41, Tyr54, Phe140 Leu141, Asn142, Ser144, Cys145, His163, His164, Met165, Glu166, His172, Asp187, Arg188, Gln189, and Thr190. There are four hydrogen bond interactions with four different amino acids, Ser144, His163, Glu166, and Arg188, at distances of 2.48, 2.64, 2.19, and 3.11 Å, respectively; one acceptor-acceptor interaction with Leu141 at a distance of 2.81 Å; and two hydrophobic Pi-alkyl interactions with Met165 at 5.07 and 5.12 Å; further enquiry indicated the presence of an Pi-Pi T-shaped interaction with His41 at a distance of 5.63 Å (Fig. 3A). The high affinity of compound 8 is also associated by Van der Waals forces created on the backbone of the amide substituents with Tyr54, Phe140, Asn142, Cys145, His164, His172, Asp187, Gln189, and Thr190, which helped create a strong cohesive environment, thereby stabilizing the complex formed.
Binding ability of compound 6 is also very favorable. A hydrogen bond formed with Cys145 of the 6LU7 at a distance of 3.60 Å. The docking score for the best interaction pose is À7.20 kcal/mol. The conformational energy of the ligand is optimized by the presence of three hydrophobic interactions (Pi-alkyl interaction with Met165 at a distance of 5.25 Å, Pi-Sigma with Thr25 at a distance of 3.66 Å, and Pi-sulfur with His41 at a distance of 4.61 Å) that largely involve charge transfers to help intercalate the ligand in the binding site of the 6LU7. The high affinity of compound 6 is associated with Van der Waals forces between the backbone of the amide substituents with Thr24, Thr26, Leu27, Met49, Leu141, Asn142, Gly143, Ser144, His164, Glu166, and Gln189, which undoubtedly created a strong cohesive environment to stabilize the complex formed (Fig. 1B).
Compound 11 also has a favorable binding affinity and is associated with three H-bonds formed with the amino acids of 6LU7 (one Pi-donor hydrogen bond with His41 at a distance of 3.15 Å and two C-H interactions with Phe140 and Glu166 at distances of 3.27 and 3.55 Å, respectively). The docking score for the best interaction is À7.20 kcal/mol. The conformational energy of the ligand was minimized by the presence of five hydrophobic interactions (one Pi-sulfur with Met49 at 4.31 Å and four Pi-alkyl interactions with Cys145, His163, Met165, and His172 at distances of 5. 10, 4.74, 4.35, and 4.81 Å, respectively), which largely involve charge transfer to intercalate the ligand in the binding site of 6LU7; further enquiry indicated the presence of an Pi-Pi Tshaped interaction with His41 at a distance of 4.21 Å. The high affinity of compound 11 is associated with Van der Waals forces created on the backbone of the amide substituents with Tyr54, Leu141, Asn142, Gly143, Ser144, His164, Asp187, Arg188, and Gln189 (Fig. 1C).
The binding ability of compound 9 has an interaction score of À7.10 kcal/mol. Two hydrogen binds formed with two amino acids (one Pi-donor hydrogen bond with Glu166 at a distance of 3.11 Å and one C-H interaction with Gln189 at a distance of 3.55 Å). The conformational energy of the ligand is minimized by the presence of two Pi-alkyl interactions with Met49 and Met165 at distances of 4.84 and 4.59 Å, respectively. As with the other compounds, Van der Waals forces are present on the backbone of the amide substituents with His41, Tyr45, Phe140, Leu141, Asn142, Gly143, Ser144, Cys145, His164, His172, Asp187, and Arg188 (Fig. 1D).
The docking score for the best interaction of compound 10 is À7.10 kcal/mol, with two hydrogen bonds formed with two amino acids (one Pi-donor hydrogen bond with Glu166 at a distance of 4.19 Å and one C-H interaction with Thr190 at a distance of 3.47 Å). The energy of the ligand is minimized by the abundance of two Pi-alkyl interactions with Cys145 and Met165 at distances of 5.25 and 5.23 Å, respectively, and one Pi-Sigma interaction with His41 at a distance of 3.72 Å.
Compound 7 firmly binds at the target site of 6LU7 with three hydrogen bonds composed with three amino acids (two Pi-donor hydrogen bonds with Ser144 and Cys145 at distances of 3.03 and 3.76 Å, respectively, and one C-H interaction with Gln189 at a distance of 3.39 Å). The stability of this compound in the binding site is attributed to the three hydrophobic interactions (one Pi-alkyl with Met165 at a distance of 4.45 Å, one Pi-anion with Glu166 at a distance of 4.12 Å, and one Pi-sulfur with His41 at a distance of 4.36 Å). The high affinity of compound 7 is also associated with the presence of Van der Waals forces with Phe140, Leu141, Asn142, Gly143, His163, His164, His172, Arg188, Thr190, and Gln192 (Fig. 1F). The docking score for the best interaction pose is À7.00 kcal/mol. Furthermore, the high binding energy of ligand 8, which bears the common (1-(3-(phenyldisulfaneyl)-5-(pyridin-n-yl)-1H-1,2,4triazol-1-yl)ethan-1-one; where n = 3,4,5) structure with ligands 9, 10 and 11 can be explained by the number of hydrogen bonds formed between amino acids and radicals inserted in the group of 1(3-(phenyldisulfaneyl). We notice that the binding energy of ligand 8 (Table 2, Fig. 3A) decreases drastically when the number of hydrogen bonds between the (R-phenyl) disulfanyl group and the predicted active sites of Mpro decrease. The nitrogen group of the ligand 8 forms three regular hydrogen bonds, while the ethyl acetate and methanol groups of the ligand 9 (Fig. 1D), 10 ( Fig. 3E) and 11 (Fig. 1C) form a single hydrogen bond. Therefore, the effect of substitution of radical nitrogen in the structure of the ligand 8 by each of the ethyl acetate and methanol radicals in the structure of ligands 9, 10, and 11, respectively, could be a decisive reason to explain the observed divergence between the binding energies of ligands' structures.
On the other hand, the substitution of the (pyridin-n-yl) group in the structure of ligands 8, 9, 10 and 11 by the (phenyl) group in ligands 6 ( Fig. 1B) and 7 (Fig. 1F) did not produce any hydrogen bonding interaction with any amino acid and also the binding energies of each ligand 6 and 7 did not improve, while the binding energy of these ligands decreased from À7.2 Kcal/mol (Table 3) to À7.0 Kcal/mol (Table 3). Thus, the substitution of pyridin-n-yl with phenyl was not favorable to improve the binding energy.

Biological interpretation
The inhibitory biological activity of unsymmetrical aromatic disulfides against the SARS-CoV Mpro has been previously evaluated in vitro (Wang et al., 2017). Due to the close similarity  between the SARS-CoV and SARS-CoV-2 viruses, we have opted for investigating in silico the behavior of disulfides regarding SARS-CoV-2 Mpro. Choosing this enzyme as a target gives important leads in antiviral treatment because it participates in the proteolytic processing of replicase polyproteins; therefore, it plays a major role in viral gene expression and replication. Thus, the inhibition of this enzyme could hamper the replication of the viral genome and then the multiplication of SARS-CoV-2. In this study, we suggest that molecules that could inhibit 3CLproof SARS-CoV may result in the same interest to inhibit ofSARS-CoV-2 Mpro, because of their high sequence similarity. This hypothesis was investigated using MD analysis.
Based on the MD results, the two critical amino acids in the active site of the target are Cys145 and His41. The inhibition of SARS-CoV 3CLpro has been reported to mainly occur via an irreversible covalent bond with Cys145 and, in other several cases, by reversible interaction with Cys145 and His41, which form a catalytic dyad, where Cys145 acts a nucleophile (Ramajayam et al., 2010;Shao et al., 2008;Jacobs et al., 2013;Thanigaimalai et al., 2013a;Liu et al., 2014;Thanigaimalai et al., 2013b) and His41 acts as a base that polarizes and deprotonates the nucleophile (Cys145) to increase its reactivity; this dyad is then attacked by the inhibitor to form an intermediate complex before releasing the enzyme. In our study, compounds 6, 7, 10, and 11 display interactions with   Cys145 and His41, suggesting that these molecules could have an inhibitory activity against SARS-CoV-2 Mpro. For compounds 8 and 9, even if they form a stable complex, they will have no pharmacological effect because there is no stabilization of Cys145 or His41. SARS-CoV-2 Mpro is a cysteine protease, in which the unsymmetrical aromatic disulfides may play the role of an inhibitor, which competes with the substrate (competitive inhibitor). We believe that compounds 6, 7, 10, and 11 could be potential inhibitors for this enzyme and, consequently, may act as drug candidates against COVID-19. To confirm this assumption, an in vitro study against SARS-CoV-2 Mpro is recommended.

Pharmacokinetics and drug-likeness prediction
Pharmacokinetics and drug-likeness prediction for the six compounds were performed using SwissADME (Daina et al., 2017) to evaluate individual ADME behaviors of target compounds.
Despite the fact that chemical compounds have activity and selectivity, they may not be viable candidates for drug development. There are a number of rules and models to follow in the classification of such compound as potential drug candidate. The most widely accepted one is Lipinski's rule of five (Lipinski et al., 1997). If two parameters are out of range, ''poor absorption or permeability is possible".
The compound may be absorbed in gastrointestinal tract if any one of the parameters is not met. Hence, Lipinski's rule suggests that active molecules with good bioavailability and oral absorption have molecular weights <500, Log P < 5, less than five hydrogen bond donor (HBD) groups, and less than 10 hydrogen bond acceptor (HBA) groups (Duchowicz and Castro, 2009). A candidate that violates more than one of these criteria is less likely to be developed as a prospective oral drug.
The studied compounds are confronted to the Lipinski criteria, the results show that only four compounds listed in Table 4 meet perfectly these criteria. Thus, there would be no problem with oral bioavailability for all proposed compounds.
The Veber rule defines drug-likeness constraints as rotatable bond (RB) count 10 and polar surface area (PSA) 140 (Veber et al., 2002). PSA is strongly related to the absorption properties of active compounds. The obtained PSA for compound 8 is greater than 140, suggesting that this compound has strong polarity and is not readily absorbed by membranes. In addition, the obtained PSAs for compounds 6, 7, 9, 10, and 11 are all less than 140; which indicates that these active compounds have good oral absorption; namely good solubility and better penetration through protein membranes. The number of RBs for the six tested compounds are less than 10; as the number of RBs increases, the molecules become more flexible and more adaptable for efficient interaction with a particular binding pocket. The screening of potential candidates using Veber's rules, allowed us to select only five compounds (6, 7, 9, 10, and 11) that meet the criteria of drug-likeness assessment, whereas, compound 8 is rejected for reason of violation of one Veber criterion (Table 3).
In drug discovery schemes, oral administration and high solubility will favor complete absorption, while poor solubility limits the absorption of the drug in the gastrointestinal tract. Data listed in Table 4 suggests that compounds 6, 7, 9, 10, and 11 have good to moderate solubility, the Log S value is between À6 and À4. Thus, the afore mentioned compounds could exhibit good oral adsorption. The lower value of Log S obtained for compound 6 indicates that is hydrophilic and is more soluble (À4 < LogS < À2), and that absorption takes place in the gastrointestinal tract. The surface area and lipophilicity of this organ influences the rate of absorption of compounds. Compounds 6, 7, and 11 display high predicted gastrointestinal absorption (Table S3).
Five isoenzymes of the CYP enzyme system are mainly involved in the metabolism of most drugs. All six compounds are predicted to be inhibitors of CYP2C19 and CYP1C2. This suggested that they may be metabolized in the liver. Moreover, they are predicted not to be inhibitors of isoenzyme CYP2D6. This is very useful, because these compounds are expected to not have CYP metabolism interactions with other drugs and would not present an increased hepatotoxicity risk.
Drug oral bioavailability is the fractional extent of the drug dosage that finally reaches the therapeutic active site; an acceptable probability score of 55% indicates that the drug candidate has met the rule of five. All candidate compounds have shown a score of 55%, indicating good bioavailability.
For the parameters of the distribution, the blood-brain barrier (BBB) protects the brain from harmful substances. The capacity of a medicine to cross into the brain is a critical criterion to consider when lowering side effects and toxicities or improving the efficacy of treatments having pharmacological action in the brain. It is believed that no compounds pass through the BBB; also, the results show that there is no P-gp substrate for all the screened compounds.
When a drug is bioavailable and administered at the proper dose rate, it reaches a steady-state level. After evaluating characteristics of compound permeability through various barriers, we predicted the excretion parameters total renal clearance (log mL/ min/kg); OCT2 substrate.  The OCT2 is a protein transporter that has a vital contribution to renal absorption, disposition and clearance of drug compounds. Evaluation of the transfer of a candidate compound by OCT2 provides useful information as regards not only its clearance but also its potential contraindications. In silico analysis revealed that all compounds except 6 and 7 were not substrates of human OCT2, which is important for the excretion of cationic molecules. So, compounds 6 and 7 are predicted to undergo a renal uptake process and stay longer in the body unlike other compounds that will not have a problem with renal absorption.
The faster the excretion process, the greater the compound's CLTOT value, From Table S4 it can be seen that the log CLTOT value of the test compounds ranges from À0.151 to 0.045 ml/min/kg and from these values can be predicted that all compounds tested are relatively rapidly excreted from the body.

Molecular dynamics simulations
The MD simulation of Mpro and associated complexes was performed for 100 ns. The RMSDs of the proteins and associated complexes were calculated with respect to their initially prepared structures ( Fig. 2A). The RMSD values of each system reached equilibrium after 40 ns and then leveled out. The average RMSD of Mpro alone is 0.232 ± 0.052 nm. The RMSDs of complexes of compounds 6, 7, and 11 are 0.206 ± 0.043, 0.265 ± 0.034, and 0.197 ± 0. 028 nm, respectively. The RMSD curves of Mpro alone and all associated complexes were similar with only minor deviations. The RMSD of the entire system is less than 0.3 nm, indicting the good stability of the complexes in the aqueous environment.
The average RMSF values of backbones C a atoms of Mpro alone and associated complexes are calculated to assess dynamic behavior of amino acids of the SARS-CoV-2 Mpro. The analysis of RMSF gives useful insights regarding the structural fluctuations of different regions of the protein. Greater fluctuations in residues signify a less stable nature of the protein. The RMSF of Mpro alone and associated complexes is shown in Fig. 5B and remained below 0.2 nm for the majority of the amino acids of Mpro, indicating the stability of the ligands with the test protein. However, several fluctuations are observed at the N-and C-termini that may be due to the hanging position of these residues. Several fluctuations are also found in the certain regions of Mpro in the presence of ligands, which may be due to the dynamic behavior of ligands within the bonding region. The average RMSF of individual atoms of ligands were also calculated to check the dynamics of ligand and showed some variation, denoting the dynamical shift from their initial position.
The R g of Mpro in the absence and presence of ligands was also calculated (Fig. 3A) as this indicates the stability of protein during MD simulation. In general, a compact protein shows less variation in R g values compared with those of an expanded protein. The average R g values of Mpro alone and in complex with compounds 6, 7, or 11 were 2.188 ± 0.013, 2.198 ± 0.018, 2.207 ± 0.012, and 2. 233 ± 0.021 nm, respectively. These values signify the stability of the ligands with Mpro. The stability was further verified by calculating the changes in SASA. The SASA of Mpro and associated complexes is shown in Fig. 3B. The SASA of all complexes are similar to Mpro alone, further validating the stability of the ligands with the protein.
The interactions of ligands with Mpro were studied by calculating the hydrogen bond profiles. The number of hydrogen bonds with respect to time is shown in Fig. 4A. Hydrogen bonds are present between the ligands and Mpro over the entire simulation period. The residues of Mpro that are involved in hydrogen bond formation with compound 6areHis41, Asn142, Gly143, Ser144, Cys145, His163, and Gln189. Similarly, Thr24, Thr25, Thr26, His41, Ser46, Met49, Leu50, Asn142, Gly143, Ser144, Cys145, His163, Glu166, and Gln189 are involved in hydrogen bond formation with compound 7, and Ser144, His163, Arg188, Gln189, and Gln192 formed hydrogen bonds with compound 11; the lower number of hydrogen bonds formed with compound 11 may be due to the interaction of water molecules at the binding sites.
The effect of binding of compounds 6, 7, and 11 on the secondary structure of Mpro was also assessed. The percentages of secondary structures in Mpro and associated complexes with compound 6, 7, and 11 are shown in Fig. 4B. The percentages of coil, bsheet, b-bridge, bend, turn, and a-helix were 25.44%, 26.35%, 2.02%, 7.05%, 16.44%, and 18.30%, respectively. The interaction of ligands with Mpro resulted in insignificant changes in the secondary structures of the protein. This data further validates the structural stability of Mpro in complex with compound 6, 7, or 11.
The binding energies and forces involved in the complexation of compounds 6, 7, and 11 with Mpro were studied using the MM-PBSA analysis. The binding energy was calculated for 100 snapshots extracted from the trajectory at equal intervals of the entire MD simulation.
The different energies observed in the principal interaction of the ligands with Mpro are presented in Table 4. In typical drugprotein interactions, the major forces observed are hydrogen bonds, electrostatic forces, hydrophobic interactions and Van der Waals forces, which positively or negatively contribute to the total binding energy . The energetic contribution of the Van der Waals forces is the major contributor and most favorable in the overall binding. Electrostatic energy and SASA energy also positively contributed to the binding of all ligands with Mpro. In contrary, altered solvation energy binding of all ligands to Mpro.

Ligand-receptor interactions with human ACE2 receptor
To further support the results of molecular docking and molecular dynamics as well as ADMET properties related to the potential candidate molecules 6, 7 and 11 to inhibit SARS CoV-2 Mpro, we evaluate the ability of these molecules to interact with the active sites of the human ACE2 receptor (PDB 1R4L) (https://www.rcsb. org/structure/1R4L). This is to predict the most important interactions likely to occur in the inhibition of human ACE2 activity by the disulfide derivatives. We perform this procedure in order to predict the modes of interaction of molecules 6, 7 and 11 with human ACE2 receptors responsible for viral binding to host cells (Mohapatra et al.; Abdalla et al.) and to identify active reference sites likely to be targeted by disulfide derivatives to disrupt the interaction of SARS-CoV-2 virus with ACE2 receptors and thus inhibit entry of this virus. For this purpose, we prepared the ACE2 protein and its ligands 6, 7 and 11 with the same tools and methodology used to prepare the main protease of SARS-CoV-2 in this work. Coordinates (x = 40.605, y = 2.495, 22.919) and size (80x80x80) containing the entire ACE2 protein structure were identified as the insertion site for the 3 ligands. The summary of the obtained results by molecular docking is presented in Table 5 and the different conformations poses of the ligands in the ACE2 receptor pocket are presented in Fig. 5.
Through the results presented in Table 5, we notice that the binding energies of ligands 6, 7 and 11 with the predicted reference sites are high and close to each other.
Therefore, it is difficult to prefer one molecule among others, and this may be an appropriate reason to suggest the threescreened molecules as potential inhibitors of ACE2 as well as Mpro towards removing Coronavirus in the future.
Finally, we suggest the disulfides 6, 7 and 11 as SARS CoV-2 Mpro inhibitors because of their good interaction ability with tested protein. The compounds have also excellent bioavailability and do not pose toxicity concerns. The complexes of these compounds are also stable when simulated in aqueous medium. Therefore, based on these observations, the selected three compounds are potential candidates of SARS CoV-2 Mpro inhibitors.

Conclusion
In this work, molecular docking was performed on 40 unsymmetrical aromatic disulfides known as anti-SARS-CoV agents. Many of these compounds have shown remarkable binding interactions with SARS-CoV-2 Mpro. The outcomes have shown that the selected three disulfides (6, 7, and 11) exhibit potent ability to fight SARS-CoV-2in-silico. Also, they have shown good stability in MD simulation with SARS-CoV-2 Mpro in aqueous conditions. In addition, many of the investigated disulfides have shown good oral bioavailability and high gastrointestinal absorption. Finally, further studies should be performed in vitro and in vivo to assert if these compounds could be drug candidates to treat SARS-CoV-2. This work could constitute an in-silico approach toward the identification of novel inhibitors of anti-SARS-CoV-2.

Funding
We would like to express our gratefulness to the ''Agence Universitaire de la Francophonie (AUF)" for funding this research project (reference: AUF-DRM6588. Title: Composé dérivés de l'eugénol: hémisynthèse ciblée via la relation structure activité et investigation de l'activité anti-SARS-CoV-2). Table 5 Modes of interaction of ligands 6, 7, and 11 to ACE2 receptors and the most important predicted reference sites.