Dynophore-Based Approach in Virtual Screening: A Case of Human DNA Topoisomerase IIα

In this study, we utilized human DNA topoisomerase IIα as a model target to outline a dynophore-based approach to catalytic inhibitor design. Based on MD simulations of a known catalytic inhibitor and the native ATP ligand analog, AMP-PNP, we derived a joint dynophore model that supplements the static structure-based-pharmacophore information with a dynamic component. Subsequently, derived pharmacophore models were employed in a virtual screening campaign of a library of natural compounds. Experimental evaluation identified flavonoid compounds with promising topoisomerase IIα catalytic inhibition and binding studies confirmed interaction with the ATPase domain. We constructed a binding model through docking and extensively investigated it with molecular dynamics MD simulations, essential dynamics, and MM-GBSA free energy calculations, thus reconnecting the new results to the initial dynophore-based screening model. We not only demonstrate a new design strategy that incorporates a dynamic component of molecular recognition, but also highlight new derivates in the established flavonoid class of topoisomerase II inhibitors.     

Identification of New KRAS G12D Inhibitors through Computer-Aided Drug Discovery Methods

Owing to several mutations, the oncogene Kirsten rat sarcoma 2 viral oncogene homolog (KRAS) is activated in the majority of cancers, and targeting it has been pharmacologically challenging. In this study, using an in silico approach comprised of pharmacophore modeling, molecular docking, and molecular dynamics simulations, potential KRAS G12D inhibitors were investigated. A ligand-based common feature pharmacophore model was generated to identify the framework necessary for effective KRAS inhibition. The chemical features in the selected pharmacophore model comprised two hydrogen bond donors, one hydrogen bond acceptor, two aromatic rings and one hydrophobic feature. This model was used for screening in excess of 214,000 compounds from InterBioScreen (IBS) and ZINC databases. Eighteen compounds from the IBS and ten from the ZINC database mapped onto the pharmacophore model and were subjected to molecular docking. Molecular docking results highlighted a higher affinity of four hit compounds towards KRAS G12D in comparison to the reference inhibitor, BI-2852. Sequential molecular dynamics (MD) simulation studies revealed all four hit compounds them possess higher KRAS G12D binding free energy and demonstrate stable polar interaction with key residues. Further, Principal Component Analysis (PCA) analysis of the hit compounds in complex with KRAS G12D also indicated stability. Overall, the research undertaken provides strong support for further in vitro testing of these newly identified KRAS G12D inhibitors, particularly Hit1 and Hit2.     

Identification of Novel Potential β-N-Acetyl-D-Hexosaminidase Inhibitors by Virtual Screening, Molecular Dynamics Simulation and MM-PBSA Calculations

Chitinolytic β-N-acetyl-d-hexosaminidases, as a class of chitin hydrolysis enzyme in insects, are a potential species-specific target for developing environmentally-friendly pesticides. Until now, pesticides targeting chitinolytic β-N-acetyl-d-hexosaminidase have not been developed. This study demonstrates a combination of different theoretical methods for investigating the key structural features of this enzyme responsible for pesticide inhibition, thus allowing for the discovery of novel small molecule inhibitors. Firstly, based on the currently reported crystal structure of this protein (OfHex1.pdb), we conducted a pre-screening of a drug-like compound database with 8 × 10(6) compounds by using the expanded pesticide-likeness criteria, followed by docking-based screening, obtaining 5 top-ranked compounds with favorable docking conformation into OfHex1. Secondly, molecular docking and molecular dynamics simulations are performed for the five complexes and demonstrate that one main hydrophobic pocket formed by residues Trp424, Trp448 and Trp524, which is significant for stabilization of the ligand-receptor complex, and key residues Asp477 and Trp490, are respectively responsible for forming hydrogen-bonding and π-π stacking interactions with the ligands. Finally, the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) analysis indicates that van der Waals interactions are the main driving force for the inhibitor binding that agrees with the fact that the binding pocket of OfHex1 is mainly composed of hydrophobic residues. These results suggest that screening the ZINC database can maximize the identification of potential OfHex1 inhibitors and the computational protocol will be valuable for screening potential inhibitors of the binding mode, which is useful for the future rational design of novel, potent OfHex1-specific pesticides.     

