Next Article in Journal
Evaluation of Morphological and Structural Skin Alterations on Diabetic Subjects by Biophysical and Imaging Techniques
Previous Article in Journal
Proteomic Analysis of Ferrochelatase Interactome in Erythroid and Non-Erythroid Cells
Previous Article in Special Issue
Dissipative Particle Dynamics Simulations for Shape Change of Growing Lipid Bilayer Vesicles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Molecular Dynamics Simulation of the Complex of PDE5 and Evodiamine

1
School of Life Sciences, Tokyo University of Pharmacy and Life Sciences, Tokyo 192-0392, Japan
2
The Institute of Statistical Mathematics, Tokyo 190-8562, Japan
3
School of Pharmacy, Tokyo University of Pharmacy and Life Sciences, Tokyo 192-0392, Japan
*
Author to whom correspondence should be addressed.
Life 2023, 13(2), 578; https://doi.org/10.3390/life13020578
Submission received: 8 June 2022 / Revised: 20 January 2023 / Accepted: 13 February 2023 / Published: 18 February 2023

Abstract

:
Alzheimer’s disease is an irreversible neurological disorder for which there are no effective small molecule therapeutics. A phosphodiesterase 5 (PDE5) inhibitor is a candidate medicine for the treatment of Alzheimer’s disease. Rutaecarpine, an indole alkaloid found in Euodiae Fructus, has inhibitory activity for PDE5. Euodiae Fructus contains more evodiamine than rutaecarpine. Therefore, we performed molecular dynamics simulations of the complex of PDE5 and evodiamine. The results showed that the PDE5 and (−)-evodiamine complexes were placed inside the reaction center compared to the case of PDE5 and (+)-evodiamine complex. The binding of (−)-evodiamine to PDE5 increased the root-mean-square deviation and radius of gyration of PDE5. In the PDE5 with (−)-evodiamine complex, the value of the root-mean-square fluctuation of the M-loop, which is thought to be important for activity, increased. This result suggests that (−)-evodiamine may have inhibitory activity.

1. Introduction

Euodiae Fructus is an important crude drug and comprises Goshuyuto, a Kampo medicine clinically used for chronic headaches and migraine attacks [1]. The active compounds of Euodiae Fructus are known as indole alkaloids and quinone alkaloids such as rutaecarpine, evodiamine, and evocarpine (Figure 1a,b) [2]. It is known that the amount of evodiamine in Euodiae Fructus is about the same as the amount of rutaecarpine [2]. Evodiamine has one chiral center. Euodiae Fructus contains more (+)-evodiamine than (−)-evodiamine [3].
Alzheimer’s disease (AD) is characterized by neuronal loss, senile cell groups, and intracellular neurofibrillary changes that cause memory loss [4]. Donepezil [5] and Galantamine [6], acetylcholinesterase inhibitors, and Memantine [7], an N-methyl-D-aspartate (NMDA) receptor antagonist, are known as drugs for AD. In general, it is suggested that AD begins as a disorder produced at least in part by amyloid-β (Aβ). The cAMP-responsive element binding protein (CREB) is an important factor involved in memory [8,9], and phosphorylation of this protein is critical for CREB activity. A typical protein that activates CREB is protein kinase A (PKA), which is known to be activated by cAMP and subsequently phosphorylates CREB [10]. Aβ is known to cause synaptic dysfunction by inhibiting CREB phosphorylation [11].
Phosphodiesterase (PDE) is an enzyme that hydrolyzes cyclic phosphate diesters of cAMP and cGMP. Therefore, when the function of PDE is inhibited, blood levels of cAMP or cGMP increase. PDE inhibitors modulate signaling pathways by increasing levels of cAMP or cGMP and activate CREB via PKA [12,13]; activation of CREB has been reported to promote neurotransmission, synaptic plasticity, and memory formation [14]. mRNA of the PDE5, one of the family proteins of PDE, has been detected in the human cortex, hippocampus, and striatum, and PDE5 inhibitors have been reported to cross the blood–brain barrier and promote vasodilation by inhibiting PDE in neurons and glial cells [15,16]. The phosphodiesterase 5 (PDE5) protein inhibitor is a candidate for the cure for AD [17,18]. Xian-Feng and coworkers synthesized a series of novel PDE5 inhibitors based on rutaecarpine and conducted studies on PDE5 inhibitory activity and in vivo memory and cognitive function [19]. The IC50 (half maximal inhibitory concentration) value of rutaecarpine (Figure 2a) showed PDE5 inhibition at approximately 1230 nM, indicating that it is a promising candidate as an inhibitor. Compound 8i, a rutaecarpine derivative (Figure 2b), was found to be a strong inhibitor with an IC50 value of about 86 nM.
The C-terminal catalytic domain of PDE5 has four subdomain active sites (Q pocket, H pocket, L region, and M site) [20]. The Q pocket is a highly conserved region in the PDE5 catalytic domain and contains the common Gln817, Phe820, Val782, and Tyr612 amino acid residues. The conserved Gln817 amide group forms a bidentate hydrogen bond with the guanine group of cGMP, stabilizing the cGMP structure of the active site. Other residues in this pocket form hydrophobic interactions with cGMP. The H pocket is composed of variable hydrophobic residues that contribute to the inhibitor selectivity of PDE5. The M site contains ions that bind to metal-binding residues common to PDE5, thus stabilizing the structure and activating the hydroxide of the active site. The active site is activated to promote hydrolysis of cGMP and interacts directly with cGMP substrates and indirectly with existing inhibitors. The L region is known to play a role in inhibitor binding by changing its conformation from closed to open [21].
Additionally, the C-terminal catalytic site of PDE5 has two characteristic loop structures. The unit containing residues 661–676 in PDE5 is called the H-loop, and the one containing residues 787–812 is called the M-loop (Figure 3a) [22]. These two loop structures are close to the active site and are important for the activity of the inhibitor. It has been shown that the selective inhibition of PDE5 by sildenafil, which is the best known PDE5 inhibitor, involves the binding between sildenafil and the M-loop or H-loop [20]. It has also been reported that the H-loop and M-loop structures are disordered when sildenafil was bound [22]. As shown in Figure 3b, sildenafil binds with a bidentate hydrogen bond between Gln817.
It has been shown that rutaecarpine, one of the active components of Euodiae Fructus, exhibits PDE5 inhibitory activity [19]. Evodiamine, also from Euodiae Fructus, is an alkaloid with a similar structure to rutaecarpine, but has some differences, such as a methyl group on the N14 side-chain [3].
We performed docking simulations of (+)-evodiamine and (−)-evodiamine with PDE5. Using each complex with smallest free energy, we carried out molecular dynamics (MD) simulations. Our results show that both evodiamines docked in the same pocket as the existing PDE5 inhibitor, and (−)-evodiamine showed the structural change in the M-loop.

