Information decay in molecular docking screens against holo, apo, and modeled conformations of enzymes

NMR structures of small molecules bound to a model of a CUG RNA repeat expansion

Trinucleotide repeat expansions fold into long, stable hairpins and cause a variety of incurable RNA gain-of-function diseases such as Huntington's disease, the myotonic dystrophies, and spinocerebellar ataxias. One approach for treating these diseases is to bind small molecules to these structured RNAs. Both Huntington's disease-like 2 (HDL2) and myotonic dystrophy type 1 (DM1) are caused by a r(CUG) repeat expansion, or r(CUG)exp. The RNA folds into a hairpin structure with a periodic array of 1 × 1 nucleotide UU loops (5'CUG/3'GUC; where the underlined nucleotides indicate the Us in the internal loop) that sequester various RNA-binding proteins (RBPs) and hence the source of its gain-of-function. Here, we report nuclear magnetic resonance (NMR)-refined structures of single 5'CUG/3'GUC motifs in complex with three different small molecules, a di-guandinobenzoate (1), a derivative of 1 where the guanidino groups have been exchanged for imidazole (2), and a quinoline with improved drug-like properties (3). These structures were determined using NMR spectroscopy and simulated annealing with restrained molecular dynamics (MD). Compounds 1, 2, and 3 formed stacking and hydrogen bonding interactions with the 5'CUG/3'GUC motif. Compound 3 also formed van der Waals interactions with the internal loop. The global structure of each RNA-small molecule complexes retains an A-form conformation, while the internal loops are still dynamic but to a lesser extent compared to the unbound form. These results aid our understanding of ligand-RNA interactions and enable structure-based design of small molecules with improved binding affinity for and biological activity against r(CUG)exp. As the first ever reported structures of a r(CUG) repeat bound to ligands, these structures can enable virtual screening campaigns combined with machine learning assisted de novo design.

Structural basis of sugar recognition by SCFFBS2 ubiquitin ligase involved in NGLY1 deficiency

The cytosolic peptide:N-glycanase (PNGase) is involved in the quality control of N-glycoproteins via the endoplasmic reticulum-associated degradation (ERAD) pathway. Mutations in the gene encoding cytosolic PNGase (NGLY1 in humans) cause NGLY1 deficiency. Recent findings indicate that the F-box protein FBS2 of the SCFFBS2 ubiquitin ligase complex can be a promising drug target for NGLY1 deficiency. Here, we determined the crystal structure of bovine FBS2 complexed with the adaptor protein SKP1 and a sugar ligand, Man3GlcNAc2, which corresponds to the core pentasaccharide of N-glycan. Our crystallographic data together with NMR data revealed the structural basis of disparate sugar-binding specificities in homologous FBS proteins and identified a potential druggable pocket for in silico docking studies. Our results provide a potential basis for the development of selective inhibitors against FBS2 in NGLY1 deficiency.

NMR structures of small molecules bound to a model of an RNA CUG repeat expansion

Trinucleotide repeat expansions fold into long, stable hairpins and cause a variety of incurable RNA gain-of-function diseases such as Huntington's disease, the myotonic dystrophies, and spinocerebellar ataxias. One approach for treating these diseases is to bind small molecules to the structured RNAs. Both Huntington's disease-like 2 (HDL2) and myotonic dystrophy type 1 (DM1) are caused by a r(CUG) repeat expansion, or r(CUG)exp. The RNA folds into a hairpin structure with a periodic array of 1×1 nucleotide UU loops (5'CUG/3'GUC; where the underlined nucleotides indicate the Us in the internal loop) that sequester various RNA-binding proteins (RBP) and hence the source of its gain-of-function. Here, we report NMR-refined structures of single 5'CUG/3'GUC motifs in complex with three different small molecules, a di-guandinobenzoate (1), a derivative of 1 where the guanidino groups have been exchanged for imidazole (2), and a quinoline with improved drug-like properties (3). These structures were determined using nuclear magnetic resonance (NMR) spectroscopy and simulated annealing with restrained molecular dynamics (MD). Compounds 1, 2, and 3 formed stacking and hydrogen bonding interactions with the 5'CUG/3'GUC motif. Compound 3 also formed van der Waals interactions with the internal loop. The global structure of each RNA-small molecule complexes retains an A-form conformation, while the internal loops are still dynamic but to a lesser extent compared to the unbound form. These results aid our understanding of ligand-RNA interactions and enable structure-based design of small molecules with improved binding affinity for and biological activity against r(CUG)exp. As the first ever reported structures of RNA r(CUG) repeats bound to ligands, these structures can enable virtual screening campaigns combined with machine learning assisted de novo design.

Computational evaluation and benchmark study of 342 crystallographic holo-structures of SARS-CoV-2 Mpro enzyme

The coronavirus disease 19 pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has led to a global health crisis with millions of confirmed cases and related deaths. The main protease (Mpro) of SARS-CoV-2 is crucial for viral replication and presents an attractive target for drug development. Despite the approval of some drugs, the search for effective treatments continues. In this study, we systematically evaluated 342 holo-crystal structures of Mpro to identify optimal conformations for structure-based virtual screening (SBVS). Our analysis revealed limited structural flexibility among the structures. Three docking programs, AutoDock Vina, rDock, and Glide were employed to assess the efficiency of virtual screening, revealing diverse performances across selected Mpro structures. We found that the structures 5RHE, 7DDC, and 7DPU (PDB Ids) consistently displayed the lowest EF, AUC, and BEDROCK scores. Furthermore, these structures demonstrated the worst pose prediction results in all docking programs. Two structural differences contribute to variations in docking performance: the absence of the S1 subsite in 7DDC and 7DPU, and the presence of a subpocket in the S2 subsite of 7DDC, 7DPU, and 5RHE. These findings underscore the importance of selecting appropriate Mpro conformations for SBVS, providing valuable insights for advancing drug discovery efforts.

DiffBindFR: an SE(3) equivariant network for flexible protein-ligand docking

Molecular docking, a key technique in structure-based drug design, plays pivotal roles in protein-ligand interaction modeling, hit identification and optimization, in which accurate prediction of protein-ligand binding mode is essential. Conventional docking approaches perform well in redocking tasks with known protein binding pocket conformation in the complex state. However, in real-world docking scenario without knowing the protein binding conformation for a new ligand, accurately modeling the binding complex structure remains challenging as flexible docking is computationally expensive and inaccurate. Typical deep learning-based docking methods do not explicitly consider protein side chain conformations and fail to ensure the physical plausibility and detailed atomic interactions. In this study, we present DiffBindFR, a full-atom diffusion-based flexible docking model that operates over the product space of ligand overall movements and flexibility and pocket side chain torsion changes. We show that DiffBindFR has higher accuracy in producing native-like binding structures with physically plausible and detailed interactions than available docking methods. Furthermore, in the Apo and AlphaFold2 modeled structures, DiffBindFR demonstrates superior advantages in accurate ligand binding pose and protein binding conformation prediction, making it suitable for Apo and AlphaFold2 structure-based drug design. DiffBindFR provides a powerful flexible docking tool for modeling accurate protein-ligand binding structures.

