Two-dimensional turbulent flow chromatography coupled on-line to liquid chromatography–mass spectrometry for solution-based ligand screening against multiple proteins

We present herein a novel bioseparation/chemical analysis strategy for protein–ligand screening and affinity ranking in compound mixtures, designed to increase screening rates and improve sensitivity and ruggedness in performance. The strategy is carried out by combining on-line two-dimensional turbulent flow chromatography (2D-TFC) with liquid chromatography–mass spectrometry (LC–MS), and accomplished through the following steps: (1) a reversed-phase TFC stage to separate the protein/ligand complex from the unbound free molecules, (2) an on-line dissociation process to release the bound ligands from the complexes, and (3) a second mixed-mode cation-exchange/reversed-phase TFC stage to trap the bound ligands and to remove the proteins and salts, followed by LC–MS analysis for identification and determination of the binding affinities. The technique can implement an ultra-fast isolation of protein/ligand complex with the retention time of a complex peak in about 5 s, and on-line prepare the “clean” sample to be directly compatible with the LC–MS analysis. The improvement in performance of this 2D-TFC/LC–MS approach over the conventional approach has been demonstrated by determining affinity-selected ligands of the target proteins acetylcholinesterase and butyrylcholinesterase from a small library with known binding affinities and a steroidal alkaloid library composed of structurally similar compounds. Our results show that 2D-TFC/LC–MS is a generic and efficient tool for high-throughput screening of ligands with low-to-high binding affinities, and structure-activity relationship evaluation.     

Tunable transition metal-ligand complexation for enhanced elucidation of flavonoid diglycosides by electrospray ionization mass spectrometry

A tunable ESI-MS/MS strategy for differentiation of flavone and flavanone diglycoside isomers based on metal complexation with auxiliary ligands is reported. The addition of a metal salt and an auxiliary ligand to a flavonoid solution results in the formation of [M(II) (flavonoid-H) auxiliary ligand](+) complexes, where M(II) is a transition metal. A series of auxiliary ligands with electron-withdrawing substituents were synthesized to tailor the relative metal binding affinities of the ligands and thus directly influence the stabilities, and consequently the dissociation pathways, of the complexes. Upon collisionally activated dissociation, the complexes yield fragmentation patterns in which the abundances of key diagnostic ions are enhanced, thus facilitating isomer differentiation.     

Estimating protein-ligand binding affinity using high-throughput screening by NMR

Many of today's drug discovery programs use high-throughput screening methods that rely on quick evaluations of protein activity to rank potential chemical leads. By monitoring biologically relevant protein-ligand interactions, NMR can provide a means to validate these discovery leads and to optimize the drug discovery process. NMR-based screens typically use a change in chemical shift or line width to detect a protein-ligand interaction. However, the relatively low throughput of current NMR screens and their high demand on sample requirements generally makes it impractical to collect complete binding curves to measure the affinity for each compound in a large and diverse chemical library. As a result, NMR ligand screens are typically limited to identifying candidates that bind to a protein and do not give any estimate of the binding affinity. To address this issue, a methodology has been developed to rank binding affinities for ligands based on NMR screens that use 1D (1)H NMR line-broadening experiments. This method was demonstrated by using it to estimate the dissociation equilibrium constants for twelve ligands with the protein human serum albumin (HSA). The results were found to give good agreement with previous affinities that have been reported for these same ligands with HSA.     

Nonequilibrium capillary electrophoresis of equilibrium mixtures--a single experiment reveals equilibrium and kinetic parameters of protein-DNA interactions

