Ligand-Based Virtual Screening and Molecular Docking of Benzimidazoles as Potential Inhibitors of Triosephosphate Isomerase Identified New Trypanocidal Agents

Trypanosoma cruzi (T. cruzi) is a parasite that affects humans and other mammals. T. cruzi depends on glycolysis as a source of adenosine triphosphate (ATP) supply, and triosephosphate isomerase (TIM) plays a key role in this metabolic pathway. This enzyme is an attractive target for the design of new trypanocidal drugs. In this study, a ligand-based virtual screening (LBVS) from the ZINC15 database using benzimidazole as a scaffold was accomplished. Later, a molecular docking on the interface of T. cruzi TIM (TcTIM) was performed and the compounds were grouped by interaction profiles. Subsequently, a selection of compounds was made based on cost and availability for in vitro evaluation against blood trypomastigotes. Finally, the compounds were analyzed by molecular dynamics simulation, and physicochemical and pharmacokinetic properties were determined using SwissADME software. A total of 1604 molecules were obtained as potential TcTIM inhibitors. BP2 and BP5 showed trypanocidal activity with half-maximal lytic concentration (LC₅₀) values of 155.86 and 226.30 µM, respectively. Molecular docking and molecular dynamics simulation analyzes showed a favorable docking score of BP5 compound on TcTIM. Additionally, BP5 showed a low docking score (−5.9 Kcal/mol) on human TIM compared to the control ligand (−7.2 Kcal/mol). Both compounds BP2 and BP5 showed good physicochemical and pharmacokinetic properties as new anti-T. cruzi agents.     

Synthesis,Biological Activity,and Molecular Dynamics Study of Novel Series of a Trimethoprim Analogs as Multi-Targeted Compounds: Dihydrofolate Reductase (DHFR) Inhibitors and DNA-Binding Agents

Eighteen previously undescribed trimethoprim (TMP) analogs containing amide bonds (1–18) were synthesized and compared with TMP, methotrexate (MTX), and netropsin (NT). These compounds were designed as potential minor groove binding agents (MGBAs) and inhibitors of human dihydrofolate reductase (hDHFR). The all-new derivatives were obtained via solid phase synthesis using 4-nitrophenyl Wang resin. Data from the ethidium displacement test confirmed their DNA-binding capacity. Compounds 13–14 (49.89% and 43.85%) and 17–18 (41.68% and 42.99%) showed a higher binding affinity to pBR322 plasmid than NT. The possibility of binding in a minor groove as well as determination of association constants were performed using calf thymus DNA, T4 coliphage DNA, poly (dA-dT)², and poly (dG-dC)². With the exception of compounds 9 (IC50 = 56.05 µM) and 11 (IC50 = 55.32 µM), all of the compounds showed better inhibitory properties against hDHFR than standard, which confirms that the addition of the amide bond into the TMP structures increases affinity towards hDHFR. Derivatives 2, 6, 13, 14, and 16 were found to be the most potent hDHFR inhibitors. This molecular modelling study shows that they interact strongly with a catalytically important residue Glu-30.     

Activation of Sirtuin 2 Inhibitors Employing Photoswitchable Geometry and Aqueous Solubility

Because isoenzymes of the experimentally and therapeutically extremely relevant sirtuin family show high similarity, addressing the unique selectivity pocket of sirtuin 2 is a promising strategy towards selective inhibitors. An unrelated approach towards selective inhibition of isoenzymes with varied tissue distribution is targeted drug delivery or spatiotemporal activation by photochemical activation. Azologization of two nicotinamide-mimicking lead structures was undertaken to combine both approaches and yielded a set of 33 azobenzenes and azopyridines that have been evaluated for their photochemical behaviour and bioactivity. For some compounds, inhibitory activity reached the sub-micromolar range in their thermodynamically favoured E form and could be decreased by photoisomerization to the metastable Z form. Besides, derivatization with long-chain fatty acids yielded potent sirtuin 2 inhibitors, featuring another intriguing aspect of azo-based photoswitches. In these compounds, switching to the Z isomer increased aqueous solubility and thereby enhanced biological activity by up to a factor of 21. The biological activity of two compounds was confirmed by hyperacetylation of sirtuin specific histone proteins in a cell-based activity assay.     

