key: cord-0982082-c7akr7g4
authors: Bui, Thanh Q.; Phuong Loan, Huynh Thi; Ai My, Tran Thi; Quang, Duong Tuan; Phuong Thuy, Bui Thi; Nhan, Vo Duy; Quy, Phan Tu; Van Tat, Pham; Dao, Duy Quang; Trung, Nguyen Tien; Huynh, Lam K.; Ai Nhung, Nguyen Thi
title: A density functional theory study on silver and bis-silver complexes with lighter tetrylene: are silver and bis-silver carbenes candidates for SARS-CoV-2 inhibition? Insight from molecular docking simulation
date: 2020-08-21
journal: RSC advances
DOI: 10.1039/d0ra05159d
sha: 46ba001b9b37cd904f3adb90d404d267af98f698
doc_id: 982082
cord_uid: c7akr7g4

Ribavirin and remdesivir have been preclinically reported as potential drugs for the treatment of SARS-CoV-2 infection, while light silver tetrylene complexes (NHE(Ph)–AgCl and (NHE(Ph)–AgCl)(2) with E = C, Si, and Ge) have gained significant interest due to their promising applicability on the cytological scale. Firstly, the structures and bonding states of silver–tetrylene complexes (NHE–Ag) and bis-silver–tetrylene complexes (NHE–Ag-bis) were investigated using density functional theory (DFT) at the BP86 level with the def2-SVP and def2-TZVPP basis sets. Secondly, the inhibitory capabilities of the carbene complexes (NHC–Ag and NHC–Ag-bis) and the two potential drugs (ribavirin and remdesivir) on human-protein ACE2 and SARS-CoV-2 protease PDB6LU7 were evaluated using molecular docking simulation. The carbene ligand NHC bonds in a head-on configuration with AgCl and (AgCl)(2), whereas, the other NHE (E = Si and Ge) tetrylene ligands bond in a side-on mode to the metal fragments. The bond dissociation energy (BDE) of the NHE–Ag bond in the complex families follows the order of NHC–Ag > NHSi–Ag > NHGe–Ag and NHSi–Ag-bis > NHGe–Ag-bis > NHC–Ag-bis. The natural bond orbital analysis implies that the [NHE(Ph)→AgCl] and [(NHE(Ph))(2)→(AgCl)(2)] donations are derived mainly from the σ- and π-contributions of the ligands. The docking results indicate that both the ACE2 and PDB6LU7 proteins are strongly inhibited by silver–carbene NHC–Ag, bis-silver–carbene NHC–Ag-bis, ribavirin, and remdesivir with the docking score energy values varying from −17.5 to −16.5 kcal mol(−1) and −16.9 to −16.6 kcal mol(−1), respectively. The root-mean-square deviation values were recorded to be less than 2 Å in all the calculated systems. Thus, the present study suggests that silver–carbene NHC–Ag and bis-silver–carbene NHC–Ag-bis complexes are potential candidates to inhibit ACE2 and PDB6LU7, and thus potentially conducive to prevent infection caused by the SARS-CoV-2 virus.

