|Home | About | Journals | Submit | Contact Us | Français|
Cyclophilin A (CypA) is an overexpressed protein in lung cancer tumors and as a result is a potential therapeutic and diagnostic target. Here we utilize an H/D exchange- and MALDI mass spectrometry-based assay, termed single-point SUPREX (Stability of Unpurified Proteins from Rates of H/D Exchange), to screen two chemical libraries, including the 1280-compound LOPAC library and the 9600 compound DIVERSet library, for binding to CypA. This work represents the first application of single-point SUPREX using a pooled ligand approach, which we demonstrate is capable of screening rates as fast as six seconds/ligand. The false positive and false negative rates determined in the current work using a set of control samples were 0% and 9%, respectively. A false positive rate of 20% was found in screening the actual libraries. Eight novel ligands to CypA were discovered including: 2-(α-naphthoyl)ethyltrimethyl-ammonium iodide, (E)-3-(4-t-Butylphenylsulfonyl)-2-propenenitrile, 3-(N-benzyl-N-isopropyl)amino-1-(naphthalen-2-yl)propan-1-one, cis-diammineplatinum (II) chloride, 1-(3,5-dichlorophenyl)-1H-pyrrole-2,5-dione, N-(3-chloro-1,4-dioxo-1,4-dihydro-2-naphthalenyl)-N-cyclohexylacetamide, 1-[2-(3,4-dimethoxyphenyl)ethyl]-1H-pyrrole-2,5-dione, and 4-(2-methoxy-4-nitrophenyl)-1-methyl-10-oxa-4-azatricyclo[18.104.22.168~2,6~]dec-8-ene-3,5-dione. These compounds, which had moderate binding affinities to CypA (i.e., Kd values in the low micromolar range), provide new molecular scaffolds that might be useful in the development of CypA targeted diagnostic imaging or therapeutic agents for lung cancer.
The ready availability of natural and synthetic compound libraries of peptides, oligonucleotides, and small molecules (MW <1000 amu) has created a new route to the discovery of novel ligands that bind to specific protein targets. This route, which typically involves the use of high-throughput screening (HTS) strategies to rapidly screen library members for protein-ligand binding, can be useful for the identification of new protein-targeted diagnostic and/or therapeutic agents. For example, the protein target in this work, CypA, is a peptidyl-prolyl cis-trans isomerase that has been shown to be overexpressed in lung cancer tumor cells,1, 2 and as such is an attractive target for potential diagnostic imaging probes. It also appears to play an important role in the lung tumor growth3 and thus could be an important therapeutic target as well.
A goal of this work was to identify novel CypA ligands that might be useful in the development of new CypA-targeted molecular imaging agents for the diagnosis of lung cancer. To this end, small molecule compounds in two chemical libraries were screened for CypA binding. One library was the Library of Pharmacologically Active Compounds, or LOPAC library, which consists of 1280 pharmacologically active compounds that span a range of cell signaling and neuroscience areas and that include a number of already marketed drugs. The other library was the 9600-member DIVERSet library, which consists of a highly diverse collection of drug-like small molecules selected through the use of a 3D pharmacophore analysis method to ensure broad coverage among the various types of interactions between target and ligand.
The high-throughput screening assay used in this work relies on an H/D exchange- and mass spectrometry-based methodology, termed single-point SUPREX, to detect the increase in a protein’s thermodynamic stability upon ligand binding. We initially developed the single-point SUPREX assay in the context of a model protein-peptide system4, and the assay was used to screen the 880-compounds in the Prestwick Chemical Library for binding to CypA.5 While no new CypA ligands were identified in the Prestwick Chemical Library screening, that work did establish a 3-minute/ligand throughput of the single-point SUPREX protocol. As part of the current work, we describe a compound-pooling strategy that increases the screening rate of the single-point SUPREX protocol to less than 10 seconds/ligand. Our results demonstrate the feasibility of the new methodology for screening 1000–10,000 member libraries, and identify eight novel ligands to CypA.
The human CypA used in this work was obtained from T7 Express E. coli transformed with pGEX-CypA, a plasmid encoding the glutathione S-transferase (GST) fusion protein of CypA. Briefly, the E. coli was grown up in LB-Amp media (sodium chloride, tryptone, yeast extract, and ampicillin) at 37°C. Isopropyl β-D-1-thiogalactopyranoside (Mediatech, Inc., Mannassas, VA) was added to a final concentration of 0.125 mM to induce overexpression when the OD600 of the culture reached 0.6–0.8 units. After 3.5 hours, the E. coli was pelleted by centrifuging at 4000 × g for 10 min. Cells were lysed with the addition of BugBuster reagent (EMD Chemicals, Inc., Gibbstown, NJ) and the insoluble material was pelleted at 17,000 × g for 20 min. CypA was purified from the supernatant using GST•Bind Resin (EMD Chemicals, Inc.), and it was concentrated with a Vivaspin 20 centrifugal concentrator (Sartorius, Bohemia, NY). Thrombin (EMD Chemicals, Inc.) was added to remove the GST tag. The GST tag was subsequently removed via the GST•Bind Resin, and the thrombin was removed with a HiTrap Benzamidine FF column (GE Healthcare Bio-Sciences Corp., Piscataway, NJ). The CypA was further concentrated to about 100 µM. All CypA concentrations were determined using absorbance measurements at 280 nm and an extinction coefficient of 8730 M−1 cm−1 as calculated through ProtParam using the protein’s amino acid sequence (http://ca.expasy.org/tools/protparam.html).
The LOPAC library (Sigma-Aldrich, St. Louis, MO) was obtained from the Small Molecule Synthesis and Screening Facility in the Center for Chemical Biology at Duke University. The library compounds were provided as 10 mM solutions in dimethyl sulfoxide (DMSO) and were pipetted into 128 10-member pools for a resulting concentration of 1 mM of each compound. The 9600-compound DIVERSet library was from ChemBridge Corp. (San Diego, CA) and was provided as 10 mM solutions in DMSO. The DIVERSet library compounds were combined into 960 10-member pools for a resulting concentration of 1 mM of each compound. 3-[5-(4-Fluorobenzylidene)-4-oxo-2-thioxo-1,3-thiazolidin-3-yl]propanoic acid, N-(3-chloro-1,4-dioxo-1,4-dihydro-2-naphthalenyl)-N-cyclohexylacetamide, 1-(3,5-dichlorophenyl)-1H-pyrrole-2,5-dione, 4-(2-methoxy-4-nitrophenyl)-1-methyl-10-oxa-4-azatricyclo[22.214.171.124~2,6~]dec-8-ene-3,5-dione, and 1-[2-(3,4-dimethoxyphenyl)ethyl]-1H-pyrrole-2,5-dione were from ChemBridge Corp. (San Diego, CA). Cyclosporin A (CsA) was purchased from LKT Laboratories (St. Paul, MN). Deuterium oxide (D2O; 99.9% atom D), sodium deuteroxide (40 wt.% in D2O, 99.9% atom D), sinapic acid (SA; 98% titration), and cis-diammineplatinum (II) chloride (cisplatin) were from Sigma-Aldrich (St. Louis, MO). (E)-3-(4-t-butylphenylsulfonyl)-2-propenenitrile (Bay 11–7085) and 2-(α-naphthoyl)ethyltrimethylammonium iodide (α-NETA) were purchased from Alexis Biochemicals (Lausen, Switzerland). 3-(N-Benzyl-N-isopropyl)amino-1-(naphthalen-2-yl)propan-1-one hydrochloride (ZM 39923) was from Tocris Bioscience (Ellisville, MO). DMSO was from Acros Organics (Geel, Belgium), and phosphoric acid-d3 was from Cambridge Isotope Laboratories, Inc. (Andover, MA). Guanidinium chloride (GdmCl) was purchased from EMD Chemicals, Inc. (Gibbstown, NJ), and trifluoroacetic acid (TFA) was from Halocarbon (River Edge, NJ). Acetonitrile (ACN) was purchased from Fisher (Fair Lawn, NJ). The chymotrypsin, N-succinyl-Ala-Ala-Pro-Phe-p-nitroanilide, and trifluoroethanol were all from Sigma-Aldrich (St. Louis, MO).
A deuterated solution of 20 mM sodium phosphate buffer, pD 7.4, containing 1.5 M GdmCl was used for the single point SUPREX experiments. A series of solutions containing 20 mM phosphate buffer, pD 7.4, containing GdmCl concentrations ranging from 0 to 4 M was prepared for the full SUPREX analyses performed on CypA in the presence of the hit compounds. The GdmCl used for buffer preparation was previously deuterated by repeated dissolution and lyophilization (4 times total) of 19 g GdmCl in 25 mL D2O. The pD of the buffers were adjusted by adding sodium deuteroxide and monitoring the pH using a Jenco microcomputer pH-Vision 6072 pH meter equipped with a Futura calomel pH electrode from Beckman (Fullerton, CA). The pD of the buffers was obtained by adding 0.4 to the pH reading to correct for isotope effects.6 The GdmCl concentration in the buffers was measured using a refractometer (Bausch & Lomb, Rochester, NY) utilizing the method reported by Nozaki.7
MALDI mass spectra were acquired on an Ultraflex II TOF/TOF from Bruker Daltonics (Billerica, MA). Spectra were the sum of 100 laser shots collected in the linear mode using a 100 Hz Nd:YAG laser. Positive ion mass spectra were manually collected using the following parameters: an ion source 1 voltage of 25 kV, an ion source 2 voltage of 23.4 kV, a lens voltage of 6.5 kV, and a delay time of 100 ns. The MALDI matrix was SA, and it was prepared as a saturated solution containing 25% water and 75% ACN with 0.1% TFA.
For the screening assay, 1 µL of each 10-member pool containing 1 mM of each library compound in DMSO was combined with 9 µL of the deuterated exchange buffer containing 1.5 M GdmCl. A 1 µL aliquot of protonated CypA solution (i.e., CypA dissolved in H2O-containing buffer), which was ~100 µM in CypA, was added to the resulting 10-µL volume of the D2O-containing compound/exchange buffer, and the H/D exchange reaction proceeded for 35 minutes before it was quenched by transferring 1 µL of the exchange reaction to 9 µL of ice-cold matrix solution. A 1–2 µL aliquot of the quenched reaction mixture was spotted onto a MALDI sample target and dried with N2. Five replicate MALDI mass spectra were collected to determine an average molecular mass of the deuterated CypA for each sample, and this was used to calculate an average ΔMass value. Average ΔMass values were obtained by subtracting the expected average molecular mass of the fully protonated CypA from the average molecular mass measured for each sample. Samples were prepared and analyzed in sets of twelve using a 12-channel pipette. Each set of twelve samples included ten pools of 10 compounds, one positive control containing 1 mM CsA (Kd ~30–200 nM8–12), and one negative control containing only DMSO. The H/D exchange reactions were staggered at 10-min intervals, so while one sample set was undergoing H/D exchange, MALDI mass spectra were being collected on the previous sample set.
The average molecular weight of CypA was determined from the MALDI-TOF mass spectra using either an Excel macro or MATLAB (The MathWorks, Inc., Natick, MA) script, both of which were developed in-house. Each takes the spectral data in text file form containing m/z and intensity values for each data point and performs a 19-point floating average smoothing of the data, a calibration using the internal standard, and a center of mass calculation for the CypA peak.
In the LOPAC library screening, the ΔMass values of the control data were averaged, and the cut-off value was determined by subtracting 2.5 standard deviations from the average negative control ΔMass value. In the DIVERSet library screening, the ΔMass values of the control data from the screen were smoothed using a 9-point central moving average. An average ΔMass value (i.e., a ΔMassav value) for each negative control was calculated using the 9-point moving average and was subsequently used to calculate a cut-off value for each set of ten pools. The average of all the 9-point standard deviations was calculated, multiplied by 2.5, and then subtracted from the ΔMassav value to determine the cut-off values. Samples with ΔMass value less than or equal to the cut-off value were classified as preliminary hits (i.e., one of the compounds in the pool of 10 binds to CypA).
In Eq. (1), σc+ is the standard deviation of the positive control, σc− is the standard deviation of the negative control, µc+ is the mean of the positive control, and µc− is the mean of the negative control.
The preliminary pools classified as hits were subjected to a “second-tier” analysis, in which the initial hit pools were re-screened using the single-point SUPREX protocol described above. We have previously shown that such a two tier screening approach reduces the occurrence of false positives in single-point SUPREX.5 Each hit-pool was re-screened under the same conditions as the initial screen. The individual compounds from the pools that were confirmed as hit-pools in the second-tier screen were individually screened using the single-point SUPREX protocol, just as described above for the compound-pools, except that the ligand stock solution only contained a single ligand at 1 mM.
Full SUPREX analyses were performed on CypA in the absence and in the presence of the individual compounds selected as hits. In these full SUPREX analyses 1 µl of a 1 mM DMSO solution of the hit ligand (or just DMSO in the case of the CypA anlayses performed in the absence of ligand) was added to 9 µl of the SUPREX buffers, which in the full SUPREX analyses contained [GdmCl] concentrations ranging from 0 to 4 M (see above). A 1 µL aliquot of protonated CypA solution, which was ~100 µM in CypA, was added to each 10-µL volume of compound/exchange buffer, and the H/D exchange reaction proceeded for 35 minutes before it was quenched by transferring 1 µL of the exchange reaction to 9 µL of ice-cold matrix solution. A 1–2 µL aliquot of the quenched reaction mixture onto a MALDI sample target and dried with N2. Ten replicate MALDI mass spectra were collected to determine an average molecular mass of the deuterated CypA for each sample, and this was used to calculate an average ΔMass value at each denaturant concentration for each ligand.
Plots of ΔMass vs. [GdmCl] were constructed for the CypA in the presence and absence of the hit compounds. The data in the plots were fit to the following four-parameter sigmoidal equation using a nonlinear regression in SigmaPlot (Systat Software Inc., San Jose, CA).
In Eq. (2), ΔM0 is the pre-transition baseline (i.e., the change in mass of CypA before global exchange), a is the amplitude of the curve in Da, [GdmCl] is the molar concentration of GdmCl, is the [GdmCl] at the transition midpoint of the curve, and b is a parameter describing the steepness of the transition. These parameters were typically allowed to float when fitting the data to Eq. (2).
The values obtained from Eq. (2) were used in Eq. (3) to determine the magnitude of the increase in thermodynamic stability (i.e., decrease in folding free energy), assuming the m-value, which is a measure of a protein’s folding cooperability, is constant in the presence and absence of ligand.
In Eq. (3), ΔΔGf is the change in the protein’s folding free energy (i.e., the protein’s free energy of binding), m is the m-value, and is the change in the value of the CypA upon binding to a ligand. Eq. (3) was derived from the general form of the SUPREX equation found in reference 14. The m-value was assigned a value of 3.7 kcal mol−1 M−1 based on an average of the values reported for CypA in reference 12.
The following equation was used to calculate dissociation constants, Kd, for the protein-ligand complexes.15
In Eq. (4), [L] is the molar concentration of the free compound, n is the number of independent binding sites (equal to 1), R is the universal gas constant, and T is the temperature in K. Because the ligand concentration was present at 10-fold excess over the CypA, the total ligand concentration was used for the free ligand concentration.
The theoretical SURPEX curves generated in this work (Figure 1) were generated exactly as described in reference 5, with the exception that the pre- and post-transition baselines were set to 45 and 65 Da (respectively) in order to more accurately represent the pre- and post-transition baselines observed in the current work.
The ability of the hit compounds to inhibit the peptidyl prolyl cis-trans isomerase activity of CypA was determined using a chymotrypsin-coupled enzyme assay that has been previously described.16, 17 In these assays, 875 µL of ice-cold buffer (50 mM HEPES, 100 mM NaCl, 1% DMSO, pH 8.0) containing 5 nM CypA and 10 µM of the hit-compound was combined with 100 µL of an ice-cold chymotrypsin solution (30 mg/mL in 0.001 HCl) in a quartz cuvette. A 25-µL aliquot of a stock solution of the peptide substrate N-succinyl-Ala-Ala-Pro-Phe-p-nitroanilide, which was prepared in trifluoroethanol at a concentration of 3 mM, was added to the cuvette. The contents of the cuvette were rapidly mixed and the absorbance at 390 nm was measured over a period of 3 minutes in order to monitor the formation of p-nitroaniline upon chymotrypsin cleavage of the substrate peptide, which only occurs when the proline residue is in the trans-configuration. CypA catalyzes the cis-trans isomerization of proline residues, thus the catalytic activity of CypA can be indirectly determined from the rate at which chymotrypsin cleaves the peptide substrate.
In Eqs. (5) and (6), [L] is the concentration of free ligand, Ptotal is the sum of the free and bound protein concentration, Ltotal is the sum of the free and bound ligand concentration, Kd is the dissociation constant of the CypA complex with the hit compound, which was determined in our binding affinity measurements (see above), and %Free is the fraction of CypA protein that is free in solution (i.e., not bound to ligand).
Single-point SUPREX is an abbreviated version of the SUPREX technique, which exploits the H/D exchange properties of globally protected amide protons in a protein to measure the increase in a protein’s thermodynamic stability upon ligand binding in solution. In a full SUPREX analysis, the extent to which a protein, or protein complexed with ligand, undergoes amide H/D exchange during a specified amount of time is determined as a function of denaturant concentration (see Figure 1). Ultimately, the denaturant concentration at the transition midpoint is used to ascertain the thermodynamic stability of the protein or protein-ligand complex. Proteins complexed with ligand have a SUPREX transition midpoint that is shifted to a higher denaturant concentration than that of the protein alone, and the magnitude of the shift is related to the binding affinity.
In single-point SUPREX the extent to which a protein, or a protein in the presence of ligand, undergoes amide H/D exchange during a specified amount of time is determined at a single denaturant concentration. The denaturant concentration is chosen such that there is a maximum difference between the ΔMass values expected for the protein and for the protein complexed with ligands having the minimum binding affinity desired in the selection. Based on the SUPREX curves expected for CypA in the presence and absence of a hypothetical ligand that binds CypA with a Kd value of 10 µM (see Figure 1), a 1.5 M denaturant concentration was chosen for the single-point SUPREX experiments in this work. This allowed for the selection of library compounds with Kd values ≤ 10 µM. We reasoned that ligands with Kd values up to ~10 µM would be useful lead compounds in our search for CypA-targeted diagnostic and imaging agents.
The 1280 and 9600 compounds in the LOPAC and DIVERSet libraries, respectively, were screened for binding to CypA using the protocol outlined in Figure 2. The protocol involved pooling ten library compounds per well of a microtiter plate prior to initiating the single-point SUPREX experiment. The library pools were prepared in 10 µl of deuterated buffer containing 1.5 M GdmCl and 10% DMSO, and the library compounds in each pool were present at a concentration of ~100 µM per compound. The single-point SUPREX protocol was initiated upon addition of CypA to each library pool. The CypA was allowed to undergo H/D exchange for 35 min in the presence of each library pool. Ultimately, the H/D exchange reactions were quenched upon the addition of a MALDI matrix solution, which also prepared the CypA protein in each sample for a MALDI-TOF analysis. The MALDI-TOF analysis was used to determine the number of amide protons in CypA that were exchanged with solvent deuterons during the 35-min H/D exchange reaction (i.e., the ΔMass value). The magnitude of the ΔMass value was used as described below to determine whether or not a given pool of library compounds contained binding ligand(s).
The pools of 10 library compounds, along with two sets of control samples, were screened for CypA binding using the single-point SUPREX protocol described above. The distributions of ΔMass values obtained in the analysis of the 1089 pools and 251 control samples are shown in Figure 3. The control samples included 126 positive controls that contained CsA, which is a known tight binding ligand to CypA, and 125 negative controls that did not contain a ligand. The data obtained on these control samples, which were analyzed with every set of 10 pools, was used to generate the cut-off values for hit selection and to help establish the false positive and false negative rates of the strategy employed in this work. Figure 4 shows the ΔMass values obtained for the positive and negative controls obtained during the analysis of each library.
In the LOPAC library screening, a cut-off value for the selection of hits in the first screening was set at 49.0 Da, which was 2.5 standard deviations below the average ΔMass value obtained for the negative controls. The ΔMass value of one negative control in the LOPAC library (see Figure 4) was determined using the Q-test at the 95% level to be an outlier18, and it was not used in the cut-off value calculation. The ΔMass value cut-offs used for the selection of hits in the DIVERSet library screening were calculated by subtracting 2.5 standard deviations from a 9-point moving average of the negative control data. Thus, a specific ΔMass value cut-off was determined for each set of ten pools that were analyzed in the DIVERSet library screening. We previously showed that the use of such a moving average in the calculation of ΔMass cut-off values is useful in single-point SUPREX screening experiments that are performed over the course of multiple days, as was done in the DIVERSet library screening.5
Ultimately, a list of hit-pools (i.e., those that yielded ΔMass values below the cut-off values) was generated. A total of 8 out of the 131 LOPAC library pools and 13 out of 958 DIVERSet Library pools were identified “hit-pools” in this first tier of screening. The hit-pools from each library were re-screened in a second tier of screening. We have previously shown that such a two-tier screening strategy significantly reduces the false positive rate using the single-point SUPREX assay.5 The two-tier screening strategy is useful because the assay is subject to random error, which is reflected in the Gaussian distributions observed for the ΔMass values recorded here (see Figure 3). The main sources of random error include the mass measurement uncertainty and the differential back-exchange (see below). After the second tier of screening, five hit-pools were identified from the LOPAC library, and five hit-pools were identified from the DIVERSet library.
The pooling strategy described here ultimately requires that separate analyses be performed on the individual compounds in the selected pools in order to determine which compound(s) yielded the “hit.” (i.e, were responsible for the low ΔMass value). Thus, the 100 compounds from the ten hit-pools identified from the LOPAC and DIVERSet libraries were individually screened for binding to CypA using the single-point SUPREX protocol. Ultimately, nine compounds were identified as CypA binding ligands in these individual screens (see Table 1). This included one compound from each of the ten hit-pools with the exception of one hit-pool that did not yield an individual hit. The binding affinity of these nine newly identified CypA ligands were then measured using the conventional SUPREX protocol.
Full SUPREX analyses were performed on CypA in the presence of the selected compounds in order to determine their affinity for CypA (Figure 5). The individual data points in a full SUPREX curve are subject to the same ΔMass value uncertainties as those in the single-point SUPREX protocol. However, protein-ligand binding analyses using the full SUPREX protocol are much less sensitive to this error because multiple ΔMass values are used to define a full SUPREX curve (see Figure 5). The full SUPREX analyses indicated that 8 of the 9 selected compounds yielded measurable binding interactions with CypA, i.e., a shift in the transition midpoint of the SUPREX curves. The SUPREX-derived Kd values of these 8 compounds ranged from 0.200–37 µM, with the majority being in the low micromolar range (Table 1).
Dissociation constants for two hit compounds, 5 and 8 from the DIVERSet library, could not be determined due to difficulties associated with obtaining MALDI ion signals for the CypA in the presence of these compounds at high [GdmCl] concentrations. At these elevated concentrations of denaturant, a MALDI ion signal for CypA was not observed and multiple peaks at m/z 200–1000 Da higher than that expected for CypA were detected. However, the ΔMass values obtained at the lower denaturant concentrations used in our full SUPREX analyses of CypA in the presence of these individual compounds were consistent with those expected for a pre-transition baseline, suggesting that binding was, in fact, detected, and that the transition midpoint of the curve was at least 1.6 and 2.3 M for 5 and 8, respectively. Binding was not detected in our full SUPREX analyses of CypA in the presence of one of the selected ligands, 6 (i.e., the SUPREX transition midpoints for CypA were similar in the presence and absence of this ligand), which appeared to be a false positive.
The nine newly identified CypA ligands were also tested for their ability to inhibit CypA’s peptidyl-prolyl isomerase activity. Initial velocities of the chymotrypsin-catalyzed hydrolysis of the peptide substrate N-succinyl-Ala-Ala-Pro-Phe-p-nitroanilide were determined in the presence and absence of CypA and then in the presence of CypA and each newly identified CypA ligand. The substrate peptide is cleaved by chymostrypsin when the Ala-Pro peptide bond is in the trans configuration. CypA increases the apparent rate that chymotrypsin hydrolyzes the scissile bond in this substrate by catalyzing the conversion of the cis isomer to the trans isomer. The degree to which CsA and the nine hit compounds identified in this study impact the rate of this hydrolysis reaction in the presence of CypA was determined from initial rate measurements (see Figure 6).
The chymotrypsin-coupled enzyme assays revealed that seven of the nine compounds (including Compounds 1–5, 7, and 9) inhibited CypA’s isomerase activity, and the degree to which each of these ligands inhibited CypA was generally consistent with the amount of free protein calculated based on each ligand’s measured Kd value, with the exception of Compound 8 that appears to bind CypA with relatively high affinity but have minimal effect on CypA’s isomerase activity. It is not surprising that such a compound was identified in this work, as the single-point SUPREX assay does not select for compounds with specific activities, rather just binding interactions. The variety of compounds that can be selected in the single-point SUPREX assay (e.g., both those that alter a specific protein function and those that do not) can be an advantage of the methodology over other HTS approaches for certain applications, such as the one described here to identify molecular imaging agents. Interestingly, Compound 6, which was not found to have a measurable CypA binding interaction in our full SUPREX analysis (Kd value > 40 µM), did not significantly inhibit the isomerase activity of CypA isomerase. Compound 6 is most likely a false positive.
The first tier of the 1280-compound LOPAC library screening was accomplished by researchers working over the course of one day for 3.5 hours, which was the time it took to perform the H/D exchange reactions and collect the mass spectral data. This equates to approximately 1.5 minutes per pool or less than 10 seconds/compound. The first tier of the 9600-compound Chembridge DIVERSet library screening was accomplished by two researchers working over the course of four days for a total of 22 hours, which again was the time it took to perform the H/D exchange reactions and collect the mass spectral data. This equates to a screening rate of approximately 1 min/pool or about 6 seconds/compound.
The overall throughput achieved in this work, 6 seconds/compound, is about 30 times faster than that previously reported for the single-point SUPREX protocol.5 In theory, the pooling strategy used here should increase the throughput of the original single-point SUPREX protocol by 10-fold. The additional 3-fold increase in throughput realized in this work is likely due to the use of a high-throughput MALDI-TOF instrument, which was equipped with a high repetition rate laser.
The false positive and false negative rates of a high-throughput screening assay can provide a measure of its efficiency. A total of 126 positive controls and 125 negative controls were analyzed during the single-point SUPREX screening experiments described here. Using the ΔMass cut-off values described above, false positive and false negative rates of 0% and 9% were determined from the control data, which are both identical to that previously observed for the control data in our Prestwick Chemical Library screening, which were 0 and 9%, respectively.5
During the screening of the pooled library samples, false negatives go undetected, and false positives are only detected after subsequent screens and/or analyses. A total of 10 hit-pools were initially identified after our analysis of the 1,089 10-compound pools. Subsequent analyses of the individual ligands from these 10 hit-pools identified hit-compounds from all but two pools, suggesting a false positive rate of 20%. Interestingly, one of the two false positive hit pools consistently produced low ΔMass values upon re-analysis, even though it failed to produce a hit-compound when the ligands from the pool were individually analyzed for CypA binding. This suggests that the higher occurrence of false positives in the pooled samples (i.e., ~20%) compared to the false positives in the controls (i.e., 0%), which were not pooled, may be a result of the pooling strategy. This is also supported by the observation that the 20% false positive rate determined from the pooled library compounds in this work was larger than the 0% false positive rate determined from the library compounds in our earlier screen of the Prestwick Chemical Library, which did not employ a pooling strategy.5 The 20% false positive rate in this work is in line with typical pharmaceutical screens, which have been estimated to have false positive rates of about 40%.19
A complicating issue in the pooling strategy is the potential for higher false positive rates due to promiscuous aggregation and non-specific binding of library compounds.20, 21 There have been reports of small molecules aggregating into particles with diameters ranging from 95 to 400 nm.21 Separate studies have shown that enzymes can be inhibited through reversible adsorption of enzyme onto the surface of such molecular aggregates.20 This phenomenon is worse at higher concentrations of small molecules relative to protein concentration. Sequestration of the protein of interest in a screening assay by aggregate particles would thus increase the protein’s stability due to non-specific binding. This could explain why some pools (with total small molecule concentrations of 1 mM in the analysis buffers) were selected as positive hits, but then none of the individual compounds from the pools would be selected as binding ligands when they were analyzed individually (with the total small molecule concentrations on the order of 100 µM).
A measure of assay robustness is the Z’-factor. The Z’-factor is maximized by increasing the difference between the ΔMass values of hits and non-hits and by minimizing the standard deviations of the ΔMass values. The difference between hit and non-hit ΔMass values is largely a property of the specific protein under study, and it is directly related to the number of amide protons that are globally protected in the protein’s three-dimensional structure. The standard deviations of the ΔMass values are limited by the mass spectrometer. The standard deviations obtained using the MALDI TOF mass spectrometer in this work averaged 3.7 Da during the screening of the LOPAC library and 3.4 Da in the DIVERSet. Generally, HTS assays with Z’-factors ≥ 0 are useful for screening large combinatorial libraries.13 The distribution of Z’ –factors determined from the positive and negative controls in this work is shown in Figure 4C.
The Z’ –factors values varied in this work mainly because the amplitude of the assay signal (i.e., ΔMass value difference between the positive and negative controls) varied over the course of the screening (see Figure 4B). This was a result of differential back-exchange (a random effect of not being able to determine the ΔMass values in the MALDI-TOF experiment at exactly the same time after they are quenched) and the fact that the standard deviations for the mass measurements of some controls were large and clearly outliers. While the moving average procedure accounted for this differential back-exchange to set the cut-off values in the screen, it is clear that the back-exchange reaction negatively impacts the Z’-value. While there was a large range of Z’-values (i.e., −1.3–0.7), the large majority of Z’-values were positive with most clustered between 0.1 and 0.5 (see Figure 4C) using the protocol described here. The distribution of Z’-factors observed here using our pooling strategy is very similar to that observed in our previous application of the single-point SUPREX assay that did not utilize a pooling strategy.5 In theory, the Z’-factors in the single-point assay should be normally distributed. The higher than expected number of low Z’-factors in here (Figure 4C), and in our previous work,5 is likely due to the back-exchange reaction, which compromise Z’-values by reducing the separation between hits and non-hits. These problems associated with the back-exchange reaction can be minimized by recording ΔMass values of the H/D exchanged protein samples in the MALDI-TOF experiment as soon after they are quenched as is possible.
Full SUPREX analyses of the nine hit-compounds selected in our screen (Figure 4) revealed that eight of the nine selected compounds bound CypA with Kd values in the low micromolar to high nanomolar range (Table 1). These eight hit-compounds are all novel CypA ligands. It is interesting to find that CypA binds to 4 (cisplatin) since this drug is a commonly used anti-tumor agent. CypA has been found to be overexpressed in many cancer types including non-small cell lung cancer1, 2, pancreatic cancer22, 23, endometrial carcinoma24, and oral cancer.25 In addition, CypA has been shown to be associated with the progression of breast cancer.26 The mode of action for cisplatin-induced apoptosis is not completely understood but is at least partly through the generation of reactive oxygen species.27 It has previously been observed that cells overexpressing CypA have been more resistant to cisplatin-induced death by minimizing stress-induced apoptosis, while offering no additional resistance to other anti-cancer drugs.28 We hypothesize that this may be the result of CypA binding to cisplatin as shown in our results, thus decreasing the amount of cisplatin available to the cell to generate reactive oxygen species.
It has also been shown that Compound 2 (Bay 11–7085) inhibits the tumor necrosis factor-α (TNFα)-induced phosphorylation of IκBα, the inhibitor form of nuclear factor kappa B (NFκB), which serves to decrease expression of the inflammatory molecules ICAM-1, VCAM-1, and E-selectin.29 The specific protein target of Bay 11–7085 was not determined in the study. However, it was speculated that the phosphorylation inhibition might be due to inhibition of a signaling protein upstream of the IκBα. It has been demonstrated that CypA has the ability to induce phosphorylation of IκBα in inflammation sites caused by rheumatoid arthritis30, and thus it is feasible to hypothesize that the mode of action of Bay 11–7085 involves CypA being inhibited by Bay 11–7085, which in turn decreases the phosphorylation of IκBα and in turn the activation of NFκB-activated expression of these inflammatory molecules. The other selected compounds from the LOPAC library have known roles. Compound 3 (ZM 39923 HCl) is a selective Janus Kinase 3 (JAK3) inhibitor, binding to JAK3’s ATP binding site.31 It has also been found to inhibit transglutamase during a HTS assay.32 Compound 1 (α-NETA) is a selective, fluorescent inhibitor of choline acetyltransferase and is used for the investigation of acetylcholine synthesis.33 To date, the hits from the Chembridge DIVERSet library have no previously known function to the best of our knowledge.
This study demonstrated the use of single-point SUPREX for small (<10,000 compounds) libraries. Using a pooling strategy, the throughput was increased more than ten-fold compared to what has been previously reported. However, one drawback of the pooling strategy was a higher false positive rate of 20%. While the pooling strategy did not specifically impact Z’-factors of the single-point SUPREX assay, it is important to recognize that the Z’-values determined in this work were generally between 0 and 0.5, with approximately one-third of the values being less than 0. As a result of these relatively low Z’ values, the assay requires at least a two-tier screening strategy. This was feasible for the relatively small (1000–10,000) compound libraries studied here. However, the relatively low Z’ values of the assay may limit applications to larger compound libraries. The results of this study confirmed the assay’s capacity to identify novel ligands to cyclophilin. In this case, eight novel CypA ligands were discovered, with one compound estimated to have a Kd of less than 200 nM. Additionally, the use of a traditional biochemical assay showed that seven of the eight novel ligands also inhibited the isomerase activity of CypA. The assay is unique among high-throughput screening assays in its generality (i.e., ligand selection is based solely on the introduction of structural stability and not necessarily on the modulation a protein’s specific biological activity).
This work was supported by an NIH grant to M.C.F. (R21-CA-127064). The authors are very grateful to Vladimir Poltoratsky for his help in pooling the DIVERSet library.