Evaluation of AlphaFold2 Structures for Hit Identification across Multiple Scenarios

The introduction of AlphaFold2 (AF2) has sparked significant enthusiasm and generated extensive discussion within the scientific community, particularly among drug discovery researchers. Although previous studies have addressed the performance of AF2 structures in virtual screening (VS), a more comprehensive investigation is still necessary considering the paramount importance of structural accuracy in drug design. In this study, we evaluate the performance of AF2 structures in VS across three common drug discovery scenarios: targets with holo, apo, and AF2 structures; targets with only apo and AF2 structures; and targets exclusively with AF2 structures. We utilized both the traditional physics-based Glide and the deep-learning-based scoring function RTMscore to rank the compounds in the DUD-E, DEKOIS 2.0, and DECOY data sets. The results demonstrate that, overall, the performance of VS on AF2 structures is comparable to that on apo structures but notably inferior to that on holo structures across diverse scenarios. Moreover, when a target has solely AF2 structure, selecting the holo structure of the target from different subtypes within the same protein family produces comparable results with the AF2 structure for VS on the data set of the AF2 structures, and significantly better results than the AF2 structures on its own data set. This indicates that utilizing AF2 structures for docking-based VS may not yield most satisfactory outcomes, even when solely AF2 structures are available. Moreover, we rule out the possibility that the variations in VS performance between the binding pockets of AF2 and holo structures arise from the differences in their biological assembly composition.

Rational Design of Drugs Targeting G-Protein-Coupled Receptors: Ligand Search and Screening

G protein-coupled receptors (GPCRs) are transmembrane proteins that participate in many physiological processes and represent major pharmacological targets. Recent advances in structural biology of GPCRs have enabled the development of drugs based on the receptor structure (structure-based drug design, SBDD). SBDD utilizes information about the receptor-ligand complex to search for suitable compounds, thus expanding the chemical space of possible receptor ligands without the need for experimental screening. The review describes the use of structure-based virtual screening (SBVS) for GPCR ligands and approaches for the functional testing of potential drug compounds, as well as discusses recent advances and successful examples in the application of SBDD for the identification of GPCR ligands.

Discovery of new antiviral agents through artificial intelligence: In vitro and in vivo results

In this research, we employed a deep reinforcement learning (RL)-based molecule design platform to generate a diverse set of compounds targeting the neuraminidase (NA) of influenza A and B viruses. A total of 60,291 compounds were generated, of which 86.5 % displayed superior physicochemical properties compared to oseltamivir. After narrowing down the selection through computational filters, nine compounds with non-sialic acid-like structures were selected for in vitro experiments. We identified two compounds, DS-22-inf-009 and DS-22-inf-021 that effectively inhibited the NAs of both influenza A and B viruses (IAV and IBV), including H275Y mutant strains at low micromolar concentrations. Molecular dynamics simulations revealed a similar pattern of interaction with amino acid residues as oseltamivir. In cell-based assays, DS-22-inf-009 and DS-22-inf-021 inhibited IAV and IBV in a dose-dependent manner with EC50 values ranging from 0.29 μM to 2.31 μM. Furthermore, animal experiments showed that both DS-22-inf-009 and DS-22-inf-021 exerted antiviral activity in mice, conferring 65 % and 85 % protection from IAV (H1N1 pdm09), and 65 % and 100 % protection from IBV (Yamagata lineage), respectively. Thus, these findings demonstrate the potential of RL to generate compounds with promising antiviral properties.

Automated benchmarking of combined protein structure and ligand conformation prediction

The prediction of protein-ligand complexes (PLC), using both experimental and predicted structures, is an active and important area of research, underscored by the inclusion of the Protein-Ligand Interaction category in the latest round of the Critical Assessment of Protein Structure Prediction experiment CASP15. The prediction task in CASP15 consisted of predicting both the three-dimensional structure of the receptor protein as well as the position and conformation of the ligand. This paper addresses the challenges and proposed solutions for devising automated benchmarking techniques for PLC prediction. The reliability of experimentally solved PLC as ground truth reference structures is assessed using various validation criteria. Similarity of PLC to previously released complexes are employed to judge PLC diversity and the difficulty of a PLC as a prediction target. We show that the commonly used PDBBind time-split test-set is inappropriate for comprehensive PLC evaluation, with state-of-the-art tools showing conflicting results on a more representative and high quality dataset constructed for benchmarking purposes. We also show that redocking on crystal structures is a much simpler task than docking into predicted protein models, demonstrated by the two PLC-prediction-specific scoring metrics created. Finally, we introduce a fully automated pipeline that predicts PLC and evaluates the accuracy of the protein structure, ligand pose, and protein-ligand interactions.

Equivariant Flexible Modeling of the Protein-Ligand Binding Pose with Geometric Deep Learning

Flexible modeling of the protein-ligand complex structure is a fundamental challenge for in silico drug development. Recent studies have improved commonly used docking tools by incorporating extra-deep learning-based steps. However, such strategies limit their accuracy and efficiency because they retain massive sampling pressure and lack consideration for flexible biomolecular changes. In this study, we propose FlexPose, a geometric graph network capable of direct flexible modeling of complex structures in Euclidean space without the following conventional sampling and scoring strategies. Our model adopts two key designs: scalar-vector dual feature representation and SE(3)-equivariant network, to manage dynamic structural changes, as well as two strategies: conformation-aware pretraining and weakly supervised learning, to boost model generalizability in unseen chemical space. Benefiting from these paradigms, our model dramatically outperforms all tested popular docking tools and recently advanced deep learning methods, especially in tasks involving protein conformation changes. We further investigate the impact of protein and ligand similarity on the model performance with two conformation-aware strategies. Moreover, FlexPose provides an affinity estimation and model confidence for postanalysis.

Discovery of a novel cardiac-specific myosin modulator using artificial intelligence-based virtual screening

