Structure-kinetic Relationships Research Articles

Pitfalls and Considerations in Determining the Potency and Mutant Selectivity of Covalent Epidermal Growth Factor Receptor Inhibitors.

Enzyme inhibitors that form covalent bonds with their targets are being increasingly pursued in drug development. Assessing their biochemical activity relies on time-dependent assays, which are distinct and more complex compared with methods commonly employed for reversible-binding inhibitors. To provide general guidance to the covalent inhibitor development community, we explored methods and reported kinetic values and experimental factors in determining the biochemical activity of various covalent epidermal growth factor receptor (EGFR) inhibitors. We showcase how liquid handling and assay reagents impact kinetic parameters and potency interpretations, which are critical for structure-kinetic relationships and covalent drug design. Additionally, we include benchmark kinetic values with reference inhibitors, which are imperative, as covalent EGFR inhibitor kinetic values are infrequently consistent in the literature. This overview seeks to inform best practices for developing new covalent inhibitors and highlight appropriate steps to address gaps in knowledge presently limiting assay reliability and reproducibility.

Insight into Thermodynamic and Kinetic Profiles in Small-Molecule Optimization.

Structure-activity relationships (SARs) and structure-property relationships (SPRs) have been considered the most important factors during the drug optimization process. For medicinal chemists, improvements in the potencies and druglike properties of small molecules are regarded as their major goals. Among them, the binding affinity and selectivity of small molecules on their targets are the most important indicators. In recent years, there has been growing interest in using thermodynamic and kinetic profiles to analyze ligand-receptor interactions, which could provide not only binding affinities but also detailed binding parameters for small-molecule optimization. In this perspective, we are trying to provide an insight into thermodynamic and kinetic profiles in small-molecule optimization. Through a highlight of strategies on the small-molecule optimization with specific cases, we aim to put forward the importance of structure-thermodynamic relationships (STRs) and structure-kinetic relationships (SKRs), which could provide more guidance to find safe and effective small-molecule drugs.

Kinetic profiling and functional characterization of 8-phenylxanthine derivatives as A2B adenosine receptor antagonists

A2B adenosine receptor (A2BAR) antagonists have therapeutic potential in inflammation-related diseases such as asthma, chronic obstructive pulmonary disease and cancer. However, no drug is currently clinically approved, creating a demand for research on novel antagonists. Over the last decade, the study of target binding kinetics, along with affinity and potency, has been proven valuable in early drug discovery stages, as it is associated with improved in vivo drug efficacy and safety. In this study, we report the synthesis and biological evaluation of a series of xanthine derivatives as A2BAR antagonists, including an isothiocyanate derivative designed to bind covalently to the receptor. All 28 final compounds were assessed in radioligand binding experiments, to evaluate their affinity and for those qualifying, kinetic binding parameters. Both structure-affinity and structure-kinetic relationships were derived, providing a clear relationship between affinity and dissociation rate constants. Two structurally similar compounds, 17 and 18, were further evaluated in a label-free assay due to their divergent kinetic profiles. An extended cellular response was associated with long A2BAR residence times. This link between a ligand’s A2BAR residence time and its functional effect highlights the importance of binding kinetics as a selection parameter in the early stages of drug discovery.

Identification of Slow-Binding Inhibitors of the BoNT/A Protease.

Botulinum neurotoxin A (BoNT/A) is a lethal toxin, which causes botulism, and is categorized as a bioterrorism threat, which causes flaccid paralysis and death. Botulinum A neurotoxicity is governed through its light chain (LC), a zinc metalloprotease. Pharmacological investigations aimed at negating BoNT/A's LC have typically looked to inhibitors that have been shown to inhibit the light chain's activity by reversible zinc chelation within its active site. This report outlines the first examples of nonpeptidic inhibitors of the BoNT/A LC that possess slow-binding kinetics, a needed logic to counteract the longevity of BoNT/A. Cyclopropane, alkyl, and alkenyl derivatives of 2,4-dichlorocinamic hydroxamic acid (DCHA) were shown to possess both one-step and two-step slow-binding kinetics. Structure-kinetic relationships (SKRs) were observed and were rationalized with the aid of docking models that predicted improved interactions with residues within a hydrophobic cleft adjacent to the active site.

Unraveling binding mechanism and kinetics of macrocyclic Gαq protein inhibitors

