Single-Point Hessian Calculations for Improved Vibrational Frequencies and Rigid-Rotor-Harmonic-Oscillator Thermodynamics

Molecular Recognition and Chiral Discrimination from NMR and Multi-Scale Simulations

Chiral molecules are particularly interesting to the pharmaceutical and agrochemical sectors due to their chemical and physical properties. Separation and identification of enantiomers are critical for a broad range of compounds, and discriminating stereoisomers in solution remains a key analytical challenge. Nuclear Magnetic Resonance (NMR) with Matrix-Assisted Diffusion-Ordered Spectroscopy (MAD), in the presence of chiral resolving agents, has emerged as a tool to explore these chiral mixtures. However, insight into the intermolecular interactions that lead to chiral recognition is still limited. Here, we combine experimental MAD studies with computational approaches to investigate the enantioselective discrimination of Mandelic Acid (MA) enantiomers using (R)-BINOL and (S)-BINOL. Molecular dynamics simulations explain the differences in diffusion coefficients for heterochiral complexes. Furthermore, quantum mechanical calculations confirmed enantioselective binding preferences due to differences in Gibbs free energies, highlighting the fundamental interactions and structural criteria that explain the NMR shielding and the diffusion trends. This integrated approach bridges experimental and theoretical studies, offering a comprehensive understanding of chiral recognition mechanisms and elucidating the observed heterochiral preference of BINOL for MA enantiomers. Our findings advance the field of chiral analysis and lay a foundation for future developments for identifying stereoisomers and recognition modes underlying enantioselective binding.

Heterophyllin B: Combining Isotropic and Anisotropic NMR for the Conformational Analysis of a Natural Occurring Cyclic Peptide

Heterophyllin B is a natural occurring cyclic peptide with diverse attributed bioactivities. NMR-based conformational analysis of cyclic peptides often poses a challenge due to limited isotropic solution-state NMR data. In this study, we combined isotropic and anisotropic NMR observables including J-coupling, NOEs, amide proton temperature coefficients, and residual dipolar couplings (RDCs), which enabled the determination of a minimal conformational ensemble of heterophyllin B in methanol at density functional theory (DFT) accuracy. For conformational sampling of a cyclic peptide with a high degree of conformational freedom, we proposed a computational strategy that combines the Conformer-Rotamer Ensemble Sampling Tool (CREST) with the Commandline Energetic SOrting (CENSO). This combined computational and NMR-based approach offers a robust framework for the conformational analysis of cyclic peptides.

Exploring the Mechanism of β-Cyclodextrin-Encased Phenolic Acids Functionalized with TPP for Antioxidant Activity and Targeting

Oxidative stress on the mitochondria in a human cell is attributed to several life-risking conditions, and as such, the importance of molecular structures packed with antioxidant properties and structural characteristics to enter the cell to help prevent such stress has been substantially relevant in recent years. In this study, we investigated the antioxidant properties of triphenylphosphonium (TPP)-conjugated phenolic acids encapsulated in β-cyclodextrin (β-CD). We synthesized TPP conjugates of caffeic, coumaric, and cinnamic acids and formed inclusion complexes with β-CD. Our results showed successful encapsulation of TPP conjugates in β-CD with high efficiency. The TPP conjugates maintained antioxidant activity, with slight reductions observed in β-CD complexes. Furthermore, cell viability studies showed low cytotoxicity of the dds. Computational analyses revealed that TPP conjugation preserved the chemical reactivity of the phenolic acids. Molecular dynamics simulations demonstrated stable inclusion complexes with β-CD and the free energy calculations indicated that TPP conjugation significantly enhanced the ability of caffeic acid to translocate across mitochondrial membranes. These results highlight the potential of TPP-conjugated phenolic acids encapsulated in β-CD as effective antioxidants with improved mitochondrial targeting capabilities.

QCxMS2 - a program for the calculation of electron ionization mass spectra via automated reaction network discovery

We present a new fully-automated computational workflow for the calculation of electron ionization mass spectra by automated reaction network discovery, transition state theory and Monte-Carlo simulations. Compared to its predecessor QCxMS [S. Grimme, Angew. Chem., Int. Ed., 52, 6306-6312] based on extensive molecular dynamics (MD) simulations, QCxMS2's more efficient approach of using stationary points on the potential energy surface (PES) enables the usage of accurate quantum chemical methods. Fragment geometries and reaction paths are optimized with fast semi-empirical quantum mechanical (SQM) methods and reaction barriers are refined at the density functional theory (DFT) level. This composite approach using GFN2-xTB geometries in combination with energies at the ωB97X-3c level proved to be an efficient combination. On a small but diverse test set of 16 organic and inorganic molecules, QCxMS2 spectra are more accurate than ones from QCxMS yielding on average a higher mass spectral matching of 0.700 compared to QCxMS with 0.622, and is more robust with a minimal matching of 0.498 versus 0.100. Further improvements were observed when both geometries and energies were computed at the ωB97X-3c level, yielding an average matching score of 0.730 and a minimal score of 0.527. Due to its higher accuracy and robustness while maintaining computational efficiency, we propose QCxMS2 as a complementary, more reliable and systematically improvable successor to QCxMS for elucidating fragmentation pathways and predicting electron ionization mass spectra of unknown chemical substances, e.g., in analytical chemistry applications. If coupled to currently developed improved SQM methods, QCxMS2 opens an efficient route to accurate, and routine mass spectra predictions. The QCxMS2 program suite is freely available on GitHub.

Investigating Interaction Dynamics of an Enantioselective Peptide-Catalyzed Acylation Reaction

Modern nuclear magnetic resonance (NMR) methods like carbon relaxation dispersion in the rotating frame (13C-R1ρ) and proton chemical exchange saturation transfer (1H-CEST) are key methods to investigate molecular recognition in biomacromolecules and to detect molecular motions on the μs to s timescale, revealing transient conformational states. Changes in kinetics can be linked to binding, folding, or catalytic events. Here, we investigated whether these methods allow detection of changes in the dynamics of a small, highly selective peptide catalyst during recognition of its enantiomeric substrates. The flexible tetrapeptide Boc-l-(π-Me)-His-AGly-l-Cha-l-Phe-OMe, used for the monoacetylation of cycloalkane-diols, is probed at natural abundance using 13C-R1ρ and 1H-CEST. Indeed, we detected differences in dynamics of the peptide upon interaction with the diol. Importantly, these differ depending on the enantiomer of the substrate used. These enantiospecific influences of the substrates on the dynamics of the peptide are rationalized using computational techniques. We find that even though one enantiomer reacts faster, as confirmed by reaction monitoring, the other is more tightly bound in DCM (as confirmed by 1H-saturation transfer difference (STD) measurements). These findings provide insights into the recognition of the substrates and explain the selectivity differences observed between the solvents toluene and DCM.

