面向虚拟筛选的GPU加速的分子对接方法

虚拟筛选是在计算机上对化合物分子进行模拟预筛选,找出容易和药物靶标结合的小分子(配体),从而降低实际实验测试次数,提高药物先导化合物的发现效率。常用的分子对接软件可以用于基于结构的虚拟筛选,寻找配体与靶标的最佳的作用模式和结合构象,并通过打分函数来筛选出潜在的配体。现有的对接软件如AutoDock Vina等在分子对接过程中需要耗费大量时间和计算资源,特别是面对大规模分子对接时,过长的筛选时间不能满足应用需求,因此,本文在最高效的QVina2对接软件基础上,提出一种基于GPU的QVina 2并行化方法QVina2-GPU,利用GPU硬件高度并行体系加速分子对接。具体包括增加初始化分子构象数量,以扩展蒙特卡罗的迭代局部搜索中线程的并行规模,增加蒙特卡罗的迭代搜索的广度以减少每次蒙特卡罗迭代搜索深度,并利用Wolfe-Powell准则改进局部搜索算法,提高了对接精度,进一步减少蒙特卡罗迭代搜索深度,最后,在NVIDIA Geforce RTX 3090平台上在公开的配体数据库上验证了QVina2-GPU的性能,实验表明在保证分子对接精度的基础上,我们提出的QVina2-GPU对Qvina2的平均加速比达到5.18倍,最大加速比达到12.28倍。     

Synthesis,molecular docking and biological evaluation of metronidazole derivatives containing piperazine skeleton as potential antibacterial agents

Metronidazole has a broad-spectrum antibacterial activity. Hereby a series of novel metronidazole derivatives were designed and synthesized based on nitroimidazole scaffold in order to find some more potent antibacterial drugs. For these compounds which were reported for the first time, their antibacterial activities against Escherichia coli, Pseudomonas aeruginosa, Bacillus subtilis and Staphylococcus aureus were tested. These compounds showed good antibacterial activities against Gram-positive strains. Compound 4m represented the most potent antibacterial activity against S. aureus ATCC 25923 with MIC of 0.003 μg/mL and it showed the most potent activity against S. aureus TyrRS with IC₅₀ of 0.0024 μM. Molecular docking of 4m into S. aureus tyrosyl-tRNA synthetase active site were also performed to determine the probable binding mode.     

Proteome based mapping and molecular docking revealed DnaA as a potential drug target against Shigella sonnei

Shigella sonnei is one of the major causes of diarrhea and remained a critical microbe responsible for higher morbidity and mortality rates resulting from dysentery every year across the world. Antibiotic therapy of Shigella diseases plays a critical role in decreasing the prevalence as well as the fatality rate of this infection. However, the management of these diseases remains challenging, owing to the overall increase in resistance against many antimicrobials. The situation necessitates the rapid development of effective and feasible S. sonnei treatments. In the present study, the subtractive genomics approach was utilized to find the potential drug targets for S. sonnei strain Ss046. Various tools of bioinformatics were implemented to remove the human-specific homologous and pathogen-specific paralogous sequences from the bacterial proteome. Then, metabolic pathway and subcellular location analysis were performed of essential bacterial proteins to describe their role in various cellular processes. Only one essential protein i-e Chromosomal replication initiator protein DnaA was found in the proteome of the pathogen that could be used as a potent target for designing new drugs. 3D structure prediction of DnaA protein was carried out using Phyre 2. Molecular docking of 5000 phytochemicals was performed against DnaA to identify four top-ranked phytochemicals (Riccionidin A, Dothistromin, Fustin, and Morin) based on scoring functions and interaction with the active site. This study suggests that these phytochemicals could be used as antibacterial drugs to treat S. sonnei infections in the future. To confirm their efficacy and evaluate their drug potency, further in vitro analyses are required.     

17(E)-Picolinylidene androstane derivatives as potential inhibitors of prostate cancer cell growth: Antiproliferative activity and molecular docking studies

