Exploring crucial structural attributes of quinolinyl methoxyphenyl sulphonyl-based hydroxamate derivatives as ADAM17 inhibitors through classification-dependent molecular modelling approaches.

A Disintegrin and Metalloproteinase 17 (ADAM17), a Zn2+-dependent metalloenzyme of the adamalysin family of the metzincin superfamily, is associated with various pathophysiological conditions including rheumatoid arthritis and cancer. However, no specific inhibitors have been marketed yet for ADAM17-related disorders. In this study, 94 quinolinyl methoxyphenyl sulphonyl-based hydroxamates as ADAM17 inhibitors were subjected to classification-based molecular modelling and binding pattern analysis to identify the significant structural attributes contributing to ADAM17 inhibition. The statistically validated classification-based models identified the importance of the P1' substituents such as the quinolinyl methoxyphenyl sulphonyl group of these compounds for occupying the S1' - S3' pocket of the enzyme. The quinolinyl function of these compounds was found to explore stable binding of the P1' substituents at the S1' - S3' pocket whereas the importance of the sulphonyl and the orientation of the P1' moiety also revealed stable binding. Based on the outcomes of the current study, four novel compounds of different classes were designed as promising ADAM17 inhibitors. These findings regarding the crucial structural aspects and binding patterns of ADAM17 inhibitors will aid the design and discovery of novel and effective ADAM17 inhibitors for therapeutic advancements of related diseases.

Halogenated bisphenol A derivatives potently inhibit human and rat 11ß-hydroxysteroid dehydrogenase 1: Structure-activity relationship and molecular docking.

Chlorinated bisphenol A (BPA) derivatives are formed during chlorination process of drinking water, whereas bisphenol S (BPS) and brominated BPA and BPS (TBBPA and TBBPS) were synthesized for many industrial uses such as fire retardants. However, the effect of halogenated BPA and BPS derivatives on glucocorticoid metabolizing enzyme 11ß-hydroxysteroid dehydrogenase 1 (11ß-HSD1) remains unclear. The inhibitory effects of 6 BPA derivatives in the inhibition of human and rat 11ß-HSD1 were investigated. The potencies for inhibition on human 11ß-HSD1 were TBBPA (IC50, 3.87 µM) = monochloro BPA (MCBPA, 4.08 µM) = trichloro BPA (TrCBPA, 4.41 µM) > tetrachloro BPA (TCBPA, 9.75 µM) > TBBPS (>100 µM) = BPS (>100 µM), and those for rat 11ß-HSD1 were TrCBPA (IC50, 2.76 µM) = MCBPA (3.75 µM) > TBBPA (39.58 µM) > TCBPA = TBBPS = BPS. All these BPA derivatives are mixed/competitive inhibitors of both human and rat enzymes. Molecular docking studies predict that MCBPA, TrCBPA, TCBPA, and TBBPA all bind to the active site of human 11ß-HSD1, forming hydrogen bonds with catalytic residue Ser170 except TCBPA. Regression of the lowest binding energy with IC50 values revealed a significant inverse linear regression. In conclusion, halogenated BPA derivatives are mostly potent inhibitors of human and rat 11ß-HSD1, and there is structure-dependent inhibition.

Discovery of novel and potent InhA direct inhibitors by ensemble docking-based virtual screening and biological assays.