Direct modulation of cardiac myosin function has emerged as a therapeutic target for both heart disease and heart failure. However, the development of myosin-based therapeutics has been hampered by the lack of targeted in vitro screening assays. In this study we use Artificial Intelligence-based virtual high throughput screening (vHTS) to identify novel small molecule effectors of human β-cardiac myosin. We test the top scoring compounds from vHTS in biochemical counter-screens and identify a novel chemical scaffold called 'F10' as a cardiac-specific low-micromolar myosin inhibitor. Biochemical and biophysical characterization in both isolated proteins and muscle fibers show that F10 stabilizes both the biochemical (i.e. super-relaxed state) and structural (i.e. interacting heads motif) OFF state of cardiac myosin, and reduces force and left ventricular pressure development in isolated myofilaments and Langendorff-perfused hearts, respectively. F10 is a tunable scaffold for the further development of a novel class of myosin modulators.

Worth the Weight: Sub-Pocket EXplorer (SubPEx), a Weighted Ensemble Method to Enhance Binding-Pocket Conformational Sampling

Structure-based virtual screening (VS) is an effective method for identifying potential small-molecule ligands, but traditional VS approaches consider only a single binding-pocket conformation. Consequently, they struggle to identify ligands that bind to alternate conformations. Ensemble docking helps address this issue by incorporating multiple conformations into the docking process, but it depends on methods that can thoroughly explore pocket flexibility. We here introduce Sub-Pocket EXplorer (SubPEx), an approach that uses weighted ensemble (WE) path sampling to accelerate binding-pocket sampling. As proof of principle, we apply SubPEx to three proteins relevant to drug discovery: heat shock protein 90, influenza neuraminidase, and yeast hexokinase 2. SubPEx is available free of charge without registration under the terms of the open-source MIT license: http://durrantlab.com/subpex/.

CHARMM-GUI-Based Induced Fit Docking Workflow to Generate Reliable Protein-Ligand Binding Modes

Molecular docking is a preferred method to predict ligand binding modes and their binding energy to target protein receptors, which is critical in early phase structure-based drug discovery. However, there is a persistent challenge in docking that can be attributed to the induced fit effect, as receptor binding sites undergo induced fit conformational changes upon ligand binding to achieve better binding modes. In this work, based on CHARMM-GUI LBS Finder& Refiner and High-Throughput Simulator, we present a straightforward CHARMM-GUI induced fit docking (CGUI-IFD) workflow to generate reliable protein-ligand binding modes. The CGUI-IFD workflow generates an ensemble of receptor binding site conformations through ligand-binding site (LBS) refinement, runs rigid receptor docking, and performs high-throughput molecular dynamics (MD) simulations of protein-ligand complex structures in explicit solvents. The results are evaluated based on the ligand root-mean-square deviation (RMSD)-based binding stability and the molecular mechanics generalized Born surface area binding energy. For a benchmark test, we used 258 cross-docking protein-ligand pairs across 41 target proteins from the Schrodinger IFD-MD data set. The application of CGUI-IFD on this data set shows 80% success rate (within 2.5 Å RMSD from the experimental structures). We expect that the CGUI-IFD workflow can be useful to generate reliable ligand binding modes for cross-docking cases.

Worth the weight: Sub-Pocket EXplorer (SubPEx), a weighted-ensemble method to enhance binding-pocket conformational sampling

Structure-based virtual screening (VS) is an effective method for identifying potential small-molecule ligands, but traditional VS approaches consider only a single binding-pocket conformation. Consequently, they struggle to identify ligands that bind to alternate conformations. Ensemble docking helps address this issue by incorporating multiple conformations into the docking process, but it depends on methods that can thoroughly explore pocket flexibility. We here introduce Sub-Pocket EXplorer (SubPEx), an approach that uses weighted ensemble (WE) path sampling to accelerate binding-pocket sampling. As proof of principle, we apply SubPEx to three proteins relevant to drug discovery: heat shock protein 90, influenza neuraminidase, and yeast hexokinase 2. SubPEx is available free of charge without registration under the terms of the open-source MIT license: http://durrantlab.com/subpex/.

Lessons Learnt from COVID-19: Computational Strategies for Facing Present and Future Pandemics

Since its outbreak in December 2019, the COVID-19 pandemic has caused the death of more than 6.5 million people around the world. The high transmissibility of its causative agent, the SARS-CoV-2 virus, coupled with its potentially lethal outcome, provoked a profound global economic and social crisis. The urgency of finding suitable pharmacological tools to tame the pandemic shed light on the ever-increasing importance of computer simulations in rationalizing and speeding up the design of new drugs, further stressing the need for developing quick and reliable methods to identify novel active molecules and characterize their mechanism of action. In the present work, we aim at providing the reader with a general overview of the COVID-19 pandemic, discussing the hallmarks in its management, from the initial attempts at drug repurposing to the commercialization of Paxlovid, the first orally available COVID-19 drug. Furthermore, we analyze and discuss the role of computer-aided drug discovery (CADD) techniques, especially those that fall in the structure-based drug design (SBDD) category, in facing present and future pandemics, by showcasing several successful examples of drug discovery campaigns where commonly used methods such as docking and molecular dynamics have been employed in the rational design of effective therapeutic entities against COVID-19.

Phytonutrient Inhibitors of SARS-CoV-2/NSP5-Encoded Main Protease (Mpro) Autocleavage Enzyme Critical for COVID-19 Pathogenesis

The genomic reshuffling, mutagenicity, and high transmission rate of the SARS-CoV-2 pathogen highlights an urgent need for effective antiviral interventions for COVID-19 control. Targeting the highly conserved viral genes and/or gene-encoded viral proteins such as main proteinase (M^pro), RNA-dependent RNA polymerase (RdRp) and helicases are plausible antiviral approaches to prevent replication and propagation of the SARS-CoV-2 infection. Coronaviruses (CoVs) are prone to extensive mutagenesis; however, any genetic alteration to its highly conserved M^pro enzyme is often detrimental to the viral pathogen. Therefore, inhibitors that target the M^pro enzyme could reduce the risk of mutation-mediated drug resistance and provide effective antiviral protection. Several existing antiviral drugs and dietary bioactives are currently repurposed to treat COVID-19. Dietary bioactives from three ayurvedic medicinal herbs, 18 β-glycyrrhetinic acid (ΔG = 8.86 kcal/mol), Solanocapsine (ΔG = 8.59 kcal/mol), and Vasicoline (ΔG = 7.34 kcal/mol), showed high-affinity binding to M^pro enzyme than the native N3 inhibitor (ΔG = 5.41 kcal/mol). Flavonoids strongly inhibited SARS-CoV-2 M^pro with comparable or higher potency than the antiviral drug, remdesivir. Several tannin hydrolysates avidly bound to the receptor-binding domain and catalytic dyad (His₄₁ and Cys₁₄₅) of SARS-CoV-2 M^pro through H-bonding forces. Quercetin binding to M^pro altered the thermostability of the viral protein through redox-based mechanism and inhibited the viral enzymatic activity. Interaction of quercetin-derivatives with the M^pro seem to be influenced by the 7-OH group and the acetoxylation of sugar moiety on the ligand molecule. Based on pharmacokinetic and ADMET profiles, several phytonutrients could serve as a promising redox nutraceutical for COVID-19 management.