2. Model and Methods

2.1. Simulation Preparation

We prepared the structure of PDE5 and its (+) or (−)-evodiamine complex. The PDE5 and evodiamine structures for our MD simulations were obtained from existing databases (PDB and Zinc20). The X-ray crystal structure of catalytic domain of human PDE5A (ID: 1XOZ) [23] was obtained as the structure of PDE5. As the structure of (+) or (−)-evodiamine, we used the three-dimensional structural data in mol2 file format (ID: ZINC898159 and ID: ZINC2031813, respectively) from the Zinc20 database [24]. The structures of each complex were created from the obtained PDE5 structure and the structure of evodiamine using SwissDock [25]. SwissDock resulted in 32 different clusters and a total of 256 structures, of which the lowest energy model was used.

2.2. Molecular Dynamics Simulation

MD simulations of protein with or without the binding of evodiamine were performed using GROMACS (ver 5.1.2 [26], double precision) for 200 ns. A cubic system was constructed so that the distance between the protein and the system boundary is at least 1 nm. The total overall energy was then minimized to 1000 kJ mol−1 nm−1 by steepest descent minimization. We performed the simulation with NVT ensemble for 100 ps and NPT ensemble for 300 ps. After energy minimization, we used velocity-rescaling [27] for NVT and Parrinello–Rahman pressure coupling [28] for NPT. The simulation is performed under periodic boundary conditions. Amber99SB-ILDN [29] was used for the force field of PDE5 and ions, and general AMBER force field (gaff) [30] was used for the force field of evodiamine. We used TIP3P as our water model. The parameters for MD simulations are shown in Table 1. Zn2+ and Mg2+ are the metal ions of PDE5, while Na+ is the counter ion to adjust the charge of the system. Since PDE5 naturally resides in the cytoplasm, it does not have the same ionic properties as saline. Therefore, we chose a total charge of zero for the simulation rather than physiological conditions. The size of the box is the size after the short simulation with NPT ensemble.

2.3. Quantum Mechanical Calculation

The structures of (+)-evodiamine at 1, 75, and 200 ns and (−)-evodiamine at 1 and 75 ns obtained by MD simulations were structurally optimized by quantum mechanical (QM) calculations with Gaussian16 Revision B.01 [31]. We used B3LYP [32] with 6-31G* basis set. Here, a very tight option was used for the structural optimization option. For the grid, Ultrafine was used. Up to 250 steps of optimization calculations were performed. Then, for each optimization structure, vibrational frequency was calculated under the same conditions.

2.4. Analysis