The interest in divalent carbon compounds has been revitalized by the discovery of highly stable nucleophilic N-heterocyclic carbenes (NHCs) by Arduengo in 1991. 1,2 Due to their s-donor capabilities, ease of preparation, functionalization, 3, 4 and strong metal-ligand interactions that allow stable multimetallic arrangements, [5] [6] [7] NHCs constitute a promising alternative to tertiary phosphines. Arduengo synthesized the rst NHC-Ag(I) complex before applying the free carbene method two years later to successfully isolate free NHC. 8 Besides, silver-NHC complexes [9] [10] [11] have been reported to possess both antibiotic [12] [13] [14] and anticancer properties [15] [16] [17] [18] and have untapped potential as drug candidates. Youngs et al. (2004) reported the rst use of silver NHCs as antimicrobial agents. 19 Also, recent studies have proposed a structureanticancer activity correlation for NHC-silver(I) complexes. 20, 21 In particular, 1-methyl-3-(p-cyanobenzyl)-benzimidazol-2ylidene silver acetate, 21 methyl caffeine-derived silver acetate 22, 23 and their structural analogues have been studied. The results implied that there is a signicant improvement in biological applicability, especially in clinical therapy, with the formation of metal-ligand coordination. 24, 25 Carbenes are electrophiles in their free state and become nucleophiles in an N-heterocyclic-coordinated system. Nheterocyclic carbenes (NHCs) are a well-known class of ligands with the capacity to form complexes with most maingroup metals 26 and transition metals, 27 including rare earth elements. 28 NHC-Ag(I) complexes have gained signicant attention among the family of NHC metal complexes considering their synthesizability and applicability. In a previous study, two silver(I) complexes, silver(I)-2,6-bis(ethanolimidazolemethyl)pyridine hydroxide and silver(I)-2,6-bis(propanolimidazolemethyl)pyridine hydroxide were shown to exhibit stronger antimicrobial activity than AgNO 3 against Escherichia coli, Staphylococcus aureus and Pseudomonas aeruginosa. 19 Another important contribution was from Ghosh's research group for the synthesis and antimicrobial evaluation of NHCsilver complexes derived from 1-benzyl-3-tert-butylimidazole. 17 Recently, Gurbuz and colleagues showed that new imidazolidin-2-ylidene silver complexes can exhibit effective antimicrobial activity against a variety of bacteria and fungi. 29 Coronaviruses (CoVs) are positive-sense, single-stranded RNA viruses infecting a wide range of animal hosts. 30 A new virus rst emerged in Wuhan, China, which has presently caused a global pandemic. It has been identied as a novel coronavirus (SARS-CoV-2), and is closely related to the severe acute respiratory syndrome from 2003 (SARS-CoV). 31 Currently, no specic treatments have been found to be effective against this new virus. Thus, a practical approach is to determine whether or not existing antiviral drugs are reasonably effective for the treatment of SARS-CoV-2 infection. Several drugs, such as ribavirin, interferon, lopinavir-ritonavir, and corticosteroids, have been administered to patients with SARS or MERS, although the total efficacy of some drugs remains controversial. 32 Wu Zhong et al. evaluated the antiviral effectiveness of ve FAD-approved drugs, including ribavirin, penciclovir, nitazoxanide, nafamostat, chloroquine and two well-known broadspectrum antiviral drugs, remdesivir (GS-5734) and favipiravir (T-705), against clinically isolated SARS-CoV-2 via in vitro experiments. 33 In fact, ribavirin can be a potential drug for the suppression of 2019 SARS-CoV-2 considering its clinical efficacy demonstrated in 2003 SARS-CoV and 2012 MERS-CoV. 34 Ribavirin (i.e. 1-b-D-ribofuranosyl-1,2,4-triazole-3carboxamide) is a guanosine analogue with a broad spectrum of antiviral activity against RNA and DNA viruses. This drug was patented in 1971 and approved for medical use in 1986. 35 It is on the World Health Organization's list of essential medicines, the safest and most effective medicines needed in the healthcare system. 36 In the treatment of hepatitis C, ribavirin is used in combination with other medications such as simeprevir, sofosbuvir, peginterferon alfa-2b and peginterferon alfa-2a. When administrated for the treatment of viral hemorrhagic fevers, it exhibits similar effectiveness for Lassa fever, Crimean-Congo hemorrhagic fever, and hantavirus infection. Remdesivir (GS-5734) is a nucleotide analog prodrug with broad-spectrum antiviral activity, which has been shown to inhibit lovirus, coronavirus, and paramyxovirus replication. 37, 38 In vitro, remdesivir shows potent antiviral activity against both the Malaysian and Bangladesh genotypes of Nipah virus. It reduces replication of Nipah virus Malaysia in primary human lung microvascular endothelial cells by more than four orders of magnitude. 38 This indicates that the further testing of the efficacy of remdesivir against Nipah virus infection in vivo is appropriate.

