Halogen Bond Structure and Dynamics from Molecular Simulations

Evaluating experimental, knowledge-based and computational cocrystal screening methods to advance drug-drug cocrystal fixed-dose combination development

Fixed-dose combinations (FDCs) offer significant advantages to patients and the pharmaceutical industry alike through improved dissolution profiles, synergistic effects and extended patent lifetimes. Identifying whether two active pharmaceutical ingredients have the potential to form a drug-drug cocrystal (DDC) or interact is an essential step in determining the most suitable type of FDC to formulate. The lack of coherent strategies to determine if two active pharmaceutical ingredients that can be co-administered can form a cocrystal, has significantly impacted DDC commercialisation. This review aims to accelerate the development of FDCs and DDCs by evaluating existing experimental, knowledge-based and computational cocrystal screening methods; the background of their development, their application in screening for cocrystals and DDCs, and their limitations are discussed. The evaluation provided in this review will act as a guide for selecting suitable screening methods to accelerate FDC development.

Interplay of halogen bonding and solvation in protein-ligand binding

Halogen bonding is increasingly utilized in efforts to achieve high affinity and selectivity of molecules designed to bind proteins, making it paramount to understand the relationship between structure, dynamics, and thermodynamic driving forces. We present a detailed analysis addressing this problem using a series of protein-ligand complexes involving single halogen substitutions - F, Cl, Br, and I - and nearly identical structures. Isothermal titration calorimetry reveals an increasingly favorable binding enthalpy from F to I that correlates with the halogen size and σ-hole electropositive character, but is partially counteracted by unfavorable entropy, which is constant from F to Cl and Br, but worse for I. Consequently, the binding free energy is roughly equal for Cl, Br, and I. QM and solvation-free-energy calculations reflect an intricate balance between halogen bonding, hydrogen bonds, and solvation. These advances have the potential to aid future drug design initiatives involving halogenated compounds.

Noncovalent interaction guided selectivity of haloaromatic isomers in a flexible porous coordination polymer

Porous, supramolecular structures exhibit preferential encapsulation of guest molecules, primarily by means of differences in the order of (noncovalent) interactions. The encapsulation preferences can be for geometry (dimension and shape) and the chemical nature of the guest. While geometry-based sorting is relatively straightforward using advanced porous materials, designing a "chemical nature" specific host is not. To introduce "chemical specificity", the host must retain an accessible and complementary recognition site. In the case of a supramolecular, porous coordination polymer (PCP) [Zn(o-phen)(ndc)] (o-phen: 1,10-phenanthroline, ndc: 2,6-naphthalenedicarboxylate) host, equipped with an adaptable recognition pocket, we have discovered that the preferential encapsulation of a haloaromatic isomer is not only for dimension and shape, but also for the "chemical nature" of the guest. This selectivity, i.e., preference for the dimension, shape and chemical nature, is not guided by any complementary recognition site, which is commonly required for "chemical specificity". Insights from crystal structures and computational studies unveil that the differences in the different types of noncovalent host-guest interaction strengths, acting in a concerted fashion, yield the unique selectivity.

Gate Opening without Volume Change Triggers Cooperative Gas Interactions, Underpins an Isotherm Step in Metal-Organic Frameworks

Three halogenated metal-organic frameworks (MOFs) reported recently exhibited a second step in their CO₂ gas adsorption isotherms. The emergence of halogen-bonding interactions beyond a threshold gas pressure between the framework halogen and the CO₂ guest was conjectured to be the underlying reason for the additional step in the isotherm. Our investigation employing periodic density functional theory calculations did not show significant interactions between the halogen and CO₂ molecules. Further, using a combination of DFT-based ab initio molecular dynamics and grand canonical Monte Carlo simulations, we find that the increased separation of framework nitrate pairs facing each other across the pore channel enables the accommodation of an additional CO₂ molecule which is further stabilized by cooperative interactions─an observation that facilely explains the second isotherm step. The increased separation between the nitrate groups can occur without any lattice expansion, consistent with experiments. The results point to a structural feature to achieve this isotherm step in MOFs that neither possess large pores nor exhibit large-scale structural changes such as breathing.

Molecular Simulation of Covalent Bond Dynamics in Liquid Silicon

Many atomic liquids can form transient covalent bonds reminiscent of those in the corresponding solid states. These directional interactions dictate many important properties of the liquid state, necessitating a quantitative, atomic-scale understanding of bonding in these complex systems. A prototypical example is liquid silicon, wherein transient covalent bonds give rise to local tetrahedral order and consequent nontrivial effects on liquid-state thermodynamics and dynamics. To further understand covalent bonding in liquid silicon, and similar liquids, we present an ab initio-simulation-based approach for quantifying the structure and dynamics of covalent bonds in condensed phases. Through the examination of structural correlations among silicon nuclei and maximally localized Wannier function centers, we develop a geometric criterion for covalent bonds in liquid Si. We use this to monitor the dynamics of transient covalent bonding in the liquid state and estimate a covalent bond lifetime. We compare covalent bond dynamics to other processes in liquid Si and similar liquids and suggest experiments to measure the covalent bond lifetime.

Lone Pair Rotational Dynamics in Solids

Traditional classifications of crystalline phases focus on nuclear degrees of freedom. Through the examination of both electronic and nuclear structure, we introduce the concept of an electronic plastic crystal. Such a material is classified by crystalline nuclear structure, while localized electronic degrees of freedom-here lone pairs-exhibit orientational motion at finite temperatures. This orientational motion is an emergent phenomenon arising from the coupling between electronic structure and polarization fluctuations generated by collective motions, such as phonons. Using ab initio molecular dynamics simulations, we predict the existence of electronic plastic crystal motion in halogen crystals and halide perovskites, and suggest that such motion may be found in a broad range of solids with lone pair electrons. Such fluctuations in the charge density should be observable, in principle, via synchrotron scattering.