Quantifying the Interactions between Biomolecules: Guidelines for

Subscriber access provided by UNIV OF TEXAS DALLAS

Perspective

Quantifying the Interactions between Biomolecules: Guidelines for Assay Design and Data Analysis Peter J Tonge ACS Infect. Dis., Just Accepted Manuscript • DOI: 10.1021/acsinfecdis.9b00012 • Publication Date (Web): 12 Mar 2019 Downloaded from http://pubs.acs.org on March 16, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Infectious Diseases

Quantifying the Interactions between Biomolecules: Guidelines for Assay Design and Data Analysis

Peter J. Tonge* Center for Advanced Study of Drug Action, Departments of Chemistry and Radiology, John S. Toll Drive, Stony Brook University, Stony Brook, NY 11794-3400, United States *Email: [email protected]

1 ACS Paragon Plus Environment

ACS Infectious Diseases 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The accurate and precise determination of binding interactions plays a central role in fields such as drug discovery where structure-activity relationships guide the selection and optimization of drug leads. Binding is often assessed by monitoring the response caused by varying one of the binding partners in a functional assay, or by using methods where the concentrations of free and/or bound ligand can be directly determined. In addition, there are also many approaches where binding leads to a change in the properties of the binding partner(s) that can be directly quantified such as an alteration in mass or in a spectroscopic signal. The analysis of data resulting from these techniques invariably relies on computer software that enable rapid fitting of the data to non-linear multiparameter equations. The objective of this perspective is to serve as a reminder of the basic assumptions that are used in deriving these equations and thus that should be considered during assay design and subsequent data analysis. The result is a set of guidelines for authors considering submitting their work to journals such as ACS Infectious Diseases.

Keywords: Binding, Concentration (dose)-response curves, IC50, Nonlinear regression, Replicates, Reproducibility


Page 2 of 30

Page 3 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


Biomolecular interactions, such as protein-protein, protein-ligand and protein-nucleic acid interactions, occur because the complex that is formed is thermodynamically more stable than the unbound species. However, although thermodynamics provide the driving force for binding, the rate of formation and breakdown of the complex is a function of the transition state barrier on the binding reaction coordinate. A full description of the binding event thus requires parameters such as IC50 or Kd values that report on the thermodynamic stability of the complex as well as the on and off rates (kon and koff) that quantify the lifetime of the complex. Many methods now exist for determining both the thermodynamics and kinetics of biomolecular interactions, and are routinely employed across biological space, for example by underpinning the generation of structure-reactivity relationships (SAR) or structure-kinetic relationships (SKR) that guide the selection and optimization of drug leads. This includes competitive radioligand binding,1 the use of mass spectrometry to quantify unbound ligand,2 methods that detect binding such as fluorescence (anisotropy, FRET),3-6 bioluminescence,7 or surface plasmon resonance (SPR),8,9 and assays that monitor change in activity as a function of ligand (agonist, antagonist, inhibitor) concentration.10 Access to user-friendly programs that can fit data to nonlinear equations has greatly facilitated the ability of investigators to obtain quantitative insight into their systems. In particular, there is no longer a reliance on linearized transformations of common equations which can distort experimental errors, such as the Scatchard plot for equilibrium binding data and the Lineweaver Burk plot for enzyme kinetics.11 However, use of programs such as GraphPad Prism or GraFit to perform nonlinear regression generally requires no knowledge of the assumptions and precepts that underlie the equations used for data analysis. There are of course many excellent, authoritative books and papers on assay design and data analysis. This includes Robert Copeland’s book “Evaluation of enzyme inhibitors in drug discovery: a guide for medicinal chemists and pharmacologists”,12 and the Assay Guidance Manual (Eli Lilly and NCATS).13 This perspective is not intended to replace these sources of information. Nor do I attempt to discuss underlying complexities such as the statistics of nonlinear regression.14 Rather I seek to raise awareness of some common misunderstandings and thus provide guidelines for authors wishing to publish in ACS ID. In this perspective, I use the interaction of small molecule chemical compounds with proteins as a paradigm for binding interactions.

Rigor and Reproducibility: Biological and Technical Replicates Many organizations now offer courses and seminars on the Responsible Conduct of Research (RCR). For example, the US National Institutes of Health (NIH) requires that all students funded on NIH grants receive training in RCR.15 One component of RCR is Data Management which includes topics such as the rigor and reproducibility of scientific research. Obviously, reproducibility is a fundamental goal in the design and development of appropriate, robust assays for quantifying biomolecular interactions. In this regard I thought it worth clarifying the difference between biological and technical replicates which are a component of assessing the reproducibility of scientific findings and often a subject of some misunderstanding. Blainey et al.16 define biological replicates as “parallel measurements of biologically distinct samples that capture random biological variation, which may itself be a subject of study or a noise source” and technical replicates as “repeated measurements of the same sample that represent independent measures of the random noise associated with protocols or equipment.” Thus, technical 3 ACS Paragon Plus Environment


replicates provide information on the precision of the measurement method, while biological replicates inform about sample to sample variation in the behavior of separate reagent preps, cell cultures, or animals.17 The number of replicates and their treatment in subsequent data analysis, for example whether or not the replicates are averaged, depends on the type of experiment and its purpose. Below I briefly introduce nonlinear regression and discuss replicates in the context of curve fitting.

Guidelines for Nonlinear Regression Nonlinear regression is a method for fitting an experimental x-y data set, such as a concentrationresponse relationship, to a mathematical equation. This is achieved by systematically varying the values of the parameters in the equation until the parameter values giving the best agreement between the data and the equation are found. The best fit is defined as the set of parameter values that minimize the squared differences between the measured and calculated y values, summed over all data points (so-called “least-squares” regression). There are several considerations in applying nonlinear regression, including the choice of model (i.e. the fitting equation), whether any parameters should be constrained (such as the Hill coefficient h in an IC50 model), choice of initial values for each parameter, how replicate data points are treated, deciding whether and how to weight the data points, and how to detect and handle outliers. Although this Perspective does not attempt to discuss each of the above topics, it is possible to give some initial guidelines especially on the treatment of replicates. Investigators may use whatever level of replication they consider appropriate for measurements that are exclusively aimed at helping make decisions on how best to proceed, and that are not intended for publication. However, minimal standards of reproducibility must be met for any result to be publishable, and high standards of reproducibility are required for results on which a major conclusion depends. For example, in general any IC50 value reported in a publication should be determined using replicate (typically triplicate) measurements at each inhibitor concentration, and the entire IC50 measurement should be repeated at least once (it being acceptable to use the same enzyme and inhibitor preps) to show that it is reproducible. In characterizing the final, optimized inhibitor compound(s), on which the manuscript’s claims of important biological activity or other major conclusions are based, higher levels of replication are typically required. These replicates will ideally include using separate preparations of enzyme and inhibitor. This higher level of rigor is required because the reproducibility of the results obtained using these key compounds is central to the validity of the entire publication, and it is therefore essential to show that the results are reproducible, and that activity does not vary in unexpected ways from one enzyme prep to another, and/or does not result from some contaminant in the inhibitor prep. Regardless of the experimental design, investigators should clearly specify how many and what kind of replicates were performed for each experiment reported. In an enzymatic assay, the duplicate or triplicate initial velocity measurements at each inhibitor concentration may come from different wells in the same multi-well plate, or from individual enzyme assays performed in a cuvette. Assuming that the reagents are stable, and the stock solutions are homogeneous (i.e. all components are fully soluble), the differences in the results obtained between the replicate wells provide a measure of the stochastic variability in diluting and dispensing the reagents, 4 ACS Paragon Plus Environment

Page 4 of 30

Page 5 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


plus any irreproducibility in the performance of the detection instrument (e.g. the plate reader). Two different approaches are commonly employed when applying nonlinear regression to such a data set. In one approach, each measurement is treated as an independent data point, and the equation is fitted simultaneously to all data points. Alternatively, some investigators will instead average the replicate measurements, and then fit the resulting average values to the equation. In most curve fitting programs, it is possible to weight the averaged values according to the spread or the standard deviation of the replicate measurements for each condition. Doing so ensures that the fitting process places greater weight on fitting the averaged values that were determined most precisely, as shown by the close agreement between the replicate measurements, while placing lower weight on fitting values that showed greater differences between the replicates. In general, fitting averaged data without individually weighting the averaged values should be avoided, because it ignores information contained in the data set about the reliability of each measurement. While arguments can be made about the relative statistical validity of averaging or not averaging replicate data points before curve fitting, if the data are of good quality (i.e. reasonably accurate and precise) then the curve fitting will return very similar parameter values regardless of whether the data are averaged and how they are weighted. Conversely, if the data quality is poor, the parameter values resulting from the curve fit will be unreliable regardless of the fitting approach used. One lesson arising from the above discussion is that it is important to check the default settings for nonlinear regression used by the curve fitting software to determine how and whether the program weights the different data points in the set. A second point is that, before averaging replicates, it is important to examine the results to check that there is no evidence of systematic error. For example, if data from a plate reader assay shows a significant difference between replicate measurements depending on which row or column of the assay plate the samples occupy, then this is indicative of a systematic rather than a random variation in results. Another common example is the observation that replicate measurements made at various times over the course of a day increase or decrease systematically, indicating that one or other of the reagents is not stable. If systematic error is evident, it is not appropriate to average the measurements. Instead, the source of the systematic error should be identified and eliminated, to allow a valid experiment to be performed. Some common errors in data fitting involve the generation of best-fit parameter values that fall outside scientifically reasonable limits. For example, when fitting a set of inhibitor concentration-response data to an inhibition curve, it is common for investigators who are inexperienced in quantitative analysis to generate curve fits that extrapolate to a reaction velocity below zero at saturating inhibitor concentrations, or to a velocity above that of the uninhibited reaction when extrapolated to zero inhibitor. These outcomes are physically impossible, and it is incumbent upon the investigator to make sure that the parameter values they obtain after curve fitting their data are indeed reasonable. One approach that can help avoid these errors is to plot the curve fit for a range of x-variable values that extend for at least 2-3 logs above and below the IC50, so that any non-physical behavior at these extremes of inhibitor concentration become apparent. It is often appropriate to constrain the values of one or more parameters in the fitting equation to fix the minimum or maximum y-variable values to known limits or control values. Thus, it is necessary for the investigator to understand the fitting equation (though not necessarily the algorithm by which the data are fit to this equation), and especially that they know which parameter in the equation corresponds to which feature of the fitted curve.



