DRUG DESIGN OF CORTICOTROPHIN RELEASING FACTOR RECEPTOR-1 ANTAGONIST FROM DERIVATIVES OF α ASARONE BY MOLECULAR DOCKING
HTML Full TextDRUG DESIGN OF CORTICOTROPHIN RELEASING FACTOR RECEPTOR-1 ANTAGONIST FROM DERIVATIVES OF α ASARONE BY MOLECULAR DOCKING
S. Suba Malani * and A. Sasi Kumar
Department of Physiology, Tagore Dental College & Hospital, Rathinamangalam, Chennai - 600127, Tamil Nadu, India.
ABSTRACT: Stress is one of the most significant health problem in modern society. This paves way for the researchers to focus and investigate into the biological pathways linking stress and health. Hence the goal of the present study is to discover novel drug derivatives from the active principle α asarone of the plant Acorus calamus Linn. (ACL) which was found to have antioxidant &antidepressor effect the basic physiological response to stress is through activation of Corticotrophin –releasing factor Receptor 1 (CRFR-1) which is predominantly found in pituitary gland is mainly responsible for signalling the adrenal gland receptor thereby increasing the Glucocorticoids as a response to stress. The 3D structure of CRFR 1 was derived from RCSB-PDB database and the α asarone was derive from Pubchem database. A total of 100 ligands were derived from α asarone using in silico method by ACD chemsketch software. Rapid virtual screening of these compounds were performed in the dock in tool iGEMDOCK v2.0. Based on the binding energy a total of seven ligands were further docked using Auto Dock 4.0 software. The selected seven ligands were then analysed for the drug relevant properties based on “Lipinski’s rule of five” and other drug like properties. Results of the current study showed that one of the derivatives of α asarone namely 4[(hyydroxyethynyl)amino]phenol-methane have shown excellent binding energy valves of -100.26 kcal/mol with a good drug likeliness and drug score. Hence this compound is a excellent drug candidate in the inhibition of CRFR-1 and can be potentially used to manage the elevated levels of glucocorticoids level during stress.
Keywords: |
CRFR-1,α asarone, Molecular docking, In-silico method
INTRODUCTION: Stress is one of the most significant health problem in modern society as it is the major factor or the stimuli for the etiology of various diseases the threatens organism’s survival. Exposure to stressful condition shows its ill effect in the cellular homeostasis.
The brain has evolved adaptive mechanism for coping with stress and responds to stressor in highly stereotyped ways by producing free radicals which result in oxidative stress. Besides the locus coeruleus – noradrenalin / sympathetic nervous system 1 a major physiological response to stress is through the Hypothalamic pituitary adrenal axis. The paraventricular nucleus (PVN) of the hypothalamus synthesize and secrete corticotrophin-releasing factor (CRF), which will act on the anterior pituitary gland to produce Adernocorticotropic hormone(ACTH) which will further act on the adrenal cortex for the secretion of glucocorticoids from Zona Fasciculate 2.
The central response to stress via HPA axis is mediated through the corticotrophin – releasing factor (CRF) and its receptor CRFR -1 which is found majorly in the anterior pituitary and certain other parts of brain 3. Binding of CRF to the type 1 receptors of anterior pituitary causes activation of adenylyl cyclase in the cell membrane and increase formation of cAMP in the cell cytoplasm in turn activates the intracellular enzymes resulting in production of ACTH 4. Studies showed that CRFR-1 mutant mice produced a statistically significant reduction in anxiety level 5, 6 and it is an evident that the binding of the antistressor or antidepressor drugs with the CRFR-1 is important in the regulation of stress. Hence the study was designed to derive an inhibitors for CRFR-1.
Free radicals are the resultant of aerobic metabolism which consist of reactive oxygen species ROS and reactive nitrogen species RNS 7. Antioxidants are the only solution to decrease oxidative stress and minimize the incidence of various pathology caused by oxidants.
During stress the antioxidants level in the body normally could not neutralize the free radicals. Administration of antioxidants as drug one way to avoid the damage induced by the free radicals or another approach is to administer drugs that either boost antioxidants or reduce free radicals generation.Some native Indian plants said to have natural antioxidant properties. Acorus calamus Linn. (ACL) “sweet flag” is a well known medicinal plant. The various extract of the plant like rhizome, roots, and aromatic oils has been used for various medicinal purposes. The ACL is rich in alkaloids, phenolics and flavonoids which are found to be responsible for its beneficial biological activities 8, 9. Isomeric forms of asarone are α and β, the α asarone has been reported to be major bioactive constituent on free radicals scavengers and lipid peroxidation.
This scavenging free radical of ACL has been found to be beneficial in overcoming excess production of oxygen free radicals generated due to continous exposure of stress 10. The plant Acorus calamus Linn. mentioned in Indian literature; but not much scientific work to show the mechanism of action in stress relevant receptors like CRFR-1. In silico is a computerized technique which is widely used in the medical research field to discover newer drugs without using animal model. The novel drug with its affinity with the receptors, the absorption, metabolism, excreation and toxic profile can be studied by this in silico method 11.The link between computational and experimental strategies has been great values in identification and development of novel promising compounds.
MATERIALS AND METHODS:
Target preparation and active site prediction: The three dimensional structure of CRF type 1 receptor was found from the database RCSB-PDB. (http://www.rcsb.org/ pdb/ explore/ explore). The PDB code is 4k5y. The active site or the possible binding sites of CRF type 1receptor (CRFR-1) was found using binding site predictor 3D LIGAND SITE which is an online tool 12.
Generation and Optimization of Ligand: The active compound alpha asarone which is present in Acorous calamus Linn. is selected for the present study and the SDF format of the structure was downloaded from Pub chem database (which is shown in Table 1). Further it is converted into MDL format in Open Babel software.
TABLE 1: THE PUBCHEM ID, IUPAC NAME AND STRUCTURE OF α ASARONE
Compound Id | IUPAC Name | Structure of α- asarone |
CID 636822 | 1,2,4-trimethoxy-5-[(1E)-prop-1-en-1-yl]benzene
(C12H16O3) |
![]() |
Using ACD Chem Sketch totally 100 ligands in 2D Format were derived with IUPAC names and the compounds were stored in MOL format. The ligands which are in MOL format was converted into PDB format for docking in Open Babel Software. ACD ChemSketch is a latest and powerful software in which we can draw desire chemical molecule, its IUPAC names and can be stored in desired formats 13. A software tool called iGEMDOCK v2.0 was used to perform rapid virtual screening of the 100 compounds 14. A population Size of 150 is set with 70 generation and one solution for quick docking. After completion of docking seven ligands were selected as they had a good binding energy. The drug relevant properties and drug likeliness of the seven ligands was analysed on “Lipinski’s rule of five” 15& OSIRIS Property Explorer. (http://www. organicchemistry.org/prog/peo/). Finally the seven new derivatives were taken for further accurate docking.
Molecular docking: The seven compounds were subjected to accurate or very slow docking by setting population size of 800 with 80 generation and 10 solutions. The post docking results of iGEMDOCK was analysed based on their electrostatic profile, hydrogen-bonding, and Van der Waals interaction.
RESULT & DISCUSSION: The corticotrophin – releasing factor receptor CRFR -1 and its 3D structure which was downloaded from RCSB-PDB is further viewed in RASMOL Structure scheme. (as shown in Fig. 1). In the Rasmol structure the alpha helices are shown in magenta, Beta sheets in yellow, turns are in paleblue colous and rest of the residues in white colour.
FIG. 1: 3D STRUCTURE OF CRFR1USINGRASMOL STRUCTURE COLOUR SCHEME
The 100 compounds were derived from the active principle α asarone using ACD Chemsketch and it was converted in to PDB format in Open Babel Software. Virtual rapis screening in iGEMDOCK software was performed for the100 compounds and finally seven best compounds were indentified on the basis of total binding energy (kcal/mol), Vanderwaals force and H bond acceptor. The results of seven best derivatives in accurate iGEMDOCK scores were displayed in Table 2. The seven selected derivatives were subjected for
Lipinski’s properties which is called as Lipinski’s Rule of five. The rule states that the compound is said to be poor in absorption or permeation if the compound has the following property: there are more than 5 H – bond donors (nHDon), more than 10 H-bond acceptors (N+O), molecular weight (MW) is over 500, Moriguchi’s log P(X LogP) is over 4.15 15, 16. The results of Lipinski’s properties were depicted in Table 3. As the final step the drug relevant properties of the seven compounds were done and the results in Table 4.
TABLE 2: THE RESULTS OF iGEMDOCK (ACCURATE) SHOWING BINDING ENERGIES OF THE SEVEN DERIVATIVES
S.no | Ligand | Total binding energy (kcal/mol) | Vanderwaals force | H bond |
1. | 1,2-diethoxy-4-methoxy-5-[(1E)-prop-1-en-1-yl]benzene | -94.7552 | -87.7552 | -7 |
2. | ({4-(chloromethoxy)-5-methoxy-2-
[(1E)-prop-1-en-1-yl]phenoxy}methyl)sodium |
-76.1521 | -76.1521
|
-10.3323 |
3. | 1-(fluoromethoxy)-2,4-dimethoxy-5-
[(1E)-prop-1-en-1-yl]benzene |
-80.3066 | -69.9743
|
-17.3323
|
4. | [{2-methoxy-5-[(1E)-3-phosphanylprop-1-en-1-yl]benzene 1,4diyl}bis(oxymethanediyl)]diphosphae | -83.592 | -76.592
|
-10.3328
|
5. | 4,5-dimethoxy-2-[(1E)-prop-1-en-1-yl]phenol | -93.816 | -77.0741 | -10.5 |
6 | [(1E)-2-chloro-1-(2,4,5-trimethoxyphenyl)
prop-1-en-1yl] potassium |
-81.1627 | -70.6627
|
-12.5
|
7. | 4[(hyydroxyethynyl)amino]phenol-methane | -100.26 | -58.8591 | -17.9436 |
TABLE 3: THE RESULT OF THE LIPINSKI’S PROPERTIES OF THE SELECTED SEVEN DERIVATIVES
S.no | Ligand | Molecular
weight |
Xlog p | H bond
donor |
H bond
Acceptor |
1. | 1,2-diethoxy-4-methoxy-5-[(1E)-prop-1-en-1-yl]benzene | 126.147 | -0.747 | 2 | 4 |
2. | ({4-(chloromethoxy)-5-methoxy-2-[(1E)
-prop-1-en-1-yl]phenoxy}methyl)sodium |
140.128 | -0.045 | 2 | 4 |
3. | 1-(fluoromethoxy)-2,4-dimethoxy-5-
[(1E)-prop-1-en-1-yl]benzene |
150.124 | -0.368 | 4 | 5 |
4. | [{2-methoxy-5-[(1E)-3-phosphanylprop-1-en-1-yl]benzene-1,4diyl}bis(oxymethanediyl)]diphosphane | 160.301 | -0.365 | 4 | 5 |
5. | 4,5-dimethoxy-2-[(1E)-prop-1-en-1-yl]phenol | 170.632 | -0.348 | 2 | 3 |
6 | [(1E)-2-chloro-1-(2,4,5-trimethoxyphenyl)
prop-1-en-1yl]potassium |
156.163 | -0.659 | 3 | 2 |
7. | 4[(hyydroxyethynyl)amino]phenol-methane | 188.145 | -1.165 | 2 | 5 |
TABLE 3: THE DRUG RELEVANT PROPERTIES OF SELECTED LIGANDS
S. no | Ligand | Drug
likeness |
Drug score | Mutagenic
|
Tumorigenic | Irritant
|
1. | 1,2-diethoxy-4-methoxy-5-[(1E)
-prop-1-en-1-yl]benzene |
0.42 | 0.7 | NO | NO | NO |
2. | ({4-(chloromethoxy)-5-methoxy-2-[(1E)
-prop-1-en 1yl]phenoxy}methyl)sodium |
0.26 | 0.5 | NO | NO | NO |
3. | 1-(fluoromethoxy)-2,4-dimethoxy-5-[(1E)
-prop-1-en-1-yl]benzene |
0.52 | 0.5 | NO | NO | NO |
4. | [{2-methoxy-5-[(1E)-3-phosphanylprop-1-en-1-yl]benzene 1,4diyl}bis(oxymethanediyl)]diphosphane | 0.24 | 0.6 | NO | NO | NO |
5. | 4,5-dimethoxy-2-[(1E)-prop-1-en-1-yl]phenol | 0.16 | 0.52 | NO | NO | NO |
6 | [(1E)-2-chloro-1-(2,4,5-trimethoxyphenyl)
prop-1-en-1yl]potassium |
0.27 | 0.5 | NO | NO | NO |
7. | 4[(hyydroxyethynyl)amino]phenol-methane | 2.03 | 0.92 | NO | NO | NO |
The results of the drug relevant properties of the seven ligands shows that the seventh ligand 4[(hydroxyethynyl) amino] phenol - methane (structure shown in Fig. 2), was found to have a good drug likeness of 2.03 and drug score of 0.92. all the ligands does not possess to have mutagenic, tumorigenic, irritant profile. The docking pose of 4[(hydroxyethynyl) amino] phenol-methane with CRFR-1 is shown in Fig. 3. The seventh ligand 4[(hydroxyethynyl) amino] phenol-methane which was derived from α-asarone using ACD chemsketch is proved to fall into the Lipinsk’s rules of five and have sufficient drug likeliness and acceptable ADME properties (absorption, distribution, metabolism, excretion) and no toxicity profile so that it can be taken into consideration for further clinical trials.
FIG. 2: STRUCTURE OF 4[(HYYDROXYEHYNYL) AMINO] PHENOL-METHANE AND ITS TOXICITY PROFILE, DRUG LIKENESS, DRUG SCORE USING OSIRIS SOFTWARE
FIG. 3: DOCKING POSE OF CRFR1 WITH (4[(HYYDROXYEHYNYL) AMINO] PHENOL-METHANE) AN EFFECTIVE DERIVATIVE FROM Α ASARONE
CONCLUSION: Stress involves activation of the Hypothalamic - Pituitary - Adrenal axis. This activation consists of release of the stress neuropeptide corticotropin-releasing hormone (CRH) from hypothalamicneurons and activation of CRH receptors (CRHRs) within the pituitary. The resulting release of adrenocorticotropic hormone leads to secretion of corticosteroids from the adrenal gland into the circulation. The release of glucocorticoids has its ill effects in cellular level producing alteration in cellular homeostasis like damage to mitochondria. Endoplasmic reticulum and further more. These changes in the cellular structure is basic pathogenis of various degenerative diseases. Indeed, it is know that blocking CRHR signaling will reduce stress as it will not increase ACTH secretion. Moreover, the current study shows that the derivate of α asarone of Acorous calamus has got an effective inhibitors activity on the of the CRH receptor type 1 (CRFR1) In the present study, with the help of molecular docking studies, a novel compound 4[(hyydroxyethynyl)amino] phenol-methane is found to have the ability to inhibit CRFR1 an important protein in regulating stress in turn gaining cellular homeostasis.
ACKNOWLEDGEMENT: The authors are thankful to the department of Bio informatics, Tagore Medical College and Hospital. Dr. Chitraa R. Chandran, Principal Tagore Dental College & Hospital, Chennai. Dr. Ramaswamy, Professor, Department of Physiology, Saveetha Medical College and Hospital.
CONFLICT OF INTEREST: Conflict of interest declared none.
REFERENCES:
- Robert M. Berne, Matthew N. Levy, Text Book of Physiology, 5th Elsevier publishers; 2004. p. 915-16.
- Kim E. Barrelt, Scott Boitano, Susan M. Barman, Heddwen L brooks, Ganong’s Review of Medical Physiology, 23rd p.353-54
- Sean M. Smith, Wylie W. Vale: The role of the hypothamic-pituitary-adrenal axis in neuroendocrine responses to stress. Dialogues Clin Neurosci. 2006; 4: 383-395.
- Guyton & Hall, Text book of Medical Physiology, 11th Elsevier publishers; p. 955-56.
- Smith GW., Aubry JM., Dellu F: Corticotropin releasing factor receptor 1 -deficient mice display decreased anxiety, impaired stress response, and aberrant neuroendocrine development. 1998; 20: 1093–1102.
- Timpl P., Spanagel R., Sillaber: Impaired stress response and reduced anxiety in mice lacking a functional corticotrophin - releasing hormone receptor 1. Nat Genet.1998; 19: 162–66.
- Chandanachoudhury, Baru, suparnasen: An comparative study of the in vitro antioxidant property of different extracts of ACL. Scholars research library 2014; 1:8-18
- Balakumbahan, K. Rajamani. Acorus calamus: An overview. Journal of Medicinal Plants Research. 2010; 4:2740-2745,
- Arasan Elaya Raja, M. Vijayalakshmi: Acorus calamus Linn: Chemistry and Biology. Research J. Pharm and Tech 2009; 2: 256-261.
- [Manikandan S, R.S. Devi, .Antioxidant property of alpha asarone against noise stress induced changes in different region of rat brain: Pharmacological research 2005; 52: 467-474.
- Kannan I, Abiramavalli R: Design of Cholera Toxin Antagonists by Molecular Docking from Polyphenolic Compounds of Tea, International Journal of PharmTech Research. 2015; 7: 54-63.
- Wass MN, Kelley LA and Sternberg MJ, 3DLigandSite: predicting ligand-binding sites using similar structures, NAR. 2010; 38: 469-473.
- ACD/Chem Sketch Freeware, version 10.00, Advanced Chemistry Development, Inc., Toronto, ON, Canada, 2012.
- Yang JM and Chen CC, "Gemdock: A generic evolutionary method for molecular
- " Proteins: Structure, Function, and Bioinfor.2004; 288-304.
- Lipinski C.A., Lombardo F., Dominy B.W and Feeney P.J: Advanced Drug Delivery Reviews 2001; 46: 3-26.
- Lipinski C.A., Lombardo F., Dominy B.W and Feeney P.J: Advanced Drug Delivery Reviews. 1997; 46: 3-26.
How to cite this article:
Malani SS and Kumar AS: Drug design of corticotrophin releasing factor receptor-1 antagonist from derivatives of α asarone by molecular docking. Int J Pharm Sci Res 2017; 8(11): 4704-09.doi: 10.13040/IJPSR.0975-8232.8(11).4704-09.
All © 2013 are reserved by International Journal of Pharmaceutical Sciences and Research. This Journal licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License.
Article Information
24
4704-4709
536
947
English
IJPSR
S. S. Malani * and A. S. Kumar
Department of Physiology, Tagore Dental College & Hospital, Rathinamangalam, Chennai, Tamil Nadu, India.
suba_malani@tagoremch.com
22 March, 2017
23 May, 2017
27 May, 2017
10.13040/IJPSR.0975-8232.8(11).4704-09
01 November, 2017