Search for Articles:
Title / Keyword
Author / Affiliation
Journal
Article Type

Section
Special Issue
Volume
Issue
Number
Page

Logical OperatorOperator
Search Text
Search Type

remove_circle_outline

Open AccessArticle

# In silico Investigation on the Inhibiting Role of Nicotine/Caffeine by Blocking the S Protein of SARS-CoV-2 Versus ACE2 Receptor

by 1,2, 3, 4, 1, 2, 5, 6 and 7,*
1
Department of Physics, Shahid Rajaee Teacher Training University, Lavizan, Tehran 16788-15811, Iran
2
Department of Physics, Azarbaijan Shahid Madani University, Tabriz 53714-161, Iran
3
Department of Microbiology, Faculty of Science, University of Maragheh, P.O. Box 55181-83111, Maragheh, Iran
4
Department of Genetic Engineering and Molecular Genetics, Zanjan University, Zanjan 45371-38791, Iran
5
Department of Chemical Engineering, School of Engineering and Applied Science, Khazar University, AZ1096 Baku, Azerbaijan
6
Department of Chemistry, University of Illinois, Urbana, IL 61801, USA
7
International Center for Materials Nanoarchitechtonics (WPI-MANA), National Institute for Materials Science (NIMS), 1-1 Namiki, Tsukuba, Ibaraki 305-0044, Japan
*
Author to whom correspondence should be addressed.
Microorganisms 2020, 8(10), 1600;
Received: 29 August 2020 / Revised: 13 October 2020 / Accepted: 13 October 2020 / Published: 17 October 2020

## Abstract

In this paper, we studied the in silico interaction of angiotensin-converting enzyme 2 (ACE2) human receptor with two bioactive compounds, i.e., nicotine and caffeine, via molecular dynamic (MD) simulations. The simulations reveal the efficient blocking of ACE2 by caffeine and nicotine in the exposure to the spike (S) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We have selected the two most important active sites of ACE2-S protein, i.e., 6LZG and 6VW1, which are critically responsible in the interaction of S protein to the receptor and thus, we investigated their interaction with nicotine and caffeine through MD simulations. Caffeine and nicotine are interesting structures for interactions because of their similar structure to the candidate antiviral drugs. Our results reveal that caffeine or nicotine in a specific molar ratio to 6LZG shows a very strong interaction and indicate that caffeine is more efficient in the interaction with 6LZG and further blocking of this site against S protein binding. Further, we investigated the interaction of ACE2 receptor- S protein with nicotine or caffeine when mixed with candidate or approved antiviral drugs for SARS-CoV-2 therapy. Our MD simulations suggest that the combination of caffeine with ribavirin shows a stronger interaction with 6VW1, while in case of favipiravir+nicotine, 6LZG shows potent efficacy of these interaction, proposing the potent efficacy of these combinations for blocking ACE2 receptor against SARS-CoV-2.

## 1. Introduction

The rapid pandemic spread of a novel coronavirus disease 2019 (COVID-19), officially known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is currently the most significant ongoing challenge to human health [1,2,3]. This new coronavirus can result in severe upper and lower respiratory infections, pneumonia, and kidney failure, and infected patients have a relatively high risk of mortality [4]. SARS-CoV-2 is the third highly pathogenic pneumonia coronavirus after SARS-CoV-1 and MERS-CoV (middle east respiratory syndrome-related coronavirus) to become a pandemic [5,6]. According to the World Health Organization (WHO), by 1 October 2020, over 34 million infections and 1,000,000 deaths have been reported [7]. COVID-19 has had a major impact on the normal daily lives of people throughout the world. Therefore, many researchers have been directed into his issue to overcome the pandemic outbreak of the COVID-19 virus. To fight this virus, three crucial and inevitable subjects are emphatically highlighted: developing an effective vaccine, discovery and/or synthesis of an effective drug to inhibit the progress of the virus in the infected body, and developing a fast and cheap technique with higher precision to diagnose infections.
To develop a new drug or repurpose one, it is essential to understand how a coronavirus infects the human body, how it attaches to cells, and how its progress in the host cells occurs. Looking at coronavirus structure, they are assembled as enveloped viruses with a positive-sense, single-stranded ribonucleic acid (RNA) genome [8,9,10]. SARS-CoV-2 contains four structural proteins including spike (S), envelope (E), membrane (M), and nucleocapsid (N) proteins [11,12]. It is well-known that SARS-CoV-2 engages its S protein to bind to the host cells through the protein-based receptor on the cells surface [13], called angiotensin converting enzyme 2 (ACE2). The structure of the complex between the SARS-CoV-2 S protein receptor binding domain (RBD) and the ACE2 receptor was determined to verify that the ACE2 receptor plays a notable role in mediating the entry of SARS-CoVs. ACE2 has an important role regulating the renin-angiotensin system (RAS), and it is confirmed that SARS-CoV infection reduces ACE2 expression and attenuates acute lung failure through blocking the reninâ€“angiotensin pathway [14]. Finding a compound that can block the formation of this complex or disrupt in RBD-ACE2 complex has been suggested as a reasonable strategy to come up with a rational drug discovery for COVID-19 [15]. Although the development and discovery of a novel drug takes several years [16] with heavy costs, releasing a repurposed drug to market significantly lowers the consumed cost and time. Therefore, repurposing a drug with the above strategy has been recommended to rapidly deal with the current pandemic disease [17,18]. However, although several efforts are underway, there is no promising antiviral agent released for SARS-CoV-2. Therefore, we still need to review existing drugs in order to discover their antiviral properties [19].
Nicotine is a bioactive natural alkaloid that can be found in plants and acts as stimulant and inhibitor in the biological system of humans. Although nicotine contributes to the development of cardiovascular and pulmonary disease in the long term, it has an effective influence on the homeostasis of RAS by upregulating the angiotensin-converting enzyme (ACE)/angiotensin (ANG)-II/ANG II type 1 receptor axis and downregulating the compensatory ACE2/ANG-(1-7)/Mas receptor axis [20]. It is found that serum ACE activity experiences a significant increase immediately after smoking, and it requires at least 20 min to return to normal and reach to control levels [21]. In China, a systematic review and meta-analysis studies on 5960 currently smoking patients demonstrated that pharmaceutical nicotine should be considered as a potential treatment option in COVID-19 [22]. Also, other investigations have revealed the positive effects of nicotine as a potential treatment for COVID-19 in 6515 patients and has shown that hospitalized current smokers, compared to non-smokers, had higher recovery rates. This suggests a hypothesis that nicotine may be protective against severe COVID-19, which is biologically plausible and should be further investigated through clinical trials on nicotine as a drug candidate [23]. These clinical trials should test the effects of the smoking and nicotine on the risk of being infected with COVID-19 (NCT04429815). Based on the findings regarding nicotine that have been reported so far, we further investigated nicotine and caffeine (thanks to the similar chemical structure backbones of caffeine and nicotine) through in silico calculations to understand their potency blocking the S protein binding to the ACE2 receptor and, consequently, in suppressing COVID-19 progress in an infected body.
Although there are many reports of repurposing common antiviral agents for COVID-19 treatment, interaction between RBD-ACE2 complex and caffeine and/or nicotine through molecular dynamics (MD) and molecular docking methods have not yet been considered. Computational methods based on MD and molecular docking offer versatile and cost-effective tools for the preliminary examination of any potential drugs [24,25]. The mixed ligands of caffeine/nicotine with the currently proposed antiviral agents are our further focus of investigation, which we believe merit consideration. We hypothesized that there could be a significant S protein-ACE2 blocking potential activity through forming a complex between nicotine or caffeine with other current antiviral drugs (e.g., favipiravir, ribavirin, remdesivir, chloroquine, hydroxychloroquine, oseltamivir, and valganciclovir), so we investigated it through MD and molecular docking simulations. We also evaluated molars ratios of nicotine/caffeine to RBD/CTD-ACE2, solely and together with the other antiviral agents.

