Role of Deprotonation Free Energies in pKa Prediction and Molecule Ranking

Theoretical insight of reactive oxygen species scavenging mechanism in lignin waste depolymerization products

Apart from natural products and synthesis, phenolic compounds can be produced from the depolymerization of lignin, a major waste in biofuel and paper production. This process yields a plethora of aryl propanoid phenolic derivatives with broad biological activities, especially antioxidant properties. Due to its versatility, our study focuses on investigating the antioxidant mechanisms of several phenolic compounds obtained from renewable and abundant resources, namely, syringol (Hs), 4-allylsyringol (HAs), 4-propenylsyringol (HPns), and 4-propylsyringol (HPs). Employing the density functional theory (DFT) approach in conjunction with the QM-ORSA protocol, we aim to explore the reactivity of these compounds in neutralizing hydroperoxyl radicals in physiological and non-polar media. Kinetic and thermodynamic parameter calculations on the antioxidant activity of these compounds were also included in this study. Additionally, our research utilizes the activation strain model (ASM) for the first time to explain the reactivity of the HT and RAF mechanisms in the peroxyl radical scavenging process. It is predicted that HPs has the best rate constant in both media (1.13 × 108 M-1 s-1 and 1.75 × 108 M-1 s-1, respectively). Through ASM analysis, it is observed that the increase in the interaction energy due to the formation of intermolecular hydrogen bonds during the reaction is an important feature for accelerating the hydrogen transfer process. Furthermore, by examining the physicochemical and toxicity parameters, only Hs is not suitable for further investigation as a therapeutic agent because of potential toxicity and mutagenicity. However, overall, all compounds are considered potent HOO˙ scavengers in lipid-rich environments compared to previously studied antioxidants.

Reliable and accurate prediction of basic pK[Formula: see text] values in nitrogen compounds: the pK[Formula: see text] shift in supramolecular systems as a case study

This article presents a quantitative structure-activity relationship (QSAR) approach for predicting the acid dissociation constant (pK[Formula: see text]) of nitrogenous compounds, including those within supramolecular complexes based on cucurbiturils. The model combines low-cost quantum mechanical calculations with QSAR methodology and linear regressions to achieve accurate predictions for a broad range of nitrogen-containing compounds. The model was developed using a diverse dataset of 130 nitrogenous compounds and exhibits excellent predictive performance, with a high coefficient of determination (R[Formula: see text]) of 0.9905, low standard error (s) of 0.3066, and high Fisher statistic (F) of 2142. The model outperforms existing methods, such as Chemaxon software and previous studies, in terms of accuracy and its ability to handle heterogeneous datasets. External validation on pharmaceutical ingredients, dyes, and supramolecular complexes based on cucurbiturils confirms the reliability of the model. To enhance usability, a script-like tool has been developed, providing a streamlined process for users to access the model. This study represents a significant advancement in pK[Formula: see text] prediction, offering valuable insights for drug design and supramolecular system optimization.

Ortho Effects of Tricarboxylate Linkers in Regulating Topologies of Rare-Earth Metal-Organic Frameworks

A linker design strategy is developed to attain novel polynuclear rare-earth (RE) metal-organic frameworks (MOFs) with unprecedented topologies. We uncover the critical role of ortho-functionalized tricarboxylate ligands in directing the construction of highly connected RE MOFs. The acidity and conformation of the tricarboxylate linkers were altered by substituting with diverse functional groups at the ortho position of the carboxyl groups. For instance, the acidity difference between carboxylate moieties resulted in forming three hexanuclear RE MOFs with novel (3,3,3,10,10)-c wxl, (3,12)-c gmx, and (3,3,3,12)-c joe topologies, respectively. In addition, when a bulky methyl group was introduced, the incompatibility between the net topology and ligand conformation guided the co-appearance of hexanuclear and tetranuclear clusters, generating a novel 3-periodic MOF with a (3,3,8,10)-c kyw net. Interestingly, a fluoro-functionalized linker prompted the formation of two unusual trinuclear clusters and produced a MOF with a fascinating (3,8,10)-c lfg topology, which could be gradually replaced by a more stable tetranuclear MOF with a new (3,12)-c lee topology with extended reaction time. This work enriches the polynuclear clusters library of RE MOFs and unveils new opportunities to construct MOFs with unprecedented structural complexity and vast application potential.

Solvation entropy, enthalpy and free energy prediction using a multi-task deep learning functional in 1D-RISM

Simultaneous calculation of entropies, enthalpies and free energies has been a long-standing challenge in computational chemistry, partly because of the difficulty in obtaining estimates of all three properties from a single consistent simulation methodology. This has been particularly true for methods from the Integral Equation Theory of Molecular Liquids such as the Reference Interaction Site Model which have traditionally given large errors in solvation thermodynamics. Recently, we presented pyRISM-CNN, a combination of the 1 Dimensional Reference Interaction Site Model (1D-RISM) solver, pyRISM, with a deep learning based free energy functional, as a method of predicting solvation free energy (SFE). With this approach, a 40-fold improvement in prediction accuracy was delivered for a multi-solvent, multi-temperature dataset when compared to the standard 1D-RISM theory [Fowles et al., Digital Discovery, 2023, 2, 177-188]. Here, we report three further developments to the pyRISM-CNN methodology. Firstly, solvation free energies have been introduced for organic molecular ions in methanol or water solvent systems at 298 K, with errors below 4 kcal mol^-1 obtained without the need for corrections or additional descriptors. Secondly, the number of solvents in the training data has been expanded from carbon tetrachloride, water and chloroform to now also include methanol. For neutral solutes, prediction errors nearing or below 1 kcal mol^-1 are obtained for each organic solvent system at 298 K and water solvent systems at 273-373 K. Lastly, pyRISM-CNN was successfully applied to the simultaneous prediction of solvation enthalpy, entropy and free energy through a multi-task learning approach, with errors of 1.04, 0.98 and 0.47 kcal mol^-1, respectively, for water solvent systems at 298 K.

