Impact of dielectric constant of solvent on the formation of transition metal-ammine complexes.

The DFT-level computational investigations into Gibbs free energies (ΔG) demonstrate that as the dielectric constant of the solvent increases, the stabilities of [M(NH3 )n ]2+/3+ (n = 4, 6; M = selected 3d transition metals) complexes decrease. However, there is no observed correlation between the stability of the complex and the solvent donor number. Analysis of the charge transfer and Wiberg bond indices indicates a dative-bond character in all the complexes. The solvent effect assessed through solvation energy is determined by the change in the solvent accessible surface area (SASA) and the change in the charge distribution that occurs during complex formation. It has been observed that the SASA and charge transfer are different in the different coordination numbers, resulting in a variation in the solvent effect on complex stability in different solvents. This ultimately leads to a change between the relative stability of complexes with different coordination numbers while increasing the solvent polarity for a few complexes. Moreover, the findings indicate a direct relationship between ΔΔG (∆Gsolvent -∆Ggas ) and ΔEsolv , which enables the computation of ΔG for the compounds in a particular solvent using only ΔGgas and ΔEsolv . This approach is less computationally expensive.

Polyparameter Linear Free Energy Relationships for Partitioning of Neutral Organic Compounds to Storage Lipids.

Storage lipids are an important compartment in the bioaccumulation of neutral organic compounds. Reliable models for predicting storage lipid-water and storage lipid-air partition coefficients (Kislip/w and Kislip/a), as well as their temperature dependence, are considered useful. Polyparameter linear free energy relationships (PP-LFERs) are accurate, general, and mechanistically clear models for predicting partitioning-related physicochemical quantities. About a decade ago, PP-LFERs were calibrated for Kislip/w at the physiological temperature of 37 °C. However, to date, a comprehensive collection and sufficiently reliable PP-LFERs for Kislip/w and Kislip/a at the most common standard temperature of 25 °C are still lacking. In this study, experimentally based Kislip/w and/or Kislip/a values at 25 °C for 278 compounds were extensively collected or converted from the literature. Subsequently, PP-LFERs were calibrated for Kislip/w and Kislip/a at 25 °C, performing well over 10 orders of magnitude with root-mean-square errors of 0.17-0.21 log units for compounds with reliable descriptors. Furthermore, standard internal energy changes of transfer from water or air to storage lipids for 158 compounds were derived and used to calibrate PP-LFERs for estimating the temperature dependence of Kislip/w or Kislip/a. Additionally, using PP-LFERs, low-density polyethylene was confirmed to be a better storage lipid analogue than silicone and polyoxymethylene in the equilibrium passive sampling of nonpolar and H-bond acceptor polar compounds.

The Influence of Ions on the Electrochemical Stability of Aqueous Electrolytes.

The electrochemical stability window of water is known to vary with the type and concentration of dissolved salts. However, the underlying influence of ions on the thermodynamic stability of aqueous solutions has not been fully understood. Here, we investigated the electrolytic behaviors of aqueous electrolytes as a function of different ions. Our findings indicate that ions with high ionic potentials, i.e., charge density, promote the formation of their respective hydration structures, enhancing electrolytic reactions via an inductive effect, particularly for small cations. Conversely, ions with lower ionic potentials increase the proportion of free water molecules-those not engaged in hydration shells or hydrogen-bonding networks-leading to greater electrolytic stability. Furthermore, we observe that the chemical environment created by bulky ions with lower ionic potentials impedes electrolytic reactions by frustrating the solvation of protons and hydroxide ions, the products of oxygen evolution reaction (OER) and hydrogen evolution reaction (HER), respectively. We found that the solvation of protons plays a more substantial role than that of hydroxide, which explains a greater shift for OER than for HER, a puzzle that cannot be rationalized by the notion of varying O-H bond strengths of water. These insights will help the design of aqueous systems.

The Impact of the Solvent Dielectric Constant on A←NH3 Dative Bond Depends on the Nature of the Lewis Electron-Pair Systems.

