Next Article in Journal
Inhibitory Effect of Osthole from Cnidium monnieri (L.) Cusson on Fusarium oxysporum, a Common Fungal Pathogen of Potato
Previous Article in Journal
Dynamics of Phloridzin and Related Compounds in Four Cultivars of Apple Trees during the Vegetation Period
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Search for Non-Protein Protease Inhibitors Constituted with an Indole and Acetylene Core

by
Marco A. Almaraz-Girón
1,
Ernesto Calderón-Jaimes
2,*,
Adrián Sánchez Carrillo
2,
Erik Díaz-Cervantes
3,*,
Edith Castañón Alonso
2,
Alejandro Islas-Jácome
1,
Armando Domínguez-Ortiz
1 and
Sandra L. Castañón-Alonso
1,*
1
Departament de Química, Universidad Autónoma Metropolitana-Iztapalapa, San Rafael Atlixco 186, Col. Vicentina, Iztapalapa, Ciudad de México C.P. 09340, Mexico
2
Laboratory de Investigación en Inmunoquímica, Unidad de Investigación en Inmunología Proteómica, Hospital Infantil de México Federico Gómez, Calle Dr. Márquez Nº 162, Col. Doctores, Delegación Cuauhtémoc, Ciudad de México C.P. 06720, Mexico
3
Centro Interdisciplinario del Noreste, Departament de Alimentos, Universidad de Guanajuato, Tierra Blanca, Guanajuato C.P. 37975, Mexico
*
Authors to whom correspondence should be addressed.
Molecules 2021, 26(13), 3817; https://doi.org/10.3390/molecules26133817
Submission received: 19 March 2021 / Revised: 22 May 2021 / Accepted: 25 May 2021 / Published: 23 June 2021
(This article belongs to the Section Organic Chemistry)

Abstract

:
A possible inhibitor of proteases, which contains an indole core and an aromatic polar acetylene, was designed and synthesized. This indole derivative has a molecular architecture kindred to biologically relevant species and was obtained through five synthetic steps with an overall yield of 37% from the 2,2′-(phenylazanediyl)di(ethan-1-ol). The indole derivative was evaluated through docking assays using the main protease (SARS-CoV-2-Mpro) as a molecular target, which plays a key role in the replication process of this virus. Additionally, the indole derivative was evaluated as an inhibitor of the enzyme kallikrein 5 (KLK5), which is a serine protease that can be considered as an anticancer drug target.

1. Introduction

Proteases are essential for replication/transmission of viruses, parasites, and bacteria that cause infectious disease in mammals. These enzymes regulate most physiological processes by controlling the activation, synthesis, and turnover of all proteins. Inhibitors of such proteases could potentially be useful in the treatment of diseases as diverse as cancer; parasitic, fungal, and viral infections (e.g., schistosomiasis, malaria, C. Albicans, human immunodeficiency virus (HIV), hepatitis, herpes); and inflammatory, immunological, respiratory, cardiovascular, and neurodegenerative disorders including Alzheimer′s disease [1]. On the other hand, proteases are classified according to their catalytic site [2]. Inhibitors can be generally classified into two large groups based on their structural dichotomy: low molecular weight peptidomimetic inhibitors and protein protease inhibitors composed of one or more peptide chains. Protease inhibitors can be further classified into five groups (serine, threonine, cysteine, aspartyl, and metalloprotease inhibitors) according to the mechanism employed at the active site of proteases they inhibit [3].
Otherwise, chymotrypsin-like serine proteases have been found in vertebrates, bacteria, and also viruses; the hepatitis C virus NS3 protease is an example of a viral protease in this family [4]. Therefore, the NS3 protease of hepatitis C virus (HCV) has been studied intensively over the past few years as a target for antiviral drugs [3,4,5]. Several research groups have explored the use of indole derivatives to be candidates for HCV protease inhibitors. Ismail [6,7] reported a series of indole derivatives, represented by 1 (Figure 1), as potent HCV NS3/4A serine protease inhibitors; molecular modeling simulation studies were used in order to predict the biological activity of the proposed compounds. Narjes et al. [8,9,10] reported a series of allosteric indole-N-acetamide inhibitors, typified by 2 (Figure 1), which are potent inhibitors of the HCV NS5B polymerase. Inhibition of NS5B enzyme has emerged as an attractive strategy for the development of antivirals to combat HCV because NS5B plays a crucial role in viral replication.
It is interesting to note that, lopinavir/ritonavir is a co-formulated human immunodeficiency virus (HIV)-specific protease inhibitor [11], interferon-α2b is a registered agent for chronic hepatitis C [11], umifenovir (3, Figure 1) has been demonstrated to inhibit in vitro infection with globally prevalent pathogenic viruses, including hepatitis C, and umifenovir has an indole core [12]. Viral proteases are well-validated drug targets that have led to various approved drugs, for example, against chronic infections with HIV or hepatitis C virus (HCV) [13], which employ aspartyl and serine proteases, respectively [1,14]. Further, several reviews have been published describing compounds containing acetylene groups and related functional groups, isolated from plants, fungi, microorganisms, and marine invertebrates, which present interesting biological activities [15,16,17]. There are acids with two groups of consecutive alkynes that act as inhibitors of HIV protease [15].
On the other hand, one of the best-characterized drug targets among coronaviruses is the main protease (Mpro, also called 3-chymotrypsin-like cysteine protease, 3CLpro); inhibiting the activity of this enzyme would block the viral replication. Because no human proteases with a similar cleavage specificity are known, such inhibitors are unlikely to be toxic [18,19]. Dai et al. [20] designed indole derivatives 4 (Figure 1), which are the strongest known SARS-CoV-2 Mpro inhibitors, with inhibition constants in the two-digit nanomolar range [14,20]. Different drugs were tested as possible therapeutic candidates against SARS-CoV-2 (severe acute respiratory syndrome–coronavirus 2) infection [11,12], such as: remdesivir, lopinavir/ritonavir, favipiravir, umifenovir, EIDD-2801, convalescent plasma collected from donors who have survived an infectious disease, interferons, corticosteroids, vitamin C, vitamin D, zinc, the non-steroidal anti-inflammatory drugs (NSAIDs) [21], antibodies (for example, REGN-COV2) [22], among others. On the other hand, studies in silico to evaluate the potential of hundreds [23] or thousands [24,25] or millions [13] of molecules against the active site of SARS-CoV-2 have been reported [26,27,28].
Moreover, kallikrein-related peptidases (known as KLKs) constitute a single family of 15 highly conserved trypsin- or chymotrypsin-like serine proteases encoded by the largest uninterrupted cluster of protease-encoding genes (KLK1–15) in the human genome [29]. Aberrant KLK expression patterns have been reported mainly in hormone-dependent malignancies, such as those in breast, ovary, and prostate, and have been widely implicated as cancer biomarkers. In ovarian cancer, overexpression of KLK5 protease in tumor tissues is associated with more advanced stages and grade of the disease [30]. On the other hand, KLK5 plays a central role in the degradation of corneodesmosomes, which are the main adhesive structures in the cornified cell layer. KLK5 is also involved in the activation of other epithelial serine proteases such as kallikrein-related peptidase. Although KLK5 is important for the maintenance of skin homeostasis, its overactivation can impair the skin barrier function and contribute to the discovery of various dermatoses [31]. Kim et al. [31] reported that the indole-3-acetonitrile-4-methoxy-2-S-β-d-glucopyranoside KLK5 protease activity. Additionally, the KLK5 protein has been employed as a molecular target for drug-transport single-wall carbon nanotubes and was reported their in silico assay by Díaz-Cervantes et al. [32].
Furthermore, lipophilicity constitutes a physicochemical property of paramount importance in medicinal chemistry [33]. Lipophilicity is usually measured by the partition coefficient, Log P, or log Kow (logarithm of n-octanol/water partition coefficient), which plays an essential role in absorption, distribution, metabolism, and excretion (ADME) characteristics of drugs, while also affecting their pharmacodynamic and toxicological profiles, such as the extent of plasma protein binding, accumulation in tissues, and unpredictable poisonous mechanism [33,34]. Regarding the activity and log P relationship, the low log P values mean less activity of the molecule in some cases [35].
Herein, we disclose our findings on the synthesis of derivative 15 (Scheme 1), which contains an indole core and aromatic polar acetylene, as a possible potential protease inhibitor. Derivative 15 was evaluated through docking assays, in which the main protease (SARS-CoV-2-Mpro) was used as a molecular target. This protein plays a key role in the replication process of this virus and has been reported in crystallographic data [19,36]. Despite there being other reports of recombinant virus proteases, the current work selected the non-recombinant protein [37]. Additionally, we decided to evaluate the indole-acetylene derivative 15 as an inhibitor of the KLK5 protein in order to see the scope of 15 as an inhibitor of different proteases.

