Submit manuscript...
MOJ
eISSN: 2379-6294

Toxicology

Research Article Volume 4 Issue 5

Computer aided drug designing for proteasome inhibitor to find out potential drug candidate for multiple myloma

Vinita Srivastava,1 Shashi Prabha Agrawal,2 Swinder Jeet Singh Kalra,3 Alka Dubey4

1Associate Professor from Department of Chemistry, D.A-V, College Kanpur INDIA
2Assistant Professor from Department of Zoology, D.A-V College Kanpur, INDIA
3Associate Professor from Department of Chemistry, D.A-V, College Kanpur INDIA
4Visiting Research Associate, in Bioinformatics Infrastructure Facility (Funded by Dept. of Biotechnology, Govt. of India) of Forest Research Institute, Dehradun, INDIA

Correspondence: Alka Dubey, Visiting Research Associate, in Bioinformatics Infrastructure Facility (Funded by Dept. of Biotechnology, Govt. of India) of Forest Research Institute, Dehradun, India

Received: September 12, 2018 | Published: October 23, 2018

Citation: Srivastava V, Agrawal SP, Kalra SJS, et al. Computer aided drug designing for proteasome inhibitor to find out potential drug candidate for multiple myloma. MOJ Toxicol. 2018;4(5):352–356. DOI: 10.15406/mojt.2018.04.00126

Download PDF

Abstract

Block the action of proteosome via Proteasome inhibitors. This inhibitors act as drug candidate for blocking the malfunctioning of proteosome. In market very few drugs available for the multiple myloma treatment in this finding we update very famous drug bortezomib for multiple myloma treatment via computational platforms and we found this drug bind via protein ligand interaction with favorable statical & structural (Tables Below) computation platforms.

Keywords: Proteasome inhibitors, bortezomib, multiple myloma, ligand

Introduction

The incorporation of the first proteasome inhibitor, bortezomib, into anti-myeloma armamentarium can be considered a major milestone in treatment of multiple myeloma, greatly improving the response rates and overall survival in front-line and relapsed/ refractory settings. Although present in all cells, proteasomes are relatively abundant in multiple myeloma cells making that disease a target for proteasome inhibitors. The proteasome Protein activity regulation by synthesis and degradation maintains cellular metabolic integrity and proliferation. The proteasome, a multimeric protease complex, is central to cellular protein regulation by degrading many proteins, thus activating some pathways and shutting down others. Proteins conjugated to multiple units of the polypeptide ubiquitin are degraded by the proteasome. Widely studied inhibitors included: lacta-cystin, a streptomyces metabolite, which is metabolized to lactacystin b-lactone, the active proteasome inhibitor peptide aldehydes, such as carbobenzoxyl-leucinyl-leucinyl-leuc-inal-H (MG-132); boronic acid peptides and others Dipeptide boronic acid derivatives, potent proteasome inhibitors, could be administered to mice, and some were orally bioavailable and relatively stable under physiological conditions. Boronic acid peptides were known to inhibit serine proteases, whereas the dipeptide boronates had a high degree of selectivity for the proteasome and were not potent inhibitors of many common proteases. Critical to the clinical development of proteasome inhibitors was establishing a pharmacodynamic assay which measured the potency and duration of inhibition of the proteasomes in normal blood cells. In animal-model studies, the dipeptide boronic acid, PS-341, rapidly disappeared from the vascular compartment. An accurate pharmacodynamic assay was developed to supplement pharmacokinetic measurements. Both the chymo-tryptic and tryptic activities of the proteasome could be detected using fluorogenic kinetic assays. These fluorogenic kinetic assays were optimized for both whole blood and blood cells.1

Treatment of myeloma has developed rapidly over the past two decades with the advent of proteasome inhibitors and immunomodulatory agents. Despite their dramatic impact on overall patient outcomes, there remains a high-risk group of patients, constituting 20–30% of all cases, who have not benefited to the same extent and continue to have poor outcomes. There is no single pathogenic mechanism that can currently be used to unify and define high-risk disease; instead, the common features are clinically aggressive behaviours, including therapy resistance, proliferation and evasion of apoptosis, which are the consequences of following several evolutionary trajectories driven by distinct gene–gene interactions. The key molecular hallmarks of high-risk disease biology include increased proliferation and the development of a high-risk ecosystem that facilitates both cancer cell survival and failure of the immune response. Genetic lesions underlie the key hallmarks of high-risk disease states and may occur as tumour initiating or progression events. No one of these lesions explains all high-risk disease, and in order to effectively target high-risk disease, we need to understand how these drivers arise, interact with each other and mediate their downstream effects. Pathologically, the high-risk biological states of multiple myeloma are the end stage of a multi-step clinical progression system typical of multiple myeloma comprising benign monoclonal gammopathy; an intermediate stage lacking clinical damage; multiple myeloma itself, which has a range of clinical behaviours; and an easily recognizable leukaemic phase wherein disease is no longer confined to the bone marrow. This progression system is driven by subclonal competition and selective pressures within the bone marrow microenvironment.2 The proteasome is a multisubunit protease complex with an apparent sedimentation coefficient of 20S. Two types of regulatory complexes, named PA700 and PA28, bind to both ends of the cylindrical 20S proteasome to form the dumbbell-like and football-like proteasomes, respectively. The former complex, named the 26S proteasome, is a eukaryotic ATP-dependent protease and appears to be well organized as a large complex of 2 MDa, consisting of approximately 40 polypeptides, to facilitate rapid proteolysis. It is assumed to be a protein “death machine”, destroying a variety of cellular proteins that have acquired a specific degradation signal(s) such as a multiubiquitin chain.3