We report a rapid and efficient synthesis of A-ring modified 17α-picolyl and 17(E)-picolinylidene androstane derivatives from dehydroepiandrosterone. Compounds were validated spectroscopically and structurally characterized by X-ray crystallography. Virtual screening by molecular docking against clinical targets of steroidal anticancer drugs (ERα, AR, Aromatase and CYP17A1) suggests that 17(E)-picolinylidene, but not 17α-picolyl androstanes could specifically interact with CYP17A1 (17α-hydroxylase) with similar geometry and affinity as Abiraterone, a 17-pyridinyl androstane drug clinically used in the treatment of prostate cancer. In addition, several 17(E)-picolinylidene androstanes demonstrated selective antiproliferative activity against PC3 prostate cancer cells, which correlates with Abiraterone antiproliferative activity and predicted CYP17A1 binding affinities. Based on these preliminary results, 17(E)-picolinylidene androstane derivatives could be a promising starting point for the development of new compounds for the treatment of prostate cancer.     

Molecular docking based screening of predicted potential inhibitors for VP40 from Ebola virus

Ebola virus is a member of Filoviridae and cause severe human disease with 90 percent mortality. The life cycle of Ebola contains an assembly stage which is mediated by VP40 proteins. VP40 subunits oligomerize and form ring-structures which are either octamers or hexamers. Prevention of VP40 matrix protein assembly prevents virus particle formation as well as virus budding. In the present study we simulated the biological condition for a single VP40 subunit. Then a library containing 120.000 drugs like chemicals was used as the virtual screening database. Top 10 successive hits were then analyzed regarding absorption, distribution, metabolism, and excretion properties. Moreover probable accessorial human protein targets and toxicity properties of successive hits were analyzed by in silico tools. We found 4 chemicals that could bind VP40 subunits in a manner that by making an interfering steric condense prevents matrix protein oligomerization. The pharmacokinetic and toxicity studies also validated the potential of 4 finlay successive hits to be considered as a new anti-Ebola drugs.     

Pharmacoinformatics and molecular docking reveal potential drug candidates against Schizophrenia to target TAAR6

Schizophrenia (SZ) is a complex disabling disorder that leads to the mental disability and afflicts 1% of the world's total population and placed in top ten medical disorders. In current work, bioinformatics analyses were carried out on Trace amine (TA)-associated receptor 6 (TAAR6) to recognize the potential drugs and compounds against SZ. Comparative modeling and threading-based approaches were utilized for the structure prediction of TAAR6. Fifty-nine predicted structures were evaluated by various model assessment techniques and final model having only eight amino acids in the outlier region and 98.5% overall quality factor was chosen for further pharmacoinformatics and molecular docking analyses. From an extensive literature review, 11 Food and Drug Administration (FDA) approved drugs were analyzed by computational techniques and Aripiprazole was found as the most effective drug against SZ by targeting TAAR6. Here, we report five novel molecules which exhibited the highest binding affinity, effective drug properties, and interestingly, observed better results than the approved selected drugs against SZ by targeting TAAR6. The docking analyses revealed that Arg-92, Trp-98, Gln-191, Thr-192, Ala-290, Cys-291, Tyr-293, and Glu-294 residues were observed as critical interacting residues in receptor-ligand interactions. Absorption, distribution, metabolism, excretion, and toxicity (ADMET) properties, Lipinski rule of five, highest binding affinity coupled with virtual screening (VS), and pharmacophore modeling approach illustrated that aripiprazole (−8.6 kcal/mol) and TAAR6_0094 (−9.3 kcal/mol) are potential inhibitors for targeting TAAR6. It is suggested that schizophrenic patients have to use Aripiprazole for the medication of SZ by targeting TAAR6 and develop effective therapies by utilizing scrutinized novel compound.     

Synthesis,molecular docking and biological evaluation of coumarin derivatives containing piperazine skeleton as potential antibacterial agents

A series of 4-hydroxycoumarin derivatives were designed and synthesized in order to find some more potent antibacterial drugs. Their antibacterial activities against Escherichia coli, Pseudomonas aeruginosa, Bacillus subtilis and Staphylococcus aureus were tested. These compounds showed good antibacterial activities against Gram-positive strains. Compound 4g represented the most potent antibacterial activity against Bacillus subtilis and S. aureus with MIC of 0.236, 0.355 μg/mL, respectively. What’s more, it showed the most potent activity against SaFabI with IC₅₀ of 0.57 μM. Molecular docking of 4g into S. aureus Enoyl-ACP-reductase active site were performed to determine the probable binding mode, while the QSAR model was built to check the previous work as well as to introduce new directions.     

Identification of influenza PA-Nter endonuclease inhibitors using pharmacophore- and docking-based virtual screening