2. Results and Discussion

2.1. Chemical Synthesis

Our synthesis strategy to access 2,2′-[(4-{[4-({4-[(1E)-3- (1H-indol-3-yl)-3-oxoprop-1-en-1-yl]phenoxy}acetamido)phenyl]ethynyl}phenyl)azanediyl]di(ethane-2,1-diyl)dipropionate (15) began with the iodination of 2,2′-(phenylazanediyl)di(ethan-1-ol) (5), which was treated with a mixture of KI, MeOH, and water, after which it was treated with Na2SO3 to give the iodophenyl-aminoethanol 6 with a yield of 81%. Then, a Sonogashira reaction between 6 and the alkyne 7 was performed to provide 8 in 87% yield [38]. An esterification was carried out between 8 and the propionic acid (11) using a propylphosphonic anhydride solution (T3P®) [39,40,41]; the ester 10 was obtained in 79% yield. The reduction of the nitro group of 10 was carried out using sodium hydrosulfite (Na2S2O4) [42,43]; the amine 11 was obtained in 82% yield. Indole derivative 14 was formed from indole 12 and acid 13 as described by Ismail et al. [6]. Finally, the condensation reaction between 11 and 14 using a propylphosphonic anhydride solution (T3P®) [39,40,41] allowed achievement of the formation of the amide bond, with an 83% yield. Compound 15 was obtained through the five synthetic steps with a total yield of 36.98% (Scheme 1).
Once the chemical structure of the indole-acetylene derivative 15 was confirmed, docking studies were conducted to find out the potential of molecule 15 as a coronavirus protease inhibitor COVID-19 and inhibitor of serine proteases.

2.2. Docking Studies with SARS-CoV-2-Mpro

Figure 2 shows the in silico molecular couplings between the SARS-CoV-2-Mpro and the four selected ligands: indole-acetylene derivative 15, co-crystallized molecule 16 (known as N3, which is a Michael acceptor inhibitor) [36], 1,5-disubstituted tetrazole-1,2,3-triazoles 17 and 18 (Figure 3) [26]. Note that, all the ligands were coupled in the right-site cavity, which contains the active site as explained below, and 1,5-disubstituted tetrazole-1,2,3-triazoles 17 and 18 were selected due to their higher selectivity to the SARS-CoV-2-Mpro, reported by Díaz-Cervantes and co-workers [28].
The derivative for the molecular coupling is depicted in Figure 2; the main interactions are shown in Table 1. The ligands are ordered from the smallest to largest interaction with the SARS-CoV-2-Mpro, with the 18 molecule being the best ligand for this target. This can be explained with the LE values, and it can be noticed that the ligand with better interaction presents higher Hbond and VdW interactions than the others.
Although 15 is not the best ligand in Table 1, this presents a more favorable value of LE than the N3 inhibitor 16 (so-called Co-Crystal in the present work) and exergonic energy, which means that this ligand can interact in a favorable way with the SARS-CoV-2-Mpro. It is important to highlight in Table 1 that Co-Crystal 16, 17, 18, 1a, 1b, 2, 3, and 4 ligands were the only molecules selected to present a comparison point of the interactions between the SARS-CoV-2-Mpro and the synthesized molecule 15.
Considering that the aim of this work was to evaluate the interactions between 15 and the protease present in the pandemic virus (SARS-CoV-2-Mpro), Figure 4 shows the position of the ligand in the active site of the selected protease. Note that the virus protein is a cysteine protease with the Cys145 (and Ser144) and that His41 is the key residue in the catalytic site, which interacts with the central site of the synthesized ligand.
Regarding the obtained pose for 15 (Figure 4), the key interactions are depicted in Figure 5. The principal residues of the SARS-CoV-2-Mpro that interact with 15 are the Ser1, His164, and Phe305 in the case of the hydrogen bond interactions and His41 and Glu166 for the electrostatic interactions.
Moreover, to understand the bioactive conformation obtained from molecular docking it is necessary to evaluate the ligand–target hydrophobic interactions. Figure 6 shows the position of 15 due to the hydrophobic behavior of the protein. Hydrophobic and hydrophilic surfaces are present in the active site and the molecule gets a pose dependent on these interactions, in addition to the HBond, VdW, and Elstat interactions. In this order, is clear that the central site of 15 is surrounded by hydrophobic surfaces of the viral target, at the same time that some hydrophilic moieties of this interact with the red (hydrophilic) surfaces of the target, promoting the bioactive conformation of this ligand. To better understand the ligand–target interactions and to evaluate the possible change of conformation of the target, a normal mode analysis of the target was conducted using the iMod server [44] which demonstrated that the active sites present lower modification.
Finally, to conclude the analysis of the interactions between the SARS-CoV-2-Mpro and 15, the bioactive conformation of the latter was inserted into the SARS-CoV-2-Mpro pharmacophore model computed by Díaz-Cervantes and coworkers [26] to generate an evaluation of the fragments that need the molecule to obtain better interactions with the target. Figure 7 shows that 15 needs one Hy fragment (on the left side) and two HA fragments and one HD fragment in the upper site to obtain a better interaction with the SARS-CoV-2-Mpro.

