Molecular docking analysis of phytochemicals from ethanolic extract of crescentia cujete with the auto inhibited parkin catalytic domain

The autoinhibited Parkin catalytic domain (PDB ID: 4BM9) receptor has been described to have a role in the ubiquitination of α-syn in Parkinson's disease. Therefore, it is of interest to discuss the molecular docking analysis data of phytochemicals from ethanolic extract of Crescentia cujete with the auto inhibited Parkin catalytic domain. We report the docking features of the phytochemical named 1, 2-Ethanediamine, N-(2-aminoethyl) with the target protein for further consideration towards the design and development of anti-Parkinson agents.


Background:
Neurodegeneration is a process, which involves neuropathological condition, and brain aging due to lack of physical movement and mental relaxations lead to all the stress related disorders [1]. The most second common neurodegenerative disorder is Parkinson disease, which is characterized by four main symptoms slowness of movement, muscle rigor, involuntary and not controllable shaking and loss of postural balance along with some secondary manifestations like a decline in memory, language, problem-solving and other thinking skills that affect a person's ability to perform everyday activities, soft speech and difficulty in swallowing due to uncoordinated movements of mouth with the throat [2]. It is characterized by the loss of pigmented dopaminergic neurons in the substantia nigral pars compacta (SNpc) of the midbrain and the presence of lewy bodies and occurs due to inhibition of mitochondrial complex-1 different mechanisms of cell damage like 190 ©Biomedical Informatics (2020) excitotoxicity, calcium homeostasis, apoptosis, protein aggregation and interaction between genetic and environmental factors [3,4]. Several  The available drugs are levodopa, carbidopa, apomorphine, amantadine, orphenadrine, benzhexol, benztropine, selegiline, pergola and many more. These drugs are effectively reversing the symptoms of Parkinson and improve the level of dopamine. The greatest disadvantage at present is available potent synthetic drugs lies in their adverse effects like constipation, ulcer, respiratory depression, hypertension, toxicity and reappearance of symptoms after discontinuation [13]. To overcome the adverse effects, a safe and effective alternate should be developed, as phytotherapeutic agents having neuroprotective effects might be necessary with minimal or null toxicity are formulated in PD treatment. The medicinal plants contain several chemical components of therapeutic value so they can be used as drugs or formulations to treat neurodegenerative diseases.
Crescentia cujeteis a small tree belongs to the family Binoniacea, which grows about 6-10m tall with a wide crown and long branches covered with clusters of tripinnate leaves and gourd-like fruit and these branches are arranged as simple elliptical leaves clustered at the anode. C.cujete leaves are used to treat cold, bronchitis, cough, asthma, urethritis and also to cure haematomas and tumours [14]. The leaves of C.cujete plants have high potential for development as a natural antioxidant and used as a central nervous system depressant in traditional medicine [15]. Therefore, it is of interest to discuss the molecular docking analysis data of phytochemicals from ethanolic extract of Crescentia cujete with the auto inhibited parkin catalytic domain.

Methods: Plant collection:
The leaves of Crescentia cujete were collected from different localities of Coimbatore District and authenticated by the Botanical Survey of India (BSI), Southern Regional Centre, Tamilnadu Agricultural University campus, Coimbatore. A voucher specimen (No: BSI/SRC/5/23/2017/Tech 2021) has been deposited at the Herbarium of the Botany department.

Preparation of plant extracts:
The leaves were cleaned and shade dried for 7 days, then grounded well to a fine powder. About 600 g of dry powder was extracted with ethanol (80%) at 70°C by continuous hot percolation using soxhlet apparatus. The extraction was continued for 2 days. The ethanolic extract was then filtered and kept in a hot air oven at 40°C for 1day to evaporate the ethanol from it. A greenish brown residue was obtained.

Gas chromatography-mass spectrometry (GC-MS) analysis:
Ethanolic extract of Crescentiacujetewas lyophilized and the sample was subjected to GC-MS equipment Thermo MS DSQII (Thermo Fisher Scientific, USA) for the analysis. The equipment has a DB 35 -MS Capillary Standard non-polar column with dimensions of 30 mm ×0.25 mm ID ×0.25 µm films. The carrier gas used is Helium with at a low of 1.0 ml/min. The injector was operated at 250°C and the oven temperature was programmed as follows: 60ºC for 15 min, then gradually increased to 280°C at 3 min. The identification of components was based on Willey and NIST libraries as well as the comparison of their retention indices [16]. Identification of Active site residues Active site residues for the protein structure 4bm9 was 191 ©Biomedical Informatics (2020) identified by using CASTp server [17]. The residues are SER145,  LEU162, ARG163, VAL164, ALA172, THR173, LEU174, THR175,  LEU176, THR177, GLU178, GLY179, HIE227, GLU426, LYS427,  ASN428, LYS431, HIS433, MET434, LYS435, TRP453, ARG455,  MET458, GLY459, TRP462, and VAL465 respectively.

