Evaluation of an Efficient 3D-RISM-SCF Implementation as a Tool for Computational Spectroscopy in Solution

Thermodynamic Study of N-Methylformamide and N,N-Dimethyl-Formamide

An extensive thermodynamic study of N-methylformamide (CAS RN: 123-39-7) and N,N-dimethylformamide (CAS RN: 68-12-2), is presented in this work. The liquid heat capacities of N-methylformamide were measured by Tian-Calvet calorimetry in the temperature interval (250-300) K. The vapor pressures for N-methylformamide and N,N-dimethylformamide were measured using static method in the temperature range 238 K to 308 K. The ideal-gas thermodynamic properties were calculated using a combination of the density functional theory (DFT) and statistical thermodynamics. A consistent thermodynamic description was developed using the method of simultaneous correlation, where the experimental and selected literature data for vapor pressures, vaporization enthalpies, and liquid phase heat capacities and the calculated ideal-gas heat capacities were treated together to ensure overall thermodynamic consistency of the results. The resulting vapor pressure equation is valid from the triple point to the normal boiling point temperature.

The accuracy limit of chemical shift predictions for species in aqueous solution

Interpreting NMR experiments benefits from first-principles predictions of chemical shifts. Reaching the accuracy limit of theory is relevant for unambiguous structural analysis and dissecting theoretical approximations. Since accurate chemical shift measurements are based on using internal reference compounds such as trimethylsilylpropanesulfonate (DSS), a detailed comparison of experimental with theoretical data requires simultaneous consideration of both target and reference species ensembles in the same solvent environment. Here we show that ab initio molecular dynamics simulations to generate liquid-state ensembles of target and reference compounds, including explicitly their short-range solvation environments and combined with quantum-mechanical solvation models, allows for predicting highly accurate 1H (∼0.1-0.5 ppm) and aliphatic 13C (∼1.5 ppm) chemical shifts for aqueous solutions of the model compounds trimethylamine N-oxide (TMAO) and N-methylacetamide (NMA), referenced to DSS without any system-specific adjustments. This encompasses the two peptide bond conformations of NMA identified by NMR. The results are used to derive a general-purpose guideline set for predictive NMR chemical shift calculations of NMA in the liquid state and to identify artifacts of force field models. Accurate predictions are only obtained if a sufficient number of explicit water molecules is included in the quantum-mechanical calculations, disproving a purely electrostatic model of the solvent effect on chemical shifts.

Recent developments and applications of reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED): A hybrid model of quantum chemistry and integral equation theory of molecular liquids

The significance of solvent effects in electronic structure calculations has long been noted, and various methods have been developed to consider this effect. The reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED) is a hybrid model that combines the integral equation theory of molecular liquids with quantum chemistry. This method can consider the statistically convergent solvent distribution at a significantly lower cost than molecular dynamics simulations. Because the RISM theory explicitly considers the solvent structure, it performs well for systems where hydrogen bonds are formed between the solute and solvent molecules, which is a challenge for continuum solvent models. Taking advantage of being founded on the variational principle, theoretical developments have been made in calculating various properties and incorporating electron correlation effects. In this review, we organize the theoretical aspects of RISM-SCF-cSED and its distinctions from other hybrid methods involving integral equation theories. Furthermore, we carefully present its progress in terms of theoretical developments and recent applications.

Dissecting the Molecular Origin of g-Tensor Heterogeneity and Strain in Nitroxide Radicals in Water: Electron Paramagnetic Resonance Experiment versus Theory

Nitroxides are common EPR sensors of microenvironmental properties such as polarity, numbers of H-bonds, pH, and so forth. Their solvation in an aqueous environment is facilitated by their high propensity to form H-bonds with the surrounding water molecules. Their g- and A-tensor elements are key parameters to extracting the properties of their microenvironment. In particular, the gxx value of nitroxides is rich in information. It is known to be characterized by discrete values representing nitroxide populations previously assigned to have different H-bonds with the surrounding waters. Additionally, there is a large g-strain, that is, a broadening of g-values associated with it, which is generally correlated with environmental and structural micro-heterogeneities. The g-strain is responsible for the frequency dependence of the apparent line width of the EPR spectra, which becomes evident at high field/frequency. Here, we address the molecular origin of the gxx heterogeneity and of the g-strain of a nitroxide moiety (HMI: 2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl, C9H19N2O) in water. To treat the solvation effect on the g-strain, we combined a multi-frequency experimental approach with ab initio molecular dynamics simulations for structural sampling and quantum chemical EPR property calculations at the highest realistically affordable level, including an explicitly micro-solvated HMI ensemble and the embedded cluster reference interaction site model. We could clearly identify the distinct populations of the H-bonded nitroxides responsible for the gxx heterogeneity experimentally observed, and we dissected the role of the solvation shell, H-bond formation, and structural deformation of the nitroxide in the creation of the g-strain associated with each nitroxide subensemble. Two contributions to the g-strain were identified in this study. The first contribution depends on the number of hydrogen bonds formed between the nitroxide and the solvent because this has a large and well-understood effect on the gxx-shift. This contribution can only be resolved at high resonance frequencies, where it leads to distinct peaks in the gxx region. The second contribution arises from configurational fluctuations of the nitroxide that necessarily lead to g-shift heterogeneity. These contributions cannot be resolved experimentally as distinct resonances but add to the line broadening. They can be quantitatively analyzed by studying the apparent line width as a function of microwave frequency. Interestingly, both theory and experiment confirm that this contribution is independent of the number of H-bonds. Perhaps even more surprisingly, the theoretical analysis suggests that the configurational fluctuation broadening is not induced by the solvent but is inherently present even in the gas phase. Moreover, the calculations predict that this broadening decreases upon solvation of the nitroxide.

