Determination Of Dissociation Constants Research Articles

Determination of Binding Constants and Gas Phase Stabilities of Artificial Carbohydrate Receptor Complexes Using Electrospray Mass Spectrometry.

In recent years, binding studies to determine complex stabilities and selectivities of artificial carbohydrate receptors with glycosides have been mainly performed using 1H NMR, isothermal titration calorimetry (ITC), and other spectroscopic titration techniques. Native electrospray ionization (ESI) mass spectrometry is used only to verify the complex stoichiometries, although determination of dissociation constants is also possible. Herein, the binding of a 1,3,5-substituted 2,4,6-triethylbenzene-based receptor (CHR) to four alkyl-β-d-glucosides with varying alkyl side chain lengths (methyl (MGP), hexyl (HGP), octyl (OGP), and dodecyl (DGP)-β-d-glucosides), which was analyzed by ESI Fourier transform ion cyclotron resonance mass spectrometry (ESI-FT-ICR-MS) under optimized spray conditions in both ion modes, is reported. The complexes of the receptor with different sugars could be detected in 1:1 and 2:1 stoichiometries. Dissociation constants calculated for the 1:1 complexes showed a stability trend depending on the length of the alkyl side chain of the sugar: CHR:DGP > CHR:OGP > CHR:HGP > CHR:MGP. Gas phase stabilities determined by CID-MS confirm this relative trend in binding affinities. These findings substantiate the validity and applicability of ESI-MS as a method for investigating noncovalent complex stabilities and thus support research in the field of molecular recognition of carbohydrates by artificial receptors.

Quantifying genome-wide transcription factor binding affinities for chromatin using BANC-seq.

Transcription factors (TFs) bind specific DNA sequences to regulate transcription. Apart from DNA sequences, local factors such as DNA accessibility and chromatin structure determine the affinity of a TF for any given locus. Including these factors when measuring TF-DNA affinities has proven difficult. To address this challenge, we recently developed a method called binding affinities in native chromatin by sequencing (BANC-seq). In BANC-seq, intact mammalian nuclei are incubated with a concentration range of epitope-tagged TF, followed by either chromatin immunoprecipitation or cleavage under target and release using nuclease with spike-in DNA. This allows determination of apparent dissociation constant (KdApp) values, defined by the concentration of TF at which half-maximum binding occurs, across the genome. Here we present a detailed stepwise protocol for BANC-seq, including downstream data analysis. In principle, any molecular biologist should be able to perform a BANC-seq experiment in as little as 1.5 d (excluding analysis). However, preprocessing and analysis of the sequencing data does require some experience in command-line shell and R programming.

A simulation-guided “swapping” protocol for NMR titrations to study protein–protein interactions

Solution nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for characterizing protein–protein interactions. NMR-monitored titration proves to be effective in determining dissociation constants, particularly for Kd values in the micromolar to millimolar range. In conventional NMR titrations, planning for optimal titration conditions requires a prior estimate of Kd. In addition, a highly concentrated ligand stock solution is often required, posing challenges when the ligand exhibits limited solubility and stability at elevated concentrations. To overcome these constraints, we propose a simulation-guided “swapping” protocol for NMR titrations. Guided by simulations of the binding curves, two samples, one with zero and the other with maximum ligand concentration, but both containing identical protein concentrations, initiate the titration. Using a “swapping” strategy, intermediate ligand concentrations in between those of the two initial samples are generated without the need of a concentrated ligand stock, while maintaining constant protein concentration. More importantly, this protocol facilitates estimation of Kd by early titration points and allows on-the-fly optimization of the titration points. The proposed approach enhances the efficiency of NMR titrations and provides a straightforward means to optimize the experimental conditions for the titrations.

Determination of dissociation constants of cephalosporin antibiotics by cellmetry method

Acid dissociation constants of three cephalosporin antibiotics (cefapirin, ceftiofur, and cefotaxime) were calculated by a newly developed methodology. Plane-wave DFT calculations were performed to determine the pKa values, and by choosing the appropriate cell sizes, accurate values could be calculated. Some characteristic points were found which helped us to find correlations among the structural and physic-chemical parameters, and correlation factors were defined as well. This present study can be a base for further approaches to determining acid dissociation constants of cephalosporin molecules.Graphical

Towards the "Eldorado" of pKa Determination: A Reliable and Rapid DFT Model.