## 2. Methods

The interactions of two main compounds including nicotine and caffeine with the human ACE2 receptor -S protein (i.e., 6VW1 and 6LZG) were investigated. The chemical structure backbones of nicotine and caffeine resemble some candidate anti-COVID-19 drugs (i.e., favipiravir and valganciclovir; see Figure 1). Therefore, we aimed to investigate the possible interactions and potential blocking activity of these two molecules with the CTD/RBD-ACE2, which are the epitopes of S protein and the corresponding receptor responsible for COVID-19 cell binding, respectively. This blocking behavior is intensified when the compounds are mixed with other antiviral agents. Our focus will be directed toward the potency of nicotine (and caffeine) as COVID-19 antiviral drugs by blocking the ACE2 receptor mechanism.

#### 2.1. Details of Structure ACE2 Receptors and SARS-CoV-2 Spike Protein Complexes

The several crystal structures of the complex between SARS-CoV-2 S protein and ACE2 receptor are recognized on the Protein Data Bank (PDB) (as PDB: 6LZG, 6VW1, 6M17) by X-ray diffraction and cryoelectron microscopy (cryo-EM). The 6LZG and 6VW1 structures as epitopes of the COVID-19 S protein were selected in this study. S protein in the 6LZG area has a non-chimeric nature, and the reported resolution of 6LZG is 2.50 $â„«$, which is suitable for the generation of docking [26]. The crystal structure of the SARS-CoV-2 S protein RBD in complex with ACE2 (6VW1) has a resolution of 2.68 $â„«$ [27]. The SARS-CoV-2 S protein has two main functional subunits, including S1 and S2 subunit [26]. S1 can be further divided into an N-terminal domain (NTD) and a C-terminal domain (CTD), both of which can function as a receptor-binding entity [28]. The related crystal structure of the CTD SARS-CoV-2 S protein in complex with ACE2 is reported as 6LZG [26].
The ACE2 receptor has proteins that can be further divided into N-terminal peptidase M2 domain and C-terminal collectrin domain [26]. Thus, the complex of the ACE2 receptor and the S protein of SARS-CoV-2 are connected to each other through the RBD of the S protein and the ACE2 peptidase domain (PD). Hence, a strong binding between the RBD of the S protein and the CTD region of ACE2 receptor is to be expected [29]. The crystal structure of the SARS-CoV-2 S protein and the ACE2 receptor complex were visualized by SWISS-PdB viewer 4.1.0.22 from Plan-les-Oiiates/Geneva, Switzerland [30].

#### 2.2. Selected Antiviral Agents

Several antiviral agents have been investigated for the treatment of COVID-19. Accordingly, we selected the seven antiviral drugs introduced by their chemical structure in Figure 2, and we selected the drug compounds shown in Table S1. We theoretically investigated the interaction of nicotine/caffeine with the RBD/CTD-ACE2 receptor (e.g., 6VW1 and 6LZG) when combined with these antiviral drugs to intensify the inhibition mechanism.

#### 2.3. Molecular Dynamics Simulations