Angiotensin-Converting Enzyme 2 (ACE2) is an integral membrane glycoprotein, which is known to be expressed in most tissues such as the kidney, endothelium, lungs, and heart. 39, 40 However, ACE2 is the host receptor of SARS-CoV-2 and SARS-CoV. 40, 48 Therefore, if the ACE2 protein can be inhibited, cells can be temporarily protected against infection from SARS-CoV-2. The structural database of ACE2 (https:// www.uniprot.org/uniprot/Q9BYF1) can be referenced from UniProtKB (Scheme 1A). 41 Moreover, protein PDB6LU7 is the main protease of SARS-CoV-2, functioning as a proteolytic enzyme that cuts polyproteins into functional pieces. Thus, if PDB6LU7 can be inhibited, the replication of SARS-CoV-2 can be restrained, temporizing its spread in order for human immunity to produce a sufficient number of antibodies to respond against this pathogen. The immune response can be neutralization, agglutination, or phagocytosis. The structure of PDB6LU7 (DOI: 10.2210/pdb6LU7/pdb) of SARS-CoV-2 has been archived recently by the Worldwide Protein Data Bank (Scheme 1B). 41 out molecular docking simulation on the inhibitability of several potential candidates against ACE2 and PDB6LU7. 42, 43 These studies included two drugs, ribavirin and remdesivir, which exhibited signicant inhibitory capability. Also, our preliminary investigation based on a similar virtual approach revealed that these two drugs interacted with the above two proteins mainly via their N-heterocyclic groups. This led to the implication of promising inhibitory efficacy of tetrylene complexes considering their resemblance of their molecular structure to N-heterocyclic rings.

Our research included theoretical calculations of the quantum properties of the complexes using density functional theory (DFT) and further predictions on their inhibitability developed by molecular docking simulation (MDS). Firstly, the nature of bonding and the extent of s and p interactions between the tetrylone ligands NHE Ph and fragment AgCl in [NHE Ph -AgCl] (NHE-Ag) and [(NHE Ph -AgCl) 2 ] (NHE-Ag-bis) (E ¼ C, Si, and Ge) complexes were investigated (Scheme 2). The geometry and energy of the compounds NHE-Ag and NHE-Agbis and the free ligands NHE were optimized and calculated to reach the equilibrium structures of the complexes before their bond dissociation energies (BDEs) were examined with gradient-corrected density functional theory (DFT). The bonding states of the complexes were also examined by calculating the energetically low-lying occupied molecular orbitals (LOMO) for the sand p-orbitals in the complexes and free ligands using natural bond orbital (NBO) analysis. Then, the inhibitory effects of the considered NHC-Ag, NHC-Ag-bis, ribavirin and remdesivir compounds (Scheme 3) on the proteins ACE2 and PDB6LU7, which are crucial to tackle the infection of virus SARS-CoV-2 in the human physiological environment, were also investigated. To the best of our knowledge, to date, there is no experimental and theoretical information reported about the prevention of SARS-CoV-2 related to carbene complexes.

Geometry optimization of the molecules was carried out without symmetry constraints using Turbomole 6.0 (ref. 44) and Gaussian 09 (ref. 45 ) at the BP86 (ref. 46 and 47)/def2-SVP 48 level of theory. The former dened the coordinates of the molecules before the latter was utilised for optimizing the calculations. The molecules were also calculated by vibrational frequencies at the same functional to conrm that the structures were in global minimum on the potential energy surface (PES). For the silver atom, small-core quasi-relativistic effective core potentials (ECPs) were used. 49 Single-point energies at the BP86/def2-SVP level optimized geometries were calculated with the application of the frozen-core approximation for non-valence-shell electrons. The same functional for geometry optimization was carried out at the larger def2-TZVPP 50 basis set. RI approximation was used for all structure optimizations using the appropriate auxiliary basis sets. The bond dissociation energy (BDE), D e (kcal mol À1 ) was determined at the BP86 level in conjunction with the def2-TZVPP basis set using the BP86/def2-SVP level optimized geometries. Considering information on the chemical bonding in the investigated compounds, the analysis on natural bond orbital (NBO) 51 was performed to propose their intermolecular interactions and to imply donor-acceptor bonds. The NBO analysis was carried out at the BP86/def2-TZVPP//BP86/def2-SVP level of theory, which was proposed for the calculation of the Wiberg bond orders (WBI) and natural partial charges (NPA) as well as localized plotting molecular orbitals and orbital energies using NBO 5.1 available in Gaussian 09. Electron density distributions were revealed by bonding analysis. The HOMO energy, E HOMO , indicates the tendency of a molecule to donate electrons, whereas the E LUMO of a molecule refers to its electron accepting ability. The energy gap, DE ¼ E LUMO À E HOMO , illustrates the tendency of inhibition efficiency of organic molecules towards the surface of metal elements. The ionization potential (I) and electron affinity (A) of the inhibitory molecules were calculated by applying Koopmans' theorem 52 related to the HOMO and LUMO energy as: I ¼ ÀE HOMO and A ¼ ÀE LUMO . The obtained ionization potential and electron affinity values were used to yield the electronegativity (c), global hardness (h), and global soness (S) of the molecules, which were calculated using the following three equations: c ¼ (I + A)/2; h ¼ (I À A)/2; and S ¼ 1/h, respectively.

All calculations for molecular docking simulation were investigated using the MOE 2015.10 soware. The structural information of the proteins (ACE2 and PDB6LU7) and the structures of 6 potential drugs were required in order to implement the molecular docking. The structures were used to simulate the interactions between the drugs and proteins, and then to evaluate the bonding and docking score energy results. The results included the potential drug congurations, docking score (DS) energy, root-mean-square deviation (RMSD), types of interactions, and respective distances between the potential drug and proteins. A molecular docking simulation procedure involves 3 steps as follows. [53] [54] [55] [56] (a) Protein and ligand preparation. The structures of the ACE2 and PDB6LU7 proteins are available at UniProtKB 41 and the Worldwide Protein Data Bank, 57 respectively. The Quickprep tool was applied to prepare their structures and 3D protonation. The active zones of the proteins were conrmed based on the ligand position within a radius of 4.5 A, and the presence of important amino acids and of the protein structures were saved in the *.pdb format. The six potential drugs were optimized via Conj Grad for the minimum energy; termination for energy change of 0.0001 kcal mol À1 ; max interactions of 1000; and using the Gasteiger-Huckel charge. Molecular dynamics on the MOE 2015.10 system were performed, and the structures of the proteins and potential drugs were saved in the *.sdf format.

(b) Docking investigation. The docking simulation parameters were set with the number of poses at 10 and kept for further analysis of the interaction, the maximum number of solutions per iteration at 1000, and the maximum number of solutions per fragmentation at 200.

(c) Docking results analysis. The docking score energy (DS) values demonstrated the binding ability between the potential drugs and the proteins (ACE2 and PDB6LU7). Evaluation of the docking score energy (DS) provided further insight into the interactions between the potential drugs and ACE2 and PDB6LU7 proteins. The performance of an interaction on the 2D and 3D planes of the drug-protein complexes was examined. The interactions formed between the potential drugs and important amino acids in the site-site distances of ACE2 and PDB6LU7 were analysed. Hydrogen bonds, ion bonds, p-p interactions, cation-p interactions, and van der Waals interactions were detected to gain information on the hydrophilic, hydrophobic, and solvent interactions. The interactions between the potential drugs and the proteins (ACE2 and PDB6LU7) were determined, allowing an inference to be reached on the inhibitory effects of the potential drugs on the host receptor ACE2 and SARS-CoV-2 main protease (PBD6LU7).

The optimized geometry with bond length, bond angle and bending angle of NHE-Ag, NHE-Ag-bis, and ligands NHE (E ¼ C, Si, and Ge) are shown in Fig. 1 , 2, and 3, respectively. The calculated Ag-E bond lengths of NHE-Ag (E ¼ C, Si, and Ge), i.e. Ag-C, Ag-Si and Ag-Ge are 2.076 A, 2.338 A, and 2.448 A, respectively. These results are consistent with the work on lessbulky N-heterocyclic carbene, silylene, and germylene complexes of AgCl investigated by Frenking et al. 58 In detail, the calculated Ag-C bond length of NHC-Ag was the shortest (2.057 A), which increased to 2.404 A for the germylene complex. The equilibrium structure of the NHC-Ag complex in Fig. 1 shows that the lighter NHC ligands are bonded in a head-on manner to the metal fragment AgCl, i.e. a 180 alignment. However, the bending angle becomes much more noticeable if the E atom becomes heavier (bending angles, a, of NHSi-Ag ¼ 174.2 and NHGe-Ag ¼ 170.6 ). Regarding the bis-complexes, the Ag1-C and Ag2-C bond lengths in the NHC-Ag-bis carbene complexes are both 2.068 A (Fig. 2) . These are also in good agreement with the results reported in another work by Frenking et al., 59 in which both Ag1-C and Ag2-C bonds are 2.084 A in the [(NHC H ) 2 -Ag] complex. Besides, the Ag1-E and Ag2-E bonds calculated for the heavier homologues, i.e. NHSi-Ag-bis (2.363 A) and NHGe-Ag-bis (2.432 A), are longer than that in the NHC-Ag-bis adduct. The equilibrium structure of NHC-Ag-bis shows that the NHC ligands bond in a head-on orientation to the silver atoms with an aligned angle a of 180.0 . The structures of NHE-Ag-bis contain a side-on bonded ligand, whose E ¼ Si and Ge exhibit the bending angle of 174.0 in NHSi-Ag-bis and of This journal is © The Royal Society of Chemistry 2020 RSC Adv., 2020, 10, 30961-30974 | 30965 170.1 in NHGe-Ag-bis, respectively. The results are rather similar to the values previously obtained from the NHE-Ag complexes. Fig. 3 shows that the E-N bond lengths in the complexes positively correlate with the mass of their homologous free ligands. In detail, the E-N lengths in the NHE complexes are shorter than that in their NHE-bis counterparts. In addition, the optimized geometries of the silver-complexes NHE-Ag and NHE-Ag-bis (E ¼ C, Si, and Ge) at the BP86/def2-SVP level with dispersive effect are presented in Fig. S1 and S2 (ESI †), respectively. The results reveal that the recorded deviation in the structural parameters, including bond length and bonding angles, is insignicant. Therefore, the geometries obtained without dispersive interactions are reliable. Fig. 4 shows the BDEs for the Ag-NHE Me bonds in NHE-Ag and NHE-Ag-bis. The energy exhibited a signicant decrease from the carbene complex NHC-Ag (D e ¼ 53.4 kcal mol À1 ) to the germylene NHGe-Ag complex (D e ¼ 32.8 kcal mol À1 ), as shown in Fig. 4A . The calculations suggest that the bonding of the NHC ligand in NHC-Ag is the strongest and that in the heavier homologues NHE-Ag (E ¼ Si and Ge) is weaker. These observations are in good agreement with the results reported in previous studies. [60] [61] [62] The BDE tendency of AgCl-carbene and its homologous is also similar to other corresponding values for the AgCl-NHE H complexes (D e ¼ 56.5-29.9 kcal mol À1 ). 58 Moreover, the BDEs calculated for the NHE-Ag-bis systems show a different pattern. Fig. 4B indicates that the BDEs exhibited the smallest value of D e ¼ 12.6 kcal mol À1 for NHC-Ag-bis, maximum for NHSi-Ag-bis (D e ¼ 22.9 kcal mol À1 ), and then smaller for NHGe-Ag-bis (D e ¼ 18.4 kcal mol À1 ). This observation is also in good agreement with a reported work on NHE Me -(AuCl) 2 (E ¼ C-Ge) transition metal tetrylene complexes. 63 Although the BDEs consistently correlate with their associated values of either bond length or bending angle, a, obtained from the DFT analysis of the NHE-Ag complexes, the expected correlations were not observed when applied for the NHE-Ag-bis compounds. This inconsistency still requires further in-depth investigations in order to yield appropriate justication.

Bonding analysis of lighter tetrylene complexes Table 1 shows the results of the NBO partitioning analysis, including Wiberg Bond Indices (WBI) and natural partial charge (NPA) of the NHE-Ag and NHE-Ag-bis complexes. The calculated partial charges show that the metal fragments AgCl and (AgCl) 2 in the complexes always carry a negative charge from NHC-Ag (À0.23 e) and NHC-Ag-bis (À0.44 e). The magnitude of charge donation to the AgCl fragment in the NHE-Ag and NHE-Ag-bis complexes is lower than the corresponding donation observed from other transition metal moieties such as W(CO) 5 and Mo(CO) 4 in similar NHE Me ligands. 60, 62 In fact, the more negative charges of À0,47 e and À0.77 e have been reported in the transition metal fragments W(CO) 5 show the occupied molecular orbitals and orbital energies of the s-type and p-type MOs calculated derived from the two types of NHE-Ag and NHE-Ag-bis complexes at the BP86/TZVPP level. This reveals that the energy levels of the p-donor orbitals of the NHE-Ag complexes are higher than that of their s-donor orbitals, except for NHC-Ag. In contrast, the NHE-Ag-bis complexes present the opposite tendency, i.e. the energy levels of their p-donor orbitals are lower than their s-donor counterparts. Also, the extensively occupied shape of the molecular orbitals indicates that NHE Ph /AgCl and (NHE Ph ) 2 /(AgCl) 2 not only perform signicant s donation, but also exhibit noticeable p donation in the tetrylene complexes.

Moreover, Fig. 6 provides information on the interactions of the separated NHE donor-acceptor fragments and AgCl adducts in each side of the NHE-Ag-bis bis-counterparts. The former still excludes the interactions between the two NHE-AgCl bulkfragments by dative bonds in NHE-Ag-bis. Therefore, HOMO-5 and HOMO-13, illustrating the donor-acceptor bonds in the NHE-Ag-bis complexes, are plotted in Fig. 7 to investigate their intramolecular interactions. The HOMO-5 distribution of NHC-Ag-bis suggests that the NHC)Ag-Cl-Cl-Ag/NHC p backdonation types are insignicant. Also, there is weak s and p bonding between the carbon central atoms and silver atoms in the metal fragments (AgCl). The profound interaction relates to (AgCl)-(AgCl). Furthermore, the distributions of HOMO-13 for NHSi-Ag-bis and NHGe-Ag-bis suggest that the NHSi)Ag-Cl-Cl-Ag/NHSi and NHGe)Ag-Cl-Cl-Ag/NHGe s, p, and p This journal is © The Royal Society of Chemistry 2020 RSC Adv., 2020, 10, 30961-30974 | 30967 back-donation types are substantial. This observation can be explained by the end-on bonded mode characterized for carbene complexes, leading to the major stabilizing strength of the fragments (Ag-Cl-Cl-Ag), and thus is unconducive for mirrorstabilization, i.e. p-back-donation contributions.