End-Point Affinity Estimation of Galectin Ligands by Classical and Semiempirical Quantum Mechanical Potentials

The use of quantum mechanical potentials in protein-ligand affinity prediction is becoming increasingly feasible with growing computational power. To move forward, validation of such potentials on real-world challenges is necessary. To this end, we have collated an extensive set of over a thousand galectin inhibitors with known affinities and docked them into galectin-3. The docked poses were then used to systematically evaluate several modern force fields and semiempirical quantum mechanical (SQM) methods up to the tight-binding level under consistent computational workflow. Implicit solvation models available with the tested methods were used to simulate solvation effects. Overall, the best methods in this study achieved a Pearson correlation of 0.7-0.8 between the computed and experimental affinities. There were differences between the tested methods in their ability to rank ligands across the entire ligand set as well as within subsets of structurally similar ligands. A major discrepancy was observed for a subset of ligands that bind to the protein via a halogen bond, which was clearly challenging for all the tested methods. The inclusion of an entropic term calculated by the rigid-rotor-harmonic-oscillator approximation at SQM level slightly worsened correlation with experiment but brought the calculated affinities closer to experimental values. We also found that the success of the prediction strongly depended on the solvation model. Furthermore, we provide an in-depth analysis of the individual energy terms and their effect on the overall prediction accuracy.

Calculation of Adsorbate Free Energy Using the Damping Function Method

Adsorbate free energies are important parameters in surface chemistry and catalysis. Because of its simplicity, the harmonic oscillator (HO) model remains the most widely used method for calculating adsorbate free energy in many fields, including microkinetic modeling. However, it is well-known that the HO method is ineffective for weak adsorption. In this study, we propose a translational model with a diffusion barrier to calculate the state functions of near free translation. Furthermore, an effective mass is introduced in this model. To address hindered translation uniformly, a diffusion barrier-based damping function (DF) is proposed that effectively links the harmonic vibration and translation limits. Adsorbates are divided into three categories according to their adsorption strength and diffusion barrier height. Adsorbed hydrogen atoms have a strong binding energy and relatively high vibrational frequency but a low diffusion barrier. The HT and our proposed DF methods predict that the adsorbed hydrogen atoms behave as translation above room temperature, while the previous DF method predicts that they behave as vibration at any temperature. At last, the dehydrogenation reaction of propane on the Pt(111) surface was taken as an example to illustrate the influence of different methods on the thermodynamic functions.

Expanding the Repertoire of Large Scaffolds with Syn and Anti Macrocyclic Metacyclophanes

Understanding how changes in structure translate to changes in molecular shape is key to catalyst optimization and molecular design in medicinal chemistry and materials. One key contributor to the molecular shape is the relative orientation of substituents on a scaffold. Macrocyclic metacyclophanes display their two arenes in a parallel or antiparallel fashion, resulting in anti or syn conformations that lead to disparate relative orientations of the aryl substituents. This work reports the synthesis of new 14- and 16-membered metacyclophanes and the elucidation of their anti/syn preferences by 1H NMR and computational conformational analysis. Most metacyclophanes studied herein display a strong anti or syn preference and, thus, have well-defined substituent orientations. We propose that anti/syn conformational preferences arise from the minimization of torsional strain along the backbone of the macrocycle, which leads to the prediction that metacyclophanes with remote aryl substituents will adopt the same conformation as their unsubstituted counterparts. Exit vector analysis also reveals that anti-metacyclophanes project their substituents into regions in three-dimensional space that are not accessed by other common large scaffolds, e.g., [2.2]paracyclophanes and ferrocenes. This work also demonstrates how ring size and functional groups, two parameters commonly optimized in macrocycle design, can be used to tune molecular shape.

Structure and intermolecular interactions in ionic liquid 1-ethyl-3-methylimidazolium bromide and its aqueous solutions investigated by vibrational spectroscopy and quantum chemical computations

The recently developed efficient protocol to explicit quantum mechanical modeling of structure and IR spectra of liquids and solutions (S. A. Katsyuba, S. Spicher, T. P. Gerasimova, S. Grimme, J. Phys. Chem. B 2020, 124, 6664) is applied to ionic liquid (IL) 1-ethyl-3-methylimidazolium bromide (EmimBr), its C2-deuterated analog [Emim-d]Br and its aqueous solutions. It is shown that the solvation strongly modifies frequencies and IR intensities of the CH/CD stretching vibrations (νCH/νCD) of the imidazolium ring. The main vibrational spectroscopic features of the neat IL are reproduced by the simulations for a cluster (EmimBr)9, in which all three imidazolium CH moieties of the solvated cation form short contacts with three Br- anions, and another two Br- anions are located on top and bottom of imidazolium ring. Cluster models of aqueous solutions reproduce the experimental vibrational frequencies of actual solutions, provided that the Br- anion of solvated contact ion pair (CIP) is situated on top of imidazolium ring, and CH/CD moieties of the latter participate in short contacts with surrounding water molecules. Both structural and spectroscopic analysis allow to interpret the short contacts CH/CD⋯Br- and CH/CD⋯OH2 as hydrogen bonds of approximately equal strength. Enthalpies of bonding of these liquid-state H-bonds, estimated with the use of empirical correlations, amount to ca. 1.4 kcal⋅mol-1, while the analogous estimates obtained for the gas-phase charged species [Emim]2Br+ increase to 5.6 kcal⋅mol-1. It is shown that formation of solvent-shared ion pair (SIP) in aqueous solution, where the counterions of IL are separated by two water molecules H-bonded to a Br- anion, produces frequency shifts ΔνCH/CD, strongly different from the case of CIP formation. This difference can be used for IR/Raman spectroscopic differentiation of the type of solvated ion pairs of EmimBr or other related ILs.

Reaction-Free Energies for Complexation of Carbohydrates by Tweezer Diboronic Acids

The accurate calculation of reaction-free energies (ΔrG°) for diboronic acids and carbohydrates is challenging due to reactant flexibility and strong solute-solvent interactions. In this study, these challenges are addressed with a semiautomatic workflow based on quantum chemistry methods to calculate conformational free energies, generate microsolvated solute structural ensembles, and compute ΔrG°. Workflow parameters were optimized for accuracy and precision while controlling computational costs. We assessed the accuracy by studying three reactions of diboronic acids with glucose and galactose, finding that the conformational entropy contributes significantly (by 3-5 kcal/mol at room temperature). Explicit solvent molecules improve the computed ΔrG° accuracy by about 4 kcal/mol compared to experimental data, though using 13 or more water molecules reduced precision and increased computational overhead. After fine-tuning, the workflow demonstrated remarkable accuracy, with an absolute error of about 2 kcal/mol compared to experimental ΔrG° and an average interquartile range of 2.4 kcal/mol. These results highlight the workflow's potential for designing and screening tweezer-like ligands with tailored selectivity for various carbohydrates.