Molecular docking studies provide useful information in drug discovery before the execution of any costly empirical research. For example, Durdagi et al. used MD and molecular docking to analyze the binding interactions between novel fullerene inhibitors and human immunodeficiency virus type I aspartic protease (HIV-1 PR); these complexes, as HIV-1 PR inhibitors, possess higher bioactivity [31]. In the present work, our calculations based on molecular dynamic (MD) simulations were performed using Materials Studio 7 (Accelrys, Inc., San Diego, CA, USA). Furthermore, molecular docking has been utilized for more accurate and complimentary calculations regarding the binding conditions of the ligand (nicotine/caffeine-drugs) and RBD-ACE2 complex [32]. A structural representation of the interaction between ACE2/SARS-CoV-2-CTD and caffeine and nicotine is shown in Figure 3.
The structures of the RBD-ACE2 and ACE2-CTD complexes in bulk water were built under the periodic boundary conditions for MD simulations. The packing density was determined by an experimental bulk density of 1.0 $( g / cm 3 )$.
The structures of energy minimization were implemented using an alternating sequence of steepest descent [33] and conjugate gradient [34] optimizations, with 0.001 $( kcal / mol )$ in the forcite module. For the intra and inter-molecular interactions and clusters calculations, the ab initio polymer consistence force fields (PCFF) was used [35,36,37]. The Andersen method [38] and the Berendsen method [39] were used to maintain the temperature () and the pressure, respectively. The simulation was performed using an integration time-step of 1.0 $fs$. All the dynamics runs were performed in a spline width of 1.0 $â„«$ and a buffer width of 0.5 $â„«$.
The non-bonding energy interactions were calculated within a cutoff distance of 15.5 $â„«$. The long-range electrostatic interactions were computed using the Ewald summation method (ESM) [40]. It is generally believed that the system reaches equilibrium when the parameters such as energy and temperature fluctuations fall within the range of 5%. The calculation examples are represented in Figure 4. In Figure 4aâ€“c, the total, non-bonding, and potential energies for CTD-ACE2+caffeine, RBD-ACE2+hydroxychloroquine+caffeine, and CTD-ACE2+favipiravir+nicotine, respectively, can be observed. Moreover, Figure 4dâ€“f represents the temperature profile versus the simulated time for CTD-ACE2+caffeine, RBD-ACE2+hydroxychloroquine+caffeine, and CTD-ACE2+favipiravir+nicotine structures.
For further optimizing the membrane structure, we used an annealing simulation. Here, we raised the temperature from 300 $K$ to 500 $K$ and 100 $ps$ for constant number of particles, pressure, and temperature (NPT) dynamics. In the next step, with a 250 $ps$ NPT ensemble, we decreased the temperature from 300 $K$ to 20 $K$. After relaxing the structures, NPT and constant number of particles, volume, and temperature (NVT) ensembles were performed for 500 $ps$.

#### 2.4. Interaction Energy Calculation

By eliminating the water molecules from the optimized structures, we computed the interaction energy between the RBD/CTD-ACE2 complex and the drug ($E Int ( RD )$). Next, by removing the drug molecule, the energy of the isolated S protein-receptor, $E receptor$, was computed. Finally, by deleting the receptor from the S protein-receptor+drugs ensemble, the remaining energetic term for drugs, $E Drug$, was obtained.
The generalized interaction energy (IE) between the RBD/CTD-ACE2 and the drug that describes the system [41,42] is given by
$âˆ† E I n t ( R D ) = E R e c e p t e r / D r u g âˆ’ ( E R e c e p t o r âˆ’ E D r u g ) .$

#### 2.5. Molecular Docking

It is suggested that a fast and economical tool (such as molecular docking) can be combined with molecular dynamics technology (precision but time consuming) and simulate more reliable and highly precise calculations in proteinâ€“ligand complexes. Moreover, for the fast screening of large libraries, docking can be efficiently utilized for exploring the conformations of the protein receptor, optimizing the structures of the final complexes, and calculating the accurate energies [43].
For targeting the interface of the S protein (in SARS-CoV-2) and the ACE2 receptor by a blocking agent, we evaluated the blocking potential of nicotine, caffeine, and their combined forms with other antiviral drugs. From the crystal structure reported by PDB for RBD-ACE2, the water molecules, 2-acetamido-2-deoxy-beta-D-glucopyranose, 1,2-ethanediol ligands, and the Zn and chloride ions were removed. Also, from crystal structure of ACE2/SARS-CoV-2-CTD, the water molecules, 2-acetamido-2-deoxy-beta-D-glucopyranose ligand, and the Zn ion were removed. The hydrogen atoms in these structures were enhanced in this study. We used the AutoDock v4.2 package (AutoDock is being developed and maintained in the Forli Laboratory, with support from the US National Institutes of Health) for the docking study [44,45]. Also, the charges of the molecules were applied. We selected a 60 Ã— 60 Ã— 60 $â„«$ grid box, and the distance between two grid points was set at 1.0 $â„«$ centering on the structures. In this paper, the rigid structure of the proteins was considered, so that, in this state, the drug is assumed to be fixed in shape. By using the Lamarckian genetic algorithm (LGA) [46], we performed molecular docking. In molecular docking through genetic algorithms (GA) [47], the particular arrangement of a ligand and a protein can be defined by a set of values describing the translation, orientation, and conformation of the ligand with respect to the protein. In the GA, there are some variable states in the ligandâ€™s state, and each variable state can correspond to a specific gene. Random pairs of individuals are mated using a process of crossover, in which new individuals inherit genes from either parent. In addition, some offspring undergo random mutation, in which one gene change to a random extent [46]. Thigh profile was achieved under the following conditions: an initial population of 150 randomly placed individuals and a maximum number of 2.5 Ã— 106 energy evaluations, a maximum number of 27,000 generations, a mutation rate of 0.02, a crossover rate of 0.80, and an elitism value of 1.

## 3. Results and Discussion

Here, we present our numerical results from calculating the interaction between RBD-ACE2 and CTD-ACE2 and antiviral drugs with nicotine or caffeine. We highlight two ligands (i.e., caffeine and nicotine) that we believe may be of special interest for experimental evaluation.

#### 3.1. Molar Ratios of Nicotine and Caffeine to RBD/CTD-ACE2 Receptor

In this section, we studied the effect of different molar ratios of caffeine/nicotine on the IE complex of 6LZG/6VW1+caffeine/nicotine. Our results show that using different molar rates of caffeine and nicotine led to a change of IE in both RBD-ACE2 and ACE2/SARS-CoV-2-CTD. Based on the MD simulations results of the ACE2-S protein, we predicted the binding between ligands and proteins with low IE. As shown in Figure 5a, the existence of 4 mol of caffeine in the ACE2/SARS-CoV-2-CTD complex showed the lowest IE complex. As shown this figure, the top compound that could target the binding interface between S protein and ACE2 is 4 mol of caffeine. In the case of nicotine, there are different IE trends in both 6VW1 and 6LZG, but the trends slightly change in comparison with caffeine, which are illustrated in Figure 5b. Overall, compared to nicotine, caffeine has a stronger interaction with the RBD-ACE2 receptor, especially with an ACE2/SARS-CoV-2-CTD complex. We can see that 4 mol of caffeine and the 6LZG complex do not behave similarly to 6VW1. It is therefore necessary to study how caffeine is connected with the amino acids of 6LZG and 6VW1. In Section 3.5, we will describe how amino acids are connected to caffeine in detail.