G proteins represent intracellular switches that transduce signals relayed from G protein-coupled receptors. The structurally related macrocyclic depsipeptides FR900359 (FR) and YM-254890 (YM) are potent, selective inhibitors of the Gαq protein family. We recently discovered that radiolabeled FR and YM display strongly divergent residence times, which translates into significantly longer antiasthmatic effects of FR. The present study is aimed at investigating the molecular basis for this observed disparity. Based on docking studies, we mutated amino acid residues of the Gαq protein predicted to interact with FR or YM, and recombinantly expressed the mutated Gαq proteins in cells in which the native Gαq proteins had been knocked out by CRISPR-Cas9. Both radioligands showed similar association kinetics, and their binding followed a conformational selection mechanism, which was rationalized by molecular dynamics simulation studies. Several mutations of amino acid residues near the putative binding site of the “lipophilic anchors” of FR, especially those predicted to interact with the isopropyl group present in FR but not in YM, led to dramatically accelerated dissociation kinetics. Our data indicate that the long residence time of FR depends on lipophilic interactions within its binding site. The observed structure-kinetic relationships point to a complex binding mechanism of FR, which likely involves snap-lock- or dowel-like conformational changes of either ligand or protein, or both. These experimental data will be useful for the design of compounds with a desired residence time, a parameter that has now been recognized to be of utmost importance in drug development.

Prediction of the Drug–Target Binding Kinetics for Flexible Proteins by Comparative Binding Energy Analysis

There is growing consensus that the optimization of the kinetic parameters for drug-protein binding leads to improved drug efficacy. Therefore, computational methods have been developed to predict kinetic rates and to derive quantitative structure-kinetic relationships (QSKRs). Many of these methods are based on crystal structures of ligand-protein complexes. However, a drawback is that each ligand-protein complex is usually treated as having a single structure. Here, we present a modification of COMparative BINding Energy (COMBINE) analysis, which uses the structures of ligand-protein complexes to predict binding parameters. We introduce the option of using multiple structures to describe each ligand-protein complex in COMBINE analysis and apply this to study the effects of protein flexibility on the derivation of dissociation rate constants (koff) for inhibitors of p38 mitogen-activated protein (MAP) kinase, which has a flexible binding site. Multiple structures were obtained for each ligand-protein complex by performing docking to an ensemble of protein configurations obtained from molecular dynamics simulations. Coefficients to scale ligand-protein interaction energies determined from energy-minimized structures of ligand-protein complexes were obtained by partial least squares regression, and they allowed for the computation of koff values. The QSKR model obtained using single, energy-minimized crystal structures for each ligand-protein complex had higher predictive power than the QSKR model obtained with multiple structures from ensemble docking. However, incorporation of ligand-protein flexibility helped to highlight additional ligand-protein interactions that lead to longer residence times, such as interactions with residues Arg67 and Asp168, which are close to the ligand in many crystal structures. These results show that COMBINE analysis is a promising method to guide the design of compounds that bind to flexible proteins with improved binding kinetics.

Machine Learning Analysis of τRAMD Trajectories to Decipher Molecular Determinants of Drug-Target Residence Times.

Drug-target residence times can impact drug efficacy and safety, and are therefore increasingly being considered during lead optimization. For this purpose, computational methods to predict residence times, τ, for drug-like compounds and to derive structure-kinetic relationships are desirable. A challenge for approaches based on molecular dynamics (MD) simulation is the fact that drug residence times are typically orders of magnitude longer than computationally feasible simulation times. Therefore, enhanced sampling methods are required. We recently reported one such approach: the τRAMD procedure for estimating relative residence times by performing a large number of random acceleration MD (RAMD) simulations in which ligand dissociation occurs in times of about a nanosecond due to the application of an additional randomly oriented force to the ligand. The length of the RAMD simulations is used to deduce τ. The RAMD simulations also provide information on ligand egress pathways and dissociation mechanisms. Here, we describe a machine learning approach to systematically analyze protein-ligand binding contacts in the RAMD trajectories in order to derive regression models for estimating τ and to decipher the molecular features leading to longer τ values. We demonstrate that the regression models built on the protein-ligand interaction fingerprints of the dissociation trajectories result in robust estimates of τ for a set of 94 drug-like inhibitors of heat shock protein 90 (HSP90), even for the compounds for which the length of the RAMD trajectories does not provide a good estimation of τ. Thus, we find that machine learning helps to overcome inaccuracies in the modeling of protein-ligand complexes due to incomplete sampling or force field deficiencies. Moreover, the approach facilitates the identification of features important for residence time. In particular, we observed that interactions of the ligand with the sidechain of F138, which is located on the border between the ATP binding pocket and a hydrophobic transient sub-pocket, play a key role in slowing compound dissociation. We expect that the combination of the τRAMD simulation procedure with machine learning analysis will be generally applicable as an aid to target-based lead optimization.