The radius of gyration (Rg), root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), and number of hydrogen bonds were calculated and analyzed for PDE5 and the complexes. These calculations were performed using GROMACS modules. Rg is defined by
R g = i = 1 m i r i 2 i = 1 m i
where m i is the mass of atom i , and r i is the position of atom i . RMSD is defined by
RMSD = i = 1 N m i r i t r i ( t 0 2 i = 1 N m i
where N is the number of atoms, m i is the mass of atom i , t 0 is the start time of the simulation, t is time, and r i is the position of atom i . This calculation was performed after fitting was done using the least squares method. RMSF is defined by
RMSF i = 1 T T 0 t j = T 0 T r i t j r i t 2
where t is time, T 0 is the reference time for the calculation, T is time at the end, r i is position of residue i , and r i t is position averaged over time of particle 𝑖. The relationship between RMSF and B-factor B i is
B i = 8 π 2 3 RMSF i 2
where RMSFi is the RMSF of residue i .
The formation of hydrogen bonds was determined from the distance and angle between the donor atom, acceptor atom, and hydrogen atom. Hydrogen bonds are counted when both of the following two conditions are met (Figure 4). First, the distance between donor and acceptor is less than 0.30 nm. Secondly, the angle θ made by hydrogen, donor and acceptor is less than 30 degrees.

3. Results

3.1. Docking Simulation

Docking simulations were performed using SwissDock shown in Figure 5. As a result, we found one docking site with a binding Gibbs free energy of ΔG = −7.568 kcal/mol, shown in orange in Figure 5a, in (+)-evodiamine. Similiarly, as shown in Figure 5b, (−)-evodiamine has a docking site with the binding free energy of ΔG = −7.568 kcal/mol, shown in orange, and a binding free energy of ΔG = −7.190 kcal/mol, shown in cyan. In both cases, the binding free energy at the orange site was the lowest, so MD simulations were performed for the orange site.

3.2. RMSD and Rg

We performed MD simulations of PDE5 and PDE5 with (+)-evodiamine and PDE5 with (−)-evodiamine. The overlapping snapshots at 0 ns and 200 ns (Figure 6) showed that the M-loop of (−)-evodiamine has changed.
We calculated Rg to ensure that the protein was not unfolding. In Figure 7, we show the Rg values at 200 ns of PDE5. The changes in these values of Rg are relatively small.
We computed the RMSD as shown in Figure 8. The structure of the reference was at 0 ns. In Figure 8a, we show the time evolution of RMSD for PDE5 and PDE5 complexes. The RMSD values of PDE5 with (−)-evodiamine complex were significantly different from those for PDE5 and PDE5 with (+)-evodiamine complex. The change in the RMSD values of the PDE5 with (−)-evodiamine complex became smaller after 50 ns. After 50 ns, PDE5 complexes show stable structure.
The RMSD was calculated to investigate the conformational changes of evodiamine. In Figure 8b, we show the time evolution of RMSD for both cases of evodiamine. This result suggests that there are two states in the structure of evodiamine. The first state is a structure with an RMSD of 0.02 nm to 0.07 nm, and the second state is a structure with an RMSD of 0.10 nm or higher.
We compared conformations of evodiamine. In Figure 9, we show five snapshots of evodiamine. From those snapshots, we can see that one state (Figure 9a1,a3,b2) is a tilted structure, and another state (Figure 9a2,b1) is a planar structure.
We calculated the dihedral angle consisting of C5-N6-C13b-C13a (Figure 10a) to confirm whether this methyl group change is caused by the change in RMSD of evodiamine. The results showed that its time evolution was consistent with the change in RMSD. This indicates that the conformation of evodiamine is changing.
The results of MD simulation show a high frequency of occurrence of structures with a dihedral angle of about −100 (100) degrees (Figure 10b). To investigate the structural stability of (+)-evodiamine, the energies of their structures at 1 ns and 75 ns were calculated by quantum mechanical method. As a result, it was found that the energy of the structure at 75 ns was larger by 0.279 kcal/mol than that at 1 ns. This indicates that the structure at 75 ns is less stable than the other structure.

3.3. Structural Changes in the Loops of the Complex

We have calculated the RMSD of PDE5 and complexes separately for each loop and non-loop structure for changes in RMSD as shown in Figure 11. The results show that in each model, the change in the H-loop was small, while the M-loop in PDE5 with (−)-evodiamine complex showed large RMSD values.

3.4. Effect of Evodiamine Binding on Fluctuations

We calculated the RMSF to study how evodiamine binding affects protein fluctuations. In Figure 12a, we show the RMSF of PDE5 and the complex of PDE5 from 50 ns to 200 ns, when the change over time of whole RMSD becomes small. RMSF results showed that the fluctuations at residues 606–614 and 760–830 were larger in the (−)-evodiamine complex than in PDE5 and its (+)-evodiamine complex. Figure 12c–e shows the average structures at 50–200 ns. In these snapshots, the color changes depending on the value of the B-factor calculated from the RMSF, with the larger value being closer to red. The residue of 607–612 had two states (Figure 12b). Therefore, the peak of 607–612 is not a fluctuation but large RMSF value (Figure 12a).

3.5. Hydrogen Bonds between PDE5 and Evodiamine

To confirm the binding of PDE5 to evodiamine, the number of intermolecular hydrogen bonds was calculated. Figure 13a shows the time evolution of the number of intermolecular hydrogen bonds formed between PDE5 and evodiamine. There are up to three hydrogen bonds between PDE5 and evodiamine (Figure 13a). A maximum of three bonds are formed from the beginning of the simulation up to 50 ns, but after that, no hydrogen bonds were formed for a long time. First, evodiamine forms a hydrogen bond with Tyr612 (H-loop) in the Q pocket (Figure 13b). It then forms a hydrogen bond with Gln817 (M-loop) and is fixed in orientation (Figure 13c). In particular, the bidentate hydrogen bond is similar to the binding of sildenafil with Gln817, as shown in Figure 3b. Finally, a multi-membered ring of (−)-evodiamine is sandwiched between the M-loops in Figure 13d.

4. Discussion

Docking simulations by SwissDock showed that (+)-evodiamine was not in the pocket, contrary to the case of tadalafil, which is an inhibitor. However, in the MD simulation, the multimembered ring moved into the pocket (Figure 14a). In this case, the multimembered ring of (+)-evodiamine was in the same position as the multimembered ring of tadalafil. There is a multi-membered ring between Phe786 and Phe820 in the pocket (Figure 14b). Therefore, a flat structure is considered to be more likely to fit into the pocket. Comparison of RMSDs showed that the complexes have a less altered M-loop structure than PDE5 alone. Then, PDE5 with (+)-evodiamine complex suggested that binding of PDE5 to the ligand constructed the M-loop.
On the other hand, (−)-evodiamine was found to have two docking sites, a reaction center pocket shown with sildenafil (PDBID: 1UDT, Figure 14a), and an allosteric pocket, which was found by T. Zhang et al. [33]. In the results, the ΔG for the reaction center pocket and allosteric pocket are -7.5681434 and −7.189539 kcal/mol, respectively, suggesting that they do not selectively bind to either pocket. In the allosteric pocket, (+)-evodiamine is not pocketed in the simulation. In the case of (−)-evodiamine, we have performed an MD simulation starting with this molecule in the reaction center pocket. In the future, MD simulation starting in the allosteric pockets should be performed.
Comparing the PDE5 and both PDE5 complexes, we found differences in the conformational change of the M-loop in the (−)-evodiamine complex. The planar structure of evodiamine is less stable than the tilted structure. The (−)-evodiamine is ring-inverted with the start of MD. This dihedral angle of the side-chain of Phe786 near the M-loop was found to be different from the case of other models (PDE5 or PDE5 with (+)-evodiamine). This is due to the fact that the (−)-evodiamine was guided by hydrogen bonds and moved to enter between the M-loop. It has been suggested that existing inhibitors can also disrupt the loop structure [20], and our results support this. The binding of nitrogen of evodiamine to Gln817 is essential for this hydrogen bond.
Tadalafil is known to have stronger inhibitory activity than sildenafil [34]. However, unlike sildenafil, it is reported to form only one hydrogen bond between Gln817 and the NH of the indole ring, without forming a bidentate hydrogen bond [35]. Relying on this strong hydrogen bond, tadalafil is stabilized in the active site. The tadalafil backbone is also stiffer than that of sildenafil, suggesting a higher affinity for tadalafil with less entropy loss [23]. Evodiamine has a similar main backbone structure to tadalafil than sildenafil, but there are differences in the formation of bidentate hydrogen bonds. The main backbone structure of evodiamine is similar to that of tadalafil. This simulation confirms that in (−)-evodiamine, the ring inversion changes the conformation of the methyl group of this N14, exposing N14 and forming a bidentate hydrogen bond with PDE5. Our quantum mechanical calculations show that the conformation forming this bidentate hydrogen bond is stable. The methyl group of (−)-evodiamine was positioned in a pseudo-equatorial position. Our results suggest that evodiamine forms a bidentate hydrogen bond similar to that of sildenafil, an existing PDE5 inhibitor. It is possible to consider using a different substituent for this methyl group to create a derivative with stronger binding strength.
Evodiamine binding increased the RMSF of the M-loop and H-loop. The M-loop is known to be involved in the selective inhibitory activity of several families of PDEs [17]. The present results did not confirm the binding of the two loop structures. However, the larger RMSF of these structures may have an effect on the selective inhibitory activity. Future studies with other proteins such as PDE4 should be interesting. A novel method for predicting ligand–protein binding affinity from conformational fluctuations of BRD4 protein by machine learning is reported [36]. It is not the magnitude of the RMSF that is important, but the residue where variation occurs. It is necessary to analyze the characteristics of the RMSF in MD simulations of various inhibitors, including evodiamine in this study. In the future, detailed calculations of coupling free energies should be performed. Therefore, MM/GBSA calculations should be considered along with new MD simulations.

5. Conclusions

PDE5 is a potential target protein for the treatment of Alzheimer’s disease. We performed MD simulations of the PDE5 complexes by docking evodiamine, one of the indole alkaloids in Euodiae Fructus, to PDE5. Binding with (+)-evodiamine increased Rg overall, but rather decreased it in the M-loop. This suggests that the M-loop would possibly be stabilized by the ligand. The results showed that the binding of (−)-evodiamine increased the RMSD and Rg of PDE5. The conformational change in evodiamine was consistent with a change in the C5-N6-C13b-C13a dihedral angle. This indicates that both evodiamine stereoisomers are ring-inverted and stable conformationally. In (−)-evodiamine, the conformation of the M-loop changed significantly. In the PDE5 complexes, an increase in the RMSF of the M-loop was observed, which was attributed to the ring inversion of (−)-evodiamine in the pocket and the associated change in the shape of the M-loop. This result suggests that (−)-evodiamine may have inhibitory activity. Moreover, the addition of functional groups to allow hydrogen bonding at stable conformations may be an indicator for designing more potent inhibitors.

Author Contributions

Conceptualization, A.K. and Y.N.; methodology, A.K. and M.N.; software, A.K. and M.N.; validation, A.K.; formal analysis, A.K., Y.N. and M.N.; investigation, M.N. and A.K.; resources, R.M. and M.T.; data curation, A.K.; writing—original draft preparation, writing—review and editing, Y.M. and M.T.; visualization, Y.N., A.K. and M.N.; supervision, R.M., Y.M. and M.T.; project administration, M.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Odaguchi, H.; Wakasugi, A.; Ito, H.; Hisakazu, S.; Gono, Y.; Sakai, F.; Hanawa, T. The efficacy of goshuyuto, a typical Kampo (Japanese herbal medicine) formula, in preventing episodes of headache. Curr. Med. Res. Opin. 2006, 22, 1587–1597. [Google Scholar] [CrossRef] [PubMed]
  2. Kano, Y.; Qine, Z.; Komatsu, K. On the evaluation of preparation of Chinese medicinal prescription. VI. The changes of the alkaloid contests by processing of Evodia fruit. Yakugaku Zasshi 1991, 111, 32–35. [Google Scholar] [CrossRef] [Green Version]
  3. Nguyen, N.V.T.; Lee, K.R.; Lee, Y.J.; Choi, S.; Kang, J.S.; Mar, W.; Kim, K.H. Chiral high-performance liquid chromatographic separation of evodiamine enantiomers and rutaecarpine, isolated from Evodiae fructus. J. Pharm. Biomed. Anal. 2013, 81–82, 151–159. [Google Scholar] [CrossRef] [PubMed]
  4. Livingston, G.; Sommerlad, A.; Orgeta, V.; Costafreda, S.G.; Huntley, J.; Ames, D.; Ballard, C.; Banerjee, S.; Burns, A.; Cohen-Mansfield, J.; et al. Dementia prevention, intervention, and care. Lancet 2017, 390, 2673–2734. [Google Scholar] [CrossRef] [Green Version]
  5. Ikeda, M.; Mori, E.; Matsuo, K.; Nakagawa, M.; Kosaka, K. Donepezil for dementia with Lewy bodies: A randomized, placebo-controlled, confirmatory phase III trial. Alzheimers Res. Ther. 2015, 7, 4. [Google Scholar] [CrossRef] [Green Version]
  6. Samochocki, M.; Höffle, A.; Fehrenbacher, A.; Jostock, R.; Ludwig, J.; Christner, C.; Radina, M.; Zerlin, M.; Ullmer, C.; Pereira, E.F.R.; et al. Galantamine is an allosterically potentiating ligand of neuronal nicotinic but not of muscarinic acetylcholine receptors. J. Pharmacol. Exp. Ther. 2003, 305, 1024–1036. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Nakamura, S.; Murayama, N.; Noshita, T.; Katsuragi, R.; Ohno, T. Cognitive dysfunction induced by sequential injection of amyloid-β and ibotenate into the bilateral hippocampus; protection by memantine and MK-801. Eur. J. Pharmacol. 2006, 548, 115–122. [Google Scholar] [CrossRef]
  8. Tully, T. Regulation of gene expression and its role in long-term memory and synaptic plasticity. Proc. Natl. Acad. Sci. USA 1997, 29, 4239–4241. [Google Scholar] [CrossRef] [Green Version]
  9. Yin, J.C.P.; Tully, T. CREB and the formation of long-term memory. Curr. Opin. Neurobiol. 1996, 6, 264–268. [Google Scholar] [CrossRef]
  10. Gonzalez, G.A.; Montminy, M.R. Cyclic AMP stimulates somatostatin gene transcription by phosphorylation of CREB at serine 133. Cell 1989, 59, 675–680. [Google Scholar] [CrossRef]
  11. Dineley, K.T.; Westerman, M.; Bui, D.; Bell, K.; Ashe, K.H.; Sweatt, J.D. β-amyloid activates the mitogen-activated protein kinase cascade via hippocampal α7 nicotinic acetylcholine receptors: In vitro and in vivo mechanisms related to Alzheimer’s disease. J. Neurosci. 2001, 21, 4125–4133. [Google Scholar] [CrossRef] [Green Version]
  12. Abel, T.; Nguyen, P.V.; Barad, M.; Deuel, T.A.S.; Kandel, E.R.; Bourtchouladze, R. Genetic demonstration of a role for PKA in the late phase of LTP and in hippocampus-based long-term memory. Cell 1997, 88, 615–626. [Google Scholar] [CrossRef] [Green Version]
  13. Bernabeu, R.; Bevilaqua, L.; Ardenghi, P.; Bromberg, E.; Schmitz, P.; Bianchin, M.; Izquierdo, I.; Medina, J.H. Involvement of hippocampal cAMP/cAMP-dependent protein kinase signaling pathways in a late memory consolidation phase of aversively motivated learning in rats. Proc. Natl. Acad. Sci. USA 1997, 94, 7041–7046. [Google Scholar] [PubMed] [Green Version]
  14. Puzzo, D.; Staniszewski, A.; Deng, S.X.; Privitera, L.; Leznik, E.; Liu, S.; Zhang, H.; Feng, Y.; Palmeri, A.; Landry, D.W.; et al. Phosphodiesterase 5 Inhibition Improves Synaptic Function, Memory, and Amyloid-β Load in an Alzheimer’s Disease Mouse Model. J. Neurosci. 2009, 29, 8075–8086. [Google Scholar] [CrossRef] [PubMed]
  15. Saavedra, A.; Giralt, A.; Arumí, H.; Alberch, J.; Pérez-Navarro, E. Regulation of hippocampal cGMP levels as a candidate to treat cognitive deficits in huntington’s disease. PLoS ONE 2013, 8, e73664. [Google Scholar] [CrossRef] [Green Version]
  16. García-Barroso, C.; Ricobaraza, A.; Pascual-Lucas, M.; Unceta, N.; Rico, A.J.; Goicolea, M.A.; Sallés, J.; Lanciego, J.L.; Oyarzabal, J.; Franco, R.; et al. Tadalafil crosses the blood-brain barrier and reverses cognitive dysfunction in a mouse model of AD. Neuropharmacology 2013, 64, 114–123. [Google Scholar] [CrossRef] [PubMed]
  17. Cuadrado-Tejedor, M.; Hervias, I.; Ricobaraza, A.; Puerta, E.; Pérez-Roldán, J.M.; García-Barroso, C.; Franco, R.; Aguirre, N.; García-Osta, A. Sildenafil restores cognitive function without affecting β-amyloid burden in a mouse model of Alzheimer’s disease. Br. J. Pharmacol. 2011, 164, 2029–2041. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Huai, Q.; Liu, Y.; Francis, S.H.; Corbin, J.D.; Ke, H. Crystal structures of phosphodiesterases 4 and 5 in complex with inhibitor 3-isobutyl-1-methylxanthine suggest a conformation determinant of inhibitor selectivity. J. Biol. Chem. 2004, 279, 13095–13101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Huang, X.; Dong, Y.; Wang, J.; Ke, H.; Song, G.; Xu, D. Novel PDE5 inhibitors derived from rutaecarpine for the treatment of Alzheimer’s disease. Bioorg. Med. Chem. Lett. 2020, 30, 127097. [Google Scholar] [CrossRef]
  20. Sung, B.; Hwang, K.Y.; Jeon, Y.H.; Lee, J.I.; Heo, Y.; Kim, J.H.; Moon, J.; Yoon, J.M.; Hyun, Y.; Kim, E.; et al. Structure of the catalytic domain of human phosphodiesterase 5 with bound drug molecules. Nature 2003, 425, 98–102. [Google Scholar]
  21. Kayık, G.; Tüzün, N.Ş.; Durdagi, S. Investigation of PDE5/PDE6 and PDE5/PDE11 selective potent tadalafil-like PDE5 inhibitors using combination of molecular modeling approaches, molecular fingerprint-based virtual screening protocols and structure-based pharmacophore development. J. Enzyme Inhib. Med. Chem. 2017, 32, 311–330. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Wang, H.; Ye, M.; Robinson, H.; Francis, S.H.; Ke, H. Conformational variations of both phosphodiesterase-5 and inhibitors provide the structural basis for the physiological effects of vardenafil and sildenafil. Mol. Pharmacol. 2008, 73, 104–110. [Google Scholar] [CrossRef]
  23. Card, G.L.; England, B.P.; Suzuki, Y.; Fong, D.; Powell, B.; Lee, B.; Luu, C.; Tabrizizad, M.; Gillette, S.; Ibrahim, P.N.; et al. Structural basis for the activity of drugs that inhibit phosphodiesterases. Structure 2004, 12, 2233–2247. [Google Scholar] [CrossRef] [Green Version]
  24. Irwin, J.J.; Tang, K.G.; Young, J.; Dandarchuluun, C.; Wong, B.R.; Khurelbaatar, M.; Moroz, Y.S.; Mayfield, J.; Sayle, R.A. ZINC20―A Free ultralarge-scale chemical database for ligand discovery. J. Chem. Inf. Model. 2020, 60, 6065–6073. [Google Scholar] [CrossRef] [PubMed]
  25. Grosdidier, A.; Zoete, V.; Michielin, O. SwissDock, a protein-small molecule docking web service based on EADock DSS. Nucleic Acids Res. 2011, 39, W270–W277. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. Software X 2015, 1–2, 19–25. [Google Scholar] [CrossRef] [Green Version]
  27. Bussi, G.; Donadio, D.; Parrinello, M. Canonical sampling through velocity-rescaling. J. Chem. Phys. 2007, 126, 14101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Parrinello, M.; Rahman, A. Polymorphic transitions in single crystals: A new molecular dynamics method. J. Appl. Phys. 1981, 52, 7182–7190. [Google Scholar] [CrossRef]
  29. Lindorff-Larsen, K.; Piana, S.; Palmo, K.; Maragakis, P.; Klepeis, J.L.; Dror, R.O.; Shaw, D.E. Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins 2010, 78, 1950–1958. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Wang, J.; Wolf, R.M.; Caldwell, J.W.; Kollman, P.A.; Case, D.A. Development and testing of a General Amber Force Field. J. Comp. Chem. 2004, 25, 1157–1174. [Google Scholar] [CrossRef]
  31. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16 Revision B.01; Gaussian Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  32. Becke, A.D. Density functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef] [Green Version]
  33. Zhang, T.; Lai, Z.; Yuan, S.; Huang, Y.; Dong, G.; Sheng, C.; Ke, H.; Luo, H. Discovery of evodiamine derivatives as highly selective PDE5 inhibitors targeting a unique allosteric pocket. J. Med. Chem. 2020, 63, 9828–9837. [Google Scholar] [CrossRef] [PubMed]
  34. Sawamura, F.; Kato, M.; Fujita, K.; Nakazawa, T.; Beardsworth, A. Tadalafil, a long-acting inhibitor of PDE5, improves pulmonary hemodynamics and survival rate of monocrotaline-induced pulmonary artery hypertension in rats. J. Pharmacol. Sci. 2009, 111, 235–243. [Google Scholar] [CrossRef] [Green Version]
  35. de Oliveira, I.P.; Lescano, C.H.; De Nucci, G. Q817G mutation in phosphodiesterase type 5: Conformational analysis and dissociation profile of the inhibitor Tadalafil. Chem. Biol. Drug. Des. 2019, 93, 419–429. [Google Scholar] [CrossRef] [PubMed]
  36. Yasuda, I.; Endo, K.; Yamamoto, E.; Hirano, Y.; Yasuoka, K. Differences in ligand-induced protein dynamics extracted from an unsupervised deep learning approach correlate with protein-ligand binding affinities. Commun. Biol. 2022, 5, 481. [Google Scholar] [CrossRef]
