Quantifying Interactions and Solvent Effects Using Molecular Balances and Model Complexes

Molecular Ratchets and Kinetic Asymmetry: Giving Chemistry Direction

Over the last two decades ratchet mechanisms have transformed the understanding and design of stochastic molecular systems-biological, chemical and physical-in a move away from the mechanical macroscopic analogies that dominated thinking regarding molecular dynamics in the 1990s and early 2000s (e.g. pistons, springs, etc), to the more scale-relevant concepts that underpin out-of-equilibrium research in the molecular sciences today. Ratcheting has established molecular nanotechnology as a research frontier for energy transduction and metabolism, and has enabled the reverse engineering of biomolecular machinery, delivering insights into how molecules 'walk' and track-based synthesisers operate, how the acceleration of chemical reactions enables energy to be transduced by catalysts (both motor proteins and synthetic catalysts), and how dynamic systems can be driven away from equilibrium through catalysis. The recognition of molecular ratchet mechanisms in biology, and their invention in synthetic systems, is proving significant in areas as diverse as supramolecular chemistry, systems chemistry, dynamic covalent chemistry, DNA nanotechnology, polymer and materials science, molecular biology, heterogeneous catalysis, endergonic synthesis, the origin of life, and many other branches of chemical science. Put simply, ratchet mechanisms give chemistry direction. Kinetic asymmetry, the key feature of ratcheting, is the dynamic counterpart of structural asymmetry (i.e. chirality). Given the ubiquity of ratchet mechanisms in endergonic chemical processes in biology, and their significance for behaviour and function from systems to synthesis, it is surely just as fundamentally important. This Review charts the recognition, invention and development of molecular ratchets, focussing particularly on the role for which they were originally envisaged in chemistry, as design elements for molecular machinery. Different kinetically asymmetric systems are compared, and the consequences of their dynamic behaviour discussed. These archetypal examples demonstrate how chemical systems can be driven inexorably away from equilibrium, rather than relax towards it.

Aromatic-Carbonyl Interactions as an Emerging Type of Non-Covalent Interactions

Aromatic-carbonyl (Ar···C═O) interactions, attractive interactions between the arene plane and the carbon atom of carbonyl, are in the infancy as one type of new supramolecular bonding forces. Here the study and functionalization of aromatic-carbonyl interactions in solution is reported. A combination of aromatic-carbonyl interactions and dynamic covalent chemistry provided a versatile avenue. The stabilizing role and mechanism of arene-aldehyde/imine interactions are elucidated through crystal structures, NMR studies, and computational evidence. The movement of imine exchange equilibria further allowed the quantification of the interplay between arene-aldehyde/imine interactions and dynamic imine chemistry, with solvent effects offering another handle and matching the electrostatic feature of the interactions. Moreover, arene-aldehyde/imine interactions enabled the reversal of kinetic and thermodynamic selectivity and sorting of dynamic covalent libraries. To show the functional utility diverse modulation of fluorescence signals is realized with arene-aldehyde/imine interactions. The results should find applications in many aspects, including molecular recognition, assemblies, catalysis, and intelligent materials.

Distinction and Quantification of Noncovalent Dispersive and Hydrophobic Effects

The possibilities of comparing computational results of noncovalent interactions with experimental data are discussed, first with respect to intramolecular interactions. For these a variety of experimental data such as heats of formation, crystal sublimation heats, comparison with energy minimized structures, and spectroscopic data are available, but until now largely have not found widespread application. Early force field and QM/MP2 calculations have already shown that the sublimation heats of hydrocarbons can be predicted with an accuracy of ±1%. Intermolecular interactions in solution or the gas phase are always accompanied by difficult to compute entropic contributions, like all associations between molecules. Experimentally observed T∆S values contribute 10% to 80% of the total ∆G, depending on interaction mechanisms within the complexes, such as, e.g., hydrogen bonding and ion pairing. Free energies ∆G derived from equilibrium measurements in solution allow us to define binding increments ∆∆G, which are additive and transferable to a variety of supramolecular complexes. Data from more than 90 equilibrium measurements of porphyrin receptors in water indicate that small alkanes do not bind to the hydrophobic flat surfaces within a measuring limit of ∆G = ±0.5 kJ/mol, and that 20 functions bearing heteroatoms show associations by dispersive interactions with up to ∆G = 8 kJ/mol, roughly as a function of their polarizability. Aromatic systems display size-dependent affinities ∆G as a linear function of the number of π-electrons.

