Mira Ahinko

Negative image-based (NIB) screening is a rigid molecular docking methodology that can also be employed in docking rescoring. During the NIB screening, a negative image is generated based on the target protein’s ligand-binding cavity by inverting its shape and electrostatics. The resulting NIB model is a drug-like entity or pseudo-ligand that is compared directly against ligand 3D conformers, as is done with a template compound in the ligand-based screening. This cavity-based rigid docking has been demonstrated to work with genuine drug targets in both benchmark testing and drug candidate/lead discovery. Firstly, the study explores in-depth the applicability of different ligand 3D conformer generation software for acquiring the best NIB screening results using cyclooxygenase-2 (COX-2) as the example system. Secondly, the entire NIB workflow from the protein structure preparation, model build-up, and ligand conformer generation to the similarity comparison is performed for COX-2. Accordingly, hands-on instructions are provided on how to employ the NIB methodology from start to finish, both with the rigid docking and docking rescoring using noncommercial software. The practical aspects of the NIB methodology, especially the effect of ligand conformers, are discussed thoroughly, thus, making the methodology accessible for new users.

1. Cytochrome P450 (CYP) enzymes constitute an essential xenobiotic metabolizing system that regulates the elimination of lipophilic compounds from the body. Convenient and affordable assays for CYP enzymes are important for assessing these metabolic pathways. 2. In this study, 10 novel profluorescent coumarin derivatives with various substitutions at carbons 3, 6 and 7 were developed. Molecular modeling indicated that 3-phenylcoumarin offers an excellent scaffold for the development of selective substrate compounds for various human CYP forms, as they could be metabolized to fluorescent 7-hydroxycoumarin derivatives. Oxidation of profluorescent coumarin derivatives to fluorescent metabolites by 13 important human liver xenobiotic-metabolizing CYP forms was determined by enzyme kinetic assays. 3. Four of the coumarin derivatives were converted to fluorescent metabolites by CYP1 family enzymes, with 6-methoxy-3-(4-trifluoromethylphenyl)coumarin being oxidized selectively by CYP1A2 in human liver microsomes. Another set of four compounds were metabolized by CYP2A6 and CYP1 enzymes. 7-Methoxy-3-(3-methoxyphenyl)coumarin was oxidized efficiently by CYP2C19 and CYP2D6 in a nonselective fashion. 4. The advantages of the novel substrates were 1) an excellent signal-to-background ratio, 2) selectivity for CYP1 forms, and 3) convenient multiwell plate measurement, allowing for precise determination of potential inhibitors of important human hepatic forms CYP1A2, CYP2C19 and CYP2D6.

The estimation of the correct binding mode and affinity of a ligand into a target protein using computational methods is challenging. However, docking can introduce poses from which the correct binding mode could be identified using other methods. Here, we analyzed the reliability of binding energy estimation using the molecular mechanics‐generalized Born surface area (MMGBSA) method without and with energy minimization to identify the likely ligand binding modes within docking results. MMGBSA workflow (a) outperformed docking in recognizing the correct binding modes of androgen receptor ligands and (b) improved the correlation coefficient of computational and experimental results of rescored docking poses to phosphodiesterase 4B. Combined with stability and atomic distance analysis, MMGBSA helped to (c) identify the binding modes and sites of metabolism of cytochrome P450 2A6 substrates. The standard deviation of estimated binding energy within one simulation was lowered by minimization in all three example cases. Minimization improved the identification of the correct binding modes of androgen receptor ligands. Although only three case studies are shown, the results are analogous and indicate that these behaviors could be generalized. Such identified binding modes could be further used, for example, with free energy perturbation methods to understand binding energetics more accurately.

Despite the large computational costs of molecular docking, the default scoring functions are often unable to recognize the active hits from the inactive molecules in large-scale virtual screening experiments. Thus, even though a correct binding pose might be sampled during the docking, the active compound or its biologically relevant pose is not necessarily given high enough score to arouse the attention. Various rescoring and post-processing approaches have emerged for improving the docking performance. Here, it is shown that the very early enrichment (number of actives scored higher than 1% of the highest ranked decoys) can be improved on average 2.5-fold or even 8.7-fold by comparing the docking-based ligand conformers directly against the target protein's cavity shape and electrostatics. The similarity comparison of the conformers is performed without geometry optimization against the negative image of the target protein's ligand-binding cavity using the negative image-based (NIB) screening protocol. The viability of the NIB rescoring or the R-NiB, pioneered in this study, was tested with 11 target proteins using benchmark libraries. By focusing on the shape/electrostatics complementarity of the ligand-receptor association, the R-NiB is able to improve the early enrichment of docking essentially without adding to the computing cost. By implementing consensus scoring, in which the R-NiB and the original docking scoring are weighted for optimal outcome, the early enrichment is improved to a level that facilitates effective drug discovery. Moreover, the use of equal weight from the original docking scoring and the R-NiB scoring improves the yield in most cases.

Monoamine oxidase B (MAO-B) catalyzes deamination of monoamines such as neurotransmitters dopamine and norepinephrine. Accordingly, small-molecule MAO-B inhibitors potentially alleviate the symptoms of dopamine-linked neuropathologies such as depression or Parkinson’s disease. Coumarin with a functionalized 3-phenyl ring system is a promising scaffold for building potent MAO-B inhibitors. Here, a vast set of 3-phenylcoumarin derivatives was designed using virtual combinatorial chemistry or rationally de novo and synthesized using microwave chemistry. The derivatives inhibited the MAO-B at 100 nM - 1 µM. The IC50 value of the most potent derivative 1 was 56 nM. A docking-based structure-activity relationship analysis summarizes the atom-level determinants of the MAO-B inhibition by the derivatives. Finally, the cross-reactivity of the derivatives was tested against monoamine oxidase A and a specific subset of enzymes linked to estradiol metabolism, known to have coumarin-based inhibitors. Overall, the results indicate that the 3-phenylcoumarins, especially derivative 1, present unique pharmacological features worth considering in future drug development.

A comprehensive set of 3-phenylcoumarin analogues with polar substituents was synthesised for blocking oestradiol synthesis by 17-β-hydroxysteroid dehydrogenase 1 (HSD1) in the latter part of the sulphatase pathway. Five analogues produced ≥62% HSD1 inhibition at 5 µM and, furthermore, three of them produced ≥68% inhibition at 1 µM. A docking-based structure-activity relationship analysis was done to determine the molecular basis of the inhibition and the cross-reactivity of the analogues was tested against oestrogen receptor, aromatase, cytochrome P450 1A2, and monoamine oxidases. Most of the analogues are only modestly active with 17-β-hydroxysteroid dehydrogenase 2 – a requirement for lowering effective oestradiol levels in vivo. Moreover, the analysis led to the synthesis and discovery of 3-imidazolecoumarin as a potent aromatase inhibitor. In short, coumarin core can be tailored with specific ring and polar moiety substitutions to block either the sulphatase pathway or the aromatase pathway for treating breast cancer and endometriosis.