Evaluation of AlphaFold2 structures as docking targets

AlphaFold2 is a promising new tool for researchers to predict protein structures and generate high-quality models, with low backbone and global root-mean-square deviation (RMSD) when compared with experimental structures. However, it is unclear if the structures predicted by AlphaFold2 will be valuable targets of docking. To address this question, we redocked ligands in the PDBbind datasets against the experimental co-crystallized receptor structures and against the AlphaFold2 structures using AutoDock-GPU. We find that the quality measure provided during structure prediction is not a good predictor of docking performance, despite accurately reflecting the quality of the alpha carbon alignment with experimental structures. Removing low-confidence regions of the predicted structure and making side chains flexible improves the docking outcomes. Overall, despite high-quality prediction of backbone conformation, fine structural details limit the naive application of AlphaFold2 models as docking targets.

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.

Holo Protein Conformation Generation from Apo Structures by Ligand Binding Site Refinement

An important part in structure-based drug design is the selection of an appropriate protein structure. It has been revealed that a holo protein structure that contains a well-defined binding site is a much better choice than an apo structure in structure-based drug discovery. Therefore, it is valuable to obtain a holo-like protein conformation from apo structures in the case where no holo structure is available. Meeting the need, we present a robust approach to generate reliable holo-like structures from apo structures by ligand binding site refinement with restraints derived from holo templates with low homology. Our method was tested on a test set of 32 proteins from the DUD-E data set and compared with other approaches. It was shown that our method successfully refined the apo structures toward the corresponding holo conformations for 23 of 32 proteins, reducing the average all-heavy-atom RMSD of binding site residues by 0.48 Å. In addition, when evaluated against all the holo structures in the protein data bank, our method can improve the binding site RMSD for 14 of 19 cases that experience significant conformational changes. Furthermore, our refined structures also demonstrate their advantages over the apo structures in ligand binding mode predictions by both rigid docking and flexible docking and in virtual screening on the database of active and decoy ligands from the DUD-E. These results indicate that our method is effective in recovering holo-like conformations and will be valuable in structure-based drug discovery.

Evaluating the accuracy of the AMBER protein force fields in modeling dihydrofolate reductase structures: misbalance in the conformational arrangements of the flexible loop domains

Protein flexible loop regions were once thought to be simple linkers between other more functional secondary structural elements. However, as it becomes clearer that these loop domains are critical players in a plethora of biological processes, accurate conformational sampling of 3D loop structures is vital to the advancement of drug design techniques and the overall growth of knowledge surrounding molecular systems. While experimental techniques provide a wealth of structural information, the resolution of flexible loop domains is sometimes low or entirely absent due to their complex and dynamic nature. This highlights an opportunity for de novo structure prediction using in silico methods with molecular dynamics (MDs). This study evaluates some of the AMBER protein force field's (ffs) ability to accurately model dihydrofolate reductase (DHFR) conformations, a protein complex characterized by specific arrangements and interactions of multiple flexible loops whose conformations are determined by the presence or absence of bound ligands and cofactors. Although the AMBER ffs, including ff19SB, studied well model most protein structures with rich secondary structure, results obtained here suggest the inability to significantly sample the expected DHFR loop-loop conformations - of the six distinct protein-ligand systems simulated, a majority lacked consistent stabilization of experimentally derived metrics definitive the three enzyme conformations. Although under-sampling and the chosen ff parameter combinations could be the cause, given past successes with these MD approaches for many protein systems, this suggests a potential misbalance in available ff parameters required to accurately predict the structure of multiple flexible loop regions present in proteins.Communicated by Ramaswamy H. Sarma.

Extensive Sampling of Molecular Dynamics Simulations to Identify Reliable Protein Structures for Optimized Virtual Screening Studies: The Case of the hTRPM8 Channel

(1) Background: Virtual screening campaigns require target structures in which the pockets are properly arranged for binding. Without these, MD simulations can be used to relax the available target structures, optimizing the fine architecture of their binding sites. Among the generated frames, the best structures can be selected based on available experimental data. Without experimental templates, the MD trajectories can be filtered by energy-based criteria or sampled by systematic analyses. (2) Methods: A blind and methodical analysis was performed on the already reported MD run of the hTRPM8 tetrameric structures; a total of 50 frames underwent docking simulations by using a set of 1000 ligands including 20 known hTRPM8 modulators. Docking runs were performed by LiGen program and involved the frames as they are and after optimization by SCRWL4.0. For each frame, all four monomers were considered. Predictive models were developed by the EFO algorithm based on the sole primary LiGen scores. (3) Results: On average, the MD simulation progressively enhances the performance of the extracted frames, and the optimized structures perform better than the non-optimized frames (EF1% mean: 21.38 vs. 23.29). There is an overall correlation between performances and volumes of the explored pockets and the combination of the best performing frames allows to develop highly performing consensus models (EF1% = 49.83). (4) Conclusions: The systematic sampling of the entire MD run provides performances roughly comparable with those previously reached by using rationally selected frames. The proposed strategy appears to be helpful when the lack of experimental data does not allow an easy selection of the optimal structures for docking simulations. Overall, the reported docking results confirm the relevance of simulating all the monomers of an oligomer structure and emphasize the efficacy of the SCRWL4.0 method to optimize the protein structures for docking calculations.

Selective CDK9 Inhibition by Natural Compound Toyocamycin in Cancer Cells

Simple Summary

By combining drug screens, transcriptomic studies, and in vitro assays, our study identified the natural product toyocamycin as a potent and selective CDK9 inhibitor. Thus, toyocamycin can be used as a new small molecule to modulate CDK9 activity in preclinical research. Through docking simulations, we identified its specific binding sites, which could spark some interest to design novel small molecule CDK9 inhibitors.

Abstract