Searching for new antiviral agents, we focused our interest on the influenza PA-Nter endonuclease. Therefore, we developed a three-dimensional pharmacophore model which contains the binding features addressed to the metal-chelating active site. The obtained hypothesis has been fruitfully employed to select three “hit compounds” through an in silico screening campaign on our in-house database of small molecules. We studied the binding poses of these hit compounds using molecular docking, and subjected them to an enzymatic assay with recombinant PA-Nter endonuclease. Compound 20 proved the most active inhibitor of the endonucleolytic cleavage reaction, with an IC₅₀ value of 12?μM.     

Identification of nonplanar small molecule for G-quadruplex grooves: molecular docking and molecular dynamic study

DNA G-quadruplex is an attractive drug target for anticancer therapy. Most G-quadruplex ligands have little selectivity, due to π-stacking interaction with common G-tetrads surface. Thanks to the varieties of G-quadruplex grooves, the groove-binding ligand is expected to create high selectivity. Therefore, developing novel molecular geometries that target G-quadruplex groove has been paid growing attention. In this work, steroid FG, a special nonplanar and nonaromatic small molecule, interacting with different conformations of G-quadruplexes has been studied by molecular docking and molecular dynamics simulations. The results showed the selectivity of the hydrophobic group of steroid FG for the wide groove of antiparallel G-quadruplex. The methyl groups on the tetracyclic ring of steroid represent the specific binding ability for the small hydrophobic cavity formed by reversed stacking of G-tetrads in antiparallel G-quadruplex groove. This work provides new insight for developing new classes of G-quadruplex groove-binding ligands.     

新型冠状病毒相关靶点小分子抑制剂虚拟筛选

为确定治疗新型冠状病毒(SARS-CoV-2)感染的候选药物,开展了针对SARS-CoV-2的药物虚拟筛选研究。以SARS-CoV-2的刺突蛋白(S蛋白)和3CL蛋白酶(主蛋白酶)作为药物靶点,以美国食品药品监督管理局(FDA)批准上市的2 726个小分子药物作为候选,通过分子对接方法,筛选出了3种(阿巴瑞克(Abarelix)、西曲瑞克(Cetrorelix)、鞣酸(Tannic acid))与S蛋白具有较强结合能力的小分子药物,1种(戈舍瑞林(Goserelin))与3CL蛋白酶具有较好结合能力的小分子药物,它们理论上都具有抑制新型冠状病毒复制的效果。将靶向3CL蛋白酶的候选药物与辉瑞公司开发的药物Paxlovid进行比较,发现其作用位点均集中于3CL蛋白酶的第130-200位的残基周围,具有相似的结合位点与相互作用。此外也对候选药物的物理与化学性质及与基因相互作用进行了分析。本研究可为开发新型冠状病毒感染的治疗药物提供参考。     

Identification of novel acetylcholinesterase inhibitors: Indolopyrazoline derivatives and molecular docking studies

The synthesis of novel indolopyrazoline derivatives (P1-P4 and Q1-Q4) has been characterized and evaluated as potential anti-Alzheimer agents through in vitro Acetylcholinesterase (AChE) inhibition and radical scavenging activity (antioxidant) studies. Specifically, Q3 shows AChE inhibition (IC₅₀: 0.68 ± 0.13 μM) with strong DPPH and ABTS radical scavenging activity (IC₅₀: 13.77 ± 0.25 μM and IC₅₀: 12.59 ± 0.21 μM), respectively. While P3 exhibited as the second most potent compound with AChE inhibition (IC₅₀: 0.74 ± 0.09 μM) and with DPPH and ABTS radical scavenging activity (IC₅₀: 13.52 ± 0.62 μM and IC₅₀: 13.13 ± 0.85 μM), respectively. Finally, molecular docking studies provided prospective evidence to identify key interactions between the active inhibitors and the AChE that furthermore led us to the identification of plausible binding mode of novel indolopyrazoline derivatives. Additionally, in-silico ADME prediction using QikProp shows that these derivatives fulfilled all the properties of CNS acting drugs. This study confirms the first time reporting of indolopyrazoline derivatives as potential anti-Alzheimer agents.     

Identification of natural compounds as potent inhibitors of SARS-CoV-2 main protease using combined docking and molecular dynamics simulations

