Activity of phytochemical constituents of Curcuma longa (turmeric) and Andrographis paniculata against coronavirus (COVID-19): an in silico approach

Background In early 2020, many scientists are rushing to discover novel drugs and vaccines against the coronavirus, and treatments for COVID-19, because coronavirus disease 2019 (COVID-19), a life-threatening viral disease, affected first in China and quickly spread throughout the world. In this article, in silico studies have been performed to explore the binding modes of chemical constituents for natural remedies like Curcuma longa (turmeric) and Andrographis paniculata against COVID-19 (PDB ID 5R82) targeting coronavirus using Schrodinger suit 2019-4. The molecular docking studies are performed by the Glide module, in silico ADMET screening was performed by the QikProp module, and binding energy of ligands was calculated using the Prime MM-GB/SA module. Results The chemical constituents from turmeric like cyclocurcumin and curcumin and from Andrographis paniculata like andrographolide and dihydroxy dimethoxy flavone are significantly binding with the active site of SARS CoV-2 main protease with Glide score more than − 6 when compared to the currently used drugs hydroxychloroquine (− 5.47) and nelfinavir (− 5.93). When compared to remdesivir (− 6.38), cyclocurcumin from turmeric is significantly more active. The docking results of the compounds exhibited similar mode of interactions with SARS CoV-2. Main protease and the residues THR24, THR25, THR26, LEU27, SER46, MET49, HIE41, GLN189, ARG188, ASP187, MET165, HIE164, PHE181, and THR54 play a crucial role in binding with ligands. Conclusion Based on in silico investigations, the chemical constituents from turmeric like cyclocurcumin and curcumin and from Andrographis paniculata like andrographolide and dihydroxy dimethoxy flavone, significantly binding with the active site of SARS CoV-2 main protease, may produce significant activity and be useful for further development.


Background
Coronavirus disease 2019 (COVID-19) is a life-threatening disease which was affected first in China and quickly spread throughout the world [1][2][3][4][5][6]. According to the WHO data, as of the second week of April 2020, there are 21.5 lakhs peoples in the world affected by COVID-19, out of these more than 1.5 lakhs peoples died. With more asymptomatic infections being found among COVID-19 cases, it is worthy of consideration the detailed current evidence and understanding of the transmission of SARS CoV, MERS-CoV, and SARS CoV-2 and discussion on pathogen inactivation methods on coronaviruses is very important [7][8][9][10][11][12].
In this emergency situation, it is very difficult to discover novel drugs with all clinical trials and also determine the side effects, adverse effects, etc. So, it is important to find some natural remedies for the prevention and treatment of COVID-19. From the literatures, the natural products like Curcuma longa (turmeric) and Andrographis  paniculata were reported for various biological activities and used traditionally for curing many diseases. Also, there is no or minimum side effects reported when compared to allopathic drugs. The dried and powdered root Curcuma longa (turmeric) is belonging to the Zingiberaceae family, which is being cultivated in many countries worldwide. It has many uses such as textile dyes, herbal medicines, or food products. The biological properties of its chemical components were reported for inhibition of platelet aggregation [13], anti-diabetic [14], anti-tumor [15][16][17], anti-inflammatory effects [18], antioxidant effects [19], anti-platelet aggregation effects [20], gastroprotective effects [21], lipid-lowering effects [22], Alzheimer's effects [23], etc.
As part of our ongoing research on searching the potent biological molecules against various diseases by in silico and wet lab methods [34][35][36][37][38][39][40][41][42][43][44], we have designed and evaluated various heterocyclic compounds for their biological activities. Using different modules (Glide, QikProp, and Prime) of Schrödinger suite LLC various computational methods like molecular docking, ADMET screening, and bindingfree energy, calculations were performed to find the interactions responsible for SARS CoV-2 main protease inhibition. These studies will provide the requirement of key structural features in the design of potential drug candidates.

Methods
The 3D crystal structure of COVID-19 protein called SARS CoV-2 main protease receptor co-crystallized with 6-(ethylamino) pyridine-3-carbonitrile (PDB ID 5R82, resolution 1.31 Å) was retrieved from the protein data bank. The protein was prepared using the protein preparation wizard of epic module of Schrödinger suite 2019-4. The protein structure retrieved from the RCSB protein data bank is a monomer with co-crystallized ligand. The protein was prepared by using the protein preparation wizard by refining bond orders, addition of hydrogens, and deleting water molecules beyond 5 Å, and missing chains are included by using the Prime module [45] of Schrödinger suite 2019-4. Protein minimization was performed using optimized potentials for liquid simulations (OPLS3) molecular force field with RMSD of crystallographic heavy atoms kept at 0.30 Å. A grid box was generated to define the centroid of the active site. All the compounds were docked into the catalytic pocket of SARS CoV-2 main protease by using the Glide module of Schrödinger suite 2019-4 in extra precision (XP) mode [46]. The ligands with significant Glide scores have more binding affinity towards SARS CoV-2 main protease enzyme. To predict the free energy of binding for the set of ligands in complex with a receptor, post-docking energy minimization studies were performed using Prime molecular mechanics-generalized Born surface area (MM-GB/SA) of Schrödinger 2019-4. The energy for minimized XP docked pose of ligand-receptor complex was calculated using the OPLS3 force field and generalized Born/surface area (GB/SA) continuum VSGB 2.0 solvent model [47,48].