#### 3.2. Interaction Energy Complex of RBD-ACE2 and CTD-ACE2 with Nicotine or Caffeine and Drugs

In this paper, we set the structures energies as follows:
$âˆ† E B = E B o n d F i n a l âˆ’ E B o n d R e f e r e n c e$
$âˆ† E N B = E N o n âˆ’ B o n d F i n a l âˆ’ E N o n âˆ’ B o n d R e f e r e n c e$
$âˆ† E v d W = E v a n d e r W a a l s F i n a l âˆ’ E v a n d e r W a a l s R e f e r e n c e$
where $âˆ† E B$ and $âˆ† E NB$ refer to the systems bond and non-bonded intra-molecular interaction energies, respectively. $E Bond Final$ and $E Non âˆ’ Bond Final$ describe the bond and non-bonded energies of the final systems, respectively, and $E Bond Reference$ and $E Non âˆ’ Bond Reference$ are the bond and non-bonded energies of the reference receptors and drugs, respectively. $âˆ† E vdW$ is the van der Waals energy, and $E vanderWaals Final$ and $E vanderWaals Reference$ refer to the reference and final energies of the systems, respectively. The interaction energies, i.e., $âˆ† E B$, $âˆ† E NB$ and $âˆ† E vdW$, of antiviral drugs with nicotine or caffeine on the RBD/CTD-ACE2 complex were evaluated (Figure 6 and Figure 7). According to the results, the system with low IE values is more stable (i.e., when the IE values are more negative, the interaction force between RBD-ACE2 or CTD-ACE2 with the ligand is noteworthy). It is observed that the IE values of the RBD/CTD-ACE2+caffeine complex and RBD/CTD-ACE2+nicotine were found to be more stable. Also, an improvement in the IE values was observed when the antiviral drugs were combined with either nicotine or caffeine, showing that the antiviral drugs enhance the tendency of a significant interaction with the ACE2 receptor and the S protein. Creating adhesion between the receptors and ligands causes the breaking of the non-bonded intramolecular bonds in the receptor protein, and as a result, this causes a free energy penalty. Therefore, the competition between the inter- and intramolecular interactions are critical in determining the successful adhesion of ACE2 with ligand(s) (e.g., caffeine or nicotine). As illustrated in Figure 6a,b, of the 28 possible interactions (see Table S1), the IEs of favipiravir+nicotine with CTD-ACE2 and ribavirin+caffeine with RBD/CTD-ACE2 are lower than those of the other interactions. There are several possible reasons for this observation: a high contact area between the S protein-ACE2 and the combined ligands; the nucleosidic nature of favipiravir and ribavirin, which allows for better adhesion; or a high amount of hydrogen bonding can be formed between the protein and the combined ligands [48]. Notably, the matchable hydrophobicity and the complementarity shape of the receptor protein and ligands play a pivotal role in determining the blocking activity of the ligand [49].
When the difference between $âˆ† E B$ and $âˆ† E NB$ decreases, by increasing the $âˆ† E B$ and decreasing the $âˆ† E NB$, the receptor has a tendency to agglomerate. These results are due to the short distance of binding force between the receptors and/or creating hydrogen bonds between the receptors and drugs. Therefore, the most effective combinations are favipiravir with nicotine in exposure to CTD-ACE2, valganciclovir with nicotine in exposure to RBD-ACE2, favipiravir with caffeine in exposure to RBD-ACE2, and ribavirin with caffeine in exposure to RBD-ACE2, as illustrated in Figure 6c,d.
Figure 7a,b reveals that, in case of the bonding and van der Waals interactions, in exploring the combination of caffeine/nicotine with antiviral drugs in exposure to ACE2-S protein, the combination of favipiravir with caffeine has a prevailing interaction with RBD-ACE2. In the case of nicotine, valganciclovir has a better combination result in interaction strength during exposure to RBD-ACE2. Therefore, as shown in Figure 7c,d, the most efficient compounds are valganciclovir with nicotine versus RBD-ACE2, remdesivir with nicotine versus RBD-ACE2, oseltamivir with nicotine versus CTD-ACE2, favipiravir and caffeine versus RBD-ACE2, and finally ribavirin and caffeine versus RBD-ACE2.

#### 3.3. Molar Ratios of Favipiravir and Ribavirin to S Protein-ACE2 Receptor

We have reported the determination of molar ratio resonance assignments for the S protein-ACE2 in complex with ligands, e.g., favipiravir and ribavirin. As shown in Figure 8, the contribution of RBD-ACE2 dynamics can be computed for different molar ratios for favipiravir and ribavirin in the absence of nicotine and caffeine. It is observed that one mol of favipiravir+CTD-ACE2, in comparison with ribavirin+RBD-ACE2, has a strong interaction.

#### 3.4. Interaction Energy of ACE2 with Ranked Compounds

Here, in the absence of an S protein, we computed the interactions of the ACE2 receptor with the candidate drugs (e.g., favipiravir, ribavirin, and remdesivir). We aimed to observe the effect of blocking ACE2 without the S protein. Thus, the best conformations were selected on the basis of interaction energy or binding affinity between two interacting systems. Therefore, we selected the top three combinations that respond well in Section 3.2, such as ribavirin, favipiravir, and remdesivir with caffeine. Our MD simulations show that the interaction energies for ACE2+ribavirin+caffeine, ACE2+favipiravir+caffeine, and ACE2+remdesivir+caffeine were found as follows: 469.61, 843.31, and âˆ’396.68 ($kcal / mol$), respectively. As a result, ACE2+remdesivir+caffeine exhibited a very strong interaction in the absence of the S protein.

#### 3.5. Molecular Docking of the S Protein-ACE2 with Compounds

