A conserved water-mediated hydrogen bond network defines bosutinib's kinase selectivity

SCN as a local probe of protein structural dynamics

The dynamics of lysozyme is probed by attaching -SCN to all alanine residues. The one-dimensional infrared spectra exhibit frequency shifts in the position of the maximum absorption of 4 cm-1, which is consistent with experiments in different solvents and indicates moderately strong interactions of the vibrational probe with its environment. Isotopic substitution 12C → 13C leads to a redshift by -47 cm-1, which agrees quantitatively with experiments for CN-substituted copper complexes in solution. The low-frequency, far-infrared part of the protein spectra contains label-specific information in the difference spectra when compared with the wild type protein. Depending on the position of the labels, local structural changes are observed. For example, introducing the -SCN label at Ala129 leads to breaking of the α-helical structure with concomitant change in the far-infrared spectrum. Finally, changes in the local hydration of SCN-labeled alanine residues as a function of time can be related to the reorientation of the label. It is concluded that -SCN is potentially useful for probing protein dynamics, both in the high-frequency part (CN-stretch) and in the far-infrared part of the spectrum.

Under the Influence of Water: Molecular Recognition of Organic Hydrophilic Molecules in Water with a Prismarene Host Driven by Hydration Effects

A water-soluble prism[5]arene host can form endo-cavity complexes with hydrophilic organic substances in water by displacing frustrated water molecules from its deep cavity. Water molecules structured at both rims of the prismarene host can mediate hydrogen bonding interactions with the guest. Water-mediated hydrogen bonding interactions were invoked here to elucidate the elevated binding affinities and selectivity of the prismarene host toward hydrophilic organic guests. We show that water at the interface of a host-guest complex can act as an extension of the host structure, facilitating the accommodation of neutral guests within the binding site. This study highlights the crucial role of water in facilitating supramolecular interactions between a deep-cavity prismarene host and organic hydrophilic guests in aqueous medium.

Network pharmacology- and cell-based assessments identify the FAK/Src pathway as a molecular target for the antimetastatic effect of momordin Ic against cholangiocarcinoma

Previous studies have indicated the efficacy of momordin Ic (MIc), a plant-derived triterpenoid, against several types of cancers, implying its potential for further development. However, comprehensive insights into the molecular mechanisms and targets of MIc in cholangiocarcinoma (CCA) are lacking. This study aimed to investigate the actions of MIc against CCA at the molecular level. Network pharmacology analysis was first employed to predict the mechanisms and targets of MIc. The results unveiled the potential involvement of MIc in apoptosis and cell migration, pinpointing Src and FAK as key targets. Subsequently, cell-based assays, in accordance with FAK/Src-associated metastasis, were conducted, demonstrating the ability of MIc to attenuate the metastatic behaviours of KKU-452 cells. The in vitro results further indicated the capability of MIc to suppress the epithelial-mesenchymal transition (EMT) process, notably by downregulating EMT regulators, including N-cadherin, vimentin, ZEB2 and FOXC1/2 expression. Furthermore, MIc suppressed the activation of the FAK/Src signalling pathway, influencing critical downstream factors such as MMP-9, VEGF, ICAM-1, and c-Myc. Molecular docking simulations also suggested that MIc could interact with FAK and Src domains and restrain kinases from being activated by hindering ATP binding. In conclusion, this study employs a comprehensive approach encompassing network pharmacology analysis, in vitro assays, and molecular docking to unveil the mechanisms and targets of MIc in CCA. MIc mitigates metastatic behaviours and suppresses key pathways, offering a promising avenue for future therapeutic strategies against this aggressive cancer.

The structures of salt-inducible kinase 3 in complex with inhibitors reveal determinants for binding and selectivity

The salt-inducible kinases (SIKs) 1 to 3, belonging to the AMPK-related kinase family, serve as master regulators orchestrating a diverse set of physiological processes such as metabolism, bone formation, immune response, oncogenesis, and cardiac rhythm. Owing to its key regulatory role, the SIK kinases have emerged as compelling targets for pharmacological intervention across a diverse set of indications. Therefore, there is interest in developing SIK inhibitors with defined selectivity profiles both to further dissect the downstream biology and for treating disease. However, despite a large pharmaceutical interest in the SIKs, experimental structures of SIK kinases are scarce. This is likely due to the challenges associated with the generation of proteins suitable for structural studies. By adopting a rational approach to construct design and protein purification, we successfully crystallized and subsequently solved the structure of SIK3 in complex with HG-9-91-01, a potent SIK inhibitor. To enable further SIK3-inhibitor complex structures we identified an antibody fragment that facilitated crystallization and enabled a robust protocol suitable for structure-based drug design. The structures reveal SIK3 in an active conformation, where the ubiquitin-associated domain is shown to provide further stabilization to this active conformation. We present four pharmacologically relevant and distinct SIK3-inhibitor complexes. These detail the key interaction for each ligand and reveal how different regions of the ATP site are engaged by the different inhibitors to achieve high affinity. Notably, the structure of SIK3 in complex with a SIK3 specific inhibitor offers insights into isoform selectivity.

Discovery of a Covalent Inhibitor Selectively Targeting the Autophosphorylation Site of c-Src Kinase

Nonreceptor tyrosine kinase c-Src plays a crucial role in cell signaling and contributes to tumor progression. However, the development of selective c-Src inhibitors turns out to be challenging. In our previous study, we performed posttranslational modification-inspired drug design (PTMI-DD) to provide a plausible way for designing selective kinase inhibitors. In this study, after identifying a unique pocket comprising a less conserved cysteine and an autophosphorylation site in c-Src as well as a promiscuous covalent inhibitor, chemical optimization was performed to obtain (R)-LW-Srci-8 with nearly 75-fold improved potency (IC50 = 35.83 ± 7.21 nM). Crystallographic studies revealed the critical C-F···C═O interactions that may contribute to tight binding. The kinact and Ki values validated the improved binding affinity and decreased warhead reactivity of (R)-LW-Srci-8 for c-Src. Notably, in vitro tyrosine kinase profiling and cellular activity-based protein profiling (ABPP) cooperatively indicated a specific inhibition of c-Src by (R)-LW-Srci-8. Intriguingly, (R)-LW-Srci-8 preferentially binds to inactive c-Src with unphosphorylated Y419 both in vitro and in cells, subsequently disrupting the autophosphorylation. Collectively, our study demonstrated the feasibility of developing selective kinase inhibitors by cotargeting a nucleophilic residue and a posttranslational modification site and providing a chemical probe for c-Src functional studies.

Classifying protein kinase conformations with machine learning

