A Protein Data Bank survey reveals shortening of intermolecular hydrogen bonds in ligand-protein complexes when a halogenated ligand is an H-bond donor

Revealing Propolis Potential Activity on Inhibiting Estrogen Receptor and Heat Shock Protein 90 Overexpressed in Breast Cancer by Bioinformatics Approaches

Breast cancer is the most commonly diagnosed cancer globally, with the highest incidence of breast cancer occurring in Asian countries including Indonesia. Among the types of breast cancer, the estrogen receptor (ER)-positive subtype which is prominent with estrogen receptor alpha (ERα) and heat shock protein 90 (HSP90) overexpression genes becomes the most prevalent than the others, approximately 75% of all breast cancer cases. ERα and HSP90 play a role in breast cancer activities including breast tumor growth, invasion, and metastasis mechanism. Propolis, a natural bee product, has been explored for its anticancer activity. However, there is lack of studies that evaluated the potential inhibitor from propolis compounds to the ERα and HSP90 proteins. Therefore, this article focuses on examining the correlation between ERα and HSP90's role in breast cancer and investigating the potential of 93 unique propolis compositions in inhibiting these genes in breast cancer using in silico approaches. This study revealed the positive correlation between ERα and HSP90 genes in breast cancer disease development. Furthermore, we also found novel potential bioactive compounds of propolis against breast cancer through binding with ERα and HSP90; they were 3',4',7-trihydroxyisoflavone and baicalein-7-O-β-D glucopyranoside, respectively. Further research on these compounds is needed to elucidate deeper mechanisms and activity in the real biological system to develop new breast cancer drug treatments.

Preferred microenvironments of halogen bonds and hydrogen bonds revealed using statistics and QM/MM calculation studies

Hydrogen bonds (HBs) and halogen bonds (XBs) are two essential non-covalent interactions for molecular recognition and drug design. As proteins are heterogeneous in structure, the microenvironments of protein structures should have effects on the formation of HBs and XBs with ligands. However, there are no systematic studies reported on this effect to date. For quantitatively describing protein microenvironments, we defined the local hydrophobicities (LHs) and local dielectric constants (LDCs) in this study. With the defined parameters, we conducted an elaborate database survey on the basis of 22 011 ligand-protein structures to explore the microenvironmental preference of HBs (91 966 in total) and XBs (1436 in total). The statistics show that XBs prefer hydrophobic microenvironments compared to HBs. The polar residues like ASP are more likely to form HBs with ligands, while nonpolar residues such as PHE and MET prefer XBs. Both the LHs and LDCs (10.69 ± 4.36 for HBs; 8.86 ± 4.00 for XBs) demonstrate that XBs are prone to hydrophobic microenvironments compared with HBs with significant differences (p < 0.001), indicating that evaluating their strengths in the corresponding environments should be necessary. Quantum Mechanics-Molecular Mechanics (QM/MM) calculations reveal that in comparison with vacuum environments, the interaction energies of HBs and XBs are decreased to varying degrees given different microenvironments. In addition, the strengths of HBs are impaired more than those of XBs when the local dielectric constant's difference between the XB microenvironments and the HB microenvironments is large.

Structure-Function Study of a Novel Inhibitor of Cyclin-Dependent Kinase C in Arabidopsis

The circadian clock, an internal time-keeping system with a period of about 24 h, coordinates many physiological processes with the day-night cycle. We previously demonstrated that BML-259 [N-(5-isopropyl-2-thiazolyl) phenylacetamide], a small molecule with mammal CYCLIN DEPENDENT KINASE 5 (CDK5)/CDK2 inhibition activity, lengthens Arabidopsis thaliana (Arabidopsis) circadian clock periods. BML-259 inhibits Arabidopsis CDKC kinase, which phosphorylates RNA polymerase II in the general transcriptional machinery. To accelerate our understanding of the inhibitory mechanism of BML-259 on CDKC, we performed structure-function studies of BML-259 using circadian period-lengthening activity as an estimation of CDKC inhibitor activity in vivo. The presence of a thiazole ring is essential for period-lengthening activity, whereas acetamide, isopropyl and phenyl groups can be modified without effect. BML-259 analog TT-539, a known mammal CDK5 inhibitor, did not lengthen the period nor did it inhibit Pol II phosphorylation. TT-361, an analog having a thiophenyl ring instead of a phenyl ring, possesses stronger period-lengthening activity and CDKC;2 inhibitory activity than BML-259. In silico ensemble docking calculations using Arabidopsis CDKC;2 obtained by a homology modeling indicated that the different binding conformations between these molecules and CDKC;2 explain the divergent activities of TT539 and TT361.

Competition between electrostatic interactions and halogen bonding in the protein-ligand system: structural and thermodynamic studies of 5,6-dibromobenzotriazole-hCK2α complexes

