Ionization pK,s for a large set of molecules were calculated using reactivity models developed in the computer program SPARC. SPARC uses relatively simple computational algorithms based on fundamental chemical structure theory to estimate ionization pK,s of organic molecules strictly from molecular structure. Molecular structures are broken at each essential single bond into functional units with intrinsic properties. Reaction centers (acid or base) are identified and the impact of appended molecular structure on ionization p& is quantified by perturbation theory. Resonance, electrostatic, solvation and H-bonding models have been developed and tested on 4338pK,s for 3685 compounds. The RMS deviation for the acids and the amino reaction center was 0.36 pK, units whereas that for the in-ring N and = N reaction centers was 0.41. Microscopic ionization constants, zwitterionic constants, isoelectric points, and molecular speciation as a function of pH can be calculated using the SPARC models. solubility, partitioning phenomena, and chemical reactivity are all highly dependent upon the state of ionization in the solution phase. The ionization pK, of an organic compound is vital to environmental exposure assessment because it can be used to define the degree of ionization and the propensity for sorption to soil and sediment by cation exchange. This, in turn, can determine mobility, reaction kinetics, bioavailability, complexation, etc. In addition to being highly significant in evaluating environmental fate and effects, acid-base ionization equilibria provide an excellent arena for testing the electrostatic effects models in the computer program SPARC. Because the gain or loss of protons results in a change in molecular charge, these processes are extremely sensitive to electric field effects within the molecule. The object of this study was to test these reactivity models utilizing the large pK, databases that are available. These databases are relatively reliable (f a few tenths of a pK, unit) for smaller more soluble molecules.Numerous investigators have attempted to predict ionization pK,s using various approaches such as ab initio [I, 21 and semiempirical [3,4] methods. The energy differences between the protonated state and the unprotonated state are small compared to the total binding energy of the reactant involved. This presents a problem for ab initio computational methods that calculate absolute energy values. Computing the relatively small energy differences needed for the analysis of chemical reactivity from absolute energies requires extremely accurate calculations. Hence, the aforementioned calculations were limited to a small subclass of molecules. A more aggressive attempt has been made by Klopman et. al., [5,6]. They estimated the p&s for some 2400 molecules (? = 0.846) based on QSAR using the Multi-CASE program. Despite the relatively large number of pK,s calculated their calculator was limited to the estimation of the first pK, [6]. Recently, we described our approach to predict numerous physical pr...
The prototype computer program SPARC has been under development for several years to estimate physical properties and chemical reactivity parameters of organic compounds strictly from molecular structure. SPARC solutesolute physical process models have been developed and tested for vapor pressure (at any temperature), heat of vaporization (at 25C and the boiling point), diffusion coefficient (at 25C) and boiling point (at any pressure) for a relatively large number of organic molecules. The RMS deviation error of the predicted the vapor pressures, heats of vaporization (at any temperature) and boiling points (at any pressure) were close to the intralaboratory experimental errors.
Solvation models, based on fundamental chemical structure theory, were developed in the SPARC mechanistic tool box to predict a large array of physical properties of organic compounds in water and in non-aqueous solvents strictly from molecular structure. The SPARC self-interaction solvation models that describe the intermolecular interaction between like molecules (solute-solute or solvent-solvent) were extended to quantify solute-solvent interaction energy in order to estimate the activity coefficient in almost any solvent. Solvation models that include dispersion, induction, dipole-dipole and hydrogen bonding interactions are used to describe the intermolecular interaction upon placing an organic solute molecule in any single or mixed solvent system. In addition to estimation of the activity coefficient for 2674 organic compounds, these solvation models were validated on solubility and liquid/liquid distribution coefficient in more than 163 solvents including water. The RMS deviations of the calculated versus observed activity coefficients, solubilities and liquid/liquid distribution coefficients were 0.272 log mole fraction, 0.487 log mole fraction and 0.44 log units, respectively.
SPARC chemical reactivity and physical processes models were coupled and extended to calculate hydration equilibrium constants for aldehydes, ketones, quinazoline and substituted quinazolines compounds from molecular structure. The energy differences between the initial (anhydrible) and the final (hydrated) states in the gas phase for a molecule of interest were calculated using SPARC mechanistic perturbation models. These perturbations models quantify the interactions of the appended perturber (P) with the carbonyl reaction center of the aldehydes or ketones or with the imine reaction center of the quinazolines. The perturbations of the reaction center were factored into mechanistic components of electrostatic, resonance and steric effects in the gas phase. The solvation energy (Henrys constant) of the anhydrible and the hydrated states on going from the gas phase to the aqueous phase were calculated using SPARC physical processes models. These physical processes models quantify the intermolecular interactions between the solute and the solvent molecules upon placing a solute molecule in the aqueous phase. The intermolecular interactions are factored into dispersion, induction, dipole-dipole and H-bonding interaction mechanisms. The RMS deviation error was 0.36 pK hydration units for 36 aldehyde and ketone compounds and 0.43 pK hydration units for quinazoline and 31 substituted quinazoline compounds.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
customersupport@researchsolutions.com
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
This site is protected by reCAPTCHA and the Google Privacy Policy and Terms of Service apply.
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.