Page 6 of 30

What to report: Standard errors should be reported for each parameter and some curve fitting programs will also calculate confidence intervals which could also be reported. If parameters such as IC50 values are determined for multiple reagent preps (inhibitor batches/protein preps), then the mean of the values can be reported together with the standard deviation and the number of replicates. In plotting fitted data for publication, the experimental data should be plotted on top of the theoretical curve predicted from the equation using the best-fit parameter values. If the data are averaged, then error bars should be included that represent the range or the standard deviation of the replicate measurements that contributed to each average data point. If the experimental data do not extend to sufficiently low and high values of the x-variable to approach zero effect and a full effect, the line that represents the curve fit should not end at the lowest and highest experimental x-variable values, but should extend beyond the data to show that the best-fit equation is consistent with the highest and lowest y-variable values (e.g. full enzyme activity versus zero or background levels of activity) expected from the experiment. It should be emphasized that these recommendations only scratch the surface of the topic; there is a great deal of additional complexity in curve fitting. For example, curve fitting programs will also calculate a 2 or R2, which are measures of how well the regression fits the actual data (the Goodness of Fit). These parameters assess the variation of the actual data from the fitted curve and in general a R2 value close to 1 or a 2 value close 0 are taken as evidence that the data are fit well by the model. However, caution should be exercised in using these values since fitting data to the wrong model can still yield what appears to be a statistically ‘good fit’. For those interested in learning more about the intricacies of curve fitting there are many excellent publications some of which are referenced here.14,1820

Enzyme Assays Steady state enzyme kinetics. Many drug targets are enzymes, and hence SAR is commonly based on assays in which the effect of compounds on the rate of substrate consumption or product formation is monitored. Consequently, it is worth reminding ourselves of approximations used in deriving the Michaelis Menten equation (Equation 1) which is based on the simple reaction scheme shown in Figure 1, viz (i) [S] >> [E] so that [S]free  [S]total and (ii) initial velocity is determined in the first few % of the reaction so that [S] is essentially constant.

E+S

k1 k2

ES

kcat

E+P

Figure 1: Kinetic scheme used to derive the Michaelis-Menten equation. k1 is the rate constant for formation of ES from E+S, k2 is the rate constant for the dissociation of ES back to E + S, and kcat is the turnover number for the enzyme.

𝑉𝑚𝑎𝑥[𝑆]

1

𝑣𝑖 = 𝐾𝑚 + [𝑆]


Page 7 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


The Michaelis-Menten equation can be derived using the steady state assumption in which the concentration of ES is assumed to be constant (cf. Briggs-Haldane derivation) and in Equation 1 Vmax = kcat [E]total, where [E]total is the total enzyme concentration and Km = (k2+kcat)/k1. Km is the MichaelisMenten constant and represents the dissociation constant of all enzyme bound species.

Basic assumptions and concentration (dose)-response curves. Similar rate expressions can be derived to account for the impact of enzyme inhibitor on the rate of the reaction. In inhibition assays it is also often assumed that [I] >> [E] again so that [I]free  [I]total. This is because equilibrium constants are determined from the concentration of reactants present at equilibrium (i.e. [S]free or [I]free) and given the above approximations data analysis can use the total concentration of substrate or inhibitor added to the reaction. Enzyme inhibition is most commonly analyzed using initial velocity measurements. Formally the initial velocity is determined from a tangent to the initial portion of the reaction progress curve, and it is most convenient to use a continuous assay format so that the initial velocity can be accurately measured. However, many assays, including those run in high throughput, often rely on single time point or end point assays from which initial velocity data are extracted. In this case it is important to recognize the underlying assumption that the reaction velocity is linear until the single time point is taken. If a continuous assay format is not available, we recommend that multiple time points are taken in order to check the linearity of the reaction. One simple control is to demonstrate that vi  [E], and thus the initial velocity should double if the enzyme concentration is doubled. The change in the initial velocity as a function of inhibitor concentration is then often fit to a concentration-response relationship to obtain the IC50 value for the inhibitor. Software programs, such as GraphPad Prism and GraFit, include several forms of the standard IC50 equation that differ depending on whether the response increases or decreases with inhibitor concentration, and also whether the data are fit to a two or four parameter equation (Equations 2 and 3, respectively).

100

% 𝐼𝑛ℎ𝑖𝑏𝑖𝑡𝑖𝑜𝑛 = 1+

𝑌 = 𝑌𝑚𝑖𝑛 +

2

𝐼𝐶50 ℎ [𝐼]

( )

𝑌𝑚𝑎𝑥 ― 𝑌𝑚𝑖𝑛 1+

3

𝐼𝐶50 ℎ [𝐼]

( )

In Equations 2 and 3, where the response is assumed to increase with [I], [I] is [I]free, and h is the Hill coefficient or slope factor. Most programs set h as a variable in addition to IC50. However, in a simple binding equilibrium, where there is no cooperativity and only 1 binding site, h is expected to be 1. It is thus recommended that h is set to 1 during data fitting, and then allowed to float only if there is clear indication that the experimental data cannot be adequately modeled using h=1, for example if there is 7 ACS Paragon Plus Environment