Multidrug-resistant tuberculosis (MDR-TB) continues to spread worldwide and remains one of the leading causes of death among infectious diseases. The enoyl-acyl carrier protein reductase (InhA) belongs to FAS-II family and is essential for the formation of the Mycobacterium tuberculosis cell wall. Recent years, InhA direct inhibitors have been extensively studied to overcome MDR-TB. However, there are still no inhibitors that have entered clinical research. Here, the ensemble docking-based virtual screening along with biological assay were used to identify potent InhA direct inhibitors from Chembridge, Chemdiv, and Specs. Ultimately, 34 compounds were purchased and first assayed for the binding affinity, of which four compounds can bind InhA well with KD values ranging from 48.4 to 56.2 µM. Among them, compound 9,222,034 has the best inhibitory activity against InhA enzyme with an IC50 value of 18.05 µM. In addition, the molecular dynamic simulation and binding free energy calculation indicate that the identified compounds bind to InhA with "extended" conformation. Residue energy decomposition shows that residues such as Tyr158, Met161, and Met191 have higher energy contributions in the binding of compounds. By analyzing the binding modes, we found that these compounds can bind to a hydrophobic sub-pocket formed by residues Tyr158, Phe149, Ile215, Leu218, etc., resulting in extensive van der Waals interactions. In summary, this study proposed an efficient strategy for discovering InhA direct inhibitors through ensemble docking-based virtual screening, and finally identified four active compounds with new skeletons, which can provide valuable information for the discovery and optimization of InhA direct inhibitors.

Chalcone derivatives as xanthine oxidase inhibitors: synthesis, binding mode investigation, biological evaluation, and ADMET prediction.

Xanthine oxidase (XO) is a crucial target for the treatment of hyperuricemia and gout. A series of derivatives based on natural 3,4-dihydroxychalcone, obtained from Carthamus tinctorious and Licorice, were designed and synthesized. Nine derivatives (9a-e, 10b,c, and 15a,b) exhibited apparent XO inhibitory activity in vitro (IC50 values varied from 0.121 to 7.086 µM), 15b presented the most potent inhibitory activity (IC50 = 0.121 µM), which was 27.47-fold higher than that of allopurinol (IC50 = 3.324 µM). The SAR analysis indicated that introducing hydroxyl groups at 3'/4'/5'-position on ring A was more beneficial to the inhibition of XO than at 2'/6'-position; the removal of 3­hydroxyl group on ring B could weaken the inhibitory potency of hydroxychalcones on XO, but it was beneficial to the XO inhibitory potency of methoxychalcones. Molecule modeling studies afforded insights into the binding mode of 15b with XO and supported the findings of SAR analysis. Additionally, kinetics studies demonstrated that 15b presented a reversible and competitive XO inhibitor, which spontaneously combined with XO through hydrophobic force, and finally changed the secondary conformation of XO. Furthermore, the acute hyperuricemia model was employed to investigate the hypouricemic effect of 15b, which could effectively reduce the serum uric acid levels of rats at an oral dose of 10 mg/kg. ADMET prediction suggested that compound 15b possessed good pharmacokinetic properties. Briefly, compound 15b emerges as an interesting XO inhibitor for the treatment of hyperuricemia and gout with beneficial effects on serum uric acid levels regulating. Meanwhile, the XO inhibitors with chalcone skeleton will deserve further attention and discussion.

Progress in programmed cell death-1/programmed cell death-ligand 1 pathway inhibitors and binding mode analysis.

Programmed cell death protein 1 (PD-1)/programmed cell death protein ligand 1 (PD-L1) plays an important role in negative regulating immunity. The search for effective PD-1/PD-L1 inhibitors has been at the cutting-edge of academic and industrial medicinal chemistry, leading to the emergence of 16 clinical candidate drugs and the launch of six monoclonal antibodies (mAbs) drugs. However, due to the unclear mechanism of the interaction between drugs and substances in vivo, the screening of preclinical drugs often takes a long time. In order to shorten the time of drug development as much as possible, the binding mode analysis that can simulate the interaction between drugs and substances in vivo at the molecular level can significantly shorten the drug development process. This paper reviews the mechanism of PD-1/PD-L1 signaling pathway at the molecular level, as well as the research progress and obstacles of inhibitors. Besides, we analyzed the binding mode of recently reported PD-1/PD-L1 inhibitors with PD-1 or PD-L1 protein in detail in order to provide ideas for the development of PD-1/PD-L1 inhibitors.