We introduce a novel electrophoretic method, nonequilibrium capillary electrophoresis of equilibrium mixtures (NECEEM), and demonstrate its use for studying protein-DNA interactions. The equilibrium mixture of protein and DNA contains three components: free protein, free DNA, and the protein-DNA complex. A short plug of such a mixture is injected into the capillary, and the three components are separated under nonequilibrium conditions. The resulting electropherograms are composed of characteristic peaks and exponential curves. An easy nonnumerical analysis of a single electropherogram reveals two parameters: the equilibrium binding constant and the monomolecular rate constant of complex decay. The bimolecular rate constant of complex formation can then be calculated as the product of the two experimentally determined constants. NECEEM was applied to study the interaction between single-stranded DNA binding protein and a fluorescently labeled 15-mer oligonucleotide. It allowed us to measure for the first time the rate constant of complex decay for this important protein-DNA pair, k-1 = 0.03 s-1. The value of the equilibrium binding constant, Kb = 3.6 x 10-6 M-1, was in good agreement with those measured by other methods. As low as 10-18 mol of the protein was sufficient for the measurements. Thus, the new method is simple, informative, and highly sensitive. Moreover, it can be equally applied to other noncovalent protein-ligand complexes. These features of NECEEM make this method an indispensable tool in studies of macromolecular interactions. They also emphasize the potential role of NECEEM in the development of extremely sensitive protein assays using nucleotide aptamers.     

Evaluation of relative DNA binding affinities of anthrapyrazoles by electrospray ionization mass spectrometry

Binding interactions of a new series of anthrapyrazoles (APs) with DNA were evaluated by electrospray ionization mass spectrometry (ESI-MS). Relative binding affinities were estimated from the ESI-MS data based on the fraction of bound DNA for DNA/anthrapyrazole mixtures, and they show a correlation to the shift in melting point of the DNA measured from a previous study. Minimal sequence specificity was observed for the series of anthrapyrazoles. Upon collisionally activated dissociation of the duplex/anthrapyrazole complexes, typically ejection of the ligand was the dominant pathway for most of the complexes. However, for complexes containing AP2 or mitoxantrone, strand separation with the ligand remaining on one of the single strands was observed, indicative of a different binding mode or stronger binding.     

Screening of protein-ligand interactions using dynamic protein-affinity chromatography solid-phase extraction-liquid chromatography-mass spectrometry

A novel methodology is shown enabling the screening of mixtures of compounds for their affinity to a receptor protein. The system presented, dynamic protein-affinity chromatography solid-phase extraction (DPAC-SPE), overcomes the limitations of the existing methods by performing an incubation of the His-tagged protein with a mixture of possible ligands, which are still in their native conditions. This is followed by a fully automated affinity trapping step, coupled on-line to an LC-MS system in order to detect and identify the bound ligands. The system has been optimized using a commercially available on-line SPE system, using the estrogen receptor alpha (ERalpha) as model protein. A representative range of ligands with sub-nanomolar to millimolar affinities has been identified successfully from a mixture. The weakest binder that can be identified is norethindrone (approximately K(d)=0.1-1mM). The same setup also provides the possibilities to measure EC50 curves of both weak and strong binders.     

Affinity Screening Using Competitive Binding with Fluorine‐19 Hyperpolarized Ligands

Fluorine‐19 NMR and hyperpolarization form a powerful combination for drug screening. Under a competitive equilibrium with a selected fluorinated reporter ligand, the dissociation constant (K_D) of other ligands of interest is measurable using a single‐scan Carr–Purcell–Meiboom–Gill (CPMG) experiment, without the need for a titration. This method is demonstrated by characterizing the binding of three ligands with different affinities for the serine protease trypsin. Monte Carlo simulations show that the highest accuracy is obtained when about one‐half of the bound reporter ligand is displaced in the binding competition. Such conditions can be achieved over a wide range of affinities, allowing for rapid screening of non‐fluorinated compounds when a single fluorinated ligand for the binding pocket of interest is known.     

In situ formation of heterobimetallic salen complexes containing titanium and/or vanadium ions