Toward the elucidation of the mechanism for passive membrane permeability of cyclic peptides.

Toward the elucidation of the mechanism for passive membrane permeability of cyclic peptides.

Tert-Butylphenylthiazoles with an oxadiazole linker: a novel orally bioavailable class of antibiotics exhibiting antibiofilm activity.

The structure–activity and structure–kinetic relationships of a new tert-butylphenylthiazole series with oxadiazole linkers were conducted with the objective of obtaining a new orally available antibacterial compounds. Twenty-two new compounds were prepared, purified and identified. Their activity against methicillin-resistant Staphylococcus aureus were examined. Compound 20 with 3-hydroxyazetidine as a nitrogenous side chain showed promising activity against twenty-four clinical isolates, including vancomycin-resistant staphylococcal and enterococcal species with MIC values ranging from 4–8 μg mL−1. Additional advantages of this compound include an ability to eradicate staphylococcal biofilm mass in a dose-dependent manner as well as high metabolic stability after an oral dose of 25 mg kg−1 with a biological half-life that exceeds 5 hours and a plasma concentration (Cmax) that exceeds the MIC values.

KinITC-One Method Supports both Thermodynamic and Kinetic SARs as Exemplified on FimH Antagonists.

Affinity data, such as dissociation constants (KD ) or inhibitory concentrations (IC50 ), are widely used in drug discovery. However, these parameters describe an equilibrium state, which is often not established in vivo due to pharmacokinetic effects and they are therefore not necessarily sufficient for evaluating drug efficacy. More accurate indicators for pharmacological activity are the kinetics of binding processes, as they shed light on the rate of formation of protein-ligand complexes and their half-life. Nonetheless, although highly desirable for medicinal chemistry programs, studies on structure-kinetic relationships (SKR) are still rare. With the recently introduced analytical tool kinITC this situation may change, since not only thermodynamic but also kinetic information of the binding process can be deduced from isothermal titration calorimetry (ITC) experiments. Using kinITC, ITC data of 29 mannosides binding to the bacterial adhesin FimH were re-analyzed to make their binding kinetics accessible. To validate these kinetic data, surface plasmon resonance (SPR) experiments were conducted. The kinetic analysis by kinITC revealed that the nanomolar affinities of the FimH antagonists arise from both (i) an optimized interaction between protein and ligand in the bound state (reduced off-rate constant koff ) and (ii) a stabilization of the transition state or a destabilization of the unbound state (increased on-rate constant kon ). Based on congeneric ligand modifications and structural input from co-crystal structures, a strong relationship between the formed hydrogen-bond network and koff could be concluded, whereas electrostatic interactions and conformational restrictions upon binding were found to have mainly an impact on kon .

Ligand Desolvation Steers On-Rate and Impacts Drug Residence Time of Heat Shock Protein 90 (Hsp90) Inhibitors

Residence time and more recently the association rate constant kon are increasingly acknowledged as important parameters for in vivo efficacy and safety of drugs. However, their broader consideration in drug development is limited by a lack of knowledge of how to optimize these parameters. In this study on a set of 176 heat shock protein 90 inhibitors, structure-kinetic relationships, X-ray crystallography, and molecular dynamics simulations were combined to retrieve a concrete scheme of how to rationally slow down on-rates. We discovered that an increased ligand desolvation barrier by introducing polar substituents resulted in a significant kon decrease. The slowdown was accomplished by introducing polar moieties to those parts of the ligand that point toward a hydrophobic cavity. We validated this scheme by increasing polarity of three Hsp90 inhibitors and observed a 9-, 13-, and 45-fold slowdown of on-rates and a 9-fold prolongation in residence time. This prolongation was driven by transition state destabilization rather than ground state stabilization.

Discovery of 7-Oxo-2,4,5,7-tetrahydro-6 H-pyrazolo[3,4- c]pyridine Derivatives as Potent, Orally Available, and Brain-Penetrating Receptor Interacting Protein 1 (RIP1) Kinase Inhibitors: Analysis of Structure-Kinetic Relationships.