Protein kinases are key actors of signaling networks and important drug targets. They cycle between active and inactive conformations, distinguished by a few elements within the catalytic domain. One is the activation loop, whose conserved DFG motif can occupy DFG-in, DFG-out, and some rarer conformations. Annotation and classification of the structural kinome are important, as different conformations can be targeted by different inhibitors and activators. Valuable resources exist; however, large-scale applications will benefit from increased automation and interpretability of structural annotation. Interpretable machine learning models are described for this purpose, based on ensembles of decision trees. To train them, a set of catalytic domain sequences and structures was collected, somewhat larger and more diverse than existing resources. The structures were clustered based on the DFG conformation and manually annotated. They were then used as training input. Two main models were constructed, which distinguished active/inactive and in/out/other DFG conformations. They considered initially 1692 structural variables, spanning the whole catalytic domain, then identified ("learned") a small subset that sufficed for accurate classification. The first model correctly labeled all but 3 of 3289 structures as active or inactive, while the second assigned the correct DFG label to all but 17 of 8826 structures. The most potent classifying variables were all related to well-known structural elements in or near the activation loop and their ranking gives insights into the conformational preferences. The models were used to automatically annotate 3850 kinase structures predicted recently with the Alphafold2 tool, showing that Alphafold2 reproduced the active/inactive but not the DFG-in proportions seen in the Protein Data Bank. We expect the models will be useful for understanding and engineering kinases.

Critical Evaluation of Polarizable and Nonpolarizable Force Fields for Proteins Using Experimentally Derived Nitrile Electric Fields

Molecular dynamics (MD) simulations are frequently carried out for proteins to investigate the role of electrostatics in their biological function. The choice of force field (FF) can significantly alter the MD results, as the simulated local electrostatic interactions lack benchmarking in the absence of appropriate experimental methods. We recently reported that the transition dipole moment (TDM) of the popular nitrile vibrational probe varies linearly with the environmental electric field, overcoming well-known hydrogen bonding (H-bonding) issues for the nitrile frequency and, thus, enabling the unambiguous measurement of electric fields in proteins (J. Am. Chem. Soc. 2022, 144 (17), 7562-7567). Herein, we utilize this new strategy to enable comparisons of experimental and simulated electric fields in protein environments. Specifically, previously determined TDM electric fields exerted onto nitrile-containing o-cyanophenylalanine residues in photoactive yellow protein are compared with MD electric fields from the fixed-charge AMBER FF and the polarizable AMOEBA FF. We observe that the electric field distributions for H-bonding nitriles are substantially affected by the choice of FF. As such, AMBER underestimates electric fields for nitriles experiencing moderate field strengths; in contrast, AMOEBA robustly recapitulates the TDM electric fields. The FF dependence of the electric fields can be partly explained by the presence of additional negative charge density along the nitrile bond axis in AMOEBA, which is due to the inclusion of higher-order multipole parameters; this, in turn, begets more head-on nitrile H-bonds. We conclude by discussing the implications of the FF dependence for the simulation of nitriles and proteins in general.

Molecular mechanisms and evolutionary robustness of a color switch in proteorhodopsins

Proteorhodopsins are widely distributed photoreceptors from marine bacteria. Their discovery revealed a high degree of evolutionary adaptation to ambient light, resulting in blue- and green-absorbing variants that correlate with a conserved glutamine/leucine at position 105. On the basis of an integrated approach combining sensitivity-enhanced solid-state nuclear magnetic resonance (ssNMR) spectroscopy and linear-scaling quantum mechanics/molecular mechanics (QM/MM) methods, this single residue is shown to be responsible for a variety of synergistically coupled structural and electrostatic changes along the retinal polyene chain, ionone ring, and within the binding pocket. They collectively explain the observed color shift. Furthermore, analysis of the differences in chemical shift between nuclei within the same residues in green and blue proteorhodopsins also reveals a correlation with the respective degree of conservation. Our data show that the highly conserved color change mainly affects other highly conserved residues, illustrating a high degree of robustness of the color phenotype to sequence variation.

Revealing curcumin therapeutic targets on SRC, PPARG, MAPK8 and HSP90 as liver cirrhosis therapy based on comprehensive bioinformatic study

Cirrhosis naturally progresses through three stages: compensated, decompensated, and late decompensated, which carry an elevated risk of death. Although curcumin's anti-cirrhosis effects have been studied, underlying mechanism in preventing cirrhosis progression and the correlation between curcumin's action with upregulated genes remains insufficiently explored. In this study, we employed network pharmacology approach to construct a drug-target-disease network through bioinformatics and validate the findings with molecular docking and dynamic simulation. The curcumin-targeted liver cirrhosis network encompassed 54 nodes with 282 edges in protein-protein interactions (PPI) network. By utilizing network centrality analysis, we identified eight crucial genes. KEGG enrichment pathway revealed that these crucial genes are involved in pathway of cancer, endocrine resistance, estrogen signaling, chemical carcinogenesis-receptor activation, lipid metabolism, and atherosclerosis. Notably, these eight genes predominantly participate in cancer-related pathways. Further investigation revealed upregulation of four genes and downregulation of four others in hepatocellular carcinoma patients. These upregulated genes-MAPK8, SRC, PPARG, and HSP90AA1-strongly correlated with reduced survival probability in liver hepatocellular carcinoma patients with survival times approximately under 4000 days (∼11 years). Molecular docking and molecular dynamic results exhibited curcumin's superior binding affinities and stability compared to native ligands of MAPK8, SRC, PPARG, and HSP90AA1 within 50 ns simulations. Moreover, MM-GBSA analysis showed stronger binding energy of curcumin to MAPK8, SRC, and HSP90AA1 than native ligand. In conclusion, this study provides valuable insights into curcumin's potential mechanisms in preventing liver cirrhosis progression, specifically in HCC. These findings offer a theoretical basis for further pharmacological research into anti-HCC effect of curcumin.Communicated by Ramaswamy H. Sarma.

Bond-selective fluorescence imaging with single-molecule sensitivity