Theoretical Studies of Leu-Pro-Arg-Asp-Ala Pentapeptide (LPRDA) Binding to Sortase A of Staphylococcus aureus.

Sortase A (SrtA) of Staphylococcus aureus is a well-defined molecular target to combat the virulence of these clinically important bacteria. However up to now no efficient drugs or even clinical candidates are known, hence the search for such drugs is still relevant and necessary. SrtA is a complex target, so many straight-forward techniques for modeling using the structure-based drug design (SBDD) fail to produce the results they used to bring for other, simpler, targets. In this work we conduct theoretical studies of the binding/activity of Leu-Pro-Arg-Asp-Ala (LPRDA) polypeptide, which was recently shown to possess antivirulence activity against S. aureus. Our investigation was aimed at establishing a framework for the estimation of the key interactions and subsequent modification of LPRDA, targeted at non-peptide molecules, with better drug-like properties than the original polypeptide. Firstly, the available PDB structures are critically analyzed and the criteria to evaluate the quality of the ligand-SrtA complex geometry are proposed. Secondly, the docking protocol was investigated to establish its applicability to the LPRDA-SrtA complex prediction. Thirdly, the molecular dynamics studies were carried out to refine the geometries and estimate the stability of the complexes, predicted by docking. The main finding is that the previously reported partially chaotic movement of the ß6/ß7 and ß7/ß8 loops of SrtA (being the intrinsically disordered parts related to the SrtA binding site) is exaggerated when SrtA is complexed with LPRDA, which in turn reveals all the signs of the flexible and structurally disordered molecule. As a result, a wealth of plausible LPRDA-SrtA complex conformations are hard to distinguish using simple modeling means, such as docking. The use of more elaborate modeling approaches may help to model the system reliably but at the cost of computational efficiency.

Discovery of selective covalent cathepsin K inhibitors containing novel 4-cyanopyrimidine warhead based on quantum chemical calculations and binding mode analysis.

Cathepsin K (Cat K), mainly expressed by osteoclasts, plays an important role in bone resorption. Covalent Cat K inhibitors will show great potential in the future treatment of osteoporosis. It has been reported that the selectivity of covalent cathepsin K inhibitors was related to the drug's safety. The type of warhead has a crucial influence on the enzyme bioactivity and selectivity of covalent inhibitors. In order to develop novel covalent inhibitors with the selective new warhead, quantum chemical calculations were performed to estimate the reactivity of the nitrile warheads. Moreover, binding mode analysis between ligands and high homology Cat K, S and B revealed differences in non-covalent interactions. Novel covalent Cat K inhibitors containing 4-cyanopyrimidine warhead (11) were determined for the first time. Among them, compound 34 significantly inhibited Cat K (IC50 = 61.9 nM) with excellent selectivity compared to Cat S (>810-fold) and Cat B (>1620-fold), respectively. Binding mode analysis of Cat K-34 complex provided the basis for further optimization. Compound 34 could be a valuable lead compound for further research on safe and effective Cat K inhibitors.

Targeting homologous recombination (HR) repair mechanism for cancer treatment: discovery of new potential UCHL-3 inhibitors via virtual screening, molecular dynamics and binding mode analysis.

UCHL3 (ubiquitin C-terminal hydrolase-L3) is a de-ubiquitinating enzyme involved in the homologous recombination repair mechanism of double-strand breaks (DBS) of the DNA. Multiple studies indicated that UCHL3 inhibitors could be used in combination therapy with high therapeutic efficacy against cancer thus highlighting the validity of directing research against UCHL3 as a druggable target in oncology. In this study, a combination of virtual screening methods was utilized to identify new potential UCHL3 inhibitors. A series of UCHL3 ligands were identified by applying a combination of cheminformatics and molecular modeling filtration techniques to a ChemBl database of over two million small molecules viz. Lipinski's Rule of Five, Veber's rule, pharmacophore model, Hierarchical molecular docking, Pan-assay Interference Compounds (PAINS) alerts, toxicity filter, and single-point Molecular mechanics Poisson/Boltzmann surface area (MM/PBSA) docking pose rescoring. This multi-layer filtration strategy led to the identification of twenty-one compounds as potential UCHL3 inhibitors that were subsequently subjected to a 50 ns molecular dynamics (MD) simulations predict the stability of their ligand-protein complexes. Furthermore, MM/PBSA calculations based on MD trajectories were performed, and the energy contribution per residue to the binding energy was calculated. Three compounds, 1, 2 and 3, were finally recognized as having the highest potential of being UCHL3 inhibitors. Therefore, those were used for binding mode analysis to the UCHL3 active site, leading to identification of four residues as key for binding viz. Pro8, Leu55, Val166, and Leu168.Communicated by Ramaswamy H. Sarma.