an obvious systematic deviation between the fitted curve and the experimental data. There are a number of mechanistic and artifactual reasons that can lead to values of h that deviate from unity,12 for example if there is positive (h>1) or negative (h>[E], have been fit to an equation in which the responses increases with [I] (Equation 2) to give an IC50 of 114 nM. The data have been converted into % inhibition where the response changes from 0 to 100% inhibition over the experiment so that only 1 parameter is needed for initial data fitting, constraining the slope factor (Hill coefficient, h) to be 1. In a 2-parameter fit, h would 8 ACS Paragon Plus Environment

Page 8 of 30

Page 9 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


be allowed to vary whilst in a 4-parameter fit the range over which the response varies (Ymax-Ymin) together with a background signal (Ymin) are also variables (Equation 3). Also shown are the calculated fits if h is constrained to 2 or 0.5 where it can be seen that there is a systematic deviation between the fitted curve and the experimental data points.

If there is no prior knowledge of compound potency, then it is not possible to choose appropriate inhibitor concentrations before any measurements have been made. In this case it is recommended that the extent of inhibition is assessed using several fixed inhibitor concentrations that vary by factors of 10 – say 0.1, 1 and 10 M. Once it is possible to estimate concentrations that span a response from ~10 to 90% inhibition, it is then important to recognize that the level of inhibition is not a linear function of [I] (Figure 2), and thus inhibitor concentrations should be chosen based on a logarithmic rather than a linear scale to ensure an even distribution of data points across the plot. On a log scale, 3.16228 is halfway between 1 and 10, not 5. In addition, there should be enough data points at high and low concentrations to clearly define the end points for curve fitting. Finally, the pIC50 (-log10(IC50)) rather than the IC50 is sometimes reported to account for the exponential nature of the relationship.23 It is important to realize that the IC50 value is a function of both the enzyme and substrate concentration. Thus, caution should be exercised when comparing IC50 values for the inhibition of a specific enzyme by different labs even when determined under ‘identical’ conditions. Ki values which do not depend on [E] or [S] offer a more rigorous basis for comparison.

Coupled assays. Continuous assays often follow the change in a spectroscopic signal, such as the absorbance or fluorescence of the substrate or product, as a function of time. However, in situations where neither substrate nor product have a convenient spectroscopic signature, the enzyme reaction of interest may be coupled to a second reaction that does lead to change in signal that can be easily monitored. In this case it is important to show that the coupling reaction is not rate limiting and that the inhibitor does not affect the activity of the coupling enzyme. In addition, it is often not realized that there may be an initial lag in reaction velocity before the initial steady-state velocity is reached. The lag depends on factors such as the kcat and Km of the coupling enzyme, and can be minimized by adjusting the concentration of the coupling enzyme.24 If it is not possible to eliminate the lag phase by adjusting the assay conditions, then initial velocities should be obtained after ensuring that the lag is complete.

Ki instead of IC50. While concentration-response curves represent the bulk of activity measurements in SAR, enzyme inhibition can also be quantified by obtaining Ki values. This can be accomplished by determining kcat and Km values at different fixed inhibitor concentrations, and using this information to generate Ki values, together with information on the mechanism of inhibition: competitive (Ki), noncompetitive/mixed (Ki and Ki’) and uncompetitive (Ki’) (Figure 3). Although this is more laborious 9 ACS Paragon Plus Environment


Page 10 of 30

than an IC50 measurement, Ki values enable a better comparison of inhibitor potency between compounds since IC50 values depend on the precise experimental conditions that are used, such as substrate concentration, and provide no mechanistic information about the mode of inhibition. Ki values can be calculated from IC50 values using the Cheng-Prusoff relationship,25 however this requires knowledge of the mechanism of inhibition as well as the ratio of [S]/Km. The mechanistic information derived from Ki measurements is useful since not all compounds are competitive inhibitors (although this is often assumed in the absence of any mechanistic data). In addition, the observation of pure noncompetitive inhibition (Ki = Ki’) may be indicative of a promiscuous inhibitor since it is highly unlikely that an inhibitor has equal affinity for both E and ES.21,26 Furthermore, uncompetitive inhibition is often considered an advantage since the increase in substrate concentration caused by target inhibition will increase the level of inhibition rather than competing off the inhibitor.

E+S Ki

I EI

k1 k2

ES Ki'

kcat

E+P

I

ESI

Figure 3: Competitive, Noncompetitive/Mixed and Uncompetitive Inhibition. A competitive inhibitor binds to free enzyme and competes with the substrate while an uncompetitive inhibitor binds to the ES complex (binds after the substrate). A mixed inhibitor binds to both E and ES, while a pure noncompetitive inhibitor has equal affinity for E and ES (Ki = Ki’).

Complexities: Tight binding inhibition. Tight binding inhibition occurs when the Ki value is similar to, or below, the enzyme concentration (Ki/[E]T < 10), leading to depletion of the free inhibitor in the assay so that the assumption [I]free  [I]T is no longer valid. It is thus an experimental definition since it is based on the lowest enzyme concentration that can be used in the assay. Clearly this is a good problem to have, since we generally want compounds that are very potent. To account for inhibitor depletion under tight binding conditions, inhibition is quantified using the Morrison equation which includes the total inhibitor concentration [I]T rather than the free concentration of inhibitor present at equilibrium ([I]free) (Equation 4).

𝑣𝑖 𝑣0

=1―

𝑎𝑝𝑝 2 ([𝐸]𝑇 + [𝐼]𝑇 + 𝐾𝑎𝑝𝑝 𝑖 ) ― ([𝐸]𝑇 + [𝐼]𝑇 + 𝐾𝑖 ) ― 4[𝐸]𝑇[𝐼]𝑇

2[𝐸]𝑇

4

For tight binding inhibitors, IC50 values can seriously underestimate potency and caution should be exercised if IC50 values start approaching [E]: note that the IC50 cannot be smaller than ½[E]. Finally, If IC50 < ½[E] is observed, then this could be because the concentration of active enzyme in the reaction 10 ACS Paragon Plus Environment

Page 11 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


has been overestimated. Indeed, it should be noted that use of the Morrison equation to determine Kiapp requires an accurate value for the enzyme concentration, unlike when [I]>>[E]. If tight-binding is suspected, then IC50 values should be determined at more than one enzyme concentration: If [I]>>[E] then the IC50 value will not change if the enzyme concentration is altered, however the IC50 value will change for a tight binding inhibitor and indeed plotting IC50 against [E] should yield a straight line with a Y-intercept where IC50 = Kiapp.

Complexities: Slow binding inhibition. In addition to assumptions such as [S] >> [E] and [I] >> [E], it is also often assumed that the system rapidly reaches equilibrium (i.e. in the mixing time of the experiment). However, some compounds are slow binding inhibitors in which the steady state is reached slowly on the timescale of the assay that is used. Slow-binding inhibitors will also dissociate slowly from their targets, which has important pharmacological consequences since the time scale for breakdown of the drug-target complex can be on the same time scale as the rate of drug elimination.2729 Slow binding inhibition can occur via several mechanisms, however the most commonly encountered are a one-step mechanism, or a two-step mechanism in which the rapid formation of the initial enzymeinhibitor complex is followed by a slow step leading to the final enzyme-inhibitor complex (Figure 4).

A E+I

k3 k4

EI

B E+I

k3 k4

EI

k5 k6

EI*

Figure 4: Slow binding inhibition mechanisms. A. A one-step mechanism. B. A two-step mechanism in which the rapid formation of the initial EI complex is followed by a slow step leading to the final EI* complex. These are mechanisms A and B from Morrison and Walsh.30 Note that by convention inhibition rate constants are numbered starting with k3 since k1 and k2 are used to describe substrate binding (Figure 1).

Slow binding inhibition can be diagnosed by observing curvature in the enzyme assay under conditions where the uninhibited reaction is linear, and data are often analyzed by forward or reverse (jump dilution) progress curves. Figure 5 contains forward progress curve data for the inhibition of an enzyme by a time-dependent inhibitor. The progress curves are fit to Equation 5 giving values for vi, vs and kobs as a function of [I], where vi is the initial velocity in the presence of the inhibitor, vs is the final steady state velocity and kobs is the rate constant for onset of inhibition. The dependence of kobs on [I] can be used to distinguish between one and two-step mechanisms: in this case the plot is hyperbolic, consistent 11 ACS Paragon Plus Environment


Page 12 of 30

with an induced-fit two-step mechanism (Mechanism B), and a fit of the data to Equation 6 gives values for k5 and k6 (Figure 4) as well as Kiapp, the dissociation constant of EI. Kiapp can also be extracted directly from a plot of vi/v0 against [I], (Equation 7), where v0 is the rate in the absence of the inhibitor, while Ki*app, the overall dissociation constant, can obtained by fitting the dependence of vs/v0 on [I] to Equation 8. In a one-step mechanism, vi does not vary with [I].31 Note that in these equations we use Kiapp and Ki*app, the apparent values for the equilibrium dissociation constants. This is because progress curve analysis is performed at a single substrate concentration, which is often very high to ensure that the velocity in the absence of inhibitor is linear. Thus, the experimentally measured dissociation constant will depend on the amount of enzyme that is present as ES. The true values for the equilibrium constants can be obtained from the apparent values using the Cheng-Prusoff equation if both the mechanism of inhibition is known as well as the value for the substrate Km. For example, if [S] > Km then Kiapp will be larger than Ki by a ratio of 1+[S]/Km for a competitive inhibitor. Analogous correction factors are also available for non-competitive and uncompetitive inhibitors.

[𝑃] = 𝑣𝑠𝑡 +

𝑣𝑖 ― 𝑣𝑠

𝑘𝑜𝑏𝑠 = 𝑘6 +(1 +

𝑣𝑖 𝑣0

𝑣𝑠 𝑣0

=

=

(1 ― 𝑒( ― 𝑘𝑜𝑏𝑠𝑡))

𝑘𝑜𝑏𝑠

𝑘5 𝐾𝑎𝑝𝑝 𝑖 /[𝐼]

5

6

)

1

7

1 + ([𝐼]/𝐾𝑎𝑝𝑝 𝐼 )

1

8

1 + ([𝐼]/𝐾𝐼∗ 𝑎𝑝𝑝)


Page 13 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


Figure 5: Progress curve analysis of a two-step slow-binding inhibitor. Under conditions where the reaction velocity is linear in the absence of inhibitor (v0), curvature in the presence of inhibitor is diagnostic of slow binding inhibition.30 The figure shows forward progress curve analysis for the inhibition of an enzyme simulated using Kintek,32 which follows a two-step induced fit mechanism in which the rapid formation of the initial enzyme-inhibitor complex (EI) is followed by the slow isomerization of EI to EI*. A. Fitting of the data to the progress curve equation (Equation 5) yields values for vi, the initial velocity, vs, the final steady-state velocity, and kobs, the rate constant for formation of the steady state. B. The hyperbolic dependence of kobs on [I] is consistent with the twostep induced-fit mechanism and fitting of the data to Equation 6 gives k5= 0.34 min-1, k6 = 0.029 min-1, and Kiapp = 0.34 M. C. Consistent with a two-step mechanism, vi varies with [I], and a fit of the data to Equation 7 also gives a value for Kiapp. D. A fit of vs/v0 against [I] to Equation 8 gives Ki*app =0.026 M. Several challenges exist in analyzing slow binding inhibitors. For example, it may be difficult, due to issues such as substrate solubility, cost or availability, to achieve reaction conditions where the initial velocity in the absence of inhibitor is linear for sufficient time to observe curvature in the presence of inhibitor. In addition, it may be difficult to distinguish one and two-step binding mechanisms. For example, if the initial binding event in a two-step mechanism has a relatively high Kiapp, then it may be difficult to reach sufficient inhibitor concentration to observe curvature in the kobs v [I] plot. In this case, a two-step binding mechanism can give the appearance of a one-step mechanism. Note that for twostep binding vi is also expected to vary with [I] (Figure 5C), however this may also be difficult to detect for the same reason. Finally, if the rate of dissociation of inhibitor from the enzyme is very slow, then it may be difficult to distinguish reversible from irreversible inhibition since the intercept on the Y axis of 13 ACS Paragon Plus Environment


the kobs v [I] plot will be close to 0. For reversible inhibitors with slow off rates it may be easier to perform a jump dilution (reverse progress curve) assay, in which enzyme and inhibitor are preincubated for several hours before dilution into a reaction mixture containing substrate. In this case the reaction velocity will increase as the inhibitor slowly dissociates from enzyme leading to a reverse progress curve. For slowly dissociating inhibitors it may also be possible to prepare and purify the enzyme-inhibitor complex and then directly monitoring the rate of inhibitor dissociation after diluting the complex. For high affinity inhibitors it will be necessary to use a very sensitive method to quantify free inhibitor concentration such as mass spectrometry or radioactivity.33 Several examples of the above scenarios are shown in Figure 6.

Figure 6: Examples of slow-binding inhibition analysis. A. kobs v [I] plot for the slow-binding inhibition of acetylcholinesterase by C547 which follows a two-step binding mechanism. Adapted with permission from Kharlamova, A. D., Lushchekina, S. V., Petrov, K. A., Kots, E. D., Nachon, F., Villard-Wandhammer, M., Zueva, I. V., Krejci, E., Reznik, V. S., Zobov, V. V., Nikolsky, E. E., and Masson, P. (2016) Slow-binding inhibition of acetylcholinesterase by an alkylammonium derivative of 6-methyluracil: mechanism and possible advantages for myasthenia gravis treatment, Biochem J 473, 1225-1236. DOI 10.1042/BCJ20160084. Copyright 2016 The Biochemical Society.34 B. kobs v [I] plot for the slow-binding inhibition of rhomboid protease GlpG by peptidyl ketoamide compound 10 which follows a one-step 14 ACS Paragon Plus Environment

Page 14 of 30

Page 15 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


binding mechanism. Adapted with permission from Ticha, A., Stanchev, S., Vinothkumar, K. R., Mikles, D. C., Pachl, P., Began, J., Skerle, J., Svehlova, K., Nguyen, M. T. N., Verhelst, S. H. L., Johnson, D. C., Bachovchin, D. A., Lepsik, M., Majer, P., and Strisovsky, K. (2017) General and Modular Strategy for Designing Potent, Selective, and Pharmacologically Compliant Inhibitors of Rhomboid Proteases, Cell Chem Biol 24, 1523-1536 e1524. DOI 10.1016/j.chembiol.2017.09.007. Available under the terms of the Creative Commons Attribution License (CC BY).35 C. kobs v [I] plot for the slow-binding inhibition of polypeptide deformylase by actinonin which follows a two-step binding mechanism. In this case the intercept on the Y-axis is close to 0, and so the data have been analyzed using a modified version of Equation 6 where k6 is set to 0. In this case k5 and Kiapp are actually kinact and KI, which are the parameters for quantifying irreversible enzyme inactivation (see below). Adapted with permission from Van Aller, G. S., Nandigama, R., Petit, C. M., DeWolf, W. E., Jr., Quinn, C. J., Aubart, K. M., Zalacain, M., Christensen, S. B., Copeland, R. A., and Lai, Z. (2005) Mechanism of time-dependent inhibition of polypeptide deformylase by actinonin, Biochemistry 44, 253-260. DOI 10.1021/bi048632b. Copyright 2005 American Chemical Society36 D. Jump-dilution progress curve analysis for the inhibition of LpxC by six inhibitors (A-F). Adapted with permission from Walkup, G. K., You, Z., Ross, P. L., Allen, E. K., Daryaee, F., Hale, M. R., O'Donnell, J., Ehmann, D. E., Schuck, V. J., Buurman, E. T., Choy, A. L., Hajec, L., MurphyBenenato, K., Marone, V., Patey, S. A., Grosser, L. A., Johnstone, M., Walker, S. G., Tonge, P. J., and Fisher, S. L. (2015) Translating slow-binding inhibition kinetics into cellular and in vivo effects, Nat Chem Biol 11, 416-423. DOI 10.1038/nchembio.1796. Copyright © 2015, Springer Nature.2

Complexities: Irreversible inhibition. Irreversible inhibitors are compounds that bind to the target and do not dissociate. As noted above, some slow binding inhibitors, including those that form a covalent complex with the target as well as those that are non-covalently bound, can appear to be irreversible if the off rate is very slow. In this situation, it is important to perform additional assays to determine whether the compounds are indeed reversible, such as jump dilution assays to check for reactivation of the enzyme. While it is tempting to use IC50 measurements to quantify the ‘potency’ of irreversible inhibitors, such experiments are flawed since they cannot account for the time-dependence of inhibition and virtually any IC50 can be achieved if the reaction is allowed to incubate for long enough. Instead, it is more appropriate to determine kinact/KI, which is the second order rate constant for formation of the irreversible enzyme-inhibitor complex. Irreversible inhibitors often operate through a two-step mechanism in which the inhibitor first binds reversibly to the enzyme followed by covalent bond formation (Figure 7).

E+I

k3 k4

EI

k5

EI*

Figure 7: Two step mechanism for irreversible inhibition. Reversible formation of the initial EI complex is followed by a second irreversible step leading to the covalent enzyme-inhibitor complex EI*. The kinetic mechanism is analogous to the reversible two step mechanism in Figure 4 except that k6 = 0. Irreversible inhibition is normally quantified by kinact/KI, the second order rate constant for the formation of EI*, where KI is the concentration of inhibitor required to reach half-maximal rate of inactivation of enzyme, and kinact is the maximum rate of inactivation at saturating inhibitor concentrations. Note that 15 ACS Paragon Plus Environment


Page 16 of 30

KI is not the same as Ki, the equilibrium constant for dissociation of EI where Ki = k4/k3. While KI can be numerically equal to Ki (e.g. when k4 >> k5), this is often not the case and kinact/KI should be used to quantify inhibitor potency.

Analysis of irreversible inhibition follows a similar approach to that used for slow binding inhibitors except that a simplified progress curve equation can be used since k6 = 0 and vs = 0 (Equation 9). Simulated data are shown in Figure 8 where it can be seen that the appearance of the plots is very similar to those for two-step slow binding reversible inhibition except that the intercept of the kobs v [I] plot passes through 0. In addition, vs values will have a finite value for reversible inhibitors whereas vs values obtained from forward progress curves should tend to 0.

Figure 8: Irreversible inhibition. The kinetics for irreversible enzyme inhibition quantified by progress curve analysis. A. Time-dependent enzyme inactivation as a function of inhibitor concentration (M) has been analyzed using a simplified version of the progress curve equation (Equation 9) since k6 = 0 which yields values of kobs and vi as described in Figure 5.12 B. A plot of kobs v [I] is hyperbolic, consistent with a two-step mechanism in which the initial non-covalent association of the inhibitor with the enzyme is followed by a second-step leading to formation of the final covalent enzyme-inhibitor complex. Fitting of the data to Equation 10 yields values for kinact = 0.033 min-1, KIapp =0.28 M, and kinact/Kiapp = 0.12 µM min-1. Again, by analogy to equations for slow binding inhibitors, Equation 10 includes KIapp, the apparent value for KI, since the presence of substrate will affect the concentration of inhibitor required to reach ½kinact.Time-dependent enzyme inactivation was simulated using Kintek.32

𝑣𝑖

[𝑃] = 𝑘𝑜𝑏𝑠 (1 ― 𝑒

𝑘𝑜𝑏𝑠 =

( ― 𝑘𝑜𝑏𝑠𝑡)

)

9

𝑘𝑖𝑛𝑎𝑐𝑡

10

1 + (𝐾𝑎𝑝𝑝 𝐼 /[𝐼])


Page 17 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


If the initial binding event is weak, then it may not be possible to saturate the enzyme in order to determine values for kinact and KIapp individually. In this case the plot of kobs v [I] will be linear with a slope of kinact/Kiapp. An interesting recent example of this behavior is the covalent inhibition of KRASG12C by compounds that bind weakly (Ki > 64 M) but react rapidly with C12 (kinact > 0.019 s-1).37

Direct Binding Assays Binding kinetics can also be determined using approaches that directly measure the formation and breakdown of the drug-target complex. These assays are used when no functional assay is available, or when binding of ligands to the free enzyme is being assessed. Numerous methods are available, including the use of radiolabeled compounds and mass spectrometry, which provide the highest levels of sensitivity, to approaches in which there is a change in spectroscopic signal on binding. In many cases, the analysis of data from direct-binding assays depends on the same fundamental principles and assumptions described for enzyme assays. Below I briefly discuss several methods that are in common use.

Radioligand binding assays. Receptor-ligand binding assays have historically relied on the use of radiolabeled ligands which permit binding to be measured for purified receptors or receptors on cell surfaces.38,39 In these assays the concentration of bound ligand is often determined directly after free ligand is removed from an immobilized receptor preparation using a wash step (often performed with cold buffer). This approach assumes that the wash step is fast relative to the rate of ligand dissociation. In addition, it is important to vary incubation times to ensure that the system is at equilibrium and also to account for non-specific binding. Data analysis often assumes that the concentration of receptor ([R]) is very low so that [L]>>[R] at every concentration of ligand ([L]) and thus that [L]free  [L]total. However, if Kd/[R] < 10 then ligand depletion will occur, and the quadratic Morrison equation must be used, with the caveat that the concentration of receptor is accurately known. Finally, a simple one-site binding event should have a Hill slope (h) of 1. Once a radioligand binding assay has been developed then the binding of additional ligands can be measured using a competitive radioligand binding assay. This can be performed under equilibrium conditions to obtain the Kd for binding but can also be adapted, as described by Motulsky and Mahan, to provide the on and off-rates for ligand binding.1,40

Fluorescence binding assays. Many direct binding assays are based on a change in spectroscopic signal upon formation of the receptor-ligand complex. This includes methods in which the fluorescence of the protein and/or ligand is altered. The change in fluorescence is determined as a function of ligand concentration, and the Kd for the ligand is obtained by fitting the data to a binding isotherm using 17 ACS Paragon Plus Environment


equations described above where the same assumptions, such as [L]>>[R], may apply. Again, it is suggested that h is set to 1 for initial data fitting. Fluorescence-based methods are generally less sensitive than techniques based on radioactivity or mass spectrometry, with a typical detection limit of ~10-8 M, although there are examples of assays where fluorescence and radiolabeling perform similarly.3 Important experimental considerations include the need to account for the inner filter effect, in which the fluorophore is sufficiently concentrated to attenuate the excitation light and may also re-absorb some of the emitted light. In addition, while binding assays with radiotracers of known specific activity can give the absolute concentration of the receptor, the relative change in fluorescence intensity at saturation often cannot be directly related to concentration. If there is no change in ligand fluorescence upon binding, the anisotropy or polarization may still be affected on complex formation due to a decrease in the rate of tumbling. Indeed, a common approach to the development of a binding assay is to append a fluorophore onto a ligand in a way that is not anticipated to perturb binding, and then to use steady-state and time-resolved anisotropy measurements to quantify binding. A second common fluorescence-based assay is to attach a second fluorophore to the protein and use Förster/fluorescence resonance energy transfer (FRET) to monitor binding either using a steady-state or time-resolved detection.41 FRET and methods such as bioluminescence resonance energy transfer (BRET) have been translated to quantifying drug-target interactions in living cells, such as by tagging the target with NanoLuc luciferase and monitoring BRET to a fluorescent ligand (NanoBRET).7,42 Using a fluorescent competitor ligand, NanoBRET and other proximity-based assays can be used to assess binding of unlabeled ligands in an analogous approach to the competitive radioligand binding assay described above. When a competition assay is used, it is advisable to accurately determine the binding of the competitor so that the kinetic parameters for association and dissociation of the competitor are not also variables in fitting the competition binding data.

Bioluminescent reporter assays In addition to bioluminescence resonance energy transfer (BRET), luciferases have found wide application in cell-based genetic reporter assays in which regulatory elements that control gene transcription are coupled to luciferase gene expression. This approach has proved particularly useful for drug targets such as GPCRs and nuclear hormone receptors which regulate gene transcription. Andruska et al. report an HTS to discover inhibitors of 17β-estradiol (E2)-estrogen receptor α (ERα) induced breast cancer cell proliferation,43 using a luciferase reporter whose expression is driven by three copies of the consensus estrogen response element (ERE)3. The authors discuss compounds that show activity in the primary screen due to a direct effect on luciferase or that are broadly cytotoxic, in addition to compounds that affect cell proliferation through inhibition of ERα

Surface plasmon resonance (SPR) and isothermal titration calorimetry (ITC). Biophysical methods play a major role in characterizing the structure, kinetics and thermodynamics of biomolecular interactions. 18 ACS Paragon Plus Environment

Page 18 of 30

Page 19 of 30

While techniques such as X-ray crystallography and NMR spectroscopy are primarily used to provide insight into the structure of protein-ligand complexes, surface plasmon resonance (SPR) and isothermal titration calorimetry (ITC) are both label free approaches to quantifying protein-protein and proteinligand interactions.44 There are many excellent papers on the application of both SPR and ITC to the analysis of biomolecular interactions and I only mention them here to draw the reader’s attention to these methods. In SPR one of the binding components is tethered to a surface and the change in SPR signal caused by the interaction with the second component yields both the thermodynamic and kinetic parameters for the interaction.9 An example of a protein-protein interaction quantified by SPR is shown in Figure 9.45 In addition to being label free, SPR can be deployed in relatively high throughput and is heavily used in the pharmaceutical industry for drug discovery campaigns. Experimental caveats include the ability to stably immobilize one of the binding partners and to then regenerate the surface after each binding reaction, and mass transfer effects that can mask the actual binding reaction, for example if the transfer of the analyte from bulk solution to the sensor surface is slower than the binding event.46,47 In addition, the tethered protein can also have different properties compared to the untethered protein in solution which may alter ligand binding. HTS libraries invariably contain compounds, such as promiscuous aggregators, that give false hits due to assay interference,26 and Giannetti et al. give clear examples of sensorgrams that arise from nonideal binding events caused by such compounds.48 One interesting recent advance is the development of a ‘chaser’ method that improves the ability of SPR to accurately quantify the kinetics of slowly dissociating ligands and thus overcome problems associated with signal drift.8

3500 Response (RU)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


A

2600 1700 800 -100

0

400

800 Time (s)

1200

1600

Figure 9: SPR Sensorgram. Biacore analysis of soluble monomeric TNF family receptor Fn14 binding to the TNF family ligand TWEAK. A. Sensorgram for 40−360 nM soluble receptor binding to surface derivatized with TWEAK. The experimental data (black lines) were globally fit to a 1:1 binding model (red lines) using BIAevaluation to determine kinetic rate constants. B. The signal observed at equilibrium, Req, plotted as a function of soluble receptor concentration, fit to a hyperbolic, single-site binding equation. Adapted from Day, E. S., Cote, S. M., and Whitty, A. (2012) Binding efficiency of protein-protein complexes, Biochemistry 51, 9124-9136. DOI 10.1021/bi301039t. Copyright © 2012, American Chemical Society.45



ITC is a solution-based method in which the heat liberated (exothermic) or absorbed (endothermic) by the binding event is used to measure thermodynamic binding parameters such as Kd, ΔH and ΔS (Figure 10).49 In contrast to SPR, ITC measures the binding reaction in solution, thus avoiding immobilization of a binding partner, but also usually requires larger amounts of sample and is lower throughput. One important aspect of ITC assay design is to minimize heat changes caused by mixing the solutions that contain the protein (cell) and analyte (syringe), and thus it is desirable to dissolve both components in identical buffer solutions. It is also important to ensure that there is sufficient time between each injection so that the system can come to equilibrium. Errors in concentration of the protein will impact the stoichiometry (n), whilst errors in the concentration of the analyte will affect n, Kd, and ΔH.50 In addition, concentrations should be chosen so that nx[protein]/Kd = 10-50, however meeting this requirement for either very high or low affinity ligands may be challenging. For example, high affinity ligands will require low protein concentrations which will result in small changes in the amount of heat absorbed or released for each injection which may be difficult to measure. Conversely, weak binders will require high concentrations of both protein and analyte which may be difficult to attain due to solubility and/or availability of reagents. Whilst ITC can quantify interactions over a range of affinities from nM to sub-mM, this range depends on the absolute heat change and can be extended by using a competitive ligand.51 Finally, although ITC is used almost exclusively for equilibrium binding measurements, there are reports that this method can also yield the kinetic parameters for binding.52,53

Figure 10: Binding of Biotin Protein Ligase A (BirA) and Biotin Analyzed by ITC. A: Raw calorimetric data obtained by the titration of Biotin (900 M) into wild type BirA (50 µM) at 25°C in 10 mM Tris-HCl buffer (pH 7.5) containing 30 mM NaCl, 200 mM KCl, and 2.5 mM MgCl2. Inset: Enlarged area of the titration curve showing that the system is at equilibrium before each injection. B: The integrated heats of injection plotted as a function of the Biotin/BirA molar ratio, fit to a single binding site model using Microcal Origin® software after correcting for the heat of dilution to give n = 0.981 ± 0.007, Ka = 3.6x105 ± 2.0x104 M-1 (Kd = 2.8 M), H = -12.2 ± 0.1 kcal/mol and S = -15.3 cal/mol/deg. Adapted from Bockman, M. R.; Engelhart, C. A.; Dawadi, S.; Larson, P.; Tiwari, D.; Ferguson, D. M.; Schnappinger, D.; Aldrich, C. C. (2018) Avoiding Antibiotic Inactivation in Mycobacterium tuberculosis by Rv3406 through Strategic Nucleoside Modification ACS Infect Dis, 4, 1102-1113. DOI: 10.1021/acsinfecdis.8b00038. Copyright © 2018, American Chemical Society.54 20 ACS Paragon Plus Environment

Page 20 of 30

Page 21 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


Both methods require specialized equipment and in each case the software that controls the instrument also includes curve fitting programs to analyze the data. Following the theme of this perspective, investigators are urged not to treat the packaged data analysis programs as black boxes for data fitting but instead to understand which equations are being fitted to the data and what, if any, assumptions are made in this procedure. In SPR, curve fitting of the sensorgrams will yield the on and off rates for ligand binding which are then used to calculate a Kd value. However, Kd values can also be calculated directly by plotting the signal at steady state (Response, RU) as a function of ligand concentration provided that the injection time is sufficiently long for the system to reach steady state.55 In ITC the binding isotherm is generated by plotting the amount of heat released or absorbed in each analyte injection as a function of analyte concentration. Curve fitting then yields n (stoichiometry), Ka (association binding constant where Ka = 1/Kd), and ΔH, from which ΔS and ΔG are calculated.56 Initial curve fitting often assumes a 1:1 binding interaction in which case a value of n=1 is expected. However, as noted above N can often deviate from 1 due to inaccuracies in either the protein or ligand concentration.

Summary: Guidelines for publication of quantitative inhibition/binding data In summary, many methods are available for quantifying biomolecular interactions. The choice of method will depend of the specific system under investigation, and in this perspective, I have only attempted to highlight some of the most common approaches in order to draw attention to key experimental considerations in developing assays and analyzing the data that are produced. Instruments such as those used for SPR and ITC measurements come with their own software for analyzing the binding data and while I have not attempted to describe the mathematical models used for these approaches, users should be aware of the underlying assumptions used here too.

A few of the key take home messages are summarized below: General points (i) (ii) (iii) (iv)

Provide full experimental details for each assay including protein and ligand concentrations, buffer conditions, reaction temperature, incubation times and number of replicates. When establishing assays, consider the methods that will be used to analyze the data. In many cases data fitting assumes that [I]free  [I]total. Provide data plots together with the fitted curve(s), and the equation(s) used for the data analysis. Report standard errors for the calculated parameters. Nonlinear regression includes the following steps: choice of model, whether to constrain any parameters, selection of initial values for each parameter, whether to use differential weighting, how to detect and handle outliers, and whether to average replicates before data fitting. In general, there should be at least two or better three replicates at each experimental 21 ACS Paragon Plus Environment


(v) (vi)

set of conditions (e.g. inhibitor concentration). The replicates can be treated as individual data points in curve fitting or the averaged data can be analyzed whilst using the standard deviation of the replicates to weight the data. Fitting averaged data without individually weighting the averaged values should be avoided. Parameters for the final, optimized inhibitor compound(s) should ideally be based on replicates determined from separate preparations of enzyme and inhibitor. Curve fitting programs enable data to be analyzed using very complex mathematical models. Generally, an increase in the number of variables used in data fitting will improve the goodness of fit. However, a valid mechanistic reason must be advanced for increasing the number of variables used to fit the data. This may include information obtained from other approaches. For example, the observation of two different enzyme-inhibitor structures (EI and EI*) by X-ray crystallography supports the two-step slow-onset mechanism for the inhibition of the enoyl-ACP reductase from Mycobacterium tuberculosis revealed by progress curve kinetics.57

Enzyme assays and concentration-response curves (vii)

(viii) (ix)

(x)

(xi)

(xii)

Single time point assays generally assume that the change in activity is linear with time, from the start of the reaction up until the time of measurement. Evidence should be provided that supports this assumption. For coupled assays check that the coupling enzyme is not rate-limiting, and that your inhibitor is not inhibiting the coupling enzyme(s). When fitting data to a concentration-response equation, use a two-parameter fit to obtain IC50 values and initially set h=1. Provide a mechanistic hypothesis if h varies significantly from 1. In addition, if the data are fit using a four-parameter model, discuss the origin of (e.g.) a background signal. In contrast to Ki values, IC50 values do not give information about the mechanism of inhibition. For instance, the observation of pure noncompetitive inhibition (Ki = Ki’) may be indicative of a promiscuous aggregator. IC50 values depend on the precise experimental conditions that are used in the assay, such as enzyme and substrate concentration, and thus cannot be used to compare inhibition of a specific enzyme by different labs. In contrast, Ki values offer a more rigorous basis for comparison. IC50 values can be converted to Ki values using the Cheng-Prusoff relationship if the mechanism of inhibition is known and also the [S]/Km ratio.

Tight binding inhibition (xiii)

(xiv)

If IC50 values start to approach [E], tight binding inhibition may be occurring. In this case, IC50 measurements can seriously underestimate the potency of the compound. A valid IC50 value cannot be less than ½[E]. Tight binding inhibition will lead to a depletion in the concentration of free inhibitor. This can be accounted for using the Morrison equation provided that Ki/[E] > 0.01.


Page 22 of 30

Page 23 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


Slow and irreversible inhibition (xv)

(xvi)

(xvii)

(xviii)

(xix)

(xx) (xxi)

Check for time-dependent effects in binding/inhibition. IC50 shifts following preincubation of enzyme and inhibitor before initiating the reaction by addition of substrate can be used to check for slow-binding inhibition, but this only works for an inhibitor that binds in the absence of the substrate. In addition, the decrease in activity following preincubation could have other explanations, for example that the inhibitor is causing the enzyme to denature. Slow binding inhibition can often be diagnosed by observing curvature in reaction progress curves under conditions where the reaction in the absence of inhibitor is linear. However, it is important to show that curvature does not result from instability of the enzyme under the assay conditions. In addition, linearity may only be accomplished by using high concentrations of substrate, and thus the apparent dissociation constants obtained under the specific assay conditions (Kiapp and Ki*app) may differ significantly from the true values as described by the Cheng-Prusoff equations. Slow binding inhibition can occur by both 1 and 2-step mechanisms but the ability to distinguish these two mechanisms may be affected by inhibitor solubility and/or the dissociation constant of the first step in a 2-step mechanism. For example, a linear increase in kobs as a function of [I] is consistent with a one-step mechanism. However, this plot can also be linear for a two-step mechanism if the highest inhibitor concentration that is used does not saturate the enzyme (i.e. [I] < Kiapp). If the off-rate (k4 or k6) is very slow, then the intercept on the Y-axis of the kobs v [I] plot may be close to 0, and thus it may be difficult to distinguish reversible slow-binding inhibition from irreversible inhibition. Additional methods may then be needed to verify that the inhibition is indeed irreversible and/or to quantify the off-rate for very slowly dissociating compounds. For compounds that bind very slowly, it may be easier to use reverse progress curve analysis (jump dilution), where enzyme and inhibitor are pre-incubated prior to dilution into the reaction solution.2 Tight-binding inhibition may occur with very potent compounds, requiring data analysis that is based on [I]total rather than [I]free as described by the Morrison equation. SKR for irreversible inhibition should use kinact/KI and not IC50 values.

Direct binding assays (xxii) (xxiii)

(xxiv)

Many direct binding assays have the same assumptions as those used for activity measurements, for example that [L]free  [L]total. Competitive radioligand or fluorescence-based assays are often used to assess the binding of other ligands. Generally, it is advisable to first quantify the binding of the competitor on its own in order to reduce the number of variables in fitting data obtained from the competition assay. Data obtained by methods such as SPR or ITC will usually be analyzed by software provided with the instruments. Whilst this is convenient, investigators should determine what assumptions have been made in deriving the mathematical models that are used for data analysis. 23 ACS Paragon Plus Environment


Common Problems (xxv)

The observation of a background rate (drift) before all the reagents have been combined could have several explanations including instrument drift, reagent instability or changes in reagent solubility. It can be useful to scan the UV-vis spectrum of the reagent since precipitation will lead to light scattering which will manifest as an increase in absorption across the whole absorbance spectrum, increasing as the wavelength decreases (since scattering intensity  -4). Proteins and other reagents can also bind to surfaces, leading to changes in the effective concentration in solution. One useful control is to double [E] and show that the observed rate doubles. This is a good check that the change in signal is due to the enzyme catalyzed reaction. (xxvi) Systematic errors in replicates could be due to reagent instability, and it is important to stagger replicates over the whole course of an experiment (e.g. at several times during the day) to check for systematic changes in activity. Proteins with essential Cys residues can lose activity due to oxidation, and reducing reagents (such as DTT) are often added to stock solutions or assays. Controls should be performed to check that the reducing reagent does not react with other reagents. The oxidized forms of some reducing reagents can have absorbance at 280 nm, so the background absorbance of a buffer containing DTT can increase slowly with time. Ligands (inhibitors) can also be unstable, especially if they contain reactive groups and/or aggregate or precipitate either before or after being added to the assay. (xxvii) IC50 values less the ½[E] could occur if the enzyme is less active than assumed, for example if it is not pure or if a fraction of the enzyme is inactive (e.g. due to oxidation – see above). (xxviii) There can be several mechanistic or artefactual explanations for concentration-response plots with Hill coefficients (h) differing from unity. For instance, compounds that form colloidal aggregates will likely inhibit enzymes through a non-specific mode of action which often manifests as a steep response (h > 1). Shoichet and co-workers have suggested a number of approaches to test for promiscuous inhibition.58 (xxix) Tight binding inhibition will also result in concentration-response plots with h > 1 and can have important consequences on the use of IC50 values for establishing SAR. Since IC50  ½[E], the dynamic range of SAR resulting from IC50 values is limited by the enzyme concentration and thus IC50 values can seriously underestimate compound potency (see discussion in Chapter 7, Tight Binding Inhibition, of Copeland12). Assuming that the assay has not bottomed-out, then repeating the IC50 measurement at a lower enzyme concentration can be informative. For instance, if tight binding conditions are present, then reducing [E] will lead to a reduction in the apparent IC50 and a decrease in h values towards 1. This limitation in IC50 values again reinforces the benefit of determining Ki values.


Page 24 of 30

Page 25 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


Abbreviations BRET: Bioluminescence resonance energy transfer FRET: Förster/fluorescence resonance energy transfer h: Hill coefficient IC50: Concentration of inhibitor (ligand) that results in 50% inhibition (effect) ITC: Isothermal titration calorimetry SKR: Structure-kinetic relationship SPR: Surface plasmon resonance

AUTHOR INFORMATION Corresponding Author *E-mail: [email protected]. Tel: (631) 632-7907.

ORCID Peter J. Tonge: 0000-0003-1606-3471

Conflict of Interest The author is a founder of Chronus Pharmaceuticals Inc.

ACKNOWLEDGMENTS This work was supported by grants from the National Institutes of Health (GM102864, EB024549, EB027050), the National Science Foundation (1817837) and a PhRMA Foundation Fellowship.



References (1) Motulsky, H. J.; Mahan, L. C. (1984) The kinetics of competitive radioligand binding predicted by the law of mass action Mol Pharmacol, 25, 1-9. DOI: (2) Walkup, G. K.; You, Z.; Ross, P. L.; Allen, E. K.; Daryaee, F.; Hale, M. R.; O'Donnell, J.; Ehmann, D. E.; Schuck, V. J.; Buurman, E. T.; Choy, A. L.; Hajec, L.; Murphy-Benenato, K.; Marone, V.; Patey, S. A.; Grosser, L. A.; Johnstone, M.; Walker, S. G.; Tonge, P. J.; Fisher, S. L. (2015) Translating slow-binding inhibition kinetics into cellular and in vivo effects Nat Chem Biol, 11, 416-23. DOI: 10.1038/nchembio.1796 (3) Handl, H. L.; Vagner, J.; Yamamura, H. I.; Hruby, V. J.; Gillies, R. J. (2005) Development of a lanthanidebased assay for detection of receptor-ligand interactions at the delta-opioid receptor Anal Biochem, 343, 299-307. DOI: 10.1016/j.ab.2005.05.040 (4) Zwier, J. M.; Roux, T.; Cottet, M.; Durroux, T.; Douzon, S.; Bdioui, S.; Gregor, N.; Bourrier, E.; Oueslati, N.; Nicolas, L.; Tinel, N.; Boisseau, C.; Yverneau, P.; Charrier-Savournin, F.; Fink, M.; Trinquet, E. (2010) A fluorescent ligand-binding alternative using Tag-lite(R) technology J Biomol Screen, 15, 1248-59. DOI: 10.1177/1087057110384611 (5) Wilson, D. J.; Shi, C.; Duckworth, B. P.; Muretta, J. M.; Manjunatha, U.; Sham, Y. Y.; Thomas, D. D.; Aldrich, C. C. (2011) A continuous fluorescence displacement assay for BioA: an enzyme involved in biotin biosynthesis Anal Biochem, 416, 27-38. DOI: 10.1016/j.ab.2011.05.003 (6) Rinken, A.; Lavogina, D.; Kopanchuk, S. (2018) Assays with Detection of Fluorescence Anisotropy: Challenges and Possibilities for Characterizing Ligand Binding to GPCRs Trends Pharmacol Sci, 39, 187-199. DOI: 10.1016/j.tips.2017.10.004 (7) Stoddart, L. A.; Kilpatrick, L. E.; Hill, S. J. (2018) NanoBRET Approaches to Study Ligand Binding to GPCRs and RTKs Trends Pharmacol Sci, 39, 136-147. DOI: 10.1016/j.tips.2017.10.006 (8) Quinn, J. G.; Pitts, K. E.; Steffek, M.; Mulvihill, M. M. (2018) Determination of Affinity and Residence Time of Potent Drug-Target Complexes by Label-free Biosensing J Med Chem, 61, 5154-5161. DOI: 10.1021/acs.jmedchem.7b01829 (9) Cooper, M. A. (2002) Optical biosensors in drug discovery Nat Rev Drug Discov, 1, 515-28. DOI: 10.1038/nrd838 (10) Yu, H.; Dranchak, P.; Li, Z.; MacArthur, R.; Munson, M. S.; Mehzabeen, N.; Baird, N. J.; Battalie, K. P.; Ross, D.; Lovell, S.; Carlow, C. K.; Suga, H.; Inglese, J. (2017) Macrocycle peptides delineate lockedopen inhibition mechanism for microorganism phosphoglycerate mutases Nat Commun, 8, 14932. DOI: 10.1038/ncomms14932 (11) Motulsky, H. J.; Ransnas, L. A. (1987) Fitting Curves to Data Using Nonlinear-Regression - a Practical and Nonmathematical Review FASEB J, 1, 365-374. DOI: 10.1096/fasebj.1.5.3315805 (12) Copeland, R. A. Evaluation of enzyme inhibitors in drug discovery: a guide for medicinal chemists and pharmacologists; 2nd ed.; Wiley: Hoboken, N.J., 2013. (13) In Assay Guidance Manual; Sittampalam, G. S., Coussens, N. P., Brimacombe, K., Grossman, A., Arkin, M., Auld, D., Austin, C., Baell, J., Bejcek, B., Chung, T. D. Y., Dahlin, J. L., Devanaryan, V., Foley, T. L., Glicksman, M., Hall, M. D., Hass, J. V., Inglese, J., Iversen, P. W., Kahl, S. D., Kales, S. C., Lal-Nag, M., Li, Z., McGee, J., McManus, O., Riss, T., Trask, O. J., Jr., Weidner, J. R., Xia, M., Xu, X., Eds.; Eli Lilly & Company and the National Center for Advancing Translational Sciences: Bethesda (MD), 2004. (14) Motulsky, H. J.; Christopoulos, A.; Oxford University Press: New York, NY USA, 2004. (15) Division of Biomedical Research Workforce National Institutes of Health. Requirement for Instruction in the Responsible Conduct of Research. https://grants.nih.gov/grants/guide/notice-files/NOT-OD10-019.html (16) Blainey, P.; Krzywinski, M.; Altman, N. (2014) Points of significance: replication. Nat Methods, 11, 879-880. DOI: 10.1038/nmeth.3091 26 ACS Paragon Plus Environment

Page 26 of 30

Page 27 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


(17) Vaux, D. L.; Fidler, F.; Cumming, G. (2012) Replicates and repeats--what is the difference and is it significant? A brief discussion of statistics and experimental design EMBO Rep, 13, 291-6. DOI: 10.1038/embor.2012.36 (18) Cornish-Bowden, A. (2014) Analysis and interpretation of enzyme kinetic data Perspectives in Science, 1, 121-125. DOI: 10.1016/j.pisc.2014.02.010 (19) Motulsky, H. J. GraphPad Statistics Guide. http://www.graphpad.com/guides/prism/7/statistics/index.htm (20) Vaux, D. L. (2012) Research methods: Know when your numbers are significant Nature, 492, 180-1. DOI: 10.1038/492180a (21) McGovern, S. L.; Caselli, E.; Grigorieff, N.; Shoichet, B. K. (2002) A common mechanism underlying promiscuous inhibitors from virtual and high-throughput screening J Med Chem, 45, 1712-22. DOI: 10.1021/jm010533y (22) Straus, O. H.; Goldstein, A. (1943) Zone behavior of enzymes: illustrated by the effect of dissociation constant and dilution on the system cholinesterase-physostigmine J Gen Physiol, 26, 559-585. DOI: (23) Navare, M. Use pIC50 Instead of IC50 to Think About Your Assays Logarithmically. https://www.collaborativedrug.com/using-pic50-instead-of-ic50-will-change-the-way-you-seeyour-drug-discovery-data/ (24) Copeland, R. A. Enzymes : a practical introduction to structure, mechanism, and data analysis; 2nd ed.; Wiley: New York, 2000. (25) Cheng, Y.; Prusoff, W. H. (1973) Relationship between the inhibition constant (K1) and the concentration of inhibitor which causes 50 per cent inhibition (I50) of an enzymatic reaction Biochem Pharmacol, 22, 3099-108. DOI: 10.1016/0006-2952(73)90196-2 (26) Aldrich, C.; Bertozzi, C.; Georg, G. I.; Kiessling, L.; Lindsley, C.; Liotta, D.; Merz, K. M., Jr.; Schepartz, A.; Wang, S. (2017) The Ecstasy and Agony of Assay Interference Compounds ACS Med Chem Lett, 8, 379-382. DOI: 10.1021/acsmedchemlett.7b00056 (27) Tonge, P. J. (2018) Drug-Target Kinetics in Drug Discovery ACS Chem Neurosci, 9, 29-39. DOI: 10.1021/acschemneuro.7b00185 (28) Copeland, R. A. (2016) The drug-target residence time model: a 10-year retrospective Nat Rev Drug Discov, 15, 87-95. DOI: 10.1038/nrd.2015.18 (29) Swinney, D. C. (2009) The role of binding kinetics in therapeutically useful drug action Curr Opin Drug Discov Devel, 12, 31-9. DOI: (30) Morrison, J. F.; Walsh, C. T. (1988) The behavior and significance of slow-binding enzyme inhibitors Adv Enzymol Relat Areas Mol Biol, 61, 201-301. DOI: 10.1002/9780470123072.ch5 (31) Morrison, J. F. (1982) The Slow-Binding and Slow, Tight-Binding Inhibition of Enzyme-Catalyzed Reactions Trends Biochem Sci, 7, 102-105. DOI: Doi 10.1016/0968-0004(82)90157-8 (32) Johnson, K. A.; Simpson, Z. B.; Blom, T. (2009) Global kinetic explorer: a new computer program for dynamic simulation and fitting of kinetic data Anal Biochem, 387, 20-9. DOI: 10.1016/j.ab.2008.12.024 (33) Yu, W.; Neckles, C.; Chang, A.; Bommineni, G. R.; Spagnuolo, L.; Zhang, Z.; Liu, N.; Lai, C.; Truglio, J.; Tonge, P. J. (2015) A [(32)P]NAD(+)-based method to identify and quantitate long residence time enoyl-acyl carrier protein reductase inhibitors Anal Biochem, 474, 40-9. DOI: 10.1016/j.ab.2014.12.022 (34) Kharlamova, A. D.; Lushchekina, S. V.; Petrov, K. A.; Kots, E. D.; Nachon, F.; Villard-Wandhammer, M.; Zueva, I. V.; Krejci, E.; Reznik, V. S.; Zobov, V. V.; Nikolsky, E. E.; Masson, P. (2016) Slow-binding inhibition of acetylcholinesterase by an alkylammonium derivative of 6-methyluracil: mechanism and possible advantages for myasthenia gravis treatment Biochem J, 473, 1225-36. DOI: 10.1042/BCJ20160084 27 ACS Paragon Plus Environment


(35) Ticha, A.; Stanchev, S.; Vinothkumar, K. R.; Mikles, D. C.; Pachl, P.; Began, J.; Skerle, J.; Svehlova, K.; Nguyen, M. T. N.; Verhelst, S. H. L.; Johnson, D. C.; Bachovchin, D. A.; Lepsik, M.; Majer, P.; Strisovsky, K. (2017) General and Modular Strategy for Designing Potent, Selective, and Pharmacologically Compliant Inhibitors of Rhomboid Proteases Cell Chem Biol, 24, 1523-1536 e4. DOI: 10.1016/j.chembiol.2017.09.007 (36) Van Aller, G. S.; Nandigama, R.; Petit, C. M.; DeWolf, W. E., Jr.; Quinn, C. J.; Aubart, K. M.; Zalacain, M.; Christensen, S. B.; Copeland, R. A.; Lai, Z. (2005) Mechanism of time-dependent inhibition of polypeptide deformylase by actinonin Biochemistry, 44, 253-60. DOI: 10.1021/bi048632b (37) Hansen, R.; Peters, U.; Babbar, A.; Chen, Y.; Feng, J.; Janes, M. R.; Li, L. S.; Ren, P.; Liu, Y.; Zarrinkar, P. P. (2018) The reactivity-driven biochemical mechanism of covalent KRAS(G12C) inhibitors Nat Struct Mol Biol. DOI: 10.1038/s41594-018-0061-5 (38) Hulme, E. C.; Trevethick, M. A. (2010) Ligand binding assays at equilibrium: validation and interpretation Br J Pharmacol, 161, 1219-37. DOI: 10.1111/j.1476-5381.2009.00604.x (39) Bigott-Hennkens, H. M.; Dannoon, S.; Lewis, M. R.; Jurisson, S. S. (2008) In vitro receptor binding assays: general methods and considerations Q J Nucl Med Mol Imaging, 52, 245-53. DOI: (40) Guo, D.; Heitman, L. H.; AP, I. J. (2017) Kinetic Aspects of the Interaction between Ligand and G Protein-Coupled Receptor: The Case of the Adenosine Receptors Chem Rev, 117, 38-66. DOI: 10.1021/acs.chemrev.6b00025 (41) Gunther, J. R.; Du, Y.; Rhoden, E.; Lewis, I.; Revennaugh, B.; Moore, T. W.; Kim, S. H.; Dingledine, R.; Fu, H.; Katzenellenbogen, J. A. (2009) A set of time-resolved fluorescence resonance energy transfer assays for the discovery of inhibitors of estrogen receptor-coactivator binding J Biomol Screen, 14, 181-93. DOI: 10.1177/1087057108329349 (42) Vasta, J. D.; Corona, C. R.; Wilkinson, J.; Zimprich, C. A.; Hartnett, J. R.; Ingold, M. R.; Zimmerman, K.; Machleidt, T.; Kirkland, T. A.; Huwiler, K. G.; Ohana, R. F.; Slater, M.; Otto, P.; Cong, M.; Wells, C. I.; Berger, B.-T.; Hanke, T.; Glas, C.; Ding, K.; Drewry, D. H.; Huber, K. V. M.; Willson, T. M.; Knapp, S.; Müller, S.; Meisenheimer, P. L.; Fan, F.; Wood, K. V.; Robers, M. B. (2017) Quantitative, WideSpectrum Kinase Profiling in Live Cells for Assessing the Effect of Cellular ATP on Target Engagement Cell Chem Biol. DOI: 10.1016/j.chembiol.2017.10.010 (43) Andruska, N.; Mao, C.; Cherian, M.; Zhang, C.; Shapiro, D. J. (2012) Evaluation of a luciferase-based reporter assay as a screen for inhibitors of estrogen-ERalpha-induced proliferation of breast cancer cells J Biomol Screen, 17, 921-32. DOI: 10.1177/1087057112442960 (44) Renaud, J. P.; Chung, C. W.; Danielson, U. H.; Egner, U.; Hennig, M.; Hubbard, R. E.; Nar, H. (2016) Biophysics in drug discovery: impact, challenges and opportunities Nat Rev Drug Discov, 15, 679-98. DOI: 10.1038/nrd.2016.123 (45) Day, E. S.; Cote, S. M.; Whitty, A. (2012) Binding efficiency of protein-protein complexes Biochemistry, 51, 9124-36. DOI: 10.1021/bi301039t (46) Rich, R. L.; Myszka, D. G. (2010) Grading the commercial optical biosensor literature-Class of 2008: 'The Mighty Binders' J Mol Recognit, 23, 1-64. DOI: 10.1002/jmr.1004 (47) Schuck, P. (1997) Use of surface plasmon resonance to probe the equilibrium and dynamic aspects of interactions between biological macromolecules Annu Rev Biophys Biomol Struct, 26, 541-566. DOI: 10.1146/annurev.biophys.26.1.541 (48) Giannetti, A. M.; Koch, B. D.; Browner, M. F. (2008) Surface plasmon resonance based assay for the detection and characterization of promiscuous inhibitors J Med Chem, 51, 574-80. DOI: 10.1021/jm700952v (49) Leavitt, S.; Freire, E. (2001) Direct measurement of protein binding energetics by isothermal titration calorimetry Curr Opin Struct Biol, 11, 560-6. DOI: 10.1016/S0959-440X(00)00248-7


Page 28 of 30

Page 29 of 30 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


(50) Gruner, S.; Neeb, M.; Barandun, L. J.; Sielaff, F.; Hohn, C.; Kojima, S.; Steinmetzer, T.; Diederich, F.; Klebe, G. (2014) Impact of protein and ligand impurities on ITC-derived protein-ligand thermodynamics Biochim Biophys Acta, 1840, 2843-50. DOI: 10.1016/j.bbagen.2014.04.018 (51) Velazquez-Campoy, A.; Freire, E. (2006) Isothermal titration calorimetry to determine association constants for high-affinity ligands Nat Protoc, 1, 186-91. DOI: 10.1038/nprot.2006.28 (52) Di Trani, J. M.; De Cesco, S.; O'Leary, R.; Plescia, J.; do Nascimento, C. J.; Moitessier, N.; Mittermaier, A. K. (2018) Rapid measurement of inhibitor binding kinetics by isothermal titration calorimetry Nat Commun, 9, 893. DOI: 10.1038/s41467-018-03263-3 (53) Hansen, L. D.; Transtrum, M. K.; Quinn, C.; Demarse, N. (2016) Enzyme-catalyzed and binding reaction kinetics determined by titration calorimetry Biochim Biophys Acta, 1860, 957-966. DOI: 10.1016/j.bbagen.2015.12.018 (54) Bockman, M. R.; Engelhart, C. A.; Dawadi, S.; Larson, P.; Tiwari, D.; Ferguson, D. M.; Schnappinger, D.; Aldrich, C. C. (2018) Avoiding Antibiotic Inactivation in Mycobacterium tuberculosis by Rv3406 through Strategic Nucleoside Modification ACS Infect Dis, 4, 1102-1113. DOI: 10.1021/acsinfecdis.8b00038 (55) O'Shannessy, D. J.; Brigham-Burke, M.; Soneson, K. K.; Hensley, P.; Brooks, I. (1993) Determination of Rate and Equilibrium Binding Constants for Macromolecular Interactions Using Surface Plasmon Resonance: Use of Nonlinear Least Squares Analysis Methods Anal Biochem, 212, 457-468. DOI: 10.1006/abio.1993.1355 (56) Damian, L. In Protein-Ligand Interactions. Methods in Molecular Biology (Methods and Protocols); Williams, M., Daviter, T., Eds.; Humana Press: Totowa, NJ, 2013; Vol. 1008. (57) Li, H. J.; Lai, C. T.; Pan, P.; Yu, W.; Liu, N.; Bommineni, G. R.; Garcia-Diaz, M.; Simmerling, C.; Tonge, P. J. (2014) A structural and energetic model for the slow-onset inhibition of the Mycobacterium tuberculosis enoyl-ACP reductase InhA ACS Chem Biol, 9, 986-93. DOI: 10.1021/cb400896g (58) Feng, B. Y.; Shoichet, B. K. (2006) A detergent-based assay for the detection of promiscuous inhibitors Nat Protoc, 1, 550-3. DOI: 10.1038/nprot.2006.77



TOC Graphic


Page 30 of 30

Quantifying the Interactions between Biomolecules: Guidelines for

Recommend Documents