Emerging Conformational-Analysis Protocols from the RTCONF55-16K Reaction Thermochemistry Conformational Benchmark Set

RTCONF55-16K is a new, reactive conformational data set based on cost-efficient methods to assess different conformational analysis protocols. Our reference calculations underpinned the accuracy of the CENSO (Grimme et al. J. Phys. Chem. A, 2021, 125, 4039) procedure and resulted in alternative recipes with different cost-accuracy compromises. Our general-purpose and economical protocols (CENSO-light and zero, respectively) were found to be 10-30 times faster than the original algorithm, adding only 0.4-0.7 kcal/mol absolute error to the relative free energy estimates.

Dissociation of HBr in Water Clusters Based on a Hybrid Density Functional Approach

The dissociation of acidic molecules within a microscopic water environment is crucial for understanding intermolecular interactions such as hydrogen bonding. This study explores the optimal configurations of HBr(H2O)n=1-7 using hybrid density functional theory. According to the different mixed cluster structures, the corresponding HBr bond lengths, single-point energies, and introduced proton-transfer parameters are computed and analyzed. The findings indicate that a minimum of three water molecules is necessary for the dissociation of HBr. Subsequently, this conclusion is reinforced through the decomposition of energy components between the acid molecule and water clusters, calculation of hydrogen bonding energies, and analysis of vibrational infrared spectroscopy.

Chemical accuracy for ligand-receptor binding Gibbs energies through multi-level SQM/QM calculations

Calculating the Gibbs energies of binding of ligand-receptor systems with a thermochemical accuracy of ± 1 kcal mol-1 is a challenge to computational approaches. After exploration of the conformational space of the host, ligand and their resulting complexes upon coordination by semi-empirical GFN2 MD and meta-MD simulations, the systematic refinement through a multi-level improvement of binding modes in terms of electronic energies and solvation is able to give Gibbs energies of binding of drug molecules to CB[8] and β-CD macrocyclic receptors with such an accuracy. The accurate treatment of a small number of structures outperforms system-specific force-matching and alchemical transfer model approaches without an extensive sampling and integration.

Molecular folding governs switchable singlet oxygen photoproduction in porphyrin-decorated bistable rotaxanes

Rotaxanes are mechanically interlocked molecules where a ring (macrocycle) is threaded onto a linear molecule (thread). The position of the macrocycle on different stations on the thread can be controlled in response to external stimuli, making rotaxanes applicable as molecular switches. Here we show that bistable rotaxanes based on the combination of a Zn(II) tetraphenylporphyrin photosensitizer, attached to the macrocycle, and a black-hole-quencher, attached to the thread, are capable of singlet oxygen production which can be switched on/off by the addition of base/acid. However, we found that only a sufficiently long linker between both stations on the thread enabled switchability, and that the direction of switching was inversed with regard to the original design. This unexpected behavior was attributed to intramolecular folding of the rotaxanes, as indicated by extensive theoretical calculations. This evidences the importance to take into account the conformational flexibility of large molecular structures when designing functional switchable systems.

Computational 3D Models of Fe2+/3+-Aβ1-42 Complexes Associated with Alzheimer's Disease

The interaction between iron and amyloid-beta (Aβ) peptides has received significant attention in Alzheimer's disease (AD) research due to its potential implications in developing this pathology. However, the coordination preferences of iron and Aβ1-42 have not been thoroughly investigated or remain unknown. This study employs a computational protocol that combines homology modeling techniques with quantum mechanics (DTF-xTB) calculations to build and evaluate several 3D models of Fe2+/3+-Aβ1-42. Our results reveal well-defined complexes for both the metal and peptide moieties, and we discuss the molecular interactions stabilizing these complexes by elucidating the coordinating environments and binding preferences. These proposed models offer valuable insights into the role of iron in Alzheimer's disease (AD) pathology.

Computational analysis of R-X oxidative addition to Pd nanoparticles

Oxidative addition (OA) is a necessary step in mechanisms of widely used synthetic methodologies such as the Heck reaction, cross-coupling reactions, and the Buchwald-Hartwig amination. This study pioneers the exploration of OA of aryl halide to palladium nanoparticles (NPs), a process previously unaddressed in contrast to the activity of well-studied Pd(0) complexes. Employing DFT modeling and semi-empirical metadynamics simulations, the oxidative addition of phenyl bromide to Pd nanoparticles was investigated in detail. Energy profiles of oxidative addition to Pd NPs were analyzed and compared to those involving Pd(0) complexes forming under both ligand-stabilized (phosphines) and ligandless (amine base) conditions. Metadynamics simulations highlighted the edges of the (1 1 1) facets of Pd NPs as the key element of oxidative addition activity. We demonstrate that OA to Pd NPs is not only kinetically facile at ambient temperatures but also thermodynamically favorable. This finding accentuates the necessity of incorporating OA to Pd NPs in future investigations, thus providing a more realistic view of the involved catalytic mechanisms. These results enhance the understanding of aryl halide (cross-)coupling reactions, reinforcing the concept of a catalytic "cocktail". This concept posits dynamic interconversions between diverse active and inactive centers, collectively affecting the outcome of the reaction. High activity of Pd NPs in direct C-X activation paves the way for novel approaches in catalysis, potentially enhancing the field and offering new catalytic pathways to consider.

NBN- and BNB-Phenalenyls: the Yin and Yang of Heteroatom-doped π Systems

NBN- and BNB-doped phenalenyls are isoelectronic to phenalenyl anions and cations, respectively. They represent a pair of complementary molecules that have essentially identical structures but opposite properties as electron donors and acceptors. The NBN-phenalenyls 1-4 considered here were prepared from N,N'-dimethyl-1,8-diaminonaphthalene and readily available boron-containing building blocks (i. e., BH3⋅SMe2 (1), p-CF3-C6H4B(OH)2 (2), C6H5B(OH)2 (3), or MesBCl2/iPr2NEt (4)). Treatment of 1 with 4-Me2N-2,6-Me2-C6H2Li gave the corresponding NBN derivative 5. The BNB-phenalenyl 6 was synthesized from 1,8-naphthalenediyl-bridged diborane(6), PhNH2, and MesMgBr. A computational study reveals that the photoemission of 1, 4, and 5 originates from locally excited (LE) states at the NBN-phenalenyl fragments, while that of 2 is dominated by charge transfer (CT) from the NBN-phenalenyl to the p-CF3-C6H4 fragment. Depending on the dihedral angle θ between its Ph and NBN planes, compound 3 emits mainly from a less polar LE (θ >55°) or more polar CT state (θ <55°). In turn, the energetic preference for either state is governed by the polarity of the solvent used. An equimolar aggregate of the NBN- and BNB-phenalenyls 3 and 6 (in THF/H2O) shows a distinct red-shifted emission compared to that of the individual components, which originates from an intermolecular CT state.

Copper-catalysed electrophilic carboamination of terminal alkynes with benzyne looked at through the computational lens