Binding mode characterization of 13b in the monomeric and dimeric states of SARS-CoV-2 main protease using molecular dynamics simulations.

The main protease, Mpro/3CLpro, plays an essential role in processing polyproteins translated from viral RNA to produce functional viral proteins and therefore serve as an attractive target for discovering COVID-19 therapeutics. The availability of both monomer and dimer crystal bound with a common ligand, '13b' (α-ketoamide inhibitor), opened up opportunities to understand the Mpro mechanism of action. A comparative analysis of both forms of Mpro was carried out to elucidate the binding site architectural differences in the presence and absence of '13b'. Molecular dynamics simulations suggest that the presence of '13b' enhances the stability of Mpro than the unbound APO form. The N- and C- terminals of both the protomers stabilize each other, and making it's interface essential for the active form of Mpro. In comparison to monomer, the relatively high affinity of '13b' is gained in dimer pocket due to the high stability of the pocket by the interaction of S1 residue of chain B with residues F140, E166 and H172 of chain A, which is absent in monomer. The comprehensive essential dynamics, protein structure network analysis and thermodynamic profiling highlight the hot-spots, pivotal in molecular recognition process at protein-ligand and protein-protein interaction levels, cross-validated through computational alanine scanning study. A comparative description of '13b' binding mechanism in both forms illustrates valuable insights into the inhibition mechanism and the selection of critical residues suitable for the structure-based approaches for the identification of more potent Mpro inhibitors.Communicated by Ramaswamy H. Sarma.

Molecular Determinants and Pharmacological Analysis for a Class of Competitive Non-transported Bicyclic Inhibitors of the Betaine/GABA Transporter BGT1.

The betaine/GABA transporter 1 (BGT1) is a member of the GABA transporter (GAT) family with still elusive function, largely due to a lack of potent and selective tool compounds. Based on modeling, we here present the design, synthesis and pharmacological evaluation of five novel conformationally restricted cyclic GABA analogs related to the previously reported highly potent and selective BGT1 inhibitor (1S,2S,5R)-5-aminobicyclo[3.1.0]hexane-2-carboxylic acid (bicyclo-GABA). Using [3H]GABA radioligand uptake assays at the four human GATs recombinantly expressed in mammalian cell lines, we identified bicyclo-GABA and its N-methylated analog (2) as the most potent and selective BGT1 inhibitors. Additional pharmacological characterization in a fluorescence-based membrane potential assay showed that bicyclo-GABA and 2 are competitive inhibitors, not substrates, at BGT1, which was validated by a Schild analysis for bicyclo-GABA (pK B value of 6.4). To further elaborate on the selectivity profile both compounds were tested at recombinant α1ß2γ2 GABAA receptors. Whereas bicyclo-GABA showed low micromolar agonistic activity, the N-methylated 2 was completely devoid of activity at GABAA receptors. To further reveal the binding mode of bicyclo-GABA and 2 binding hypotheses of the compounds were obtained from in silico-guided mutagenesis studies followed by pharmacological evaluation at selected BGT1 mutants. This identified the non-conserved BGT1 residues Q299 and E52 as the molecular determinants driving BGT1 activity and selectivity. The binding mode of bicyclo-GABA was further validated by the introduction of activity into the corresponding GAT3 mutant L314Q (38 times potency increase cf. wildtype). Altogether, our data reveal the molecular determinants for the activity of bicyclic GABA analogs, that despite their small size act as competitive inhibitors of BGT1. These compounds may serve as valuable tools to selectively and potently target BGT1 in order to decipher its elusive pharmacological role in the brain and periphery such as the liver and kidneys.