CK2 is a member of the CMGC group of eukaryotic protein kinases and a cancer drug target. It can be efficiently inhibited by halogenated benzotriazoles and benzimidazoles. Depending on the scaffold, substitution pattern, and pH, these compounds are either neutral or anionic. Their binding poses are dictated by a hydrophobic effect (desolvation) and a tug of war between a salt bridge/hydrogen bond (to K68) and halogen bonding (to E114 and V116 backbone oxygens). Here, we test the idea that binding poses might be controllable by pH for ligands with near-neutral pK_a, using the conditionally anionic 5,6-DBBt and constitutively anionic TBBt as our models. We characterize the binding by low-volume Differential Scanning Fluorimetry (nanoDSF), Isothermal Calorimetry (ITC), Hydrogen/Deuterium eXchange (HDX), and X-ray crystallography (MX). The data indicate that the ligand pose away from the hinge dominates for the entire tested pH range (5.5-8.5). The insensitivity of the binding mode to pH is attributed to the perturbation of ligand pK_a upon binding that keeps it anionic in the ligand binding pocket at all tested pH values. However, a minor population of the ligand, detectable only by HDX, shifts towards the hinge in acidic conditions. Our findings demonstrate that electrostatic (ionic) interactions predominate over halogen bonding.

Leveraging structural and 2D-QSAR to investigate the role of functional group substitutions, conserved surface residues and desolvation in triggering the small molecule-induced dimerization of hPD-L1

Small molecules are rising as a new generation of immune checkpoints’ inhibitors, with compounds targeting the human Programmed death-ligand 1 (hPD-L1) protein are pioneering this area of research. Promising examples include the recently disclosed compounds from Bristol-Myers-Squibb (BMS). These molecules bind specifically to hPD-L1 through a unique mode of action. They induce dimerization between two hPD-L1 monomers through the hPD-1 binding interface in each monomer, thereby inhibiting the PD-1/PD-L1 axis. While the recently reported crystal structures of such small molecules bound to hPD-L1 reveal valuable insights regarding their molecular interactions, there is still limited information about the dynamics driving this unusual complex formation. The current study provides an in-depth computational structural analysis to study the interactions of five small molecule compounds in complex with hPD-L1. By employing a combination of molecular dynamic simulations, binding energy calculations and computational solvent mapping techniques, our analyses quantified the dynamic roles of different hydrophilic and lipophilic residues at the surface of hPD-L1 in mediating these interactions. Furthermore, ligand-based analyses, including Free-Wilson 2D-QSAR was conducted to quantify the impact of R-group substitutions at different sites of the phenoxy-methyl biphenyl core. Our results emphasize the importance of a terminal phenyl ring that must be present in any hPD-L1 small molecule inhibitor. This phenyl moiety overlaps with a very unfavorable hydration site, which can explain the ability of such small molecules to trigger hPD-L1 dimerization.

Synthesis of Novel Halogenated Heterocycles Based on o-Phenylenediamine and Their Interactions with the Catalytic Subunit of Protein Kinase CK2

Protein kinase CK2 is a highly pleiotropic protein kinase capable of phosphorylating hundreds of protein substrates. It is involved in numerous cellular functions, including cell viability, apoptosis, cell proliferation and survival, angiogenesis, or ER-stress response. As CK2 activity is found perturbed in many pathological states, including cancers, it becomes an attractive target for the pharma. A large number of low-mass ATP-competitive inhibitors have already been developed, the majority of them halogenated. We tested the binding of six series of halogenated heterocyclic ligands derived from the commercially available 4,5-dihalo-benzene-1,2-diamines. These ligand series were selected to enable the separation of the scaffold effect from the hydrophobic interactions attributed directly to the presence of halogen atoms. In silico molecular docking was initially applied to test the capability of each ligand for binding at the ATP-binding site of CK2. HPLC-derived ligand hydrophobicity data are compared with the binding affinity assessed by low-volume differential scanning fluorimetry (nanoDSF). We identified three promising ligand scaffolds, two of which have not yet been described as CK2 inhibitors but may lead to potent CK2 kinase inhibitors. The inhibitory activity against CK2α and toxicity against four reference cell lines have been determined for eight compounds identified as the most promising in nanoDSF assay.

Statistical physics modeling and interpretation of the adsorption of enantiomeric terpenes onto the human olfactory receptor OR1A1