A detailed computational mechanistic study of the copper-catalysed three-component-type electrophilic carboamination of terminal alkynes with benzyne and an archetypal O-benzoylhydroxylamine electrophile is presented. Probing various plausible pathways for relevant elementary steps and scrutinising performance degradation pathways, with the aid of a reliable computational protocol applied to a realistic catalyst model combined with kinetic analysis, identified the pathways preferably traversed in productive catalysis. It entails rapid alkynylcupration of in situ generated benzyne to deliver the arylcopper nucleophile that undergoes amination with the O-benzoylhydroxylamine electrophile to afford copper benzoate. Umpolung-enabled electrophilic amination favours a multistep SN2-type oxidative addition/N-C bond-forming reductive elimination sequence involving a short-lived formal {P^P}CuIII carboxylate amido aryl intermediate. SN2-type displacement of the benzoate leaving group at the arylcopper nucleophile, which represents the catalyst resting state, is predicted to be the turnover limiting step. Alkynolysis transforms copper benzoate back to catalytically competent alkynylcopper. The computational probe of a wider range of substrates reveals that only severely ring-strained C6-arynes, C6-cycloalkynes and electron-deficient cyclopropenes featuring a highly reactive C≡C linkage could replace benzyne. Moreover, strict control of stationary benzyne concentration is indispensable for electrophilic carboamination to ever become achievable.

Solvate Suite: A Command-Line Interface for Molecular Simulations and Multiscale Microsolvation Modeling

In this work, we introduce the Solvate Suite, a comprehensive and modular command-line interface designed for molecular simulation and microsolvation modeling. The suite interfaces with widely used scientific software, streamlining computational experiments for liquid systems through the automated creation of simulation boxes and topology with adjustable simulation parameters. Furthermore, it has features for graphical and statistical analysis of simulated properties and extraction of trajectory configurations with various filters. Additionally, it introduces innovative strategies for microsolvation modeling with a multiscale approach, employing equilibrated dynamics to identify favorable solute-solvent interactions and enabling full cluster optimization for free-energy calculations without imaginary frequency contamination.

Structural flexibility of favipiravir and its structural analogues in solutions: experimental and computational insight

Combined UV-vis and quantum chemical studies of the structural flexibility and tautomerism of 6-R-3-hydroxy-2-pyrazine carboxamides in solutions revealed that their keto-enol transformations are accompanied by the deprotonation of enol tautomers and the formation of the corresponding anionic species. Both the solvent and the 6-R substituent strongly influence the relative abundance of the above forms in solutions. Anions are not formed in 1,2-dichloroethane (DCE), but the probability of deprotonation in neutral water and N,N-dimethylformamide (DMF) increases in the order R = H < F < NO2. Only enol tautomers of all solutes are found in DCE. DMF stabilizes keto forms only moderately and assists much strongly in the deprotonation of all three compounds. Water tends to stabilize both keto tautomers and deprotonated anions: the keto form dominates in the case of R = H (antiviral drug T-1105), the anions are found exclusively for R = NO2, and the aqueous solution of another antiviral drug, favipiravir (R = F), contains both the keto tautomer and the anionic form. The results of quantum chemical free energy calculations are in agreement with the experimental observations.

Acidity constants and protonation sites of cyclic dinucleotides determined by capillary electrophoresis, quantum chemical calculations, and NMR spectroscopy

Cyclic dinucleotides (CDNs) are important second messengers in bacteria and eukaryotes. Detailed characterization of their physicochemical properties is a prerequisite for understanding their biological functions. Herein, we examine acid-base and electromigration properties of selected CDNs employing capillary electrophoresis (CE), density functional theory (DFT), and nuclear magnetic resonance (NMR) spectroscopy to provide benchmark pKa values, as well as to unambiguously determine the protonation sites. Acidity constants (pKa) of the NH+ moieties of adenine and guanine bases and actual and limiting ionic mobilities of CDNs were determined by nonlinear regression analysis of the pH dependence of their effective electrophoretic mobilities measured by CE in aqueous background electrolytes in a wide pH range (0.98-11.48), at constant temperature (25°C), and constant ionic strength (25 mM). The thermodynamic pKa values were found to be in the range 3.31-4.56 for adenine and 2.28-3.61 for guanine bases, whereas the pKa of enol group of guanine base was in the range 10.21-10.40. Except for systematic shifts of ∼2 pKa, the pKa values calculated by the DFT-D3//COSMO-RS composite protocol that included large-scale conformational sampling and "cross-morphing" were in a relatively good agreement with the pKas determined by CE and predict N1 atom of adenine and N7 atom of guanine as the protonation sites. The protonation of the N1 atom of adenine and N7 atom of guanine in acidic background electrolytes (BGEs) and the dissociation of the enol group of guanine in alkaline BGEs was confirmed also by NMR spectroscopy.

Computationally assisted vibrational spectroscopy of nucleic acid bases. 2. Thymine

As in the case of cytosine [Phys. Chem. Chem. Phys. 2023, 25, 24121-24128], Raman and infrared (IR) spectra of aqueous thymine and its N-deuterated derivative, thymine-d2 have been computationally reproduced and interpreted with the use of the recently developed efficient protocol to explicit quantum mechanical modeling of structure and IR spectra of liquids and solutions [J. Phys. Chem. B, 2020, 124, 6664-6670]. A cluster model of a solute surrounded by 30 water molecules is shown to be sufficient to reproduce experimental vibrational frequencies and relative Raman intensities with the use of B3LYP-D3/def2-TZVP or B3LYP-D3/aug-cc-pVDZ simulations. Analogous PBE-D3 computations provided a less good, but still reasonably accurate, modeling of Raman spectra. It is shown that strong changes of frequencies and relative intensities of the Raman bands of thymine, caused by its hydration, can be interpreted mainly as a result of hydrogen bonding with 6 nearest water molecules. Non-negligible improvement of the quality of simulations for larger clusters comprising water molecules that do not have direct contacts with the solute, suggests that spectroscopic effects of hydration should be ascribed to the joined action of solute-solvent and solvent-solvent interactions. Nevertheless, the moderate number of water molecules required for successful simulations of the Raman spectra of aqueous thymine, suggests that the vibrational modes and derivatives of the polarizability of the solute are mainly locally influenced, while the effect of bulk water is rather modest.

Machine learning-based correction for spin-orbit coupling effects in NMR chemical shift calculations