Discovery of Tyrosine Kinase Inhibitors by Docking into an Inactive Kinase Conformation Generated by Molecular Dynamics

Several small molecules that bind to the inactive DFG‐out conformation of tyrosine kinases (called type II inhibitors) have shown a good selectivity profile over other kinase targets. To obtain a set of DFG‐out structures, we performed an explicit solvent molecular dynamics (MD) simulation of the complex of the catalytic domain of a tyrosine kinase receptor, ephrin type‐A receptor 3 (EphA3), and a manually docked type II inhibitor. Automatic docking of four previously reported type II inhibitors was used to select a single snapshot from the MD trajectory for virtual screening. High‐throughput docking of a pharmacophore‐tailored library of 175 000 molecules resulted in about 4 million poses, which were further filtered by van der Waals efficiency and ranked according to a force‐field‐based energy function. Notably, around 20 % of the compounds with predicted binding energy smaller than ?10 kcal mol^?1 are known type II inhibitors. Moreover, a series of 5‐(piperazine‐1‐yl)isoquinoline derivatives was identified as a novel class of low‐micromolar inhibitors of EphA3 and unphosphorylated Abelson tyrosine kinase (Abl1). The in silico predicted binding mode of the new inhibitors suggested a similar affinity to the gatekeeper mutant T315I of Abl1, which was verified in vitro by using a competition binding assay. Additional evidence for the type II binding mode was obtained by two 300 ns MD simulations of the complex between N‐(3‐chloro‐4‐(difluoromethoxy)phenyl)‐2‐(4‐(8‐nitroisoquinolin‐5‐yl)piperazin‐1‐yl)acetamide and EphA3.     

3D QSAR pharmacophore modeling, in silico screening, and density functional theory (DFT) approaches for identification of human chymase inhibitors

Human chymase is a very important target for the treatment of cardiovascular diseases. Using a series of theoretical methods like pharmacophore modeling, database screening, molecular docking and Density Functional Theory (DFT) calculations, an investigation for identification of novel chymase inhibitors, and to specify the key factors crucial for the binding and interaction between chymase and inhibitors is performed. A highly correlating (r = 0.942) pharmacophore model (Hypo1) with two hydrogen bond acceptors, and three hydrophobic aromatic features is generated. After successfully validating "Hypo1", it is further applied in database screening. Hit compounds are subjected to various drug-like filtrations and molecular docking studies. Finally, three structurally diverse compounds with high GOLD fitness scores and interactions with key active site amino acids are identified as potent chymase hits. Moreover, DFT study is performed which confirms very clear trends between electronic properties and inhibitory activity (IC(50)) data thus successfully validating "Hypo1" by DFT method. Therefore, this research exertion can be helpful in the development of new potent hits for chymase. In addition, the combinational use of docking, orbital energies and molecular electrostatic potential analysis is also demonstrated as a good endeavor to gain an insight into the interaction between chymase and inhibitors.     

In Silico Evaluation of the Binding Site of Fucosyltransferase 8 and First Attempts to Synthesize an Inhibitor with Drug-Like Properties

Core fucosylation of N-glycans is catalyzed by fucosyltransferase 8 and is associated with various types of cancer. Most reported fucosyltransferase inhibitors contain non-drug-like features, such as charged groups. New starting points for the development of inhibitors of fucosyltransferase 8 using a fragment-based strategy are presented. Firstly, we discuss the potential of a new putative binding site of fucosyltransferase 8 that, according to a molecular dynamics (MD) simulation, is made accessible by a significant motion of the SH3 domain. This might enable the design of completely new inhibitor types for fucosyltransferase 8. Secondly, we have performed a docking study targeting the donor binding site of fucosyltransferase 8, and this yielded two fragments that were linked and trimmed in silico. The resulting ligand was synthesized. Saturation transfer difference (STD) NMR confirmed binding of the ligand featuring a pyrazole core that mimics the guanine moiety. This ligand represents the first low-molecular-weight compound for the development of inhibitors of fucosyltransferase 8 with drug-like properties.     

