A recursive molecular docking coupled with energy-based pose-rescoring and MD simulations to identify <i>h</i>sGC βH-NOX allosteric modulators for cardiovascular dysfunctions
Creators: name: Khalid, Rana Rehan
name: Maryam, Arooma
name: Çınaroğlu, Süleyman Selim
name: Siddiqi, Abdul Rauf
name: Sezerman, Osman Ugur
Publisher: name: Taylor & Francis
Year: 2021
Subjects: subject: Biochemistry
subject: Cell Biology
subject: Genetics
subject: FOS: Biological sciences
subject: Pharmacology
subject: 39999 Chemical Sciences not elsewhere classified
subject: FOS: Chemical sciences
subject: 69999 Biological Sciences not elsewhere classified

description: Modulating the activity of human soluble guanylate cyclase (<i>h</i>sGC) through allosteric regulation of the βH-NOX domain has been considered as an immediate treatment for cardiovascular disorder (CVDs). Currently available βH-NOX domain-specific agonists including cinaciguat are unable to deal with the conundrum raised due to oxidative stress in the case of CVDs and their associated comorbidities. Therefore, the idea of investigating novel compounds for allosteric regulation of <i>h</i>sGC activation has been rekindled to circumvent CVDs. Current study aims to identify novel βH-NOX domain-specific compounds that can selectively turn on sGC functions by modulating the conformational dynamics of the target protein. Through a comprehensive computational drug-discovery approach, we first executed a target-based performance assessment of multiple docking (PLANTS, QVina, LeDock, Vinardo, Smina) scoring functions based on multiple performance metrices. QVina showed the highest capability of selecting true-positive ligands over false positives thus, used to screen 4.8 million ZINC15 compounds against βH-NOX domain. The docked ligands were further probed in terms of contact footprint and pose reassessment through clustering analysis and PLANTS docking, respectively. Subsequently, energy-based AMBER rescoring of top 100 low-energy complexes, per-residue energy decomposition analysis, and ADME-Tox analysis yielded the top three compounds i.e. ZINC000098973660, ZINC001354120371, and ZINC000096022607. The impact of three selected ligands on the internal structural dynamics of the βH-NOX domain was also investigated through molecular dynamics simulations. The study revealed potential electrostatic interactions for better conformational dialogue between βH-NOX domain and allosteric ligands that are critical for the activation of <i>h</i>sGC as compared to the reference compound. HighlightsScoring capability of QVina, Vinardo, Smina, PLANTS, and LeDock were assessed against <i>the h</i>sGC βH-NOX specific active and decoy dataset.Target biasness was achieved whilst screening by retaining two out of three substructures in screened compound that are critical for heme-independent <i>h</i>sGC βH-NOX activation.Target-specific biasness was further ensured by preserving the interactions of functionally important residues in the top three hits.MM-GB/PBSA yields energetically favored docked complexes.Principal Component Analysis (PCA) exhibited noticeable fluctuations in ligand-binding regions of the active-site that are responsible for heme and NO- independent activation. Scoring capability of QVina, Vinardo, Smina, PLANTS, and LeDock were assessed against <i>the h</i>sGC βH-NOX specific active and decoy dataset. Target biasness was achieved whilst screening by retaining two out of three substructures in screened compound that are critical for heme-independent <i>h</i>sGC βH-NOX activation. Target-specific biasness was further ensured by preserving the interactions of functionally important residues in the top three hits. MM-GB/PBSA yields energetically favored docked complexes. Principal Component Analysis (PCA) exhibited noticeable fluctuations in ligand-binding regions of the active-site that are responsible for heme and NO- independent activation.