Differential Hydration of Phenol and Phenoxy Radical and the Energetics of the Phenol O−H Bond in Solution

Unraveling the acid-base characterization and solvent effects on the structural and electronic properties of a bis-bidentate bridging ligand

Understanding the interactions and the solvent effects on the distribution of several species in equilibrium and how it can influence the ¹H-NMR properties, spectroscopy (UV-vis absorption), and the acid-base equilibria can be especially challenging. This is the case of a bis-bidentate bridging ligand bis(2-pyridyl)-benzo-bis(imidazole), where the two pyridyl and four imidazolyl nitrogen atoms can be protonated in different ways, depending on the solvent, generating many isomeric/tautomeric species. Herein, we report a combined theoretical-experimental approach based on a sequential quantum mechanics/molecular mechanics procedure that was successfully applied to describe in detail the acid-base characterization and its effects on the electronic properties of such a molecule in solution. The calculated free-energies allowed the identification of the main species present in solution as a function of the solvent polarity, and its effects on the magnetic shielding of protons (¹H-NMR chemical shifts), the UV-vis absorption spectra, and the acid-base equilibrium constants (pK_as) in aqueous solution. Three acid-base equilibrium constants were experimentally/theoretically determined (pK_a1 = 1.3/1.2, pK_a2 = 2.1/2.2 and pK_a5 = 10.1/11.3) involving mono-deprotonated and mono-protonated cis and trans species. Interestingly, other processes with pK_a3 = 3.7 and pK_a4 = 6.0 were also experimentally determined and assigned to the protonation and deprotonation of dimeric species. The dimerization of the most stable neutral species was investigated by Monte Carlo simulations and its electronic effects were considered for the elucidation of the UV-vis absorption bands, revealing transitions mainly with the charge-transfer characteristic and involving both the monomeric species and the dimeric species. The good matching of the theoretical and experimental results provides an atomistic insight into the solvent effects on the electronic properties of this bis-bidentate bridging ligand.

A QM/MM study of the conformation stability and electronic structure of the photochromic switches derivatives of DHA/VHF in acetonitrile solution

We present a detailed theoretical study of the electronic absorption spectra and thermochemistry of molecular photoswitches composed of one and two photochromic units of dihydroazulene (DHA)/vinylheptafulvene (VHF) molecules. Six different isomers are considered depending on the ring opening/closure forms of the DHA units. The solvent effect of acetonitrile is investigated using a sequential Molecular Mechanics/Quantum Mechanics approach. The thermochemical investigations of these photochromic molecules were performed using the Free Energy Perturbation method, and the simulations were performed using Configurational Bias Monte Carlo. We show that to open the 5-member ring of the DHA, there is no significant gain in thermal release of energy for the back reaction when a unit or two DHA units are considered. Overall, we found agreement between the solvation free energy based on Monte Carlo simulations and the continuum solvent model. However, the cavitation term in the continuum model is shown to be a source of disagreement when the non-electrostatic terms are compared. The electronic absorption spectra are calculated using TDDFT CAM-B3LYP/cc-pVDZ. Agreement with experiment is obtained within 0.1 eV, considering statistically uncorrelated configurations from the simulations. Inhomogeneous broadening is also considered and found to be well described in all cases.

Energetics of Radical Formation in Eumelanin Building Blocks: Implications for Understanding Photoprotection Mechanisms in Eumelanin

The supramolecular structure of melanin pigments is characterized by a high concentration of radical species. Therefore, the energetics of the radical formation in melanin building blocks is key for understanding the structure and the electronic properties of the pigments at the molecular level. Nevertheless, the radical energetics of even the simplest melanin building blocks are largely unknown. In order to address this fundamental issue, the bond dissociation enthalpies (BDEs) for the melanin monomers 5,6-dihydroxy-1H-indole-2-carboxylic acid (DHICA), 1H-indole-5,6-diol (DHI), and 1H-indole-5,6-dione (IQ) were determined through high-accuracy ab initio quantum chemistry methods. Our results provide strong evidence of the importance on BDEs for explaining the experimentally observed dependence of the antioxidant properties of eumelanin pigments on the DHICA/DHI ratio, and the role that these two species play on the photoprotection mechanism.

Photoionization access to cyclodextrin-encapsulated resveratrol phenoxy radicals and their repair by ascorbate across the phase boundary

Photoionization with UV-A generates phenoxy radical ResO˙ of the famous red-wine antioxidant resveratrol inside cyclodextrin cavities without the participation of auxiliary reactants. This permits the isolated exploration of the ResO˙ repair by ascorbate at the water–cyclodextrin interface.