The selection of a "perfect tool" for the theoretical determination of acid-base dissociation constants (Ka) is still puzzling. Recently, we developed a user-friendly model exploiting CAM-B3LYP for determining pKa with impressive reliability. Herein, a new challenge is faced, examining a panel of functionals belonging to different rungs of the "Jacob's ladder" organization, which classifies functionals according to their level of theory. Specifically, meta-generalized gradient approximations (GGAs), hybrid-GGAs, and the more complex range-separated hybrid (RSH)-GGAs were investigated in predicting the pKa of differently substituted carboxylic acids. Therefore, CAM-B3LYP, WB97XD, B3PW91, PBE1PBE, PBEPBE and TPSSTPSS were used, with 6-311G+(d,p) as the basis set and the solvation model based on density (SMD). CAM-B3LYP showed the lowest mean absolute error value (MAE = 0.23) with relatively high processing time. PBE1PBE and B3PW91 provided satisfactory predictions (MAE = 0.34 and 0.38, respectively) with moderate computational time cost, while PBEPBE, TPSSTPSS and WB97XD led to unreliable results (MAE > 1). These findings validate the reliability of our model in predicting carboxylic acids pKa, with MAE well below 0.5 units, using a simplistic theoretical level and a low-cost computational approach.

Speeding-up the Determination of Protein-Ligand Affinities by STD NMR: The Reduced Data Set STD NMR Approach (rd-STD NMR).

STD NMR spectroscopy is a powerful ligand-observed NMR tool for screening and characterizing the interactions of small molecules and low molecular weight fragments with a given macromolecule, identifying the main intermolecular contacts in the bound state. It is also a powerful analytical technique for the accurate determination of protein-ligand dissociation constants (KD) of medium-to-weak affinity, of interest in the pharmaceutical industry. However, accurate KD determination and epitope mapping requires a long series of experiments at increasing saturation times to carry out a full analysis using the so-called STD NMR build-up curve approach and apply the "initial slopes approximation". Here, we have developed a new protocol to bypass this important limitation, which allows us to obtain initial slopes by using just two saturation times and, hence, to very quickly determine precise protein-ligand dissociation constants by STD NMR.

The effect of TG2-inhibitory monoclonal antibody zampilimab on tissue fibrosis in human in vitro and primate in vivo models of chronic kidney disease.

Fibrotic remodeling is the primary driver of functional loss in chronic kidney disease, with no specific anti-fibrotic agent available for clinical use. Transglutaminase 2 (TG2), a wound response enzyme that irreversibly crosslinks extracellular matrix proteins causing dysregulation of extracellular matrix turnover, is a well-characterized anti-fibrotic target in the kidney. We describe the humanization and characterization of two anti-TG2 monoclonal antibodies (zampilimab [hDC1/UCB7858] and BB7) that inhibit crosslinking by TG2 in human in vitro and rabbit/cynomolgus monkey in vivo models of chronic kidney disease. Determination of zampilimab half-maximal inhibitory concentration (IC50) against recombinant human TG2 was undertaken using the KxD assay and determination of dissociation constant (Kd) by surface plasmon resonance. Efficacy in vitro was established using a primary human renal epithelial cell model of tubulointerstitial fibrosis, to assess mature deposited extracellular matrix proteins. Proof of concept in vivo used a cynomolgus monkey unilateral ureteral obstruction model of chronic kidney disease. Zampilimab inhibited TG2 crosslinking transamidation activity with an IC50 of 0.25 nM and Kd of <50 pM. In cell culture, zampilimab inhibited extracellular TG2 activity (IC50 119 nM) and dramatically reduced transforming growth factor-β1-driven accumulation of multiple extracellular matrix proteins including collagens I, III, IV, V, and fibronectin. Intravenous administration of BB7 in rabbits resulted in a 68% reduction in fibrotic index at Day 25 post-unilateral ureteral obstruction. Weekly intravenous administration of zampilimab in cynomolgus monkeys with unilateral ureteral obstruction reduced fibrosis at 4 weeks by >50%, with no safety signals. Our data support the clinical investigation of zampilimab for the treatment of kidney fibrosis.

Molecular Mechanism of STIL Coiled-Coil Domain Oligomerization.