As one of the most powerful analytical methods for molecular and solid-state structure elucidation, NMR spectroscopy is an integral part of chemical laboratories associated with a great research interest in its computational simulation. Particularly when heavy atoms are present, a relativistic treatment is essential in the calculations as these influence also the nearby light atoms. In this work, we present a Δ-machine learning method that approximates the contribution to 13C and 1H NMR chemical shifts that stems from spin-orbit (SO) coupling effects. It is built on computed reference data at the spin-orbit zeroth-order regular approximation (ZORA) DFT level for a set of 6388 structures with 38 740 13C and 64 436 1H NMR chemical shifts. The scope of the methods covers the 17 most important heavy p-block elements that exhibit heavy atom on the light atom (HALA) effects to covalently bound carbon or hydrogen atoms. Evaluated on the test data set, the approach is able to recover roughly 85% of the SO contribution for 13C and 70% for 1H from a scalar-relativistic PBE0/ZORA-def2-TZVP calculation at virtually no extra computational costs. Moreover, the method is transferable to other baseline DFT methods even without retraining the model and performs well for realistic organotin and -lead compounds. Finally, we show that using a combination of the new approach with our previous Δ-ML method for correlation contributions to NMR chemical shifts, the mean absolute NMR shift deviations from non-relativistic DFT calculations to experimental values can be halved.

Affinity of Telluronium Chalcogen Bond Donors for Lewis Bases in Solution: A Critical Experimental-Theoretical Joint Study

Telluronium salts [Ar2 MeTe]X were synthesized, and their Lewis acidic properties towards a number of Lewis bases were addressed in solution by physical and theoretical means. Structural X-ray diffraction analysis of 21 different salts revealed the electrophilicity of the Te centers in their interactions with anions. Telluroniums' propensity to form Lewis pairs was investigated with OPPh3 . Diffusion-ordered NMR spectroscopy suggested that telluroniums can bind up to three OPPh3 molecules. Isotherm titration calorimetry showed that the related heats of association in 1,2-dichloroethane depend on the electronic properties of the substituents of the aryl moiety and on the nature of the counterion. The enthalpies of first association of OPPh3 span -0.5 to -5 kcal mol-1 . Study of the affinity of telluroniums for OPPh3 by state-of-the-art DFT and ab-initio methods revealed the dominant Coulombic and dispersion interactions as well as an entropic effect favoring association in solution. Intermolecular orbital interactions between [Ar2 MeTe]+ cations and OPPh3 are deemed insufficient on their own to ensure the cohesion of [Ar2 MeTe ⋅ Bn ]+ complexes in solution (B=Lewis base). Comparison of Grimme's and Tkatchenko's DFT-D4/MBD-vdW thermodynamics of formation of higher [Ar2 MeTe ⋅ Bn ]+ complexes revealed significant molecular size-dependent divergence of the two methodologies, with MBD yielding better agreement with experiment.

Minimal Peptoid Dynamics Inform Self-Assembly Propensity

Peptoids are structural isomers of natural peptides, with side chain attachment at the amide nitrogen, conferring this class of compounds with the ability to access both cis and trans ω torsions as well as an increased diversity of ψ/φ states with respect to peptides. Sampling within these dimensions is controlled through side chain selection, and an expansive set of viable peptoid residues exists. It has been shown recently that "minimal" di- and tripeptoids with aromatic side chains can self-assemble into highly ordered structures, with size and morphological definition varying as a function of sequence pattern (e.g., XFF and FXF, where X = a nonaromatic peptoid monomer). Aromatic groups, such as phenylalanine, are regularly used in the design of minimal peptide assemblers. In recognition of this, and to draw parallels between these compounds classes, we have developed a series of descriptors for intramolecular dynamics of aromatic side chains to discern whether these dynamics, in a preassembly condition, can be related to experimentally observed nanoscale assemblies. To do this, we have built on the atomistic peptoid force field reported by Weiser and Santiso (CGenFF-WS) through the rigorous fitting of partial charges and the collation of Charmm General Force Field (CGenFF) parameters relevant to these systems. Our study finds that the intramolecular dynamics of side chains, for a given sequence, is dependent on the specific combination of backbone ω torsions and that homogeneity of sampling across these states correlates well with the experimentally observed ability to assemble into nanomorphologies with long-range order. Sequence patterning is also shown to affect sampling, in a manner consistent for both tripeptoids and tripeptides. Additionally, sampling similarities between the nanofiber forming tripeptoid, Nf-Nke-Nf in the cc state, and the nanotube forming dipeptide FF, highlight a structural motif which may be relevant to the emergence of extended linear assemblies. To assess these properties, a variety of computational approaches have been employed.

ConfSolv: Prediction of Solute Conformer-Free Energies across a Range of Solvents

Predicting Gibbs free energy of solution is key to understanding the solvent effects on thermodynamics and reaction rates for kinetic modeling. Accurately computing solution free energies requires the enumeration and evaluation of relevant solute conformers in solution. However, even after generation of relevant conformers, determining their free energy of solution requires an expensive workflow consisting of several ab initio computational chemistry calculations. To help address this challenge, we generate a large data set of solution free energies for nearly 44,000 solutes with almost 9 million conformers calculated in 41 different solvents using density functional theory and COSMO-RS and quantify the impact of solute conformers on the solution free energy. We then train a message passing neural network to predict the relative solution free energies of a set of solute conformers, enabling the identification of a small subset of thermodynamically relevant conformers. The model offers substantial computational time savings with predictions usually substantially within 1 kcal/mol of the free energy of the solution calculated by using computational chemical methods.

Revisiting Homochiral versus Heterochiral Interactions through a Long Detective Story of a Useful Azobis-Nitrile and Puzzling Racemate

This paper documents and reinvestigates the solid-state and crystal structures of 4,4'-azobis-4-cyanopentanoic acid (ACPA), a water-soluble azobis-nitrile of immense utility as a radical initiator in living polymerizations and a labile mechanophore that can be embedded within long polymer chains to undergo selective scission under mechanical activation. Surprisingly, for such applications, both the commercially available reagent and their derivatives are used as "single initiators" when this azonitrile is actually a mixture of stereoisomers. Although the racemate and meso compounds were identified more than half a century ago and their enantiomers were separated by classical resolution, there have been confusing narratives dealing with their characterization, the existence of a conglomeratic phase, and fractional crystallization. Our results report on the X-ray crystal structures of all stereoisomers for the first time, along with further details on enantiodiscrimination and the always intriguing arguments accounting for the stability of homochiral versus heterochiral crystal aggregates. To this end, metadynamic (MTD) simulations on stereoisomer molecular aggregates were performed to capture the incipient nucleation events at the picosecond time scale. This analysis sheds light on the driving homochiral aggregation of ACPA enantiomers.

Comprehensive Mechanistic View of the Hydrolysis of Oxadiazole-Based Inhibitors by Histone Deacetylase 6 (HDAC6)