Two Methods, One Goal: Structural Differences between Cocrystallization and Crystal Soaking to Discover Ligand Binding Poses.

In lead optimization, protein crystallography is an indispensable tool to analyze drug binding. Binding modes and non-covalent interaction inventories are essential to design follow-up synthesis candidates. Two protocols are commonly applied to produce protein-ligand complexes: cocrystallization and soaking. Because of its time and cost effectiveness, soaking is the more popular method. Taking eight ligand hinge binders of protein kinase A, we demonstrate that cocrystallization is superior. Particularly for flexible proteins, such as kinases, and larger ligands cocrystallization captures more reliable the correct binding pose and induced protein adaptations. The geometrical discrepancies between soaking and cocrystallization appear smaller for fragment-sized ligands. For larger flexible ligands that trigger conformational changes of the protein, soaking can be misleading and underestimates the number of possible polar interactions due to inadequate, highly impaired positions of protein amino-acid side and main chain atoms. Thus, if applicable cocrystallization should be the gold standard to study protein-ligand complexes.

Structural evidence for pheromone discrimination by the pheromone binding protein 3 from Plutella xylostella.

Insect pheromone binding proteins (PBPs) are believed to have a high degree of pheromone selectivity, acting as the first filter to discriminate specific pheromones from other volatile compounds. Herein, we provide evidence using homology-based model for the pheromone discrimination of Plutella xylostella pheromone binding protein 3 (PxPBP3). Combining molecular dynamics simulations and in vitro binding assays, two dominant sites are determined to be essential for the PxPBP3 to discriminate (Z)-11-hexadecenyl acetate (Hexadecenyl) from (Z)-11-hexadecenal (Hexadecenal). As the first key site for pheromone discrimination, Arg111 is indispensable to the PxPBP3-Hexadecenyl interaction. However, its importance in the binding of Hexadecenal to PxPBP3 is greatly reduced. A second site where pheromone discrimination occurs is a small loop (residues 34-38) in PxPBP3. It is shown that the hydrophobic strength provided by three hydrophobic residues (Phe34, Tyr37, and Trp38) in the small loop is significantly biased in the two complexes formed by PxPBP3 and the two pheromones. The discrimination capacity of PxPBP3 indicates that the P. xylostella pheromones may not share the same peri-receptor pathway, although they both show high affinity to PxPBP3.

Structure guided design of potent indole-based ATX inhibitors bearing hydrazone moiety with tumor suppression effects.

ATX was capable of catalyzing the hydrolysis of LPC to the lipid mediator LPA which attracted considerable attention on the development of potent ATX inhibitors. Herein, driven by the HTS product indole-based lead 1, a hybridization strategy was utilized to construct the trifluoroacetyl hydrazone moiety through assembling the phenyl thiazole fragment to the indole skeleton of lead 1. After a systematic structure guided optimization, by cycling the phenyl thiazole to the compacted benzothiazole or decreasing the lipophilicity, two promising ATX inhibitors (9j and 25a) were identified with IC50 values of 2.1 nM and 19.0 nM, respectively. All compounds were tested a panel of cancer cell lines and a preliminary affinity on breast cancer cell lines (SI > 16.5) were observed which shed a light on their potential application of breast cancer relevant cases. Through a dedicated docking study, the intramolecular pseudo-ring within the trifluoroacetylhydrazone moiety played a significant role in constraining the binding poses of 9j and 25a. Finally, a binding free energy calculation was conducted to examine the contribution of different interactions in binding affinity.