Factors Impacting σ- and π-Hole Regions as Revealed by the Electrostatic Potential and Its Source Function Reconstruction: The Case of 4,4'-Bipyridine Derivatives

Recent developments in molecular modeling tools and applications related to pharmaceutical and biomedical research

In modern pharmaceutical and biomedical research, molecular modeling represents a useful tool to explore processes and their mechanistic bases at the molecular level. Integrating experimental and virtual analysis is a fruitful approach to study ligand-receptor interaction in chemical, biochemical and biological environments. In these fields, molecular docking and molecular dynamics are considered privileged techniques for modeling (bio)macromolecules and related complexes. This review aims to present the current landscape of molecular modeling in pharmaceutical and biomedical research by examining selected representative applications published in the last years and highlighting current topics and trends of this field. Thus, a systematic compilation of all published literature has not been attempted herein. After a brief overview of the main theoretical and computational tools used to investigate mechanisms at molecular level, recent applications of molecular modeling in drug discovery, ligand binding and for studying protein conformation and function will be discussed. Furthermore, specific sections will be devoted to the application of molecular modeling for unravelling enantioselective mechanisms underlying the enantioseparation of chiral compounds of pharmaceutical and biomedical interest as well as for studying new forms of noncovalent interactivity identified in biochemical and biological environments. The general aim of this review is to provide the reader with a modern overview of the topic, highlighting advancements and outlooks as well as drawbacks and pitfalls still affecting the applicability of theoretical and computational methods in the field of pharmaceutical and biomedical research.

Unravelling dispersion forces in liquid-phase enantioseparation. Part I: Impact of ferrocenyl versus phenyl groups

BACKGROUND

Highly ordered chiral secondary structures as well as multiple (tunable) recognition sites are the keys to success of polysaccharide carbamate-based chiral selectors in enantioseparation science. Hydrogen bonds (HBs), dipole-dipole, and π-π interactions are classically considered the most frequent noncovalent interactions underlying enantioselective recognition with these chiral selectors. Very recently, halogen, chalcogen and π-hole bonds were also identified as interactions working in polysaccharide carbamate-based selectors to promote enantiomer distinction. On the contrary, the function of dispersion interactions in this field was not explored so far.

RESULTS

The enantioseparation of chiral ferrocenes featuring chiral axis or chiral plane as stereogenic elements was performed by comparing five polysaccharide carbamate-based chiral columns, with the aim to identify enantioseparation outcomes that could be reasonably determined by dispersion forces, making available a reliable experimental data set for future theoretical studies to confirm the heuristic hypothesis. The effects of mobile phase polarity and temperature on the enantioseparation were considered, and potential recognition sites on analytes and selectors were evaluated by electrostatic potential (V) analysis and molecular dynamics (MD). In this first part, the enantioseparation of 3,3'-dibromo-5,5'-bis-ferrocenylethynyl-4,4'-bipyridine bearing two ferrocenylethynyl units linked to an axially chiral core was performed and compared to that of the analyte featuring the same structural motif with two phenyl groups in place of the ferrocenyl moieties. The results of this study showed the superiority of the ferrocenyl compared to the phenyl group, as a structural element favouring enantiodifferentiation.

SIGNIFICANCE AND NOVELTY

Even if dispersion (London) forces have been envisaged acting in liquid-phase enantioseparations, focused studies to explore possible contributions of dispersion forces with polysaccharide carbamate-based selectors are practically missing. This study allowed us to collect experimental information that support the involvement of dispersion forces as contributors to liquid-phase enantioseparation, paving the way to a new picture in this field.

Methylammonium Tetrel Halide Perovskite Ion Pairs and Their Dimers: The Interplay between the Hydrogen-, Pnictogen- and Tetrel-Bonding Interactions

The structural stability of the extensively studied organic-inorganic hybrid methylammonium tetrel halide perovskite semiconductors, MATtX₃ (MA = CH₃NH₃⁺; Tt = Ge, Sn, Pb; X = Cl, Br, I), arises as a result of non-covalent interactions between an organic cation (CH₃NH₃⁺) and an inorganic anion (TtX₃^-). However, the basic understanding of the underlying chemical bonding interactions in these systems that link the ionic moieties together in complex configurations is still limited. In this study, ion pair models constituting the organic and inorganic ions were regarded as the repeating units of periodic crystal systems and density functional theory simulations were performed to elucidate the nature of the non-covalent interactions between them. It is demonstrated that not only the charge-assisted N-H···X and C-H···X hydrogen bonds but also the C-N···X pnictogen bonds interact to stabilize the ion pairs and to define their geometries in the gas phase. Similar interactions are also responsible for the formation of crystalline MATtX₃ in the low-temperature phase, some of which have been delineated in previous studies. In contrast, the Tt···X tetrel bonding interactions, which are hidden as coordinate bonds in the crystals, play a vital role in holding the inorganic anionic moieties (TtX₃^-) together. We have demonstrated that each Tt in each [CH₃NH₃⁺•TtX₃^-] ion pair has the capacity to donate three tetrel (σ-hole) bonds to the halides of three nearest neighbor TtX₃^- units, thus causing the emergence of an infinite array of 3D TtX₆^4- octahedra in the crystalline phase. The TtX₄^4- octahedra are corner-shared to form cage-like inorganic frameworks that host the organic cation, leading to the formation of functional tetrel halide perovskite materials that have outstanding optoelectronic properties in the solid state. We harnessed the results using the quantum theory of atoms in molecules, natural bond orbital, molecular electrostatic surface potential and independent gradient models to validate these conclusions.

Recognition in the Domain of Molecular Chirality: From Noncovalent Interactions to Separation of Enantiomers