Bortezomib is an antitumor drug that competitively inhibits proteasome beta-1 and beta-5 subunits. While the impact of bortezomib on protein stability is known, the effect of this drug on intracellular peptides has not been previously explored. A quantitative peptidomics technique was used to examine the effect of treating human embryonic kidney 293T (HEK293T) cells with 5–500nM bortezomib for various lengths of time (30minutes to 16 hours), and human neuroblastoma SH-SY5Y cells with 500nM bortezomib for 1 hour. Although bortezomib treatment decreased the levels of some intracellular peptides, the majority of peptides were increased by 50–500 nM bortezomib.4

Material & methods

Database NCBI (National Center for Biotechnology Information)5 PDB (Protein Data Bank).6 Drug Bank.7 Tools: BLAST (Basic Local Alignment Search Tool).8 Model validation: SAVES (Structural Analysis and Verification Server)9 Model visualization: Chimera, Rasmol, Pymol, discovery studio, Binding site analysis: Qsite Finder, Pocket Finder,10 Dockingtool: Auto Dock, hex, PATCHDOCK, (hexserver.loria.fr/)Automated Docking Server: Online different type of docking server, The first step in methodology is collection of sequences data from NCBI. Sequence alignment: The protein sequences of multiple myloma (>sp|Q96EZ4.2|MYEOV_HUMAN Rec Name: Full=Myeloma-over expressed gene protein; Alt Name: Full=Oncogene in multiple myeloma) were obtained from NCBI/PDB after that the homology modeling of sequence is done then selection of the best model is done with the help of core region and model validation a binding site is also predicted via online tools then go for docking for identification of potential ligand with minimum energy.11

Results and discussion

The results analysis base on Sequence of Rec Name: Full=Myeloma-overexpressed gene protein; Alt Name: Full=Oncogene in multiple myeloma

>sp|Q96EZ4.2|MYEOV_HUMAN Rec Name: Full=Myeloma-overexpressed gene protein; Alt Name: Full=Oncogene in multiple myeloma.12

MALRICVTYTPALPIGLCTRCCLCLEQSPSWCHCLRGVSFLTFHLHQSVPLGDRDSLLMFTRQAGHFV
EGSKAGRSRGRLCLSQALRVAVRGAFVSLWFAAGAGDRERNKGDKGAQTGAGLSQEAEDVDVSRARR
VTDAPQGTLCGTGNRNSGSQSARVVGVAHLGEAFRVGVEQAISSCPEEVHGRHGLSMEIMWARMDV
ALRSPGRGLLAGAGALCMTLAESSCPDYERGRRACLTLHRHPTPHCSTWGLPLRVAGSWLTVVTVEA
LGGWRMGVRRTGQVGPTMHPPPVSGASPLLLHHLLLLLLIIILTC 

Expectation Value=0.001, Search Tool = blast, Mask Low Complexity=yes) via BLASTP .this blast mainly use for protein the results of Computer Aided Drug Designing of Proteosome Inhibitors to find the potential drug candidate for Multiple myloma based on different type potential parameters Homology modeling, also known as comparative modeling of protein, refers to constructing anatomic-resolution model of the "target" protein formats amino acid sequence and an experimental three-dimensional structure of a related homologous protein. In this project homology modeling completed with Geno3D & Phyre it is an automatic web server for protein molecular

Model validation

Model validation completed with the help of SAVES server. It is type of online web server for model validation and for analyzing protein structure for validity and assessing how correct they are it’s based on six programs Results: the Model 1 is pass by the SAVES server (Figure 1).

Figure 1 Ramchandran plot & Free residues structures.

Visualization

Selected model analysis based on 2D and 3D structure visualization and the model visualization completed with the help of different type of software’s for example: Chimera, Rasmol, and discovery studio computational predicted models via bioinformatics approaches and the method of homology modelling is based on the observations because that’s protein tertiary structure is better conserved than amino acid sequence (Figures 2-4).13

Figure 2 3D Structure of protein in a form of back bone & ball sticks.

Figure 3 3D Structure of protein in a form of space fill & strand structure.

Figure 4 3D Structure of protein in a form of Ribbon & labeled.

