

Mechanism Study of Aliskiren and Its Analogues by Molecular Dynamic Simulation
English
Mechanism Study of Aliskiren and Its Analogues by Molecular Dynamic Simulation
-
Key words:
- renin inhibitors
- / mechanism
- / molecular docking
- / molecular dynamics (MD) simulations
-
1. INTRODUCTION
RAS (renin-angiotensin system) is an important endocrine system for regulating blood pressure and fluid homeostasis, playing a vital role in the internal environment[1-3]. The overactivation of RAS induces excessive generation of angiotensin II, which leads to the cardiovascular disease[4, 5]. The common drugs for the regulation of RAS are angiotensin I converting enzyme inhibitors (ACEIs), angiotensin receptor blockers (ARBs) and renin inhibitors (RIs)[6-8]. However, ACEIs and ARBs will not totally block the RAS. RIs are considered more effective to prevent and reverse the damage of the target organ, inhibiting renin directly at the rate-limiting step[9, 10]. Based on the structure-based drug design, the first approved nonpeptide renin inhibitor aliskiren was discovered, which possessed the ability to interact with the hydrophobic S1/S3-binding pocket of renin[11]. With the success of aliskiren, the study of renin inhibitors attracted widespread attention[12, 13]. Many researchers have recognized the significance of structure-activity relationship and the action mechanism in renin inhibitors[14-16], which is the key to design new drugs. As some side effects exist in aliskiren[17, 18], it is still necessary to find better renin inhibitors.
Computer-aided drug design (CADD) as a brand new technique for drug discovery is widely used in various pharmaceutical researches[19]. The advancement in CADD enables novel techniques for renin inhibitors to explore the action mechanism[20]. Many minuscule mechanism details of pharmacology were revealed by the method[21-24]. Molecular docking and molecular dynamics (MD) simulation are the principal approaches to calculate the interaction between the drug and targets[25-27]. In recent years, the methods have been widely applied for analyzing rational conformation between drug and its receptor, which promotes the theory discussion of internal mechanism[28-31]. As the renin inhibitors were reported with five major active pockets for bounding to the renin[11], it was essential to understand the details between renin inhibitors and active pockets. Herein, by means of molecular docking and molecular dynamics simulation, two homologous renin inhibitors in different activities were applied to observe the combinative disparity with renin on the molecule level.
2. COMPUTATIONAL METHODS
2.1 Molecular docking
The molecules of high activity (aliskiren) and low activity (compound 1) were chosen as the ligands, which were constructed by Chem3D module of Chemoffice, and their information is exhibited in Fig. 1. Energies of the ligands were auto minimized in the MMFF94 force field for the next docking. The X-ray crystal of renin (ID: 2v10, resembling the common skeleton in different active renin inhibitors) was used as the receptor. Extra crystalline water molecule and unrelated ions such as Na+ and CL- of renin were removed by Pymol 2.3.0 software. AutoDock 4.2 was applied for calculating the binding mode of renin and renin inhibitor. The ligands and receptors added with hydrogen and charges were saved as the format pdbqt for preparation. And the grid box was set as a cube with the length of 1.5 Å. The docking was executed in the form of semi-flexibility with the rigid receptor and flexible ligand. And the optimal binding modes of ligand and receptor were screen out by Lamarckian genetic algorithm.
Figure 1
2.2 Molecular dynamics simulation
MD simulations were adopted to study the mechanism between renin and the different inhibitors by GROMACS 5.1.2 package. Gromacs96 force filed was used for the receptor and ligands[32]. The simulation was carried out in a cubic type of box filled with SPC water, and Na+ and CL- were added for neutralizing the system[33]. Canonical NVT (isothermal-isochoric ensemble) and NPT (isothermal-isobaric ensemble) with periodic boundary conditions were set for the equilibration. The temperature and pressure were kept constant at 300 and 101.3 kPa, respectively. Both NPT and NVT of equilibration were performed for 100 ps. Cut-off distances for the calculation of Coulomb and van der Waals interactions were set as 1.0 nm. Finally, 25 ns simulation was continued with the result of 5 ns simulation as the initial conformation. The simulation conditions for both two complexes (renin-compound 1 and aliskiren) and renin in the absence of ligand were the same.
3. RESULTS AND DISCUSSION
3.1 Result of molecular docking
Optimal conformations for complexes of receptor and ligand (renin-aliskiren and renin-compound 1) were screened out. The highly active ligand aliskiren was obtained with a docking energy value of –8.30 kJ /mol. And the lowly active compound 1 was obtained by the same operation with an energy value of –6.08 kJ/mol. Conformations of the two inhibitors in comparison with the prime ligand of renin are shown in Fig. 2. The docking results denoted the objects corresponding to the way of original ligand bound in the renin crystal. And the interactions of the docked ligand with renin are explicated in Fig. 3.
Figure 2
Figure 3
3.2 Result of molecular dynamics simulation
3.2.1 Landscape of the combination
In the molecular dynamics simulation, the root-mean-square deviation (RMSD) trajectory was performed to investigate the stability of the complexes (renin-aliskiren and renin-compound 1). As Fig. 4 shows, the RMSD of renin combined with aliskiren or compound 1 reached equilibrium faster than the renin in the absence of ligands. And the variation in RMSD of renin binding with ligand is more stable than renin without ligand, indicating that the stabilization of renin would be enhanced in the presence of renin inhibitors.
Figure 4
To find out more information of the interaction and combination between renin and renin inhibitors, the binding energy was calculated by g_mmpbsa, a plugin of Gromacs software, with the results performed in Table 2. In solution, renin inhibitors bound to renin were mainly mediated by van der Waals force, while the role of electrostatic force was very faint. As described in Table 2, the total binding energies of renin-aliskiren and renin-compound 1 were respectively at –271.78 and –221.22 kJ/mol, demonstrating that the affinity of renin and aliskiren was stronger, in accordance with the high experimental pIC50 of aliskiren.
Table 1
Pocket S3sp S1 S3 S1' S2' Residues in renin
(code: 2v10)Gln 19 Gly 228 Val36 Phe 119 Thr18 Leu 121 Ser84 Pro 306 Gly40 Tyr83 Tyr 20 Ala 229 Asp 38 Phe 124 Gln19 Ala 122 Leu224 Thr309 Ser41 Glu135 Val 36 Ser 230 Tyr 83 Val 127 Pro118 Phe 124 Asp226 Thr312 Leu81 Ile137 Val 127 Ala 317 Thr 85 Asp 226 Phe119 Ser 230 Gly228 Ala314 Arg82 Tyr83 Tyr 162 Ile305 Pro 306 Table 2
Complex Polar contributions (kJ/mol) Nonpolar contributions (kJ/mol) ΔGbind (kJ/mol) ΔEele ΔGele, sol ΔEvdw ΔGnonpol, sol Renin-aliskiren –21.88 105.79 –327.39 –28.30 –271.78 Renin-compound 1 –44.97 86.21 –242.01 –20.45 –221.22 The b-factor (temperature factor) in Fig. 5 was introduced to describe the details of overall combination. The mean structures of renin, renin-aliskiren and renin-compound 1 were calculated to investigate the b-factor. In Fig. 5, the stabilization was described by colors, and the extents of stability in proper order were blue, green, yellow and red. Blue represents the most stable conformation and red refers to the most active conformation. According to the b-factor, most residues of the complexes renin-aliskiren and renin-compound 1 are blue, representing that the structures are stable (Figs. 5a and 5b). However, in the case of the renin without ligand, there are many other colors except blue (Fig. 5c). And in the active pockets, this phenomenon was prominent. Apparently, contrasting with the renin absent of ligands, the stability of renin was promoted by the renin inhibitors.
Figure 5
The covariance matrix of α-C atom in renin residues based on 2-dimension principal component analyses (PCA) was performed in Fig. 6. It was applied for describing the intensive motion of particles. In the covariance matrix, red indicated that two atoms were moving towards the same direction, while blue was on the opposite. By comparing with the renin without ligand in Fig. 6a, the color depth of the blue decreased to some extent in Fig. 6b and Fig. 6c. This suggested the renin inhibitors would reduce irrelevance between renin residues and promote stability of renin due to the presence of van der Waals forces generated by ligands.
Figure 6
3.2.2 Details of the binding sites
In order to detect the specific details of the combination between renin and renin inhibitors, the interactions of renin inhibitors in active pockets of renin were investigated. The important binding pockets (S1, S1', S2', S3, S3sp) for the renin inhibitors had been described by Rahuel et al, and the components of the pockets are listed in Table 1[11]. The hydrophobic interactions had influence on the inhibitors binding to the active sites. As depicted in Fig. 7a, aliskiren was restricted to a large hydrophobic pocket consist of the key binding site S1, S1', S2', S3 and S3sp. Compound 1 was in the hydrophobic pocket formed by S1, S1', S3 and S3sp in Fig. 7b. Distinctly, aliskiren occupied more binding sites than compound 1 due to the larger size. Comparison for the interactions of the two different ligands is shown in Fig. 8. The residues surrounded with the ligand at range of 6 Å were described as ribbons, and the ligands are shown as sticks in Figs. 8a and 8c. According to Fig. 8, the conformation of aliskiren was more extended than that of compound 1. Therefore, there was more opportunity for aliskiren to contact with the key residues of renin. And compared with compound 1, the orientation of aliskiren appeared better to fit the surface of the hydrophobic pockets.
Figure 7
Figure 8
A higher resolution scan of the interaction between renin and renin inhibitors was demonstrated in Figs. 8b and 8d. As shown in the images, the binding occurred mostly through weak non-covalent interactions, such as van der Waals force and hydrogen bonding. Among the forces, van der Waals force was the major driving forces. In Figs. 8b and 8d, more interactions were found between renin and aliskiren than that between renin with compound 1. According to Rahuel, the interaction with S1/S3-binding pocket in renin promised the vitro activity. Besides, the additional hydrogen bonding or hydrophobic effects with binding site S3sp of inhibitors were essential for the higher bioactivity[11]. Exactly, the result of MD simulation was matched with the rule. Both the highly and lowly active renin inhibitors could interact with S1, but only aliskiren formed hydrogen bonding with the S3sp. As Fig. 8b shows, aliskiren possessed electrostatic interactions with pocket S3sp (Gly228), S2' (Gln135, Gly40) and S1 (Asp226, Asp38), conjugation effect with S1 (Phe124, Tyr83, Val127), S1' (Leu224), S2' (Ile137) and S3sp (Val127) and van der Waals force with all the pockets above. However, compound 1 possessed electrostatic and conjugate interactions with S1 (Tyr83 and Thr85) and S3 (Pro118, Phe119, Ala122 and Phe124) merely (Fig. 8d), and the van der Waals interactions of compound 1 were less than that of aliskiren. Apparently, 1 was not optimistic to interact with the important pocket, thus resulting in the undesired inhibitory activity.
In addition, Asp32 and Asp 215 were imperative for catalysis of renin. Occupying the sites was favor for inhibiting activity of renin inhibitors. The binding energy variation of the two inhibitors with Asp38 and Asp 226 was extracted in Fig. 9. As it depicted, aliskiren possessed more stable and powerful binding energy with the two key residues than compound 1. This reveals the fact that compound 1 is not effective enough to meet the interactive needs of renin.
Figure 9
The binding energy contribution in Fig. 10 revealed the role of residues to the combination. In the complex of renin-aliskiren, residues Gln19, Tyr83, Thr85 and Phe124 made the significant contribution to the binding. And in the complex of renin-compound 1, Arg82, Tyr83 and Phe124 were vital for the combination. The contribution of residue Tyr83 and Phe124 in pocket S1 was prominent in both the two complexes. It proved that pocket S1 was important for the bioactivity of renin inhibitors as well.
Figure 10
4. CONCLUSIONS
The match manners of renin and renin inhibitors were detected by molecular docking in the study. For more insight, the highly active inhibitor aliskiren and the lowest active compound 1 were compared in aqueous solution through the RMSD, binding energy, b-factor, PCA and internal interactions by MD simulations. The mechanism of different inhibitors exhibiting diverse activity has been revealed. The results of RMSD, b-factor and PCA confirmed that the stabilization of renin could be promoted in the presence of the ligands. The binding energy revealed that the combination of renin and renin inhibitors was mainly mediated by van der Waals force. And the higher active inhibitor had stronger binding ability with key residues Asp38 and Asp226. The calculated interactions demonstrated that both the renin inhibitors interacted with a large hydrophobic pocket of renin, and the structure of higher active inhibitor was more extended to access to interact with the residues in the binding sites. The interaction of renin inhibitor with S1 and S3 binding site guaranteed the bioactivity, and the additional hydrogen-bonding interactions interaction to S3sp could enhance the activity. The results help us comprehend the mechanism of renin inhibitor more profoundly.
-
-
[1]
Carey, R. M. The intrarenal renin-angiotensin system in hypertension. Adv. Chronic Kidney Dis. 2015, 22, 204-210. doi: 10.1053/j.ackd.2014.11.004
-
[2]
Kangussu, L. M.; Marzano, L. A. S.; Souza, C. F.; Dantas, C. C.; Miranda, A. S.; Simoes, E. S. The renin-angiotensin system and the cerebrovascular diseases: experimental and clinical evidence. Protein Pept. Lett. 2020, 27, 463-475.
-
[3]
Te-Riet, L.; van-Esch, J. H.; Roks, A. J.; van-den-Meiracker, A. H.; Danser, A. H. Hypertension: renin-angiotensin-aldosterone system alterations. Circ. Res. 2015, 116, 960-975. doi: 10.1161/CIRCRESAHA.116.303587
-
[4]
Nguyen, G.; Delarue, F.; Burckle, C.; Bouzhir, L.; Giller, T.; Sraer, J. D. Pivotal role of the renin/prorenin receptor in angiotensin II production and cellular responses to renin. J. Clin. Invest. 2002, 109, 1417-1427. doi: 10.1172/JCI0214276
-
[5]
Giani, J. F.; Janjulia, T.; Taylor, B.; Bernstein, E. A.; Shah, K.; Shen, X. Z.; McDonough, A. A.; Bernstein, K. E.; Gonzalez-Villalobos, R. A. Renal generation of angiotensin II and the pathogenesis of hypertension. Curr. Hypertens Rep. 2014, 16, 477-5. doi: 10.1007/s11906-014-0477-1
-
[6]
Ramya, K.; Suresh, R.; Kumar, H. Y.; Kumar, B. R. P.; Murthy, N. B. S. Decades-old renin inhibitors are still struggling to find a niche in antihypertensive therapy. A fleeting look at the old and the promising new molecules. Bioorg. Med. Chem. 2020, 28, 115466-12. doi: 10.1016/j.bmc.2020.115466
-
[7]
Wright, J. M.; Musini, V. M.; Gill, R. First-line drugs for hypertension. Cochrane Database Syst. Rev. 2010, 128, 47-52.
-
[8]
Neal, B.; MacMahon, S.; Chapman, N. Blood pressure lowering treatment trialists effects of ACE inhibitors, calcium antagonists, and other blood-pressure-lowering drugs: results of prospectively designed overviews of randomised trials. Blood pressure lowering treatment trialists' collaboration. Lancet 2000, 356, 1955-1964. doi: 10.1016/S0140-6736(00)03307-9
-
[9]
Hamet, P. Direct renin inhibition: mechanistic advantages and disadvantages compared with angiotensin-converting enzyme inhibitors and angiotensin II receptor blockers. Can. J. Cardiol. 2008, 24, 44-49. doi: 10.1016/S0828-282X(08)71038-9
-
[10]
Angeli, F.; Reboldi, G.; Mazzotta, G.; Poltronieri, C.; Garofoli, M.; Ramundo, E.; Biadetti, A.; Verdecchia, P. Safety and efficacy of aliskiren in the treatment of hypertension and associated clinical conditions. Curr. Drug Saf. 2012, 7, 76-85. doi: 10.2174/157488612800492771
-
[11]
Rahuel, J.; Rasetti, V.; Maibaum, J.; Rueger, H.; Goschke, R.; Cohen, N. C.; Stutz, S.; Cumin, F.; Fuhrer, W.; Wood, J. M.; Grutter, M. G. Structure-based drug design: the discovery of novel nonpeptide orally active inhibitors of human renin. Chem. Biol. 2000, 7, 493-504. doi: 10.1016/S1074-5521(00)00134-4
-
[12]
Pantzaris, N. D.; Karanikolas, E.; Tsiotsios, K.; Velissaris, D. Renin inhibition with aliskiren: a decade of clinical experience. J. Clin. Med. 2017, 6, 61-80. doi: 10.3390/jcm6060061
-
[13]
Ichihara, A.; Sakoda, M.; Kurauchi-Mito, A.; Narita, T.; Kinouchi, K.; Bokuda, K.; Itoh, H. New approaches to blockade of the renin-angiotensin-aldosterone system: characteristics and usefulness of the direct renin inhibitor aliskiren. J. Pharmacol. Sci. 2010, 113, 296-300. doi: 10.1254/jphs.10R04FM
-
[14]
Loganathan, L.; Muthusamy, K. Investigation of drug interaction potentials and binding modes on direct renin inhibitors: a computational modeling studies. Lett. Drug Des. Discov. 2019, 16, 919-938. doi: 10.2174/1570180815666180827113622
-
[15]
Politi, A.; Durdagi, S.; Moutevelis-Minakakis, P.; Kokotos, G.; Papadopoulos, M. G.; Mavromoustakos, T. Application of 3D QSAR comfa/comsia and in silico docking studies on novel renin inhibitors against cardiovascular diseases. Eur. J. Med. Chem. 2009, 44, 3703-3711. doi: 10.1016/j.ejmech.2009.03.040
-
[16]
Al-Nadaf, A. H.; Taha, M. O. Discovery of new renin inhibitory leads via sequential pharmacophore modeling, QSAR analysis, in silico screening and in vitro evaluation. J. Mol. Graph. Model 2011, 29, 843-864. doi: 10.1016/j.jmgm.2011.02.001
-
[17]
Sanoski, C. A. Aliskiren: an oral direct renin inhibitor for the treatment of hypertension. Pharmacotherapy 2009, 29, 193-212. doi: 10.1592/phco.29.2.193
-
[18]
Morisawa, N.; Sugano, N.; Yamakawa, T.; Kuriyama, S.; Yokoo, T. Successful long-term effects of direct renin inhibitor aliskiren in a patient with atherosclerotic renovascular hypertension. CEN Case Rep. 2017, 6, 66-73. doi: 10.1007/s13730-016-0246-x
-
[19]
Kapetanovic, I. M. Computer-aided drug discovery and development (CADDD): in silico-chemico-biological approach. Chem. Biol. Interact. 2008, 171, 165-176. doi: 10.1016/j.cbi.2006.12.006
-
[20]
Subramanian, G. Computational modeling and design of renin inhibitors. Bioorg. Med. Chem. Lett. 2013, 23, 460-465. doi: 10.1016/j.bmcl.2012.11.059
-
[21]
Quezada, G. R.; Piceros, E.; Robles, P; Moraga, C.; Jeldres, R. I. Polyacrylic acid to improve flotation tailings management: understanding the chemical interactions through molecular dynamics. Metals. 2021, 11, 987-993. doi: 10.3390/met11060987
-
[22]
Dimitropoulos, N.; Papakyriakou, A.; Dalkas, G. A.; Sturrock, E. D.; Spyroulias G. A. A computational approach to the study of the binding mode of dual ACE/NEP inhibitors. J. Chem. Inf. Model 2010, 50, 388-396. doi: 10.1021/ci9005047
-
[23]
Lam, P. Y.; Jadhav, P. K.; Eyermann, C. J.; Hodge, C. N.; Ru, Y.; Bacheler, L. T.; Meek, J. L.; Otto, M. J.; Rayner, M. M.; Wong, Y. N. Rational design of potent, bioavailable, nonpeptide cyclic ureas as HIV protease inhibitors. Science 1994, 263, 380-4. doi: 10.1126/science.8278812
-
[24]
Singh, A.; Paliwal, S. K.; Sharma, M.; Mittal, A.; Sharma, S.; Sharma, J. P. In silico and in vitro screening to identify structurally diverse non-azole CYP51 inhibitors as potent antifungal agent. J. Mol. Graph. Model 2016, 63, 1-7. doi: 10.1016/j.jmgm.2015.10.014
-
[25]
Kairys, V.; Baranauskiene, L.; Kazlauskiene, M.; Matulis, D.; Kazlauskas, E. Binding affinity in drug design: experimental and computational techniques. Expert. Opin. Drug Discov. 2019, 14, 755-768. doi: 10.1080/17460441.2019.1623202
-
[26]
Zhu, G. F.; Wang, Y.; Liu, J.; Wang, H.; Xi, L.; Du, L. F. Interaction between ginkgolic acid and human serum albumin by spectroscopy and molecular modeling methods. J. Solution Chem. 2014, 43, 1232-1249. doi: 10.1007/s10953-014-0200-5
-
[27]
Tang, J. H.; Liang, G. B.; Zheng, C. Z.; Lian, N. Investigation on the binding behavior of ellagic acid to human serum albumin in aqueous solution. J. Solution Chem. 2013, 42, 226-238. doi: 10.1007/s10953-012-9938-9
-
[28]
Wang, Z. Z.; Ma, C. Y.; Yang, J.; Gao, Q. B.; Sun, X. D.; Ding, L. N.; Liu, H. M. Investigating the binding mechanism of (4-cyanophenyl) glycine derivatives as reversible LSD1 by 3D-QSAR, molecular docking and molecular dynamics simulations. J. Mol. Struct. 2019, 1175, 698-707. doi: 10.1016/j.molstruc.2018.08.029
-
[29]
Yan, W.; Lin, G.; Zhang, R.; Liang, Z.; Wu, L.; Wu, W. Studies on molecular mechanism between ACE and inhibitory peptides in different bioactivities by 3D-QSAR and MD simulations. J. Mol. Liquids 2020, 304, 112702-11. doi: 10.1016/j.molliq.2020.112702
-
[30]
Geng, S.; Jiang, Z.; Ma, H.; Wang, Y.; Liu, B.; Liang, G. Interaction mechanism of flavonoids and bovine beta-lactoglobulin: experimental and molecular modelling studies. Food Chem. 2020, 312, 126066-8. doi: 10.1016/j.foodchem.2019.126066
-
[31]
Lv, Y.; Wang, Y.; Zheng, X.; Liang, G. Reveal the interaction mechanism of five old drugs targeting VEGFR2 through computational simulations. J. Mol. Graph. Model 2020. 96, 107538-9. doi: 10.1016/j.jmgm.2020.107538
-
[32]
Kouza, M.; Hu, C. K.; Zung, H.; Li, M. S. Protein mechanical unfolding: importance of non-native interactions. J. Chem. Physics 2009, 131, 215103-6. doi: 10.1063/1.3272275
-
[33]
Abraham, M. J.; Gready, J. E. Optimization of parameters for molecular dynamics simulation using smooth particle-mesh Ewald in GROMACS 4.5. J. Comput. Chem. 2011, 32, 2031-2040. doi: 10.1002/jcc.21773
-
[1]
-
Table 1. Component of the Pockets of Renin
Pocket S3sp S1 S3 S1' S2' Residues in renin
(code: 2v10)Gln 19 Gly 228 Val36 Phe 119 Thr18 Leu 121 Ser84 Pro 306 Gly40 Tyr83 Tyr 20 Ala 229 Asp 38 Phe 124 Gln19 Ala 122 Leu224 Thr309 Ser41 Glu135 Val 36 Ser 230 Tyr 83 Val 127 Pro118 Phe 124 Asp226 Thr312 Leu81 Ile137 Val 127 Ala 317 Thr 85 Asp 226 Phe119 Ser 230 Gly228 Ala314 Arg82 Tyr83 Tyr 162 Ile305 Pro 306 Table 2. Binding Energy of Renin-aliskiren and Renin-compound 1
Complex Polar contributions (kJ/mol) Nonpolar contributions (kJ/mol) ΔGbind (kJ/mol) ΔEele ΔGele, sol ΔEvdw ΔGnonpol, sol Renin-aliskiren –21.88 105.79 –327.39 –28.30 –271.78 Renin-compound 1 –44.97 86.21 –242.01 –20.45 –221.22 -

计量
- PDF下载量: 3
- 文章访问数: 553
- HTML全文浏览量: 68