Comparing Ligand Interactions with Multiple Receptors via Serial Docking

Molecular docking-based virtual screening: Challenges in hits identification for Anti-SARS-Cov-2 activity

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-Cov-2) requires finding new drugs or repurposing drugs for clinical use. Molecular docking belongs to structure-based drug design providing a fast method for identifying the hit compounds with antiviral activity against SARS-Cov-2. However, the weakness of the docking method is compounded by the limited crystallographic information and comparison drugs due to the novelty of this virus can present challenges in identifying hits of anti-SARS-Cov-2. In the current review, we highlighted several aspects, especially those related to the target structure, docking validation, and virtual hit selection, that need to be considered to obtain reliable docking results. Here, we discussed several cases pertaining to the issue highlighted and approaches that could be used to solve them.

Structural Elucidation of Rift Valley Fever Virus L Protein towards the Discovery of Its Potential Inhibitors

Rift valley fever virus (RVFV) is the causative agent of a viral zoonosis that causes a significant clinical burden in domestic and wild ruminants. Major outbreaks of the virus occur in livestock, and contaminated animal products or arthropod vectors can transmit the virus to humans. The viral RNA-dependent RNA polymerase (RdRp; L protein) of the RVFV is responsible for viral replication and is thus an appealing drug target because no effective and specific vaccine against this virus is available. The current study reported the structural elucidation of the RVFV-L protein by in-depth homology modeling since no crystal structure is available yet. The inhibitory binding modes of known potent L protein inhibitors were analyzed. Based on the results, further molecular docking-based virtual screening of Selleckchem Nucleoside Analogue Library (156 compounds) was performed to find potential new inhibitors against the RVFV L protein. ADME (Absorption, Distribution, Metabolism, and Excretion) and toxicity analysis of these compounds was also performed. Besides, the binding mechanism and stability of identified compounds were confirmed by a 50 ns molecular dynamic (MD) simulation followed by MM/PBSA binding free energy calculations. Homology modeling determined a stable multi-domain structure of L protein. An analysis of known L protein inhibitors, including Monensin, Mycophenolic acid, and Ribavirin, provide insights into the binding mechanism and reveals key residues of the L protein binding pocket. The screening results revealed that the top three compounds, A-317491, Khasianine, and VER155008, exhibited a high affinity at the L protein binding pocket. ADME analysis revealed good pharmacodynamics and pharmacokinetic profiles of these compounds. Furthermore, MD simulation and binding free energy analysis endorsed the binding stability of potential compounds with L protein. In a nutshell, the present study determined potential compounds that may aid in the rational design of novel inhibitors of the RVFV L protein as anti-RVFV drugs.

In silico structural elucidation of RNA-dependent RNA polymerase towards the identification of potential Crimean-Congo Hemorrhagic Fever Virus inhibitors

The Crimean-Congo Hemorrhagic Fever virus (CCHFV) is a segmented negative single-stranded RNA virus (-ssRNA) which causes severe hemorrhagic fever in humans with a mortality rate of ~50%. To date, no vaccine has been approved. Treatment is limited to supportive care with few investigational drugs in practice. Previous studies have identified viral RNA dependent RNA Polymerase (RdRp) as a potential drug target due to its significant role in viral replication and transcription. Since no crystal structure is available yet, we report the structural elucidation of CCHFV-RdRp by in-depth homology modeling. Even with low sequence identity, the generated model suggests a similar overall structure as previously reported RdRps. More specifically, the model suggests the presence of structural/functional conserved RdRp motifs for polymerase function, the configuration of uniform spatial arrangement of core RdRp sub-domains, and predicted positively charged entry/exit tunnels, as seen in sNSV polymerases. Extensive pharmacophore modeling based on per-residue energy contribution with investigational drugs allowed the concise mapping of pharmacophoric features and identified potential hits. The combination of pharmacophoric features with interaction energy analysis revealed functionally important residues in the conserved motifs together with in silico predicted common inhibitory binding modes with highly potent reference compounds.

Modelling three-dimensional protein structures for applications in drug design

A structural perspective of drug target and anti-target proteins, and their molecular interactions with biologically active molecules, largely advances many areas of drug discovery, including target validation, hit and lead finding and lead optimisation. In the absence of experimental 3D structures, protein structure prediction often offers a suitable alternative to facilitate structure-based studies. This review outlines recent methodical advances in homology modelling, with a focus on those techniques that necessitate consideration of ligand binding. In this context, model quality estimation deserves special attention because the accuracy and reliability of different structure prediction techniques vary considerably, and the quality of a model ultimately determines its usefulness for structure-based drug discovery. Examples of G-protein-coupled receptors (GPCRs) and ADMET-related proteins were selected to illustrate recent progress and current limitations of protein structure prediction. Basic guidelines for good modelling practice are also provided.

