Molecular docking analysis of Cianidanol fromGinkgo biloba with HER2+ breast cancer target

HER2 is a known therapeutic target for about 30% of breast cancer patients where HER2 is over expressed and this is referred to as HER2 positive breast cancer. This subtype is characterized by a clinical behavior know to be especially aggressive. Improved HER2 targeting agents such as trastuzumab, pertuzumb, lapatinib and ado-trastuzumab emtansine are available. Some patients have shown no response to treatment while others show progress to these agents. Therefore, it is of interest to screen HER2+ with phyto-chemical lead compound from Ginkgo biloba using molecular docking techniques. We screened 25 phyto-chemicals from literature with HER2+. Results show that cianidanol have an acceptable binding energy of (-8.2kcal/mol). Thus, we report the binding properties of cianidanol with HER2+.


Background:
Human Epidermal Growth Factor Receptor type 2 (HER2) belongs to the family of human epidermal growth factor receptors (HER/EGFR/ERBB), which also includes HER1, HER3, and HER4. They are a special family of oncogenic proteins whose amplification has been shown to play important roles in the development and progression of the certain aggressive type of breast cancer. Recently, HER2 has become an important biomarker and therapeutic target for about 30% of cases of breast cancer in patients. [1] Despite recent breakthroughs, breast cancer remains the most prevalent type of cancer in women and the second most deadly disease in advanced countries [2]. HER2 is of particular interest in breast cancer because, in one of its subtype (amounting to about 15 -20 % of all cases), it is over-expressed, giving it the name HER2+ breast cancer. This subtype is particularly characterized by a clinical behavior known to be especially aggressive [3]. In about 50 % of cases where this receptor is over-expressed, estrogen receptor (ER) and/or progesterone receptor (PR) is also over-expressed [4]. While considering this type of breast carcinoma, the development of targeted therapies has considerably improved its prognosis. Although advanced cases are still mostly considered incurable and there is a wide variation in survival among patients [5] Hitherto, there are no known activating ligands of HER2 directly implicated in triggering its signaling cascade, but there are reports that the receptor is activated by homo-and heterodimerization with other known HER family receptors. These include the HER3, HER4 and HER1 [6,7]. Dimerization of HER2 and HER3 has been shown to be the most potent oncogenic pair, leading to activation of their kinases with subsequent activation of downstream signaling pathways which includes the phosphatidylinositol-3-kinase (PI3K)/Akt and Mitogen-Activated Protein Kinase (MAPK) [6,7,8,9]. Consequently, the activated signaling cascade promotes expression of more oncogenes that are involved in orchestrating the survival of tumor cells, their proliferation and differentiation, thus fostering tumor progression [9, 10]. Accurately assessing HER2 status is crucial in decision making involving treatment for patients with breast cancer. A false-negative status for HER2 may lead to decisions where anti-HER2-directed therapy is omitted as much as a falsepositive status may bring about unnecessary administration of cost-ineffective and extended treatment with no known benefits [11].
Research involving HER2 targeting agents has progressed over the years with discoveries showing this to be one of the most productive research interests in oncogenic drug development. Trastuzumab targeted therapy as well as agents such as pertuzumab, lapatinib and ado-trastuzumab emtansine (T-DM1) is the standard treatment for HER2+ breast cancer patients. Although there have been several insights into improved HER2 targeting and considerable efficacy shown by target agents, some patients have shown no response to treatment with others eventually progressing [12]. These insights include inhibition of the HER family dimerization; which is an established mechanism of anti-HER2 therapy resistance by patients [13], delivering HER2-targeted chemotherapy where antibody-drug conjugate T-DM1 (a combination of trastuzumab and emtansine) is used to deliver cytotoxic therapy directly to breast cancer cells [14], targeting HER2 / ER crosstalk which has been shown to promote tumor resistance in HER2+ / ER+ tumors [15] and the use of mutated HER2 as a target in HER2 non-amplified breast cancer in patients where somatic mutation in HER 2 gene is responsible for the activation of HER2 signaling pathway independent of dimerization of receptors [16]. Ginkgo biloba nut has been used in traditional Chinese medicine for the management and treatment of varied medical conditions such as asthma and cough. It is also being used in the negation and treatment of brain, systemic circulatory disorders, and Alzheimer's disease [17]. In addendum, pharmacological properties exhibited by Ginkgo biloba phytochemicalsinclude cell cycle regulatory, antioxidant, antiproliferative, anti-angiogenic and antiestrogenic activities [18].
Despite the great efficacy shown by drugs in the clinic against HER2 positive breast cancer, the reported resistance in patients with long-term trastuzumab treatment must be overcome and there is a need to identify and develop a novel therapeutic agent that can decrease the amount of HER2 in breast cancer cells which at the same time will possess a molecular mechanism of pharmacological activity that will subdue the resistance mechanisms employed against known agents.