Computational Investigation Identified Potential Chemical Scaffolds for Heparanase as Anticancer Therapeutics

Heparanase (Hpse) is an endo-β-D-glucuronidase capable of cleaving heparan sulfate side chains. Its upregulated expression is implicated in tumor growth, metastasis and angiogenesis, thus making it an attractive target in cancer therapeutics. Currently, a few small molecule inhibitors have been reported to inhibit Hpse, with promising oral administration and pharmacokinetic (PK) properties. In the present study, a ligand-based pharmacophore model was generated from a dataset of well-known active small molecule Hpse inhibitors which were observed to display favorable PK properties. The compounds from the InterBioScreen database of natural (69,034) and synthetic (195,469) molecules were first filtered for their drug-likeness and the pharmacophore model was used to screen the drug-like database. The compounds acquired from screening were subjected to molecular docking with Heparanase, where two molecules used in pharmacophore generation were used as reference. From the docking analysis, 33 compounds displayed higher docking scores than the reference and favorable interactions with the catalytic residues. Complex interactions were further evaluated by molecular dynamics simulations to assess their stability over a period of 50 ns. Furthermore, the binding free energies of the 33 compounds revealed 2 natural and 2 synthetic compounds, with better binding affinities than reference molecules, and were, therefore, deemed as hits. The hit compounds presented from this in silico investigation could act as potent Heparanase inhibitors and further serve as lead scaffolds to develop compounds targeting Heparanase upregulation in cancer.     

Novel Hybrid Virtual Screening Protocol Based on Molecular Docking and Structure-Based Pharmacophore for Discovery of Methionyl-tRNA Synthetase Inhibitors as Antibacterial Agents

Methione tRNA synthetase (MetRS) is an essential enzyme involved in protein biosynthesis in all living organisms and is a potential antibacterial target. In the current study, the structure-based pharmacophore (SBP)-guided method has been suggested to generate a comprehensive pharmacophore of MetRS based on fourteen crystal structures of MetRS-inhibitor complexes. In this investigation, a hybrid protocol of a virtual screening method, comprised of pharmacophore model-based virtual screening (PBVS), rigid and flexible docking-based virtual screenings (DBVS), is used for retrieving new MetRS inhibitors from commercially available chemical databases. This hybrid virtual screening approach was then applied to screen the Specs (202,408 compounds) database, a structurally diverse chemical database. Fifteen hit compounds were selected from the final hits and shifted to experimental studies. These results may provide important information for further research of novel MetRS inhibitors as antibacterial agents.     

Targeting the SARS-CoV-2 HR1 with Small Molecules as Inhibitors of the Fusion Process

The rapid and global propagation of the novel human coronavirus that causes severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has produced an immediate urgency to discover promising targets for the treatment of this virus. In this paper, we studied the spike protein S2 domain of SARS-CoV-2 as it is the most conserved component and controls the crucial fusion process of SARS-CoV-2 as a target for different databases of small organic compounds. Our in silico methodology, based on pharmacophore modeling, docking simulation and molecular dynamics simulations, was first validated with ADS-J1, a potent small-molecule HIV fusion inhibitor that has already proved effective in binding the HR1 domain and inhibiting the fusion core of SARS-CoV-1. It then focused on finding novel small molecules and new peptides as fusion inhibitors. Our methodology identified several small molecules and peptides as potential inhibitors of the fusion process. Among these, NF 023 hydrate (MolPort-006-822-583) is one of the best-scored compounds. Other compounds of interest are ZINC00097961973, Salvianolic acid, Thalassiolin A and marine_160925_88_2. Two interesting active peptides were also identified: AP00094 (Temporin A) and AVP1227 (GBVA5). The inhibition of the spike protein of SARS-CoV-2 is a valid target to inhibit the virus entry in human cells. The discussed compounds reported in this paper led to encouraging results for future in vitro tests against SARS-CoV-2.     

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.     