Advances and Prospects in Understanding London Dispersion Interactions in Molecular Chemistry

London dispersion (LD) interactions are the main contribution of the attractive part of the van der Waals potential. Even though LD effects are the driving force for molecular aggregation and recognition, the role of these omnipresent interactions in structure and reactivity had been largely underappreciated over decades. However, in the recent years considerable efforts have been made to thoroughly study LD interactions and their potential as a chemical design element for structures and catalysis. This was made possible through a fruitful interplay of theory and experiment. This review highlights recent results and advances in utilizing LD interactions as a structural motif to understand and utilize intra- and intermolecularly LD-stabilized systems. Additionally, we focus on the quantification of LD interactions and their fundamental role in chemical reactions.

Probing Non-Covalent Interactions through Molecular Balances: A REG-IQA Study

The interaction energies of two series of molecular balances (1-X with X = H, Me, OMe, NMe2 and 2-Y with Y = H, CN, NO2, OMe, NMe2) designed to probe carbonyl…carbonyl interactions were analysed at the B3LYP/6-311++G(d,p)-D3 level of theory using the energy partitioning method of Interacting Quantum Atoms/Fragments (IQA/IQF). The partitioned energies are analysed by the Relative Energy Gradient (REG) method, which calculates the correlation between these energies and the total energy of a system, thereby explaining the role atoms have in the energetic behaviour of the total system. The traditional "back-of-the-envelope" open and closed conformations of molecular balances do not correspond to those of the lowest energy. Hence, more care needs to be taken when considering which geometries to use for comparison with the experiment. The REG-IQA method shows that the 1-H and 1-OMe balances behave differently to the 1-Me and 1-NMe2 balances because the latter show more prominent electrostatics between carbonyl groups and undergoes a larger dihedral rotation due to the bulkiness of the functional groups. For the 2-Y balance, REG-IQA shows the same behaviour across the series as the 1-H and 1-OMe balances. From an atomistic point of view, the formation of the closed conformer is favoured by polarisation and charge-transfer effects on the amide bond across all balances and is counterbalanced by a de-pyramidalisation of the amide nitrogen. Moreover, focusing on the oxygen of the amide carbonyl and the α-carbon of the remaining carbonyl group, electrostatics have a major role in the formation of the closed conformer, which goes against the well-known n-π* interaction orbital overlap concept. However, REG-IQF shows that exchange-correlation energies overtake electrostatics for all the 2-Y balances when working with fragments around the carbonyl groups, while they act on par with electrostatics for the 1-OMe and 1-NMe2. REG-IQF also shows that exchange-correlation energies in the 2-Y balance are correlated to the inductive electron-donating and -withdrawing trends on aromatic groups. We demonstrate that methods such as REG-IQA/IQF can help with the fine-tuning of molecular balances prior to the experiment and that the energies that govern the probed interactions are highly dependent on the atoms and functional groups involved.

Illuminating a Solvent-Dependent Hierarchy for Aromatic CH/π Complexes with Dynamic Covalent Glyco-Balances

CH/π interactions are prevalent among aromatic complexes and represent invaluable tools for stabilizing well-defined molecular architectures. Their energy contributions are exceptionally sensitive to various structural and environmental factors, resulting in a context-dependent nature that has led to conflicting findings in the scientific literature. Consequently, a universally accepted hierarchy for aromatic CH/π interactions has remained elusive. Herein, we present a comprehensive experimental investigation of aromatic CH/π complexes, employing a novel approach that involves isotopically labeled glyco-balances generated in situ. This innovative strategy not only allows us to uncover thermodynamic insights but also delves into the often less-accessible domain of kinetic information. Our analyses have yielded more than 180 new free energy values while considering key factors such as solvent properties, the interaction geometry, and the presence and nature of accompanying counterions. Remarkably, the obtained results challenge conventional wisdom regarding the stability order of common aromatic complexes. While it was believed that cationic CH/π interactions held the highest strength, followed by polarized CH/π, nonpolarized CH/π, and finally anionic CH/π interactions, our study reveals that this hierarchy can be subverted depending on the environment. Indeed, the performance of polarized CH/π interactions can match or even outcompete that of cationic CH/π interactions making them a more reliable stabilization strategy across the entire spectrum of solvent polarity. Overall, our results provide valuable guidelines for the selection of optimal interacting partners in every chemical environment, allowing the design of tailored aromatic complexes with applications in supramolecular chemistry, organocatalysis, and/or material sciences.