In this section, we present our molecular docking results from calculating the interaction of RBD-ACE2 and CTD-ACE2 with drug compounds and nicotine/caffeine. Basically, AutoDock computes the best binding mode using binding free energy evaluation. Also, AutoDock calculates the binding free energy using inter-molecular energy (which consists of van der Waals energy, hydrogen bonding energy, desolation energy, and electrostatic energy), internal energy of ligand, and torsional free energy [50]. The binding energy is due to the energy contributed by different amino acids or residues around the cavity of target protein on interaction with the screened molecule. These residues contribute energy is due to different interactions like hydrogen bonding, van der Waals, electrostatic interactions, Ï€-Ï€ stacking, etc. [51].
As the binding of the S protein to ACE2 is undesirable, it is preferable to diminish the ligandâ€“interface interactions that may bridge, and therefore stabilize, the interaction between the S protein and the ACE2 receptor. The strongest binding affinity, initial repurposing may be better suited to nicotine+favipiravir+ CTD-ACE2 and caffeine+ribavirin+6VW1 (Figure S1). The molecular docking calculations illustrates that the nicotine+favipiravir+ CTD-ACE2 and caffeine+ribavirin+ RBD-ACE2, equal to âˆ’7.13 and âˆ’6.76 ($kcal / mol$), respectively, have the lowest binding free energy and, hence, the strongest interactions. It seems that variation in the binding free energies occurs due to the difference in the hydrophobic interactions and hydrogen bonding formation between RBD/CTD-ACE2â€²s amino acid residues with drugs and nicotine/caffeine. The system stabilization can occur by lowering the binding free energy for the most stable conformation. The nicotine+favipiravir+ CTD-ACE2 forms hydrophobic interactions with nine amino acids, which are Ala386, Ala387, Arg393, Asn33, Gln388, Glu37, His34, Phe390, and Pro389, from the target ACE2 with nicotine and Arg403 and Tyr505 of the S protein. The hydrogen bonds created the interactions between the N groups of the nicotine with amino acid residues Arg403 and Tyr505, whereas residues Ala396, Asn394, Asn397, Asp206, Glu398, Gly395, and Lys362 of ACE2 receptor have hydrophobic interactions with favipiravir (see Table 1 and Figures S2 and S3). Also, for caffeine+ribavirin+RBD-ACE2, there are hydrophobic interactions between six amino acids (e.g., Ala386, Ala387, Arg393, Asp355, Glu37, Gly352 of ACE2) and ribavirin. There are other residues, such as Asp405, Gly502, Gly504, Tyr505, and Val503, of SARS-CoV-2 chimeric RBD that are connected with ribavirin. There are interactions between the hydrogen bonding groups of the ribavirin with amino acid residues (Lys353, Gly354, Asp350) that bind to ACE2 and amino acids, like Tyr505 or Gly504 of the S protein. The two hydrogen bonds formed with the amino acid residues are (a) the hydrogen atom of the O-H group with amino acids, including Gly504 and Lys353; and (b) the nitrogen atom of N-terminal of ribavirin with amino acids such as Asp350, Tyr505, and Gly354 (see Table 2 and Figures S4 and S5).

## 4. Conclusions

In summary, we theoretically evaluated (using MD simulations) the interactions of two crucial active sites of the S protein (i.e., 6LZG, 6VW1) when it is complexed with the ACE2 receptor using two accessible natural bioactive alkaloids, i.e., nicotine and caffeine. Then, the combination of nicotine and caffeine with antiviral drugsâ€”including favipiravir, ribavirin, remdesivir, chloroquine, hydroxychloroquine, oseltamivir, and valganciclovirâ€”were investigated as potential complimentary agents of the ACE2 receptor. The results of the MD simulations revealed a promising binding tendency in caffeine/nicotine with the ACE2 receptor; consequently, the blocking of the ACE2 receptor against SARS-CoV-2 can occur. The results indicate that the molar ratio of caffeine and nicotine with 6LZG and 6VW1 also has a major impact on the blocking activities of the ACE2 receptor. We have shown that, in the case of the 6VW1 complex, caffeine with favipiravir and ribavirin form a more permeant structure against SARS-CoV-2 in term of non-bonding interaction energy. Also, we show that the combination of nicotine and favipiravir in blocking 6LZG and the combination of caffeine and ribavirin in blocking 6VW1 were more successful. In conclusion, our results suggest that nicotine and caffeine compounds have the capability to interact with the S protein and ACE2 and interfere with their binding through blocking the active sites. The results might have significant applications in considering them as drug candidates in therapeutic treatment of SARS-CoV-2 infection.

## Supplementary Materials

The following are available online at /2076-2607/8/10/1600/s1, Figure S1: Interaction between ACE2-RBD/CTD (6VW1, 6LZG) complex with (a) 6LZG+favipiravir+nicotine and (b) 6VW1+ribavirin+caffeine that green region and magnified inlays represent active sites; Figure S2: 2D representation of favipiravir binding mode with receptor binding site of SARS-CoV-2 S protein (6LZG+favipiravir+nicotine); Figure S3: 2D representation of nicotine binding mode with receptor binding site of SARS-CoV-2 S protein (6LZG+favipiravir+nicotine); Figure S4: 2D representation of ribavirin binding mode with receptor binding site of SARS-CoV-2 S protein (6VW1+ribavirin+caffeine); Figure S5: 2D representation of caffeine binding mode with receptor binding site of SARS-CoV-2 S protein (6VW1+ribavirin+caffeine); Table S1: Selected drug compounds with S protein (RBD) and ACE2 receptor.

## Author Contributions

Conceptualization, M.H. and R.M.; Data curation, S.M., N.S. and E.D.; Investigation, S.M., A.E., M.H. and E.D.; Software, S.M.; Supervision, E.D.; Writingâ€”original draft, S.M.; Writingâ€”review & editing, M.E. (Mehrnaz Entesari), A.E., M.E. (Mohammad Esmailpour), R.M., N.S. and E.D. All authors have read and agreed to the published version of the manuscript.

## Funding

This research was supported by Shahid Rajaee Teacher Training University.

## Acknowledgments

We would like to thank Shahid Rajaee Teacher Training University for giving our group the opportunity to work. Furthermore, we would like to thank Shahid Rajaee Teacher Training University for providing cluster facility for simulation work, and other enough facilities for our group.

## Conflicts of Interest

The authors declare no conflict of interest.

## References