We report the discovery of 7-oxo-2,4,5,7-tetrahydro-6 H-pyrazolo[3,4- c]pyridine derivatives as a novel class of receptor interacting protein 1 (RIP1) kinase inhibitors. On the basis of the overlay study between HTS hit 10 and GSK2982772 (6) in RIP1 kinase, we designed and synthesized a novel class of RIP1 kinase inhibitor 11 possessing moderate RIP1 kinase inhibitory activity and P-gp mediated efflux. The optimization of the core structure and the exploration of appropriate substituents utilizing SBDD approach led to the discovery of 22, a highly potent, orally available, and brain-penetrating RIP1 kinase inhibitor with excellent PK profiles. Compound 22 significantly suppressed necroptotic cell death both in mouse and human cells. Oral administration of 22 (10 mg/kg, bid) attenuated disease progression in the mouse experimental autoimmune encephalomyelitis (EAE) model of multiple sclerosis (MS). Moreover, analysis of structure-kinetic relationship (SKR) for our novel chemical series was also discussed.

Structure-Affinity Relationships and Structure-Kinetic Relationships of 1,2-Diarylimidazol-4-carboxamide Derivatives as Human Cannabinoid 1 Receptor Antagonists.

We report on the synthesis and biological evaluation of a series of 1,2-diarylimidazol-4-carboxamide derivatives developed as CB1 receptor antagonists. These were evaluated in a radioligand displacement binding assay, a [35S]GTPγS binding assay, and in a competition association assay that enables the relatively fast kinetic screening of multiple compounds. The compounds show high affinities and a diverse range of kinetic profiles at the CB1 receptor and their structure–kinetic relationships (SKRs) were established. Using the recently resolved hCB1 receptor crystal structures, we also performed a modeling study that sheds light on the crucial interactions for both the affinity and dissociation kinetics of this family of ligands. We provide evidence that, next to affinity, additional knowledge of binding kinetics is useful for selecting new hCB1 receptor antagonists in the early phases of drug discovery.

Drug-Target Kinetics in Drug Discovery.

The development of therapies for the treatment of neurological cancer faces a number of major challenges including the synthesis of small molecule agents that can penetrate the blood-brain barrier (BBB). Given the likelihood that in many cases drug exposure will be lower in the CNS than in systemic circulation, it follows that strategies should be employed that can sustain target engagement at low drug concentration. Time dependent target occupancy is a function of both the drug and target concentration as well as the thermodynamic and kinetic parameters that describe the binding reaction coordinate, and sustained target occupancy can be achieved through structural modifications that increase target (re)binding and/or that decrease the rate of drug dissociation. The discovery and deployment of compounds with optimized kinetic effects requires information on the structure–kinetic relationships that modulate the kinetics of binding, and the molecular factors that control the translation of drug–target kinetics to time-dependent drug activity in the disease state. This Review first introduces the potential benefits of drug-target kinetics, such as the ability to delineate both thermodynamic and kinetic selectivity, and then describes factors, such as target vulnerability, that impact the utility of kinetic selectivity. The Review concludes with a description of a mechanistic PK/PD model that integrates drug–target kinetics into predictions of drug activity.

Rationalizing the Binding Kinetics for the Inhibition of the Burkholderia pseudomallei FabI1 Enoyl-ACP Reductase.

There is growing awareness of the link between drug-target residence time and in vivo drug activity, and there are increasing efforts to determine the molecular factors that control the lifetime of a drug-target complex. Rational alterations in the drug-target residence time require knowledge of both the ground and transition states on the inhibition reaction coordinate, and we have determined the structure-kinetic relationship for 22 ethyl- or hexyl-substituted diphenyl ethers that are slow-binding inhibitors of bpFabI1, the enoyl-ACP reductase FabI1 from Burkholderia pseudomallei. Analysis of enzyme inhibition using a two-dimensional kinetic map demonstrates that the ethyl and hexyl diphenyl ethers fall into two distinct clusters. Modifications to the ethyl diphenyl ether B ring result in changes to both on and off rates, where residence times of up to ∼700 min (∼11 h) are achieved by either ground state stabilization (PT444) or transition state destabilization (slower on rate) (PT404). By contrast, modifications to the hexyl diphenyl ether B ring result in residence times of 300 min (∼5 h) through changes in only ground state stabilization (PT119). Structural analysis of nine enzyme:inhibitor complexes reveals that the variation in structure-kinetic relationships can be rationalized by structural rearrangements of bpFabI1 and subtle changes to the orientation of the inhibitor in the binding pocket. Finally, we demonstrate that three compounds with residence times on bpFabI1 from 118 min (∼2 h) to 670 min (∼11 h) have in vivo efficacy in an acute B. pseudomallei murine infection model using the virulent B. pseudomallei strain Bp400.

Binding kinetics in drug discovery - A current perspective.