Dispersion Interactions in Condensed Phases and inside Molecular Containers

ConspectusThe past decade has seen significant progress in the understanding and appreciation of the importance of London dispersion interactions (LDIs) in supramolecular systems and solutions. The Slater-Kirkwood formula relates LDIs to the molecular polarizabilities of the two interacting molecular species (α) and their interaction distance (a dependence of R-6). When advancing arguments related to intermolecular interactions, it is frequently assumed that molecules with larger molecular polarizabilities are more amenable to larger LDIs. However, arguments related to molecular polarizabilities are not always transferable to the condensed phase. In fact, the underlying bulk and molecular polarizabilities of common solvents show opposing trends. The intuitive concept that aromatic molecules are more polarizable than saturated hydrocarbons and that perfluorinated molecules are less polarizable than saturated hydrocarbons applies to the condensed phase only. When treating association phenomena in solution, where LDIs are generally very attenuated, the use of bulk polarizabilities is recommended, which are experimentally accessible through either refractive index measurements or suitable solvatochromic probes. Such probes can also be used to assess polarizabilities inside molecular container compounds, such as cucurbit[n]urils (CBn), cyclodextrins, calixarenes, and hemicarcerands. These macrocyclic cavities can have extreme microenvironments. For example, the inner concave phase of CB7 has been shown to be weakly polarizable, falling in between the gas phase and perfluorohexane; those of β-cyclodextrin and p-sulfonatocalix[4]arene have been found to be similarly polarizable as water and alkanes, respectively, and the inside of hemicarcerands displays a very large bulk polarizability, exceeding that of diiodomethane. CBn compounds are privileged molecular container compounds, which we exemplify in this Account through case studies. (1) CBn macrocycles are prime water-soluble receptors for hydrocarbons, allowing for the reduction of the binding free energies to two components: the hydrophobic effect and dispersion interactions. To understand hydrocarbon binding, we initiated the HYDROPHOBE challenge, which revealed the shortcomings of both quantum-chemical and molecular dynamics approaches. (2) The smallest CBn receptor, CB5, is uniquely suited to bind the entire noble gas series, where hydrophobic effects and dispersion interactions operate in opposite directions. CB5 was revaled to be a unique synthetic receptor for noble gases, with the dominant driving force being the recovery of the cavitation energies for the hydration of noble gases in aqueous solution. Computational methods that encounter challenges in predicting hydrocarbon affinities and trends for CB6 and CB7 perform well for noble gases binding to CB5. (3) The larger homologue, CB8, allows one to set up intermolecular interaction chambers by the encapsulation of a (first) aromatic guest, thereby tuning LDIs inside the receptor cavity. In this manner, CB8 can be modulated to preferentially bind unsaturated and aromatic rather than saturated hydrocarbons, while the unmodified cavities of the smaller macrocycles CB6 and CB7 show selective binding of saturated hydrocarbons. (4) The (charged) host-guest complexes of CBn hosts are sufficiently stable in the gas phase, allowing for the study of the influence of LDIs on inner-phase chemical reactions. These studies are particularly interesting for the theoretical analysis of isolated host-guest LDIs, as experimental and computational data are directly comparable in the gas phase due to the absence of the solvation effect.

The Monomer-Dimer Equilibrium of Triscatechol Titanium(IV)-Based Hierarchical Helicates as a Tool for the Development of Molecular Balances and Molecular Switches

