Intermolecular Interaction in the H2S–H2 Complex: Molecular Beam Scattering Experiments and Ab-Inito Calculations

Situating the phosphonated calixarene-cytochrome C association by molecular dynamics simulations

Protein-calixarenes binding plays an increasingly central role in many applications, spanning from molecular recognition to drug delivery strategies and protein inhibition. These ligands obey a specific bio-supramolecular chemistry, which can be revealed by computational approaches, such as molecular dynamics simulations. In this paper, we rely on all-atom, explicit-solvent molecular dynamics simulations to capture the electrostatically driven association of a phosphonated calix-[4]-arene with cytochome-C, which critically relies on surface-exposed paired lysines. Beyond two binding sites identified in direct agreement with the x-ray structure, the association has a larger structural impact on the protein dynamics. Then, our simulations allow a direct comparison to analogous calixarenes, namely, sulfonato, similarly reported as "molecular glue." Our work can contribute to a robust in silico predictive tool to assess binding sites for any given protein of interest for crystallization, with the specificity of a macromolecular cage whose endo/exo orientation plays a role in the binding.

The dawn of hydrogen and halogen bonds and their crucial role in collisional processes probing long-range intermolecular interactions

This perspective review focuses on the results of an internally consistent study developed in the Perugia laboratory, centered on the fundamental interaction components that, at large intermolecular distances, determine the formation of weak intermolecular hydrogen (HB) and halogen (XB) bonds. This investigation exploits old and novel molecular beam scattering experiments involving several gaseous prototypical systems. In particular, we focus on the kinetic energy dependence of the total (elastic + inelastic) integral cross-sections. Of particular interest is the measure of quantum interference patterns in the energy dependence of cross-sections of targeted systems and their shift compared to that of known reference systems. We interpreted these findings as interaction energy stabilization components, such as charge transfer, σ-hole, and polar flattening, that emerge at intermediate separation distance ranges and selectively manifest for specific geometries of collision complexes. Another significant observable we discuss is the absolute value of the cross-section and its dependence on permanent multipole moments of the collisional partners. Specifically, we show how the spontaneous orientation of rotationally cold and polar molecules, due to the electric field gradient associated with the interaction between permanent multipole moments, can significantly modify the magnitude of the total cross-section, even at high values of the impact parameter. We are confident that the present results can help extend the force field formulation to various interacting systems and carry out molecular dynamics simulations under conditions of application interest.

Interaction of the H2S molecule with molecular hydrogen: Ab initio potential energy surface and scattering calculations

The potential energy surface describing the interaction of H₂S with molecular hydrogen has been computed through explicitly correlated coupled cluster calculations including single, double, and (perturbative) triple excitations [CCSD(T)-F12a], with the assumption of fixed molecular geometries. The computed points were fit to an analytical form suitable for time-independent quantum scattering calculations of rotationally inelastic cross sections and rate coefficients. Representative sets of energy dependent state-to-state rotationally inelastic cross sections and rate coefficients are presented and discussed.

Insight into the halogen-bond nature of noble gas-chlorine systems by molecular beam scattering experiments, ab initio calculations and charge displacement analysis

We have carried out molecular-beam scattering experiments and high-level ab initio investigations on the potential energy surfaces of a series of noble-gas-Cl2 adducts. This effort has permitted the construction of a simple, reliable and easily generalizable analytical model potential formulation, which is based on a few physically meaningful parameters of the interacting partners and transparently shows the origin, strength, and stereospecificity of the various interaction components. The results demonstrate quantitatively beyond doubt that the interaction between a noble-gas (Ng) atom - even He - and Cl2 in a collinear configuration is characterized by weak halogen bond (XB) formation, accompanied by charge transfer (CT) from the Ng to chlorine. This characteristic, which stabilizes the adduct, rapidly disappears on going towards the T-shaped configuration, dominated by pure van der Waals (vdW) forces. Similarly, a pure vdW interaction takes place - with no CT component in any configuration - if Cl2 is present in the lowest πg* → σu* excited state, because the change in electron density that accompanies the excitation eliminates the Cl2 polar flattening and σ hole, making the XB interaction inaccessible.

