Sunil Kumar1*, Jaya Gupta1, Awanish Kumar2*, Ankur Mohan3
1Faculty of Biotechnology, Institute of Bioscience and Technology, Shri Ramswaroop Memorial University, India
2Department of Biotechnology, National Institute of Technology, Raipur, Chhattisgarh, India
3MRD-LifeSciencesTM Lucknow, Uttar Pradesh
*Corresponding authors: Sunil Kumar, Faculty of Biotechnology, Institute of Bioscience and Technology, Shri Ramswaroop Memorial University, Barabanki (UP), India. Tel: +917379893670; Email: email@example.com; firstname.lastname@example.org
Received Date: 09 August, 2017; Accepted Date: 01 September, 2017; Published Date: 08 September, 2017
recent years, it has become increasingly evident that pathological neoangiogenesis
processes occur in chronic inflammation and Chronic Liver Diseases (CLD). The
challenge for the upcoming years is the characterization of the molecular basis
and pathways of angiogenic disorders in an integrated manner. If the currently
uncertain Patho-genetic role of angiogenesis in CLD can be demonstrated, it
could have a prognostic value in the evaluation of disease progression. The current work involves the
identification of the protein CPEB4, which is the main cause of CLD&
We have analyzed its role in the disease by using bioinformatics tools. The sequence structure and functional analysis of the concerned protein has been studied through BLAST and Phyre2.0. Since the precise structure of CPEB4 protein is not yet available, therefore, theoretical model for the query protein was generated using Easy Modeler 4.0. The protein-ligand docking was done with use of Auto dock software by drugs downloaded from Drug Bank. Ligand N-[4-(6-Carbamoyl-7-methoxy-4-quinolyl) oxy-2-trifluoromethylphenyl]-N'-methylurea showed an impressive binding energy of -14.3 kcal/mol, in comparison with the already proved drugs Nateglinide & Lenvatinib (binding energy: -8.8 kcal/mol& -8 kcal/mol respectively). Overall result suggests that this inhibitor could represent a promising molecule for the effective treatment of CRC.
Keywords: Angiogenesis; Chronic Liver Diseases; CPEB4 protein Inhibitor Therapy
Chronic Liver Diseases (CLDs) are one of the major causes of death, which are still increasing year-on-year . Therefore, knowledge about the pathophysiology of CLDs and its complications is of uttermost importance. There are clear role of angiogenesis in the disease progression of various liver diseases [2,3]. Looking closer at the Pathophysiology of Portal Hypertension (PH) , fibrosis [5,6], cirrhosis , Non-Alcoholic Steato Hepatitis (NASH)  and Hepatocellular Carcinoma (HCC) . We find that angiogenesis is a recurring factor in the disease progression. The close relationship between the progression of CLDs and angiogenesis emphasizes the need for anti-angiogenic therapy as a tool for blocking or slowing down the disease progression . The fact that angiogenesis plays pivotal role in CLD give rise to new opportunities for treating CLDs and its complications . Angiogenesis in liver is characterized by intrahepatic vascular remodeling with capillarization of the sinusoids and development of intrahepatic shunts, which lead to increased hepatic resistance and decreased effective hepatocytes perfusion . It has become increasingly clear that angiogenesis occurring during chronic wound healing and fibro genesis provides a key contribution to disease progression and complications .
The progression from fibrosis to cirrhosis, the end-point of CLDs, is distinguished by a prolonged inflammatory and fibrogenic process that leads to an abnormal angioarchitecture distinctive for cirrhosis . Several mechanisms are responsible for the angiogenic switch during the pathogenesis of CLDs [10,11]. First, CLDs are characterized by chronic inflammation . Increased intrahepatic vascular resistance is primarily caused by anatomical changes, such as fibrotic scar tissue compressing portal and central venues . In addition, the formation of fibrotic septa, as well as sinusoidal capillarization, can result in an increased resistance to blood flow and oxygen delivery . This results in hypoxia and the transcription of hypoxia-sensitive pro-angiogenic genes, usually modulated trough HIF . Also, an increased contribution of the hepatic artery to the sinusoidal blood flow leads to sensitization of the hepatocytes to abnormal high-oxygen conditions .
Hepatic Stellate Cell motility and migration promote coverage of HSC around sinusoids, causing sinusoidal constriction and contributing to the hepatic resistance in cirrhosis . The characteristic fenestrated phenotype of the sinusoidal ECs is lost and an organized basement membrane is established, which leads to an impairment of oxygen diffusion even though the increased arterial flow provides a high supply of oxygenated blood . The deposition of collagen in the space of Disseaccentuates the narrowing and distortion of the sinusoidal lumen, further restricting micro vascular blood flow . This is aggravated by leucocytes, either mechanically trapped in the narrowed sinusoids or adhering to the endothelium, as a result of activation of a hepatic micro vascular inflammatory response . The hypoxic liver tissue causes an up-regulation of VEGF, Ang and their receptors in HSCs, not only enhancing the hypoxia-induced angiogenesis but also stimulating activation and migration of HSC . As stated previously, activated HSCs induce an inflammatory response and enhance angiogenesis . New vessels themselves can significantly contribute to the perpetuation of the inflammatory response by expressing chemokines and adhesion molecules promoting the recruitment of inflammatory cells . Furthermore, angiogenesis, early in the course of a CLD, may contribute to the transition from acute to chronic inflammation . This work explores the potential usefulness of angiogenesis inhibitor N-[4-(6-Carbamoyl-7-methoxy-4-quinolyl)-oxy-2-trifluoromethylphenyl]-N'-methylureain the treatment of CLD using in silico platform.
2.1. Sequence Alignment
The FASTA sequence of CPEB4 query protein was obtained from online source of National Centre for Biotechnology Information, NCBI (http://www.ncbi.nlm.nih.gov) and aligned by using BLASTp sequence alignment tool (http://blast.ncbi.nlm.nih.gov/Blast.cgi). Chain a of H1N1 influenza A virus (PDB ID: 4F15) has showed maximum identity with the query protein. Hence, it has been selected as a template protein for structure prediction of the targeted protein.
2.2 Structure Prediction
The three-dimensional secondary structure of CPEB4 protein has been predicted with the help of online server 3DPSSM and Phyre2. Ten protein templates have been selected showing maximum identity with our query protein structure. The PDB file format of the selected templates have been obtained from Protein Data Bank (PDB), a computational database, having resolution <3.0, R-Value <0.5 and X-Ray crystallographic structures. With the help of the software EASY MODELLER (version 4.0), the 3-dimensional model of query protein has been created by submitting PDB file formats of these template protein to the software program.
2.3. Validation of Predicted Model (Loop Modelling)
Structural Analysis and Verification Server, SAVES (Prochek) is an online bioinformatics tool which has been used in our work for validation of the query protein model (Ramachandran Plot Analysis) of the query protein. The model query protein has been rebuilt into its secondary structure by further submitting it into the RAMPAGE. The final secondary structure obtained after validation has been viewed using the Swiss PDB Viewer (SPDB 4.10).
2.4. Prediction of the Active Site
The query protein, having 704 amino acid residues, has been submitted to online Dog Site Scorer, an Active site prediction and Analysis server. Prediction of active site is necessary to identify the ligand binding site, present on the query protein.
2.5. Drug Library Search
Drug Bank database has been used to obtain the structure of drug molecules (in MOL format) already available for the treatment of swine influenza. The MOL files of these structures were then converted into PDB format with the help of software Open-Babel. These molecules were used as ligand during molecular docking .
2.6. Molecular Docking
Molecular docking has been performed using Auto Dock software against the active site of the query protein with the available ligand . The ligand molecule (that showed minimum free energy) has been used as a base molecule for designing in silico structure of drug molecule against Swine Influenza. A computational pipeline was developed for large-scale molecular docking of drugs to protein targets. Briefly, we collected all 3D structures available for each drug target, determined binding pockets in the structures, and docked drugs to each pocket. Results were collected and thresholds were applied to select the top predicted interactions, which were then visually inspected.
2.7. Finding Similar Drug Molecule of Best Docked Ligand
With the help of Pubchem server, the in-silico drug molecules were obtained. The molecule that showed minimum free energy when binding with the active site of the query protein has been selected as the final drug molecule.
2.8. Toxicity Prediction
The toxicity of the selected ligand has been checked with the help of an online server ACD Labs I Labs. This online server gives the information regarding the physical and chemical properties of any molecule by following the five rules of Lipinsky, Health effects and ADME properties. Quantitative Structure-Activity Relationships (QSAR) is quantitative models which relate the variation in measures of activity in a series of chemical compounds to the variation in chemical structure between compounds in the series. The use of QSAR in the study and identification of potentially toxic chemicals is receiving increased attention, as demonstrated by several recent reviews (for example, McKinney, 1985; Hansch, 1985; Enslein, 1984; Birge and Cassidy, 1983).
2.9. Energy Minimization of Protein Drug Complex
Energy was minimized by Gromacs software on Linux operating system of the protein-drug complex structure, so that our complex will be strong enough to be stable. Minimization provides information that is complementary to that obtained from molecular dynamics or Monte Carlo. Ensembles of dynamics or Monte Carlo structures are useful for calculating thermodynamic averages and estimating entropy, but the large number of structures involved makes detailed microscopic analysis cumbersome. Minimized structures in some sense represent the underlying configurations about which fluctuations occur during dynamics and, as such, provide a convenient and meaningful basis for structural analysis (Hagler, 1985; Struthers et al., 1984).
The FASTA sequence of CPEB4 (Accession No.: AAI43960) obtained from NCBI has showed maximum identity of 32% with Chain A of Crystal Structure of a Heterogeneous Nuclear Ribo nucleoprotein A/b (hnrpab) From Homo Sapiens At 2.15 A Resolution (PDB: 3S7R) (Figure 1).
By using the secondary structure, 3-dimensional tertiary structure of the query protein has been created by using software Modeler. The obtained three-dimensional computational model of query protein has been allowed for validation by following Ramachandran Plot analysis, using RAMPAGE, an online server. It has been viewed under Swiss PDB Viewer (Figure 3-5).
The active site of the query protein obtained from Dog Site Scorer (Figure 6), has been allowed to molecular docking and a minimum free energy of (binding energy: -8.8 kcal/mol& -8 kcal/mol, Nateglinide Lenvatinib simultaneously) (Figure 7) selected as a base molecule for further synonyms drugs. The novel drug molecule obtained by us from pub chem server, possess a free energy of -14.3 kcal/mol. IUPAC name of the novel drug stands N-[4-(6-Carbamoyl-7-methoxy-4-quinolyl)oxy-2-trifluoromethylphenyl]-N'-methylurea (Figure 8).
The obtained molecule is the most suitable with consideration of the Lipinsky’s rule for the toxicity of the chemical compounds in vivo (Table 1).
Therefore, this molecule has been chosen for the best possible drug against angiogenesis and chronic liver disease. Energy was minimized by Gromacs software on Linux operating system of the protein-drug complex structure (Figure 9) so that our complex will be strong enough to be stable.
The in-silico identification of the biological targets of small organic molecules has attracted substantial attention from biologists, chemists and pharmacologists. In the present work we have found that angiogenesis and chronic liver disease is caused by CPEB4 protein. Through homology modeling and loop modeling we have got the final structure of our protein and after its docking with the best result, from Pubchem server, we have obtained a suitable ligand having minimum free energy when docked with the active site of the query protein. In vivo toxicity of the novel drug has been checked. The de-novo structure based drug designing approach has led to a potential ligand (drug) molecule that blocks the active site of CPEB4 protein preventing the responsible gene of angiogenesis and chronic liver disease [18,19]. Energy was minimized by Gromacs software on Linux operating system of the protein-drug complex structure so that our complex will be strong enough to be stable. Our in-silico analysis shows that N-[4-(6-Carbamoyl-7-methoxy-4-quinolyl) oxy-2-trifluoromethylphenyl]-N'-methylurea could offer an improved treatment CLD patients. Further prospective pre-clinical and clinical studies are still needed to confirm our findings. The active site of the query protein obtained from Dog Site Scorer (Figure 6), has been allowed to molecular docking and a minimum free energy of (binding energy: -8.8 kcal/mol& -8 kcal/mol, Nateglinide& Lenvatinib simultaneously) (Figure 7) selected as a base molecule for further synonyms drugs. The novel drug molecule obtained by us from Pubchem server, possess a free energy of -14.3 kcal/mol. IUPAC name of the novel drug stands N-[4-(6-Carbamoyl-7-methoxy-4-quinolyl)oxy-2-trifluoromethylphenyl]-N'-methylurea.
Authors thank Shri Ramswaroop Memorial University, Barabanki, Uttar Pradesh, India; MRD-Life SciencesTM, Luck now, Uttar Pradesh, India and National Institute of Technology, Raipur (CG), India for their continuous support and assistance during the course of research work and scientific writing.
Figure 1: Maximum Identity value of CPEB4 after running BLASTp.
Figure 2: Minimum Dope Score showed by the template protein molecules after Homology Modeling and Running Easy Modeler.
Figure 3: Validation with RAMPAGE Server.
Figure4: Ramachandran Plot Analysis with RAMPAGE.
Figure 5: Three-dimensional model of CPEB4 protein viewed under Swiss PDB Viewer.
Figure 6: Active Site of Query Protein obtained by Dog Site Scorer (in yellow colour).
Figure 7: Chemical Structures of Nateglinide & Lenvatinib simultaneously.
Figure 8: Chemical Structure of N-[4-(6-Carbamoyl-7-methoxy-4-quinolyl) oxy-2-trifluoromethylphenyl]-N'-methylurea.
Figure 9: GROMACS version in the 5.x series.
123.90 0.3 cm3
No. of Hydrogen Bond Donors
No. of Hydrogen Bond Acceptors
Table 1: Various chemical parameters of N-[4-(6-Carbamoyl-7-methoxy-4-quinolyl) oxy-2-trifluoromethylphenyl]-N'-methylurea.
13. Lopes Baitello ME, Tenani GD, Ferreira RF, Victor Nogueira, Souza Pinhel MA, et al. (2016) VEGF Polymorphisms Related to Higher Serum Levels of Protein Identify Patients with Hepatocellular Carcinoma. Can J Gastroenterol Hepatol 2016: 9607054.
14. Han X, Zhu X, Hong Z, Wei L, Ren Y, et al. (2017) Structure-based Rational Design of Novel Inhibitors against Fructose-1,6-Bisphosphate Aldolase from Candida albicans. J Chem Inf Model 57: 1426-1438.
Citation: Kumar S, Gupta J, Kumar A, Mohan A (2017) Structure Prediction and Drug Designing of a Target Protein CPEB4 on In-Silico Platform for Angiogenesis and Chronic Liver Disease. J Pharmacovigil Pharm Ther: JPPT-119. DOI: 10.29011/JPPT-119. 100119