2.3. Docking Studies with KLK5 Protease

The interaction energies of 15 with another therapeutic target, the KLK5 protease, are shown in Table 2. In this case, a nanostructured drug-transport (a single-wall carbon nanotube functionalized [32] with altretamine, a so-called Tube-altre) has been used as a comparative point (Figure 8A).
Table 2 shows that the ligand efficiency of 15 is better than that obtained by the Tube-altre ligand, apparently due to the higher hydrogen bond interactions of 15 than with another molecule. Note that KLK5 is a serine protease, His57 and the Ser195 being two of the key residues in the active site, which interact with the central site of 15 (Figure 8B).
Additionally, 15 interacts with KLK5 protease in terms of hydrogen bonds through Lys60, Gly193, Asp194, Ser195, and Ser214 (see Figure 9A). At the same time, Figure 9B shows the electrostatic interaction between 15 and His57, Lys60, Asp217, and Arg224. Note that the hydrogen bond is the key interaction in the ligand–target coupling.
Moreover, the hydrophobic interactions were computed, and it is clear that 15 obtains a bioactive pose due to the hydrophobic surfaces present in the active site of the KLK5 protease, as Figure 10 depicts, which is in accordance with the hydrophobic surface study of MPro. The hydrophilic surfaces of 15 are on the tips of the molecule and the hydrophobic moieties (the center of 15) and interact with the blue surfaces of the target.
Finally, to obtain the first approach with regard to the possible biodistribution of our synthesized molecule, we conducted a logP analysis of it through the server ALOGPS 2.1 [45]. The computed logP value was compared with the reference molecules. Table 3 shows that the logP for 15 is higher than that of other molecules, which can mean a better distribution of this molecule through the cellular membranes and perhaps into the virus capsid.

3. Materials and Methods

3.1. Materials and Instrumentation

Commercial reagents were purchased from Sigma Aldrich (St. Louis, MO, USA) and were used as received. Proton nuclear magnetic resonance (1H NMR) spectra and carbon nuclear magnetic resonance (13C NMR) spectra were recorded on a Bruker Advance 400 MHz spectrometer and a Varian Unity Inova 300 MHz and 75 MHz; both are reported in parts per million and are referenced to tetramethylsilane or residual solvent. Data are reported as follows: Chemical shift, integration, multiplicity (br = broad, s = singlet, d = doublet, t = triplet, q = quartet, m = multiplet), coupling constant in Hertz (Hz). Infrared (IR) spectra were obtained using a (in total attenuated reflection: FTIR-TAR mode) Perkin Elmer Spectrum 400 and Perkin Elmer 1605 spectrometer by KBr. Data are represented in frequency of adsorption (cm−1). Mass spectrometry was obtained with a JEOL SMX-102 A instrument (FAB+).

3.2. Chemical Synthesis

3.2.1. 2,2′-[(4-iodophenyl)azanediyl]di(ethan-1-ol) (6)

2,2′-[(4-iodophenyl)azanediyl]di(ethan-1-ol) (6) was prepared according to the procedure described in the reference [39].

3.2.2. 2,2′-({4-[(4-nitrophenyl)ethynyl]phenyl}azanediyl)di(ethan-1-ol) (8)

2,2′-[(4-iodophenyl)azanediyl]di(ethan-1-ol) (6) (2 g, 6.52 mmol) was dissolved in TEA (20 mL) and tetrahydrofuran (THF 20 mL); then 1-ethynyl-4-nitrobenzene (7) (1 g, 6.79 mmol) was added to this solution and stirred. This was followed by the addition of triphenylphosphine (PPh3, 0.02 g, 7.61 × 10−4 mol), copper iodide (CuI, 0.2 g, 1.05 mmol), and chloride bis triphenylphosphine palladium (II) ([Pd(PPh3)2Cl2], 0.2 g, 2.85 × 10−4 mol). The mixture was heated at 50 °C under nitrogen for 48 h, changing its color from clear yellow to orange. Afterwards, the amine salt was separated from the reaction mixture by filtration, and the filtrate was concentrated by evaporation. The obtained solid was then washed with THF, hexane, and finally acetone in order to remove the amine salt and catalysts. The product was purified by column chromatography on silica gel (1:1 hexane/AcOEt, added 5% of MeOH) to afford the desired product in 87% yield (from 6) as an orange powder. 1H NMR (400 MHz, Acetone d-6) δ (ppm): 8.24 (d, 2H, J = 8.8 Hz), 7.70 (d, 2H, J = 8.8 Hz), 7.39 (d, 2H, J = 9.2 Hz), 6.80 (d, 2H, J = 9.2 Hz), 3.34 (S, 2H) 3.78 (t, 4H, J = 5.4 Hz), 3.63 (t, 4H, J = 5.8 Hz), 13C NMR (100 MHz, Acetone d-6) δ (ppm): 149.2, 146.4, 133.14, 131.6, 131.2, 123.7, 111.7, 107.7, 96.9, 86.05, 59.1, 53.9. FT-IR (KBr), ν (cm−1): 3229, 3098, 2958, 2206, 1584, 1511, 1399, 1349, 1044, 856, 799. (FAB+) (m/z): 326. m.p. = 151–153 °C.

3.2.3. 2,2′-({4-[(4-nitrophenyl)ethynyl]phenyl}azanediyl)di(ethane-2,1-diyl)dipropianate (10)

