Screening of Natural Products Targeting SARS-CoV-2-ACE2 Receptor Interface - A MixMD Based HTVS Pipeline

The COVID-19 pandemic, caused by novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a severe global health crisis now. SARS-CoV-2 utilizes its Spike protein receptor-binding domain (S-protein) to invade human cell through binding to Angiotensin-Converting Enzyme 2 receptor (ACE2). S-protein is the key target for many therapeutics and vaccines. Potential S-protein-ACE2 fusion inhibitor is expected to block the virus entry into the host cell. In many countries, traditional practices, based on natural products (NPs) have been in use to slow down COVID-19 infection. In this study, a protocol was applied that combines mixed solvent molecular dynamics simulations (MixMD) with high-throughput virtual screening (HTVS) to search NPs to block SARS-CoV-2 entry into the human cell. MixMD simulations were employed to discover the most promising stable binding conformations of drug-like probes in the S-protein-ACE2 interface. Detected stable sites were used for HTVs of 612093 NPs to identify molecules that could interfere with the S-protein-ACE2 interaction. In total, 19 NPs were selected with rescoring model. These top-ranked NP-S-protein complexes were subjected to classical MD simulations for 300 ns (3 replicates of 100 ns) to estimate the stability and affinity of binding. Three compounds, ZINC000002128789, ZINC000002159944 and SN00059335, showed better stability in all MD runs, of which ZINC000002128789 was predicted to have the highest binding affinity, suggesting that it could be effective modulator in RBD-ACE2 interface to prevent SARS-CoV-2 infection. Our results support that NPs may provide tools to fight COVID-19

Substrate Selectivity of Coumarin Derivatives by Human CYP1 Enzymes: In Vitro Enzyme Kinetics and In Silico Modeling

Of the three enzymes in the human cytochrome P450 family 1, CYP1A2 is an important enzyme mediating metabolism of xenobiotics including drugs in the liver, while CYP1A1 and CYP1B1 are expressed in extrahepatic tissues. Currently used CYP substrates, such as 7-ethoxycoumarin and 7-ethoxyresorufin, are oxidized by all individual CYP1 forms. The main aim of this study was to find profluorescent coumarin substrates that are more selective for the individual CYP1 forms. Eleven 3-phenylcoumarin derivatives were synthetized, their enzyme kinetic parameters were determined, and their interactions in the active sites of CYP1 enzymes were analyzed by docking and molecular dynamic simulations. All coumarin derivatives and 7-ethoxyresorufin and 7-pentoxyresorufin were oxidized by at least one CYP1 enzyme. 3-(3-Methoxyphenyl)-6-methoxycoumarin (19) was 7-O-demethylated by similar high efficiency [21-30 ML/(min.mol CYP)] by all CYP1 forms and displayed similar binding in the enzyme active sites. 3-(3-Fluoro-4-acetoxyphenyl)coumarin (14) was selectively 7-O-demethylated by CYP1A1, but with low efficiency [0.16 ML/(min mol)]. This was explained by better orientation and stronger H-bond interactions in the active site of CYP1A1 than that of CYP1A2 and CYP1B1. 3-(4-Acetoxyphenyl)-6-chlorocoumarin (20) was 7-O-demethylated most efficiently by CYP1B1 [53 ML/(min.mol CYP)], followed by CYP1A1 [16 ML/(min.mol CYP)] and CYP1A2 [0.6 ML/(min.mol CYP)]. Variations in stabilities of complexes between 20 and the individual CYP enzymes explained these differences. Compounds 14, 19, and 20 are candidates to replace traditional substrates in measuring activity of human CYP1 enzymes

Fragment‐ and Negative Image‐Based Screening of Phosphodiesterase 10A Inhibitors