Figure 1. The active compounds of Euodiae Fructus. (a) Rutaecarpine and (b) evodiamine, indole alkaloids isolated from Euodiae Fructus. (c) Quinone alkaloid isolated from Euodiae Fructus.
Figure 1. The active compounds of Euodiae Fructus. (a) Rutaecarpine and (b) evodiamine, indole alkaloids isolated from Euodiae Fructus. (c) Quinone alkaloid isolated from Euodiae Fructus.
Life 13 00578 g001
Figure 2. Compounds that have shown PDE5 inhibitory activity. (a) Rutaecarpine and (b) 8i, which is an analog.
Figure 2. Compounds that have shown PDE5 inhibitory activity. (a) Rutaecarpine and (b) 8i, which is an analog.
Life 13 00578 g002
Figure 3. (a) Structure of PDE5 and (b) binding with sildenafil. In (a), the red area is the H-loop, the blue area is the M-loop, and orange is the residues of Q pocket. In (b), the bidentate hydrogen bond is shown between Gln817 in PDE5 and sildenafil.
Figure 3. (a) Structure of PDE5 and (b) binding with sildenafil. In (a), the red area is the H-loop, the blue area is the M-loop, and orange is the residues of Q pocket. In (b), the bidentate hydrogen bond is shown between Gln817 in PDE5 and sildenafil.
Life 13 00578 g003
Figure 4. Conditions for hydrogen bonding.
Figure 4. Conditions for hydrogen bonding.
Life 13 00578 g004
Figure 5. Results of docking simulation. (a) (+)-evodiamine, (b) (−)-evodiamine.
Figure 5. Results of docking simulation. (a) (+)-evodiamine, (b) (−)-evodiamine.
Life 13 00578 g005
Figure 6. Snapshots of PDE and evodiamine complex structure for MD simulation. (a) Superimposed snapshots of the structure of PDE5 at 0 ns and 200 ns. (b) PDE5 with (+)-evodiamine complex, (c) PDE5 with (−)-evodiamine complex. Green: 0 ns, orange: 200 ns.
Figure 6. Snapshots of PDE and evodiamine complex structure for MD simulation. (a) Superimposed snapshots of the structure of PDE5 at 0 ns and 200 ns. (b) PDE5 with (+)-evodiamine complex, (c) PDE5 with (−)-evodiamine complex. Green: 0 ns, orange: 200 ns.
Life 13 00578 g006
Figure 7. Time evolution of Rg.
Figure 7. Time evolution of Rg.
Life 13 00578 g007
Figure 8. Time evolution of RMSD. (a) RMSD of PDE5 in each model, (b) evodiamine.
Figure 8. Time evolution of RMSD. (a) RMSD of PDE5 in each model, (b) evodiamine.
Life 13 00578 g008
Figure 9. Structures of each evodiamine extracted from MD simulations (green) were optimized for QM calculations (cyan) of (+)-evodiamine: (a1) 1 ns, (a2) 75 ns, (a3) 200 ns. (−)-evodiamine: (b1) 1 ns, (b2) 200 ns.
Figure 9. Structures of each evodiamine extracted from MD simulations (green) were optimized for QM calculations (cyan) of (+)-evodiamine: (a1) 1 ns, (a2) 75 ns, (a3) 200 ns. (−)-evodiamine: (b1) 1 ns, (b2) 200 ns.
Life 13 00578 g009
Figure 10. Time evolution of conformational change of (+)-evodiamine (blue) and (−)-evodiamine (red): (a) evodiamine C5-N6-C13b-C13a, (b) dihedral angle of C5-N6-C13b-C13a.
Figure 10. Time evolution of conformational change of (+)-evodiamine (blue) and (−)-evodiamine (red): (a) evodiamine C5-N6-C13b-C13a, (b) dihedral angle of C5-N6-C13b-C13a.
Life 13 00578 g010
Figure 11. Time evolution of RMSD for H-loop and M-loop. (a) PDE5, (b) PDE5-(+)-evodiamine complex, (c) PDE5-(−)-evodiamine complex. (d) Position of H-loop and M-loop of PDE5.
Figure 11. Time evolution of RMSD for H-loop and M-loop. (a) PDE5, (b) PDE5-(+)-evodiamine complex, (c) PDE5-(−)-evodiamine complex. (d) Position of H-loop and M-loop of PDE5.
Life 13 00578 g011
Figure 12. (a) RMSF at 50–200 ns. (b) Snapshot of PDE5 with (−)-evodiamine complex at 50 ns and 200 ns. The structure of PDE5 (c) and PDE5 with (+)-evodiamine complex (d) and PDE5 with (−)-evodiamine (e). In (ce), the closer the color is to red, the higher the value of the B-factor is.
Figure 12. (a) RMSF at 50–200 ns. (b) Snapshot of PDE5 with (−)-evodiamine complex at 50 ns and 200 ns. The structure of PDE5 (c) and PDE5 with (+)-evodiamine complex (d) and PDE5 with (−)-evodiamine (e). In (ce), the closer the color is to red, the higher the value of the B-factor is.
Life 13 00578 g012
Figure 13. Intermolecular hydrogen bonds between PDE5 and evodiamine. (a) Number of hydrogen bonds between PDE5 and evodiamine. (b) Hydrogen bond between the oxygen atom of (−)-evodiamine and the hydroxy group of Tyr612. (c) Hydrogen bonds between the two nitrogen atoms of evodiamine and the side of Gln817. (d) Hydrogen bond between the nitrogen atom of evodiamine and the carbonyl group of Gln817.
Figure 13. Intermolecular hydrogen bonds between PDE5 and evodiamine. (a) Number of hydrogen bonds between PDE5 and evodiamine. (b) Hydrogen bond between the oxygen atom of (−)-evodiamine and the hydroxy group of Tyr612. (c) Hydrogen bonds between the two nitrogen atoms of evodiamine and the side of Gln817. (d) Hydrogen bond between the nitrogen atom of evodiamine and the carbonyl group of Gln817.
Life 13 00578 g013
Figure 14. (a) (+)-evodiamine (magenta), (−)-evodiamine (yellow), and sildenafil (cyan). (b) PDE5 and (+)-evodiamine at 148 ns. Ligand is sandwiched between two phenylalanines.
Figure 14. (a) (+)-evodiamine (magenta), (−)-evodiamine (yellow), and sildenafil (cyan). (b) PDE5 and (+)-evodiamine at 148 ns. Ligand is sandwiched between two phenylalanines.
Life 13 00578 g014
Table 1. The systems of MD simulation for PDE5 and complexes.
Table 1. The systems of MD simulation for PDE5 and complexes.
PDE5PDE5 + (+)-EvodiaminePDE5 + (−)-Evodiamine
Water60,89870,09770,067
Zn2+111
Mg2+111
Na+555
Box size (nm3)1892.42174.32175.8
Temperature (K)300300300
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kobayashi, A.; Nakajima, M.; Noguchi, Y.; Morikawa, R.; Matsuo, Y.; Takasu, M. Molecular Dynamics Simulation of the Complex of PDE5 and Evodiamine. Life 2023, 13, 578. https://doi.org/10.3390/life13020578

AMA Style

Kobayashi A, Nakajima M, Noguchi Y, Morikawa R, Matsuo Y, Takasu M. Molecular Dynamics Simulation of the Complex of PDE5 and Evodiamine. Life. 2023; 13(2):578. https://doi.org/10.3390/life13020578

Chicago/Turabian Style

Kobayashi, Ayame, Motokuni Nakajima, Yoh Noguchi, Ryota Morikawa, Yukiko Matsuo, and Masako Takasu. 2023. "Molecular Dynamics Simulation of the Complex of PDE5 and Evodiamine" Life 13, no. 2: 578. https://doi.org/10.3390/life13020578

APA Style

Kobayashi, A., Nakajima, M., Noguchi, Y., Morikawa, R., Matsuo, Y., & Takasu, M. (2023). Molecular Dynamics Simulation of the Complex of PDE5 and Evodiamine. Life, 13(2), 578. https://doi.org/10.3390/life13020578

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop