An official website of the United States government
Abstract To implement equilibrium hard‐modeling of spectroscopic titration data, the analyst must make a variety of crucial data processing choices that address negative absorbance and molar absorptivity values. The efficacy of three such methodological options is evaluated via high‐throughput Monte Carlo simulations, root‐mean‐square error surface mapping, and two mathematical theorems. Accuracy of the calculated binding constant values constitutes the key figure of merit used to compare different data analysis approaches. First, using singular value decomposition to filter the raw absorbance data prior to modeling often reduces the number of negative values involved but has little effect on the calculated binding constant despite its ability to address spectrometer noise. Second, both truncation of negative molar absorptivity values and the fast nonnegative least squares algorithms are superior to unconstrained regression because they avoid local minima; however, they introduce bias into the calculated binding constants in the presence of negative baseline offsets. Finally, we establish two theorems showing that negative values are best addressed when all the chemical solutions leading to the raw absorbance data are the result of mixing exactly two distinct stock solutions. This allows the raw absorbance data to be shifted up, eliminating negative baseline offsets, without affecting the concentration matrix, residual matrix, or calculated binding constants. Otherwise, the data cannot be safely upshifted. A comprehensive protocol for analyzing experimental absorbance datasets with is included.
more »
« less
A reliable algorithm for calculating stoichiometry parameters in the hard modeling of spectrophotometric titration data
Abstract Computation of binding constants from spectrophotometric titration data is a very popular application of chemometric hard modeling. However, the calculated values are misleading if the correct binding model is not used. Given that many supramolecular systems of interest feature unknown speciation, a priori determination of binding stoichiometry constitutes an important unsolved problem in chemometrics. We present a new and reliable algorithm for accomplishing this task, implemented using a hybrid particle swarm optimization technique. Simultaneous optimization of stoichiometry ratios and binding constants allows the optimal binding model to be calculated in just a few minutes for systems with up to four reactions. Simulated data studies demonstrate that the algorithm finds the correct stoichiometry with up to nine reactions in the absence of noise, including accurately determining species with unusual stoichiometry, such as H2G5. Application to four experimental datasets shows the algorithm is robust to experimental errors for a variety of chemical systems and binding models. This algorithm will facilitate the discovery of complex binding models, increase efficiency in titration analysis, and avert incorrect stoichiometry models, thereby improving the reliability of binding constant information in spectrophotometric titrations.
more »
« less
- Award ID(s):
- 2004005
- PAR ID:
- 10446403
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Journal of Chemometrics
- Volume:
- 36
- Issue:
- 6
- ISSN:
- 0886-9383
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Binding studies are ubiquitous in chemistry, but their extensive usefulness is undermined by false positive and false negative results. Centering on the G-protein mini-Gs, we present a thorough study with both simulated and experimental spectrophotometric titration data to diagnose the validity of both binding and non-binding models. Without the use of statistical tests like Bayesian Information Criterion (BIC) and data reconstruction fractions, spurious binding models may go undetected. Furthermore, if the signal change upon binding is too minute, false negatives can also result. Delineating such issues is paramount to effective science.more » « less
-
ABSTRACT Characterizing complicated solution phase systems in situ requires advanced modeling techniques to capture the intricate balances between the many chemical species. Due to the error inherent in any scientific measurement, a spectrophotometric titration experiment with nickel(II) and ethylenediamine (en) was repeated six times using an autotitrator to test the replicability of the data and the consistency of the resulting thermodynamic model. All six datasets could be modeled very tightly (R2 > 99.9999%) with the following eight complexes: [Ni]2+, [Ni2en]4+, [Nien]2+, [Ni2en3]4+, [Nien2]2+, [Ni2en5]4+, [Nien3]2+, and [Nien6]2+. The logK values for the stepwise associative reactions agree with existing literature values for the majority species ([Nienn = 1–3]2+) and matched expectations for the minority species; 95% confidence intervals for each logK value were determined via bootstrapping, which quantifies the variability in the binding constant value that is supported by a given dataset. The repeated experiments, which could not be successfully concatenated together, demonstrate that replication is crucial to capturing all the variability in the logK values. Conversely, bootstrapped confidence intervals across multiple experiments can be readily combined to generate an appropriate range for an experimentally determined binding constant.more » « less
-
Abstract Spacious M4L6tetrahedra can act as catalytic inhibitors for base‐mediated reactions. Upon adding only 5 % of a self‐assembled Fe4L6cage complex, the conversion of the conjugate addition between ethylcyanoacetate and β‐nitrostyrene catalyzed by proton sponge can be reduced from 83 % after 75 mins at ambient temperature to <1 % under identical conditions. The mechanism of the catalytic inhibition is unusual: the octacationic Fe4L6cage increases the acidity of exogenous water in the acetonitrile reaction solvent by favorably binding the conjugate acid of the basic catalyst. The inhibition only occurs for Fe4L6hosts with spacious internal cavities: minimal inhibition is seen with smaller tetrahedra or Fe2L3helicates. The surprising tendency of the cationic cage to preferentially bind protonated, cationic ammonium guests is quantified via the comprehensive modeling of spectrophotometric titration datasets.more » « less
-
Keskin, Ozlem (Ed.)Multistep protein-protein interactions underlie most biological processes, but their characterization through methods such as isothermal titration calorimetry (ITC) is largely confined to simple models that provide little information on the intermediate, individual steps. In this study, we primarily examine the essential hub protein LC8, a small dimer that binds disordered regions of 100+ client proteins in two symmetrical grooves at the dimer interface. Mechanistic details of LC8 binding have remained elusive, hampered in part by ITC data analyses employing simple models that treat bivalent binding as a single event with a single binding affinity. We build on existing Bayesian ITC approaches to quantify thermodynamic parameters for multi-site binding interactions impacted by significant uncertainty in protein concentration. Using a two-site binding model, we identify positive cooperativity with high confidence for LC8 binding to multiple client peptides. In contrast, application of an identical model to the two-site binding between the coiled-coil NudE dimer and the intermediate chain of dynein reveals little evidence of cooperativity. We propose that cooperativity in the LC8 system drives the formation of saturated induced-dimer structures, the functional units of most LC8 complexes. In addition to these system-specific findings, our work advances general ITC analysis in two ways. First, we describe a previously unrecognized mathematical ambiguity in concentrations in standard binding models and clarify how it impacts the precision with which binding parameters are determinable in cases of high uncertainty in analyte concentrations. Second, building on observations in the LC8 system, we develop a system-agnostic heat map of practical parameter identifiability calculated from synthetic data which demonstrates that the ability to determine microscopic binding parameters is strongly dependent on both the parameters themselves and experimental conditions. The work serves as a foundation for determination of multi-step binding interactions, and we outline best practices for Bayesian analysis of ITC experiments.more » « less