The present work aims to determine to what extent the value of the dielectric constant of the solvent can influence the dative bond in Lewis electron pair bonding systems. For this purpose, two different systems, namely H3 B←NH3 and {Zn←(NH3 )}2+ , were studied in selected solvents with significantly different dielectric constants. Based on the results from state-of-the-art computational methods using DFT, constrained DFT, energy decomposition analyses, solvent accessible surface area, and charge transfer calculations, we found that the stability of the neutral H3 B←NH3 system increases with increasing solvent polarity. In contrast, the opposite trend is observed for the positively charged {Zn←(NH3 )}2+ . The observed changes are attributed to different charge redistributions in neutral and charged complexes, which are reflected by a different response to the solvent and are quantified by changes in solvation energies.

Estimating the Gibbs Hydration Energies of Actinium and Trans-Plutonium Actinides.

The use of actinides for medical, scientific and technological purposes has gained momentum in the recent years. This creates a need to understand their interactions with biomolecules, both at the interface and as they become complexed. Calculation of the Gibbs binding energies of the ions to biomolecules, i. e., the Gibbs energy change associated with a transfer of an ion from the water phase to its binding site, could help to understand the actinides' toxicities and to design agents that bind them with high affinities. To this end, there is a need to obtain accurate reference values for actinide hydration, that for most actinides are not available from experiment. In this study, a set of ionic radii is developed that enables future calculations of binding energies for Pu3+ and five actinides with renewed scientific and technological interest: Ac3+ , Am3+ , Cm3+ , Bk3+ and Cf3+ . Reference hydration energies were calculated using quantum chemistry and ion solvation theory and agree well for all ions except Ac3+ , where ion solvation theory seems to underestimate the magnitude of the Gibbs hydration energy. The set of radii and reference energies that are presented here provide means to calculate binding energies for actinides and biomolecules.

Intermolecular Interactions, Solute Descriptors, and Partition Properties of Neutral Per- and Polyfluoroalkyl Substances (PFAS).

The environmental partition properties of perfluoroalkyl and polyfluoroalkyl substances (PFAS) must be understood for their transport and fate analysis. In this study, isothermal gas chromatographic (GC) retention times of 60 neutral PFAS were measured using four columns with different stationary phase polarities, which indicated varying polar interactions exerted by these substances. The GC data were combined with new octanol/water partition coefficient data from this study and existing partition coefficient data from the literature and used to determine the polyparameter linear free energy relationship (PP-LFER) solute descriptors. A complete set of the solute descriptors was obtained for 47 PFAS, demonstrating the characteristic intermolecular interaction properties, such as hydrogen bonding capabilities influenced by the electron-withdrawing perfluoroalkyl group. The partition coefficients between octanol and water, air and water, and octanol and air predicted by the PP-LFER models agreed with those predicted by the quantum chemically based model COSMOtherm, suggesting that both models are highly accurate for neutral PFAS and can fill the current large data gaps in partition property data. A chemical partitioning space plot was generated by using the PP-LFER-predicted partition coefficients, showing the primary importance of the air phase for the environmental distribution of nonpolar and weakly polar PFAS and the increasing significance of organic phases with increasing PFAS polarity.

Strengthening Aqueous Electrolytes without Strengthening Water.

Aqueous electrolytes typically suffer from poor electrochemical stability; however, eutectic aqueous solutions-25 wt.% LiCl and 62 wt.% H3 PO4 -cooled to -78 °C exhibit a significantly widened stability window. Integrated experimental and simulation results reveal that, upon cooling, Li+ ions become less hydrated and pair up with Cl- , ice-like water clusters form, and H⋅⋅⋅Cl- bonding strengthens. Surprisingly, this low-temperature solvation structure does not strengthen water molecules' O-H bond, bucking the conventional wisdom that increasing water's stability requires stiffening the O-H covalent bond. We propose a more general mechanism for water's low temperature inertness in the electrolyte: less favorable solvation of OH- and H+ , the byproducts of hydrogen and oxygen evolution reactions. To showcase this stability, we demonstrate an aqueous Li-ion battery using LiMn2 O4 cathode and CuSe anode with a high energy density of 109 Wh/kg. These results highlight the potential of aqueous batteries for polar and extraterrestrial missions.