Preparation of Protein:
The crystal structure of the autoinhibited Parkin catalytic domain was downloaded from Protein Data Bank (PDB) with the PDB ID: 4BM9. The protein preparation was carried out using the protein preparation wizard of Schrodinger Suite 2018. In general, the protein preparation consists of fixing structures, deleting unwanted chains and waters, fixing hetero groups and finally optimizing the fixed structure. The 3D structure of the autoinhibited Parkin catalytic domain (PDB ID: 4BM9) was shown in (Figure 1).  (Figure 2 and 3).

Results and Discussion:
Our study is focusing on the inhibition of Parkin catalytic domain. We want to enhance inhibitory activity of autoinhibited parkin catalytic domain, so we are not considering the active form of protein. GlideScore is an empirical scoring function that approximates the ligand binding free energy. It has many terms, including force field (electrostatic, van der Waals) contributions and terms rewarding or penalizing interactions known to influence ligand binding. It has been optimized for docking accuracy, database enrichment, and binding affinity prediction and be used to rank poses of different ligands, for example in virtual screening.
As it simulates a binding free energy, more negative values represent tighter binders. The Crescentia cujete compounds that were used for the In-silico study are listed (Figure 3). Molecular interactions of the ligands to 4bm9 were ranked based on various parameters such as Gscore, DockScore and Hbonds. From the docking results, six phytocompounds from Crescentia cujete showed better interaction with 4bm9 (Table 2). Docking Analysis of phytocompounds identified from ethanolic extract of Crescentia cujete are listed in (Table 3). The lead compound has obtained the highest Gscore of -7.41with hydrogen bonds. The DockScore of lead compound was also found to be -7.37. The 4bm9 and lead ligand complex showed in (Figure 4 A and B) (2020) atoms. The interaction of lead ligand with 4bm9 suggests that it is enhancing the inhibitory activity of parkin catalytic domain for overcoming the disease and makes it beneficiary for identifying novel lead molecule to cure the disease. The Gscore for rest of the five compounds 2 5-dimethyloxazolidine, Xycaine (3676), 1-(4-(2methoxyethyl) phenoxy)-3-(N-methyl-N-isopropylamino)propan-2ol, Azidophenylaceto amide and 5,8-Dibromo-7-methoxy-3methoxycarbonylpyrimido[1,6-a]indole were -4.5, -4.96, -5.24, -6.63 and -4.28 respectively.Similarly for the complex 4bm9 with 5,8-Dibromo-7-methoxy-3methoxycarbonylpyrimido[1,6-a]indole, the interacting atom in the protein was found to be oxygen (two atoms) and that of the ligand molecule was Hydrogen (two atoms) and the interacting residues were LYS 299, HIE 227 with the bond length of 2.08, 2.08 respectively. The complex 4bm9 with Azidophenylacetoamide showed interaction with three amino acid residues LEU 228, VAL 465 and HIE227 and the bond lengths were 2.14, 1.50 and 1.85 respectively. The interacting atom of protein was oxygen (three atoms) and that of ligand molecule was Hydrogen (three atoms). 4bm9 and 1-(4-(2-methoxyethyl) phenoxy) showed interaction with three amino acid residues GLN 252, VAL 465 and GLN 252 with the bond lengths 1.94, 1.74 and 2.25 where the interacting atoms were 3 oxygens in protein and 3 hydrogens in the ligand. Two interacting residues at the VAL 465 and GLN 252 were found with the protein complexed with Xycaine. The interacting atom in protein was oxygen (two atoms) and that of the ligand here is also hydrogen (two atoms) and the bond lengths were 1.59 and 2.09 respectively. In 4bm9 and 2, 5-dimethyloxazolidine complex, two oxygen atoms in the protein and two Hydrogen atoms in the ligand showed interaction with the amino acid VAL 465 with the bond length 2.08 for both the atoms.

Conclusion:
We report the docking features of the phytochemical named 1,2-Ethanediamine, N-(2-aminoethyl) with the target protein for further consideration towards the design and development of antiparkinson agents. The lead compound satisfying the Lipinski rule of 5 and ADMET properties, to be act as a drug.