Novel Sulfonamide-Based Carbamates as Selective Inhibitors of BChE

A series of 14 target benzyl [2-(arylsulfamoyl)-1-substituted-ethyl]carbamates was prepared by multi-step synthesis and characterized. All the final compounds were tested for their ability to inhibit acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) in vitro, and the selectivity index (SI) was determined. Except for three compounds, all compounds showed strong preferential inhibition of BChE, and nine compounds were even more active than the clinically used rivastigmine. Benzyl {(2S)-1-[(2-methoxybenzyl)sulfamoyl]-4-methylpentan-2-yl}carbamate (5k), benzyl {(2S)-1-[(4-chlorobenzyl)sulfamoyl]-4-methylpentan-2-yl}carbamate (5j), and benzyl [(2S)-1-(benzylsulfamoyl)-4-methylpentan-2-yl]carbamate (5c) showed the highest BChE inhibition (IC₅₀ = 4.33, 6.57, and 8.52 µM, respectively), indicating that derivatives 5c and 5j had approximately 5-fold higher inhibitory activity against BChE than rivastigmine, and 5k was even 9-fold more effective than rivastigmine. In addition, the selectivity index of 5c and 5j was approx. 10 and that of 5k was even 34. The process of carbamylation and reactivation of BChE was studied for the most active derivatives 5k, 5j. The detailed information about the mode of binding of these compounds to the active site of both BChE and AChE was obtained in a molecular modeling study. In this study, combined techniques (docking, molecular dynamic simulations, and QTAIM (quantum theory of atoms in molecules) calculations) were employed.     

新型乙酰胆碱酯酶抑制剂的设计与虚拟筛选

运用分子对接软件GOLD,研究了延胡索类生物碱corydaline与乙酰胆碱酯酶的结合模型,以coryaaline为先导化合物,设计与虚拟筛选一系列四氢异喹啉衍生物。预测活性做好的化合物1A环与Trp84发生π-π堆积;质子化的氮原子与Phe330发生阳离子-π作用;B环与Trp279发生π-π作用;1位的羟甲基、4位甲氧基、C环上的甲氧基分别与Tyr334、Tyr21、Tyrl30形成氢键作用。     

Investigation into the Feasibility of Thioditaloside as a Novel Scaffold for Galectin‐3‐Specific Inhibitors

Galectin‐3 is extensively involved in metabolic and disease processes, such as cancer metastasis, thus giving impetus for the design of specific inhibitors targeting this β‐galactose‐binding protein. Thiodigalactoside (TDG) presents a scaffold for construction of galectin inhibitors, and its inhibition of galectin‐1 has already demonstrated beneficial effects as an adjuvant with vaccine immunotherapy, thereby improving the survival outcome of tumour‐challenged mice. A novel approach—replacing galactose with its C2 epimer, talose—offers an alternative framework, as extensions at C2 permit exploitation of a galectin‐3‐specific binding groove, thereby facilitating the design of selective inhibitors. We report the synthesis of thioditaloside (TDT) and crystal structures of the galectin‐3 carbohydrate recognition domain in complexes with TDT and TDG. The different abilities of galactose and talose to anchor to the protein correlate with molecular dynamics studies, likely explaining the relative disaccharide binding affinities. The feasibility of a TDT scaffold to enable access to a particular galectin‐3 binding groove and the need for modifications to optimise such a scaffold for use in the design of potent and selective inhibitors are assessed.     

Structural Determinants of CX-4945 Derivatives as Protein Kinase CK2 Inhibitors: A Computational Study