Coiled-coil domains (CCDs) play key roles in regulating both healthy cellular processes and the pathogenesis of various diseases by controlling protein self-association and protein-protein interactions. Here, we probe the mechanism of oligomerization of a peptide representing the CCD of the STIL protein, a tetrameric multi-domain protein that is over-expressed in several cancers and associated with metastatic spread. STIL tetramerization is mediated both by an intrinsically disordered domain (STIL400-700) and a structured CCD (STIL CCD718-749). Disrupting STIL oligomerization via the CCD inhibits its activity in vivo. We describe a comprehensive biophysical and structural characterization of the concentration-dependent oligomerization of STIL CCD peptide. We combine analytical ultracentrifugation, fluorescence and circular dichroism spectroscopy to probe the STIL CCD peptide assembly in solution and determine dissociation constants of both the dimerization, (KD = 8 ± 2 µM) and tetramerization (KD = 68 ± 2 µM) of the WT STIL CCD peptide. The higher-order oligomers result in increased thermal stability and cooperativity of association. We suggest that this complex oligomerization mechanism regulates the activated levels of STIL in the cell and during centriole duplication. In addition, we present X-ray crystal structures for the CCD containing destabilising (L736E) and stabilising (Q729L) mutations, which reveal dimeric and tetrameric antiparallel coiled-coil structures, respectively. Overall, this study offers a basis for understanding the structural molecular biology of the STIL protein, and how it might be targeted to discover anti-cancer reagents.

Targeting SARS-CoV-2 Macrodomain-1 to Restore the Innate Immune Response Using In Silico Screening of Medicinal Compounds and Free Energy Calculation Approaches.

Among the different drug targets of SARS-CoV-2, a multi-domain protein known as NSP3 is a critical element of the translational and replication machinery. The macrodomain-I, in particular, has been reported to have an essential role in the viral attack on the innate immune response. In this study, we explore natural medicinal compounds and identify potential inhibitors to target the SARS-CoV-2-NSP3 macrodomain-I. Computational modeling and simulation tools were utilized to investigate the structural-dynamic properties using triplicates of 100 ns MD simulations. In addition, the MM/GBSA method was used to calculate the total binding free energy of each inhibitor bound to macrodomain-I. Two significant hits were identified: 3,5,7,4'-tetrahydroxyflavanone 3'-(4-hydroxybenzoic acid) and 2-hydroxy-3-O-beta-glucopyranosyl-benzoic acid. The structural-dynamic investigation of both compounds with macrodomain-I revealed stable dynamics and compact behavior. In addition, the total binding free energy for each complex demonstrated a robust binding affinity, of ΔG -61.98 ± 0.9 kcal/mol for Compound A, while for Compound B, the ΔG was -45.125 ± 2.8 kcal/mol, indicating the inhibitory potential of these compounds. In silico bioactivity and dissociation constant (KD) determination for both complexes further validated the inhibitory potency of each compound. In conclusion, the aforementioned natural products have the potential to inhibit NSP3, to directly rescue the host immune response. The current study provides the basis for novel drug development against SARS-CoV-2 and its variants.

Synthesis and Determination of Acid Dissociation Constants of Bis-Acyl Thiourea Derivatives

N,N'-((dodecane-1,12-diylbis(azanediyl))bis(carbonothioyl))bis(4-nitrobenzamide) 5 and N,N'-((dodecane-1,12-diylbis(azanediyl))bis(carbonothioyl))bis(3-nitrobenzamide) 6 as bis acyl thiourea derivatives were synthesized and their molecular structures were characterized using 1H NMR, 13C NMR, COSY, DEPT, HMQC, FT-IR, and HRMS techniques. The acid dissociation constants (pKa) of the bis-acyl thiourea derivatives 5, 6 were determined potentiometrically and spectrophotometrically. The pKa values of products 5, 6 were determined in 50% (v/v) dimethyl sulfoxide–water hydro-organic solvent in the presence of 0.1 molL1 ionic strength of NaCl and in the acidic medium at 25±0.1 °C, and two pKa values were calculated for each compound with the HYPERQUAD computer program using the data obtained from the potentiometric titrations performed. In addition, three pKa values for each compound were determined in the calculations using the HypSpec program from the data obtained from the spectrophotometric titrations performed under the conditions where the potentiometric titrations were performed. For compounds 5 and 6, spectrophotometrically, pKa1 was 3.56±0.08 and 3.87±0.01, respectively, pKa2 was 7.11±0.08 and 7.05±0.01, respectively, and pKa3 was 12.30±0.08 and 11.82±0.02, respectively. It can be said that pKa1, pKa2, and pKa3 values may belong to enol, enthiol, and NH, respectively. Moreover, for compounds 5 and 6, potentiometrically, pKa2 was 7.06±0.13 and 6.94±0.11, respectively, and pKa3 was 12.11±0.06 and 11.17±0.06, respectively, and it can be said that pKa2 and pKa3 values may belong to enthiol and NH, respectively. It is seen that the pKa values determined spectrophotometrically and potentiometrically are compatible with each other.