Towards a transferable nonelectrostatic model for continuum solvation: The electrostatic and nonelectrostatic energy correction model.

In this work, we introduce an electrostatic and non-electrostatic (ENE) correction to the solvation energy based on the Solvent-Accessible Surface Area (SASA) of the solute and the solvent static dielectric constant. The proposed correction was developed for neutral solutes in non-aqueous solvents, considering three different implicit solvation models based on a Self-Consistent Reaction Field treatment of solute-solvent mutual polarization using an Apparent Surface Charge formalism, namely the Integral Equation Formalism of the Polarizable Continuum Model using a continuous surface charge scheme (PCM), the Solvation Model based on solute electron density (SMD), and the generalized Finite-Difference Poisson-Boltzmann (FDPB) model. The proposed correction was parametrized on a diverse training set of 4980 solvation data from the Solv@tum database of experimental solvation energies, and validated on the non-aqueous subset of the MNSOL database comprising 2140 solvation energies. The performances of the proposed ENE models with minimal and extended parameters formulations have been analyzed and the latter variant has been further compared to the widely used Cavity, Dispersion, and Solvent structural effects (CDS) non-electrostatic model originally developed for the SMx family of implicit solvation models. Overall, a very good agreement between the computed solvation energies with the ENE correction and the reference experimental data has been found on both the training and test sets for all continuum solvation models considered. Furthermore, results for the ENE correction are on par with the reference CDS non-electrostatic model for both SMD and FDPB electrostatics, but with the advantage of using a lower number of parameters and thus an improved transferability between different electrostatics treatments.

Evaluation of Two Cosolvency Models to Predict Solute Partitioning between Polymers (LDPE) and Water - Ethanol Simulating Solvent Mixtures.

PURPOSE: Binary water - ethanol mixtures, by mimicking a clinically relevant medium's polarity-driven extraction strength, facilitate experimental modeling of patient exposure to chemicals which can potentially leach from a plastic material for pharmaceutical applications. Estimates of patient exposure could consequently benefit from a quantitative concept for tailoring the extraction strength of the simulating solvent mixture towards the one of the clinically relevant medium. METHODS: The hypothetical partition coefficient based upon the differential solubility between water-ethanol mixtures and water, [Formula: see text], has been calculated by the log-linear model from Yalkowsky and coworkers and a cosolvency model based on Abraham-type linear solvation energy relationships (LSERs). Then, by applying a thermodynamic cycle using the partition coefficient LDPE/water, [Formula: see text], partitioning between LDPE and the ethanol in water mixture was calculated and experimentally verified for a wide array of chemically diverse solutes. RESULTS: The partition coefficients between LDPE and volume fractions of 0.1, 0.2, 0.35 and 0.5 of ethanol in water calculated by this approach correlated well with experimentally obtained values. The LSER based model was found slightly superior over the log-linear cosolvency model. CONCLUSIONS: Solubilization strength projection by means of cosolvency models in combination with LSER predicted partition coefficients LDPE/water enable the tailored preparation of water-ethanol simulating solvent mixtures when input parameters from the clinically relevant medium are available. This approach can increase the reliability of patient exposure estimations and avoid overly complex extraction profiles, thus minimizing time and resources for chemical safety risk assessments on plastic materials used in pharmaceutical applications.

Applicability Domain of Polyparameter Linear Free Energy Relationship Models Evaluated by Leverage and Prediction Interval Calculation.

