Pharmacophore-based virtual screening approach for identification of potent natural modulatory compounds of human Toll-like receptor 7

Toll-like receptor 7 (TLR7) is a transmembrane glycoprotein playing very crucial role in the signaling pathways involved in innate immunity and has been demonstrated to be useful in fighting against infectious disease by recognizing viral ssRNA & specific small molecule agonists. In order to find novel human TLR7 (hTLR7) modulators, computational ligand-based pharmacophore modeling approach was used to identify the molecular chemical features required for the modulation of hTLR7 protein. A training set of 20 TLR7 agonists with their known experimental activity was used to create pharmacophore model using 3D-QSAR pharmacophore generation (HypoGen algorithm) module in Discovery Studio. The best developed hypothesis consists of four pharmacophoric features namely, one hydrogen bond donor (HBD), one ring aromatic (RA), and two hydrophobic (HY) character. The developed hypothesis was then validated by different methods such as cost analysis, test set method, and Fischer’s test method for consistency. Hence, this validated model was further employed for screening of natural hit compounds from InterBioScreen Natural product database, consisting of more than 60,000 natural compounds and derivatives. The screened hit compounds were subsequently filtered by Lipinski’s rule of 5, ADME and toxicity parameters and molecular docking studies to remove the false positive rates. Finally, molecular docking analysis led to identification of the (3a′S,6a′R)-3′-(3,4-dihydroxybenzyl)-5′-(3,4-dimethoxyphenethyl)-5-ethyl-3′,3a′-dihydro-2′H-spiro[indoline-3,1′-pyrrolo[3,4-c]pyrrole]-2,4′,6′(5′H,6a′H)-trione (Compound ID: STOCK1N-65837) as potent hTLR7 modulator due to its better docking score and molecular interactions compared to other compounds. The result of virtual screening was further validated using molecular dynamics (MD) simulation analysis. Thus, a 30 ns MD simulation analysis revealed high stability and effective binding of STOCK1N-65837 within the binding site of hTLR7. Therefore, the present study provides confidence for the utility of the selected chemical feature based pharmacophore model to design novel TLR7 modulators with desired biological activity.