To explain the end-on and side-on structures of the carbenes, silylenes, and germylenes, the highest occupying molecular orbitals according to the sand p-orbital of the free NHE (E ¼ C, Si, and Ge) ligands were plotted, as shown in Fig. 8. Fig. 8 displays the shape of HOMO, HOMO-2, and HOMO-3 for the free ligands considering their symmetrical states. Most of the HOMOs of the NHE ligands emerge with p symmetry except for NHC, while all the HOMOs representing s bonds are symmetrical, including HOMO (NHC), HOMO-2 (NHSi), and HOMO-3 (NHGe). Fig. 8 also shows that the energy level of the p orbitals increases with an increase in the mass of the ligand, while the corresponding gures for the s orbitals show the reverse tendency. This can be explained by the fact that s-orbital energy decreases due to the shrinking phenomenon of the sorbital resulting from the relativistic effect for the heavier homologues. 60,61 These energy level tendencies of the energetically highest-lying sand p-orbitals of the NHE ligand verify the preference of the heavier NHSi and NHGe ligands for side-on coordination to the metal, in which the s-donation takes place through the p orbital of the ligand. 60 Table 2 summarizes the quantum chemical parameters related to the molecular electronic structures of the studied NHE-Ag, NHE-Ag-bis (E ¼ C, Si, and Ge) complexes. In principle, the E HOMO value represents the electron-accepting capability and E LUMO provides information on the capability of electron donation. Theoretically, either a high HOMO energy or low LUMO energy will be conducive to the binding ability of the complexes to protein since polypeptide molecules have been proposed and well-proven to exhibit electrical conductivity. 64, 65 The modern explanation for the electron transfer includes the super-exchange theory (or electron tunneling) and the electron hopping model. 66 The E HOMO values of the silver-tetrylene complexes follow the decreasing order of NHC-Ag > NHSi-Ag > NHGe-Ag and the E HOMO values of the bis-silver-tetrylene complexes follow a similar trend, i.e. NHC-Ag-bis > NHSi-Ag-bis > NHGe-Ag-bis. These observations are also consistent with the tendencies observed from E LUMO , but with a narrower differential. In detail, the HOMO energy varies from 4.879 eV to 5.617 eV and the LUMO energy is in the range of 2.188-3.178 eV. Overall, the energy gap (DE GAP ) values of the NHE-Ag-bis complexes are slightly lower than that of their mono counterpart silylenes (2.767 eV to 2.798 eV) and germylene (2.359 eV to 2.439 eV). Conversely, the carbene complexes show the opposite pattern, in which the corresponding value for NHC-Ag (2.653 eV) is smaller than that of NHC-Ag-bis (2.710 eV). In particular, NHGe-Ag-bis exhibits the smallest energy gap (2.359 eV), while NHSi-Ag possesses the largest energy gap (2.798 eV). Table 2 also shows that the electronegativity (c) value exhibited the decreasing order of NHC-Ag > NHSi-Ag > NHGe-Ag and NHC-Ag-bis > NHSi-Ag-bis > NHGe-Ag-bis and the largest value belongs to NHGe-Ag (4.398). Thus, the noticeable quantum chemical parameters, including energy gap (DE GAP ), ionization potential (I), electron affinity (A), and electronegativity (c), imply that all the studied complexes can be considered as versatile ligands and can to create polar interactions with highly polarized amino acids in the protein structure.