Spin-State Splittings in 3d Transition-Metal Complexes Revisited: Toward a Reliable Theory Benchmark

A new composite method for the calculation of spin-crossover energies in 3d transition-metal complexes based on multireference methods is presented. The method reduces to MRCISD+Q at the complete-basis-set (CBS) level for atomic ions, for which it gives excitation energies with a mean absolute error of only ca. 0.01 eV. For molecular complexes, the CASPT2+δMRCI composite approach corresponds to a CASPT2/CBS calculation augmented by a high-level MRCISD+Q-CASPT2 correction with a smaller ligand basis set. For a set of [Fe(He)6]n+ test complexes, the approach reproduces full MRCISD+Q/CBS results to within better than 0.04 eV, without depending on any arbitrary IPEA shifts. The high-quality CASPT2+δMRCI method has then been applied to a series of 3d transition-metal hexaqua complexes in aqueous solution, augmented by an elaborate 3D-RISM-SCF solvent treatment of the underlying structures. It provides unprecedented agreement with experiment for the lowest-lying vertical spin-flip excitation energies, except for the Fe³⁺ system. Closer examination of the latter case provides strong evidence that the observed lowest-energy excitation at 1.56 eV, which has been used frequently for evaluating quantum-chemical methods, does not arise from the iron(III) hexaqua complex in solution, but from its singly deprotonated counterpart, [Fe(H2O)5OH]2+.

Reference interaction site model self-consistent field with constrained spatial electron density approach for nuclear magnetic shielding in solution

We propose a new hybrid approach combining quantum chemistry and statistical mechanics of liquids for calculating the nuclear magnetic resonance (NMR) chemical shifts of solvated molecules. Based on the reference interaction site model self-consistent field with constrained spatial electron density distribution (RISM–SCF–cSED) method, the electronic structure of molecules in solution is obtained, and the expression for the nuclear magnetic shielding tensor is derived as the second-order derivative of the Helmholtz energy of the solution system. We implemented a method for calculating chemical shifts and applied it to an adenine molecule in water, where hydrogen bonding plays a crucial role in electronic and solvation structures. We also performed the calculations of 17O chemical shifts, which showed remarkable solvent dependence. While converged results could not be sometimes obtained using the conventional method, in the present framework with RISM–SCF–cSED, an adequate representation of electron density is guaranteed, making it possible to obtain an NMR shielding constant stably. This introduction of cSED is key to extending the method’s applicability to obtain the chemical shift of various chemical species. The present demonstration illustrates our approach’s superiority in terms of numerical robustness and accuracy.

Unusually Large Effects of Charge-assisted C-H⋅⋅⋅F Hydrogen Bonds to Anionic Fluorine in Organic Solvents: Computational Study of 19 F NMR Shifts versus Thermochemistry

A comparison of computed ¹⁹ F NMR chemical shifts and experiment provides evidence for large specific solvent effects for fluoride-type anions interacting with the σ*(C-H) orbitals in organic solvents like MeCN or CH₂ Cl₂ . We show this for systems ranging from the fluoride ion and the bifluoride ion [FHF]^- to polyhalogen anions [ClF_x ]^- . Discrepancies between computed and experimental shifts when using continuum solvent models like COSMO or force-field-based descriptions like the 3D-RISM-SCF model show specific orbital interactions that require a quantum-mechanical treatment of the solvent molecules. This is confirmed by orbital analyses of the shielding constants, while less negatively charged fluorine atoms (e. g., in [EF₄ ]^- ) do not require such quantum-mechanical treatments to achieve reasonable accuracy. The larger ¹⁹ F solvent shift of fluoride in MeCN compared to water is due to the larger coordination number in the former. These observations are due to unusually strong charge-assisted C-H⋅⋅⋅F^- hydrogen bonds, which manifest beyond some threshold negative natural charge on fluorine of ca. < -0.6 e. The interactions are accompanied by sizable free energies of solvation, in the order F^- ≫[FHF]^- >[ClF₂ ]^- >[ClF₄ ]^- . COSMO-RS solvation free energies tend to moderately underestimate those from the micro-solvated cluster treatment. Red-shifted and intense vibrational C-H stretching bands, potentially accessible in bulk solution, are further spectroscopic finger prints.