The statistical physics approach has been well studied by our research team for liquid and gaseous adsorption systems. This treatment is based on the grand canonical partition function to give new interpretations of the adsorption process at molecular level for chemical senses: olfaction and taste. This work represents a contribution to understand the olfaction mechanism of four of enantiomeric terpenes by applying a statistical physics treatment that allows giving a physico-chemical meaning to parameters involved in the analytical model. It is possible to estimate the number of adsorbed molecules per site, the anchorage number, the receptor density, the concentration at half saturation and the molar adsorption energy. Through this selection of the best fitting model and through fitted values of these parameters, we showed that the adsorption of carvone and limonene enantiomers is not a multilayer process but a monolayer monosite process (monolayer adsorption model with identical and independent sites (n ≠ 1)). The physico-chemical model parameters can be used for the energetic characterization of the interactions between the carvone and the limonene enantiomers and the human olfactory receptor OR1A1 and the determination of an olfactory band of order of 14 kJ/mol, 7 kJ/mol, 9 kJ/mol, 8 kJ/mol for (R)-(-)-carvone, (S)-(+)-carvone, (R)-(+)-limonene and (S)-(-)-limonene, respectively, through the determination of the adsorption energy values and the adsorption energy distributions (AEDs). Thanks to the grand canonical formalism in statistical physics, the negative values of the Gibbs free enthalpy indicate that the adsorption process of the four enantiomeric terpenes onto the human olfactory receptor OR1A1 was spontaneous. The exothermic adsorption mechanism involved in the olfactory perception was explained via the negative values of the internal energy.

Mitigate the cytokine storm due to the severe COVID-19: A computational investigation of possible allosteric inhibitory actions on IL-6R and IL-1R using selected phytochemicals

The novel corona virus 2019 (COVID 19) is growing at an increasing rate with high mortality. Meanwhile, the cytokine storm is the most dangerous and potentially life-threatening event related to COVID 19. Phyto-compounds found in existing Ayurveda drugs have the ability to inhibit the Interleukin 6 (IL-6R) and Interleukin 1 (IL-1R) receptors. IL-6R and IL-1R receptors involve in cytokine storm and recognition of phytochemicals with proven safety profiles could open a pathway to the development of the most effective drugs against cytokine storm. In this study, we intend to perform an in silico investigation of effective phyto compounds, which can be isolated from selected medicinal herbs to avoid cytokine storm, inhibiting the IL-6 and IL-1 receptor binding process. An extensive literature survey followed by virtual screening was carried out to identify phytochemicals with potential anti-hyper-inflammatory action. Flexible docking was conducted for validated models of IL-1R and IL-6R-α with the most promising phytochemicals at possible allosteric sites using AutoDock Vina. Molecular dynamics (MD) studies were conducted for selected protein-ligand complexes using LARMD server and conformational changes were evaluated. According to the results, taepeenin J had Gibbs energy (ΔG) of -10.85 kcal/mol towards IL-1R but had limited oral bioavailability. MD analysis revealed that taepeenin J can cause significant conformational movements in IL-1R. Nortaepeenin B showed a ΔG of -8.5 kcal/mol towards IL-6R-α with an excellent oral bioavailability. MD analysis predicted that it can cause significant conformational movements in IL-6R-α. Hence, the evaluated phytochemicals are potential candidates for further in vitro studies for the development of medicine against cytokine storm on behalf of SARS-COV-2 infected patients.

A competition between hydrophobic and electrostatic interactions in protein-ligand systems. Binding of heterogeneously halogenated benzotriazoles by the catalytic subunit of human protein kinase CK2

A series of chlorine-substituted benzotriazole derivatives, representing all possible substitution patterns of halogen atoms attached to the benzotriazole benzene ring, were synthetized as potential inhibitors of human protein kinase CK2. Basic ADME parameters for the free solutes (hydrophobicity, electronic properties) together with their binding affinity to the catalytic subunit of protein kinase CK2 were determined with reverse-phase HPLC, spectrophotometric titration, and Thermal Shift Assay Method, respectively. The analysis of position-dependent thermodynamic contribution of a chlorine atom attached to the benzotriazole ring confirmed the previous observation for brominated benzotriazoles, in which substitution at positions 5 and 6 with bromine was found crucial for ligand binding. In all tested halogenated benzotriazoles the replacement of Br with Cl decreases the hydrophobicity, while the electronic properties remain virtually unaffected. Supramolecular architecture identified in the just resolved crystal structures of three of the four possible dichloro-benzotriazoles shows how substitution distant from the triazole ring affects the pattern of intermolecular interactions. Summarizing, the benzotriazole benzene ring substitution pattern has been identified as the main driver of ligand binding, predominating the non-specific hydrophobic effect.

Halogen Bonding in Nucleic Acid Complexes