Bioimaging harnessing optical contrasts and chemical specificity is of vital importance in probing complex biology. Vibrational spectroscopy based on mid-infrared (mid-IR) excitation can reveal rich chemical information about molecular distributions. However, its full potential for bioimaging is hindered by the achievable sensitivity. Here, we report bond selective fluorescence-detected infrared-excited (BonFIRE) spectral microscopy. BonFIRE employs two-photon excitation in the mid-IR and near-IR to upconvert vibrational excitations to electronic states for fluorescence detection, thus encoding vibrational information into fluorescence. The system utilizes tuneable narrowband picosecond pulses to ensure high sensitivity, biocompatibility, and robustness for bond-selective biological interrogations over a wide spectrum of reporter molecules. We demonstrate BonFIRE spectral imaging in both fingerprint and cell-silent spectroscopic windows with single-molecule sensitivity for common fluorescent dyes. We then demonstrate BonFIRE imaging on various intracellular targets in fixed and live cells, neurons, and tissues, with promises for further vibrational multiplexing. For dynamic bioanalysis in living systems, we implement a high-frequency modulation scheme and demonstrate time-lapse BonFIRE microscopy of live HeLa cells. We expect BonFIRE to expand the bioimaging toolbox by providing a new level of bond-specific vibrational information and facilitate functional imaging and sensing for biological investigations.

HINT, a code for understanding the interaction between biomolecules: a tribute to Donald J. Abraham

A long-lasting goal of computational biochemists, medicinal chemists, and structural biologists has been the development of tools capable of deciphering the molecule-molecule interaction code that produces a rich variety of complex biomolecular assemblies comprised of the many different simple and biological molecules of life: water, small metabolites, cofactors, substrates, proteins, DNAs, and RNAs. Software applications that can mimic the interactions amongst all of these species, taking account of the laws of thermodynamics, would help gain information for understanding qualitatively and quantitatively key determinants contributing to the energetics of the bimolecular recognition process. This, in turn, would allow the design of novel compounds that might bind at the intermolecular interface by either preventing or reinforcing the recognition. HINT, hydropathic interaction, was a model and software code developed from a deceptively simple idea of Donald Abraham with the close collaboration with Glen Kellogg at Virginia Commonwealth University. HINT is based on a function that scores atom-atom interaction using LogP, the partition coefficient of any molecule between two phases; here, the solvents are water that mimics the cytoplasm milieu and octanol that mimics the protein internal hydropathic environment. This review summarizes the results of the extensive and successful collaboration between Abraham and Kellogg at VCU and the group at the University of Parma for testing HINT in a variety of different biomolecular interactions, from proteins with ligands to proteins with DNA.

Accounting for Solvation Correlation Effects on the Thermodynamics of Water Networks in Protein Cavities

Macromolecular recognition and ligand binding are at the core of biological function and drug discovery efforts. Water molecules play a significant role in mediating the protein-ligand interaction, acting as more than just the surrounding medium by affecting the thermodynamics and thus the outcome of the binding process. As individual water contributions are impossible to measure experimentally, a range of computational methods have emerged to identify hydration sites in protein pockets and characterize their energetic contributions for drug discovery applications. Even though several methods model solvation effects explicitly, they focus on determining the stability of specific water sites independently and neglect solvation correlation effects upon replacement of clusters of water molecules, which typically happens in hit-to-lead optimization. In this work, we rigorously determine the conjoint effects of replacing all combinations of water molecules in protein binding pockets through the use of the RE-EDS multistate free-energy method, which combines Hamiltonian replica exchange (RE) and enveloping distribution sampling (EDS). Applications on the small bovine pancreatic trypsin inhibitor and four proteins of the bromodomain family illustrate the extent of solvation correlation effects on water thermodynamics, with the favorability of replacement of the water sites by pharmacophore probes highly dependent on the composition of the water network and the pocket environment. Given the ubiquity of water networks in biologically relevant protein targets, we believe our approach can be helpful for computer-aided drug discovery by providing a pocket-specific and a priori systematic consideration of solvation effects on ligand binding and selectivity.

In Silico Evaluation of the Thr58-Associated Conserved Water with KRAS Switch-II Pocket Binders

The KRAS switch-II pocket (SII-P) has proven to be one of the most successful tools for targeting KRAS with small molecules to date. This has been demonstrated with several KRAS(G12C)-targeting covalent inhibitors, already resulting in two FDA-approved drugs. Several earlier-stage compounds have also been reported to engage KRAS SII-P with other position 12 mutants, including G12D, G12S, and G12R. A highly conserved water molecule exists in the KRAS SII-P, linking Thr58 of switch-II and Gly10 of β1 sheet. This conserved water is also present in the cocrystal structures of most of the disclosed small-molecule inhibitors but is only displaced by a handful of SII-P binders. Here, we evaluated the conserved water molecule energetics by the WaterMap for the SII-P binders with publicly disclosed structures and studied the water behavior in the presence of selected inhibitors by microsecond timescale molecular dynamics (MD) simulations using two water models (total simulation time of 120 μs). Our data revealed the high-energy nature of this hydration site when coexisting with an SII-P binder and that there is a preference for a single isolated hydration site in this location within the most advanced compounds. Furthermore, water displacement was only achieved with a few disclosed compounds and was suboptimal, as for instance a cyanomethyl group as a water displacer appears to introduce repulsion with the native conformation of Thr58. These results suggested that this conserved water should be considered more central when designing new inhibitors, especially in the design of noncovalent inhibitors targeting the SII-P.

Development of a Covalent Inhibitor of c-Jun N-Terminal Protein Kinase (JNK) 2/3 with Selectivity over JNK1

The c-Jun N-terminal kinases (JNKs) are members of the mitogen-activated protein kinase (MAPK) family, which includes JNK1-JNK3. Interestingly, JNK1 and JNK2 show opposing functions, with JNK2 activity favoring cell survival and JNK1 stimulating apoptosis. Isoform-selective small molecule inhibitors of JNK1 or JNK2 would be useful as pharmacological probes but have been difficult to develop due to the similarity of their ATP binding pockets. Here, we describe the discovery of a covalent inhibitor YL5084, the first such inhibitor that displays selectivity for JNK2 over JNK1. We demonstrated that YL5084 forms a covalent bond with Cys116 of JNK2, exhibits a 20-fold higher Kinact/KI compared to that of JNK1, and engages JNK2 in cells. However, YL5084 exhibited JNK2-independent antiproliferative effects in multiple myeloma cells, suggesting the existence of additional targets relevant in this context. Thus, although not fully optimized, YL5084 represents a useful chemical starting point for the future development of JNK2-selective chemical probes.

Utilization of Supervised Machine Learning to Understand Kinase Inhibitor Toxophore Profiles