2,2′-({4-[(4-nitrophenyl)ethynyl]phenyl}azanediyl)di(ethan-1-ol) (8) (1 g, 3.67 mmol, 1 equiv) was dissolved in THF (5 mL); then propionic acid (9) (1 mL, 13.5 mmol) was added to this solution and was stirred. The resulting solution was treated with pyridine (0.5 mL, 1 equiv) and this was followed by the addition of propylphosphonic anhydride solution (T3P®) (2 mL, 3 equiv). The mixture was heated at 50 °C under nitrogen for 48 h; the solvent was concentrated under reduced pressure and after was neutralized with sodium carbonate solution. The product was purified by column chromatography on silica gel (7:3 hexane/AcOEt) to afford the desired product in 79% yield (from 8) as orange oil. 1H NMR (300 MHz, CDCl3) δ (ppm): 8.17, (d, 2H, J = 9 Hz), 7.60 (d, 2H, J = 9 Hz), 7.42 (d, 2H, J = 9 Hz), 6.74 (d, 2H, J = 9 Hz), 4.27 (t, 4H, J = 6 Hz), 3.67 (t, 4H, J = 6 Hz), 2.33 (q, 4H, J = 9 Hz), 1.13 (t, 3H, J = 6 Hz). 13C NMR (75 MHz, CDCl3) δ (ppm): 174.4, 148, 146.3, 133.5, 131.7, 123.6, 111.7, 109.5, 96.5, 86.6, 61.1, 49.5, 27.5, 9.1. FT-IR (KBr), ν (cm−1): 3297, 2927, 2857, 2209, 1719, 1607, 1518, 1380, 1339, 1231, 861, 809.

3.2.4. 2,2′-({4-[(4-aminophenyl)ethynyl]phenyl}azanediyl) di(ethane-2,1-diyl)dipropianate (11)

2,2′-({4-[(4-nitrophenyl)ethynyl]phenyl}azanediyl)di(ethane-2,1-diyl)dipropianate (10) (0.5 g, 1.629 mmol) was dissolved in MeOH (30 mL), then an excess of sodium hydrosulfite salt (NaHSO3, 7 g, 40 mmol) in H2O (40 mL) was added to this solution and stirred. The resulting solution was stirred at room temperature under nitrogen for 24 h. After, the solvent was concentrated under reduced pressure, and the residue was taken up in water (40 mL) and was extracted with AcOEt (3 × 150 mL). The ethyl acetate was dried over anhydrous MgSO4 and was concentrated under reduced pressure. The product was purified by column chromatography on silica gel (1:1 hexane/AcOEt) to afford the desired product in 82% yield (from 10) as yellow oil. 1H NMR (300 MHz, CDCl3) δ (ppm): 7.29, (d, 2H, J = 9 Hz), 7.21 (d, 2H, J = 9 Hz), 6.63 (d, 2H, J = 9 Hz), 6.56 (d, 2H, J = 9 Hz), 4.18 (t, 4H, J = 6 Hz), 3.56 (t, 4H, J = 9 Hz), 2.25 (q, 4H, J = 9 Hz), 2.02 (s, 2H), 1.05 (t, 6H, J = 9 Hz). FT-IR (KBr), ν (cm−1): 3317, 3061, 2958, 2928, 2210, 1722, 1682, 1583, 1465, 1355, 1266, 1165, 735.

3.2.5. {4-[(1E)-3-(1H-indol-3-yl)-3-oxoprop-1-en-1-yl]phenoxy} acetic acid (14)

{4-[(1E)-3-(1H-indol-3-yl)-3-oxoprop-1-en-1-yl]phenoxy} acetic acid (14) was prepared according to the procedure described in reference [6].

3.2.6. 2,2′-[(4-{[4-({4-[(1. E)-3-(1H-indol-3-yl)-3-oxoprop-1-en-1-yl]phenoxy}acetamido)phe-nyl]ethynyl}phenyl)azanediyl]di(ethane-2,1-diyl)dipropionate (15)

A mixture of compound {4-[(1E)-3-(1H-indol-3-yl)-3-oxoprop-1-en-1-yl]phenoxy}- acetic acid (14) (0.5 g, 1.557 mmol) and 2,2′-({4-[(4-aminophenyl)ethynyl]phenyl}azanediyl) bis(ethane-2,1-diyl)dipropianate (11) (0.5 g, 1.225 mmol) was dissolved in THF (10 mL) and was stirred. Then, the solution was treated with pyridine (1.5 mL, 1 equiv) and was followed by the addition of propylphosphonic anhydride solution (T3P®) (2.5 mL, 3 equiv). The mixture was heated at 50 °C under nitrogen for 24 h and was concentrated under reduced pressure after. The product was purified by column chromatography on silica gel (1:1 hexane/AcOEt, added 5% of MeOH, followed by elution 6:4 hexane/AcOEt,) to afford the desired product in 83% yield (from 11) as yellow oil. 1H NMR (300 MHz, CDCl3) δ (ppm): 9.87 (s, 1H), 8.83 (s, 1H), 8.33-8.30 (m, 1H), 8.18 (s, 1H), 7.87–7.82 (m, 2H), 7.79 (d, 1H, J = 3 Hz), 7.52–7.49 (m, 2H), 7.44–7.41 (m, 2H), 7.34–7.30 (m, 3H), 7.23–7.20 (m, 3H), 7.08–7.03 (m, 2H), 6.67–6.62 (m, 2H), 4.63 (s, 2H), 4.19 (t, 4H, J = 6 Hz), 3.58 (t, 4H, J = 9 Hz), 2.26 (q, 4H, J = 6 Hz), 1.05 (t, 6H, J = 6 Hz). 13C NMR (75 MHz, CDCl3) δ (ppm): 190.8, 174.6, 165.3, 148.2, 148.1, 147.4, 136.1, 132.5, 132.4, 131.8, 123.9, 123.7, 122.9, 122.6, 120.9, 120.1, 115.4, 111.9, 111.6, 109.6, 91.4, 87.3, 67.7, 61.4, 49.7, 29.9, 27.7, 9.7. FT-IR (KBr), ν (cm−1): 3350, 3041, 2977, 2940, 2202, 2124, 1736, 1680, 1602, 1520, 1439, 1383, 1038, 1155, 827, 812.

3.3. Computational Assays

To model and optimize the studied ligands, the Gaussian package (G09) [46] at the PM6 [47] level was employed. The crystallographic data of the targets, SARS-CoV-2-Mpro and KLK5, were obtained via the protein data bank resource, with the PDB codes 6LU7 and 2PSY, respectively. At the same time, the proteins were corrected, and a pH adjustment using the Chimera software was performed [48].
Once the ligands and targets were modeled, optimized, and corrected, the in silico molecular couplings were performed with the Molegro Virtual Docker package (MVD) [49] using the MolDock [50] scoring function. To complete the theoretical results, a pharmacophore model [26] for the SARS-CoV-2-Mpro was employed to overlap the synthesized ligand using the ZINCPharmer server [51].

4. Conclusions