Protein kinase CK2, also known as casein kinase-2, is involved in a broad range of physiological events including cell growth, proliferation and suppression of apoptosis which are related to human cancers. A series of compounds were identified as CK2 inhibitors and their inhibitory activities varied depending on their structures. In order to explore the structure-activity correlation of CX-4945 derivatives as inhibitors of CK2, in the present study, a set of ligand- and receptor-based 3D-QSAR models were developed employing Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Index Analysis (CoMSIA). The optimum CoMFA (R(cv) (2) = 0.618, R(pred) (2) = 0.892) and CoMSIA (R(cv) (2) = 0.681, R(pred) (2) = 0.843) models exhibited reasonable statistical characteristics for CX-4945 derivatives. The results indicated that electrostatic effects contributed the most to both CoMFA and CoMSIA models. The combination of docking analysis and molecular dynamics (MD) simulation showed that Leu45, Lys68, Glu81, Val116, Asp175 and Trp176 of CK2 which formed several direct or water-bridged H-bonds with CX-4945 are crucial for CX-4945 derivatives recognition to CK2. These results can offer useful theoretical references for designing more potent CK2 inhibitors.     

Docking and Molecular Dynamic of Microalgae Compounds as Potential Inhibitors of Beta-Lactamase

Bacterial resistance is responsible for a wide variety of health problems, both in children and adults. The persistence of symptoms and infections are mainly treated with β-lactam antibiotics. The increasing resistance to those antibiotics by bacterial pathogens generated the emergence of extended-spectrum β-lactamases (ESBLs), an actual public health problem. This is due to rapid mutations of bacteria when exposed to antibiotics. In this case, β-lactamases are enzymes used by bacteria to hydrolyze the beta-lactam rings present in the antibiotics. Therefore, it was necessary to explore novel molecules as potential β-lactamases inhibitors to find antibacterial compounds against infection caused by ESBLs. A computational methodology based on molecular docking and molecular dynamic simulations was used to find new microalgae metabolites inhibitors of β-lactamase. Six 3D β-lactamase proteins were selected, and the molecular docking revealed that the metabolites belonging to the same structural families, such as phenylacridine (4-Ph), quercetin (Qn), and cryptophycin (Cryp), exhibit a better binding score and binding energy than commercial clinical medicine β-lactamase inhibitors, such as clavulanic acid, sulbactam, and tazobactam. These results indicate that 4-Ph, Qn, and Cryp molecules, homologous from microalgae metabolites, could be used, likely as novel β-lactamase inhibitors or as structural templates for new in-silico pharmaceutical designs, with the possibility of combatting β-lactam resistance     

A Computational Study on Thiourea Analogs as Potent MK-2 Inhibitors

Mitogen-activated protein kinase-activated protein kinase 2 (MK-2) has been identified as a drug target for the treatment of inflammatory diseases. Currently, a series of thiourea analogs as potent MK-2 inhibitors were studied using comprehensive computational methods by 3D-QSAR, molecular docking and molecular dynamics simulations for a further improvement in activities. The optimal 3D models exhibit high statistical significance of the results, especially for the CoMFA results with r(2) (ncv), q(2) values of 0.974, 0.536 for the internal validation, and r(2) (pred), r(2) (m) values of 0.910, 0.723 for the external validation and Roy's index, respectively. In addition, more rigorous validation criteria suggested by Tropsha were also employed to check the built models. Graphic representation of the results, as contoured 3D coefficient plots, also provides a clue to the reasonable modification of molecules: (i) The substituent with a bulky size and electron-rich group at the C5 position of the pyrazine ring is required to enhance the potency; (ii) The H-bond acceptor group in the C3 position of the pyrazine ring is likely to be helpful to increase MK-2 inhibition; (iii) The small and electropositive substituent as a hydrogen bond donor of the C2 position in the oxazolone ring is favored; In addition, several important amino acid residues were also identified as playing an important role in MK-2 inhibition. The agreement between 3D-QSAR, molecular docking and molecular dynamics simulations also proves the rationality of the developed models. These results, we hope, may be helpful in designing novel and potential MK-2 inhibitors.     

Computational Design of Miniproteins as SARS-CoV-2 Therapeutic Inhibitors