Aberrant transcription in cancer cells involves the silencing of tumor suppressor genes (TSGs) and activation of oncogenes. Transcriptomic changes are associated with epigenomic alterations such as DNA-hypermethylation, histone deacetylation, and chromatin condensation in promoter regions of silenced TSGs. To discover novel drugs that trigger TSG reactivation in cancer cells, we used a GFP-reporter system whose expression is silenced by promoter DNA hypermethylation and histone deacetylation. After screening a natural product drug library, we identified that toyocamycin, an adenosine-analog, induces potent GFP reactivation and loss of clonogenicity in human colon cancer cells. Connectivity-mapping analysis revealed that toyocamycin produces a pharmacological signature mimicking cyclin-dependent kinase (CDK) inhibitors. RNA-sequencing revealed that the toyocamycin transcriptomic signature resembles that of a specific CDK9 inhibitor (HH1). Specific inhibition of RNA Pol II phosphorylation level and kinase assays confirmed that toyocamycin specifically inhibits CDK9 (IC₅₀ = 79 nM) with a greater efficacy than other CDKs (IC₅₀ values between 0.67 and 15 µM). Molecular docking showed that toyocamycin efficiently binds the CDK9 catalytic site in a conformation that differs from other CDKs, explained by the binding contribution of specific amino acids within the catalytic pocket and protein backbone. Altogether, we demonstrated that toyocamycin exhibits specific CDK9 inhibition in cancer cells, highlighting its potential for cancer chemotherapy.

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.

Lockdown, a selective small-molecule inhibitor of the integrin phosphatase PPM1F, blocks cancer cell invasion

Phosphatase PPM1F is a regulator of cell adhesion by fine-tuning integrin activity and actin cytoskeleton structures. Elevated expression of this enzyme in human tumors is associated with high invasiveness, enhanced metastasis, and poor prognosis. Thus, PPM1F is a target for pharmacological intervention, yet inhibitors of this enzyme are lacking. Here, we use high-throughput screening to identify Lockdown, a reversible and non-competitive PPM1F inhibitor. Lockdown is selective for PPM1F, because this compound does not inhibit other protein phosphatases in vitro and does not induce additional phenotypes in PPM1F knockout cells. Importantly, Lockdown-treated glioblastoma cells fully re-capitulate the phenotype of PPM1F-deficient cells as assessed by increased phosphorylation of PPM1F substrates and corruption of integrin-dependent cellular processes. Ester modification yields Lockdown^Pro with increased membrane permeability and prodrug-like properties. Lockdown^Pro suppresses tissue invasion by PPM1F-overexpressing human cancer cells, validating PPM1F as a therapeutic target and providing an access point to control tumor cell dissemination.

Benchmarking ensemble docking methods in D3R Grand Challenge 4

The discovery of new drugs is a time consuming and expensive process. Methods such as virtual screening, which can filter out ineffective compounds from drug libraries prior to expensive experimental study, have become popular research topics. As the computational drug discovery community has grown, in order to benchmark the various advances in methodology, organizations such as the Drug Design Data Resource have begun hosting blinded grand challenges seeking to identify the best methods for ligand pose-prediction, ligand affinity ranking, and free energy calculations. Such open challenges offer a unique opportunity for researchers to partner with junior students (e.g., high school and undergraduate) to validate basic yet fundamental hypotheses considered to be uninteresting to domain experts. Here, we, a group of high school-aged students and their mentors, present the results of our participation in Grand Challenge 4 where we predicted ligand affinity rankings for the Cathepsin S protease, an important protein target for autoimmune diseases. To investigate the effect of incorporating receptor dynamics on ligand affinity rankings, we employed the Relaxed Complex Scheme, a molecular docking method paired with molecular dynamics-generated receptor conformations. We found that Cathepsin S is a difficult target for molecular docking and we explore some advanced methods such as distance-restrained docking to try to improve the correlation with experiments. This project has exemplified the capabilities of high school students when supported with a rigorous curriculum, and demonstrates the value of community-driven competitions for beginners in computational drug discovery.

A practical guide to large-scale docking

Structure-based docking screens of large compound libraries have become common in early drug and probe discovery. As computer efficiency has improved and compound libraries have grown, the ability to screen hundreds of millions, and even billions, of compounds has become feasible for modest-sized computer clusters. This allows the rapid and cost-effective exploration and categorization of vast chemical space into a subset enriched with potential hits for a given target. To accomplish this goal at speed, approximations are used that result in undersampling of possible configurations and inaccurate predictions of absolute binding energies. Accordingly, it is important to establish controls, as are common in other fields, to enhance the likelihood of success in spite of these challenges. Here we outline best practices and control docking calculations that help evaluate docking parameters for a given target prior to undertaking a large-scale prospective screen, with exemplification in one particular target, the melatonin receptor, where following this procedure led to direct docking hits with activities in the subnanomolar range. Additional controls are suggested to ensure specific activity for experimentally validated hit compounds. These guidelines should be useful regardless of the docking software used. Docking software described in the outlined protocol (DOCK3.7) is made freely available for academic research to explore new hits for a range of targets.

CHARMM-GUI LBS Finder & Refiner for Ligand Binding Site Prediction and Refinement