ConspectusHierarchical helicates are formed by noncovalent connection of two or more monomeric metal complex units, e.g., by bridging metal cations. A unique kind of hierarchical helicate is obtained from 3-carbonyl substituted catechol ligands with titanium(IV) ions in the presence of lithium cations. This kind of supramolecular complex shows in solution a "monomer-dimer" equilibrium. There are different possibilities (solvent, countercation, substituents at carbonyl unit, etc.) to shift this equilibrium to either the monomer or the dimer side. Thus, the lithium-bridged catecholate-based hierarchical helicates resemble a molecular switch. In this Account, different aspects are discussed of how this unique behavior of the dimeric titanium catecholates can be used for application.Thorough investigation of the energetics of the monomer-dimer equilibrium leads to a deeper understanding of the thermodynamic and kinetic effects of the dimerization (or dissociation) process. In this context, even weak interaction of substituents in the periphery of the complexes can be observed. Hereby on the one hand, solvent effects have an important influence and can be easily evaluated. The thorough understanding of the behavior of the monomer-dimer equilibrium allows one to develop some novel applications. In this respect, the use of the hierarchical helicate-based switch as a platform for reaction control and catalysis is described. Decent enantioselectivities up to ee = 58% can be found in Diels-Alder reactions in the periphery of the dimers, while switching to the monomer as a reaction platform still allows the cycloaddition reaction but turns the selectivity off. Additionally, it is described that catalytically important units can be introduced and hydrogenation reactions as well as Michael-type reactions are catalyzed at the helicates.Covalent connection of two catechol ester units leads to classical helicates. Depending on the alkaline metal cation, those can be switched from a compressed to an expanded form or vice versa. Hereby the monomer-dimer equilibrium is transformed into a structural switch. The switching process can be initiated by removal or addition of lithium cations (e.g., by addition of [2.1.1]cryptand). Alternative switching possibilities are based in the case of glycol bridged helicates on cation translocation isomerism and with thioester derivatives it occurs spontaneously in DMSO. Introduction of chiral tethers results in a three state switch allowing expansion/compression as well as switching of the helicity.

Experimental Quantification of Halogen⋅⋅⋅Arene van der Waals Contacts

Crystallographic and computational studies suggest the occurrence of favourable interactions between polarizable arenes and halogen atoms. However, the systematic experimental quantification of halogen⋅⋅⋅arene interactions in solution has been hindered by the large variance in the steric demands of the halogens. Here we have synthesized molecular balances to quantify halogen⋅⋅⋅arene contacts in 17 solvents and solvent mixtures using 1 H NMR spectroscopy. Calculations indicate that favourable halogen⋅⋅⋅arene interactions arise from London dispersion in the gas phase. In contrast, comparison of our experimental measurements with partitioned SAPT0 energies indicate that dispersion is sufficiently attenuated by the solvent that the halogen⋅⋅⋅arene interaction trend was instead aligned with increasing exchange repulsion as the halogen increased in size (ΔGX ⋅⋅⋅Ph =0 to +1.5 kJ mol-1 ). Halogen⋅⋅⋅arene contacts were slightly less disfavoured in solvents with higher solvophobicities and lower polarizabilities, but strikingly, were always less favoured than CH3 ⋅⋅⋅arene contacts (ΔGMe ⋅⋅⋅Ph =0 to -1.4 kJ mol-1 ).

Systematic Study of Heteroarene Stacking Using a Congeneric Set of Molecular Glues for Procaspase-6

Heteroaromatic stacking interactions are important in drug binding, supramolecular chemistry, and materials science, making protein-ligand model systems of these interactions of considerable interest. Here we studied 30 congeneric ligands that each present a distinct heteroarene for stacking between tyrosine residues at the dimer interface of procaspase-6. Complex X-ray crystal structures of 10 analogs showed that stacking geometries were well conserved, while high-accuracy computations showed that heteroarene stacking energy was well correlated with predicted overall ligand binding energies. Empirically determined KD values in this system thus provide a useful measure of heteroarene stacking with tyrosine. Stacking energies are discussed in the context of torsional strain, the number and positioning of heteroatoms, tautomeric state, and coaxial orientation of heteroarene in the stack. Overall, this study provides an extensive data set of empirical and high-level computed binding energies in a versatile new protein-ligand system amenable to studies of other intermolecular interactions.

Chalcogen Bonds: How to Characterize Them in Solution?

Chalcogen bonds (ChBs) occur between molecules containing Lewis acidic chalcogen substituents and Lewis bases. Recently, ChB emerged as a pivotal interaction in solution-based applications such as anion recognition, anion transport and catalysis. However, before moving to applications, the involvement of ChB must be established in solution. In this Concept article, we provide a brief review of the currently available experimental investigations of ChB in solution.

Probing the Size Limit of Dispersion Energy Donors with a Bifluorenylidene Balance: Magic Cyclohexyl