The impact of target binding kinetics (BK) on the clinical performance of therapeutic agents is presently a topic of intense debate in drug discovery. While retrospective studies suggest that BK is a differentiating parameter in marketed medicines, it is yet unclear how this information could be used to prioritize drug candidates during lead optimization. Motivated by the question whether BK can be understood and rationally optimized, we review the most relevant literature in the field, with special focus on selected examples from our organization. First we discuss structure-kinetic relationships (SKR), and how they can be influenced by factors such as conformational changes, molecular flexibility, hydrogen bonds, hydrophobicity, water molecules and (reversible-) covalent bonds. We then introduce the methodologies currently used for the investigation of BK parameters, briefly commenting on their strengths, weaknesses and future trends. Finally, we present our current perspective on the integration of BK in the drug discovery process, aiming to stimulate further thoughts on this important subject.

Structure-Kinetic Relationships of Passive Membrane Permeation from Multiscale Modeling.

Passive membrane permeation of small molecules is essential to achieve the required absorption, distribution, metabolism, and excretion (ADME) profiles of drug candidates, in particular intestinal absorption and transport across the blood-brain barrier. Computational investigations of this process typically involve either building QSAR models or performing free energy calculations of the permeation event. Although insightful, these methods rarely bridge the gap between computation and experiment in a quantitative manner, and identifying structural insights to apply toward the design of compounds with improved permeability can be difficult. In this work, we combine molecular dynamics simulations capturing the kinetic steps of permeation at the atomistic level with a dynamic mechanistic model describing permeation at the in vitro level, finding a high level of agreement with experimental permeation measurements. Calculation of the kinetic rate constants determining each step in the permeation event allows derivation of structure-kinetic relationships of permeation. We use these relationships to probe the structural determinants of membrane permeation, finding that the desolvation/loss of hydrogen bonding required to leave the membrane partitioned position controls the membrane flip-flop rate, whereas membrane partitioning determines the rate of leaving the membrane.

Kinetic and structural insights into the binding of histone deacetylase 1 and 2 (HDAC1, 2) inhibitors

The structure–activity and structure–kinetic relationships of a series of novel and selective ortho-aminoanilide inhibitors of histone deacetylases (HDACs) 1 and 2 are described. Different kinetic and thermodynamic selectivity profiles were obtained by varying the moiety occupying an 11Å channel leading to the Zn2+ catalytic pocket of HDACs 1 and 2, two paralogs with a high degree of structural similarity. The design of these novel inhibitors was informed by two ligand-bound crystal structures of truncated hHDAC2. BRD4884 and BRD7232 possess kinetic selectivity for HDAC1 versus HDAC2. We demonstrate that the binding kinetics of HDAC inhibitors can be tuned for individual isoforms in order to modulate target residence time while retaining functional activity and increased histone H4K12 and H3K9 acetylation in primary mouse neuronal cell culture assays. These chromatin modifiers, with tuned binding kinetic profiles, can be used to define the relation between target engagement requirements and the pharmacodynamic response of HDACs in different disease applications.

Structure-activity relationships (SAR) and structure-kinetic relationships (SKR) of sulphone-based CRTh2 antagonists

Monocyclic and bicyclic ring systems were investigated as the “core” section of a series of diphenylsulphone-containing acetic acid CRTh2 receptor antagonists. A range of potencies were observed and single-digit nanomolar potencies were obtained in both the monocyclic and bicyclic cores. Residence times for the monocyclic compounds were very short. Some of the bicyclic cores displayed better residence times. A methyl group in the northern part of the core, between the head and tail was a necessary requirement for the beginnings of long residence times. Variations of the tail substitution maximised potencies and residence times.

Design strategies to address kinetics of drug binding and residence time.

The kinetics of drug binding and drug residence time are recognized to be important in the clinical effectiveness of drug candidates. In most cases a long residence time of the drug-target complex results in an extended duration of pharmacodynamic activity, even when systemic concentrations of drug have been notably reduced through elimination routes. Hence, if selective for target, long residence times can increase the duration of drug efficacy in vivo and can significantly diminish the potential for off-target-mediated toxicities. Furthermore, a compound with a slower dissociation rate may allow a reduced dosing schedule relative to a compound with a rapid dissociation rate. Factors contributing to long residence time that could be useful to medicinal chemists in the prospective design of compounds with long residence times will be discussed in this perspective. Particular emphasis will be on case studies highlighting how kinetics can be measured, modulated based on supporting structure kinetic relationships and whether these effects are translatable into man.