Ab initio machine learning of phase space averages

Equilibrium structures determine material properties and biochemical functions. We here propose to machine learn phase space averages, conventionally obtained by ab initio or force-field-based molecular dynamics (MD) or Monte Carlo (MC) simulations. In analogy to ab initio MD, our ab initio machine learning (AIML) model does not require bond topologies and, therefore, enables a general machine learning pathway to obtain ensemble properties throughout the chemical compound space. We demonstrate AIML for predicting Boltzmann averaged structures after training on hundreds of MD trajectories. The AIML output is subsequently used to train machine learning models of free energies of solvation using experimental data and to reach competitive prediction errors (mean absolute error ∼ 0.8 kcal/mol) for out-of-sample molecules-within milliseconds. As such, AIML effectively bypasses the need for MD or MC-based phase space sampling, enabling exploration campaigns of Boltzmann averages throughout the chemical compound space at a much accelerated pace. We contextualize our findings by comparison to state-of-the-art methods resulting in a Pareto plot for the free energy of solvation predictions in terms of accuracy and time.

Reaction Entropies in Solution from Analytical Three-Dimensional Reference Interaction Site Model Derivatives with Application to Redox and Spin-Crossover Processes

An analytical approach to compute the excess entropy of solvation at constant pressure in three-dimensional reference interaction site model (3D-RISM) calculations is presented. It includes the changes in the macroscopic dielectric constant of the solvent upon variation of temperature and density. The approach is exact within the framework of force-field descriptions of the solute and gives reasonable results for self-consistently determined electrostatics as used in the 3D-RISM-self-consistent field approach, particularly for entropy differences. The new method is applied to simple examples of reaction entropies of iron complexes in aqueous solution, for which simple gas-phase calculations and many other approaches give unreliable estimates. For both redox half-reactions and spin-crossover processes, (semi)quantitative agreement with experimental reaction entropies can be achieved out of the box.

Enabling Cysteine-Free Native Chemical Ligation at Challenging Junctions with a Ligation Auxiliary Capable of Base Catalysis

Ligation auxiliaries are used in chemical protein synthesis to extend the scope of native chemical ligation (NCL) beyond cysteine. However, auxiliary-mediated ligations at sterically demanding junctions have been difficult. Often the thioester intermediate formed in the thiol exchange step of NCL accumulates because the subsequent S→N acyl transfer is extremely slow. Here we introduce the 2-mercapto-2-(pyridin-2-yl)ethyl (MPyE) group as the first auxiliary designed to aid the ligation reaction by catalysis. Notably, the MPyE auxiliary provides useful rates even for junctions containing proline or a β-branched amino acid. Quantum chemical calculations suggest that the pyridine nitrogen acts as an intramolecular base in a rate-determining proton transfer step. The auxiliary is prepared in two steps and conveniently introduced by reductive alkylation. Auxiliary cleavage is induced upon treatment with TCEP/morpholine in presence of a Mn^II complex as radical starter. The synthesis of a de novo designed 99mer peptide and an 80 aa long MUC1 peptide demonstrates the usefulness of the MPyE auxiliary.

Biomolecular Simulations with the Three-Dimensional Reference Interaction Site Model with the Kovalenko-Hirata Closure Molecular Solvation Theory

The statistical mechanics-based 3-dimensional reference interaction site model with the Kovalenko-Hirata closure (3D-RISM-KH) molecular solvation theory has proven to be an essential part of a multiscale modeling framework, covering a vast region of molecular simulation techniques. The successful application ranges from the small molecule solvation energy to the bulk phase behavior of polymers, macromolecules, etc. The 3D-RISM-KH successfully predicts and explains the molecular mechanisms of self-assembly and aggregation of proteins and peptides related to neurodegeneration, protein-ligand binding, and structure-function related solvation properties. Upon coupling the 3D-RISM-KH theory with a novel multiple time-step molecular dynamic (MD) of the solute biomolecule stabilized by the optimized isokinetic Nosé-Hoover chain thermostat driven by effective solvation forces obtained from 3D-RISM-KH and extrapolated forward by generalized solvation force extrapolation (GSFE), gigantic outer time-steps up to picoseconds to accurately calculate equilibrium properties were obtained in this new quasidynamics protocol. The multiscale OIN/GSFE/3D-RISM-KH algorithm was implemented in the Amber package and well documented for fully flexible model of alanine dipeptide, miniprotein 1L2Y, and protein G in aqueous solution, with a solvent sampling rate ~150 times faster than a standard MD simulation in explicit water. Further acceleration in computation can be achieved by modifying the extent of solvation layers considered in the calculation, as well as by modifying existing closure relations. This enhanced simulation technique has proven applications in protein-ligand binding energy calculations, ligand/solvent binding site prediction, molecular solvation energy calculations, etc. Applications of the RISM-KH theory in molecular simulation are discussed in this work.