A protein performs its task by binding a variety of ligands in its local region that is also known as the ligand-binding-site (LBS). Therefore, accurate prediction, characterization, and refinement of LBS can facilitate protein functional annotations and structure-based drug design. In this work, we present CHARMM-GUI LBS Finder & Refiner (https://www.charmm-gui.org/input/lbsfinder) that predicts potential LBS, offers interactive features for local LBS structure analysis, and prepares various molecular dynamics (MD) systems and inputs by setting up distance restraint potentials for LBS structure refinement. LBS Finder & Refiner supports 5 different commonly used simulation programs, such as NAMD, AMBER, GROMACS, GENESIS, and OpenMM, for LBS structure refinement together with hydrogen mass repartitioning. The capability of LBS Finder & Refiner is illustrated through LBS structure predictions and refinements of 48 modeled and 20 apo benchmark target proteins. Overall, successful LBS structure predictions and refinements are seen in our benchmark tests. We hope that LBS Finder & Refiner is useful to predict, characterize, and refine potential LBS on any given protein of interest.

A Brief Overview of Potential Treatments for Viral Diseases Using Natural Plant Compounds: The Case of SARS-Cov

The COVID-19 pandemic, as well as the more general global increase in viral diseases, has led researchers to look to the plant kingdom as a potential source for antiviral compounds. Since ancient times, herbal medicines have been extensively applied in the treatment and prevention of various infectious diseases in different traditional systems. The purpose of this review is to highlight the potential antiviral activity of plant compounds as effective and reliable agents against viral infections, especially by viruses from the coronavirus group. Various antiviral mechanisms shown by crude plant extracts and plant-derived bioactive compounds are discussed. The understanding of the action mechanisms of complex plant extract and isolated plant-derived compounds will help pave the way towards the combat of this life-threatening disease. Further, molecular docking studies, in silico analyses of extracted compounds, and future prospects are included. The in vitro production of antiviral chemical compounds from plants using molecular pharming is also considered. Notably, hairy root cultures represent a promising and sustainable way to obtain a range of biologically active compounds that may be applied in the development of novel antiviral agents.

A High-Resolution View of the Coordination Environment in a Paramagnetic Metalloprotein from its Magnetic Properties

Metalloproteins constitute a significant fraction of the proteome of all organisms and their characterization is critical for both basic sciences and biomedical applications. A large portion of metalloproteins bind paramagnetic metal ions, and paramagnetic NMR spectroscopy has been widely used in their structural characterization. However, the signals of nuclei in the immediate vicinity of the metal center are often broadened beyond detection. In this work, we show that it is possible to determine the coordination environment of the paramagnetic metal in the protein at a resolution inaccessible to other techniques. Taking the structure of a diamagnetic analogue as a starting point, a geometry optimization is carried out by fitting the pseudocontact shifts obtained from first principles quantum chemical calculations to the experimental ones.

Impact of Traditional Plants and their Secondary Metabolites in the Discovery of COVID-19 Treatment

BACKGROUND

Coronavirus Disease-2019 belongs to the family of viruses which cause serious pneumonia along with fever, breathing issues and infection of lungs, and was first reported in China and later spread worldwide.

OBJECTIVE

Several studies and clinical trials have been conducted to identify potential drugs and vaccines for Coronavirus Disease-2019. The present study listed natural secondary metabolites identified from plant sources with antiviral properties and could be a safer and tolerable treatment for Coronavirus Disease-2019.

METHODS

A comprehensive search on the reported studies was conducted using different search engines such as Google Scholar, SciFinder, Sciencedirect, Medline PubMed, and Scopus for the collection of research articles based on plant-derived secondary metabolites, herbal extracts, and traditional medicine for coronavirus infections.

RESULTS

Status of COVID-19 worldwide and information of important molecular targets involved in COVID- 19 are described, and through literature search, it is highlighted that numerous plant species and their extracts possess antiviral properties and are studied with respect to coronavirus treatments. Chemical information, plant source, test system type with a mechanism of action for each secondary metabolite are also mentioned in this review paper.

CONCLUSION

The present review has listed plants that have presented antiviral potential in the previous coronavirus pandemics and their secondary metabolites, which could be significant for the development of novel and a safer drug which could prevent and cure coronavirus infection worldwide.

Identification of Potential Binders of Mtb Universal Stress Protein (Rv1636) Through an in silico Approach and Insights Into Compound Selection for Experimental Validation

Millions of deaths caused by Mycobacterium tuberculosis (Mtb) are reported worldwide every year. Treatment of tuberculosis (TB) involves the use of multiple antibiotics over a prolonged period. However, the emergence of resistance leading to multidrug-resistant TB (MDR-TB) and extensively drug-resistant TB (XDR-TB) is the most challenging aspect of TB treatment. Therefore, there is a constant need to search for novel therapeutic strategies that could tackle the growing problem of drug resistance. One such strategy could be perturbing the functions of novel targets in Mtb, such as universal stress protein (USP, Rv1636), which binds to cAMP with a higher affinity than ATP. Orthologs of these proteins are conserved in all mycobacteria and act as “sink” for cAMP, facilitating the availability of this second messenger for signaling when required. Here, we have used the cAMP-bound crystal structure of USP from Mycobacterium smegmatis, a closely related homolog of Mtb, to conduct a structure-guided hunt for potential binders of Rv1636, primarily employing molecular docking approach. A library of 1.9 million compounds was subjected to virtual screening to obtain an initial set of ~2,000 hits. An integrative strategy that uses the available experimental data and consensus indications from other computational analyses has been employed to prioritize 22 potential binders of Rv1636 for experimental validations. Binding affinities of a few compounds among the 22 prioritized compounds were tested through microscale thermophoresis assays, and two compounds of natural origin showed promising binding affinities with Rv1636. We believe that this study provides an important initial guidance to medicinal chemists and biochemists to synthesize and test an enriched set of compounds that have the potential to inhibit Mtb USP (Rv1636), thereby aiding the development of novel antitubercular lead candidates.

De novo 3D models of SARS-CoV-2 RNA elements from consensus experimental secondary structures

The rapid spread of COVID-19 is motivating development of antivirals targeting conserved SARS-CoV-2 molecular machinery. The SARS-CoV-2 genome includes conserved RNA elements that offer potential small-molecule drug targets, but most of their 3D structures have not been experimentally characterized. Here, we provide a compilation of chemical mapping data from our and other labs, secondary structure models, and 3D model ensembles based on Rosetta's FARFAR2 algorithm for SARS-CoV-2 RNA regions including the individual stems SL1-8 in the extended 5' UTR; the reverse complement of the 5' UTR SL1-4; the frameshift stimulating element (FSE); and the extended pseudoknot, hypervariable region, and s2m of the 3' UTR. For eleven of these elements (the stems in SL1-8, reverse complement of SL1-4, FSE, s2m and 3' UTR pseudoknot), modeling convergence supports the accuracy of predicted low energy states; subsequent cryo-EM characterization of the FSE confirms modeling accuracy. To aid efforts to discover small molecule RNA binders guided by computational models, we provide a second set of similarly prepared models for RNA riboswitches that bind small molecules. Both datasets ('FARFAR2-SARS-CoV-2', https://github.com/DasLab/FARFAR2-SARS-CoV-2; and 'FARFAR2-Apo-Riboswitch', at https://github.com/DasLab/FARFAR2-Apo-Riboswitch') include up to 400 models for each RNA element, which may facilitate drug discovery approaches targeting dynamic ensembles of RNA molecules.

Ligand-Binding-Site Refinement to Generate Reliable Holo Protein Structure Conformations from Apo Structures

The first important step in a structure-based virtual screening is the judicious selection of a receptor protein. In cases where the holo protein receptor structure is unavailable, significant reduction in virtual screening performance has been reported. In this work, we present a robust method to generate reliable holo protein structure conformations from apo structures using molecular dynamics (MD) simulation with restraints derived from holo structure binding-site templates. We perform benchmark tests on two different datasets: 40 structures from a directory of useful decoy-enhanced (DUD-E) and 84 structures from the Gunasekaran dataset. Our results show successful refinement of apo binding-site structures toward holo conformations in 82% of the test cases. In addition, virtual screening performance of 40 DUD-E structures is significantly improved using our MD-refined structures as receptors with an average enrichment factor (EF), an EF_1% value of 6.2 compared to apo structures with 3.5. Docking of native ligands to the refined structures shows an average ligand root mean square deviation (RMSD) of 1.97 Å (DUD-E dataset and Gunasekaran dataset) relative to ligands in the holo crystal structures, which is comparable to the self-docking (i.e., docking of the native ligand back to its crystal structure receptor) average, 1.34 Å (DUD-E dataset) and 1.36 Å (Gunasekaran dataset). On the other hand, docking to the apo structures yields an average ligand RMSD of 3.65 Å (DUD-E) and 2.90 Å (Gunasekaran). These results indicate that our method is robust and can be useful to improve virtual screening performance of apo structures.

Alchemical Binding Free Energy Calculations in AMBER20: Advances and Best Practices for Drug Discovery

Predicting protein-ligand binding affinities and the associated thermodynamics of biomolecular recognition is a primary objective of structure-based drug design. Alchemical free energy simulations offer a highly accurate and computationally efficient route to achieving this goal. While the AMBER molecular dynamics package has successfully been used for alchemical free energy simulations in academic research groups for decades, widespread impact in industrial drug discovery settings has been minimal because of the previous limitations within the AMBER alchemical code, coupled with challenges in system setup and postprocessing workflows. Through a close academia-industry collaboration we have addressed many of the previous limitations with an aim to improve accuracy, efficiency, and robustness of alchemical binding free energy simulations in industrial drug discovery applications. Here, we highlight some of the recent advances in AMBER20 with a focus on alchemical binding free energy (BFE) calculations, which are less computationally intensive than alternative binding free energy methods where full binding/unbinding paths are explored. In addition to scientific and technical advances in AMBER20, we also describe the essential practical aspects associated with running relative alchemical BFE calculations, along with recommendations for best practices, highlighting the importance not only of the alchemical simulation code but also the auxiliary functionalities and expertise required to obtain accurate and reliable results. This work is intended to provide a contemporary overview of the scientific, technical, and practical issues associated with running relative BFE simulations in AMBER20, with a focus on real-world drug discovery applications.

Combining fragment docking with graph theory to improve ligand docking for homology model structures

Computational protein–ligand docking is well-known to be prone to inaccuracies in input receptor structures, and it is challenging to obtain good docking results with computationally predicted receptor structures (e.g. through homology modeling). Here we introduce a fragment-based docking method and test if it reduces requirements on the accuracy of an input receptor structures relative to non-fragment docking approaches. In this method, small rigid fragments are docked first using AutoDock Vina to generate a large number of favorably docked poses spanning the receptor binding pocket. Then a graph theory maximum clique algorithm is applied to find combined sets of docked poses of different fragment types onto which the complete ligand can be properly aligned. On the basis of these alignments, possible binding poses of complete ligand are determined. This docking method is first tested for bound docking on a series of Cytochrome P450 (CYP450) enzyme–substrate complexes, in which experimentally determined receptor structures are used. For all complexes tested, ligand poses of less than 1 Å root mean square deviations (RMSD) from the actual binding positions can be recovered. Then the method is tested for unbound docking with modeled receptor structures for a number of protein–ligand complexes from different families including the very recent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) protease. For all complexes, poses with RMSD less than 3 Å from actual binding positions can be recovered. Our results suggest that for docking with approximately modeled receptor structures, fragment-based methods can be more effective than common complete ligand docking approaches.