Comparative Modeling of Drug Target Proteins☆

In this perspective, we begin by describing the comparative protein structure modeling technique and the accuracy of the corresponding models. We then discuss the significant role that comparative prediction plays in drug discovery. We focus on virtual ligand screening against comparative models and illustrate the state-of-the-art by a number of specific examples.

Are homology models sufficiently good for free-energy simulations?

In this paper, I evaluate the usefulness of protein homology models in rigorous free-energy simulations to determine ligand affinities. Two templates were used to create models of the factor Xa protein and one template was used for dihydrofolate reductase from Plasmodium falciparum. Then, the relative free energies for several pairs of ligands were estimated using thermodynamic integration with the homology models as starting point of the simulation. These binding affinities were compared to affinities obtained when using published crystal structures as starting point of the simulations. Encouragingly, the differences between the affinities obtained when starting from either homology models or crystal structure were not statistical significant for a majority of the considered pairs of ligands. Differences between 1 and 2 kJ/mol were observed for the dihydrofolate reductase ligands and differences between 0 and 8 kJ/mol were observed for the factor Xa ligands. The largest difference for factor Xa was caused by an erroneous modeling of a loop region close to two of the ligands, and it was only observed when using one of the templates. Therefore, it is advisible to always use more than one template when creating homology models if they should be used in free-energy simulations.

Virtual ligand screening against comparative protein structure models

Virtual ligand screening uses computation to discover new ligands of a protein by screening one or more of its structural models against a database of potential ligands. Comparative protein structure modeling extends the applicability of virtual screening beyond the atomic structures determined by X-ray crystallography or NMR spectroscopy. Here, we describe an integrated modeling and docking protocol, combining comparative modeling by MODELLER and virtual ligand screening by DOCK.

Identification of the Schistosoma mansoni molecular target for the antimalarial drug artemether

Plasmodium falciparum and Schistosoma mansonii are the parasites responsible for most of the malaria and schistosomiasis cases in the world. Notwithstanding their many differences, the two agents have striking similarities in that they both are blood feeders and are targets of an overlapping set of drugs, including the well-known artemether molecule. Here we explore the possibility of using the known information about the mode of action of artemether in Plasmodium to identify the molecular target of the drug in Schistosoma and provide evidence that artemether binds to SmSERCA, a putative Ca²⁺-ATPase of Schistosoma . We also predict the putative binding mode of the molecule for both its Plasmodium and Schistosoma targets. Our analysis of the mode of binding of artemether to Ca²⁺-ATPases also provides an explanation for the apparent paradox that, although the molecule has no side effect in humans, it has been shown to possess antitumoral activity.

Molecular docking screens using comparative models of proteins

Two orders of magnitude more protein sequences can be modeled by comparative modeling than have been determined by X-ray crystallography and NMR spectroscopy. Investigators have nevertheless been cautious about using comparative models for ligand discovery because of concerns about model errors. We suggest how to exploit comparative models for molecular screens, based on docking against a wide range of crystallographic structures and comparative models with known ligands. To account for the variation in the ligand-binding pocket as it binds different ligands, we calculate "consensus" enrichment by ranking each library compound by its best docking score against all available comparative models and/or modeling templates. For the majority of the targets, the consensus enrichment for multiple models was better than or comparable to that of the holo and apo X-ray structures. Even for single models, the models are significantly more enriching than the template structure if the template is paralogous and shares more than 25% sequence identity with the target.

Modeling of peroxide activation in artemisinin derivatives by serial docking

Serial docking of artemisinin derivatives. A serial docking study was undertaken with the purpose to improve the understanding of the mechanism of production of biological activity in Artemisinin derivatives. The Heme molecule receptors were primarily chosen to represent the changing binding and oxidation states of this molecule, which are postulated to occur during the activation of the drug, in order to relate these results to the observed biological activity. The results of the docking runs were classified according to similarity to a standard orientation by a combination of automated and "by hand" procedures. One- and two-dimensional QSAR equations were used in an exploratory sense in order to study the influence of the type of receptor. Principal component and partial least squares regression techniques were used in the case of the multivariate (3D-QSAR) descriptors. The results obtained corroborate the postulated mechanism of production of biological activity as well as providing evidence that, at the moment of activation, the electronic structure of the Heme molecule approximates that of the oxygenated Heme, the drug molecule adopts a preferred orientation, and that, in addition to the important positive contribution of the O(1) - Fe interaction, there is as well a significant negative effect on the biological activity by carbons 4 - 6 of the artemisinin ring system.