Coronavirus disease 2019 (COVID-19) has emerged from China and globally affected the entire population through the human-to-human transmission of a newly emerged virus called severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The genome of SARS-CoV-2 encodes several proteins that are essential for multiplication and pathogenesis. The main protease (M^pro or 3CL^pro) of SARS-CoV-2 plays a central role in its pathogenesis and thus is considered as an attractive drug target for the drug design and development of small-molecule inhibitors. We have employed an extensive structure-based high-throughput virtual screening to discover potential natural compounds from the ZINC database which could inhibit the M^pro of SARS-CoV-2. Initially, the hits were selected on the basis of their physicochemical and drug-like properties. Subsequently, the PAINS filter, estimation of binding affinities using molecular docking, and interaction analyses were performed to find safe and potential inhibitors of SARS-CoV-2 M^pro. We have identified ZINC02123811 (1-(3-(2,5,9-trimethyl-7-oxo-3-phenyl-7H-furo[3,2-g]chromen-6-yl)propanoyl)piperidine-4-carboxamide), a natural compound bearing appreciable affinity, efficiency, and specificity towards the binding pocket of SARS-CoV-2 M^pro. The identified compound showed a set of drug-like properties and preferentially binds to the active site of SARS-CoV-2 M^pro. All-atom molecular dynamics (MD) simulations were performed to evaluate the conformational dynamics, stability and interaction mechanism of M^pro with ZINC02123811. MD simulation results indicated that M^pro with ZINC02123811 forms a stable complex throughout the trajectory of 100 ns. These findings suggest that ZINC02123811 may be further exploited as a promising scaffold for the development of potential inhibitors of SARS-CoV-2 M^pro to address COVID-19.     

Synthesis,biological evaluation and molecular docking analysis of novel benzopyrimidinone derivatives as potential anti-tyrosinase agents

2-substitued-benzopyrimidinones 2 were synthesized in high to excellent yields in a single step via condensation of 2-aminobenzamide 1 with some aryl-aldehydes in the presence of iodine. Cyclocondensation reaction of hydrazides 3 which were obtained in two steps from benzopyrimidinones 2, with some electrophilic species such as 2,4-pentandione, 2,5-hexandione, 1-phenylbutan-1,3-dione and cyclic anyhdrides provided the new compounds 4a–c, 5a–c, 6a–c, 7a–c, 8a–c and 9a–c. The synthesized compounds were characterized by spectroscopic means. They were also evaluated for their anti-tyrosinase potential. The structure-activity relationship (SAR) was discussed on the basis of the molecular docking analysis.     

Discovery of the first SecA inhibitors using structure-based virtual screening

Bacterial protein secretion is a critical and complex process. The Sec machinery provides a major pathway for protein translocation across and integration into the cellular membrane in bacteria. Small molecule probes that perturb the functions of individual member proteins within the Sec machinery will be very important research tools as well as leads for future antimicrobial agent development. Herein we describe the discovery of inhibitors, through virtual screening, that specifically act on SecA ATPase, which is a critical member of the Sec system. These are the very first inhibitors reported for intrinsic SecA ATPase.     

In silico study of M18 aspartyl amino peptidase (M18AAP) of Plasmodium vivax as an antimalarial drug target

Plasmodium vivax (Pv) is the second most malaria causing pathogen among Plasmodium species. M18 aspartic aminopeptidase (M18AAP) protein is a single gene copy present in Plasmodium. This protein is functional at the terminal stage of hemoglobin degradation of host and completes the hydrolysis process which makes it an important target for new chemotherapeutics. No experimental and structural study on M18AAP protein of P. vivax is reported till today. This paper advocates the application of multiple computational approaches like protein model prediction, ligand-based 3D QSAR study, pharmacophore, structure-based virtual screening and molecular docking simulation for identification of potent lead molecules against the enzyme. The 3D QSAR model was developed using known bioactive compounds against the PvM18AAP protein which statistically signify the k-NN model with q^2 = 0.7654. The study reports a lead molecule from ligand-centric approach with good binding affinity and possessing lowest docking score. The findings will be helpful for in-vivo and in-vitro validations and development of potent anti-malarial molecules against the drug resistant strains of malaria parasite.     

Synthesis,in vitro urease inhibitory potential and molecular docking study of Benzimidazole analogues

Despite of many diverse biological activities exhibited by benzimidazole scaffold, it is rarely explored for the urease inhibitory potential. For that purpose, benzimidazole analogues 1–19 were synthesized and screened for in vitro urease inhibitory potential. Structures of all synthetic analogues were deduced by different spectroscopic techniques. All analogues revealed inhibition potential with IC₅₀ values of 0.90 ± 0.01 to 35.20 ± 1.10 μM, when compared with the standard thiourea (IC₅₀ = 21.40 ± 0.21 μM). Limited SAR suggested that the variations in the inhibitory potentials of the analogues are the result of different substitutions on phenyl ring. In order to rationalize the binding interactions of most active compounds with the active site of urease enzyme, molecular docking study was conducted.     