Histone deacetylase (HDAC) inhibitors used in the clinic typically contain a hydroxamate zinc-binding group (ZBG). However, more recent work has shown that the use of alternative ZBGs, and, in particular, the heterocyclic oxadiazoles, can confer higher isoenzyme selectivity and more favorable ADMET profiles. Herein, we report on the synthesis and biochemical, crystallographic, and computational characterization of a series of oxadiazole-based inhibitors selectively targeting the HDAC6 isoform. Surprisingly, but in line with a very recent finding reported in the literature, a crystal structure of the HDAC6/inhibitor complex revealed that hydrolysis of the oxadiazole ring transforms the parent oxadiazole into an acylhydrazide through a sequence of two hydrolytic steps. An identical cleavage pattern was also observed both in vitro using the purified HDAC6 enzyme as well as in cellular systems. By employing advanced quantum and molecular mechanics (QM/MM) and QM calculations, we elucidated the mechanistic details of the two hydrolytic steps to obtain a comprehensive mechanistic view of the double hydrolysis of the oxadiazole ring. This was achieved by fully characterizing the reaction coordinate, including identification of the structures of all intermediates and transition states, together with calculations of their respective activation (free) energies. In addition, we ruled out several (intuitively) competing pathways. The computed data (ΔG^‡ ≈ 21 kcal·mol^-1 for the rate-determining step of the overall dual hydrolysis) are in very good agreement with the experimentally determined rate constants, which a posteriori supports the proposed reaction mechanism. We also clearly (and quantitatively) explain the role of the -CF₃ or -CHF₂ substituent on the oxadiazole ring, which is a prerequisite for hydrolysis to occur. Overall, our data provide compelling evidence that the oxadiazole warheads can be efficiently transformed within the active sites of target metallohydrolases to afford reaction products possessing distinct selectivity and inhibition profiles.

Compilation of Ionic Clusters with the Rock Salt Structure: Accurate Benchmark Thermochemical Data, Assessment of Quantum Chemistry Methods, and the Convergence Behavior of Lattice Energies

In the present study, computational quantum chemistry is used to obtain lattice energies (LEs) for a range of ionic clusters with the NaCl structure. Specifically, the compounds include NaF, NaCl, MgO, MgS, KF, CaO, and CaS clusters, (MX)_n, with n = 1, 2, 4, 6, 8, 12, 16, 24, 32, 40, 50, 60, 75, 90, and 108. The highest-level W2 and W1X-2 methods are applied to the small clusters with n = 1 to 8 (the MX35 data set). The assessment with MX35 shows that, for the calculation of geometries and vibrational frequencies, the PBE0-D3(BJ) and PBE-D3(BJ) DFT methods are reasonable, but the calculation of atomization energies is more challenging. This is a result of different systematic deviations for clusters of different species. Thus, species-specific adjustments are applied for larger clusters, which are calculated with the DuT-D3 double-hybrid DFT method, the MN15 DFT method, and the PM7 semi-empirical method. They yield smoothly converging LEs to the bulk values. It is also found that, for the alkali-metal species, the LEs for a single molecule are ∼70% of the bulk values, while for the alkali-earth species, they are ∼80%. This has enabled a straightforward means to the first-principles estimation of LEs for similarly structured ionic compounds.

Dissecting Reaction Mechanisms and Catalytic Contributions in Flavoprotein Fumarate Reductases

The interconversion between fumarate and succinate is fundamental to the energy metabolism of nearly all organisms. This redox reaction is catalyzed by a large family of enzymes, fumarate reductases and succinate dehydrogenases, using hydride and proton transfers from a flavin cofactor and a conserved Arg side-chain. These flavoenzymes also have substantial biomedical and biotechnological importance. Therefore, a detailed understanding of their catalytic mechanisms is valuable. Here, calibrated electronic structure calculations in a cluster model of the active site of the Fcc₃ fumarate reductase were employed to investigate various reaction pathways and possible intermediates in the enzymatic environment and to dissect interactions that contribute to catalysis of fumarate reduction. Carbanion, covalent adduct, carbocation, and radical intermediates were examined. Significantly lower barriers were obtained for mechanisms via carbanion intermediates, with similar activation energies for hydride and proton transfers. Interestingly, the carbanion bound to the active site is best described as an enolate. Hydride transfer is stabilized by a preorganized charge dipole in the active site and by the restriction of the C1-C2 bond in a twisted conformation of the otherwise planar fumarate dianion. But, protonation of a fumarate carboxylate and quantum tunneling effects are not critical for catalysis of the hydride transfer. Calculations also suggest that the driving force for enzyme turnover is provided by regeneration of the catalytic Arg, either coupled with flavin reduction and decomposition of a proposed transient state or directly from the solvent. The detailed mechanistic description of enzymatic reduction of fumarate provided here clarifies previous contradictory views and provides new insights into catalysis by essential flavoenzyme reductases and dehydrogenases.

NDIPhos as a platform for chiral supramolecular ligands in rhodium-catalyzed enantioselective hydrogenation

Chiral naphthalene diimide ligands (NDIPhos) were exploited in rhodium-catalyzed enantioselective hydrogenation. The key feature of these ligands is their ability to self-assemble via π-π interactions to mimic bidentate ligands, offering a complementary method to traditional supramolecular strategies. This concept was further substantiated by computations with the composite electronic-structure method r²SCAN-3c.

The Interplay of Weakly Coordinating Anions and the Mechanical Bond: A Systematic Study of the Explicit Influence of Counterions on the Properties of (Pseudo)rotaxanes

Weakly coordinating anions (WCAs) have attracted much attention in recent years due to their ability to stabilise highly reactive cations. It may well be argued, however, that a profound understanding of what truly defines a WCA is still lacking, and systematic studies to unravel counterion effects are scarce. In this work, we investigate a supramolecular pseudorotaxane formation reaction, subject to a selection of anions, ranging from strongly to weakly coordinating, which not only aids in fostering our knowledge about anion coordination properties, but also provides valuable theoretical insight into the nature of the mechanical bond. We employ state-of-the-art DFT-based methods and tools, combined with isothermal calorimetry and 1H NMR experiments, to compute anion-dependent Gibbs free association energies ΔGa, as well as to evaluate intermolecular interactions. We find correlations between ΔGa and the anions’ solvation energies, which are exploited to calculate physico-chemical reaction parameters in the context of coordinating anions. Furthermore, we show that the binding situation within the (pseudo)rotaxanes can be mostly understood by straight-forward electrostatic considerations. However, quantum-chemical effects such as dispersion and charge-transfer interactions become more and more relevant when WCAs are employed.

Probing a Hydrogen-π Interaction Involving a Trapped Water Molecule in the Solid State

The detection and characterization of trapped water molecules in chemical entities and biomacromolecules remains a challenging task for solid materials. We herein present proton-detected solid-state Nuclear Magnetic Resonance (NMR) experiments at 100 kHz magic-angle spinning and at high static magnetic-field strengths (28.2 T) enabling the detection of a single water molecule fixed in the calix[4]arene cavity of a lanthanide complex by a combination of three types of non-covalent interactions. The water proton resonances are detected at a chemical-shift value close to zero ppm, which we further confirm by quantum-chemical calculations. Density Functional Theory calculations pinpoint to the sensitivity of the proton chemical-shift value for hydrogen-π interactions. Our study highlights how proton-detected solid-state NMR is turning into the method-of-choice in probing weak non-covalent interactions driving a whole branch of molecular-recognition events in chemistry and biology.

