Design, synthesis and the structure-activity relationship of agonists targeting on the ALDH2 catalytic tunnel

Repurposing Drugs for Inhibition against ALDH2 via a 2D/3D Ligand-Based Similarity Search and Molecular Simulation

Aldehyde dehydrogenase-2 (ALDH2) is a crucial enzyme participating in intracellular aldehyde metabolism and is acknowledged as a potential therapeutic target for the treatment of alcohol use disorder and other addictive behaviors. Using previously reported ALDH2 inhibitors of Daidzin, CVT-10216, and CHEMBL114083 as reference molecules, here we perform a ligand-based virtual screening of world-approved drugs via 2D/3D similarity search methods, followed by the assessments of molecular docking, toxicity prediction, molecular simulation, and the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) analysis. The 2D molecular fingerprinting of ECFP4 and FCFP4 and 3D molecule-shape-based USRCAT methods show good performances in selecting compounds with a strong binding behavior with ALDH2. Three compounds of Zeaxanthin (q = 0), Troglitazone (q = 0), and Sequinavir (q = +1 e) are singled out as potential inhibitors; Zeaxanthin can only be hit via USRCAT. These drugs displayed a stronger binding strength compared to the reported potent inhibitor CVT-10216. Sarizotan (q = +1 e) and Netarsudil (q = 0/+1 e) displayed a strong binding strength with ALDH2 as well, whereas they displayed a shallow penetration into the substrate-binding tunnel of ALDH2 and could not fully occupy it. This likely left a space for substrate binding, and thus they were not ideal inhibitors. The MM-PBSA results indicate that the selected negatively charged compounds from the similarity search and Vina scoring are thermodynamically unfavorable, mainly due to electrostatic repulsion with the receptor (q = -6 e for ALDH2). The electrostatic attraction with positively charged compounds, however, yielded very strong binding results with ALDH2. These findings reveal a deficiency in the modeling of electrostatic interactions (in particular, between charged moieties) in the virtual screening via the 2D/3D similarity search and molecular docking with the Vina scoring system.

Computational Investigation of Structural Basis for Enhanced Binding of Isoflavone Analogues with Mitochondrial Aldehyde Dehydrogenase

Isoflavone compounds are potent inhibitors against mitochondrial aldehyde dehydrogenase (ALDH2) for the treatment of alcoholism and drug addiction, and an in-depth understanding of the underlying structural basis helps design new inhibitors for enhanced binding. Here, we investigated the binding poses and strengths of eight isoflavone analogues (including CVT-10216 and daidzin) with ALDH2 via computational methods of molecular docking, molecular dynamics (MD) simulation, molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), steered MD, and umbrella sampling. Neither the Vina scoring of docked and MD-sampled complexes nor the nonbonded protein-inhibitor interaction energy from MD simulations is able to reproduce the relative binding strength of the inhibitors compared to experimental IC₅₀ values. Considering the solvation contribution, MM-PBSA and relatively expensive umbrella sampling yield good performance for the relative binding (free) energies. The isoflavone skeleton prefers to form π-π stacking, π-sulfur, and π-alkyl interactions with planar (Phe and Trp) or sulfur-containing (Cys and Met) residues. The enhanced inhibition of CVT-10216 originates from both end groups of the isoflavone skeleton offering strong van der Waals contacts and from the methylsulfonamide group at the 4' position by hydrogen bonding (HB) with neighboring receptor residues. These results indicate that the hydrophobic binding tunnel of ALDH2 is larger than the isoflavone skeleton in length and thus an extended hydrophobic core is likely a premise for potent inhibitors.