The computational study showed that the indole and acetylene fragments in our molecule is a molecular motif that could be explored to further the design of plausible pharmacophores of varied interest for new pharmaceuticals that could be used to combat existing diseases or future pandemics that are a priority for humankind. One of the strategies used in this exploration was the search for new and better protease inhibitors. We carried out the synthesis of the hybrid compound 15 in our search for molecules that may have a synergy with an indole and acetylene core. Through in silico assays, it has been demonstrated that 15 is a molecule with a plausible inhibition activity against the SARS-CoV-2-Mpro, presenting a higher interaction than the N3 inhibitor reported in other studies. At the same time, 15 can be considered as an anticancer agent, due to the favorable interactions with the KLK5, a protein present in ovarian cancer. Finally, 15 can be proposed as a protease inhibitor, considering the present computational assays.

Author Contributions

M.A.A.-G. writing, review-reviser, conceptualization, visualization, draft writing of the manuscript and investigation; E.C.-J. funding acquisition, conceptualization and investigation; A.S.C. methodology and investigation; E.D.-C. Software, investigation and draft writing; E.C.A. methodology and investigation; A.I.-J. visualization, review and editing and investigation; A.D.-O. visualization and investigation; S.L.C.-A. investigation, project administration, conceptualization, formal analysis, methodology, supervision, visualization, review and editing. All authors contributed equally to the present work. All authors have read and agreed to the published version of the manuscript.

Funding

We would like to thank PRODEP–SEP–México and FONDOS FEDERALES–México (Grant HIM/2012 /003 SSA. 999) for financial support.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in this article.

Acknowledgments