Supramolecular Optimization of Sensory Function of a Hemicurcuminoid through Its Incorporation into Phospholipid and Polymeric Polydiacetylenic Vesicles: Experimental and Computational Insight

This work presents the synthesis of a new representative of hemicurcuminoids with a nonyloxy substituent (HCur) as a fluorescent amphiphilic structural element of vesicular aggregates based on phosphatidylcholine (PC), phosphatidylserine (PS), and 10,12-pentacosadiynoic acid (PCDA). Both X-ray diffraction analysis of the single crystal and ¹H NMR spectra of HCur in organic solvents indicate the predominance of the enol-tautomer of HCur. DFT calculations show the predominance of the enol tautomer HCur in supramolecular assemblies with PC, PS, and PCDA molecules. The results of the molecular modeling show that HCur molecules are surrounded by PC and PS with a rather weak exposure to water molecules, while an exposure of HCur molecules to water is enhanced under its supramolecular assembly with PCDA molecules. This is in good agreement with the higher loading of HCur into PC(PS) vesicles compared to PCDA vesicles converted into polydiacetylene (PDA) ones by photopolymerization. HCur molecules incorporated into HCur-PDA vesicles exhibit greater planarity distortion and hydration effect in comparison with HCur-PC(PS) ones. HCur-PDA is presented as a dual fluorescence-chromatic nanosensor responsive to a change in pH within 7.5-9.5, heavy metal ions and polylysine, and the concentration-dependent fluorescent response is more sensitive than the chromatic one. Thus, the fluorescent response of HCur-PDA allows for the distinguishing between Cd²⁺ and Pb²⁺ ions in the concentration range 0-0.01 mM, while the chromatic response allows for the selective sensing of Pb²⁺ over Cd²⁺ ions at their concentrations above 0.03 mM.

Quantum chemical calculation of the vapor pressure of volatile and semi volatile organic compounds

The vapor pressure is a specific and temperature-dependent parameter that describes the volatility of a substance and thus its driving force for evaporation or sublimation into the gas phase. Depending on the magnitude of the vapor pressure, there are different methods for experimental determination. However, these are usually associated with a corresponding amount of effort and become less accurate as the vapor pressure decreases. For purposes of vapor pressure prediction, algorithms were developed that are usually based on quantitative structure-activity relationships (QSAR). The quantum mechanical (QM) approach followed here applies an alternative, much less empirical strategy, where the change in Gibbs free energy for the transition from the condensed to the gas phase is obtained from conformer ensembles computed for each phase separately. The results of this automatic, so-called CRENSO workflow are compared with experimentally determined vapor pressures for a large set of environmentally relevant compounds. In addition, comparisons are made with the single structure-based COSMO-RS QM approach, linear-free-energy relationships (LFER) as well as results from the SPARC program. We show that our CRENSO workflow is superior to conventional prediction models and provides reliable vapor pressures for liquids and sub-cooled liquids over a wide pressure range.

Electrochemical recycling of homogeneous catalysts

Homogeneous catalysts have rapid kinetics and keen reaction selectivity. However, their widespread use for industrial catalysis has remained limited because of challenges in reusability. Here, we propose a redox-mediated electrochemical approach for catalyst recycling using metallopolymer-functionalized electrodes for binding and release. The redox platform was investigated for the separation of key platinum and palladium homogeneous catalysts used in organic synthesis and industrial chemical manufacturing. Noble metal catalysts for hydrosilylation, silane etherification, Suzuki cross-coupling, and Wacker oxidation were recycled electrochemically. The redox electrodes demonstrated high sorption uptake for platinum-based catalysts (Q_max up to 200 milligrams of platinum per gram of adsorbent) from product mixtures, with up to 99.5% recovery, while retaining full catalytic activity over multiple cycles. The combination of mechanistic studies and electronic structure calculations indicate that selective interactions with anionic intermediates during the catalytic cycle played a key role in the separations. Last, continuous flow cell studies support the scalability and favorable technoeconomics of electrochemical recycling.

Best-Practice DFT Protocols for Basic Molecular Computational Chemistry

Nowadays, many chemical investigations are supported by routine calculations of molecular structures, reaction energies, barrier heights, and spectroscopic properties. The lion's share of these quantum-chemical calculations applies density functional theory (DFT) evaluated in atomic-orbital basis sets. This work provides best-practice guidance on the numerous methodological and technical aspects of DFT calculations in three parts: Firstly, we set the stage and introduce a step-by-step decision tree to choose a computational protocol that models the experiment as closely as possible. Secondly, we present a recommendation matrix to guide the choice of functional and basis set depending on the task at hand. A particular focus is on achieving an optimal balance between accuracy, robustness, and efficiency through multi-level approaches. Finally, we discuss selected representative examples to illustrate the recommended protocols and the effect of methodological choices.

"Vermellogens" and the Development of CB[8]-Based Supramolecular Switches Using pH-Responsive and Non-Toxic Viologen Analogues

We present herein the "vermellogens", a new class of pH-responsive viologen analogues, which replace the direct linking between para-substituted pyridinium moieties within those by a hydrazone functional group. A series of such compounds have been efficiently synthesized in aqueous media by hydrazone exchange reactions, displaying a marked pH-responsivity. Furthermore, the parent N,N'-dimethylated "vermellogen": the "red thread", an analogue of the herbicide paraquat and used herein as a representative model of the series, showed anion-recognition abilities, non-reversible electrochemical behavior, and non-toxicity of the modified bis-pyridinium core. The host-guest chemistry for the "red thread" with the CB[7,8] macrocyclic receptors has been extensively studied experimentally and by dispersion corrected density functional theory methods, showing a parallel behavior to that previously described for the herbicide but, crucially, swapping the well-known redox reactive capabilities of the viologen-based inclusion complexes by acid-base supramolecular responsiveness.

Portable Models for Entropy Effects on Kinetic Selectivity

Differences in entropies of competing transition states can direct kinetic selectivity. Understanding and modeling such entropy differences at the molecular level is complicated by the fact that entropy is statistical in nature; i.e., it depends on multiple vibrational states of transition structures, the existence of multiple dynamically accessible pathways past these transition structures, and contributions from multiple transition structures differing in conformation/configuration. The difficulties associated with modeling each of these contributors are discussed here, along with possible solutions, all with an eye toward the development of portable qualitative models of use to experimentalists aiming to design reactions that make use of entropy to control kinetic selectivity.

The Role of Packing, Dispersion, Electrostatics, and Solvation in High-Affinity Complexes of Cucurbit[n]urils with Uncharged Polar Guests