Benchmarking Data Sets from PubChem BioAssay Data: Current Scenario and Room for Improvement

Developing realistic data sets for evaluating virtual screening methods is a task that has been tackled by the cheminformatics community for many years. Numerous artificially constructed data collections were developed, such as DUD, DUD-E, or DEKOIS. However, they all suffer from multiple drawbacks, one of which is the absence of experimental results confirming the impotence of presumably inactive molecules, leading to possible false negatives in the ligand sets. In light of this problem, the PubChem BioAssay database, an open-access repository providing the bioactivity information of compounds that were already tested on a biological target, is now a recommended source for data set construction. Nevertheless, there exist several issues with the use of such data that need to be properly addressed. In this article, an overview of benchmarking data collections built upon experimental PubChem BioAssay input is provided, along with a thorough discussion of noteworthy issues that one must consider during the design of new ligand sets from this database. The points raised in this review are expected to guide future developments in this regard, in hopes of offering better evaluation tools for novel in silico screening procedures.

Fragment Binding Pose Predictions Using Unbiased Simulations and Markov-State Models

Predicting the costructure of small-molecule ligands and their respective target proteins has been a long-standing problem in drug discovery. For weak binding compounds typically identified in fragment-based screening (FBS) campaigns, determination of the correct binding site and correct binding mode is usually done experimentally via X-ray crystallography. For many targets of pharmaceutical interest, however, establishing an X-ray system which allows for sufficient throughput to support a drug discovery project is not possible. In this case, exploration of fragment hits becomes a very laborious and consequently slow process with the generation of protein/ligand cocrystal structures as the bottleneck of the entire process. In this work, we introduce a computational method which is able to reliably predict binding sites and binding modes of fragment-like small molecules using solely the structure of the apoprotein and the ligand's chemical structure as input information. The method is based on molecular dynamics simulations and Markov-state models and can be run as a fully automated protocol requiring minimal human intervention. We describe the application of the method to a representative subset of different target classes and fragments from historical FBS efforts at Boehringer Ingelheim and discuss its potential integration into the overall fragment-based drug discovery workflow.

Modeling the binding of diverse ligands within the Ah receptor ligand binding domain

The Ah receptor (AhR) is a ligand-dependent transcription factor belonging to the basic helix-loop-helix Per-Arnt-Sim (bHLH-PAS) superfamily. Binding to and activation of the AhR by a variety of chemicals results in the induction of expression of diverse genes and production of a broad spectrum of biological and toxic effects. The AhR also plays important roles in several physiological responses, which has led it to become a novel target for the development of therapeutic drugs. Differences in the interactions of various ligands within the AhR ligand binding domain (LBD) may contribute to differential modulation of AhR functionality. We combined computational and experimental analyses to investigate the binding modes of a group of chemicals representative of major classes of AhR ligands. On the basis of a novel computational approach for molecular docking to the homology model of the AhR LBD that includes the receptor flexibility, we predicted specific residues within the AhR binding cavity that play a critical role in binding of three distinct groups of chemicals. The prediction was validated by site-directed mutagenesis and evaluation of the relative ligand binding affinities for the mutant AhRs. These results provide an avenue for understanding ligand modulation of the AhR functionality and for rational drug design.

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.

In Need of Bias Control: Evaluating Chemical Data for Machine Learning in Structure-Based Virtual Screening

Reports of successful applications of machine learning (ML) methods in structure-based virtual screening (SBVS) are increasing. ML methods such as convolutional neural networks show promising results and often outperform traditional methods such as empirical scoring functions in retrospective validation. However, trained ML models are often treated as black boxes and are not straightforwardly interpretable. In most cases, it is unknown which features in the data are decisive and whether a model's predictions are right for the right reason. Hence, we re-evaluated three widely used benchmark data sets in the context of ML methods and came to the conclusion that not every benchmark data set is suitable. Moreover, we demonstrate on two examples from current literature that bias is learned implicitly and unnoticed from standard benchmarks. On the basis of these results, we conclude that there is a need for eligible validation experiments and benchmark data sets suited to ML for more bias-controlled validation in ML-based SBVS. Therefore, we provide guidelines for setting up validation experiments and give a perspective on how new data sets could be generated.