Design and Synthesis of Dihydroxamic Acids as HDAC6/8/10 Inhibitors

We report the synthesis and evaluation of a class of selective multitarget agents for the inhibition of HDAC6, HDAC8, and HDAC10. The concept for this study grew out of a structural analysis of the two selective inhibitors Tubastatin A (HDAC6/10) and PCI-34051 (HDAC8), which we recognized share the same N-benzylindole core. Hybridization of the two inhibitor structures resulted in dihydroxamic acids with benzyl-indole and -indazole core motifs. These substances exhibit potent activity against HDAC6, HDAC8, and HDAC10, while retaining selectivity over HDAC1, HDAC2, and HDAC3. The best substance inhibited the viability of the SK-N-BE(2)C neuroblastoma cell line with an IC₅₀ value similar to a combination treatment with Tubastatin A and PCI-34051. This compound class establishes a proof of concept for such hybrid molecules and could serve as a starting point for the further development of enhanced HDAC6/8/10 inhibitors.     

Structure-Based Design,Docking and Binding Free Energy Calculations of A366 Derivatives as Spindlin1 Inhibitors

The chromatin reader protein Spindlin1 plays an important role in epigenetic regulation, through which it has been linked to several types of malignant tumors. In the current work, we report on the development of novel analogs of the previously published lead inhibitor A366. In an effort to improve the activity and explore the structure–activity relationship (SAR), a series of 21 derivatives was synthesized, tested in vitro, and investigated by means of molecular modeling tools. Docking studies and molecular dynamics (MD) simulations were performed to analyze and rationalize the structural differences responsible for the Spindlin1 activity. The analysis of MD simulations shed light on the important interactions. Our study highlighted the main structural features that are required for Spindlin1 inhibitory activity, which include a positively charged pyrrolidine moiety embedded into the aromatic cage connected via a propyloxy linker to the 2-aminoindole core. Of the latter, the amidine group anchor the compounds into the pocket through salt bridge interactions with Asp184. Different protocols were tested to identify a fast in silico method that could help to discriminate between active and inactive compounds within the A366 series. Rescoring the docking poses with MM-GBSA calculations was successful in this regard. Because A366 is known to be a G9a inhibitor, the most active developed Spindlin1 inhibitors were also tested over G9a and GLP to verify the selectivity profile of the A366 analogs. This resulted in the discovery of diverse selective compounds, among which 1s and 1t showed Spindlin1 activity in the nanomolar range and selectivity over G9a and GLP. Finally, future design hypotheses were suggested based on our findings.     

Potent and Selective Non‐hydroxamate Histone Deacetylase 8 Inhibitors

Specific inhibition of histone deacetylase 8 (HDAC8) has been suggested as a promising option for the treatment of neuroblastoma and T‐cell malignancies. A novel class of highly potent and selective HDAC8 inhibitors with a pyrimido[1,2‐c][1,3]benzothiazin‐6‐imine scaffold was studied that is completely different from the traditional concept of HDAC inhibitors comprising a zinc binding group (ZBG), in most cases a hydroxamate group, a spacer, and a capping group that may interact with the surface of the target protein. Although lacking a ZBG, some of the new compounds were shown to have outstanding potency against HDAC8 in the single‐digit nanomolar range. The pyrimido[1,2‐c][1,3]benzothiazin‐6‐imines also inhibited the growth of solid and hematological tumor cells. The small size and beneficial physicochemical properties of the novel HDAC inhibitor class underline the high degree of drug likeness. This and the broad structure–activity relationship suggest great potential for the further development of compounds with the pyrimido[1,2‐c][1,3]benzothiazin‐6‐imine scaffold into innovative and highly effective therapeutic drugs against cancer.     