The rationalization of non-covalent binding trends is both of fundamental interest and provides new design concepts for biomimetic molecular systems. Cucurbit[n]urils (CBn) are known for a long time as the strongest synthetic binders for a wide range of (bio)organic compounds in water. However, their host-guest binding mechanism remains ambiguous despite their symmetric and simple macrocyclic structure and the wealth of literature reports. We herein report experimental thermodynamic binding parameters (ΔG, ΔH, TΔS) for CB7 and CB8 with a set of hydroxylated adamantanes, di-, and triamantanes as uncharged, rigid, and spherical/ellipsoidal guests. Binding geometries and binding energy decomposition were obtained from high-level theory computations. This study reveals that neither London dispersion interactions, nor electronic energies or entropic factors are decisive, selectivity-controlling factors for CBn complexes. In contrast, peculiar host-related solvation effects were identified as the major factor for rationalizing the unique behavior and record-affinity characteristics of cucurbit[n]urils.

GEOM, energy-annotated molecular conformations for property prediction and molecular generation

Machine learning (ML) outperforms traditional approaches in many molecular design tasks. ML models usually predict molecular properties from a 2D chemical graph or a single 3D structure, but neither of these representations accounts for the ensemble of 3D conformers that are accessible to a molecule. Property prediction could be improved by using conformer ensembles as input, but there is no large-scale dataset that contains graphs annotated with accurate conformers and experimental data. Here we use advanced sampling and semi-empirical density functional theory (DFT) to generate 37 million molecular conformations for over 450,000 molecules. The Geometric Ensemble Of Molecules (GEOM) dataset contains conformers for 133,000 species from QM9, and 317,000 species with experimental data related to biophysics, physiology, and physical chemistry. Ensembles of 1,511 species with BACE-1 inhibition data are also labeled with high-quality DFT free energies in an implicit water solvent, and 534 ensembles are further optimized with DFT. GEOM will assist in the development of models that predict properties from conformer ensembles, and generative models that sample 3D conformations.

Measurement(s)	Conformer geometries and properties
Technology Type(s)	Computational Chemistry

Hydrocarbon Macrocycle Conformer Ensembles and 13 C-NMR Spectra

NMR as a routine analytical method provides important three-dimensional structure information of compounds in solution. Here we apply the recently released CRENSO computational workflow for the automated generation of conformer ensembles to the quantum mechanical calculation of ¹³ C-NMR spectra of a series of flexible cycloalkanes up to C₂₀ H₄₀ . We evaluate the computed chemical shifts in comparison with corresponding experimental data in chloroform. It is shown that accurate and properly averaged theoretical NMR data can be obtained in about a day of computation time on a standard workstation computer. The excellent agreement between theory and experiment enables one to deduce preferred conformations of large, non-rigid macrocycles under ambient conditions from our automated procedure.

Temperature-sensitive emission of dialkylaminostyrylhetarene dyes and their incorporation into phospholipid aggregates: Applicability for thermal sensing and cellular uptake behavior

A series of dialkylaminostyrylhetarene dyes constructed from electron-rich and electron-deficient moieties of various structures connected via vinylene π-bridges are introduced as temperature-sensitive luminophores. The temperature dependent emission of the dyes in the acidified dichloromethane solutions derives from temperature-induced shift of the equilibrium between neutral and protonated forms of the dyes. The heating-induced blue shift and intensification of emission of neutral form of the dyes make them a promising basis for development of nanoparticles exhibiting temperature-sensitivity in aqueous solutions at pH typical of biological liquids. Hydrophobicity-driven incorporation of the water insoluble dyes into L-α-phosphatidylcholine(PC)-based bilayers allows to obtain water dispersible dye-PC aggregates, and to follow their emission in the aqueous solutions. Structure of the dyes has strong impact on the efficacy of the dyes incorporation into the PC-based bilayers, temperature sensitivity of emission of the dye-PC aggregates and its reversibility under the heating/cooling cycles. This enables structural optimization of the dyes in order to obtain the dye-PC species demonstrating maximal temperature dependence and reversibility of their luminescence in aqueous solutions. The selected leader exhibits low cytotoxicity exemplified for M-HeLa and Chang Liver cell lines, while the efficient cell internalization of the dye, manifested in the staining of the cell cytoplasm, opens further opportunities for biosensing applications.

Computer-aided simulation of infrared spectra of ethanol conformations in gas, liquid and in CCl4 solution

The recently developed efficient protocol combining implicit and explicit, accurate quantum-mechanical modeling of the condensed state (Katsyuba et al., J. Chem. Phys. 155, 024507 [2021]) is used to describe the IR spectra of liquid ethanol and its solutions in CCl₄ . The relative abundance of the anti and gauche conformers of ethanol is shown to increase from ~40:60 in the gas phase to ~55:45 in the liquid phase. In spite of a moderate impact of media effects on the conformational composition of the liquid, the solvent strongly influences vibrational frequencies, IR intensities, and normal modes of each conformer, producing qualitatively different spectra compared to the gas phase and CCl₄ solution. Further, these solvent effects affecting IR frequencies and intensities depend not only on the conformation of the solvated molecule but also on the solvating species. Nevertheless, vibrational frequencies of anti and gauche conformers of liquid ethanol and its several isotopomers practically coincide with each other. Convenient liquid-state conformational markers in the fingerprint region of IR spectra are revealed for the hydroxyl-deuterated species: CH₃ CH₂ OD, CH₃ CHDOD, CH₃ CD₂ OD, and CD₃ CD₂ OD.

Quantum Chemical Calculation and Evaluation of Partition Coefficients for Classical and Emerging Environmentally Relevant Organic Compounds

Octanol/water (K_OW), octanol/air (K_OA), and hexadecane/air (K_HdA) partition coefficients are calculated for 67 organic compounds of environmental concern using computational chemistry. The extended CRENSO workflow applied here includes the calculation of extensive conformer ensembles with semiempirical methods and refinement through density functional theory, taking into account solvation models, especially COSMO-RS, and thermostatistical contributions. This approach is particularly advantageous for describing large and nonrigid molecules. With regard to K_OW and K_HdA, one can refer to many experimental data from direct and indirect measurement methods, and very good matches with results from our quantum chemical workflow are evident. In the case of the K_OA values, however, good matches are only obtained for the experimentally determined values. Larger systematic deviations between data computed here and available, nonexperimental quantitative structure-activity relationship literature data occur in particular for phthalic acid esters and organophosphate esters. From a critical analysis of the coefficients calculated in this work and comparison with available literature data, we conclude that the presented quantum chemical composite approach is the most powerful so far for calculating reliable partition coefficients because all physical contributions to the conformational free energy are considered and the structure ensembles for the two phases are generated independently and consistently.

Efficient calculation of protein-ligand binding free energy with GFN methods: the power of cluster model

The protein-ligand interactions are crucial in many biochemical processes and biomedical applications, yet it still remains challenging to accurately calculating the binding free energy of their interactions. In this work,...