Methodology: Ligand selection and preparation
The chemical structures of twenty-five phytochemicals were obtained from the database of PubChem compounds (https://pubchem.ncbi.nlm.nih.gov).The downloaded MOL SDF format of these ligands was converted to PDBQT file using PyRx tool to generate atomic coordinates and energy was minimized by using the optimization algorithm at force field set at mff (required) on PyRx.

Accession and preparation of the target protein
The HER2 protein was prepared by recovering the threedimensional crystal structure of HER2 (PDB: 5o4g) in a complex with a cocrystallized ligand from RCSB PDB (Protein Data Base). (http://www.rcsb.org/pdb/home/home.do). Removing the bound complex molecule, non-essential water molecules and all heteroatoms using the Pymol toolkit then cleaned the protein.
The co-crystallized ligand was extracted (not removed) from the active site to reveal the coordinate of the grid around the binding pocket when viewed on the pymol. . This image was generated using Discovery studio software.

Molecular docking using PyRx
After preparation of the receptor and ligands, molecular anchor (docking) analysis was done by PyRx, AutoDockVina option based on the notation functions. For our analysis, we used the exhaustive search anchoring function of PyRx, AutoDock Vina. After the minimization process, the grid resolution was centred at 76.6051 × 90.6294 × 55.7857 along x, y, and z-axes, respectively, for a grid size of 25x25x25 Å to defined the binding site. The cocrystallized binder (ligand), which serves as the default was first anchored on the linking (binding) site of HER2, and the resulting interaction was compared to that cianidanol in similar active sites using the same grid box dimension.

Validation of docking results
The results obtained were validated during the blasting of the FASTA sequence of the HER2 crystalline structure (ID: 5o4g) obtained from the protein database of ChEMBL (Www.ebi .ac.uk / chembl /). The bioactivity generated by the database, with an activity of 64, an IC50 value of 1206 and the KI value of 178 was downloaded in txt format. Missing or lost data was removed only 30 of the 1206 drug-related compounds were recovered. The compiled compounds were split and converted to 2D (in sdf format) by the DataWarrior software (version 2) and converted to pdbqt format by the PyRx tool. The binders were anchored in the HER2 binding domain using the PyRx AutoDockVina logging function. A correlation coefficient was plotted between the coupling scores of 30 generated compounds and their corresponding PCHEMBL_VALUE values (determined experimentally). The graph of the correlation coefficient of Spearman Rank was plotted to obtain the correlation (R²) between the ChEMBL compounds and their corresponding results generated experimentally.

Result and Discussion:
Human Epidermal Growth Factor Receptor type 2 (HER2), belongs to the family of oncogenic proteins whose inhibition of amplification or over-expression has been shown to be cancer target (Figure 1).
The reliability of our docking scores was validated using the online ChEMBL database, the Fasta sequence of the HER2 crystalline structure (ID: 5o4g) obtained was BLAST on www.ebi.ac.uk/chembl/. The binding site of HER was docked with the compounds obtained from the search, and a correlation coefficient graph was generated by plotting the CHEMBL's Pchem values (experimentally determined IC50) and the docking scores of the compounds obtained from the search. From the plot a strong correlation coefficient of (R 2 = 0.75) was obtained (Figure  3), this gives certitude to the verity and reliability of the computational experiment and that PyRx Auto DockVina algorithm is dependable ( Table 1).

Conclusion:
Docking studies and ADMET evaluation of cianidanol showed that this ligand is drug-gable and plays critical role in the inhibition of HER2. It could be deduced that cianidanol could service as a potential antagonistic agent against HER2 + , which is overexpressed in aggressive female breast cancer.