A rational therapeutic strategy is urgently needed for combating SARS-CoV-2 infection. Viral infection initiates when the SARS-CoV-2 receptor-binding domain (RBD) binds to the ACE2 receptor, and thus, inhibiting RBD is a promising therapeutic for blocking viral entry. In this study, the structure of lead antiviral candidate binder (LCB1), which has three alpha-helices (H1, H2, and H3), is used as a template to design and simulate several miniprotein RBD inhibitors. LCB1 undergoes two modifications: structural modification by truncation of the H3 to reduce its size, followed by single and double amino acid substitutions to enhance its binding with RBD. We use molecular dynamics (MD) simulations supported by ab initio density functional theory (DFT) calculations. Complete binding profiles of all miniproteins with RBD have been determined. The MD investigations reveal that the H3 truncation results in a small inhibitor with a −1.5 kcal/mol tighter binding to RBD than original LCB1, while the best miniprotein with higher binding affinity involves D17R or E11V + D17R mutation. DFT calculations provide atomic-scale details on the role of hydrogen bonding and partial charge distribution in stabilizing the minibinder:RBD complex. This study provides insights into general principles for designing potential therapeutics for SARS-CoV-2.     

A Hamiltonian Replica Exchange Molecular Dynamics (MD) Method for the Study of Folding,Based on the Analysis of the Stabilization Determinants of Proteins

Herein, we present a novel Hamiltonian replica exchange protocol for classical molecular dynamics simulations of protein folding/unfolding. The scheme starts from the analysis of the energy-networks responsible for the stabilization of the folded conformation, by means of the energy-decomposition approach. In this framework, the compact energetic map of the native state is generated by a preliminary short molecular dynamics (MD) simulation of the protein in explicit solvent. This map is simplified by means of an eigenvalue decomposition. The highest components of the eigenvector associated with the lowest eigenvalue indicate which sites, named “hot spots”, are likely to be responsible for the stability and correct folding of the protein. In the Hamiltonian replica exchange protocol, we use modified force-field parameters to treat the interparticle non-bonded potentials of the hot spots within the protein and between protein and solvent atoms, leaving unperturbed those relative to all other residues, as well as solvent-solvent interactions. We show that it is possible to reversibly simulate the folding/unfolding behavior of two test proteins, namely Villin HeadPiece HP35 (35 residues) and Protein A (62 residues), using a limited number of replicas. We next discuss possible implications for the study of folding mechanisms via all atom simulations.     

Computer-Aided and AILDE Approaches to Design Novel 4-Hydroxyphenylpyruvate Dioxygenase Inhibitors

4-Hydroxyphenylpyruvate dioxygenase (HPPD) is a pivotal enzyme in tocopherol and plastoquinone synthesis and a potential target for novel herbicides. Thirty-five pyridine derivatives were selected to establish a Topomer comparative molecular field analysis (Topomer CoMFA) model to obtain correlation information between HPPD inhibitory activity and the molecular structure. A credible and predictive Topomer CoMFA model was established by “split in two R-groups” cutting methods and fragment combinations (q² = 0.703, r² = 0.957, ONC = 6). The established model was used to screen out more active compounds and was optimized through the auto in silico ligand directing evolution (AILDE) platform to obtain potential HPPD inhibitors. Twenty-two new compounds with theoretically good HPPD inhibition were obtained by combining the high-activity contribution substituents in the existing molecules with the R-group search via Topomer search. Molecular docking results revealed that most of the 22 fresh compounds could form stable π-π interactions. The absorption, distribution, metabolism, excretion and toxicity (ADMET) prediction and drug-like properties made 9 compounds potential HPPD inhibitors. Molecular dynamics simulation indicated that Compounds Y12 and Y14 showed good root mean square deviation (RMSD) and root mean square fluctuation (RMSF) values and stability. According to the AILDE online verification, 5 new compounds with potential HPPD inhibition were discovered as HPPD inhibitor candidates. This study provides beneficial insights for subsequent HPPD inhibitor design.     