Exploring Experimental Sources of Multiple Protein Conformations in Structure-Based Drug Design

Receptor flexibility must be incorporated into structure-based drug design in order to portray a more accurate representation of a protein in solution. Our approach is to generate pharmacophore models based on multiple conformations of a protein and is very similar to solvent mapping of hot spots. Previously, we had success using computer-generated conformations of apo human immunodeficiency virus-1 protease (HIV-1p). Here, we examine the use of an NMR ensemble versus a collection of crystal structures, and we compare back to our previous study based on computer-generated conformations. To our knowledge, this is the first direct comparison of an NMR ensemble and a collection of crystal structures to incorporate protein flexibility in structure-based drug design. To provide an accurate comparison between the experimental sources, we used bound structures for our multiple protein structure (MPS) pharmacophore models. The models from an NMR ensemble and a collection of crystal structures were both able to discriminate known HIV-1p inhibitors from decoy molecules and displayed superior performance over models created from single conformations of the protein. Although the active-site conformations were already predefined by bound ligands, the use of MPS allows us to overcome the cross-docking problem and generate a model that does not simply reproduce the chemical characteristics of a specific ligand class. We show that there is more structural variation between 28 structures in an NMR ensemble than 90 crystal structures bound to a variety of ligands. MPS models from both sources performed well, but the model determined using the NMR ensemble appeared to be the most general yet accurate representation of the active site. This work encourages the use of NMR models in structure-based design.

Analysis of HIV wild-type and mutant structures via in silico docking against diverse ligand libraries

The FightAIDS@Home distributed computing project uses AutoDock for an initial virtual screen of HIV protease structures against a broad range of 1771 ligands including both known protease inhibitors and a diverse library of other ligands. The volume of results allows novel large-scale analyses of binding energy "profiles" for HIV structures. Beyond identifying potential lead compounds, these characterizations provide methods for choosing representative wild-type and mutant protein structures from the larger set. From the binding energy profiles of the PDB structures, a principal component analysis based analysis identifies seven "spanning" proteases. A complementary analysis finds that the wild-type protease structure 2BPZ best captures the central tendency of the protease set. Using a comparison of known protease inhibitors against the diverse ligand set yields an AutoDock binding energy "significance" threshold of -7.0 kcal/mol between significant, strongly binding ligands and other weak/nonspecific binding energies. This threshold captures nearly 98% of known inhibitor interactions while rejecting more than 95% of suspected noninhibitor interactions. These methods should be of general use in virtual screening projects and will be used to improve further FightAIDS@Home experiments.

Predicting substrates by docking high-energy intermediates to enzyme structures

With the emergence of sequences and even structures for proteins of unknown function, structure-based prediction of enzyme activity has become a pragmatic as well as an interesting question. Here we investigate a method to predict substrates for enzymes of known structure by docking high-energy intermediate forms of the potential substrates. A database of such high-energy transition-state analogues was created from the KEGG metabolites. To reduce the number of possible reactions to consider, we restricted ourselves to enzymes of the amidohydrolase superfamily. We docked each metabolite into seven different amidohydrolases in both the ground-state and the high-energy intermediate forms. Docking the high-energy intermediates improved the discrimination between decoys and substrates significantly over the corresponding standard ground-state database, both by enrichment of the true substrates and by geometric fidelity. To test this method prospectively, we attempted to predict the enantioselectivity of a set of chiral substrates for phosphotriesterase, for both wild-type and mutant forms of this enzyme. The stereoselectivity ratios of the six enzymes considered for those four substrate enantiomer pairs differed over a range of 10- to 10,000-fold and underwent 20 switches in stereoselectivities for favored enantiomers, compared to the wild type. The docking of the high-energy intermediates correctly predicted the stereoselectivities for 18 of the 20 substrate/enzyme combinations when compared to subsequent experimental synthesis and testing. The possible applications of this approach to other enzymes are considered.

Progressive Docking: A Hybrid QSAR/Docking Approach for Accelerating In Silico High Throughput Screening