M. A. Almaraz-Girón thanks PRODEP–SEP–México for the postdoctoral fellowship (Grant 511-6/2019-201). We are grateful to the Laboratorio Nacional de Caracterización de Propiedades Fisicoquímicas y Estructura Molecular (UG-UAA-CONACYT, Project: 123732) for the computing time provided. A.I.J. acknowledges “Proyecto Apoyado por el Fondo Sectorial de Investigación para la Educación CONACyT-SEP CB-2017-2018 (A1-S-32582)”. Language revision by Oscar González-Antonio and his fruitful comments are also acknowledged. Special thanks to our reviewers.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Abbenante, G.; Fairlie, D.P. Proteases Inhibitors in the clinic. Med. Chem. 2005, 1, 71–104. [Google Scholar] [CrossRef]
  2. Grzonka, Z.; Jankowska, E.; Kasprzykowski, F.; Kasprzykowska, R.; Lankiewicz, L.; Wiczk, W.; Wieczerzak, E.; Ciarkowski, J.; Drabik, P.; Janowski, R.; et al. Structural studies of cysteine proteases and their inhibitors. Acta Biochim. Pol. 2001, 48, 1–20. [Google Scholar] [CrossRef] [Green Version]
  3. Fear, G.; Komarnytsky, S.; Raskin, I. Protease inhibitors and their peptidomimetic derivatives as potential drugs. Pharmacol. Ther. 2007, 113, 354–368. [Google Scholar] [CrossRef] [PubMed]
  4. Tong, L. Viral Proteases. Chem. Rev. 2002, 102, 4609–4626. [Google Scholar] [CrossRef] [PubMed]
  5. Colarusso, S.; Koch, U.; Gerlach, B.; Steinkuhler, C.; Francesco, R.D.; Altamura, S.; Matassa, V.G.; Narjes, F. Phenethyl Amides as Novel Noncovalent Inhibitors of Hepatitis C Virus NS3/4A Protease: Discovery, Initial SAR, and Molecular Modeling. J. Med. Chem. 2003, 46, 345–348. [Google Scholar] [CrossRef]
  6. Ismail, N.S.M.; Hatori, M. Molecular Medeling based approach, Synthesis and in vitro assay to new indole inhibitors of hepatitis C NS3/4A Serine protease. Bioorg. Med. Chem. 2011, 19, 374–383. [Google Scholar] [CrossRef] [PubMed]
  7. Ismail, N.S.M.; Dine, R.S.E.; Hattori, M.; Takahashi, K.; Ihara, M. Computer based design, synthesis and biological evaluation of novel indole derivatives as HCV NS3-4A serine protease inhibitors. Bioorg. Med. Chem. 2008, 16, 7877–7887. [Google Scholar] [CrossRef] [PubMed]
  8. Harper, S.; Pacini, B.; Avolio, S.; Filippo, M.; Migliaccio, G.; Laufer, R.; Francesco, R.; Rowley, M.; Narjes, F. Development and Preliminary Optimization of Indole-N-Acetamide Inhibitors of Hepatitis C Virus NS5B Polymerase. J. Med. Chem. 2005, 48, 1314–1317. [Google Scholar] [CrossRef] [PubMed]
  9. Harper, S.; Avolio, S.; Pacini, B.; Di Filippo, M.; Altamura, S.; Tomei, L.; Paonessa, G.; Di Marco, S.; Carfi, A.; Giuliano, C.; et al. Potent Inhibitors of Subgenomic Hepatitis C Virus RNA Replication through Optimization of Indole-N-Acetamide Allosteric Inhibitors of the Viral NS5B Polymerase. J. Med. Chem. 2005, 48, 4547–4557. [Google Scholar] [CrossRef] [PubMed]
  10. Stansfield, I.; Pompei, M.; Conte, I.; Ercolani, C.; Migliaccio, G.; Jairaj, M.; Giuliano, C.; Rowley, M.; Narjes, F. Development of carboxylic acid replacements in indole-N-acetamide inhibitors of hepatitis C virus NS5B polymerase. Bioorg. Med. Chem. Lett. 2007, 17, 5143–5149. [Google Scholar] [CrossRef]
  11. Zhang, J.; Xie, B.; Hashimoto, K. Current status of potential therapeutic candidates for the COVID-19 crisis. Brain Behav. Immun. 2020, 87, 59–73. [Google Scholar] [CrossRef] [PubMed]
  12. McKee, D.L.; Sternberg, A.; Stange, U.; Laufer, S.; Naujokat, C. Candidate drugs against SARS-CoV-2 and COVID-19. Pharmacol. Res. 2020, 157, 104859. [Google Scholar] [CrossRef]
  13. Fischer, A.; Sellner, M.; Neranjan, S.; Smieško, M.; Lill, M.A. Potential Inhibitors for Novel Coronavirus Protease Identified by Virtual Screening of 606 Million Compounds. Int. J. Mol. Sci. 2020, 21, 3626. [Google Scholar] [CrossRef] [PubMed]
  14. Ullrich, S.; Nitsche, C. The SARS-CoV-2 main protease as drug target. Bioorganic Med. Chem. Lett. 2020, 30, 127377. [Google Scholar] [CrossRef]
  15. Zhou, Z.-F.; Menna, M.; Cai, Y.-S.; Guo, Y.-W. Polyacetylenes of Marine Origin: Chemistry and Bioactivity. Chem. Rev. 2015, 115, 1543–1596. [Google Scholar] [CrossRef]
  16. Minto, R.E.; Blacklock, B.J. Biosynthesis and function of polyacetylenes and allied natural products. Prog. Lipid Res. 2008, 47, 233–306. [Google Scholar] [CrossRef] [Green Version]
  17. Dembitsky, V.M. Anticancer Activity of Natural and Synthetic Acetylenic Lipids. Lipids 2006, 41, 883–924. [Google Scholar] [CrossRef]
  18. Zhang, L.; Lin, D.; Sun, X.; Curth, U.; Drosten, C.; Sauerhering, L.; Becker, S.; Rox, K.; Hilgenfeld, R. Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved α-ketoamide inhibitors. Science 2020, 368, 409–412. [Google Scholar] [CrossRef] [Green Version]
  19. Jin, Z.; Du, X.; Xu, Y.; Deng, Y.; Liu, M.; Zhao, Y.; Zhang, B.; Li, X.; Zhang, L.; Peng, C.; et al. Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors. Nature 2020, 582, 289–293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Dai, W.; Zhang, B.; Jiang, X.-M.; Su, H.; Li, J.; Zhao, Y.; Xie, X.; Jin, Z.; Peng, J.; Liu, F.; et al. Structure-based design of antiviral drug candidates targeting the SARS-CoV-2 main protease. Science 2020, 368, 1331–1335. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Bianconi, V.; Violi, F.; Fallarino, F.; Pignatelli, P.; Sahebkar, A.; Pirro, M. Is Acetylsalicylic Acid a Safe and Potentially Useful Choice for Adult Patients with COVID-19? Drugs 2020, 80, 1383–1396. [Google Scholar] [CrossRef]
  22. Yang, L.; Liu, W.; Yu, X.; Wu, M.; Reichert, J.M.; Ho, M. COVID-19 antibody therapeutics tracker: A global online database of antibody therapeutics for the prevention and treatment of COVID-19. Antibody Ther. 2020, 3, 205–212. [Google Scholar] [CrossRef]
  23. Shivanika, C.; Deepak, K.S.; Ragunathan, V.; Tiwari, P.; Sumitha, A.; Brindha, D.P. Molecular docking, validation, dynamics simulations, and pharmacokinetic prediction of natural compounds against the SARS-CoV-2 main-protease. J. Biomol. Struct. Dyn. 2020, 8, 1–27. [Google Scholar] [CrossRef]
  24. Joshi, R.S.; Jagdale, S.S.; Bansode, S.B.; Shankar, S.S.; Tellis, M.B.; Pandya, V.K.; Chugh, A.; Giri, A.P.; Kulkarni, M.J. Discovery of potential multi-target-directed ligands by targeting host-specific SARS-CoV-2 structurally conserved main protease. J. Biomol. Struct. Dyn. 2020, 5, 1–16. [Google Scholar] [CrossRef] [Green Version]
  25. Rahman, M.; Saha, T.; Islam, K.J.; Suman, R.H.; Biswas, S.; Rahat, E.U.; Hossen, R.; Islam, R.; Hossain, N.; Al Mamun, A.; et al. Virtual screening, molecular dynamics and structure–activity relationship studies to identify potent approved drugs for Covid-19 treatment. J. Biomol. Struct. Dyn. 2020, 21, 1–11. [Google Scholar] [CrossRef]
  26. Umesh; Kundu, D.; Selvaraj, C.; Singh, S.K.; Dubey, V.K. Identification of new anti-nCoV drug chemical compounds from Indian spices exploiting SARS-CoV-2 main protease as target. J. Biomol. Struct. Dyn. 2020, 2, 1–9. [Google Scholar] [CrossRef]
  27. Mazzini, S.; Musso, L.; Dallavalle, S.; Artali, R. Putative SARS-CoV-2 M Inhibitors from an In-House Library of Natural and Nature-Inspired Products: A Virtual Screening and Molecular Docking Study. Molecules 2020, 25, 3745. [Google Scholar] [CrossRef] [PubMed]
  28. Cortés-García, C.J.; Chacón-García, L.; Mejia-Benavides, J.E.; Díaz-Cervantes, E. Tackling the SARS-CoV-2 main protease using hybrid derivatives of 1,5-disubstituted tetrazole-1,2,3-triazoles: An in silico assay. PeerJ Phys. Chem. 2020, 2, e10. [Google Scholar] [CrossRef]
  29. Sotiropoulou, G.; Pampalakis, G.; Diamandis, E.P. Functional Roles of Human Kallikrein-related Peptidases. J. Biol. Chem. 2009, 284, 32989–32994. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Avgeris, M.; Papachristopoulou, G.; Polychronis, A.; Scorilas, A. Down-regulation of kallikrein-related peptidase 5 (KLK5) expression in breast cancer patients: A biomarker for the differential diagnosis of breast lesions. Clin. Proteom. 2011, 8, 5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Kim, Y.M.; Bae, K.S. Inhibitory Effects of Indole-3-Acetonitrile-4- Methoxy-2-S-Β-D-Glucopyranoside From Nasturtium officinale R. Br on Human Kallikrein 5 and 7 Protease. Nat. Prod. Commun. 2020, 15, 1–4. [Google Scholar] [CrossRef]
  32. Díaz-Cervantes, E.; Robles, J.; Aguilera-Granja, F. Understanding the structure, electronic properties, solubility in water, and protein interactions of three novel nano-devices against ovarian cancer: A computational study. J. Nanopart. Res. 2018, 20, 266. [Google Scholar] [CrossRef]
  33. Han, S.-Y.; Qiao, J.-Q.; Zhang, Y.-Y.; Lian, H.-Z.; Ge, X. Determination of n-octanol/water partition coefficients of weak ionizable solutes by RP-HPLC with neutral model compounds. Talanta 2012, 97, 355–361. [Google Scholar] [CrossRef] [PubMed]
  34. Liu, X.; Testa, B.; Fahr, A. Lipophilicity and Its Relationship with Passive Drug Permeation. Pharm Res. 2011, 28, 962–977. [Google Scholar] [CrossRef] [PubMed]
  35. Berkowitz, A.J.; Franson, A.D.; Cassals, A.G.; Donald, K.A.; Yu, A.J.; Garimallaprabhakaran, A.K.; Lynda, A.; Morrison, L.A.; Murelli, R.P. Importance of lipophilicity for potent anti-herpes simplex virus-1 activity of α-hydroxytropolones. Med. Chem. Commun. 2019, 10, 1173–1176. [Google Scholar] [CrossRef] [PubMed]
  36. Liu, X.; Zhang, B.; Jin, Z.; Yang, H.; Rao, Z.; Liu, X.; Zhang, B.; Jin, Z.; Yang, H.; Rao, Z. The Crystal Structure of COVID-19 Main Protease in Complex with an Inhibitor N3; Worldwide Protein Data Bank: London, UK, 2020. [Google Scholar] [CrossRef]
  37. Mody, V.; Ho, J.; Wills, S.; Mawri, A.; Lawson, L.; Ebert, M.C.C.J.C.; Fortin, G.M.; Rayalam, S.; Taval, S. Identification of 3-chymotrypsin like protease (3CLPro) inhibitors as potential anti-SARS-CoV-2 agents. Commun. Biol. 2021, 4, 93. [Google Scholar] [CrossRef]
  38. Castanon-Alonso, S.L.; Morales-Saavedra, O.G.; Baez-Pimiento, S.; Ortega-Martinez, R.; Antonio, A.; Rodriguez-Rosales, A.A.; Hernandez-Rojas, M.E. Synthesis and overall nonlinear optical characterization of poly(hexa-2,4-diynylen-1,6-dioxydibenzoate) containing 2,2-(4-((4-nitrophenyl)ethynyl)phenylazanediyl)diethanol. Mater. Chem. Phys. 2012, 133, 528–540. [Google Scholar] [CrossRef]
  39. Jablonkai, E.; Henyecz, R.; Milen, M.; Kóti, J.; Keglevich, G. T3P®-assisted esterification and amidation of phosphinic acids. Tetrahedron 2014, 70, 8280–8285. [Google Scholar] [CrossRef]
  40. Jablonkai, E.; Milen, M.; Drahos, L.; Keglevich, G. Esterification of five-membered cyclic phosphinic acids under mild conditions using propylphosphonic anhydride (T3P®). Tetrahedron Lett. 2013, 54, 5873–5875. [Google Scholar] [CrossRef]
  41. Dunetz, J.R.; Xiang, Y.; Baldwin, A.; Ringling, J. General and Scalable Amide Bond Formation with Epimerization-Prone Substrates Using T3P and Pyridine. Org. Lett. 2011, 13, 5048–5051. [Google Scholar] [CrossRef]
  42. Randall, A.; Scheuerman, R.; Tumelty, D. The reduction of aromatic nitro groups on solid supports using sodium hydrosulfite (Na2S2O4). Tetrahedron Lett. 2000, 41, 6531–6535. [Google Scholar] [CrossRef]
  43. Williamson, K.L. Reduction of Indigo: Sodium Hydrosulfite as a Reducing Agent. J. Chem. Educ. 1989, 66, 359. [Google Scholar] [CrossRef] [Green Version]
  44. iMODS. Available online: http://imods.chaconlab.org/ (accessed on 2 March 2021).
  45. Tetko, I.V.; Tanchuk, V.Y. Application of associative neural networks for prediction of lipophilicity in ALOGPS 2.1 program. J. Chem. Inf. Comput. Sci. 2002, 42, 1136–1145. [Google Scholar] [CrossRef]
  46. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Wallingford CT, 2009. Available online: https://gaussian.com/g09citation/ (accessed on 2 March 2021).
  47. Stewart, J.J.P. Optimization of parameters for semiempirical methods V: Modification of NDDO approximations and application to 70 elements. J. Mol. Mod. 2007, 13, 1173–1213. [Google Scholar] [CrossRef] [Green Version]
  48. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera—A visualization system for exploratory research and analysis. J. Comput. Chem. 2004, 25, 1605–1612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Yang, J.; Chen, C. GEMDOCK: A generic evolutionary method for molecular docking. Proteins 2004, 55, 288–304. [Google Scholar] [CrossRef] [PubMed]
  50. Thomsen, R.; Christensen, M.H. MolDock:  A New Technique for High-Accuracy Molecular Docking. J. Med. Chem. 2006, 49, 3315–3321. [Google Scholar] [CrossRef] [PubMed]
  51. Koes, D.R.; Camacho, C.J. ZINCPharmer: Pharmacophore search of the ZINC database. Nucleic Acids Res. 2012, 40, W409–W414. [Google Scholar] [CrossRef]