Discovery of Novel Dihydrolipoamide S-Succinyltransferase Inhibitors Based on Fragment Virtual Screening

A series of new oxadiazole sulfone derivatives containing an amide moiety was synthesized based on fragment virtual screening to screen high-efficiency antibacterial agents for rice bacterial diseases. All target compounds showed greater bactericidal activity than commercial bactericides. 3-(4-fluorophenyl)-N-((5-(methylsulfonyl)-1,3,4-oxadiazol-2-yl)methyl)acrylamide (10) showed excellent antibacterial activity against Xanthomonas oryzae pv. oryzae and Xanthomonas oryzae pv. oryzicola, with EC₅₀ values of 0.36 and 0.53 mg/L, respectively, which were superior to thiodiazole copper (113.38 and 131.54 mg/L) and bismerthiazol (83.07 and 105.90 mg/L). The protective activity of compound 10 against rice bacterial leaf blight and rice bacterial leaf streak was 43.2% and 53.6%, respectively, which was superior to that of JHXJZ (34.1% and 26.4%) and thiodiazole copper (33.0% and 30.2%). The curative activity of compound 10 against rice bacterial leaf blight and rice bacterial leaf streak was 44.5% and 51.7%, respectively, which was superior to that of JHXJZ (32.6% and 24.4%) and thiodiazole copper (27.1% and 28.6%). Moreover, compound 10 might inhibit the growth of Xanthomonas oryzae pv. oryzae and Xanthomonas oryzae pv. oryzicola by affecting the extracellular polysaccharides, destroying cell membranes, and inhibiting the enzyme activity of dihydrolipoamide S-succinyltransferase.     

Nontoxic and Naturally Occurring Active Compounds as Potential Inhibitors of Biological Targets in Liriomyza trifolii

In recent years, novel strategies to control insects have been based on protease inhibitors (PIs). In this regard, molecular docking and molecular dynamics simulations have been extensively used to investigate insect gut proteases and the interactions of PIs for the development of resistance against insects. We, herein, report an in silico study of (disodium 5′-inosinate and petunidin 3-glucoside), (calcium 5′-guanylate and chlorogenic acid), chlorogenic acid alone, (kaempferol-3,7-di-O-glucoside with hyperoside and delphinidin 3-glucoside), and (myricetin 3′-glucoside and hyperoside) as potential inhibitors of acetylcholinesterase receptors, actin, α-tubulin, arginine kinase, and histone receptor III subtypes, respectively. The study demonstrated that the inhibitors are capable of forming stable complexes with the corresponding proteins while also showing great potential for inhibitory activity in the proposed protein-inhibitor combinations.     

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.     

Analysis of Binding Determinants for Different Classes of Competitive and Noncompetitive Inhibitors of Glycine Transporters

Glycine transporters are interesting therapeutic targets as they play significant roles in glycinergic and glutamatergic systems. The search for new selective inhibitors of particular types of glycine transporters (GlyT-1 and GlyT-2) with beneficial kinetics is hampered by limited knowledge about the spatial structure of these proteins. In this study, a pool of homology models of GlyT-1 and GlyT-2 in different conformational states was constructed using the crystal structures of related transporters from the SLC6 family and the recently revealed structure of GlyT-1 in the inward-open state, in order to investigate their binding sites. The binding mode of the known GlyT-1 and GlyT-2 inhibitors was determined using molecular docking studies, molecular dynamics simulations, and MM-GBSA free energy calculations. The results of this study indicate that two amino acids, Gly373 and Leu476 in GlyT-1 and the corresponding Ser479 and Thr582 in GlyT-2, are mainly responsible for the selective binding of ligands within the S1 site. Apart from these, one pocket of the S2 site, which lies between TM3 and TM10, may also be important. Moreover, selective binding of noncompetitive GlyT-1 inhibitors in the intracellular release pathway is affected by hydrophobic interactions with Ile399, Met382, and Leu158. These results can be useful in the rational design of new glycine transporter inhibitors with desired selectivity and properties in the future.