1. Holmes, L., Jr.; Enwere, M.; Williams, J.; Ogundele, B.; Chavan, P.; Piccoli, T.; Chinaka, C.; Comeaux, C.; Pelaez, L.; Okundaye, O.; et al. Black-White Risk Dierentials in COVID-19 (SARS-COV-2) Transmission, Mortality and Case Fatality in the United States: Translational Epidemiologic Perspective and Challenges. Int. J. Environ. Res. Public Health 2020, 17, 4322. [] []
2. Meng, L.; Hua, F.; Bian, Z. Coronavirus Disease 2019 (COVID-19): Emerging and Future Challenges for Dental and Oral Medicine. J. Dent. Res. 2020, 99, 481â€“487. [] []
3. Khalilov, R.; Hosainzadegan, M.; Eftekhari, A.; Nasibova, A.; Hasanzadeh, A.; Vahedi, P.; Hosain Zadegan, H. Overview of the environmental distribution, resistance, mortality, and genetic diversity of new coronavirus (COVID-19): Review. Adv. Biol. Earth Sci. 2020, 5, 7â€“12. []
4. Tang, N.; Bai, H.; Chen, X.; Gong, J.; Li, D.; Sun, Z. Anticoagulant treatment is associated with decreased mortality in severe coronavirus disease 2019 patients with coagulopathy. J. Thromb. Haemost. 2020, 18, 1094â€“1099. [] []
5. Hoffmann, M.; Kleine-Weber, H.; Schroeder, S.; Kruger, N.; Herrler, T.; Erichsen, S.; Schiergens, T.S.; Herrler, G.; Wu, N.H.; Nitsche, A.; et al. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell 2020, 181, 271â€“280. [] [] []
6. Letko, M.; Marzi, A.; Munster, V. Functional assessment of cell entry and receptor usage for SARS-CoV-2 and other lineage B betacoronaviruses. Nat. Microbiol. 2020, 5, 562â€“569. [] [] []
7. World Health Organization. 2020. Available online: (accessed on 27 September 2020).
8. Lu, G.; Liu, D. SARS-like virus in the Middle East: A truly bat-related coronavirus causing human diseases. Protein Cell 2012, 3, 803â€“805. [] []
9. Chen, S.J.; Wang, S.C.; Chen, Y.C. Novel Antiviral Strategies in the Treatment of COVID-19: A Review. Microorganisms 2020, 8, 1259. [] []
10. Contini, C.; Caselli, E.; Martini, F.; Maritati, M.; Torreggiani, E.; Seraceni, S.; Vesce, F.; Perri, P.; Rizzo, L.; Tognon, M. COVID-19 Is a Multifaceted Challenging Pandemic Which Needs Urgent Public Health Interventions. Microorganisms 2020, 8, 1228. [] []
11. Bosch, B.J.; van der Zee, R.; de Haan, C.A.; Rottier, P.J. The coronavirus spike 15 protein is a class I virus fusion protein: Structural and functional characterization of the fusion core complex. J. Virol. 2003, 77, 8801â€“8811. [] []
12. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; Huang, B.; Shi, W.; Lu, R.; et al. A novel coronavirus from patients with pneumonia in China, 2019. N. Engl. J. Med. 2020, 382, 727â€“733. [] []
13. Li, Q.; Kang, C.B. Progress in Developing Inhibitors of SARS-CoV-2 3C-Like Protease. Microorganisms 2020, 8, 1250. [] [] []
14. Kuba, K.; Imai, Y.; Rao, S.; Gao, H.; Guo, F.; Guan, B.; Huan, Y.; Yang, P.; Zhang, Y.; Deng, W.; et al. A crucial role of angiotensin converting enzyme 2 (ACE2) in SARS coronavirus-induced lung injury. Nat. Med. 2005, 11, 875â€“879. [] [] []
15. Smith, M.D.; Smith, J.C. Repurposing Therapeutics for COVID-19: Supercomputer-Based Docking to the SARS-CoV-2 Viral Spike Protein and Viral Spike Protein-Human ACE2 Interface. Chem Rxiv. 2020. [] []
16. Van Norman, G.A. Drugs, devices, and the FDA: Part 1: An overview of approval processes for drugs. JACC Basic Transl. Sci. 2016, 1, 170â€“179. [] []
17. Wang, C.; Horby, P.W.; Hayden, F.G.; Gao, G.F. A novel coronavirus outbreak of global health concern. Lancet 2020, 395, 370â€“373. [] []
18. Wu, C.; Liu, Y.; Yang, Y.; Zhang, P.; Zhong, W.; Wang, Y.; Wang, Q.; Xu, Y.; Li, M.; Li, X.; et al. Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods. Acta Pharm. Sin. B 2020, 10, 766â€“788. [] []
19. Sethi, A.; Bach, H. Evaluation of Current Therapies for COVID-19 Treatment. Microorganisms 2020, 8, 1097. [] []
20. Oakes, J.M.; Fuchs, R.M.; Gardner, J.D.; Lazartigues, E.; Yue, X. Nicotine and the renin-angiotensin system. Am. J. Physiol. Regul. Integr. Comp. Physiol. 2018, 315, R895â€“R906. [] []
21. Kitamura, S. Effects of Cigarette Smoking on Metabolic Events in the Lung. Environ. Health. Perspect. 1987, 72, 283â€“296. [] []
22. Farsalinos, K.; Barbouni, A.; Niaura, R. Systematic review of the prevalence of current smoking among hospitalized COVID-19 patients in China: Could nicotine be a therapeutic option? Int. Emerg. Med. 2020, 15, 845â€“852. [] []
23. Farsalinos, K.; Barbouni, A.; Poulas, K.; Polosa, R.; Caponnetto, P.; Niaura, R. Current smoking, former smoking, and adverse outcome among hospitalized COVID-19 patients: A systematic review and meta-analysis. Ther. Adv. Chronic. Dis. 2020, 11, 1â€“14. [] []
24. Ferreira, L.G.; dos Santos, R.N.; Oliva, G.; Andricopulo, A.D. Molecular Docking and Structure-Based Drug Design Strategies. Molecules 2015, 20, 13384â€“13421. [] [] []
25. Kalyaanamoorthy, S.; Phoebe Chen, Y.P. Structure-based drug design to augment hit discovery. Drug Discov. Today 2011, 16, 831â€“839. [] []
26. Hussain, M.; Jabeen, N.; Raza, F.; Shabbir, S.; Baig, A.A.; Amanullah, A.; Aziz, B. Structural variations in human ACE2 may influence its binding with SARS-CoV-2 spike protein. J. Med. Virol. 2020, 1â€“7. [] []
27. Shang, J.; Ye, G.; Shi, K.; Wan, Y.; Luo, C.; Aihara, H.; Geng, Q.; Auerbach, A.; Li, F. Structural basis of receptor recognition by SARS-CoV-2. Nature 2020, 581, 221â€“224. [] [] []
28. Li, F.; Li, W.; Farzan, M.; Harrison, S.C. Structure of SARS coronavirus spike receptor-binding domain complexed with receptor. Science 2005, 309, 1864â€“1868. [] [] []
29. Wang, Q.; Zhang, Y.; Wu, L.; Niu, S.; Song, C.; Zhang, Z.; Lu, G.; Qiao, C.; Hu, Y.; Yuen, K.-K.; et al. Structural and Functional Basis of SARS-CoV-2 Entry by Using Human ACE2. Cell 2020, 181, 894â€“904. [] []
30. Johansson, M.U.; Zoete, V.; Michielin, O.; Guex, N. Defining and searching for 70 structural motifs using DeepView/Swiss-PdbViewer. BMC Bioinform. 2012, 13, 173. [] []
31. Durdagi, S.; Mavromoustakos, T.; Chronakis, N.; Papadopoulos, M.G. Computational design of novel fullerene analogues as potential HIV-1 PR inhibitors: Analysis of the binding interactions between fullerene inhibitors and HIV-1 PR residues using 3D QSAR, molecular docking and molecular dynamics simulations. Bioorg. Med. Chem. 2008, 16, 9957â€“9974. [] [] []
32. Amaro, R.E.; Baudry, J.; Chodera, J.; Demir, O.; McCammon, J.A.; Miao, Y.; Smith, J.C. Ensemble docking in drug discovery. Biophys. J. 2018, 114, 2271â€“2278. [] []
33. Forrey, C.; Saylor, D.M.; Silverstein, J.S.; Douglas, J.F.; Davis, E.M.; Elabd, Y.A. Prediction and Validation of Diffusion Coefficients in a Model Drug Delivery System Using Microsecond Atomistic Molecular Dynamics Simulation and Vapour Sorption Analysis. Soft Matter 2014, 10, 7480â€“7494. [] []
34. Eslami, M.; Nikkhah, S.J.; Hashemianzadeh, S.M.; Seyed Sajadi, S.A. The compatibility of Tacrine molecule with poly (n-butylcyanoacrylate) and Chitosan as efficient carriers for drug delivery: A Molecular Dynamics Study. Eur. J. Pharm. Sci. 2016, 82, 79â€“85. [] []
35. Sun, H.; Mumby, S.J.; Maple, J.R.; Hagler, A.T. An ab initio CFF93 all-atom force field for polycarbonates. J. Am. Chem. Soc. 1994, 116, 2978â€“2987. [] []
36. Hwang, M.J.; Stockfisch, T.P.; Hagler, A.T. Derivation of class II force fields. 2. Derivation and characterization of a class II force field, CFF93, for the alkyl functional group and alkane molecules. J. Am. Chem. Soc. 1994, 116, 2515â€“2525. [] []
37. Sun, H. Ab initio calculations and force field development for computer simulation of polysilanes. Macromolecules 1995, 28, 701â€“712. [] []
38. Andrea, T.A.; Swope, W.C.; Andersen, H.C. The role of long ranged forces in determining the structure and properties of liquid water. J. Chem. Phys. 1983, 79, 4576â€“4584. [] []
39. Berendsen, H.J.C.; Postma, J.P.M.; van Gunsteren, W.F.; DiNola, A.; Haak, J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684â€“3690. [] []
40. Ewald, P.P. Die Berechnung optischer und elektrostatischer Gitterpotentiale. Ann. Phys. 1921, 369, 253â€“287. [] []
41. Chang, K.S.; Huang, Y.H.; Lee, K.R.; Tung, K.L. Free volume and polymeric structure analyses of aromatic polyamide membranes: A molecular simulation and experimental study. J. Membr. Sci. 2010, 354, 93â€“100. [] []
42. Li, B.; Pan, F.S.; Fang, Z.P.; Liu, L.; Jiang, Z.Y. Molecular Dynamics Simulation of Diffusion Behavior of Benzene/Water in PDMS-Calix [4] arene Hybrid Pervaporation Membranes. Ind. Eng. Chem. Res. 2008, 47, 4440â€“4447. [] []
43. Alonso, H.A.; Bliznyuk, A.A.; Gready, J.E. Combining docking and molecular dynamic simulations in drug design. Med. Res. Rev. 2006, 26, 531â€“568. [] [] []
44. Morris, G.M.; Huey, R.; Lindstrom, W.; Sanner, M.F.; Belew, R.K.; Goodsell, D.S.; Olson, A.J. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J. Comput. Chem. 2009, 30, 2785â€“2791. [] [] []
45. Morris, G.M.; Goodsell, D.S.; Halliday, R.S.; Huey, R.; Hart, W.E.; Belew, R.K.; Olson, A.J. Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function. J. Comput. Chem. 1998, 19, 1639â€“1662. [] []
46. Audipudi, A.V.; Badri, R.R.; Bhaskar, C.V.S. GC-MS and In Silico Molecular Docking Analysis of Secondary Metabolites Present in Leaf Extract of Cassia occidentalis Linn; Medicinal Plants: Biodiversity, Sustainable Utilization and Conservation; Springer: Berlin/Heidelberg, Germany, 2020; pp. 501â€“508. []
47. Holland, J.H. Adaptation in Natural and Artificial Systems; University of Michigan Press: Ann Arbor, MI, USA, 1975. []
48. Khan, R.J.; Jha, R.K.; Amera, G.M.; Jain, M.; Singh, E.; Pathak, A.; Singh, R.P.; Muthukumaran, J.; Singh, A.K. Targeting SARS-CoV-2: A systematic drug repurposing approach to identify promising inhibitors against 3C-like proteinase and 2â€²-Oribose methyltransferase. J. Biomol. Struct. Dyn. 2020. [] []
49. MacRaild, C.A.; Daranas, A.H.; Bronowska, A.; Homans, S.W. Global Changes in Local Protein Dynamics Reduce the Entropic Cost of Carbohydrate Binding in the Arabinose-binding Protein. J. Mol. Biol. 2007, 368, 822â€“832. [] []
50. Chandler, D. Interfaces and the driving force of hydrophobic assembly. Nature 2005, 437, 640â€“647. [] []
51. Kumara, D.; Kumari, K.; Jayaraj, A.; Kumar, V.; Kumar, R.V.; Dass, S.K.; Chandra, R.; Singh, P. Understanding the binding affinity of noscapines with protease of SARS-CoV-2 for COVID-19 using MD simulations at different temperatures. J. Biomol. Struct. Dyn. 2020. [] []
Figure 1. Comparing the chemical structure of (a) favipiravir with nicotine and (b) caffeine with valganciclovir; their similarities from chemical structure point of view are highlighted in yellow.
Figure 1. Comparing the chemical structure of (a) favipiravir with nicotine and (b) caffeine with valganciclovir; their similarities from chemical structure point of view are highlighted in yellow.
Figure 2. The phytochemical structures of antiviral drugs studied.
Figure 2. The phytochemical structures of antiviral drugs studied.
Figure 3. Interaction between ACE2/SARS-CoV-2-CTD (6LZG) complex with (a) nicotine and (b) caffeine. The green regions (and the corresponding magnified inlays) represent the active sites.
Figure 3. Interaction between ACE2/SARS-CoV-2-CTD (6LZG) complex with (a) nicotine and (b) caffeine. The green regions (and the corresponding magnified inlays) represent the active sites.
Figure 4. Energies and temperatures versus simulation time for (a) total, (b) non-bonded, and (c) potential energies. The plots show the temperatures for (d) CTD-ACE2+caffeine (6LZG+C), (e) RBD-ACE2+hydroxychloroquine+caffeine (6VW1+HC), and (f) CTD-ACE2+favipiravir+nicotine (6LZG+FN).
Figure 4. Energies and temperatures versus simulation time for (a) total, (b) non-bonded, and (c) potential energies. The plots show the temperatures for (d) CTD-ACE2+caffeine (6LZG+C), (e) RBD-ACE2+hydroxychloroquine+caffeine (6VW1+HC), and (f) CTD-ACE2+favipiravir+nicotine (6LZG+FN).
Figure 5. The interaction energies (IEs) for different ratios of (a) caffeine and (b) nicotine structures with 6LZG and 6VW1.
Figure 5. The interaction energies (IEs) for different ratios of (a) caffeine and (b) nicotine structures with 6LZG and 6VW1.
Figure 6. The interaction and non-bonded energies for drug and caffeine/nicotine structures with 6VW1 (RBD-ACE2) and 6LZG (CTD-ACE2). (a,b) The interaction energies caffeine and nicotine, (c,d) the non-bonded energies caffeine and nicotine, respectively.
Figure 6. The interaction and non-bonded energies for drug and caffeine/nicotine structures with 6VW1 (RBD-ACE2) and 6LZG (CTD-ACE2). (a,b) The interaction energies caffeine and nicotine, (c,d) the non-bonded energies caffeine and nicotine, respectively.
Figure 7. The bond and van der Waals energies for drug and caffeine/nicotine structures with 6VW1 (RBD-ACE2) and 6LZG (CTD-ACE2). (a,b) The bond energies caffeine and nicotine, (c,d) the van der Waals energies caffeine and nicotine, respectively.
Figure 7. The bond and van der Waals energies for drug and caffeine/nicotine structures with 6VW1 (RBD-ACE2) and 6LZG (CTD-ACE2). (a,b) The bond energies caffeine and nicotine, (c,d) the van der Waals energies caffeine and nicotine, respectively.
Figure 8. The interaction energies for different ratios of the structures of favipiravir and ribavirin with ACE2/SARS-CoV-2-CTD (6LZG) and RBD-ACE2 (6VW1). The 6LZG+F and 6VW1+R display 6LZG+favipiravir and 6VW1+ribavirin, respectively.
Figure 8. The interaction energies for different ratios of the structures of favipiravir and ribavirin with ACE2/SARS-CoV-2-CTD (6LZG) and RBD-ACE2 (6VW1). The 6LZG+F and 6VW1+R display 6LZG+favipiravir and 6VW1+ribavirin, respectively.
Table 1. Contact residues of the SARS-CoV-2 CTD-ACE2 interfacing with favipiravir and nicotine. The residues involved in the hydrogen and hydrophobic interactions between S protein and ACE2 with favipiravir+nicotine.
Table 1. Contact residues of the SARS-CoV-2 CTD-ACE2 interfacing with favipiravir and nicotine. The residues involved in the hydrogen and hydrophobic interactions between S protein and ACE2 with favipiravir+nicotine.
FavipiravirNicotine
S proteinACE2S proteinACE2
Hydrogen bonding -Arg403(N)-Tyr505(N)-
Hydrophoic interactions Ala396-Asn394-Asn397-Asp206-Glu398-Gly395-Lys362 Ala386-Ala387-Arg393-Asn33-Gln388-Glu37-His34-Phe390-Pro389
Table 2. Contact residues of the SARS-CoV-2 RBD-ACE2 interfacing with ribavirin and caffeine. Residues involved in the hydrogen and hydrophobic interactions between S protein and ACE2+ribavirin+caffeine.
Table 2. Contact residues of the SARS-CoV-2 RBD-ACE2 interfacing with ribavirin and caffeine. Residues involved in the hydrogen and hydrophobic interactions between S protein and ACE2+ribavirin+caffeine.
RibavirinCaffeine
S proteinACE2S proteinACE2
Hydrogen bondingGly504(O, O-H)-Tyr505(N)Lys353(O)-Asp350(NH2)-Gly354(N)--
Hydrophoic interactionsAsp405-Gly502-Val503Ala386-Ala387-Arg393-Asp355-Glu37-Gly352-Leu351-Met383-Phe356-Thr324 Arg393-Asn394- Asp350-Glu37-Gly352-Leu351-Leu391-Leu392-Phe40-Phe390
 Publisherâ€™s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.