There have been more than 70 FDA-approved drugs to target the ATP binding site of kinases, mainly in the field of oncology. These compounds are usually developed to target specific kinases, but in practice, most of these drugs are multi-kinase inhibitors that leverage the conserved nature of the ATP pocket across multiple kinases to increase their clinical efficacy. To utilize kinase inhibitors in targeted therapy and outside of oncology, a narrower kinome profile and an understanding of the toxicity profile is imperative. This is essential when considering treating chronic diseases with kinase targets, including neurodegeneration and inflammation. This will require the exploration of inhibitor chemical space and an in-depth understanding of off-target interactions. We have developed an early pipeline toxicity screening platform that uses supervised machine learning (ML) to classify test compounds’ cell stress phenotypes relative to a training set of on-market and withdrawn drugs. Here, we apply it to better understand the toxophores of some literature kinase inhibitor scaffolds, looking specifically at a series of 4-anilinoquinoline and 4-anilinoquinazoline model libraries.

Photoenhancement of the C≡N Stretching Vibration Intensity of Aromatic Nitriles

The C≡N stretching vibration is a versatile infrared (IR) reporter that is useful for a wide range of applications. Aiming to further expand its spectroscopic utility, herein, we show that, using 4-cyanoindole and 4-cyano-7-azaindole as examples, photoexcitation can significantly shift the frequency (ν_CN) and enhance the molar extinction coefficient (ε_CN) of this vibrational mode of aromatic nitriles and that, for these indole derivatives, the enhancement factor can reach 13. Moreover, we find that while solvent relaxation at the excited electronic state(s) always leads to an increase in ε_CN, its effect on ν_CN depends on the solute and the solvent. Taken together, these results demonstrate that solvent relaxation can differently affect the local environment of the nitrile group and its conjugation with the indole ring and, more importantly, that the C≡N stretching vibration can serve as a sensitive IR probe of charge and electron transfer processes in which an aromatic nitrile is involved.

Isolation and anticancer activity evaluation of rare Bisaryl anthraquinone antibiotics from novel Streptomyces sp. strain of NW Himalayan region