Figure 1. Indole derivatives with biological activity: indole derivative 1 is a potent HCV NS3/4A serine protease inhibitor; indole-N-acetamide inhibitor 2 is a potent inhibitor of the HCV NS5B Polymerase; umifenovir (3) has been demonstrated to inhibit in vitro infection with globally prevalent pathogenic viruses, including hepatitis C; indole derivatives 4 are the strongest known SARS-CoV-2-Mpro inhibitors with inhibition constants in the two-digit nanomolar range.
Figure 1. Indole derivatives with biological activity: indole derivative 1 is a potent HCV NS3/4A serine protease inhibitor; indole-N-acetamide inhibitor 2 is a potent inhibitor of the HCV NS5B Polymerase; umifenovir (3) has been demonstrated to inhibit in vitro infection with globally prevalent pathogenic viruses, including hepatitis C; indole derivatives 4 are the strongest known SARS-CoV-2-Mpro inhibitors with inhibition constants in the two-digit nanomolar range.
Molecules 26 03817 g001
Scheme 1. Synthesis of 2,2′-[(4-{[4-({4-[(1E)-3-(1H-indol-3-yl)-3-oxoprop-1-en-1-yl]phenoxy}acetamido)phenyl]ethynyl} phenyl) azanediyl]di(ethane-2,1-diyl)dipropionate (15).
Scheme 1. Synthesis of 2,2′-[(4-{[4-({4-[(1E)-3-(1H-indol-3-yl)-3-oxoprop-1-en-1-yl]phenoxy}acetamido)phenyl]ethynyl} phenyl) azanediyl]di(ethane-2,1-diyl)dipropionate (15).
Molecules 26 03817 sch001
Figure 2. Molecular docking between the SARS-CoV-2-Mpro and the studied ligands.
Figure 2. Molecular docking between the SARS-CoV-2-Mpro and the studied ligands.
Molecules 26 03817 g002
Figure 3. Structures of N3 inhibitor 16, 1,5-disubstituted tetrazole-1,2,3-triazoles 17 and 18.
Figure 3. Structures of N3 inhibitor 16, 1,5-disubstituted tetrazole-1,2,3-triazoles 17 and 18.
Molecules 26 03817 g003
Figure 4. Interactions between 15 and the catalytic side of the SARS-CoV-2-Mpro.
Figure 4. Interactions between 15 and the catalytic side of the SARS-CoV-2-Mpro.
Molecules 26 03817 g004
Figure 5. (A) Hydrogen bond and (B) electrostatic interactions between the SARS-CoV-2-Mpro and 15.
Figure 5. (A) Hydrogen bond and (B) electrostatic interactions between the SARS-CoV-2-Mpro and 15.
Molecules 26 03817 g005
Figure 6. Hydrophobic interactions between 15 and the SARS-CoV-2-Mpro. Blue and red surfaces depict a hydrophobic and hydrophilic behavior, respectively.
Figure 6. Hydrophobic interactions between 15 and the SARS-CoV-2-Mpro. Blue and red surfaces depict a hydrophobic and hydrophilic behavior, respectively.
Molecules 26 03817 g006
Figure 7. Indole-acetylene derivative 15 into the SARS-CoV-2-Mpro pharmacophore model. Blue spheres are the aromatic segments (Ar), green spheres depict the hydrophobic segments (Hy), gray spheres depict the hydrogen donor segments (HD), orange spheres represent the hydrogen acceptor segments (HA), and purple spheres represent the hydrophobic and aromatic segments (Hy–Ar).
Figure 7. Indole-acetylene derivative 15 into the SARS-CoV-2-Mpro pharmacophore model. Blue spheres are the aromatic segments (Ar), green spheres depict the hydrophobic segments (Hy), gray spheres depict the hydrogen donor segments (HD), orange spheres represent the hydrogen acceptor segments (HA), and purple spheres represent the hydrophobic and aromatic segments (Hy–Ar).
Molecules 26 03817 g007
Figure 8. (A) Molecular docking between the KLK5 protease and the synthesized ligand 15. (B) Main locations of KLK5 protease interacting with 15.
Figure 8. (A) Molecular docking between the KLK5 protease and the synthesized ligand 15. (B) Main locations of KLK5 protease interacting with 15.
Molecules 26 03817 g008
Figure 9. (A) Hydrogen bond and (B) electrostatic interactions between 15 and KLK5 protease.
Figure 9. (A) Hydrogen bond and (B) electrostatic interactions between 15 and KLK5 protease.
Molecules 26 03817 g009
Figure 10. Hydrophobic interactions between KLK5 protease and 15. Blue and red surfaces depict a hydrophobic and hydrophilic behavior, respectively.
Figure 10. Hydrophobic interactions between KLK5 protease and 15. Blue and red surfaces depict a hydrophobic and hydrophilic behavior, respectively.
Molecules 26 03817 g010
Table 1. Principal interaction energies between the studied ligands and the SARS-CoV-2-Mpro. E means the total interaction energy, HBond is the hydrogen bond interactions, Elstat is the electrostatic interactions, VdW represents the Van der Waals interactions, and LE is the Ligand Efficiency (LE = E/#heavy atoms), all in kcal/mol.
Table 1. Principal interaction energies between the studied ligands and the SARS-CoV-2-Mpro. E means the total interaction energy, HBond is the hydrogen bond interactions, Elstat is the electrostatic interactions, VdW represents the Van der Waals interactions, and LE is the Ligand Efficiency (LE = E/#heavy atoms), all in kcal/mol.
MoleculeEHBondElstatVdWLE
1−191.25−6.11−0.64−52.36−5.63
2−186.18−2.72−0.80−49.59−6.21
3−177.44−5.56−0.7320.92−6.12
4a−209.02−6.92−1.00−51.58−6.33
4b−221.35−3.080.24−27.07−6.51
15 (this work)−224.29−5.48−0.87−54.97−4.23
Co-Crystal (16)−198.93−4.05−1.30−53.08−4.06
17−181.75−2.52−0.65−52.43−5.19
18−255.79−7.99−0.11−60.73−5.44
Table 2. Principal interaction energies between the studied ligands and the KLK5 protein. E means the total interaction energy, HBond is the hydrogen bond interactions, Elstat is the electrostatic interaction, VdW represents the Van der Waals interactions and LE is the Ligand Efficiency (LE = E/#heavy atoms), all in kcal/mol.
Table 2. Principal interaction energies between the studied ligands and the KLK5 protein. E means the total interaction energy, HBond is the hydrogen bond interactions, Elstat is the electrostatic interaction, VdW represents the Van der Waals interactions and LE is the Ligand Efficiency (LE = E/#heavy atoms), all in kcal/mol.
MoleculeEHBondElectroVdWLE
15−168.87−4.97−1.3288.53−3.19
Tube-altre [32]−518.38−0.58−2.5521.72−2.41
Table 3. Computed LogP values and molecular weight (MW) for the studied molecules.
Table 3. Computed LogP values and molecular weight (MW) for the studied molecules.
MoleculeMWLogP a
1440.464.97
2404.514.91
3477.425.14
4a449.604.08
4b461.542.93
15711.826.70
Co-Crystal (16)680.802.11
17469.552.79
18623.683.57
a LogP calculation was performed considering the optimized molecules at level PM6 through the ALOGPS server [45].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Almaraz-Girón, M.A.; Calderón-Jaimes, E.; Carrillo, A.S.; Díaz-Cervantes, E.; Alonso, E.C.; Islas-Jácome, A.; Domínguez-Ortiz, A.; Castañón-Alonso, S.L. Search for Non-Protein Protease Inhibitors Constituted with an Indole and Acetylene Core. Molecules 2021, 26, 3817. https://doi.org/10.3390/molecules26133817

AMA Style

Almaraz-Girón MA, Calderón-Jaimes E, Carrillo AS, Díaz-Cervantes E, Alonso EC, Islas-Jácome A, Domínguez-Ortiz A, Castañón-Alonso SL. Search for Non-Protein Protease Inhibitors Constituted with an Indole and Acetylene Core. Molecules. 2021; 26(13):3817. https://doi.org/10.3390/molecules26133817

Chicago/Turabian Style

Almaraz-Girón, Marco A., Ernesto Calderón-Jaimes, Adrián Sánchez Carrillo, Erik Díaz-Cervantes, Edith Castañón Alonso, Alejandro Islas-Jácome, Armando Domínguez-Ortiz, and Sandra L. Castañón-Alonso. 2021. "Search for Non-Protein Protease Inhibitors Constituted with an Indole and Acetylene Core" Molecules 26, no. 13: 3817. https://doi.org/10.3390/molecules26133817

APA Style

Almaraz-Girón, M. A., Calderón-Jaimes, E., Carrillo, A. S., Díaz-Cervantes, E., Alonso, E. C., Islas-Jácome, A., Domínguez-Ortiz, A., & Castañón-Alonso, S. L. (2021). Search for Non-Protein Protease Inhibitors Constituted with an Indole and Acetylene Core. Molecules, 26(13), 3817. https://doi.org/10.3390/molecules26133817

Article Metrics

Back to TopTop