Imaging Saturation Transfer Difference (STD) NMR: Affinity and Specificity of Protein-Ligand Interactions from a Single NMR Sample.

We have combined saturation transfer difference NMR (STD NMR) with chemical shift imaging (CSI) and controlled concentration gradients of small molecule ligands to develop imaging STD NMR, a new tool for the assessment of protein-ligand interactions. Our methodology allows the determination of protein-ligand dissociation constants (KD) and assessment of the binding specificity in a single NMR tube, avoiding time-consuming titrations. We demonstrate the formation of suitable and reproducible concentration gradients of ligand along the vertical axis of the tube, against homogeneous protein concentration, and present a CSI pulse sequence for the acquisition of STD NMR experiments at different positions along the sample tube. Compared to the conventional methodology in which the [ligand]/[protein] ratio is increased manually, we can perform STD NMR experiments at a greater number of ratios and construct binding epitopes in a fraction (∼20%) of the experimental time. Second, imaging STD NMR also allows us to screen for non-specific binders, by monitoring any variation of the binding epitope map at increasing [ligand]/[protein] ratios. Hence, the proposed method does carry the potential to speed up and smooth out the drug discovery process.

Two‐ and Three‐Dimensional Strategies for Determining Consecutive Acidity Constants of Ciprofloxacin and Enrofloxacin Using Fluorescent Excitation‐Emission Measurements

AbstractMultidimensional analysis is getting more popular by the day thanks to technical improvements of analytical instruments and development of chemometric methods. Here, we describe the application of two different chemometric methods to the analysis of three‐dimensional spectrofluorimetric dataset for green determination of dissociation constants of ciprofloxacin and enrofloxacin with minimal use of consumables. Three‐dimensional spectrofluorimetric datasets of both drugs were obtained by recording excitation‐emission matrices of only six sample solutions having different pH values. Parallel factor analysis (PARAFAC) and multivariate curve resolution (MCR) methods with the use of three components were used to decompose three‐dimensional excitation‐emission‐pH arrays. The number of components support the presence of two pKa values. The individual profiles of H2A, HA− and A2− species in excitation, emission and pH modes were obtained by the application of PARAFAC and MCR. The intersection points in the pH modes were used to predict pKa values. pKa1 and pKa2 of ciprofloxacin was reported as 6.35, 9.15 and 6.42, 9.40 by PARAFAC and MCR‐ALS, respectively. For enrofloxacin, consecutive pKa values were predicted as 6.11 and 9.06 by PARAFAC and 6.10 and 9.10 by MCR. A detailed comparison to literature was done for the estimated pKa values.

Riboflavin kinase and pyridoxine 5'-phosphate oxidase complex formation envisages transient interactions for FMN cofactor delivery.

Enzymes catalysing sequential reactions have developed different mechanisms to control the transport and flux of reactants and intermediates along metabolic pathways, which usually involve direct transfer of metabolites from an enzyme to the next one in a cascade reaction. Despite the fact that metabolite or substrate channelling has been widely studied for reactant molecules, such information is seldom available for cofactors in general, and for flavins in particular. Flavin adenine dinucleotide (FAD) and flavin mononucleotide (FMN) act as cofactors in flavoproteins and flavoenzymes involved in a wide range of physiologically relevant processes in all type of organisms. Homo sapiens riboflavin kinase (RFK) catalyses the biosynthesis of the flavin mononucleotide cofactor, and might directly interplay with its flavin client apo-proteins prior to the cofactor transfer. Non-etheless, none of such complexes has been characterized at molecular or atomic level so far. Here, we particularly evaluate the interaction of riboflavin kinase with one of its potential FMN clients, pyridoxine-5'-phosphate oxidase (PNPOx). The interaction capacity of both proteins is assessed by using isothermal titration calorimetry, a methodology that allows to determine dissociation constants for interaction in the micromolar range (in agreement with the expected transient nature of the interaction). Moreover, we show that; i) both proteins become thermally stabilized upon mutual interaction, ii) the tightly bound FMN product can be transferred from RFK to the apo-form of PNPOx producing an efficient enzyme, and iii) the presence of the apo-form of PNPOx slightly enhances RFK catalytic efficiency. Finally, we also show a computational study to predict likely RFK-PNPOx binding modes that can envisage coupling between the FMN binding cavities of both proteins for the potential transfer of FMN.