Synthesis of graphene oxide grafted by diazanyl groups and its application in recovery of lead from lead-acid wastewater

Heavy metal ion (HMI) in wastewater is a kind of resource that is wrongly placed. Recovery of heavy metal from lead-acid wastewater desires efficient and reusable functional materials. In this paper, graphene oxide-like with diazanyl groups (GOLA) was synthesized by a nucleophilic substitution reaction of graphene oxide-like with hydroxyl groups (GOLH) with diazane. GOLA exhibited good stability and recyclability in wastewater. The maximal adsorption capacity (q_max) values of GOLA for Pb(II), Cd(II), Cu(II), Ni(II), and Cr(III) ions were 505.80, 401.99, 83.48, 82.29, and 147.77 mg/g, respectively. The equilibrium time of GOLA adsorbing HMIs was 20 min. GOLA was employed to recover lead ions from lead-acid wastewater to give Pb(OH)₂ and reusable water. Therefore, this paper provides a useful method of recycling lead from lead-acid wastewater.

Uncovering Differences in Hydration Free Energies and Structures for Model Compound Mimics of Charged Side Chains of Amino Acids

Free energies of hydration are of fundamental interest for modeling and understanding conformational and phase equilibria of macromolecular solutes in aqueous phases. Of particular relevance to systems such as intrinsically disordered proteins are the free energies of hydration and hydration structures of model compounds that mimic charged side chains of Arg, Lys, Asp, and Glu. Here, we deploy a Thermodynamic Cycle-based Proton Dissociation (TCPD) approach in conjunction with data from direct measurements to obtain estimates for the free energies of hydration for model compounds that mimic the side chains of Arg+, Lys+, Asp-, and Glu-. Irrespective of the choice made for the hydration free energy of the proton, the TCPD approach reveals clear trends regarding the free energies of hydration for Arg+, Lys+, Asp-, and Glu-. These trends include asymmetries between the hydration free energies of acidic (Asp- and Glu-) and basic (Arg+ and Lys+) residues. Further, the TCPD analysis, which relies on a combination of experimental data, shows that the free energy of hydration of Arg+ is less favorable than that of Lys+. We sought a physical explanation for the TCPD-derived trends in free energies of hydration. To this end, we performed temperature-dependent calculations of free energies of hydration and analyzed hydration structures from simulations that use the polarizable Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field and water model. At 298 K, the AMOEBA model generates estimates of free energies of hydration that are consistent with TCPD values with a free energy of hydration for the proton of ca. -259 kcal/mol. Analysis of temperature-dependent simulations leads to a structural explanation for the observed differences in free energies of hydration of ionizable residues and reveals that the heat capacity of hydration is positive for Arg+ and Lys+ and negative for Asp- and Glu-.

Polyelectrolyte pKa from experiment and molecular dynamics simulation

The pK_a of a polyelectrolyte has been determined experimentally by potentiometric titration and computed using Molecular Dynamics (MD) constant pH (CpH) methodology, which allows the pK_a of each titratable site along the polymer backbone.

Interacting quantum fragments-rooted preorganized-interacting fragments attributed relative molecular stability of the Be(II) complexes of nitrilotriacetic acid and nitrilotri-3-propionic acid

A method designed to investigate, on a fundamental level, the origin of relative stability of molecular systems using Be(II) complexes with nitrilotriacetic acid (NTA) and nitrilotri-3-propionic acid (NTPA) is described. It makes use of the primary and molecular fragment energy terms as defined in the IQA/F (Interacting Quantum Atoms/Fragments) framework. An extensive classical-type investigation, focused on single descriptors (bond length, density at critical point, the size of metal ion or coordination ring, interaction energy between Be(II) and a donor atom, etc.) showed that it is not possible to explain the experimental trend. The proposed methodology is fundamentally different in that it accounts for the total energy contributions coming from all atoms of selected molecular fragments, and monitors changes in defined energy terms (e.g., fragment deformation, inter- and intra-fragment interaction) on complex formation. By decomposing combined energy terms we identified the origin of relative stability of Be(II) (NTA) and Be(II) (NTPA) complexes. We found that the sum of coordination bonds' strength, as measured by interaction energies between Be(II) ion and donor atoms, favours Be(II) (NTA) but the binding energy of Be(II) ion to the entire ligand correlates well with experimental trend. Surprisingly, the origin of Be(II) (NTPA) being more stable is due to less severe repulsive interactions with the backbone of NTPA (C and H-atoms). This general purpose protocol can be employed not only to investigate the origin of relative stability of any molecular system (e.g., metal complexes) but, in principle, can be used as a predictive tool for, e.g., explaining reaction mechanism. © 2016 Wiley Periodicals, Inc.

Mechanistic study of allopurinol oxidation using aldehyde oxidase, xanthine oxidase and cytochrome P450 enzymes

The oxidation reaction of allopurinol to its active metabolite (oxypurinol) is investigated using the AO and P450 enzymes. To the contrary of AO (and XO), the P450 enzyme can metabolize the allopurinol with a not self-inhibitory mechanism.