Energy Decomposition Analysis with a Stable Charge-Transfer Term for Interpreting Intermolecular Interactions

Many schemes for decomposing quantum-chemical calculations of intermolecular interaction energies into physically meaningful components can be found in the literature, but the definition of the charge-transfer (CT) contribution has proven particularly vexing to define in a satisfactory way and typically depends strongly on the choice of basis set. This is problematic, especially in cases of dative bonding and for open-shell complexes involving cation radicals, for which one might expect significant CT. Here, we analyze CT interactions predicted by several popular energy decomposition analyses and ultimately recommend the definition afforded by constrained density functional theory (cDFT), as it is scarcely dependent on basis set and provides results that are in accord with chemical intuition in simple cases, and in quantitative agreement with experimental estimates of the CT energy, where available. For open-shell complexes, the cDFT approach affords CT energies that are in line with trends expected based on ionization potentials and electron affinities whereas some other definitions afford unreasonably large CT energies in large-gap systems, which are sometimes artificially offset by underestimation of van der Waals interactions by density functional theory. Our recommended energy decomposition analysis is a composite approach, in which cDFT is used to define the CT component of the interaction energy and symmetry-adapted perturbation theory (SAPT) defines the electrostatic, polarization, Pauli repulsion, and van der Waals contributions. SAPT/cDFT provides a stable and physically motivated energy decomposition that, when combined with a new implementation of open-shell SAPT, can be applied to supramolecular complexes involving molecules, ions, and/or radicals.

Interaction of O2 with CH4, CF4, and CCl4 by Molecular Beam Scattering Experiments and Theoretical Calculations

Gas phase collisions of O2 by CH4, CF4, and CCl4 have been investigated with the molecular beam technique by measuring both the integral cross section value, Q, and its dependence on the collision velocity, v. The adopted experimental conditions have been appropriate to resolve the oscillating "glory" pattern, a quantum interference effect controlled by the features of the intermolecular interaction, for all the three case studies. The analysis of the Q(v) data, performed by adopting a suitable representation of the intermolecular potential function, provided the basic features of the anisotropic potential energy surfaces at intermediate and large separation distances and information on the relative role of the physically relevant types of contributions to the global interaction. The present work demonstrates that while O2-CH4 and O2-CF4 are basically bound through the balance between size (Pauli) repulsion and dispersion attraction, an appreaciable intermolecular bond stabilization by charge transfer is operative in O2-CCl4. Ab initio calculations of the strength of the interaction, coupled with detailed analysis of electronic charge displacement promoted by the formation of the dimer, fully rationalizes the experimental findings. This investigation indicates that the interactions of O2, when averaged over its relative orientations, are similar to that of a noble gas (Ng), specifically Ar. We also show that the binding energy in the basic configurations of the prototypical Ng-CF4,CCl4 systems [ Cappelletti , D. ; Chem. Eur. J. 2015 , 21 , 6234 - 6240 ] can be reconstructed by using the interactions in Ng-F and Ng-Cl systems, previously characterized by molecular beam scattering experiments of state-selected halogen atom beams. This information is fundamental to approach the modeling of the weak intermolecular halogen bond. On the basis of the electronic polarizability, this also confirms [ Aquilanti , V. ; Angew. Chem., Int. Ed. 2005 , 44 , 2356 - 2360 ] that O2 can be taken as a proper reference partner for simulating the behavior of some basic noncovalent components of the interactions involving water. Present results are of fundamental relevance to build up the force field controlling the hydrophobic behavior of prototypical apolar CX4 (X = H, F, Cl) molecules.

H2O–CH4and H2S–CH4complexes: a direct comparison through molecular beam experiments and ab initio calculations

Electron density redistribution upon the formation of the water–methane complex arises from polarisation and charge transfer effects.