Probing the Electrostatic and Steric Requirements for Substrate Binding in Human Platelet-Type 12-Lipoxygenase

Human platelet ALOX12 (hALOX12 or h12-LOX) has been implicated in a variety of human diseases. The present study investigates the active site of hALOX12 to more thoroughly understand how it positions the substrate and achieves nearly perfect regio- and stereospecificities (i.e., 100 ± 5% of the 12(S)-hydroperoxide product), utilizing site-directed mutagenesis. Specifically, we have determined that Arg402 is not as important in substrate binding as previously seen for hALOX15 but that His596 may play a role in anchoring the carboxy terminal of the arachidonic acid during catalysis. In addition, Phe414 creates a π-stacking interaction with a double bond of arachidonic acid (Δ¹¹), and Ala417/Val418 define the bottom of the cavity. However, the influence of Ala417/Val418 on the profile is markedly less for hALOX12 than that seen in hALOX15. Mutating these two residues to larger amino acids (Ala417Ile/Val418Met) only increased the generation of 15-HpETE by 24 ± 2%, but conversely, smaller residues at these positions converted hALOX15 to almost 100% hALOX12 reactivity [Gan et al. (1996) J. Biol. Chem. 271, 25412-25418]. However, we were able to increase 15-HpETE to 46 ± 3% by restricting the width of the active site with the Ala417Ile/Val418Met/Ser594Thr mutation, indicating both depth and width of the active site are important. Finally, residue Leu407 is shown to play a critical role in positioning the substrate correctly, as seen by the increase of 15-HpETE to 21 ± 1% for the single Leu407Gly mutant. These results outline critical differences between the active site requirements of hALOX12 relative to hALOX15 and explain both their product specificity and inhibitory differences.

Binding Analysis of the Inositol-Requiring Enzyme 1 Kinase Domain

Inositol-requiring enzyme 1 (IRE1) is an orchestrator of the unfolded protein response (UPR), the cellular response to endoplasmic reticulum (ER) stress that plays a crucial role in tumor development. IRE1 signaling is the most evolutionary conserved branch of the UPR. Under ER stress, the IRE1 luminal domain undergoes a conformational change to multimerize, resulting in trans-autophosphorylation and activation of the cytosolic kinase and endoribonuclease domain. Adenosine triphosphate-competitive inhibitors that bind to the IRE1 kinase site can modulate the activity of the RNase domain through an allosteric relationship between the IRE1 kinase and RNase domains. The current study aims at the investigation of available structural data of the IRE1 kinase domain and provides insights into the design of novel kinase inhibitors. To this end, a detailed analysis of IRE1 kinase active site and investigation of suitable structures for virtual screening studies were performed. The results indicate in silico target fishing as an appropriate strategy for the identification of novel IRE1 kinase binders, further validating the robustness of the in silico protocol. Importantly, the study highlights the kinase-inhibiting RNase attenuator (KIRA)-bound protein data bank 4U6R structure as the best protein structure to perform virtual screening to develop diverse and more potent KIRA-like IRE1 kinase inhibitors that are capable of allosterically affecting the RNase activity.

Computational Strategy for Bound State Structure Prediction in Structure-Based Virtual Screening: A Case Study of Protein Tyrosine Phosphatase Receptor Type O Inhibitors

Accurate protein structure in the ligand-bound state is a prerequisite for successful structure-based virtual screening (SBVS). Therefore, applications of SBVS against targets for which only an apo structure is available may be severely limited. To address this constraint, we developed a computational strategy to explore the ligand-bound state of a target protein, by combined use of molecular dynamics simulation, MM/GBSA binding energy calculation, and fragment-centric topographical mapping. Our computational strategy is validated against low-molecular weight protein tyrosine phosphatase (LMW-PTP) and then successfully employed in the SBVS against protein tyrosine phosphatase receptor type O (PTPRO), a potential therapeutic target for various diseases. The most potent hit compound GP03 showed an IC₅₀ value of 2.89 μM for PTPRO and possessed a certain degree of selectivity toward other protein phosphatases. Importantly, we also found that neglecting the ligand energy penalty upon binding partially accounts for the false positive SBVS hits. The preliminary structure-activity relationships of GP03 analogs are also reported.

Exploring the conformational landscapes of HIV protease structural ensembles using principal component analysis

HIV protease, an essential enzyme for viral particle maturation, is an important drug target of HIV. Its structural conformation is a key determinant of both biological function as well as efficient binding of protease inhibitor molecules. In the present study we analyzed 471 crystal structures of HIV-1 protease to understand the conformational changes induced by mutations or binding of various ligands and substrates. We performed principal component analysis on the ensembles of the HIV-1 protease structures to explore the conformational landscapes. The study identified structural differences between drug resistant and drug sensitive protease structures. Conformational changes were identified in the A and B chains of homo-dimeric HIV protease structures having different combinations of mutations, and also rigidity in the binding conformation of HIV drugs within the active site of the protein.© 2018 Wiley Periodicals, Inc.

Online structure-based screening of purchasable approved drugs and natural compounds: retrospective examples of drug repositioning on cancer targets

Drug discovery is a long and difficult process that benefits from the integration of virtual screening methods in experimental screening campaigns such as to generate testable hypotheses, accelerate and/or reduce the cost of drug development. Current drug attrition rate is still a major issue in all therapeutic areas and especially in the field of cancer. Drug repositioning as well as the screening of natural compounds constitute promising approaches to accelerate and improve the success rate of drug discovery. We developed three compounds libraries of purchasable compounds: Drugs-lib, FOOD-lib and NP-lib that contain approved drugs, food constituents and natural products, respectively, that are optimized for structure-based virtual screening studies. The three compounds libraries are implemented in the MTiOpenScreen web server that allows users to perform structure-based virtual screening computations on their selected protein targets. The server outputs a list of 1,500 molecules with predicted binding scores that can then be processed further by the users and purchased for experimental validation. To illustrate the potential of our service for drug repositioning endeavours, we selected five recently published drugs that have been repositioned in vitro and/or in vivo on cancer targets. For each drug, we used the MTiOpenScreen service to screen the Drugs-lib collection against the corresponding anti-cancer target and we show that our protocol is able to rank these drugs within the top ranked compounds. This web server should assist the discovery of promising molecules that could benefit patients, with faster development times, and reduced costs and risk.