A combination of protein-ligand docking and ligand-based QSAR approaches has been elaborated, aiming to speed-up the process of virtual screening. In particular, this approach utilizes docking scores generated for already processed compounds to build predictive QSAR models that, in turn, assess hypothetical target binding affinities for yet undocked entries. The "progressive docking" has been tested on drug-like substances from the NCI database that have been docked into several unrelated targets, including human sex hormone binding globulin (SHBG), carbonic anhydrase, corticosteroid-binding globulin, SARS 3C-like protease, and HIV1 reverse transcriptase. We demonstrate that progressive docking can reduce the amount of computations 1.2- to 2.6-fold (when compared to traditional docking), while maintaining 80-99% hit recovery rates. This progressive-docking procedure, therefore, substantially accelerates high throughput screening, especially when using high accuracy (slower) docking approaches and large-sized datasets, and has allowed us to identify several novel potent nonsteroidal SHBG ligands.

InSilico-Guided Target Identification of a Scaffold-Focused Library: 1,3,5-Triazepan-2,6-diones as Novel Phospholipase A2 Inhibitors

A collection of 2150 druggable active sites from the Protein Data Bank was screened by high-throughput docking to identify putative targets for five representative molecules of a combinatorial library sharing a 1,3,5-triazepan-2,6-dione scaffold. Five targets were prioritized for experimental evaluation by computing enrichment in individual protein entries among the top 2% scoring targets. Out of the five proposed proteins, secreted phospholipase A2 (sPLA2) was shown to be a true target for a panel of 1,3,5-triazepan-2,6-diones which exhibited micromolar affinities toward two human sPLA2 members.

Screening Drug-Like Compounds by Docking to Homology Models: A Systematic Study

In the absence of an experimentally solved structure, a homology model of a protein target can be used instead for virtual screening of drug candidates by docking and scoring. This approach poses a number of questions regarding the choice of the template to use in constructing the model, the accuracy of the screening results, and the importance of allowing for protein flexibility. The present study addresses such questions with compound screening calculations for multiple homology models of five drug targets. A central result is that docking to homology models frequently yields enrichments of known ligands as good as that obtained by docking to a crystal structure of the actual target protein. Interestingly, however, standard measures of the similarity of the template used to build the homology model to the targeted protein show little correlation with the effectiveness of the screening calculations, and docking to the template itself often is as successful as docking to the corresponding homology model. Treating key side chains as mobile produces a modest improvement in the results. The reasons for these sometimes unexpected results, and their implications for future methodologic development, are discussed.

Rapid computational identification of the targets of protein kinase inhibitors

We describe a method for rapidly computing the relative affinities of an inhibitor for all individual members of a family of homologous receptors. The approach, implemented in a new program, SCR, models inhibitor-receptor interactions in full atomic detail with an empirical energy function and includes an explicit account of flexibility in homology-modeled receptors through sampling of libraries of side chain rotamers. SCR's general utility was demonstrated by application to seven different protein kinase inhibitors: for each inhibitor, relative binding affinities with panels of approximately 20 protein kinases were computed and compared with experimental data. For five of the inhibitors (SB203580, purvalanol B, imatinib, H89, and hymenialdisine), SCR provided excellent reproduction of the experimental trends and, importantly, was capable of identifying the targets of inhibitors even when they belonged to different kinase families. The method's performance in a predictive setting was demonstrated by performing separate training and testing applications, and its key assumptions were tested by comparison with a number of alternative approaches employing the ligand-docking program AutoDock (Morris et al. J. Comput. Chem. 1998, 19, 1639-1662). These comparison tests included using AutoDock in nondocking and docking modes and performing energy minimizations of inhibitor-kinase complexes with the molecular mechanics code GROMACS (Berendsen et al. Comput. Phys. Commun. 1995, 91, 43-56). It was found that a surprisingly important aspect of SCR's approach is its assumption that the inhibitor be modeled in the same orientation for each kinase: although this assumption is in some respects unrealistic, calculations that used apparently more realistic approaches produced clearly inferior results. Finally, as a large-scale application of the method, SB203580, purvalanol B, and imatinib were screened against an almost full complement of 493 human protein kinases using SCR in order to identify potential new targets; the predicted targets of SB203580 were compared with those identified in recent proteomics-based experiments. These kinome-wide screens, performed within a day on a small cluster of PCs, indicate that explicit computation of inhibitor-receptor binding affinities has the potential to promote rapid discovery of new therapeutic targets for existing inhibitors.