A combination of high-resolution electrospray mass spectrometry and (1)H NMR spectroscopy has been used to prove that when a mixture of [(salen)TiO]2 complexes containing two different salen ligands (salen and salen') is formed, an equilibrium is established between the homodimers and the heterodimer [(salen)TiO2Ti(salen')]. Depending upon the structure and stereochemistry of the two salen ligands, the equilibrium may favor either the homodimers or the heterodimer. Extension of this process to mixtures of titanium(salen) complexes [(salen)TiO]2 and vanadium (V)(salen') complexes [(salen')VO] (+)Cl (-) allowed the in situ formation of the heterobimetallic complex [(salen)TiO2V(salen')] (+)X (-) to be confirmed for all combinations of salen ligands studied except when the salen ligand attached to titanium contained highly electron-withdrawing nitro-groups. The rate of equilibration between heterobimetallic complexes is faster than that between two titanium complexes as determined by line broadening in the (1)H NMR spectra. These structural results explain the strong rate-inhibiting effect of vanadium (V)(salen) complexes in asymmetric cyanohydrin synthesis catalyzed by [(salen)TiO]2 complexes. It has also been demonstrated for the first time that the titanium and vanadium complexes can undergo exchange of salen ligands and that this is catalyzed by protic solvents. However, the ligand exchange is relatively slow (occurring on a time scale of days at room temperature) and so does not complicate studies aimed at using heterobimetallic titanium and vanadium salen complexes as asymmetric catalysts. Attempts to obtain a crystal structure of a heterobimetallic salen complex led instead to the isolation of a trinuclear titanium(salen) complex, the formation of which is also consistent with the catalytic results obtained previously.     

Quantitative analysis of carbohydrate-protein interactions using glycan microarrays: determination of surface and solution dissociation constants

Carbohydrate-protein interactions on surface and in solution were quantitatively measured by a glycan microarray. Assessing carbohydrate affinities is typically difficult due to weak affinities and limited sources of structurally complex glycans. We described here a sensitive, high-throughput, and convenient glycan microarray technology for the simultaneous determination of a wide variety of parameters in a single experiment using small amounts of materials. Assay systems based on this technology were developed to analyze multivalent interactions and determine the surface dissociation constant (KD,surf) for surface-coated mannose derivatives with mannose binding lectins and antibodies. Competition experiments that employed monovalent ligands in solution yielded KD and Ki values in solution similar to equilibrium binding constants obtained in titration microcalorimetry and surface plasmon resonance experiments.     

Asymmetric Hydrogenation Using Rhodium Complexes Generated from Mixtures of Monodentate Neutral and Anionic Phosphorus Ligands

A series of monodentate neutral and anionic phosphorus ligands was synthesized and evaluated in the asymmetric rhodium‐catalyzed hydrogenation of functionalized olefins by using either catalysts containing identical ligands or catalysts generated from mixtures of two different ligands. We expected that the combination of an anionic ligand with a neutral ligand would favor the formation of hetero over homo bis‐ligand complexes due to charge repulsion. NMR spectroscopic studies confirmed that charge effects can indeed shift the equilibrium toward the hetero bis‐ligand complexes. In several cases, the combination of a neutral phosphane with an anionic phosphane, one chiral and the other achiral, furnished significantly higher enantioselectivities than analogous mixtures of two neutral ligands. The best results were obtained with a mixture of an anionic phosphoramidite and a neutral phosphoric acid diester. It is supposed that in this case a hydrogen bond between the two ligands additionally stabilizes the hetero ligand combination.     

Assays of ligand-human serum albumin binding using pulsed ultrafiltration and liquid chromatography-mass spectrometry

Two approaches were utilized to increase the throughput of pulsed ultrafiltration assays of ligand binding to human serum albumin, reducing the volume of the ultrafiltration chamber and combining pulsed ultrafiltration with high performance liquid chromatography-atmospheric pressure chemical ionization mass spectrometry (LC-MS). Affinity constants for binding of ligands to human serum albumin were determined using pulsed ultrafiltration with ultraviolet absorbance detection. The first affinity constants (Ka1) were measured for the binding of dansylsarcosine, dansylamide, 7-anilinocoumarin-4-acetic acid and warfarin, and were determined to be 1.8 x 105, 5 x 104, 8 x 104, and 2.0 x 105 M-1, respectively. The throughput of pulsed ultrafiltration analyses was tripled compared to previous pulsed ultrafiltration measurements by reducing the volume of the chamber. In addition, the use of LC-MS with pulsed ultrafiltration permitted the simultaneous comparison and rank ordering of ligand mixtures for binding to serum albumin. For example, the throughput of these pulsed ultrafiltration measurements was tripled by analyzing three ligands as a mixture.     

Mass spectrometry of protein-ligand complexes: Enhanced gas-phase stability of ribonuclease-nucleotide complexes

Noncovalent protein-ligand complexes are readily detected by electrospray ionization mass spectrometry (ESI-MS). Ligand binding stoichiometry can be determined easily by the ESI-MS method. The ability to detect noncovalent protein-ligand complexes depends, however, on the stability of the complexes in the gas-phase environment. Solution binding affinities may or may not be accurate predictors of their stability in vacuo. Complexes composed of cytidine nucleotides bound to ribonuclease A (RNase A) and ribonuclease S (RNase S) were detected by ESI-MS and were further analyzed by MS/MS. RNase A and RNase S share similar structures and biological activity. Subtilisin-cleavage of RNase A yields an S-peptide and an S-protein; the S-peptide and S-protein interact through hydrophobic interactions with a solution binding constant in the nanomolar range to generate an active RNase S. Cytidine nucleotides bind to the ribonucleases through electrostatic interactions with a solution binding constant in the micromolar range. Collisionally activated dissociation (CAD) of the 1:1 RNase A-CDP and CTP complexes yields cleavage of the covalent phosphate bonds of the nucleotide ligands, releasing CMP from the complex. CAD of the RNase S-CDP and CTP complexes dissociates the S-peptide from the remaining S-protein/nucleotide complex; further dissociation of the S-protein/nucleotide complex fragments a covalent phosphate bond of the nucleotide with subsequent release of CMP. Despite a solution binding constant favoring the S-protein/S-peptide complex, CDP/CTP remains electrostatically bound to the S-protein in the gas-phase dissociation experiment. This study highlights the intrinsic stability of electrostatic interactions in the gas phase and the significant differences in solution and gas-phase stabilities of noncovalent complexes that can result.     

Investigations of cyclophilin interactions with oligopeptides containing proline by affinity capillary electrophoresis

Affinity capillary electrophoresis using mobility-shift analysis was utilized to characterize the binding of peptide ligands to cyclophilins, which are members of the enzyme family of peptidyl-prolyl cis/trans isomerases. Peptides derived from the human immunodeficiency virus capsid protein p24 exhibited different affinities to the isoenzymes cyclophilin18 and cyclophilin20. For the interaction of the peptide hormone bradykinin with cyclophilin18, a dissociation constant of 1.4 +/- 0.1 mM was determined. Finally, the affinity of cyclophilin20 to peptides from a cellulose-bound peptide library scanning the sequence of Drosophila melanogaster protein cappuccino was investigated. The affinities of selected peptides to cyclophilin20 and a green fluorescent fusion protein with cyclophilin20 were compared.     

Identification of small molecule binding sites within proteins using phage display technology

Affinity selection of peptides displayed on phage particles was used as the basis for mapping molecular contacts between small molecule ligands and their protein targets. Analysis of the crystal structures of complexes between proteins and small molecule ligands revealed that virtually all ligands of molecular weight 300 Da or greater have a continuous binding epitope of 5 residues or more. This observation led to the development of a technique for binding site identification which involves statistical analysis of an affinity-selected set of peptides obtained by screening of libraries of random, phage-displayed peptides against small molecules attached to solid surfaces. A random sample of the selected peptides is sequenced and used as input for a similarity scanning program which calculates cumulative similarity scores along the length of the putative receptor. Regions of the protein sequence exhibiting the highest similarity with the selected peptides proved to have a high probability of being involved in ligand binding. This technique has been employed successfully to map the contact residues in multiple known targets of the anticancer drugs paclitaxel (Taxol), docetaxel (Taxotere) and 2-methoxyestradiol and the glycosaminoglycan hyaluronan, and to identify a novel paclitaxel receptor [1]. These data corroborate the observation that the binding properties of peptides displayed on the surface of phage particles can mimic the binding properties of peptides in naturally occurring proteins. It follows directly that structural context is relatively unimportant for determining the binding properties of these disordered peptides. This technique represents a novel, rapid, high resolution method for identifying potential ligand binding sites in the absence of three-dimensional information and has the potential to greatly enhance the speed of development of novel small molecule pharmaceuticals.     

Direct zonal liquid chromatographic method for the kinetic study of actinomycin-DNA binding

The binding of an anticancer drug (actinomycin D or ACTD) to double-stranded DNA (dsDNA) was studied by means of high-performance liquid chromatography (HPLC). ACTD is an antitumor antibiotic containing one chromophore group and two pentapeptidic lactone cycles that binds dsDNA. Incubations of ACTD with DNA were performed at physiological pH. The complexed and free ligand concentrations of the mixture were quantified at 440 nm from their separation on a size-exclusion chromatographic (SEC) column using the same buffer for the elution and the sample incubation. The DNA and the ACTD-DNA complexes were eluted at the column exclusion volume while the ligand was retained on the support. An apparent binding curve was obtained by plotting the amount emerging at the exclusion column volume against that eluted at free ACTD retention volume. A dissociating effect was evidenced and the binding parameters were significantly different from those obtained at equilibrium by visible absorbance titration. The equilibrium binding parameters determined by absorption spectroscopy were used as starting data in the numerical simulations of the chromatographic process. The results showed a strong dependency of the apparent binding parameters on the reaction kinetics. Finally the comparison of the apparent binding curve obtained from the HPLC experiments and from the numerical simulations permitted an evaluation of the dissociation rate constant (kd = 0.004 s(-1)).     

Experimental validation of theoretical potassium and sodium cation affinities of amides by mass spectrometric kinetic method measurements

In this study the theoretical Gaussian-2 K(+)/Na(+) binding affinities (enthalpies) at 0 K (in kJ mol(-1)) for six amides in the order: formamide (109.2/138.5) < N-methylformamide (117.7/148.6) < acetamide (118.7/149.5) < N,N-dimethylformamide (123.9/156.4) < N-methylacetamide (125.6/157.7) < N,N-dimethylacetamide (129.2/162.6), reported previously (Siu et al., J. Chem. Phys. 2001; 114: 7045-7051), were validated experimentally by mass spectrometric kinetic method measurements. By monitoring the collision-induced dissociation (CID) of K(+)/Na(+)-bound heterodimers of the amides, the relative affinities were shown to be accurate to within +/-2 kJ mol(-1). With these six theoretical K(+)/Na(+) binding affinities as reference values, the absolute K(+)/Na(+) affinities of imidazole, 1-methylimidazole, pyridazine and 1,2-dimethoxyethane were determined by the extended kinetic method, and found to be consistent (to within +/-9 kJ mol(-1)) with literature experimental values obtained by threshold-CID, equilibrium high-pressure mass spectrometry, and Fourier transform ion cyclotron resonance/ligand-exchange equilibrium methods. A self-consistent resolution is proposed for the inconsistencies in the relative order of K(+)/Na(+) affinities of amides reported in the literature. These two sets of validated K(+) and Na(+) affinity values are useful as reference values in kinetic method measurements of K(+)/Na(+) affinity of model biological ligands, such as the K(+) affinities of aliphatic amino acids.     

谷胱甘肽与D型氨基酸非共价复合物的质谱

为了研究谷胱甘肽和D型氨基酸的非共价相互作用, 将一定化学剂量比的还原型谷胱甘肽与D-苯丙氨酸、D-组氨酸或D-谷氨酰胺在室温下混合后, 温育1 h, 使反应达到平衡. 电喷雾质谱测量结果表明, 在生理pH条件下, 谷胱甘肽可以和D 型氨基酸经反应生成非共价复合物. 串级质谱中的碰撞诱导解离(CID)以及紫外光谱进一步确认了非共价复合物的生成. 为了避免严重的离子化效率和质谱信号相互抑制作用, 对谷胱甘肽和D-谷氨酰胺的相互作用作了定量的评估. 配制一系列不同初始浓度的谷胱甘肽和D型氨基酸的混合溶液, 并用电喷雾质谱测定混合溶液中不同质点的质谱峰强度, 计算了谷胱甘肽与三个D型氨基酸结合形成的复合物的解离常数. 计算结果表明, 谷胱甘肽可以和D型氨基酸结合形成不同键合强度的非共价复合物, 其稳定性按照D-谷氨酰胺、D-苯丙氨酸、D-组氨酸的次序逐渐增大.     

Toward quantitative estimates of binding affinities for protein–ligand systems involving large inhibitor compounds: A steered molecular dynamics simulation route

Understanding binding mechanisms between enzymes and potential inhibitors and quantifying protein – ligand affinities in terms of binding free energy is of primary importance in drug design studies. In this respect, several approaches based on molecular dynamics simulations, often combined with docking techniques, have been exploited to investigate the physicochemical properties of complexes of pharmaceutical interest. Even if the geometric properties of a modeled protein – ligand complex can be well predicted by computational methods, it is still challenging to rank with chemical accuracy a series of ligand analogues in a consistent way. In this article, we face this issue calculating relative binding free energies of a focal adhesion kinase, an important target for the development of anticancer drugs, with pyrrolopyrimidine‐based ligands having different inhibitory power. To this aim, we employ steered molecular dynamics simulations combined with nonequilibrium work theorems for free energy calculations. This technique proves very powerful when a series of ligand analogues is considered, allowing one to tackle estimation of protein – ligand relative binding free energies in a reasonable time. In our cases, the calculated binding affinities are comparable with those recovered from experiments by exploiting the Michaelis – Menten mechanism with a competitive inhibitor.     

Noncovalent interactions of Cu+ with N-donor ligands (pyridine, 4,4-dipyridyl, 2,2-dipyridyl, and 1,10-phenanthroline): collision-induced dissociation and theoretical studies

Collision-induced dissociation of complexes of Cu+ bound to a variety of N-donor ligands (N-L) with Xe is studied using guided ion beam tandem mass spectrometry. The N-L ligands examined include pyridine, 4,4-dipyridyl, 2,2-dipyridyl, and 1,10-phenanthroline. In all cases, the primary and lowest-energy dissociation channel observed corresponds to the endothermic loss of a single intact N-L ligand. Sequential dissociation of additional N-L ligands is observed at elevated energies for the pyridine and 4,4-dipyridyl complexes containing more than one ligand. Ligand exchange processes to produce Cu+Xe are also observed as minor reaction pathways in several systems. The primary cross section thresholds are interpreted to yield 0 and 298 K bond dissociation energies (BDEs) after accounting for the effects of multiple ion-neutral collisions, the kinetic and internal energy distributions of the reactants, and dissociation lifetimes. Density functional theory calculations at the B3LYP/6-31G* level are performed to obtain model structures, vibrational frequencies, and rotational constants for the neutral N-L ligands and the Cu+(N-L)x complexes. The relative stabilities of the various conformations of these N-L ligands and Cu+(N-L)x complexes as well as theoretical BDEs are determined from single point energy calculations at the B3LYP/6-311+G(2d,2p) level of theory using B3LYP/6-31G* optimized geometries. Excellent agreement between theory and experiment is observed for all complexes containing one or two N-L ligands, while theory systematically underestimates the strength of binding for complexes containing more than two N-L ligands. The ground-state structures of these complexes and the trends in the sequential BDEs are explained in terms of stabilization gained from sd-hybridization and repulsive ligand-ligand interactions. The nature of the binding interactions in the Cu+(N-L)x complexes are examined via natural bond orbital analyses.     

Electrospray ionization of nucleic acid aptamer/small molecule complexes for screening aptamer selectivity

Molecular recognition of small molecule ligands by the nucleic acid aptamers for tobramycin, ATP, and FMN has been examined using electrospray ionization mass spectrometry (ESI-MS). Mass spectrometric data for binding stoichiometry and relative binding affinity correlated well with solution data for tobramycin aptamer complexes, in which aptamer/ligand interactions are mediated by hydrogen bonds. For the ATP and FMN aptamers, where ligand interactions involve both hydrogen bonding and significant pi-stacking, the relative binding affinities determined by MS did not fully correlate with results obtained from solution experiments. Some high-affinity aptamer/ligand complexes appeared to be destabilized in the gas phase by internal Coulombic repulsion. In CAD experiments, complexes with a greater number of intermolecular hydrogen bonds exhibited greater gas-phase stability even in cases when solution binding affinities were equivalent. These results indicate that in at least some cases, mass spectrometric data on aptamer/ligand binding affinities should be used in conjunction with complementary techniques to fully assess aptamer molecular recognition properties.