Determination of dissociation constants via quantitative mass spectrometry

The interplay of biomolecules governs all cellular processes. Qualitative analysis of such interactions between biomolecules as well as the quantitative assessment of their binding affinities are essential for the understanding of biochemical mechanisms. As scientific interest therefore moves beyond pure structural investigation, methods that allow for the investigation of such interactions become increasingly relevant. In this perspective we outline classical methods that are applicable for the determination of binding constants and highlight specifically mass spectrometry based methods. The use of mass spectrometry to gain quantitative information about binding affinities however is a still developing field. Here, we discuss different approaches, which emerged over the last years to determine dissociation constants (KD) with mass spectrometry based methods. Specifically, we highlight the recent development of quantitative Laser Induced Liquid Bead Ion Desorption (qLILBID) mass spectrometry for the example of double stranded deoxyribonucleic acids as well as for different RNA—RNA binding protein systems. We show that quantitative laser induced liquid bead ion desorption can successfully be used for the top down investigation of complexes and their dissociation constants values ranging from low nM to low µM affinities.

Streamlined Data Processing for Determination of Equilibrium Dissociation Constants with Accurate Constant via Transient Incomplete Separation (ACTIS)

The determination of accurate equilibrium dissociation constants, Kd, of protein-small molecule complexes is important but challenging as all established methods have inherent sources of inaccuracy. Accurate Constant via Transient Incomplete Separation (ACTIS) is a new method for Kd determination using transient incomplete separation of the complex from the unbound small molecule in a pressure-driven flow inside a capillary. ACTIS is accurate, and its accuracy is invariant to variations in geometries of both the fluidic system and the flow. Furthermore, ACTIS is implemented using a simple fluidic system supporting its accuracy and providing a simple-to-follow/copy template for instrumentation. Despite the simple and robust instrumentation/acquisition, the current data processing workflow is cumbersome, time consuming, and prone to hard-to-trace human errors therefore hindering ACTIS' ability to become a practical reference method for Kd determination. This technical note describes a streamlined workflow for processing ACTIS data; the workflow is implemented as a set of open-source software tools called prACTISed (https://github.com/prACTISedProgram/prACTISed). These tools allow all steps of data processing to be performed in a fast and straightforward fashion. These practical software tools complement the simple instrumentation serving both developers and users of ACTIS.

Deconvoluting binding sites in amyloid nanofibrils using time-resolved spectroscopy.

Steady-state fluorescence spectroscopy has a central role not only for sensing applications, but also in biophysics and imaging. Light switching probes, such as ruthenium dipyridophenazine complexes, have been used to study complex systems such as DNA, RNA, and amyloid fibrils. Nonetheless, steady-state spectroscopy is limited in the kind of information it can provide. In this paper, we use time-resolved spectroscopy for studying binding interactions between amyloid-β fibrillar structures and photoluminescent ligands. Using time-resolved spectroscopy, we demonstrate that ruthenium complexes with a pyrazino phenanthroline derivative can bind to two distinct binding sites on the surface of fibrillar amyloid-β, in contrast with previous studies using steady-state photoluminescence spectroscopy, which only identified one binding site for similar compounds. The second elusive binding site is revealed when deconvoluting the signals from the time-resolved decay traces, allowing the determination of dissociation constants of 3 and 2.2 μM. Molecular dynamic simulations agree with two binding sites on the surface of amyloid-β fibrils. Time-resolved spectroscopy was also used to monitor the aggregation of amyloid-β in real-time. In addition, we show that common polypyridine complexes can bind to amyloid-β also at two different binding sites. Information on how molecules bind to amyloid proteins is important to understand their toxicity and to design potential drugs that bind and quench their deleterious effects. The additional information contained in time-resolved spectroscopy provides a powerful tool not only for studying excited state dynamics but also for sensing and revealing important information about the system including hidden binding sites.

