Quantitative Characterization of Three Carbonic Anhydrase Inhibitors by LESA Mass Spectrometry

Liquid extraction surface analysis (LESA) coupled to native mass spectrometry (MS) presents unique analytical opportunities due to its sensitivity, speed, and automation. Here, we examine whether this tool can be used to quantitatively probe protein–ligand interactions through calculation of equilibrium dissociation constants (Kd values). We performed native LESA MS analyses for a well-characterized system comprising bovine carbonic anhydrase II and the ligands chlorothiazide, dansylamide, and sulfanilamide, and compared the results with those obtained from direct infusion mass spectrometry and surface plasmon resonance measurements. Two LESA approaches were considered: In one approach, the protein and ligand were premixed in solution before being deposited and dried onto a solid substrate for LESA sampling, and in the second, the protein alone was dried onto the substrate and the ligand was included in the LESA sampling solvent. Good agreement was found between the Kd values derived from direct infusion MS and LESA MS when the protein and ligand were premixed; however, Kd values determined from LESA MS measurements where the ligand was in the sampling solvent were inconsistent. Our results suggest that LESA MS is a suitable tool for quantitative analysis of protein–ligand interactions when the dried sample comprises both protein and ligand.


■ INTRODUCTION
Protein−ligand binding is of key importance in all living organisms for the maintenance of protein structure and function, enzymatic activity, molecular interactions, signaling, and recognition. Therapeutic interventions exploit noncovalent interactions between the drug ligand and the target protein.
Understanding the binding characteristics, e.g., stoichiometry, equilibrium dissociation constants, binding sites, and nature of binding (covalent vs noncovalent interactions), is therefore of paramount importance in the discovery and development of drugs. Detailed characterization of these essential parameters can be performed by an array of analytical techniques including mass spectrometry (MS), 1 isothermal titration calorimetry (ITC), 2 surface plasmon resonance (SPR), 3 fluorescence spectroscopy, fluorescence correlation spectroscopy (FCS), 4 and nuclear magnetic resonance (NMR). 5,6 Native MS in conjunction with nanoelectrospray ionization (nESI), herein referred to as direct infusion for simplicity, has been widely exploited for the quantitative analysis of protein− ligand and protein−protein interactions, 7 including for the characterization of protein−small molecule, 8−10 protein− phosphopeptide, 11 protein−carbohydrate, 12,13 protein−protein, 14 protein−DNA, 15 and protein−RNA 16 complexes. Broadly, the capabilities of direct infusion MS for the study of protein−ligand binding are twofold. First, the resultant mass shift, as detected in the mass spectrum, is indicative of binding of a ligand to a protein. Second, the relative intensities of unbound protein ion peaks and ligand-bound protein ion peaks can be used to quantitatively monitor the extent of binding. Typically, the concentration of the protein of interest is fixed, and the ligand concentration is varied in order to perform a titration. With increasing ligand concentrations, a proportional increase in the intensity of the ligand-bound peak is observed. 17 The ratio of the protein−ligand complex and the protein, referred to as the response (R), is plotted as a function of initial ligand concentration. 15,18−20 This plot serves as the basis for determination of the equilibrium dissociation constant (K d ) via curve fitting with a linear 19 or nonlinear function. 12 The relative intensities or peak areas of the ligand-free and ligand-bound peaks are expressed for all observed charge states rather than for individual charge states in order to obtain a more comprehensive picture: individual charge states may disproportionally represent single, distinct conformational forms. 9,21 Protein−ligand titrations monitored in the gas phase assume that the ionization efficiency of unbound protein and the protein−ligand complex is similar in the case of a small ligand that does not substantially increase the mass and that the initial protein concentration does not change during the ionization and detection process. 17,22 Nevertheless, source conditions and instrumental parameters can affect the stability of complexes, necessitating cautious interpretation of the derived equilibrium dissociation constants. 11,12,14,19,23 Minimization of in-source dissociation of low-affinity, labile complexes is a critical step of gas-phase titrations and can be experimentally achieved through optimal selection of tuning conditions, reducing the diameter of the electrospray emitter, and/or application of a postacquisition correction to account for partial in-source dissociation. 24,25 The use of complementary in-solution techniques has proven to be helpful for further validation. 26,27 Here, we sought to evaluate the use of liquid extraction surface analysis (LESA) 28 MS for the measurement of binding affinities under native-like conditions. Native LESA MS appears to be an attractive strategy for the study of protein− ligand or protein−protein interactions as it affords automation prior to MS, allows direct sampling from a solid surface, and has the potential to lend itself to measurement of binding affinities in ex vivo tissue samples. Indeed, we have demonstrated analyses of endogenous protein assemblies and protein ligands from thin tissue sections by native LESA MS. 29−31 As an initial step in exploring the utility of LESA for quantitative analysis of protein−ligand binding, we examine the utility of native LESA MS for the measurement of K d values by use of a well-characterized model system comprising bovine carbonic anhydrase (CAH) and three ligands, chlorothiazide (CTZ), dansylamide (DNSA), and sulfanilamide (SLFA) ( Figure S1, Supporting Information). In numerous previous biophysical studies, CTZ,9,32 DNSA, 32,33 and SLFA 27,34,35 have shown to be potent inhibitors of various isoforms and derivatives of CAH. The carbonic anhydrases are a large group of isoenzymes, 34,36 some of which are membrane bound while others are cytosolic or secreted, that catalyze the interconversion of CO 2 and bicarbonate. 37 They represent key pharmaceutical targets owing to their involvement in infection, cancer, glaucoma, obesity, and hypertension. 36,38 Bovine CAH is a particularly well-described model protein, a 259-residue enzyme 39 that catalyzes the hydration of carbon dioxide and dehydration of bicarbonate. 37 In its hydrophobic, active center it contains a Zn 2+ ion coordinated by three imidazole moieties and a hydroxide group to which many of its inhibitors bind. 40−42 Sulfonamides and their derivatives, e.g., SLFA and CTZ, bind to CAH via coordination to the Zn 2+ in the center of the binding pocket. 36,40,43 The aromatic rings of Zn 2+ binders can occupy both the hydrophobic and the hydrophilic parts of the binding cavity, whereas the tails tend to orient toward the exit of the binding pocket. 36 DNSA favors the hydrophobic region of the binding pocket with its naphthalene ring rotated in a 54°angle. 44 We compared K d values measured by direct infusion MS with those measured by LESA MS. Two LESA strategies were trialled. In the first, the protein and ligand were mixed in solution prior to deposition and drying on a foil-covered glass slide. The resulting dried protein/ligand spots were sampled by LESA using a solvent comprising 25 mM ammonium acetate. We refer to this approach as "LESA premix ". In the second, solutions of the pure protein were deposited and dried onto foil-covered glass slides. The dried protein spot was sampled with a solvent comprising 25 mM ammonium acetate which contained the ligand. We refer to this approach as "LESA ligand ". To corroborate our data, K d values determined by MS were validated against SPR measurements. SPR is a label-free, surface-based, optical technique commonly deployed for the characterization of macromolecular interactions. 45,46 ■ EXPERIMENTAL SECTION Sample Preparation. Ten micromolar bovine carbonic anhydrase II (CAH) was mixed with various concentrations of CTZ, DNSA, and SLFA in 25 mM ammonium acetate (pH 7.0), incubated for 10 min at room temperature (21−22°C), and processed by direct infusion nESI MS or by native LESA MS. LESA was performed by use of the Triversa Nanomate platform (Advion, Ithaca, NY). Full details of the sampling procedures are described in the Supporting Information and summarized below.
LESA premix . A 1.5 μL amount of 10 μM CAH only or [CAH + ligand] was spotted onto a microscope glass slide covered with a layer of Al foil and allowed to air dry at room temperature. The dried spots were sampled with 25 mM ammonium acetate (pH 7.0) with a dwell time of 60 s following aspiration of 5.0 μL of solvent and deposition of 3.0 μL. A 3.5 μL amount was reaspirated to ensure complete sample collection.
LESA ligand . A 1.5 μL amount of 10 μM CAH was spotted onto a microscope glass slide covered with Al foil and allowed to air dry at room temperature. The dried spots were sampled with 25 mM ammonium acetate (pH 7.0) containing various concentrations of CTZ, DNSA, or SLFA at room temperature. The sampling protocol was comprised of 10 mix and repeat cycles at the settings above to control for the 10 min incubation used for LESA premix .
Mass Spectrometry and Data Analyses. All protein− ligand titration experiments were performed on an Orbitrap Elite mass spectrometer (Thermo Fisher Scientific, Bremen, Germany) with an automatic gain control target of 1 × 10 6 charges in full-scan mode, 250°C source temperature, 25− 30% RF voltage, 5 or 10 ms injection time, 1000 microscans, and 120 000 resolution. Data were acquired in triplicate in the 200−4000 m/z mass range for 3−5 min. The ion intensities of each acquisition were summed in Xcalibur (Foundation 3.0 SP2), written to a raw file, and further processed in MATLAB (version R2017a and statistics toolbox, Mathworks, Inc., Natick, MA) and Python (version 3.7.6) using in-house software (see Supplemental File 1, Supporting Information). Using MATLAB, the peak areas of protonated unbound and ligand-bound protein peaks in all observed charge states (9+, 10+, and 11+) were determined using the trapezium rule. The average ratio, R, of bound to unbound protein across the three charge states was determined. No baseline correction was performed, and only the main protonated ion was considered. Titration curves were then generated in Python. The mean fraction bound (R/(R + 1)) from experimental replicates was plotted as a function of the initial ligand concentration L 0 (μM).
Equilibrium dissociation constants were then calculated by fitting the titration curves using a nonlinear sum of least squares. In the case of direct infusion measurements, equilibrium dissociation constants were calculated by fitting the titration curves to eq 1 (described in detail in ref 25), where P 0 denotes the initial protein concentration, K d is the equilibrium dissociation constant, and f sat is a correction factor to account for gas-phase dissociation and is the experimental bound fraction at saturation.
For the LESA premix and LESA ligand measurements, corrected values of L 0 and/or P 0 were used to account for extraction efficiency and dilution effects associated with the sampling process as described below. Determination of P 0 Correction for Combined Protein Extraction Efficiency and Dilution Effects. LESA premix . A paper grid print of a 1584-well plate was placed into the plate holder of the stage directly under the glass slide covered with Al foil containing the samples in order to guide manual selection of xy coordinates. In addition, a halved 96well plate was placed beside the glass slide for collection of extracted aliquots. The advanced user interface (AUI) was used to aspirate 5 μL of 25 mM ammonium acetate and dispense that onto a 1.5 μL dried droplet of 10 μM CAH at a 3 μL/min flow rate, 1.4−1.6 mm height, mimicking the formation of a liquid junction at a 60 s dwell time. A 3.5 μL volume was reaspirated for collection into a half 96-well plate. For each location sampled (n = 15), two 1.2 μL samples were analyzed on the DeNovix DS-II spectrophotometer. The UV spectra at 280 nm (A 280 ) were used to calculate the corrected concentration (in micromolar) based on the extinction coefficient of CAH (50 070 M −1 cm −1 47 ). The absorption of 25 mM ammonium acetate was subtracted from that of CAH. Data were exported as .csv files, and the correction factor incorporating both the dilution effect and the extraction efficiency of LESA was calculated in Prism GraphPad 6.01. The mean concentration of 15 sample spots (micromolar) was used as the corrected value for P 0 .
LESA ligand . The corrected value for P 0 was determined as above except that the sampling procedure was adjusted to mimic the LESA ligand sampling procedure in which 10 mix and repeat steps were included. To aid sample collection, the blowout tab was checked in the software and the volume was released at a pressure of 3 psi and duration of 10 s. In total, 10 spots were sampled and analyzed as above. The mean Figure 1. Native mass spectra of CAH following direct infusion (a), LESA premix (b), and LESA ligand (c) sampling. (Inset) CSD of CAH in its Zn 2+bound form with the 10 + ion being the most abundant. Native mass spectra of CAH in the presence of 0.1 μM CTZ following direct infusion (d), LESA premix (e), and LESA ligand (f). Mass spectra of CAH in the presence of 1.5 μM DNSA obtained in direct infusion mode (g) and subsequently LESA premix (h) and LESA ligand (i) sampling. Native mass spectra of CAH in the presence of 5 μM SLFA following direct infusion (j), LESA premix (k), and LESA ligand (l) sampling. Dots represent the protein−ligand complexes.
Journal of the American Society for Mass Spectrometry pubs.acs.org/jasms Research Article concentration (micromolar) was used directly as the corrected value for P 0 . LESA Correction Factor for the Total Free Ligand Extraction Efficiency and Dilution Effects. The absorption spectra of different concentrations of CTZ, SLFA, and DNSA were obtained by spectrophotometry (wavelength 220−500 nm) for a solution reference curve without any sample manipulation, n = 3. The wavelength of maximum absorption was determined from these samples. A 1.5 μL amount of ligand was deposited on the surface of a glass slide covered with Al foil, air dried, placed on top of a 1584-well paper grid at different concentrations in order to determine its coordinates, and extracted in 25 mM ammonium acetate in triplicate using the aforementioned LESA AUI sampling conditions. The absorption spectra of the samples following LESA AUI were collected and exported as .csv files, after confirming no significant red or blue shift, in order to calculate the correction factor accounting for both dilution effects and extraction efficiency and any evaporation occurring during extended LESA sampling in the case of the LESA ligand for the free ligands at different concentrations. The mean and standard deviation values were calculated in Prism GraphPad 6.01 based on the solution reference values and the AUI collected samples. The mean value ± STD of the entire concentration range measured yielded the correction factor used in the K d calculations. The mean correction factor was used to multiply the initial reference L 0 values to obtain its paired corrected L 0 .
Traveling Wave Ion Mobility Spectrometry MS, Surface Plasmon Resonance, and Circular Dichroism Spectroscopy. Details are given in the Supporting Information.
■ RESULTS AND DISCUSSION Figure 1 shows the mass spectra obtained following direct infusion, LESA premix , and LESA ligand for CAH and the three ligands at various ligand concentrations (0, 0.1, 1.5, and 5 μM). For CAH, a typical narrow charge state distribution (CSD) was observed ranging from 9 + to 11 + with the 10 + charge state being the most dominant (Figure 1a−c), indicative of retention of its native structure. The deconvoluted mass, 29 088.86 Da, revealed that CAH was N-terminally acetylated with a Zn 2+ ion bound to the protein ( Figure S2, Supporting Information) and is in good agreement (<5 ppm) with the theoretically calculated mass of the protein with 29 088.78 Da (accession ID P00921).
In the presence of CTZ, a 295 Da mass shift was observed in the direct infusion ( Figure 1d) and LESA premix (Figure 1e) samples at concentrations of 0.1 μM and above. No obvious ligand binding was observed at 0.1 μM CTZ with LESA ligand MS (Figure 1f) but was observed at concentrations above 1 μM ( Figure S3, Supporting Information). For DNSA, ligand binding was observed in the mass spectra (Δ 250 Da) at 1.5 μM and above in both direct infusion ( Figure 1g) and LESA premix (Figure 1h). No ligand binding was observed following LESA ligand (Figure 1i) at this concentration; however, it could be observed at concentrations above 3 μM ( Figure S4, Supporting Information). In the presence of SLFA, the mass spectra revealed that a 172 Da mass shift, corresponding to ligand binding, occurred from a 5 μM concentration onward in direct infusion mode (Figure 1j) and in LESA premix mode (Figure 1k), but a well-resolved peak for the protein−ligand complex could not be detected in LESA ligand (Figure 1l) mode. For LESA ligand mode, ligand binding was observed at concentrations of 10 μM and above ( Figure S5, Supporting Information). There was no evidence for binding of multiple ligands at any concentration. No binding was observed for any of the ligands in the presence of 1% (v/v) formic acid ( Figures  S3−S5, Supporting Information).
The discrepancies in the concentrations at which the protein−ligand complexes were observed in the various sampling modes (direct infusion, LESA premix , and LESA ligand ) together with the deviations in the observed relative abundance suggest that the extraction efficiency of LESA and the dilution inherent in the LESA sampling process result in changes to P 0 and/or L 0 , thus shifting the equilibria described in eq 1. In order to calculate the K d values from LESA data, it is necessary to determine the P 0 and L 0 values corrected for these factors. Literature reference values are displayed with an indication of CAH isoform, temperature, and analytical technique if provided. K d values determined from LESA experiments are shown before and after corrections of L 0 and/or P 0 . a P 0 and L 0 corrected for LESA extraction efficiency and dilution. b P 0 corrected for extraction efficiency and dilution. c Uncorrected for LESA extraction efficiency and dilution.  Figure S6, Supporting Information, along with the derived binding parameters; see Table S1, Supporting Information.) The K d values calculated for LESA premix and LESA ligand MS are given with (i) P 0 and L 0 corrected for the LESA extraction efficiency and dilution, (ii) only P 0 corrected for the LESA extraction efficiency and dilution, and (iii) uncorrected P 0 and L 0 .
For CAH and CTZ, there is good agreement between the K d values calculated from direct infusion MS, P 0 -and L 0 -corrected LESA premix MS, and SPR measurements. The measured values fall within the range described in the literature. For CAH and DNSA, there is good agreement between the K d values measured by direct infusion and P 0 -and L 0 -corrected LESA premix MS. For CAH and SLFA, both the uncorrected and the corrected K d values calculated from LESA premix MS are in reasonable agreement with those calculated by direct infusion MS. Differences between the reported K d values from the SPR and MS measurements have been noted previously and the limitations discussed. 27 For the K d values calculated from LESA ligand MS measurements, the best agreement with the K d calculated from direct infusion MS for CTZ was achieved with the uncorrected value. For DNSA, neither corrected nor uncorrected LESA ligand values were in agreement with that obtained following direct infusion MS, although all were closer to those obtained by SPR. For SLFA, corrected values were in the best agreement with that obtained by direct infusion.
Calculation of K d Values Following Direct Infusion and LESA premix MS. Conventionally, the K d values are calculated from the direct infusion MS measurements after careful tuning, assuming that no major ligand dissociations occur in source or during detection and that the stoichiometry of a given protein−ligand (PL) complex formed in solution is maintained in the gas phase; 12 however, improved accuracy in the K d values has been demonstrated when in-source dissociation is taken into account. 24,25 During LESA, the airdried protein and/or protein−ligand spots are resolubilized and extracted from the substrate for MS. The extraction efficiency is below 100% if not all of the protein or ligand deposited is successfully extracted into the LESA droplet. Moreover, in the extraction process, the dried droplet is extracted into a larger volume (1.5 μL deposited sample volume extracted in 5 μL, i.e., ∼3.33-fold dilution). To address these factors, the extraction efficiencies of the protein and free ligand were determined, enabling us to calculate the corrected P 0 and L 0 values which were then used to determine the K d values using eq 1.
First, we assumed that the extraction efficiencies of the free, unbound protein and the ligand-bound protein were the same. To account for the combined dilution and extraction efficiency effects for the protein, we deposited and extracted CAH mimicking the LESA premix conditions ( Figure S7, Supporting Information) but collected the extracted sample rather than introducing it to the mass spectrometer. The protein concentration was determined by absorbance at 280 nm. An example of a representative absorption spectrum is shown in Figure S8, Supporting Information, for 10 μM CAH and 10 μM CAH following LESA premix sampling. The mean concen-tration (n = 15) after LESA premix sampling of 10 μM CAH protein spots was 1.51 ± 0.59 μM (Tables S2a and S2b, Supporting Information), indicating that the extraction efficiency was approximately 50%. These data also demonstrate the variation in protein sampling by LESA which can arise due to variable spreading of the sampling droplet, and therefore area sampled, and variation in the proportion of sampling solvent recovered on reaspiration. 52 The corrected protein concentration was used in subsequent K d value calculations for the LESA premix conditions.
To estimate a correction factor for the combined effects of the LESA premix -associated dilution and extraction efficiency for each ligand, we conducted LESA extraction, sample collection, and UV−vis spectrophotometry of the three compounds at various concentrations. A range of concentrations for each ligand without LESA sampling served as reference values. For all three ligands, the wavelength with maximum absorbance (λ max ) was determined experimentally (see Figure S9, Supporting Information, λ max = 277 (CTZ), 244 (DNSA), and 258 nm (SLFA)) The absorbance at λ max was plotted versus the ligand concentration, where it was evident that the mean absorbance values of the reference data points were higher than those following LESA premix sampling. A correction factor was calculated for each ligand across all points based on the ratio of the absorbance of the solution reference sample and the equivalent LESA sample. The mean correction factor was determined to be 0.32 ± 0.03 for CTZ, 0.33 ± 0.05 for DNSA, and 0.39 ± 0.08 for SLFA (Tables S3−S5, Supporting  Information), indicating the extraction efficiencies were ∼100% for all concentrations after accounting for the dilution upon LESA premix sampling. These mean correction factors ± 1 STD were used to determine adjusted L 0 values for the LESA premix conditions. The data used for K d value determination by direct infusion and LESA premix are shown in Figure S10, Supporting Information. The results are summarized in Table 1. For CTZ and SLFA, the corrected and uncorrected K d values do not differ significantly. For DNSA, however, the uncorrected K d value is over 2-fold greater than that for direct infusion but correction of P 0 and L 0 results in good agreement between the K d values determined by direct infusion MS and LESA premix MS.
Calculation of K d Values Following LESA ligand MS. Our initial assumption for the LESA ligand mode was that only P 0 would require correction for LESA extraction efficiency and dilution, as the ligand is in the LESA solvent. We determined the corrected P 0 concentration by mimicking the LESA ligand sampling and collecting the sample for UV−vis spectrophotometry. A representative UV−vis spectrum is shown in Figure  S11, Supporting Information. The mean extracted protein concentration (n = 10) was calculated to be 2.60 ± 0.62 μM (Table S6, Supporting Information). This value was used in subsequent K d value calculations as the corrected P 0 . Notably, this value is ∼2-fold higher than the value obtained for LESA premix , likely due to the repeated mix steps in this LESA mode. (In addition, we could not rule out some solvent evaporation during the extended length of sampling, although the sample stage was kept covered at a constant 21°C.) Individual titration plots used for determination of K d values, with corrected P 0 , can be found in Figure S12, Supporting Information.
The general lack of agreement between K d values determined by direct infusion MS and LESA ligand MS despite Journal of the American Society for Mass Spectrometry pubs.acs.org/jasms Research Article correction of P 0 suggested that our assumption about L 0 was incorrect. We next considered whether the apparent L 0 after LESA ligand may increase as a consequence of evaporation during the 10 min sampling process. The sampling process, in which ligand solutions comprising a range of concentrations were deposited, aspirated, mixed, and reaspirated from the foilcovered slide, was recreated, and the resulting samples were collected. The absorbance spectra of these samples were compared with reference solutions that had not been subjected to any manipulation ( Figure S13, Supporting Information). For DNSA and SLFA, the pairwise comparison of absorbances at the maximum wavelength (Figures S13f and S13i, Supporting Information) revealed that at higher concentrations L 0 was higher following LESA ligand sampling than that for the reference solution. Nevertheless, the mean correction factor across the concentration range was ∼1 in both cases (1.063 and 0.9904 for DNSA and SLFA, respectively, see Tables S7 and S8, Supporting Information). For CTZ, the pairwise comparison revealed higher absorbances for the LESA ligand samples at all concentrations ( Figure S13c, Supporting Information). Moreover, a blue shift in the absorption spectrum was observed. The calculated mean correction factor was 1.5-fold higher than the initial ligand concentrations (Table S9, Supporting Information). The K d values were calculated using corrected values for both P 0 and L 0 (see Table  1 and Figure S14, Supporting Information); however, the fits and agreement between values from direct infusion and LESA premix did not improve. For completeness, we also calculated the K d values with no correction for P 0 and L 0 ( Figure S14, Supporting Information). For CTZ, the uncorrected K d value is in best agreement with the value obtained by direct infusion. For DNSA, correction of L 0 and/ or P 0 did not significantly change the calculated K d values. For SLFA, the corrected values were in better agreement than the uncorrected values. It is worth noting that the binding kinetics may also result in differences in the K d values determined via LESA premix and LESA ligand . For the protein−ligand interactions studied here, the calculated time to 99% occupancy was less than 10 min in all cases. For protein−ligand interactions with slower kinetics, it may be necessary to increase the incubation times (LESA premix ) or number of mix and repeat samples (LESA ligand ). Structural Insights. We also considered whether the discrepancies seen for LESA ligand sampling might be due to conformational changes induced during the drying of the protein spots, that is, if structural changes and additional conformers are produced on drying which are maintained following LESA extraction then there is the possibility that a population of protein exists that is unable to bind ligand. To test that, we coupled native LESA MS with traveling wave ion mobility spectrometry for a single ligand concentration (30 μM) for each compound in the presence of controls. All mass spectra ( Figure S15, Supporting Information) displayed a narrow charge state distribution from 8 + to 11 + . The 8 + ions showed multiple conformers (Figure S16, Supporting Information); however, the 8 + charge state was not observed in high-resolution native MS and was not part of the titration plots described above. For the 9 + charge state, two conformations were observed for CAH following LESA whereas a single peak was observed following direct infusion, as observed previously. 53 It is also possible that the 9 + ion shows some heterogeneity in the gas phase, supported by the observation of two conformers with different CCSs reported by Harrison et al. 54 The drift time profiles for the ligand-bound CAH suggest binding to both the major and the minor conformers. The 10 + and 11 + ions were characterized by one main, dominant conformation in both direct infusion and LESA modes. Table S10, Supporting Information, summarizes the estimated mean TW CCSs N2→He of the ligand-bound and free charge states and the CCSs of CAH only. Calculation of the TW CCSs N2→He in the presence of the three CAH inhibitors indicates that the coordination of these ligands in the active center of CAH do not result in a significant increase in TW CCSs N2→He . Experimental CCS values yielded similar values for direct infusion, LESA premix , and LESA ligand samples.
Finally, we assessed the near-UV secondary structure of CAH in the presence and absence of CTZ, DNSA, and SLFA in a range of concentrations. The circular dichroism spectra did not reveal any major conformational rearrangements in the secondary structure of the protein upon interactions with these ligands ( Figure S17, Supporting Information).