Binding site prediction

According Steps selected protein model submitted in a Reptox Binding site finder for best receptor ligang binding site.

Binding residues A100 G104 D105 R106 E107 R108 G118 A119 G120 (Figure 5).

Figure 5 3D Structure of protein in a form of protein refined structure.

Suitable Ligand Selection for Receptor Ligand Binding

In this project the ligand selection for potential receptor based on ligand binding site and also available list favorable drugs for multiple myloma. Then select the potential ligand as drug candidate and change some confirmation on drug structure computationally then prepare final ligand for docking results analysis (Figure 6).

Figure 6 Structure of ligand with some conformational changes in main structure.

Docking

Prediction of the optimal physical configuration and energy between two molecules and the phenomenon which enables the interaction between receptor molecules and the ligand molecule, and mainly (Figure 7) (Figure 8).

Figure 7 Different 3 D model structure of docking model number 1 with ligand molecule.

Figure 8 3D Structure of Docking model number 1 in a form of space fill & ball sticks with ligand molecule.

Conclusion

Finally in this study of Proteasome protein structure that can be concluded10 number of model and 10 number of potential ligand with 3Dstructure which is predicted through the Insilico approaches and homology modeling and the docking of Proteasome protein with various ligand determine and the interaction between protein and ligand that’s bind on active site of the Proteasome enzyme, although docking process. It is very complicated because its depends on various parameters the main resultant obtained by different type of docking tools and docking completed with the help of HEX,PATHADOCK for identify the suitable Proteasome inhibitors via different ligands which are docked with the Proteasome protein. Only 1 numbers of ligand given the minimum energy and out of these 10 ligands. Bortezomib is playing an important role in Multiple myloma and bortezomib ligand have minimum binding energy and its work as a potential Proteasome inhibitor as per the Table 1 parameters. Perhaps the ultimate solution is to develop a potential drug candidate against this devastating epidemic.

Docking Model
Number

Binding score

Area

ACE

Transformation information

Docking Model 1

3682

468.70

-192.32    

0.94-0.44-3.0720.88-147.77145.00 (Clustering RMSD: 4.0)

Docking Model 2

3636

430.80

-178.17    

-2.71 -0.81 1.80 143.35 -66.44 68.01

Docking Model 3

3490

384.30

-144.34    

0.43 -0.85 -0.34 26.23 92.40 94.25

Docking Model 4

3480

437.50

-183.08    

2.21 -0.35 -0.59 -50.43 -95.75 176.70

Docking Model 5

3380

423.60

-128.38    

-2.34 1.44 -0.64 -17.03 -178.54 -38.76

Docking Model 6

3316

397.30

-55.96    

-0.29 0.02 0.70 -118.58 9.17 -45.83

Docking Model 7

3270

414.40

-147.82    

2.61 -0.20 -2.42 -4.51 131.65 147.67

Docking Model 8

3256

394.60

-186.25    

-0.13 -0.76 -0.68 47.04 183.08 40.81

Docking Model 9

3232

403.30

-118.56    

-1.63 -1.09 0.30 -130.45 -137.71 28.88

Docking Model 10

3218

407.60

-91.55    

3.01 0.01 2.65 108.48 40.28 94.74

Docking Model 11

3152

452.60

-103.98    

-0.64 1.20 -1.65 72.28 -132.61 -62.16

Docking Model 12

3132

340.70

-148.59    

-1.98 1.27 -1.71 -97.76 -120.81 -57.48

Docking Model 13

3100

394.40

-135.21    

-1.59 -0.27 1.32 12.58 -148.67 -78.65

Docking Model 14

3086

437.40

-97.04    

-1.29 0.21 -2.28 17.51 -27.40 -113.35

Docking Model 15

3084

403.10

-203.82    

-1.66 0.33 0.72 48.43 -16.98 -105.44

Docking Model 16

3070

380.50

-172.17    

-0.22 0.83 1.14 -103.61 23.33 -51.27

Docking Model 17

3064

404.50

-121.96    

0.57 0.31 1.12 -164.56 -66.88 56.60

Docking Model 18

3056

417.90

-81.88    

0.90 -1.02 2.40 -115.45 -85.90 144.86

Docking Model 19

3016

380.60

-121.62    

-2.61 0.29 1.39 123.67 -107.94-12.27

Docking Model 20

2996

373.30

-49.59    

-3.09 1.05 1.34 123.04 -128.71-10.05

Table 1 Amplitude in mill volts of the Lead-1 of electrocardiography in sheep

*Significant (P≤0.05); NSNot significant (P>0.05)

Acknowledgements

It is my proud privilege to express deepest sense of reverence and heart full thanks DBT, for giving me this platform I am extremely great full to Er. Neelesh Yadav from BIFC–FRI, Dehradun for giving me valuable guidance related to advance computer technology.

Conflict of interest

None.

References

Creative Commons Attribution License

©2018 Srivastava, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.