Results
Results are summarized in Tables 1, 2, and 3 and Figs. 1,  2, 3, 4, 5, 6, and 7. The results revealed that the SARS CoV-2 main protease inhibitory property of the compounds isolated from some natural products like Curcuma longa (turmeric) and Andrographis paniculata greatly depended on the chemical nature of the substituents. The chemical structures of selected major bioactive constituents of Curcuma longa (turmeric) and Andrographis paniculata are given in Fig. 1a and b. The anti-malarial drug which was currently recommended in many countries like the USA, India, etc. [49] for the treatment of COVID-19 is hydroxychloroquine (Fig. 1c).
The docking studies of the ligands to protein active sites were performed by an advanced molecular docking program Glide module of Schrodinger suite 2019 Maestro-12.2 version for determining the binding affinities of the compounds. The designed analogues were docked towards the SARS CoV-2 main protease (PDB ID 5R82) in order to ascertain their inhibitory activity. The analogues show best fit root mean square difference (RMSD) value of 0.2.
The results are summarized in Table 1. Almost all the compounds are docked in the same binding pocket.
From the molecular docking study, it was revealed that the ligands have shown agreeable Glide G score values from − 6.13 kcal/mol (T1_Curcumin) to − 6.77 kcal/mol (T4_Cyclocurcumin) when compared to the currently recommended drugs for COVID-19 hydroxychloroquine (G score − 5.47) and nelfinavir (− 5.93). When compared to remdesivir (− 6.38), cyclocurcumin from turmeric is significantly more active. From the obtained binding modes, it was illustrated that the ligands formed hydrophobic interactions and hydrogen bonding interactions with different residues THR24 to GLN192 surrounding the active pocket which was shown in Fig. 4. The ligand N1_Andrographolide exhibited hydrogen bonding interaction with some amino acid residues and with some water molecules which are shown in Fig. 5. The presence of aromatic features and different heterocyclic rings majorly contributed towards lipophilic factors (Fig. 6).
The Glide score of the standard compound hydroxychloroquine was decreased because of the rotational penalty of the side alkyl chain which was shown in Fig. 7.
Molecular docking was additionally assessed with MM-GBSA free restricting vitality which is identified with the post-scoring approach for SARS CoV-2 main protease (PDB ID 5R82) target and the values are shown in Table 3.

Discussion
From the docking results, as shown in Table 1, it is clearly demonstrated that some of the chemical constituents from turmeric like cyclocurcumin and curcumin and from Andrographis paniculata like andrographolide and dihydroxy dimethoxy flavone significantly bind with the active site of COVID-19 main protease with Glide score more than − 6 when compared to the currently recommended drug hydroxychloroquine (G score − 5.47) and significantly inhibit SARS CoV-2 main protease and may be active against COVID-19 on further process. The above compounds have good affinity to the receptor due to more lipophilic character and also due to hydrogen bonding. From the 2D-ligand interaction diagrams, almost all the compounds exhibited similar mode of interactions with SARS CoV-2 main protease and the residues THR24, THR25, THR26, LEU27, SER46, MET49, HIE41, GLN189, ARG188, ASP187, MET165, HIE164, PHE181, and THR54 play a crucial role in binding with ligands. From Fig. 5, the docking score of the ligand N1_Andrographolide is increased due to hydrogen bonding interaction with SER46 (H-bond length 3.45 Å), GLY143 (Hbond length 2.72 Å), and THR25 (H-bond length 1.90 Å) residues and with some water molecules. From Fig. 6, it is clearly demonstrated that most of the aromatic features are covered in the lipophilic region (red color) which contributed towards lipophilic factors.
From Fig. 7, the Glide score of the standard hydroxychloroquine is decreased because of the rotational

Conclusion
From the results of the docking study, the chemical constituents of Curcuma longa (turmeric) and Andrographis paniculata demonstrated better arrangement at a dynamic site. The in silico structuring strategy embraced in the present investigation helped for recognizing some lead molecules and furthermore may somewhat clarify their useful impact for further determinations like in vitro and in vivo assessments. Results from the in silico study exhibited that many of the chemical constituents from Curcuma longa (turmeric) and Andrographis paniculata family may be useful against COVID-19 by inhibiting SARS CoV-2 main protease enzyme. Based on in silico studies, the chemical constituents such as cyclocurcumin and curcumin from turmeric and andrographolide and dihydroxy dimethoxy flavone from Andrographis paniculata are significantly active against COVID-19 by inhibiting SARS CoV-2 main protease enzyme with remedial possibilities and are probably going to be helpful after further refinement. In conclusion, consuming turmeric in our diet regularly may be a useful remedy in the prevention of the coronavirus.