In summary, although the DFT results suggest that the tetrylene complexes containing Si and Ge exhibit desirable quantum properties for practical applications, their potential applicability in physiological systems is limited due to the possible biological toxicity of their elements. Therefore, NHC- Table 2 Quantum chemical parameters of the tetrylene NHE-Ag and NHE-Ag-bis (E ¼ C, Si, and Ge) complexes obtained from the NBO data calculated at the BP86/def2-TZVPP level, including the energy gap (DE GAP ), ionization potential (I), electron affinity (A) and electronegativity (c). Energy in eV Ag and NHC-Ag-bis were selected for further study using molecular docking simulation in an attempt to investigate the deactivation of the SARS-CoV-2 PDB6LU7 enzymes and protection of the host ACE2 receptors. Also, the inhibitory activities of the studied silver complexes will be compared with that of ribavirin and remdesivir, which were considered as references.

Molecular docking simulation was used to investigate the interactions of the silver-carbene NHC-Ag and bis-silver-carbene NHC-Ag-bis with two proteins, including ACE2 in the human body and PDB6LU7 of SARS-CoV-2. In detail, we determined whether or not the silver-carbene NHC-Ag and bissilver-carbene NHC-Ag-bis interact with the ACE2 and PDBLU7 proteins and inhibit their actions, and weaken them compared to their original state. The docking score energy (DS) and rootmean-square deviation (RMSD) values as well as various interaction congurations via hydrogen bonds, cation-p, p-p bonds, and ionic interactions, site-site binding, and van der Waals interactions between the carbene complexes, two control This journal is © The Royal Society of Chemistry 2020 RSC Adv., 2020, 10, 30961-30974 | 30969 drugs and ACE2 and PDB6LU7 proteins are presented in Fig. 9 , 10, Tables 3 and 4 . Firstly, the docking simulations on the ACE2 protein were investigated, and the obtained results are presented in Fig. 9 and Tables 3 and 4 . It can be seen that all the root-mean-square deviation values are smaller than 2 A. The NHC-Ag-bis complex docked with the ACE2 protein exhibits a stronger inhibitory effect in comparison to the that derived from NHC-Ag with the docking score energy values of À17.5 kcal mol À1 and À14.1 kcal mol À1 , respectively. This can be explained by the fact that NHC-Ag-bis possesses a larger volume and a higher molecular mass, leading to a signicant polarizability and strong binding capacity with amino acids. Furthermore, Tables  3 and 4 show that the NHC-Ag complex forms 10 van der Waals interactions with the protein and 3 hydrogen bonds with its p-H bonds between the phenyl ring and N-heterocyclic of the carbene complex and -Cof different amino acids, including Glu 208 (3. Table 3 Docking simulation results with docking score energy (DS), root-mean-square deviation (RMSD), and van der Waals (VDW) interactions of NHC-Ag, NHC-Ag-bis, ribavirin, and remdesivir with the amino acids of the ACE2 and PDB6LU7 proteins Table 5 summarizes the docking parameters including DS average (kcal mol À1 ), polarizability ( A 3 ) and total hydrogen Table 4 Molecular docking simulation results with critical interactions between the complexes and the two proteins (ACE2 and PDB6LU7), including interaction and distance, site-site binding, energy, cation-p, p-p bonds, ionic interactions, and total hydrogen bonds bonds when each complex/drug is docked with the ACE2 and PDB6LU7 proteins. These parameters were calculated using the QSARIS system with the Gasteiger-Marsili method. 67 The ACE2 and PDB6LU7 inhibition ability of NHC-Ag, NHC-Ag-bis, ribavirin, and remdesivir is expressed by the average DS values of À14.1, À17.2, À16.6, and À16.9 kcal mol À1 , respectively. These indicate that all the compounds can be considered to exhibit strong inhibitory effects on ACE2 and PDB6LU7. The justication is based on their relative signicance in comparison to the DS values representing the inhibitory capability of other duosystems previously reported in the literature. 42, 43, 68, 69 In addition, the polarizability values of NHC-Ag-bis and remdesivir are signicantly higher (i.e. 71.3 and 57.1 A 3 ) than the corresponding values for NHC-Ag and ribavirin (35.7 and 21.1 A 3 ).

The reason for this may be related to the fact that they possess the largest molecular mass and volume, which lead to high binding capacity with amino acids. Although NHC-Ag and ribavirin possess a smaller molecular mass and polarizability values, they also exhibit relatively strong site-site binding due to their hydrogen bond interactions, i.e. 6 and 10 hydrogen bonds, respectively. In summary, the results imply a feasible correlation between the structures of the silver-carbene and bis-silvercarbene complexes, and anti-SARS-CoV-2 activities considering their docking capability towards both the virus main protease PDB6LU7 and the host complementary receptor ACE2.

The calculated equilibrium structures of the NHC-Ag and NHC-Ag-bis complexes reveal that the NHC tetrylene ligands bond in a head-on mode to the AgCl and (AgCl) 2 bis-silver chloride fragments, while the NHE (E ¼ Si and Ge) ligands bond in a side-on conguration to AgCl and (AgCl) 2 . The calculated BDEs suggest that the NHE-AgCl bond strength decreases in the order of NHC-Ag > NHSi-Ag > NHGe-Ag. Meanwhile, the corresponding calculation implies that the bis-silver chloride-NHE bond strength decreases in the order of NHC-Si-bis > NHGe-Agbis > NHC-Ag-bis. The NBO analysis indicates that [NHE Ph / AgCl] and [(NHE Ph ) 2 /(AgCl) 2 ] donations of the complexes are mainly derived from the sand p-contributions of the ligands. The docking simulation results suggest that the studied carbene complexes (NHC-Ag and NHC-Ag-bis) are active for inhibiting the host receptor ACE2 and PDB6LU7 protein of SARS-CoV-2. However, NHC-Ag-bis has stronger inhibitory effects on both the ACE2 and PDB6LU7 proteins. The NHC-Agbis bis-silver-carbene complex, ribavirin, and remdesivir exhibit similar efficacy in inhibiting either ACE2 or PDB6LU7. Their corresponding docking score energies are À17.5 to À16.5 kcal mol À1 and À16.9 to À16.6 kcal mol À1 , respectively. The root-mean-square deviation values are always less than 2 A in all the calculated systems. The structures of the potential drugs t well with the site-site binding of the host receptor ACE2 and the viral protease PDB6LU7 based on the hydrogen bond interactions. Interestingly, the results also indicate that the inhibitory ability of potential drugs on the ACE2 and PDB6LU7 proteins seems to be correlated among the average of docking score energy, site-site active interactions of the potential complexes/drugs-proteins, and polarizability of the potential carbene and bis-carbene complexes and drugs. Thus, the obtained results in this study suggest that the NHC-Ag carbene complex and NHC-Ag-bis bis-carbene adduct are promising complexes, which can serve as a reference source of data for further research on the development of new agents to inhibit the host receptor ACE2 and the main protease PDB6LU7 of SARS-CoV-2.

The results in this study are highly conducive for the preparation of therapeutic drugs for SARS-CoV-2 considering the theoretical demonstration of the stability of the potential drugs based on the DFT results and their inhibitory effectiveness by molecular docking simulation.

There are no conicts to declare.

Analogue-based Drug Discovery III

World Health Organization model list of essential medicines: 21st list 2019, World Health Organization

The crystal structure of COVID-19 main protease in complex with an inhibitor N3