Halogen bonding (X-bonding) has attracted notable attention among noncovalent interactions. This highly directional attraction between a halogen atom and an electron donor has been exploited in knowledge-based drug design. A great deal of information has been gathered about X-bonds in protein-ligand complexes, as opposed to nucleic acid complexes. Here we provide a thorough analysis of nucleic acid complexes containing either halogenated building blocks or halogenated ligands. We analyzed close contacts between halogens and electron-rich moieties. The phosphate backbone oxygen is clearly the most common halogen acceptor. We identified 21 X-bonds within known structures of nucleic acid complexes. A vast majority of the X-bonds is formed by halogenated nucleobases, such as bromouridine, and feature excellent geometries. Noncovalent ligands have been found to form only interactions with suboptimal interaction geometries. Hence, the first X-bonded nucleic acid binder remains to be discovered.

Hydration effects on the efficacy of the Epidermal growth factor receptor kinase inhibitor afatinib

Small molecules targeting the EGFR tyrosine kinase domain have been used with some success at treating patients with non-small cell lung cancer driven by activating mutations in the kinase domain. The initial class of inhibitors displaced ATP noncovalently but were rendered ineffective due to the development of resistance mutations in the kinase domain. These were overcome by the development of covalent inhibitors such as afatinib which also bind in the ATP pocket. However pooled analysis of two recent clinical trials LUX-3 and LUX-6 demonstrated an unprecedented overall survival benefit of afatinib over chemotherapy for the EGFR^19del, but not the EGFR^L858R. In the current study we use modelling and simulations to show that structural constraints in EGFR^19del deletion result in significantly attenuated flexibilities in the binding pocket resulting in strong hydrogen and halogen bonds with afatinib in the EGFR^19del; these constraints are modulated by buried water and result in the differential affinities of afatinib for the different mutants. SNP analysis of residues surrounding the buried water points to the likelihood of further differential effects of afatinib and provides a compelling case for investigating the effects of the SNPs towards further stratification of patients for ensuring the most effective use of afatinib.

Thermodynamics parameters for binding of halogenated benzotriazole inhibitors of human protein kinase CK2α

The interaction of human CK2α (hCK2α) with nine halogenated benzotriazoles, TBBt and its analogues representing all possible patterns of halogenation on the benzene ring of benzotriazole, was studied by biophysical methods. Thermal stability of protein-ligand complexes, monitored by calorimetric (DSC) and optical (DSF) methods, showed that the increase in the mid-point temperature for unfolding of protein-ligand complexes (i.e. potency of ligand binding to hCK2α) follow the inhibitory activities determined by biochemical assays. The dissociation constant for the ATP-hCK2α complex was estimated with the aid of microscale thermophoresis (MST) as 4.3±1.8 μM, and MST-derived dissociation constants determined for halogenated benzotriazoles, when converted according to known ATP concentrations, perfectly reconstruct IC50 values determined by the biochemical assays. Ligand-dependent quenching of tyrosine fluorescence, together with molecular modeling and DSC-derived heats of unfolding, support the hypothesis that halogenated benzotriazoles bind in at least two alternative orientations, and those that are efficient hCK2α inhibitors bind in the orientation which TBBt adopts in its complex with maize CK2α. DSC-derived apparent heat for ligand binding (ΔΔHbind) is driven by intermolecular electrostatic interactions between Lys68 and the triazole ring of the ligand, as indicated by a good correlation between ΔΔHbind and ligand pKa. Overall results, additionally supported by molecular modeling, confirm that a balance of hydrophobic and electrostatic interactions contribute predominantly (~40 kJ/mol), relative to possible intermolecular halogen/hydrogen bonding (less than 10 kJ/mol), in binding of halogenated benzotriazoles to the ATP-binding site of hCK2α. This article is part of a Special Issue entitled: Inhibitors of Protein Kinases.

Thermodynamic parameters for binding of some halogenated inhibitors of human protein kinase CK2

The interaction of human CK2α with a series of tetrabromobenzotriazole (TBBt) and tetrabromobenzimidazole (TBBz) analogs, in which one of the bromine atoms proximal to the triazole/imidazole ring is replaced by a methyl group, was studied by biochemical (IC50) and biophysical methods (thermal stability of protein-ligand complex monitored by DSC and fluorescence). Two newly synthesized tri-bromo derivatives display inhibitory activity comparable to that of the reference compounds, TBBt and TBBz, respectively. DSC analysis of the stability of protein-ligand complexes shows that the heat of ligand binding (Hbind) is driven by intermolecular electrostatic interactions involving the triazole/imidazole ring, as indicated by a strong correlation between Hbind and ligand pKa. Screening, based on fluorescence-monitored thermal unfolding of protein-ligand complexes, gave comparable results, clearly identifying ligands that most strongly bind to the protein. Overall results, additionally supported by molecular modeling, confirm that a balance of hydrophobic and electrostatic interactions contribute predominantly, relative to possible intermolecular halogen bonding, in binding of the ligands to the CK2α ATP-binding site.