skip to main content


Title: Properly handling negative values in the calculation of binding constants by physicochemical modeling of spectroscopic titration data
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
NSF-PAR ID:
10459182
Author(s) / Creator(s):
 ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Journal of Chemometrics
Volume:
33
Issue:
11
ISSN:
0886-9383
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The molar absorptivity curves for [NiLn(MeOH)6-n]2+ L = pyridine, 3-methylpyridine, 4-methylpyridine, n = 1 – 4, have been simultaneously deduced by modeling composite absorbance data of a series of equilibrium solutions in dry methanol using equilibrium-restricted factor analysis, a technique for obtaining spectral and thermodynamic information for component species involved in solution equilibria. Furthermore, the stepwise formation constants at 296 K have been determined with a high degree of accuracy. For pyridine, logK1-4 = 1.272(6), 0.669(9), 0.14(2), -0.32(2), respectively. For 3-methylpyridine, logK1-4 = 1.802(9), 1.16 (1), 0.32(1), -0.46(1), respectively. For 4-methylpyridine, logK1-4 = 2.808(9), 1.114(4), 0.411(4), -0.421(9), respectively. Unrestricted factor analysis was used to confirm the precise number of unique complexes in each case. The only additional complex for which some evidence was found was the pentakis version of the pyridine complex. 
    more » « less
  2. 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
  3. Abstract. Photochemical reactions of contaminants in snow and ice can be importantsinks for organic and inorganic compounds deposited onto snow from theatmosphere and sources for photoproducts released from snowpacks into theatmosphere. Snow contaminants can be found in the bulk ice matrix, ininternal liquid-like regions (LLRs), or in quasi-liquid layers (QLLs) at theair–ice interface, where they can readily exchange with the firn air. Somestudies have reported that direct photochemical reactions occur faster inLLRs and QLLs than in aqueous solution, while others have found similarrates. Here, we measure the photodegradation rate constants for loss of thethree dimethoxybenzene isomers under varying experimental conditions,including in aqueous solution, in LLRs, and at the air–ice interface ofnature-identical snow. Relative to aqueous solution, we find modestphotodegradation enhancements (3- and 6-fold) in LLRs for two of theisomers and larger enhancements (15- to 30-fold) at the air–ice interfacefor all three isomers. We use computational modeling to assess the impact oflight absorbance changes on photodegradation rate enhancements at theinterface. We find small (2–5 nm) bathochromic (red) absorbance shifts atthe interface relative to in solution, which increases light absorption, butthis factor only accounts for less than 50 % of the measured rate constantenhancements. The major factor responsible for photodegradation rateenhancements at the air–ice interface appears to be more efficientphotodecay: estimated dimethoxybenzene quantum yields are 6- to 24-foldlarger at the interface compared to in aqueous solution and account for themajority (51 %–96 %) of the observed enhancements. Using a hypotheticalmodel compound with an assumed Gaussian-shaped absorbance peak, we find thata shift in the peak to higher or lower wavelengths can have a minor tosubstantial impact on photodecay rate constants, depending on the originallocation of the peak and the magnitude of the shift. Changes in other peakproperties at the air–ice interface, such as peak width and height (i.e.,molar absorption coefficient), can also impact rates of light absorption anddirect photodecay. Our results suggest our current understanding ofphotodegradation processes underestimates the rate at which some compoundsare broken down, as well as the release of photoproducts into theatmosphere. 
    more » « less
  4. Abstract

    Long/multi‐pathlength spectrometers are being employed increasingly in ocean science, but detailed calibration/validation protocols for these instruments remain unavailable. The lack of such protocols has led to the collection of absorption data from different instruments that are not always comparable. The goal of this work was to develop procedures that improve inter‐laboratory data reproducibility and accuracy. A World Precision Instrument (WPI) multi‐pathlength liquid capillary waveguide spectrometer along with a Shimadzu dual‐beam spectrometer were employed to evaluate the proposed protocols, which address the following issues: (1) sample cross‐contamination due to improper cell cleaning, (2) nonstandardized sample preparation, (3) errors in effective pathlength determination, and (4) offsets due to refractive index effects (i.e., high‐salinity samples). Effective pathlength calculation is critical given that the measured and manufacturer provided pathlengths can differ significantly, with the provided pathlengths resulting in data with large variation (up to 21%) and percent error (up to 33%) when employing different pathlengths. When these proposed protocols were applied, we found that the WPI spectrometer's precision (measured as a coefficient of variation) was < 0.6% for intra‐day measurements and < 1.5% for inter‐day measurements, and there was < 2% error in absorbance measurements. These procedures allowed for the correction of the significant baseline offsets caused by differences in refractive index between reference and sample, by either subtracting a collected or a generated reference spectrum that matches the offset. Our results indicate that employing the procedures described herein will provide accurate, reproducible, and pathlength‐independent absorption data.

     
    more » « less
  5. Stable isotope proxies measured in the proteinaceous fraction of archaeological mollusc shell represents an increasingly important archive for reconstructing past ecological and biogeochemical conditions of nearshore environments. A major issue, however, is understanding the impact of diagenetic alteration in sub-fossil shell isotope values. “Bulk” stable isotope values of nitrogen (δ15N), and especially carbon (δ13C) often shift strongly with increasing C/N ratios in degraded shell, resulting in unreliable data. Here, we examine preservation of an entirely new set of shell paleo-proxies, compound-specific isotopes of amino acids (CSI-AA). We examine carbon (δ13CAA) and nitrogen (δ15NAA) patterns and values from the organic fraction of California mussel (Mytilus californianus) shells from the California Channel Islands. Archaeological shell samples ranging in age from ca. 6,100 to 250 cal BP exhibiting a wide range of degradation states were collected from varied depositional environments (e.g., exposed coastal bluff, buried strata, etc.), and were directly compared to modern shells of the same species and region. Our results indicate organic matter C/N ratios as the best bulk diagnostic indicator of the relative degradation state of shell organic fraction, including changes at the molecular level. Modern shell organic C/N ratios ranged from 2.8 to 3.5, while those in archaeological shell were substantially elevated (3.4–9.5), exhibiting strong and significant negative correlations with bulk δ13C values, weight %C, and weight %N, and a significant but weaker correlation with δ15 N values. An additional “cleaning” step using weak NaOH helped to remove possible exogenous contaminants and improved bulk values of some samples. However, relative molar AA abundances revealed that some AAs, especially the two most abundant, Glycine and Alanine, progressively decreased with increasing C/N ratio. The loss of these amino acids permanently alters bulk isotope values regardless of removal of contaminants. Modeling the bulk isotope change expected due to amino acid molar composition showed major and predictable shifts in bulk δ13C values from selected AA loss, and similarly large but far more variable impacts from exogenous contaminants. In contrast to bulk data, key CSI-AA values and patterns remained almost entirely unaltered, even in the most degraded shell samples, closely matching expected biosynthetic isotope patterns in modern mussel shell. AA isotope proxies for “baseline” (δ15N-Phenylalanine and average δ13C-Essential AAs) and planktonic trophic structure (δ15N-Glutamic Acid and δ15N-Phenylalanine) were not statistically altered with degradation in any sample. Overall, we conclude that while bulk isotopes, particularly δ13C, are very likely to be unreliable in archaeological or subfossil shell with C/N ratios higher than ∼4.0, CSI-AA proxies can still be used to reconstruct past climatic and ecological conditions of the nearshore marine environment. 
    more » « less