Biosynthesis of bisaryl preanthraquinone antibiotics by various microorganisms differs in monomeric subunits as well as their dimerization positions leading to different configurations. The present study relates to the production of rare bisaryl anthraquinone antibiotics by a new Streptomyces strain isolated from Shivalik region of NW Himalayas. In vitro anticancer and anti-migratory effects of Setomimycin (9,9' bisanthraquinone antibiotic) was seen with a significant reduction in the expression of both MEK as well as ERK pathways in a dose dependent manner at 6.5 μM & 8 μM concentration in HCT-116 and 5.5 μM & 7 μM concentration in MCF-7 cells. In vivo studies in aggressive orthotopic mouse mammary carcinoma model (4T1) demonstrated about 76% reduction of primary tumor weight and 90.5% reduction in the tumor volume within two weeks. In vivo pharmacokinetics study of setomimycin revealed that it can be rapidly absorbed with an adequate plasma exposure and half-life which can be linked to its in vivo efficacy.

Intrinsically Fluorescent Anti-Cancer Drugs

At present, about one-third of the total protein targets in the pharmaceutical research sector are kinase-based. While kinases have been attractive targets to combat many diseases, including cancer, selective kinase inhibition has been challenging, because of the high degree of structural homology in the active site where many kinase inhibitors bind. Despite efficacy as cancer drugs, kinase inhibitors can exhibit limited target specificity and rationalizing their target profiles in the context of precise molecular mechanisms or rearrangements is a major challenge for the field. Spectroscopic approaches such as infrared, Raman, NMR and fluorescence have the potential to provide significant insights into drug-target and drug-non-target interactions because of sensitivity to molecular environment. This review places a spotlight on the significance of fluorescence for extracting information related to structural properties, discovery of hidden conformers in solution and in target-bound state, binding properties (e.g., location of binding sites, hydrogen-bonding, hydrophobicity), kinetics as well as dynamics of kinase inhibitors. It is concluded that the information gleaned from an understanding of the intrinsic fluorescence from these classes of drugs may aid in the development of future drugs with improved side-effects and less disease resistance.

Discovery of Potent and Orally Bioavailable Pyridine N-Oxide-Based Factor XIa Inhibitors through Exploiting Nonclassical Interactions

Activated factor XI (FXIa) inhibitors are promising novel anticoagulants with low bleeding risk compared with current anticoagulants. The discovery of potent FXIa inhibitors with good oral bioavailability has been challenging. Herein, we describe our discovery effort, utilizing nonclassical interactions to improve potency, cellular permeability, and oral bioavailability by enhancing the binding while reducing polar atoms. Beginning with literature-inspired pyridine N-oxide-based FXIa inhibitor 1, the imidazole linker was first replaced with a pyrazole moiety to establish a polar C-H···water hydrogen-bonding interaction. Then, structure-based drug design was employed to modify lead molecule 2d in the P1' and P2' regions, with substituents interacting with key residues through various nonclassical interactions. As a result, a potent FXIa inhibitor 3f (K_i = 0.17 nM) was discovered. This compound demonstrated oral bioavailability in preclinical species (rat 36.4%, dog 80.5%, and monkey 43.0%) and displayed a dose-dependent antithrombotic effect in a rabbit arteriovenous shunt model of thrombosis.

Efficient expression, purification, and visualization by cryo-EM of unliganded near full-length HER3

Biochemical analyses of membrane receptor kinases have been limited by challenges in obtaining sufficient homogeneous receptor samples for downstream structural and biophysical characterization. Here, we report a suite of methods for the efficient expression, purification, and visualization by cryo-electron microscopy (cryo-EM) of near full-length Human Epidermal Growth Factor Receptor 3 (HER3), a receptor tyrosine pseudokinase, in the unliganded state. Through transient mammalian cell expression, a two-step purification with detergent exchange into lauryl maltose neopentyl glycol (LMNG), and freezing devoid of background detergent micelle, we obtained ~6Å reconstructions of the ~60kDa fully-glycosylated unliganded extracellular domain of HER3 from just 30mL of suspension culture. The reconstructions reveal previously unappreciated extracellular domain dynamics and glycosylation sites.

Targeting the HER3 pseudokinase domain with small molecule inhibitors

HER3 is a potent oncogenic growth factor receptor belonging to the human epidermal growth factor (HER/EGFR) family of receptor tyrosine kinases. In contrast to other EGFR family members, HER3 is a pseudokinase, lacking functional kinase activity. As such, efforts to develop small molecule tyrosine kinase inhibitors against this family member have been limited. In response to HER3-specific growth factors such as neuregulin (NRG, also known as heregulin or HRG), HER3 must couple with catalytically active family members, including its preferred partner HER2. Dimerization of the intracellular HER2:HER3 kinase domains is a critical part of the activation mechanism and HER3 plays a specialized role as an allosteric activator of the active HER2 kinase partner. Intriguingly, many pseudokinases retain functionally important nucleotide binding capacity, despite loss of kinase activity. We demonstrated that occupation of the nucleotide pocket of the pseudokinase HER3 retains functional importance for growth factor signaling through oncogenic HER2:HER3 heterodimers. Mutation of the HER3 nucleotide pocket both disrupts signaling and disrupts HER2:HER3 dimerization. Conversely, ATP competitive drugs which bind to HER3, but not HER2, can stabilize HER2:HER3 dimers, induce signaling and promote cell growth in breast cancer models. This indicates a nucleotide-dependent conformational role for the HER3 kinase domain. Critically, our recent proof-of-concept work demonstrated that HER3-directed small molecule inhibitors can also disrupt HER2:HER3 dimerization and signaling, supporting the prospect that HER3 can be a direct drug target despite its lack of intrinsic activity. In this chapter we will describe methods for identifying and validating small molecule inhibitors against the HER3 pseudokinase.

Molecular Dynamics Simulations of the Conformational Plasticity in the Active Pocket of Salt-Inducible Kinase 2 (SIK2) Multi-State Binding with Bosutinib

The kinase domain is highly conserved among protein kinases 'in terms of both sequence and structure. Conformational rearrangements of the kinase domain are affected by the phosphorylation of residues and the binding of kinase inhibitors. Interestingly, the conformational rearrangement of the active pocket plays an important role in kinase activity and can be used to design novel kinase inhibitors. We characterized the conformational plasticity of the active pocket when bosutinib was bound to salt-inducible kinase 2 (SIK2) using homology modeling and molecular dynamics simulations. Ten different initial complex models were constructed using the Morph server, ranging from open to closed conformations of SIK2 binding with bosutinib. Our simulation showed that bosutinib binds SIK2 with up or down conformations of the P-loop and with all the conformations of the activation loop. In addition, the αC-helix conformation was induced by the conformation of the activation loop, and the salt bridge formed only with its open conformation. The binding affinity of the models was also determined using the molecular mechanics generalized Born surface area method. Bosutinib was found to form a strong binding model with SIK2 and hydrophobic interactions were the dominant factor. This discovery may help guide the design of novel SIK2 inhibitors.

Fumarato and Phthalato Bridged Dinuclear Metal-Organic Cu(II) and Mn(II) Compounds involving Infinite Fumarate-water Assemblies and Unusual Structure-guiding H-bonded Synthons: Antiproliferative Evaluation and Theoretical Studies

Two new dinuclear coordination compounds viz. [Cu2(µ-fum)(phen)2(H2O)6](fum)•6H2O (1) and [Mn2(µ-phth)2(phen)4]•2H2O (2) (phen = 1,10-phenanthroline, fum = fumarate and phth = phthalate) have been synthesized and characterized by elemental analysis, single...

Anti-metastatic Potential of Natural Triterpenoid Cucurbitacin B Against Cholangiocarcinoma Cells by Targeting Src Protein

Owing to the crucial role of Src in cancer metastasis, interruption of Src and its signaling has been considered a promising strategy for cancer metastasis treatment. Cucurbitacin B, a dietary triterpenoid, has been shown to possess anti-proliferative and apoptosis-inducing activities in cholangiocarcinoma (CCA) cells via suppressing the activation of FAK which is a main downstream Src effector. We hypothesized that cucurbitacin B might act as a Src suppressant which conferring anti-metastasis effect against CCA cells. To investigate this, the role of Src in regulating metastasis behavior of CCA cells and the effect of cucurbitacin B on Src-mediated metastatic phenotype of these cells were determined. The results showed that activation of Src significantly enhanced the migratory and invasive abilities of CCA cells. Molecular analysis revealed that Src-facilitated metastasis behavior of CCA cells occurred by modifying expression of a wide range of metastasis-related genes in the cells. Consistent with gene expression results, activation of Src significantly induced the protein expression of 2 important metastasis-associated molecules, MMP-9 and VEGF. Cucurbitacin B markedly suppressed activation of Src and its key effector, FAK. As a consequence, the alteration of expression profiles of metastasis-associated genes induced by Src activator in CCA cells was diminished by cucurbitacin B treatment. The compound also down-regulated Src-induced expression of MMP-9 and VEGF proteins in the cells. Moreover, molecular docking analysis revealed that cucurbitacin B could interact with Src kinase domain and possibly restrain the kinase from being activated by hindering the ATP binding. In conclusion, cucurbitacin B exhibited anti-metastatic property in CCA cells via negatively influencing Src and Src-related oncogenic signaling. This compound may therefore be a potential therapeutic drug for further development as an anti-Src agent for treatment of metastatic CCA.

Small Molecule Kinase Inhibitor Drugs (1995-2021): Medical Indication, Pharmacology, and Synthesis

The central role of dysregulated kinase activity in the etiology of progressive disorders, including cancer, has fostered incremental efforts on drug discovery programs over the past 40 years. As a result, kinase inhibitors are today one of the most important classes of drugs. The FDA approved 73 small molecule kinase inhibitor drugs until September 2021, and additional inhibitors were approved by other regulatory agencies during that time. To complement the published literature on clinical kinase inhibitors, we have prepared a review that recaps this large data set into an accessible format for the medicinal chemistry community. Along with the therapeutic and pharmacological properties of each kinase inhibitor approved across the world until 2020, we provide the synthesis routes originally used during the discovery phase, many of which were only available in patent applications. In the last section, we also provide an update on kinase inhibitor drugs approved in 2021.

Small molecule ERK5 kinase inhibitors paradoxically activate ERK5 signalling: be careful what you wish for…

ERK5 is a protein kinase that also contains a nuclear localisation signal and a transcriptional transactivation domain. Inhibition of ERK5 has therapeutic potential in cancer and inflammation and this has prompted the development of ERK5 kinase inhibitors (ERK5i). However, few ERK5i programmes have taken account of the ERK5 transactivation domain. We have recently shown that the binding of small molecule ERK5i to the ERK5 kinase domain stimulates nuclear localisation and paradoxical activation of its transactivation domain. Other kinase inhibitors paradoxically activate their intended kinase target, in some cases leading to severe physiological consequences highlighting the importance of mitigating these effects. Here, we review the assays used to monitor ERK5 activities (kinase and transcriptional) in cells, the challenges faced in development of small molecule inhibitors to the ERK5 pathway, and classify the molecular mechanisms of paradoxical activation of protein kinases by kinase inhibitors.

Addressing a Trapped High-Energy Water: Design and Synthesis of Highly Potent Pyrimidoindole-Based Glycogen Synthase Kinase-3β Inhibitors

In small molecule binding, water is not a passive bystander but rather takes an active role in the binding site, which may be decisive for the potency of the inhibitor. Here, by addressing a high-energy water, we improved the IC₅₀ value of our co-crystallized glycogen synthase kinase-3β (GSK-3β) inhibitor by nearly two orders of magnitude. Surprisingly, our results demonstrate that this high-energy water was not displaced by our potent inhibitor (S)-3-(3-((7-ethynyl-9H-pyrimido[4,5-b]indol-4-yl)(methyl)amino)piperidin-1-yl)propanenitrile ((S)-15, IC₅₀ value of 6 nM). Instead, only a subtle shift in the location of this water molecule resulted in a dramatic decrease in the energy of this high-energy hydration site, as shown by the WaterMap analysis combined with microsecond timescale molecular dynamics simulations. (S)-15 demonstrated both a favorable kinome selectivity profile and target engagement in a cellular environment and reduced GSK-3 autophosphorylation in neuronal SH-SY5Y cells. Overall, our findings highlight that even a slight adjustment in the location of a high-energy water can be decisive for ligand binding.

Water-Mediated Oxidation of Guanine by a Repair Enzyme: Simulation Using the ABEEM Polarizable Force Field

The recognition mechanism of oxidative damage in organisms has long been a research hotspot. Water is an important medium in the recognition process, but its specific role remains unknown. There is a need to develop a suitable force field that can adequately describe the electrostatic, hydrogen bond, and other interactions among the molecules in the complex system of the repair enzyme and oxidized base. The developing ABEEM polarizable force field (PFF) has been used to simulate the repaired enzyme hOGG1 and oxidized DNA (PDB ID: 1EBM) in a biological environment, and the corresponding results are better than those of the fixed-charge force fields OPLS/AA and AMBER OL15. 8-Oxo-G is recognized by Gln315 of hOGG1 mainly through hydrogen bonds mediated by continuous exchange of 2 water molecules. Phe319 and Cys253 are stacked on both sides of the π planes of bases to form sandwich structures. The charge polarization effect gives an important signal to drive the exchange of water molecules and maintains the recognition of oxidation bases by enzymes. The mediated main water molecule A and mediated auxiliary water molecule B together pull Gln315 to recognize 8-oxo-G by hydrogen bond interactions. Then, the charge polarization signal of solvent water molecule C with a large absolute charge causes the absolute charge of O atoms in water molecule A or B to increase by approximately 0.2 e, and water molecule A or B leaves Gln315 and 8-oxo-G. The other water molecule and water molecule C synergistically recognize 8-oxo-G with Gln315. Even though the water molecules between Gln315 and 8-oxo-G are removed, the MD simulation results show that water molecules appear between Gln315 and 8-oxo-G in a very short time (<2 ps). The dwell time of each water molecule is approximately 60 ps. The radial distribution function and dwell time support the correctness of the above mechanism. These polarization effects and hydrogen bonding interactions cannot be simulated by a fixed-charge force field.

Allostery governs Cdk2 activation and differential recognition of CDK inhibitors

Cyclin-dependent kinases (CDKs) are the master regulators of the eukaryotic cell cycle. To become activated, CDKs require both regulatory phosphorylation and binding of a cognate cyclin subunit. We studied the activation process of the G1/S kinase Cdk2 in solution and developed a thermodynamic model that describes the allosteric coupling between regulatory phosphorylation, cyclin binding and inhibitor binding. The results explain why monomeric Cdk2 lacks activity despite sampling an active-like state, reveal that regulatory phosphorylation enhances allosteric coupling with the cyclin subunit and show that this coupling underlies differential recognition of Cdk2 and Cdk4 inhibitors. We identify an allosteric hub that has diverged between Cdk2 and Cdk4 and show that this hub controls the strength of allosteric coupling. The altered allosteric wiring of Cdk4 leads to compromised activity toward generic peptide substrates and comparative specialization toward its primary substrate retinoblastoma (RB).

Paradoxical activation of c-Src as a drug-resistant mechanism

ATP-competitive inhibitors have been developed as promising anti-cancer agents. However, drug-resistance frequently occurs, and the underlying mechanisms are not fully understood. Here, we show that the activation of c-Src and its downstream phosphorylation cascade can be paradoxically induced by Src-targeted and RTK-targeted kinase inhibitors. We reveal that inhibitor binding induces a conformational change in c-Src, leading to the association of the active form c-Src with focal adhesion kinase (FAK). Reduction of the inhibitor concentration results in the dissociation of inhibitors from the c-Src-FAK complex, which allows c-Src to phosphorylate FAK and initiate FAK-Grb2-mediated Erk signaling. Furthermore, a drug-resistant mutation in c-Src, which reduces the affinity of inhibitors for c-Src, converts Src inhibitors into facilitators of cell proliferation by enhancing the phosphorylation of FAK and Erk in c-Src-mutated cells. Our data thus reveal paradoxical enhancement of cell growth evoked by target-based kinase inhibitors, providing potentially important clues for the future development of effective and safe cancer treatment.

How Far Are We from the Rapid Prediction of Drug Resistance Arising Due to Kinase Mutations?

In all living organisms, protein kinases regulate various cell signaling events through phosphorylation. The phosphorylation occurs upon transferring an ATP's terminal phosphate to a target residue. Because of the central role of protein kinases in several proliferative pathways, point mutations occurring within the kinase's ATP-binding site can lead to a constitutively active enzyme, and ultimately, to cancer. A select set of these point mutations can also make the enzyme drug resistant toward the available kinase inhibitors. Because of technical and economical limitations, rapid experimental exploration of the impact of these mutations remains to be a challenge. This underscores the importance of kinase-ligand binding affinity prediction tools that are poised to measure the efficacy of inhibitors in the presence of kinase mutations. To this end, here, we compare the performances of six web-based scoring tools (DSX-ONLINE, KDEEP, HADDOCK2.2, PDBePISA, Pose&Rank, and PRODIGY-LIG) in assessing the impact of kinase mutations on their interactions with their inhibitors. This assessment is carried out on a new structure-based BINDKIN benchmark we compiled. BINDKIN contains wild-type and mutant structure pairs of kinase-inhibitor complexes, together with their corresponding experimental binding affinities (in the form of IC₅₀, K _d, and K _i). The performance of various web servers over BINDKIN shows that they cannot predict the binding affinities (ΔGs) of wild-type and mutant cases directly. Still, they could catch whether a mutation improves or worsens the ligand binding (ΔΔGs) where the highest Pearson's R correlation coefficient is reached by DSX-ONLINE over the K _i dataset. When homology models are used instead of K _i-associated crystal structures, DSX-ONLINE loses its predictive capacity. These results highlight that there is room to improve the available scoring functions to estimate the impact of protein kinase point mutations on inhibitor binding. The BINDKIN benchmark with all related results is freely accessible online (https://github.com/CSB-KaracaLab/BINDKIN).

Tuning Local Hydration Enables a Deeper Understanding of Protein-Ligand Binding: The PP1-Src Kinase Case

Water plays a key role in biomolecular recognition and binding. Despite the development of several computational and experimental approaches, it is still challenging to comprehensively characterize water-mediated effects on the binding process. Here, we investigate how water affects the binding of Src kinase to one of its inhibitors, PP1. Src kinase is a target for treating several diseases, including cancer. We use biased molecular dynamics simulations, where the hydration of predetermined regions is tuned at will. This computational technique efficiently accelerates the SRC-PP1 binding simulation and allows us to identify several key and yet unexplored aspects of the solvent's role. This study provides a further perspective on the binding phenomenon, which may advance the current drug design approaches for the development of new kinase inhibitors.

Impact of structural biologists and the Protein Data Bank on small-molecule drug discovery and development

The Protein Data Bank (PDB) is an international core data resource central to fundamental biology, biomedicine, bioenergy, and biotechnology/bioengineering. Now celebrating its 50th anniversary, the PDB houses >175,000 experimentally determined atomic structures of proteins, nucleic acids, and their complexes with one another and small molecules and drugs. The importance of three-dimensional (3D) biostructure information for research and education obtains from the intimate link between molecular form and function evident throughout biology. Among the most prolific consumers of PDB data are biomedical researchers, who rely on the open access resource as the authoritative source of well-validated, expertly curated biostructures. This review recounts how the PDB grew from just seven protein structures to contain more than 49,000 structures of human proteins that have proven critical for understanding their roles in human health and disease. It then describes how these structures are used in academe and industry to validate drug targets, assess target druggability, characterize how tool compounds and other small-molecules bind to drug targets, guide medicinal chemistry optimization of binding affinity and selectivity, and overcome challenges during preclinical drug development. Three case studies drawn from oncology exemplify how structural biologists and open access to PDB structures impacted recent regulatory approvals of antineoplastic drugs.

Biologically relevant and energetically significant cooperative ternary (π–π)2/(π–π)1/(π–π)2 assemblies and fascinating discrete (H2O)21 clusters in isostructural 2,5-pyridine dicarboxylato Co(ii) and Zn(ii) phenanthroline compounds: antiproliferative evaluation and theoretical studies

Cytotoxicity in cancer cells with structure activity relationship has been explored in isostructural Co(ii) and Zn(ii) compounds involving energetically significant cooperative (π–π)₂/(π–π)₁/(π–π)₂ assemblies and fascinating (H₂O)₂₁ clusters.

Three-Dimensional Interactions Analysis of the Anticancer Target c-Src Kinase with Its Inhibitors

Src family kinases (SFKs) constitute the biggest family of non-receptor tyrosine kinases considered as therapeutic targets for cancer therapy. An aberrant expression and/or activation of the proto-oncogene c-Src kinase, which is the oldest and most studied member of the family, has long been demonstrated to play a major role in the development, growth, progression and metastasis of numerous human cancers, including colon, breast, gastric, pancreatic, lung and brain carcinomas. For these reasons, the pharmacological inhibition of c-Src activity represents an effective anticancer strategy and a few compounds targeting c-Src, together with other kinases, have been approved as drugs for cancer therapy, while others are currently undergoing preclinical studies. Nevertheless, the development of potent and selective inhibitors of c-Src aimed at properly exploiting this biological target for the treatment of cancer still represents a growing field of study. In this review, the co-crystal structures of c-Src kinase in complex with inhibitors discovered in the past two decades have been described, highlighting the key ligand-protein interactions necessary to obtain high potency and the features to be exploited for addressing selectivity and drug resistance issues, thus providing useful information for the design of new and potent c-Src kinase inhibitors.

Allosteric Regulation at the Crossroads of New Technologies: Multiscale Modeling, Networks, and Machine Learning

Allosteric regulation is a common mechanism employed by complex biomolecular systems for regulation of activity and adaptability in the cellular environment, serving as an effective molecular tool for cellular communication. As an intrinsic but elusive property, allostery is a ubiquitous phenomenon where binding or disturbing of a distal site in a protein can functionally control its activity and is considered as the "second secret of life." The fundamental biological importance and complexity of these processes require a multi-faceted platform of synergistically integrated approaches for prediction and characterization of allosteric functional states, atomistic reconstruction of allosteric regulatory mechanisms and discovery of allosteric modulators. The unifying theme and overarching goal of allosteric regulation studies in recent years have been integration between emerging experiment and computational approaches and technologies to advance quantitative characterization of allosteric mechanisms in proteins. Despite significant advances, the quantitative characterization and reliable prediction of functional allosteric states, interactions, and mechanisms continue to present highly challenging problems in the field. In this review, we discuss simulation-based multiscale approaches, experiment-informed Markovian models, and network modeling of allostery and information-theoretical approaches that can describe the thermodynamics and hierarchy allosteric states and the molecular basis of allosteric mechanisms. The wealth of structural and functional information along with diversity and complexity of allosteric mechanisms in therapeutically important protein families have provided a well-suited platform for development of data-driven research strategies. Data-centric integration of chemistry, biology and computer science using artificial intelligence technologies has gained a significant momentum and at the forefront of many cross-disciplinary efforts. We discuss new developments in the machine learning field and the emergence of deep learning and deep reinforcement learning applications in modeling of molecular mechanisms and allosteric proteins. The experiment-guided integrated approaches empowered by recent advances in multiscale modeling, network science, and machine learning can lead to more reliable prediction of allosteric regulatory mechanisms and discovery of allosteric modulators for therapeutically important protein targets.

Impact of the Protein Data Bank on antineoplastic approvals

Open access to 3D structure information from the Protein Data Bank (PDB) facilitated discovery and development of >90% of the 79 new antineoplastic agents (54 small molecules, 25 biologics) with known molecular targets approved by the FDA 2010-2018. Analyses of PDB holdings, the scientific literature and related documents for each drug-target combination revealed that the impact of public-domain 3D structure data was broad and substantial, ranging from understanding target biology (∼95% of all targets) to identifying a given target as probably druggable (∼95% of all targets) to structure-guided lead optimization (>70% of all small-molecule drugs). In addition to aggregate impact assessments, illustrative case studies are presented for three protein kinase inhibitors, an allosteric enzyme inhibitor and seven advanced-stage melanoma therapeutics.

Paradoxical activation of the protein kinase-transcription factor ERK5 by ERK5 kinase inhibitors

The dual protein kinase-transcription factor, ERK5, is an emerging drug target in cancer and inflammation, and small-molecule ERK5 kinase inhibitors have been developed. However, selective ERK5 kinase inhibitors fail to recapitulate ERK5 genetic ablation phenotypes, suggesting kinase-independent functions for ERK5. Here we show that ERK5 kinase inhibitors cause paradoxical activation of ERK5 transcriptional activity mediated through its unique C-terminal transcriptional activation domain (TAD). Using the ERK5 kinase inhibitor, Compound 26 (ERK5-IN-1), as a paradigm, we have developed kinase-active, drug-resistant mutants of ERK5. With these mutants, we show that induction of ERK5 transcriptional activity requires direct binding of the inhibitor to the kinase domain. This in turn promotes conformational changes in the kinase domain that result in nuclear translocation of ERK5 and stimulation of gene transcription. This shows that both the ERK5 kinase and TAD must be considered when assessing the role of ERK5 and the effectiveness of anti-ERK5 therapeutics.

Nicotinamide-Ponatinib Analogues as Potent Anti-CML and Anti-AML Compounds

Ponatinib is a multikinase inhibitor that is used to treat chronic myeloid leukemia patients harboring mutated ABL1(T315I) kinase. Due to the potent inhibition of FLT3, RET, and fibroblast growth factor receptors (FGFRs), it is also being evaluated against acute myeloid leukemia (AML), biliary, and lung cancers. The multikinase inhibition profile of ponatinib may also account for its toxicity, thus analogs with improved kinase selectivity or different kinase inhibition profiles could be better tolerated. The introduction of nitrogen into drug compounds can enhance efficacy and drug properties (a concept called "necessary nitrogen"). Here, we introduce additional nitrogen into the benzamide moiety of ponatinib to arrive at nicotinamide analogs. A nicotinamide analogue of ponatinib, HSN748, retains activity against FLT3, ABL1, RET, and PDGFRα/β but loses activity against c-Src and P38α. MNK1 and 2 are key kinases that phosphorylate eIF4E to regulate the protein translation complex. MNK also modulates mTORC1 signaling and contributes to rapamycin resistance. Inhibitors of MNK1 and 2 are being evaluated for anticancer therapy. Ponatinib is not a potent inhibitor of MNK1 or 2, but the nicotinamide analogs are potent inhibitors of MNKs. This illustrates a powerful demonstration of the necessary nitrogen concept to alter both the potency and selectivity of drugs.

Adipato bridged novel hexanuclear Cu(ii) and polymeric Co(ii) coordination compounds involving cooperative supramolecular assemblies and encapsulated guest water clusters in a square grid host: antiproliferative evaluation and theoretical studies

Novel adipato bridged Cu(ii) and Co(ii) complexes synthesized by considering cytotoxicity, apoptosis, ROS generation, molecular docking and pharmacophore features.

Design and Analysis of the 4-Anilinoquin(az)oline Kinase Inhibition Profiles of GAK/SLK/STK10 Using Quantitative Structure-Activity Relationships

The 4-anilinoquinoline and 4-anilinoquinazoline ring systems have been the focus of significant efforts in prior kinase drug discovery programs, which have led to approved medicines. Broad kinome profiles of these compounds have now been assessed with the advent of advanced screening technologies. These ring systems, while originally designed for specific targets including epidermal growth factor receptor (EGFR), but actually display a number of potent collateral kinase targets, some of which have been associated with negative clinical outcomes. We have designed and synthesized a series of 4-anilinoquin(az)olines in order to better understand the structure-activity relationships of three main collateral kinase targets of quin(az)oline-based kinase inhibitors: cyclin G associated kinase (GAK), STE20-like serine/threonine-protein kinase (SLK) and serine/threonine-protein kinase 10 (STK10). This was achieved through a series of quantitative structure-activity relationship (QSAR) analysis, water mapping of the kinase ATP binding sites and extensive small-molecule X-ray structural analysis.

The role of hydration effects in 5-fluorouridine binding to SOD1: insight from a new 3D-RISM-KH based protocol for including structural water in docking simulations

Misfolded Cu/Zn superoxide dismutase enzyme (SOD1) shows prion-like propagation in neuronal cells leading to neurotoxic aggregates that are implicated in amyotrophic lateral sclerosis (ALS). Tryptophan-32 (W32) in SOD1 is part of a potential site for templated conversion of wild type SOD1. This W32 binding site is located on a convex, solvent exposed surface of the SOD1 suggesting that hydration effects can play an important role in ligand recognition and binding. A recent X-ray crystal structure has revealed that 5-Fluorouridine (5-FUrd) binds at the W32 binding site and can act as a pharmacophore scaffold for the development of anti-ALS drugs. In this study, a new protocol is developed to account for structural (non-displaceable) water molecules in docking simulations and successfully applied to predict the correct docked conformation binding modes of 5-FUrd at the W32 binding site. The docked configuration is within 0.58 Å (RMSD) of the observed configuration. The docking protocol involved calculating a hydration structure around SOD1 using molecular theory of solvation (3D-RISM-KH, 3D-Reference Interaction Site Model-Kovalenko-Hirata) whereby, non-displaceable water molecules are identified for docking simulations. This protocol was also used to analyze the hydrated structure of the W32 binding site and to explain the role of solvation in ligand recognition and binding to SOD1. Structural water molecules mediate hydrogen bonds between 5-FUrd and the receptor, and create an environment favoring optimal placement of 5-FUrd in the W32 binding site.

How To Design Selective Ligands for Highly Conserved Binding Sites: A Case Study Using N-Myristoyltransferases as a Model System

A model system of two related enzymes with conserved binding sites, namely N-myristoyltransferase from two different organisms, was studied to decipher the driving forces that lead to selective inhibition in such cases. Using a combination of computational and experimental tools, two different selectivity-determining features were identified. For some ligands, a change in side-chain flexibility appears to be responsible for selective inhibition. Remarkably, this was observed for residues orienting their side chains away from the ligands. For other ligands, selectivity is caused by interfering with a water molecule that binds more strongly to the off-target than to the target. On the basis of this finding, a virtual screen for selective compounds was conducted, resulting in three hit compounds with the desired selectivity profile. This study delivers a guideline on how to assess selectivity-determining features in proteins with conserved binding sites and to translate this knowledge into the design of selective inhibitors.