A novel virtual screening methodology called fragment‐ and negative image‐based (F‐NiB) screening is introduced and tested experimentally using phosphodiesterase 10A (PDE10A) as a case study. Potent PDE10A‐specific small‐molecule inhibitors are actively sought after for their antipsychotic and neuroprotective effects. The F‐NiB combines features from both fragment‐based drug discovery and negative image‐based (NIB) screening methodologies to facilitate rational drug discovery. The selected structural parts of protein‐bound ligand(s) are seamlessly combined with the negative image of the target's ligand‐binding cavity. This cavity‐ and fragment‐based hybrid model, namely its shape and electrostatics, is used directly in the rigid docking of ab initio generated ligand 3D conformers. In total, 14 compounds were acquired using the F‐NiB methodology, 3D quantitative structure‐activity relationship modeling and pharmacophore modeling. Three of the small‐molecules inhibited PDE10A at ~27 μM to ~67 μM range in a radiometric assay. In a larger context, the study shows that the F‐NiB provides a flexible way to incorporate small‐molecule fragments into the drug discovery.peerReviewe

Substrate Selectivity of Coumarin Derivatives by Human CYP1 Enzymes : In Vitro Enzyme Kinetics and In Silico Modeling

Of the three enzymes in the human cytochrome P450 family 1, CYP1A2 is an important enzyme mediating metabolism of xenobiotics including drugs in the liver, while CYP1A1 and CYP1B1 are expressed in extrahepatic tissues. Currently used CYP substrates, such as 7-ethoxycoumarin and 7-ethoxyresorufin, are oxidized by all individual CYP1 forms. The main aim of this study was to find profluorescent coumarin substrates that are more selective for the individual CYP1 forms. Eleven 3-phenylcoumarin derivatives were synthetized, their enzyme kinetic parameters were determined, and their interactions in the active sites of CYP1 enzymes were analyzed by docking and molecular dynamic simulations. All coumarin derivatives and 7-ethoxyresorufin and 7-pentoxyresorufin were oxidized by at least one CYP1 enzyme. 3-(3-Methoxyphenyl)-6-methoxycoumarin (19) was 7-O-demethylated by similar high efficiency [21–30 ML/(min·mol CYP)] by all CYP1 forms and displayed similar binding in the enzyme active sites. 3-(3-Fluoro-4-acetoxyphenyl)coumarin (14) was selectively 7-O-demethylated by CYP1A1, but with low efficiency [0.16 ML/(min mol)]. This was explained by better orientation and stronger H-bond interactions in the active site of CYP1A1 than that of CYP1A2 and CYP1B1. 3-(4-Acetoxyphenyl)-6-chlorocoumarin (20) was 7-O-demethylated most efficiently by CYP1B1 [53 ML/(min·mol CYP)], followed by CYP1A1 [16 ML/(min·mol CYP)] and CYP1A2 [0.6 ML/(min·mol CYP)]. Variations in stabilities of complexes between 20 and the individual CYP enzymes explained these differences. Compounds 14, 19, and 20 are candidates to replace traditional substrates in measuring activity of human CYP1 enzymes.peerReviewe

Virtual Screening Strategy to Identify Retinoic Acid-Related Orphan Receptor γt Modulators

Molecular docking is a key method used in virtual screening (VS) campaigns to identify small-molecule ligands for drug discovery targets. While docking provides a tangible way to understand and predict the protein-ligand complex formation, the docking algorithms are often unable to separate active ligands from inactive molecules in practical VS usage. Here, a novel docking and shape-focused pharmacophore VS protocol is demonstrated for facilitating effective hit discovery using retinoic acid receptor-related orphan receptor gamma t (RORγt) as a case study. RORγt is a prospective target for treating inflammatory diseases such as psoriasis and multiple sclerosis. First, a commercial molecular database was flexibly docked. Second, the alternative docking poses were rescored against the shape/electrostatic potential of negative image-based (NIB) models that mirror the target’s binding cavity. The compositions of the NIB models were optimized via iterative trimming and benchmarking using a greedy search-driven algorithm or brute force NIB optimization. Third, a pharmacophore point-based filtering was performed to focus the hit identification on the known RORγt activity hotspots. Fourth, free energy binding affinity evaluation was performed on the remaining molecules. Finally, twenty-eight compounds were selected for in vitro testing and eight compounds were determined to be low μM range RORγt inhibitors, thereby showing that the introduced VS protocol generated an effective hit rate of ~29%