It is not a coincidence that both chirality and noncovalent interactions are ubiquitous in nature and synthetic molecular systems. Noncovalent interactivity between chiral molecules underlies enantioselective recognition as a fundamental phenomenon regulating life and human activities. Thus, noncovalent interactions represent the narrative thread of a fascinating story which goes across several disciplines of medical, chemical, physical, biological, and other natural sciences. This review has been conceived with the awareness that a modern attitude toward molecular chirality and its consequences needs to be founded on multidisciplinary approaches to disclose the molecular basis of essential enantioselective phenomena in the domain of chemical, physical, and life sciences. With the primary aim of discussing this topic in an integrated way, a comprehensive pool of rational and systematic multidisciplinary information is provided, which concerns the fundamentals of chirality, a description of noncovalent interactions, and their implications in enantioselective processes occurring in different contexts. A specific focus is devoted to enantioselection in chromatography and electromigration techniques because of their unique feature as "multistep" processes. A second motivation for writing this review is to make a clear statement about the state of the art, the tools we have at our disposal, and what is still missing to fully understand the mechanisms underlying enantioselective recognition.

Stereoselective Processes Based on σ-Hole Interactions

The σ-hole interaction represents a noncovalent interaction between atoms with σ-hole(s) on their surface (such as halogens and chalcogens) and negative sites. Over the last decade, significant developments have emerged in applications where the σ-hole interaction was demonstrated to play a key role in the control over chirality. The aim of this review is to give a comprehensive overview of the current advancements in the use of σ-hole interactions in stereoselective processes, such as formation of chiral supramolecular assemblies, separation of enantiomers, enantioselective complexation and asymmetric catalysis.

Anatomy of π-hole bonds: Linear systems

The list of σ-hole bonds is long and growing, encompassing both H-bonds and its closely related halogen, chalcogen, etc., sisters. These bonds rely on the asymmetric distribution of electron density, whose depletion along the extension of a covalent bond leaves a positive region of electrostatic potential from which these bonds derive their name. However, the density distributions of other molecules contain analogous positive regions that lie out of the molecular plane known as π-holes, which are likewise capable of engaging in noncovalent bonds. Quantum calculations are applied to study such π-hole bonds that involve linear molecules, whose positive region is a circular belt surrounding the molecule, rather than the more restricted area of a σ-hole. These bonds are examined in terms of their most fundamental elements arising from the spatial dispositions of their relevant molecular orbitals and the π-holes in both the total electron density and the electrostatic potential to which they lead. Systems examined comprise tetrel, chalcogen, aerogen, and triel bonds, as well as those involving group II elements, with atoms drawn from various rows of the Periodic Table. The π-hole bonds established by linear molecules tend to be weaker than those of comparable planar systems.

Dissection of the Origin of π-Holes and the Noncovalent Bonds in Which They Engage

Accompanying the rapidly growing list of σ-hole bonds has come the acknowledgment of parallel sorts of noncovalent bonds which owe their stability in large part to a deficiency of electron density in the area above the molecular plane, known as a π-hole. The origins of these π-holes are probed for a wide series of molecules, comprising halogen, chalcogen, pnicogen, tetrel, aerogen, and spodium bonds. Much like in the case of their σ-hole counterparts, formation of the internal covalent π-bond in the Lewis acid molecule pulls density toward the bond midpoint and away from its extremities. This depletion of density above the central atom is amplified by an electron-withdrawing substituent. At the same time, the amplitude of the π*-orbital is enhanced in the region of the density-depleted π-hole, facilitating a better overlap with the nucleophile's lone pair orbital and a stabilizing n → π* charge transfer. The presence of lone pairs on the central atom acts to attenuate the π-hole and shift its position somewhat, resulting in an overall weakening of the π-hole bond. There is a tendency for π-hole bonds to include a higher fraction of induction energy than σ-bonds with proportionately smaller electrostatic and dispersion components, but this distinction is less a product of the σ- or π-character and more a function of the overall bond strength.

Enantioseparation of 5,5'-Dibromo-2,2'-Dichloro-3-Selanyl-4,4'-Bipyridines on Polysaccharide-Based Chiral Stationary Phases: Exploring Chalcogen Bonds in Liquid-Phase Chromatography

The chalcogen bond (ChB) is a noncovalent interaction based on electrophilic features of regions of electron charge density depletion (σ-holes) located on bound atoms of group VI. The σ-holes of sulfur and heavy chalcogen atoms (Se, Te) (donors) can interact through their positive electrostatic potential (V) with nucleophilic partners such as lone pairs, π-clouds, and anions (acceptors). In the last few years, promising applications of ChBs in catalysis, crystal engineering, molecular biology, and supramolecular chemistry have been reported. Recently, we explored the high-performance liquid chromatography (HPLC) enantioseparation of fluorinated 3-arylthio-4,4′-bipyridines containing sulfur atoms as ChB donors. Following this study, herein we describe the comparative enantioseparation of three 5,5′-dibromo-2,2′-dichloro-3-selanyl-4,4′-bipyridines on polysaccharide-based chiral stationary phases (CSPs) aiming to understand function and potentialities of selenium σ-holes in the enantiodiscrimination process. The impact of the chalcogen substituent on enantioseparation was explored by using sulfur and non-chalcogen derivatives as reference substances for comparison. Our investigation also focused on the function of the perfluorinated aromatic ring as a π-hole donor recognition site. Thermodynamic quantities associated with the enantioseparation were derived from van’t Hoff plots and local electron charge density of specific molecular regions of the interacting partners were inspected in terms of calculated V. On this basis, by correlating theoretical data and experimental results, the participation of ChBs and π-hole bonds in the enantiodiscrimination process was reasonably confirmed.