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
-
Abstract We present a method to use long‐range CH coupling constants to derive the correct diastereoisomer from the molecular constitution of small molecules. A set of 792JCHand3JCHvalues collected from a single HSQMBC experiment on a sample of strychnine were used in the CASE‐3D (computer‐assisted 3D structure elucidation) protocol. In addition to the most commonly used3JCHcoupling constants, the subset of 322JCHvalues alone showed an excellent degree of configuration selection. The study is mainly based on comparison of DFT‐calculated2,3JCHvalues with experimental ones, critical for the case of2JCH. But the configuration selection also works well using3JCHvalues predicted from a semi‐empirical Karplus‐based equation limited to H−C−C−C fragments. The robustness, shown using strychnine as a proof of concept, makes theJ‐based CASE‐3D analysis a viable option for the application in fields such as peptide and carbohydrate research, organic synthesis, natural‐product identification and analysis, as well as medicinal chemistry.more » « less