Fluorescence Anisotropy and Polarization in the Characterization of Biomolecular Association Processes and Their Application to Study SH2 Domain Binding Affinity.

Fluorescence anisotropy (or polarization) is a powerful technique to study biomolecular association processes, by following the rotational motions of one of the two partners in the interaction, labeled with a fluorophore. It can be used to determine dissociation constants in solution, down to nM values, and unlabeled ligands can be characterized, too, by using competition experiments. In this chapter, we introduce the basic principles of the technique, compare it with other experimental approaches, and discuss the experimental details with specific examples regarding SH2 domain/phosphopeptide association processes. The experimental protocols to be used in binding experiments and displacement studies are described, as well as the caveats to be considered in performing accurate measurements.

Transient Incomplete Separation of Species with Close Diffusivity to Study the Stability of Affinity Complexes

Large molecules can be generically separated from small ones, though partially and temporarily, in a pressure-driven flow inside a capillary. This transient incomplete separation has been only applied to species with diffusion coefficients different by at least an order of magnitude. Here, we demonstrate, for the first time, the analytical utility of transient incomplete separation for species with close diffusion coefficients. First, we prove in silico that even a small difference in diffusivity can lead to detectable transient incomplete separation of species. Second, we use computer simulation to prove that such a separation can be used for the reliable determination of equilibrium dissociation constant (Kd) of complexes composed of similar-sized molecules. Finally, we demonstrate experimentally the use of this separation for the accurate determination of Kd value for a protein-aptamer complex. We conclude that "accurate constant via transient incomplete separation" (ACTIS) can serve as a reference method for affinity characterization of protein-aptamer binding in solution.

Label-Free Determination of the Kinetic Parameters of Protein-Aptamer Interaction by Surface Plasmon Resonance.

Due to their high specificity and affinity to target molecules, aptamers can be used as powerful tools in diagnostics, therapeutics, and environmental or food analytics. For the use in various applications, the detailed characterization of their binding behavior is an important step after selection to determine the interaction strength between the aptamer and its target and to find the best kinetics depending on the field of application. The surface plasmon resonance (SPR) spectroscopy is a powerful technology to investigate important parameters in molecular interaction, for example, kinetics, affinity, and specificity. The most-used system is the Biacore™ SPR system which comprises an optical biosensor for label-free monitoring of binding events in real time based on SPR. This biophysical phenomenon describes the changes in refractive index on a sensor surface which can be used to measure binding events and to determine kinetic constants. In this chapter, a detailed protocol for the determination of kinetic constants for protein-aptamer interaction is provided. An 82-nt long ssDNA aptamer which are targeted against human urokinase is used as a model system for determination of binding and dissociation constants using Biacore™ SPR technology. A detailed note section provides useful tips and pitfalls at the end of this chapter.

In Situ Assay of Interfacial Interaction between ZnO Nanoparticles and Live Cell Disturbed by Surfactants.

The interfacial interaction between pollutants and organisms is a critical process in controlling the environmental fates of pollutants; however, in situ assay of the interaction is still a great challenge. Here, in situ determination of dissociation constants (Kd) for ZnO nanoparticles (ZnO NPs) from live algal cells disturbed by different-charged surfactants was established using microscale thermophoresis (MST). Moreover, in situ measurement of the adhesion force between the ZnO NPs probe and live single cell was performed using an atomic force microscope (AFM). Results showed that the cationic cetyltrimethylammonium chloride (CTAC) and anionic sodium dodecylbenzenesulfonate (SDBS) increased but nonionic Triton X-100 (TX-100) decreased the adhesion of ZnO NPs on cells. However, the force signature exhibited a smooth single retracted peak at short distances in the SDBS- and TX-100-treated groups, distinguished from the "see-saw" pattern peak in the CTAC-treated groups. The extended Derjaguin-Landau-Verway-Overbeek (XDLVO) calculation further confirmed that SDBS and TX-100 mainly disturbed the short-range hydration on the NP-cell interface, while CTAC reduced the long-range electrostatic repulsion. Furthermore, an excellent linear correlation between Zn bioaccumulation and two parameters (Kd and adhesion force) indicated that NP-cell interfacial interactions affected Zn bioaccumulation. Thus, in situ assay provides a quantitative basis for the pollutant-organism interfacial interaction to evaluate the environmental fate and ecological risk of pollutants.