Polyparameter linear free energy relationships (PP-LFERs) are accurate and robust models employed to predict equilibrium partition coefficients (K) of organic chemicals. The accuracy of predictions by a PP-LFER depends on the composition of the respective calibration data set. Generally, extrapolation outside the domain defined by the calibration data is likely to be less accurate than interpolation. In this study, the applicability domain (AD) of PP-LFERs was systematically evaluated by calculating the leverage (h) and prediction interval (PI). Repeated simulations with experimental data showed that the root mean squared error of predictions increased with h. However, the analysis also showed that PP-LFERs calibrated with a large number (e.g., 100) of training data were highly robust against extrapolation error. For such PP-LFERs, the common definition of extrapolation (h > 3 hmean, where hmean is the mean h of all training compounds) may be excessively strict. Alternatively, the PI is proposed as a metric to define the AD of PP-LFERs, as it provides a concrete estimate of the error range that agrees well with the observed errors, even for extreme extrapolations. Additionally, published PP-LFERs were evaluated in terms of their AD using the new concept of AD probes, which indicated the varying predictive performance of PP-LFERs in the existing literature for environmentally relevant compounds.

Predicting Leachables Solubilization in Polysorbate 80 Solutions by a Linear Solvation Energy Relationship (LSER).

PURPOSE: A linear solvation energy relationship (LSER) was developed to predict the partitioning of neutral chemicals from polysorbate 80 (PS 80) micelles to water. Predicted partition coefficients were converted to a concentration dependent solubilization strength of aqueous PS 80 solutions. This solubilization strength represents a key parameter to project equilibrium levels of leaching from pharmaceutical plastic materials. METHODS: To construct the LSER model equation, partition coefficients between PS 80 micelles and water were measured via a reference phase method or collected from the literature. Multiple linear regression of partition coefficients against five publicly available solute parameters was used to obtain the LSER system parameters. RESULTS: 112 chemically diverse compounds were incorporated for LSER model regression. The model equation shows a very good fit (R2 = 0.969, SD = 0.219) for the entire dataset. The accuracy of the multi-parameter LSER model was proven to be substantially better in comparison to a single-parameter log-linear model based on the octanol-water partition coefficient. CONCLUSION: PS 80 solubilization strength in water can expediently and accurately be calculated for neutral organic compounds with the proposed LSER model. LSER system parameters provide insightful chemical information with respect to solubilization in aqueous solutions of PS 80.

Standard and transformed values of gibbs energy formation for some radicals and ions involved in biochemical reactions.

The standard and transformed values of the Gibbs formation function of a number of radicals and ions are calculated: H, OH, HO2,O, SH, NH2, CH3, H+,O-,OH-HO2-, O2-,SH-NH2-, Q-, NAD*, FMN-, FAD-. These data can be used in consideration of the thermodynamics of biochemical reactions involving free radicals.

Preparation and comparison of three zwitterionic stationary phases for hydrophilic interaction liquid chromatography.

Surface-bonded zwitterionic stationary phases have shown highlighted performances in separation of polar and hydrophilic compounds under hydrophilic interaction chromatography mode. So, it would be helpful to evaluate the characteristics of zwitterionic stationary phases with different arranged charged groups. The present work involved the preparation and comparison of three zwitterionic stationary phases. An imidazolium ionic liquid was designed and synthesized, and the cationic and anionic moieties respectively possessed positively charged imidazolium ring and negatively charged sulfonic groups. Then, the prepared ionic liquid, phosphorylcholine and an imidazolium-based zwitterionic selector were bonded on the surface of silica to obtain three zwitterionic stationary phases. The selectivity properties were characterized and compared through the relative retention of selected solute pairs, and different kinds of hydrophilic solutes mixtures were used to evaluate the chromatographic performances. Moreover, the zwitterionic stationary phases were further characterized by the modified linear solvation energy relationship model to probe the multiple interactions. All the results indicated that the types and arrangement of charged groups in zwitterionic stationary phases mainly affect the retention and separation of ionic or ionizable compounds, and for interaction characteristics the contribution from n and π electrons and electrostatic interactions displayed certain differences.

Fragment-based quantum mechanical calculation of protein-protein binding affinities.

The electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method has been successfully utilized for efficient linear-scaling quantum mechanical (QM) calculation of protein energies. In this work, we applied the EE-GMFCC method for calculation of binding affinity of Endonuclease colicin-immunity protein complex. The binding free energy changes between the wild-type and mutants of the complex calculated by EE-GMFCC are in good agreement with experimental results. The correlation coefficient (R) between the predicted binding energy changes and experimental values is 0.906 at the B3LYP/6-31G*-D level, based on the snapshot whose binding affinity is closest to the average result from the molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) calculation. The inclusion of the QM effects is important for accurate prediction of protein-protein binding affinities. Moreover, the self-consistent calculation of PB solvation energy is required for accurate calculations of protein-protein binding free energies. This study demonstrates that the EE-GMFCC method is capable of providing reliable prediction of relative binding affinities for protein-protein complexes. © 2018 Wiley Periodicals, Inc.