We report the synthesis of 14 2,2'-disubstituted 9,9'-bifluorenylidenes as molecular balances for the quantification of London dispersion interactions between various dispersion energy donors. For all balances, we measured ΔG_Z/E at 333 K using ¹H NMR in seven organic solvents. For various alkyl and aryl substituents, we generally observe a preference for the "folded" Z-isomer due to attractive London dispersion interactions. The cyclohexyl-substituted system shows the largest Z-preference in this study with ΔG_Z/E = -0.6 ± 0.05 kcal mol^-1 in all solvents, owing to the rotational freedom of cyclohexyl groups paired with their large polarizability that maximizes London dispersion interactions. On the other hand, rigid and sterically more demanding substituents like tert-butyl unexpectedly favor the unfolded E-isomer. This is a result of the close relative position in which the functional groups are positioned in this molecular balance. This close proximity is the reason for the increase of Pauli repulsion in the Z-isomers with large rigid substituents (tert-butyl, adamantyl, and diamantyl) which leads to an equilibrium shift toward the unfolded E-form. While we were able to reproduce most of our experimental trends qualitatively using contemporary computational chemistry methods, quantitative accuracy of the employed methods still needs further improvement.

Assessing the Experimental Hydrogen Bonding Energy of the Cyclic Water Dimer Transition State with a Cyclooctatetraene-Based Molecular Balance

We have conducted an experimental and computational study of cyclooctatetraene-1,4/1,6-dimethanol (1,4 and 1,6) as a molecular balance with the goal in mind to determine the otherwise inaccessible hydrogen bonding energy (HBE) of the cyclic water dimer, which constitutes a transition state. The 1,4/1,6 folding equilibrium is governed by an intramolecular hydrogen bond in the folded 1,6-isomer, in which the OH groups adopt a cyclic planar geometry, akin to the structure of the cyclic water dimer transition state. We characterized hydrogen bonding in 1,6 and reference complexes utilizing SAPT2 + (3)δMP2/aug-cc-pVTZ and selected quantum theory of atoms in molecule descriptors at M06-2XD3(0)/ma-def2-TZVPP. Additionally, we computed HBEs at the DLPNO-CCSD(T)/aug-cc-pVQZ level of theory. We find that hydrogen bonding in 1,6 is very similar to the interaction in the C_i symmetric cyclic water dimer TS, both in magnitude and character. We experimentally determined the Gibbs free energy of the folding process (ΔG_eq) in a variety of organic solvents via nuclear magnetic resonance spectroscopy measurements at room temperature. By combining experimentally obtained ΔG_eq values with corrections derived from accurate computational methods, we provide estimates for the HBE of cyclic water dimers and the cyclic water dimer TS, as the most stable cyclic water dimer.

Dissecting Solvent Effects on Hydrogen Bonding

The experimental isolation of H‐bond energetics from the typically dominant influence of the solvent remains challenging. Here we use synthetic molecular balances to quantify amine/amide H‐bonds in competitive solvents. Over 200 conformational free energy differences were determined using 24 H‐bonding balances in 9 solvents spanning a wide polarity range. The correlations between experimental interaction energies and gas‐phase computed energies exhibited wild solvent‐dependent variation. However, excellent correlations were found between the same computed energies and the experimental data following empirical dissection of solvent effects using Hunter's α/β solvation model. In addition to facilitating the direct comparison of experimental and computational data, changes in the fitted donor and acceptor constants reveal the energetics of secondary local interactions such as competing H‐bonds.

Amphiphilic Polymers for Color Dispersion: Toward Stable and Low-Viscosity Inkjet Ink

Amphiphilic random and block copolymers were synthesized as potential inkjet inks. This study evaluated the potential of these polymers for color dispersion by examining the following factors: surface tension, zeta potential, viscosity, and particle size. Acrylic acid and (ethoxyethoxy)ethyl acrylate were used as the hydrophilic molecular units. Styrene, butyl acrylate, and phenoxyethyl acrylate were used as hydrophobic units. Color dispersions were prepared by using organic dye and these amphiphilic polymers. The color dispersions containing random copolymers exhibited low viscosity, which is preferable for jetting, but the dye particles tended to sediment after the thermal aging test. In contrast, those containing block copolymers showed high viscosity, which was unsuitable for jetting. However, they retained their initial dispersion state after the aging test. The advantages and disadvantages of each monomer arrangement (random or block) were demonstrated, providing a future outlook on the molecular design of polymer dispersants for color dispersions.

Characterizing the Properties of Anion-Binding Bis(cyclopeptides) with Solvent-Independent Energy Increments