■ CONCLUSIONS
The results show that native LESA MS is a potentially useful tool for quantitative determination of protein−ligand interactions using CAH and three inhibitors covering submicromolar to low micromolar K d values as examples of monovalent interactions. To estimate the broader utility of LESA-MS, it would be particularly useful to examine a variety of larger proteins including those with a number of subunits capable of binding multiple ligands, proteins with multivalent binding sites, and proteins that bind numerous ligands with positive or negative cooperativity. When using LESA it is necessary to determine corrections for P 0 and L 0 to better describe the binding equilibrium following sampling from a substrate taking into account extraction efficiency and dilution effects. We found a good match between the K d values derived from data obtained following direct infusion MS and LESA premix MS for all three ligands. In contrast, the K d values from the LESA ligand sampling were inconsistent, and this could not be addressed by the corrections described here. We found that although LESA ligand MS is not suitable for accurate K d value calculations, it can still confirm whether a ligand can bind or not, and CCS measurements can be observed with comparable values to those from direct infusion MS and LESA premix MS. Our findings suggest that native LESA MS has the potential to be developed as an automated, high-throughput platform to support drug development.
Extended methods, software for calculation of AUCs and fitting of titration curves, chemical structures of ligands, mass spectra, SPR sensograms and calculations, photographs of the advanced user interface, UV absorption spectra, calculation of corrected P 0 and L 0 values, titration plots, TWIMS mass spectra and drift time profiles, calculated collision cross sections, CD spectra (PDF)