A Multiscale Treatment of Angeli's Salt Decomposition

Sodium trioxodinitrate's (Na2N2O3, Angeli's salt) unique cardiovascular effects have been associated with its ability to yield HNO upon dissociation under physiological conditions. Due to its potential applications in new therapies for heart failure, the dissociation of Angeli's salt has recently received increased attention. The decomposition mechanism has been previously studied by quantum mechanical methods using a continuum approximation (PCM) for the solvent effects. In this work we use our recently developed interface of the Amber and Gaussian packages via the PUPIL package to study Angeli's salt dissociation in a hybrid QM/MM scheme where the water solvent molecules are treated explicitly with classical mechanics while the solute is treated with full quantum mechanics (UB3LYP/6-31+G(d) and UMP2/6-31+G(d)) level. Multiple steered molecular dynamics was used with the Jarzynski relationship to extract the free energy profile for the process. We obtain 4.8 kcal mol(-1) and 6.4 kcal mol(-1) free energy barriers for the N-N bond breaking for UB3LYP and UMP2, respectively. The geometries and Mulliken charges for reactant, transition state, and products have been characterized through a number of hybrid QM/MM molecular dynamics runs with the N-N distance restrained to representative values of each species. The results highlight the role of individual solvent molecules for the reaction energetics and provide a comparison point against implicit solvation methods.

Kinetic and equilibrium studies of phenol adsorption by natural and modified forms of the clinoptilolite

The contribution presents results of investigation of phenol adsorption by the raw, HDTMA- and NaOH-modified clinoptilolite. The experimental data on the sorption kinetic process are modelled using the particle and the film diffusion models with calculation of a rate constant of intraparticle mass transfer. Two main stages may be divided in the kinetic sorption process: (1) initial rapid sorption (amount of phenol sorbed 85-90%) extending over the first 60 min and (2) stage of slow approach to equilibrium covering about 8h. Equilibrium isotherms for phenol adsorption are modelled using the Langmuir, Freundlich and Dubinin-Radushkevich equations. The Freundlich isotherm describes in a greater degree phenol sorption from aqueous solutions of low phenol concentrations, while the Langmuir isotherm fits better to sorption at high initial concentrations.

Theoretical insights, in the liquid phase, into the antioxidant mechanism-related parameters in the 2-monosubstituted phenols

The paper describes a DFT/B3LYP study, in the liquid phase, [using the PCM continuum model] on the O-H bond dissociation enthalpy (BDE) and ionization energy (IE) parameter values of the 2-monosubstituted phenols (2-X-ArOH), related to the H-atom transfer (HAT) and single-electron transfer (SET) mechanisms. The solvent and substituent effects on the conformers, the BDEs, and the IEs were studied using four electron-donating (EDG) and five electron-withdrawing (EWG) groups, in seven different solvents. In both the EDG- and/or EWG-substituted species of the parent compounds, radicals, and/or cation radicals, the most stable conformer is varied, depending on the medium and the substitution. The EWG-substituents increase IEs, resulting in a weaker antioxidant activity than the EDG ones; the effect appears stronger on the IEs than on BDEs. However, although the liquid-phase IEs, which are related to solution-phase oxidation potentials, decrease with the polarity and/or the hydrogen-bonding ability of the solvent, the opposite holds true for the BDEs, exhibiting a weaker effect. The gas-phase-calculated IE for benzene is among the most accurate ones in the field, compared to the experiment, that for phenol being the most accurate. In addition, calculated IEs for the 2-X-ArOH are in close agreement with the very few existing experimental ones. It is shown that the oxidation potentials are (a) highly correlated with the gas-phase ones, and (b) strongly solvent dependent. The stabilization/destabilization of the cation radical (SPC) contribution, in all media, is the decisive factor in the DeltaIE calculation. The reasonable correlations found between the DeltaBDE and DeltaIE could account well for the assumption of the simultaneous action of both mechanisms in the 2-X-ArOH, in both the gas and the liquid phase. It seems, however, that the presence of a particular solvent by itself is not sufficient enough for the HAT to SET transition. The involvement of specific ED and/or EW groups in the 2-X-ArOH seems also necessary. It appears that our theoretical approach is not only generally applicable to the set of substituents important to antioxidant activity but also useful in (a) the rational design of phenolic antioxidants and (b) affording accurate BDE and IE parameter values related to both possible antioxidant mechanisms.