The binding energies of 121 complexes between anions and bis(cyclopeptides) differing in the structure and the number of linking units between the two cyclopeptide rings were analyzed. These Gibbs free energies were obtained in earlier work for different anions, under different conditions, and with different methods. The multiparametric analysis of a subset of 42 binding energies afforded linear relationships that allowed the relatively reliable estimation of the iodide and sulfate affinity of three structurally related bis(cyclopeptides) in water/methanol and water/acetonitrile mixtures at different solvent compositions. Three parameters were required to achieve a satisfactory correlation, namely, the Gibbs free energy of transferring the respective anion from water into the solvent mixture in which complex stability was determined, and the Kamlet–Taft parameters α and β. Based on these relationships, the anion affinities of the other bis(cyclopeptides) were evaluated, giving rise to a set of energy increments that allow quantifying the effects of the linker structure or the nature of the anion on binding affinity relative to the reference system.

Assessing the persistence of chalcogen bonds in solution with neural network potentials

Non-covalent bonding patterns are commonly harvested as a design principle in the field of catalysis, supramolecular chemistry and functional materials to name a few. Yet, their computational description generally neglects finite temperature and environment effects, which promote competing interactions and alter their static gas-phase properties. Recently, neural network potentials (NNPs) trained on Density Functional Theory (DFT) data have become increasingly popular to simulate molecular phenomena in condensed phase with an accuracy comparable to ab initio methods. To date, most applications have centered on solid-state materials or fairly simple molecules made of a limited number of elements. Herein, we focus on the persistence and strength of chalcogen bonds involving a benzotelluradiazole in condensed phase. While the tellurium-containing heteroaromatic molecules are known to exhibit pronounced interactions with anions and lone pairs of different atoms, the relevance of competing intermolecular interactions, notably with the solvent, is complicated to monitor experimentally but also challenging to model at an accurate electronic structure level. Here, we train direct and baselined NNPs to reproduce hybrid DFT energies and forces in order to identify what are the most prevalent non-covalent interactions occurring in a solute-Cl$^-$-THF mixture. The simulations in explicit solvent highlight competition with chalcogen bonds formed with the solvent and the short-range directionality of the interaction with direct consequences for the molecular properties in the solution. The comparison with other potentials (e.g., AMOEBA, direct NNP and continuum solvent model) also demonstrates that baselined NNPs offer a reliable picture of the non-covalent interaction interplay occurring in solution.

An Incremental System To Predict the Effect of Different London Dispersion Donors in All‐ meta ‐Substituted Azobenzenes

Predictive models based on incremental systems exist for many chemical phenomena, thus allowing easy estimates. Despite their low magnitude in isolated systems London dispersion interactions are ubiquitous in manifold situations ranging from solvation to catalysis or in biological systems. Based on our azobenzene system, we systematically determined the London dispersion donor strength of the alkyl substituents Me, Et, iPr up to tBu. Based on this data, we were able to implement an incremental system for London dispersion for the azobenzene scheme. We propose an equation that allows the prediction of the effect of change of substituents on London dispersion interactions in azobenzenes, which has to be validated in similar molecular arrangements in the future.

Quantifying Solvophobic Effects in Organic Solvents Using a Hydrocarbon Molecular Balance

We evaluate the use of the cohesive energy density (ced) as a quantitative descriptor for solvophobic effects in organic solvents by measuring ΔG_Z/E of the rigid Z- and E-2,2'-diethynyl-9,9'-bifluorenylidene. In line with previously employed balances, solvent-dependent changes in ΔG_Z/E are predominantly induced by solvophobic effects, leading to a strong correlation with the solvent's ced. We re-emphasize the role of ceds as quantitative descriptors of solvophobic effects of organic solvents. Our experimental findings are well supported by B3LYP-D3/def2TZVP computations.

Quantum Chemical Microsolvation by Automated Water Placement

We developed a quantitative approach to quantum chemical microsolvation. Key in our methodology is the automatic placement of individual solvent molecules based on the free energy solvation thermodynamics derived from molecular dynamics (MD) simulations and grid inhomogeneous solvation theory (GIST). This protocol enabled us to rigorously define the number, position, and orientation of individual solvent molecules and to determine their interaction with the solute based on physical quantities. The generated solute-solvent clusters served as an input for subsequent quantum chemical investigations. We showcased the applicability, scope, and limitations of this computational approach for a number of small molecules, including urea, 2-aminobenzothiazole, (+)-syn-benzotriborneol, benzoic acid, and helicene. Our results show excellent agreement with the available ab initio molecular dynamics data and experimental results.