Determination of lipid bilayer affinities and solvation characteristics by electrokinetic chromatography using polymer-bound lipid bilayer nanodiscs.

Styrene-maleic acid polymer-bound lipid bilayer nanodiscs have been investigated and characterized by electrokinetic chromatography. Linear solvation energy relationship analysis was employed to characterize the changes in solvation environment of nanodiscs of varied belt to lipid ratio, belt polymer chemistry and molecular weight, and lipid composition. Increases in the lipid to belt polymer ratio resulted in smaller, more cohesive nanodiscs with greater electrophoretic mobility. Nanodisc structures with belt polymers of different chemistry and molecular weight were compared and showed only minor changes in solvent characteristics and selectivity consistent with changes in structure of the lipid bilayer. Seven phospholipid and sphingomyelin nanodiscs of different lipid composition were characterized. Changes in lipid head group structure had a significant effect on bilayer-solute interactions. In most cases, changes in alkyl tail structure had no discernible effect on solvation environment aside from those explained by changes in the gel-liquid transition temperature. Comparison to vesicles of similar lipid composition show only minor differences in solvation environment, likely due to differences in lipid composition and bilayer curvature. Together these results provide evidence that the dominant solute-nanodisc interactions are with the lipid bilayer and that head group chemistry has a greater impact on bilayer-solute interactions than alkyl tail or belt polymer structure. Nanodisc electrokinetic chromatography is demonstrated to allow characterization of solute interactions with lipid bilayers of varied composition.

Comparison of a polymeric pseudostationary phase in EKC with ODS stationary phase in RP-HPLC.

Poly(stearyl methacrylate-co-methacrylic acid) (P(SMA-co-MAA)) was induced as pseudostationary phase (PSP) in electrokinetic chromatography (EKC). The n-octadecyl groups in SMA were the same as that in octadecylsilane (ODS) C18 column. Thus, the present work focused on the comparison of selectivity between polymeric PSP and ODS stationary phase (SP), and the effect of organic modifiers on the selectivity of polymeric PSP and ODS SP. 1-butanol could directly interacted with PSP as a Class I modifier, and improved both of the methylene selectivity and polar group selectivity. When the analysis times were similar, the polymeric PSP exhibited better methylene selectivity and polar group selectivity. Although the hydrophobic groups were similar, the substituted benzenes elution order was different between polymeric PSP and ODS SP. Linear solvation energy relationships (LSER) model analysis found that polymeric PSP and ODS SP exhibited two same key factors in selectivity: hydrophobic interaction and hydrogen bonding acidity. But polymeric PSP exhibited relatively strong n- and π-electrons interaction to the analytes.

Protonation of polyaniline-coated silica stationary phase affects the retention behavior of neutral hydrophobic solutes in reversed-phase capillary liquid chromatography.

Because of its high conductivity when acid doped, polyaniline is known as a synthetic metal and is used in a wide range of applications, such as supercapacitors, biosensors, electrochromic devices, or solar and fuel cells. Emeraldine is the partly oxidized, stable form of polyaniline, consisting of alternating diaminobenzenoid and iminoquinoid segments. When acidified, the nitrogen atoms of emeraldine become protonated. Due to electrostatic repulsion between positive charges, the polarity and morphology of emeraldine chains presumably change; however, the protonation effects on emeraldine have not yet been clarified. Thus, we investigated these changes by reversed-phase capillary liquid chromatography using a linear solvation energy relationship approach to assess differences in dominant retention interactions under a significantly varied mobile phase pH. We observed that hydrophobicity dominates the intermolecular interactions under both acidic and alkaline eluent conditions, albeit to different extents. Therefore, by tuning the mobile phase pH, we can even modulate the retention of neutral hydrophobic solutes, such as aromatic hydrocarbons, because the pH-dependent charge and structure of polymer chains of the emeraldine-coated silica stationary phase show a mixed-mode separation mechanism.

