Study on the hydrolysis mechanism of phosphodiesterase 4 using molecular dynamics simulations

The molecular basis for the inhibition of phosphodiesterase-4D by three natural resveratrol analogs. Isolation, molecular docking, molecular dynamics simulations, binding free energy, and bioassay

The phosphodiesterase-4 (PDE4) enzyme is a promising therapeutic target for several diseases. Our previous studies found resveratrol and moracin M to be natural PDE4 inhibitors. In the present study, three natural resveratrol analogs [pterostilbene, (E)-2',3,5',5-tetrahydroxystilbene (THSB), and oxyresveratrol] are structurally related to resveratrol and moracin M, but their inhibition and mechanism against PDE4 are still unclear. A combined method consisting of molecular docking, molecular dynamics (MD) simulations, binding free energy, and bioassay was performed to better understand their inhibitory mechanism. The binding pattern of pterostilbene demonstrates that it involves hydrophobic/aromatic interactions with Phe340 and Phe372, and forms hydrogen bond(s) with His160 and Gln369 in the active site pocket. The present work also reveals that oxyresveratrol and THSB can bind to PDE4D and exhibits less negative predicted binding free energies than pterostilbene, which was qualitatively validated by bioassay (IC50=96.6, 36.1, and 27.0μM, respectively). Additionally, a linear correlation (R(2)=0.953) is achieved for five PDE4D/ligand complexes between the predicted binding free energies and the experimental counterparts approximately estimated from their IC50 values (≈RT ln IC50). Our results imply that hydrophobic/aromatic forces are the primary factors in explaining the mechanism of inhibition by the three products. Results of the study help to understand the inhibitory mechanism of the three natural products, and thus help the discovery of novel PDE4 inhibitors from resveratrol, moracin M, and other natural products.

Insight into the phosphodiesterase mechanism from combined QM/MM free energy simulations

Molecular dynamics simulations employing a combined quantum mechanical and molecular mechanical potential have been carried out to elucidate the reaction mechanism of the hydrolysis of a cyclic nucleotide cAMP substrate by phosphodiesterase 4B (PDE4B). PDE4B is a member of the PDE superfamily of enzymes that play crucial roles in cellular signal transduction. We have determined a two-dimensional potential of mean force (PMF) for the coupled phosphoryl bond cleavage and proton transfer through a general acid catalysis mechanism in PDE4B. The results indicate that the ring-opening process takes place through an S(N)2 reaction mechanism, followed by a proton transfer to stabilize the leaving group. The computed free energy of activation for the PDE4B-catalyzed cAMP hydrolysis is about 13 kcal·mol(-1) and an overall reaction free energy is about -17 kcal·mol(-1), both in accord with experimental results. In comparison with the uncatalyzed reaction in water, the enzyme PDE4B provides a strong stabilization of the transition state, lowering the free energy barrier by 14 kcal·mol(-1). We found that the proton transfer from the general acid residue His234 to the O3' oxyanion of the ribosyl leaving group lags behind the nucleophilic attack, resulting in a shallow minimum on the free energy surface. A key contributing factor to transition state stabilization is the elongation of the distance between the divalent metal ions Zn(2+) and Mg(2+) in the active site as the reaction proceeds from the Michaelis complex to the transition state.