Molecular Dynamics Simulation of Tryptophan Hydroxylase-1: Binding Modes and Free Energy Analysis to Phenylalanine Derivative Inhibitors

Serotonin is a neurotransmitter that modulates many central and peripheral functions. Tryptophan hydroxylase-1 (TPH1) is a key enzyme of serotonin synthesis. In the current study, the interaction mechanism of phenylalanine derivative TPH1 inhibitors was investigated using molecular dynamics (MD) simulations, free energy calculations, free energy decomposition analysis and computational alanine scanning. The predicted binding free energies of these complexes are consistent with the experimental data. The analysis of the individual energy terms indicates that although the van der Waals and electrostatics interaction contributions are important in distinguishing the binding affinities of these inhibitors, the electrostatic contribution plays a more crucial role in that. Moreover, it is observed that different configurations of the naphthalene substituent could form different binding patterns with protein, yet lead to similar inhibitory potency. The combination of different molecular modeling techniques is an efficient way to interpret the interaction mechanism of inhibitors and our work could provide valuable information for the TPH1 inhibitor design in the future.     

Combined 3D-QSAR, Molecular Docking and Molecular Dynamics Study on Derivatives of Peptide Epoxyketone and Tyropeptin-Boronic Acid as Inhibitors Against the β5 Subunit of Human 20S Proteasome

An abnormal ubiquitin-proteasome is found in many human diseases, especially in cancer, and has received extensive attention as a promising therapeutic target in recent years. In this work, several in silico models have been built with two classes of proteasome inhibitors (PIs) by using 3D-QSAR, homology modeling, molecular docking and molecular dynamics (MD) simulations. The study resulted in two types of satisfactory 3D-QSAR models, i.e., the CoMFA model (Q(2) = 0.462, R(2) (pred) = 0.820) for epoxyketone inhibitors (EPK) and the CoMSIA model (Q(2) = 0.622, R(2) (pred) = 0.821) for tyropeptin-boronic acid derivatives (TBA). From the contour maps, some key structural factors responsible for the activity of these two series of PIs are revealed. For EPK inhibitors, the N-cap part should have higher electropositivity; a large substituent such as a benzene ring is favored at the C6-position. In terms of TBA inhibitors, hydrophobic substituents with a larger size anisole group are preferential at the C8-position; higher electropositive substituents like a naphthalene group at the C3-position can enhance the activity of the drug by providing hydrogen bond interaction with the protein target. Molecular docking disclosed that residues Thr60, Thr80, Gly106 and Ser189 play a pivotal role in maintaining the drug-target interactions, which are consistent with the contour maps. MD simulations further indicated that the binding modes of each conformation derived from docking is stable and in accord with the corresponding structure extracted from MD simulation overall. These results can offer useful theoretical references for designing more potent PIs.     

Exploration of N-Arylsulfonyl-indole-2-carboxamide Derivatives as Novel Fructose-1,6-bisphosphatase Inhibitors by Molecular Simulation

A series of N-arylsulfonyl-indole-2-carboxamide derivatives have been identified as potent fructose-1,6-bisphosphatase (FBPase) inhibitors (FBPIs) with excellent selectivity for the potential therapy of type II diabetes mellitus. To explore the structure–activity relationships (SARs) and the mechanisms of action of these FBPIs, a systematic computational study was performed in the present study, including three-dimensional quantitative structure–activity relationship (3D-QSAR) modeling, pharmacophore modeling, molecular dynamics (MD), and virtual screening. The constructed 3D-QSAR models exhibited good predictive ability with reasonable parameters using comparative molecular field analysis (q² = 0.709, R² = 0.979, r_pre² = 0.932) and comparative molecular similarity indices analysis (q² = 0.716, R² = 0.978, r_pre² = 0.890). Twelve hit compounds were obtained by virtual screening using the best pharmacophore model in combination with molecular dockings. Three compounds with relatively higher docking scores and better ADME properties were then selected for further studies by docking and MD analyses. The docking results revealed that the amino acid residues Met18, Gly21, Gly26, Leu30, and Thr31 at the binding site were of great importance for the effective bindings of these FBPIs. The MD results indicated that the screened compounds VS01 and VS02 could bind with FBPase stably as its cognate ligand in dynamic conditions. This work identified several potential FBPIs by modeling studies and might provide important insights into developing novel FBPIs.     