Synthesis,molecular docking study and in vitro thymidine phosphorylase inhibitory potential of oxadiazole derivatives

We have synthesized oxadiazole derivatives (1–16), characterized by ¹H NMR, ¹³C NMR and HREI-MS and screened for thymidine phosphorylase inhibitory potential. All derivatives display varied degree of thymidine phosphorylase inhibition in the range of 1.10 ± 0.05 to 49.60 ± 1.30 μM when compared with the standard inhibitor 7-Deazaxanthine having an IC₅₀ value 38.68 ± 1.12 μM. Structure activity relationships (SAR) has been established for all compounds to explore the role of substitution and nature of functional group attached to the phenyl ring which applies imperious effect on thymidine phosphorylase activity. Molecular docking study was performed to understand the binding interaction of the most active derivatives with enzyme active site.     

Discover potential inhibitors for PFKFB3 using 3D-QSAR,virtual screening,molecular docking and molecular dynamics simulation

Abstract

The 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase-3 (PFKFB3) is a master regulator of glycolysis in cancer cells by synthesizing fructose-2,6-bisphosphate (F-2,6-BP), a potent allosteric activator of phosphofructokinase-1 (PFK-1), which is a rate-limiting enzyme of glycolysis. PFKFB3 is an attractive target for cancer treatment. It is valuable to discover promising inhibitors by using 3D-QSAR pharmacophore modeling, virtual screening, molecular docking and molecular dynamics simulation. Twenty molecules with known activity were used to build 3D-QSAR pharmacophore models. The best pharmacophore model was ADHR called Hypo1, which had the highest correlation value of 0.98 and the lowest RMSD of 0.82. Then, the Hypo1 was validated by cost value method, test set method and decoy set validation method. Next, the Hypo1 combined with Lipinski's rule of five and ADMET properties were employed to screen databases including Asinex and Specs, total of 1,048,159 molecules. The hits retrieved from screening were docked into protein by different procedures including HTVS, SP and XP. Finally, nine molecules were picked out as potential PFKFB3 inhibitors. The stability of PFKFB3-lead complexes was verified by 40?ns molecular dynamics simulation. The binding free energy and the energy contribution of per residue to the binding energy were calculated by MM-PBSA based on molecular dynamics simulation.     

Anti-SARS drug screening by molecular docking

Summary. Starting from a collection of 1386 druggable compounds obtained from the 3D pharmacophore search, we performed a similarity search to narrow down the scope of docking studies. The template molecule is KZ7088 (Chou et al., 2003, Biochem Biophys Res Commun 308: 148–151). The MDL MACCS keys were used to fingerprint the molecules. The Tanimoto coefficient is taken as the metric to compare fingerprints. If the similarity threshold was 0.8, a set of 50 unique hits and 103 conformers were retrieved as a result of similarity search. The AutoDock 3.011 was used to carry out molecular docking of 50 ligands to their macromolecular protein receptors. Three compounds, i.e., C₂₈H₃₄O₄N₇Cl, C₂₁H₃₆O₅N₆, and C₂₁H₃₆O₅N₆, were found that may be promising candidates for further investigation. The main feature shared by these three potential inhibitors as well as the information of the involved side chains of SARS Cov Mpro may provide useful insights for the development of potent inhibitors against SARS enzyme.     

A hierarchical method for molecular docking using cloud computing

Discovering small molecules that interact with protein targets will be a key part of future drug discovery efforts. Molecular docking of drug-like molecules is likely to be valuable in this field; however, the great number of such molecules makes the potential size of this task enormous. In this paper, a method to screen small molecular databases using cloud computing is proposed. This method is called the hierarchical method for molecular docking and can be completed in a relatively short period of time. In this method, the optimization of molecular docking is divided into two subproblems based on the different effects on the protein–ligand interaction energy. An adaptive genetic algorithm is developed to solve the optimization problem and a new docking program (FlexGAsDock) based on the hierarchical docking method has been developed. The implementation of docking on a cloud computing platform is then discussed. The docking results show that this method can be conveniently used for the efficient molecular design of drugs.