Computational Study on the Properties and Structure of Methyl Lactate

A theoretical study on the properties and molecular level structure of the very important green solvent methyl lactate is carried out in the gas phase and methanol and water solutions, with the solvent treated both explicitly and as a continuum. Torsional barriers giving rise to different conformers by rotation of the hydroxyl and methyl groups were analyzed using density functional theory (DFT) to establish the most stable conformer both in gas phase and solution. DFT computations on lactate dimers were also done to study short-range features, and the effect of the surrounding solvent on intra- and intermolecular hydrogen bonding was analyzed according to the polarizable continuum model approach. We have also studied lactate/water and lactate/methanol small clusters together with the corresponding binding energies. Moreover, classical molecular dynamics simulations (MD) were carried out to study medium- and large-range effects at lower computational cost. MD simulations at different pressure and temperature conditions on pure lactate were carried out, and mixtures with water and methanol of different compositions were also studied. Structural information, analyzed through the radial distribution functions, together with dynamic aspects of pure and mixed fluids were considered. The intramolecular hydrogen bonding ability of methyl lactate together with the possibility of homo- and hetero-intermolecular association determines the behavior of this molecule in pure fluids or in mixed.

Enthalpy of Solvation Correlations for Gaseous Solutes Dissolved in Water and in 1-Octanol Based on the Abraham Model

Data have been assembled on the enthalpies of solvation of 373 compounds in water and 138 compounds in 1-octanol. It is shown that an Abraham solvation equation with five descriptors can be used to correlate the experimental solvation enthalpies to within standard deviations of 3.68 kJ/mol (water) and 2.66 kJ/mol (1-octanol). The derived correlations provide very accurate mathematical descriptions of the observed enthalpies of solvation, which in the case of water span a range of 150 kJ/mol. Division of the experimental values into a training set and a test set shows that there is no bias in predictions and that the predictive capability of the correlations is better than 4 kJ/mol.

Supramolecular Structure of N,N-Bis(2-hydroxybenzyl)alkylamine: Flexible Molecular Assembly Framework for Host without Guest and Host with Guest

N,N-Bis(2-hydroxybenzyl)alkylamine derivatives form a cage-like assembly consisting of two molecules via inter- and intramolecular hydrogen bonds. The derivatives exhibit themselves as host to accept copper-ion guests under the double-oxygen-bridged dimeric system. Quantum chemical calculation suggested that the host-guest interaction is based on a charge-transfer coordination. Comparison of the crystal structures before and after complexation clarifies a rare example of a host-guest compound where the hosts maintain their cage framework through the change of hydrogen bonds to coordination bonds.

Reaction Mechanism and Tautomeric Equilibrium of 2-Mercaptopyrimidine in the Gas Phase and in Aqueous Solution: A Combined Monte Carlo and Quantum Mechanics Study

A combined Monte Carlo and quantum mechanical study was carried out to analyze the tautomeric equilibrium of 2-mercaptopyrimidine in the gas phase and in aqueous solution. Second- and fourth-order Møller-Plesset perturbation theory calculations indicate that in the gas phase thiol (Pym-SH) is more stable than the thione (Pym-NH) by ca. 8 kcal/mol. In aqueous solution, thermodynamic perturbation theory implemented on a Monte Carlo NpT simulation indicates that both the differential enthalpy and Gibbs free energy favor the thione form. The calculated differential enthalpy is DeltaH(SH)(-->)(NH)(solv) = -1.7 kcal/mol and the differential Gibbs free energy is DeltaG(SH)(-->)(NH)(solv) = -1.9 kcal/mol. Analysis is made of the contribution of the solute-solvent hydrogen bonds and it is noted that the SH group in the thiol and NH group in the thione tautomers act exclusively as a hydrogen bond donor in aqueous solution. The proton transfer reaction between the tautomeric forms was also investigated in the gas phase and in aqueous solution. Two distinct mechanisms were considered: a direct intramolecular transfer and a water-assisted mechanism. In the gas phase, the intramolecular transfer leads to a large energy barrier of 34.4 kcal/mol, passing through a three-center transition state. The proton transfer with the assistance of one water molecule decreases the energy barrier to 17.2 kcal/mol. In solution, these calculated activation barriers are, respectively, 32.0 and 14.8 kcal/mol. The solvent effect is found to be sizable but it is considerably more important as a participant in the water-assisted mechanism than the solvent field of the solute-solvent interaction. Finally, the calculated total Gibbs free energy is used to estimate the equilibrium constant.