Valproic Acid as a Potential Inhibitor of Plasmodium falciparum Histone Deacetylase 1 (PfHDAC1): An in Silico Approach

A new Plasmodium falciparum histone deacetylase1 (PfHDAC1) homology model was built based on the highest sequence identity available template human histone deacetylase 2 structure. The generated model was carefully evaluated for stereochemical accuracy, folding correctness and overall structure quality. All evaluations were acceptable and consistent. Docking a group of hydroxamic acid histone deacetylase inhibitors and valproic acid has shown binding poses that agree well with inhibitor-bound histone deacetylase-solved structural interactions. Docking affinity dG scores were in agreement with available experimental binding affinities. Further, enzyme-ligand complex stability and reliability were investigated by running 5-nanosecond molecular dynamics simulations. Thorough analysis of the simulation trajectories has shown that enzyme-ligand complexes were stable during the simulation period. Interestingly, the calculated theoretical binding energies of the docked hydroxamic acid inhibitors have shown that the model can discriminate between strong and weaker inhibitors and agrees well with the experimental affinities reported in the literature. The model and the docking methodology can be used in screening virtual libraries for PfHDAC1 inhibitors, since the docking scores have ranked ligands in accordance with experimental binding affinities. Valproic acid calculated theoretical binding energy suggests that it may inhibit PfHDAC1.     

Molecular dynamics simulations of the docking of substituted N5- deazapterins to dihydrofolate reductase

Orientations of the deazapterin ring and the conformational preferences of groups appended to the deazapterin ring in a set of 8-substituted deazapterin cations docked into the dihydrofolate reductase (DHFR) binding site have been investigated using a methodology based on the simulated annealing technique within molecular dynamics (MD) simulations. Of five possible binding pockets for the 8-substituents, identified from a preliminary manual docking study, one has been definitively eliminated after an analysis of MD trajectories, while another remains uncertain. Using a new method based on standard thermodynamic cycles and a linear approximation of polar and non-polar free energy contributions from MD averages, binding affinities of the different ligands in each binding site have been correlated with experimental dissociation constants. The study has provided insights into structure-activity relationships for use in the design of modified inhibitors of DHFR.      

Synthesis,Crystallization Studies,and in vitro Characterization of Cinnamic Acid Derivatives as SmHDAC8 Inhibitors for the Treatment of Schistosomiasis

Schistosomiasis is a neglected parasitic disease that affects more than 265 million people worldwide and for which the control strategy relies on mass treatment with only one drug: praziquantel. Based on the 3‐chlorobenzothiophene‐2‐hydroxamic acid J1075 , a series of hydroxamic acids with different scaffolds were prepared as potential inhibitors of Schistosoma mansoni histone deacetylase 8 (SmHDAC8). The crystal structures of SmHDAC8 with four inhibitors provided insight into the binding mode and orientation of molecules in the binding pocket as well as the orientation of its flexible amino acid residues. The compounds were evaluated in screens for inhibitory activity against schistosome and human HDACs. The most promising compounds were further investigated for their activity toward the major human HDAC isotypes. The most potent inhibitors were additionally screened for lethality against the schistosome larval stage using a fluorescence‐based assay. Two of the compounds showed significant, dose‐dependent killing of the schistosome larvae and markedly impaired egg laying of adult worm pairs maintained in culture.