Development and evaluation of MTLSER and QSAR models for predicting polyethylene-water partition coefficients.

Current study was aimed to make further improvements in measuring low density polyethylene (LDPE) -water partition coefficient (KPE-w) for organic chemicals. Modified theoretical linear solvation energy relationship (MTLSER) model and quantitative structure activity relationship (QSAR) model were developed for predicting KPE-w values from chemical descriptors. With the MTLSER model, α (average molecular polarizability), µ (dipole moment) and q- (net charge of the most negative atoms) as significant variables were screened. With the QSAR model, main control factors of KPE-w values, such as CrippenLogP (Crippen octanol-water partition coefficient), CIC0 (neighborhood symmetry of 0-order) and GATS2p (Geary autocorrelation-lag2/weighted by polarizabilities) were studied. As per our best knowledge, this is the first attempt to predict polymer-water partition coefficient using the MTLSER model. Statistical parameters, correlation coefficient (R2) and cross-validation coefficients (Q2) were ranging from 0.811 to 0.951 and 0.761 to 0.949, respectively, which indicated that the models appropriately fit the results, and also showed robustness and predictive capacity. Mechanism interpretation suggested that the main factors governing the partition process between LDPE and water were the molecular polarizability and hydrophobicity. The results of this study provide an excellent tool for predicting log KPE-w values of most common hydrophobic organic compounds, within the applicability domains to reduce experimental cost and time for innovation.

Phospholipid bilayer affinities and solvation characteristics by electrokinetic chromatography with a nanodisc pseudostationary phase.

Phospholipid bilayer nanodiscs composed of 1,2-dimyristoyl-sn-glycero-3-phosphocholine and synthetic maleic acid-styrene copolymer belts have been introduced as a pseudostationary phase (PSP) in electrokinetic chromatography and demonstrated good performance. The nanodiscs provide a suitable migration range and high theoretical plate counts. Using this nanodisc pseudostationary phase, the affinity of the bilayer structure for probe solutes was determined and characterized. Good correlation is observed between retention factors and octanol water partition coefficients for particular categories of solutes, but the general correlation is weak primarily because the nanodiscs show stronger affinity than octanol for hydrogen bond donors. This suggests that a more appropriate application of this technology is to measure and characterize interactions between solutes and lipid bilayers directly. Linear solvation energy relationship analysis of the nanodisc-solute interactions in this study demonstrates that the nanodiscs provide a solvation environment with low cohesivity and weak hydrogen bond donating ability, and provide relatively strong hydrogen bond acceptor strength.

Characterization of polyaniline-coated stationary phases by using the linear solvation energy relationship in the hydrophilic interaction liquid chromatography mode using capillary liquid chromatography.

A polyaniline coating was used to modify the surface of bare silica gel and octadecyl silica stationary phases to characterize the properties of altered materials. It was assumed that the mixed-mode retention was established on the basis of the polyaniline chemical structure and its combination with the original sorbents. Polyaniline was deposited onto the original surfaces during the chemical polymerization of aniline hydrochloride. The prepared materials were slurry packed into capillary columns and systematic chromatographic characterization was performed using the linear solvation energy relationship, also employing descriptors that allow inclusion of ionic interactions in the proposed retention mechanism. The retention times of 80 solutes with various chemical structures were measured in the hydrophilic interaction liquid chromatography mode. The obtained results demonstrated the significant contribution of the polyaniline coating to the retention mechanism under the given conditions; the assumed mixed-mode retention was confirmed. The dominant retention interaction for both modified stationary phases was based on the protonation of nitrogen atoms in the polyaniline structure, leading to suitable retention and selectivity for the hydrophilic analytes, especially anionic and zwitterionic species. Thus, especially, the polyaniline-coated bare silica gel sorbent seems to be promising for potential applications related to the separation of polar compounds.