Evidence for CCl/CBr⋅⋅⋅π Interactions as an Important Contribution to Protein-Ligand Binding Affinity

Photoregulated Supramolecular Polymerization through Halogen Bonding

Supramolecular polymers are able to change their structure, morphology and function in response to external stimuli. However, controlling the independence of stimuli-responses in these systems is challenging. Herein, we exploit halogen bonding (XB) as a reversible network element to regulate the photoresponsive and adaptive behavior of supramolecular polymers. To this end, we have designed a two-component system comprising an amphiphilic XB acceptor with the ability to self-assemble in aqueous media (OPE-Py) and a molecule with a dual photoresponsive and XB donor function [(E)-Azo-I]. OPE-Py self-assembles in aqueous media into supramolecular polymers, which transform into nanoparticle assemblies upon co-assembly with (E)-Azo-I. Interestingly, a third type of assembly (2D sheets) is obtained if OPE-Py is treated with (E)-Azo-I and exposed to photoirradiation. At ambient conditions, both nanoparticles and 2D sheets remain invariant over time. However, heating dissociates the XB interactions present in both assemblies, resulting in their transformation to the original fiber-like morphology of OPE-Py. Thus, breaking the communication between self-assembly and the stimuli-responses upon heating restores the original state of the system, drawing parallels to feedback loops in programming language. This work broadens the still limited scope of XB in solution assemblies and paves the way for multifunctional adaptive supramolecular systems.

Tetramerization of deoxyadenosine kinase meets the demands of a DNA replication substrate challenge in Giardia intestinalis

The protozoan parasite Giardia intestinalis is one of only a few organisms lacking de novo synthesis of DNA building blocks (deoxyribonucleotides). Instead, the parasite relies exclusively on salvaging deoxyadenosine and other deoxyribonucleosides from its host environment. Here, we report that G. intestinalis has a deoxyribonucleoside kinase with a 1000-fold higher catalytic efficiency (kcat/KM) for deoxyadenosine than the corresponding mammalian kinases and can thereby provide sufficient deoxyadenosine triphosphate levels for DNA synthesis despite the lack of de novo synthesis. Several deoxyadenosine analogs were also potent substrates and showed comparable EC50 values on cultured G. intestinalis cells as metronidazole, the current first-line treatment, with the additional advantage of being effective against metronidazole-resistant parasites. Structural analysis using cryo-EM and X-ray crystallography showed that the enzyme is unique within its family of deoxyribonucleoside kinases by forming a tetramer stabilized by extended N- and C-termini in a novel dimer-dimer interaction. Removal of the two termini resulted in lost ability to form tetramers and a markedly reduced affinity for the deoxyribonucleoside substrate. The development of highly efficient deoxyribonucleoside kinases via oligomerization may represent a critical evolutionary adaptation in organisms that rely solely on deoxyribonucleoside salvage.

Halogen Bonding on Water─A Drop in the Ocean?

Halogen bonding is a valuable interaction in drug design, offering an unconventional way to influence affinity and selectivity by leveraging the halogen atoms' ability to form directional bonds. The present study evaluates halogen-water interactions within protein binding sites, demonstrating that targeting a water molecule via halogen bonding can in specific cases contribute beneficially to ligand binding. In solving and examining the crystal structure of 2-cyclopentyl-7-iodo-1H-indole-3-carbonitrile bound to DYRK1a kinase, we identified a notable iodine-water interaction, where water accepts a halogen bond with good geometric and energetic features. This starting point triggered further investigations into the prevalence of such interactions across various halogen-bearing ligands (chlorine, bromine, iodine) in the PDB. Using QM calculations (MP2/TZVPP), we highlight the versatility and potential benefits of such halogen-water interactions, particularly when the water molecule is a stable part of the binding site's structured environment. While the interaction energies with water are lower compared to other typical halogen bond acceptors, we deem this different binding strength essential for reducing desolvation costs. We suggest that "interstitial" water molecules, as stable parts of the binding site engaging in multiple strong interactions, could be prime targets for halogen bonding. Further systematic studies, combining high-resolution crystal structures and quantum chemistry, are required to scrutinize whether halogen bonding on water is more than a "drop in the ocean".

Arylboronic Acid Pinacol Esters as Stable Boron Sources for Dihydrodibenzoborepin Derivatives and a Dibenzoborole

The general synthesis of boron-containing cyclic compounds (boracycles) necessitates toxic organotin precursors or highly reactive boron halides. Here, we report the synthesis of seven- and five-membered boracycles utilizing arylboronic acid pinacol esters (ArBpins) as stable boron sources. Grignard reagents generated from 2,2'-dibromodibenzyl or 2,2'-dibromobiphenyl reacted with ArBpins, where Ar = 9-anthryl (Anth), 2,4,6-trimethylphenyl (Mes), or 2,4,6-triisopropylphenyl (Tip), to give 10,11-dihydro-5H-dibenzo[b,f]borepins or dibenzoborole derivatives. This Bpin-based method was successfully applied to a one-shot double boracycle formation, providing a dihydrodibenzoborepin-anthracene-dihydrodibenzoborepin triad molecule in a good yield. The dihydrodibenzoborepin bearing the Anth group was directly converted to the unsaturated borepin by NBS/AIBN. All products were characterized by NMR, HRMS, and in some cases, single-crystal X-ray diffraction analysis. Additionally, the photophysical properties of the products are also reported.

Discovery of a Novel and Potent Cyclin-Dependent Kinase 8/19 (CDK8/19) Inhibitor for the Treatment of Cancer

The mediator kinases CDK8 and CDK19 control the dynamic transcription of selected genes in response to various signals and have been shown to be hijacked to sustain hyperproliferation by various solid and liquid tumors. CDK8/19 is emerging as a promising anticancer therapeutic target. Here, we report the discovery of compound 12, a novel small molecule CDK8/19 inhibitor. This molecule demonstrated not only decent enzymatic and cellular activities but also remarkable selectivity in CDK and kinome panels. Besides, compound 12 also displayed favorable ADME profiles including low CYP1A2 inhibition, acceptable clearance, and high oral bioavailability in multiple preclinical species. Robust in vivo PD and efficacy studies in mice models further demonstrated its potential use as mono- and combination therapy for the treatment of cancers.

Identification of an inhibitory pocket in falcilysin provides a new avenue for malaria drug development

Identification of new druggable protein targets remains the key challenge in the current antimalarial development efforts. Here we used mass-spectrometry-based cellular thermal shift assay (MS-CETSA) to identify potential targets of several antimalarials and drug candidates. We found that falcilysin (FLN) is a common binding partner for several drug candidates such as MK-4815, MMV000848, and MMV665806 but also interacts with quinoline drugs such as chloroquine and mefloquine. Enzymatic assays showed that these compounds can inhibit FLN proteolytic activity. Their interaction with FLN was explored systematically by isothermal titration calorimetry and X-ray crystallography, revealing a shared hydrophobic pocket in the catalytic chamber of the enzyme. Characterization of transgenic cell lines with lowered FLN expression demonstrated statistically significant increases in susceptibility toward MK-4815, MMV000848, and several quinolines. Importantly, the hydrophobic pocket of FLN appears amenable to inhibition and the structures reported here can guide the development of novel drugs against malaria.

Halogen Bonding in Brominated BODIPY Crystals: a Crystallographic and Computational Study

The study of halogen bonds (XBs) has been a subject of great interest in recent years due to its clear application in catalysis, liquid crystals, and crystal engineering. In this study, we analyzed the intermolecular interactions, in particular halogen bonds in BODIPYs with an increasing number of bromine atoms. The computational study included analyses through three different methods: the first approach of close contacts provided by mercury, then the expanded approach of the electron density partition of the molecules in the crystals provided by the analysis of Hirshfeld surfaces, and finally, the approach of the Quantum Theory of Atoms in Molecules (QT-AIM) to characterize the non-covalent interactions through finding electron density critical points between atoms and between neighboring molecules. The use of different computational methods allowed to gain insight into the interactions directing the crystal packing as the number of bromine atoms increased in the BODIPY moiety. Monocoordinated and bifurcated halogen bonds involving halide/halide were found. The penta-brominated BODIPY showed four-center cyclic nodes where each node is linked via XBs. This kind of motif can be useful in supramolecular chemistry and self-assembly.

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 ).

Abnormalities of the Halogen Bonds in the Complexes between Y2CTe (Y = H, F, CH3) and XF (X = F, Cl, Br, I)

In this work, the hydrogen bonds and halogen bonds in the complexes between Y2CTe (Y = H, F, CH3) and XF (X = F, Cl, Br, I) have been studied by quantum chemical calculations. We found three interesting abnormalities regarding the interactions. Firstly, the strength of halogen bonds increases in the order of IF < BrF < ClF < F2. Secondly, the halogen bonds formed by F2 are very strong, with an interaction energy in the range between −199.8 and −233.1 kJ/mol. Thirdly, all the halogen bonds are stronger than the hydrogen bonds in the systems we examined. All these results are against the general understanding of halogen bonds. These apparent abnormal properties are reconciled with the high polarizability of the Te atom and the strong inducing effect of F on the Te atom of Y2CTe. These findings provide a new perspective on halogen bonds. Additionally, we also proposed bonding distance-based methods to compare the strength of halogen/hydrogen bonds formed between different donor atoms and the same acceptor atom.

Hydration dynamics and IR spectroscopy of 4-fluorophenol

Halogenated groups are relevant in pharmaceutical applications and potentially useful spectroscopic probes for infrared spectroscopy. In this work, the structural dynamics and infrared spectroscopy of para-fluorophenol (F-PhOH) and phenol (PhOH) is investigated in the gas phase and in water using a combination of experiment and molecular dynamics (MD) simulations. The gas phase and solvent dynamics around F-PhOH and PhOH is characterized from atomistic simulations using empirical energy functions with point charges or multipoles for the electrostatics, Machine Learning (ML) based parametrizations and with full ab initio (QM) and mixed Quantum Mechanical/Molecular Mechanics (QM/MM) simulations with a particular focus on the CF- and OH-stretch region. The CF-stretch band is heavily mixed with other modes whereas the OH-stretch in solution displays a characteristic high-frequency peak around 3600 cm-1 most likely associated with the -OH group of PhOH and F-PhOH together with a characteristic progression below 3000 cm-1 due to coupling with water modes which is also reproduced by several of the simulations. Solvent and radial distribution functions indicate that the CF-site is largely hydrophobic except for simulations using point charges which renders them unsuited for correctly describing hydration and dynamics around fluorinated sites. The hydrophobic character of the CF-group is particularly relevant for applications in pharmaceutical chemistry with a focus on local hydration and interaction with the surrounding protein.

Probing Factor Xa Protein-Ligand Interactions: Accurate Free Energy Calculations and Experimental Validations of Two Series of High-Affinity Ligands

The accurate prediction of protein-ligand binding affinity belongs to one of the central goals in computer-based drug design. Molecular dynamics (MD)-based free energy calculations have become increasingly popular in this respect due to their accuracy and solid theoretical basis. Here, we present a combined study which encompasses experimental and computational studies on two series of factor Xa ligands, which enclose a broad chemical space including large modifications of the central scaffold. Using this integrated approach, we identified several new ligands with different heterocyclic scaffolds different from the previously identified indole-2-carboxamides that show superior or similar affinity. Furthermore, the so far underexplored terminal alkyne moiety proved to be a suitable non-classical bioisosteric replacement for the higher halogen-π aryl interactions. With this challenging example, we demonstrated the ability of the MD-based non-equilibrium free energy calculation approach for guiding crucial modifications in the lead optimization process, such as scaffold replacement and single-site modifications at molecular interaction hot spots.

A bibliographic survey of the supramolecular architectures sustained by delocalised C–I⋯π(arene) interactions in metal-organic crystals

A survey of the crystallographic literature of metal-organic crystal structures for the presence of C–I···π(arene) interactions where the iodide atom occupies a position close to plumb to the ring centroid, corresponding to a delocalised interaction, and is within the assumed sum of the van der Waals radii, i.e. 3.88 Å, has been undertaken. The majority of the 26 identified examples feature supramolecular chains of varying topology whereby C–I···π(arene) contacts are readily identified and apparently operating independently of other obvious supramolecular synthons. The next most prevalent supramolecular aggregate was zero-dimensional, containing up to a maximum of three molecules. While there were three examples of two-dimensional arrays among a series of isostructural crystal structures, no examples of three-dimensional structures largely sustained by C–I···π(arene) interactions were noted. This distribution of supramolecular aggregation patterns matched that noted for all-organic systems. In terms of the overall adoption rate, delocalised C–I···π(arene) interactions were found in 3% of crystals of metal-organic species where they could form, a percentage lower than 4% noted for all-organic crystals.

Multi-Substituted Quinolines as HIV-1 Integrase Allosteric Inhibitors

Allosteric HIV-1 integrase (IN) inhibitors, or ALLINIs, are a new class of antiviral agents that bind at the dimer interface of the IN, away from the enzymatic catalytic site and block viral replication by triggering an aberrant multimerization of the viral enzyme. To further our understanding of the important binding features of multi-substituted quinoline-based ALLINIs, we have examined the IN multimerization and antiviral properties of substitution patterns at the 6 or 8 position. We found that the binding properties of these ALLINIs are negatively impacted by the presence of bulky substitutions at these positions. In addition, we have observed that the addition of bromine at either the 6 (6-bromo) or 8 (8-bromo) position conferred better antiviral properties. Finally, we found a significant loss of potency with the 6-bromo when tested with the ALLINI-resistant IN A128T mutant virus, while the 8-bromo analog retained full effectiveness.

The involvement of CYP1A2 in biodegradation of dioxins in pigs

2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is one of the most harmful chemicals showing resistance to biodegradation. The majority of TCDD effects is mediated by the aryl hydrocarbon receptor (AhR) pathway. TCDD binding to AhR results in the activation of cytochrome P450 enzymes (CYP1A1, CYP1A2, CYP1B1) involved in dioxin biodegradation. The goal of the study was to explore the potentialrole of CYP1A2 in the metabolism of TCDD. We investigated a molecular structure of CYP1A2 and the binding selectivity and affinity between the pig CYP1A2 and: 1/ DiCDD or TCDD (dioxins differing in toxicity and biodegradability) or 2/ their selected metabolites. pCYP1A2 demonstrated higher affinity towards DiCDD and TCDD than other pCYP1 enzymes. All dioxin-pCYP1A2 complexes were found to be stabilized by hydrophobic interactions. The calculated distances between the heme oxygen and the dioxin carbon nearest to the oxygen, reflecting the hydroxylating potential of CYP1A2, were higher than in other pCYP1 enzymes. The distances between the heme iron and the nearest dioxin carbon exceeded 5 Å, a distance sufficient to allow the metabolites to leave the active site. However, the molecular dynamics simulations revealed that two access channels of CYP1A2 were closed upon binding the majority of the examined dioxins. Moreover, the binding of dioxin metabolites did not promote opening of channel S–an exit for hydroxylated products. It appears that the undesired changes in the behavior of access channels prevail over the hydroxylating potential of CYP1A2 towards TCDD and the favorable distances, ultimately trapping the metabolites at the enzyme’s active site.

Interpretation of Structure-Activity Relationships in Real-World Drug Design Data Sets Using Explainable Artificial Intelligence

In silico models based on Deep Neural Networks (DNNs) are promising for predicting activities and properties of new molecules. Unfortunately, their inherent black-box character hinders our understanding, as to which structural features are important for activity. However, this information is crucial for capturing the underlying structure-activity relationships (SARs) to guide further optimization. To address this interpretation gap, "Explainable Artificial Intelligence" (XAI) methods recently became popular. Herein, we apply and compare multiple XAI methods to projects of lead optimization data sets with well-established SARs and available X-ray crystal structures. As we can show, easily understandable and comprehensive interpretations are obtained by combining DNN models with some powerful interpretation methods. In particular, SHAP-based methods are promising for this task. A novel visualization scheme using atom-based heatmaps provides useful insights into the underlying SAR. It is important to note that all interpretations are only meaningful in the context of the underlying models and associated data.

Structural systematics in the isomorphous binary co-crystal solvates comprising 2,2'-dithiodibenzoic acid, 4-halobenzoic acid and dimethylformamide (1:1:1), for halide = chloride, bromide and iodide

The 1:1:1 binary co-crystal solvates formulated as 2,2'-dithiodibenzoic acid (DTBA), 4-halobenzoic acid (4-XBA) and dimethylformamide (DMF) for X = Cl (1), Br (2) and I (3) are isomorphous and the...

N-chlorobenzimidazoles as efficient and structurally diverse amphoteric halogen bond donors in crystal engineering

The ability of amphoteric N-chlorobenzimidazoles to self-associate into 1D chains through strong and linear N−Cl•••N halogen bond interactions is demonstrated. Less polarisable Cl atom is strongly activated thanks to the...

Green by Design: Convergent Synthesis, Computational Analyses, and Activity Evaluation of New FXa Inhibitors Bearing Peptide Triazole Linking Units

Green chemistry implementation has led to promising results in waste reduction in the pharmaceutical industry. However, the early sustainable development of pharmaceutically active compounds and ingredients remains a considerable challenge. Herein, we wish to report a green synthesis of new pharmaceutically active peptide triazoles as potent factor Xa inhibitors, an important drug target associated with the treatment of diverse cardiovascular diseases. The new inhibitors were synthesized in three steps, featuring cycloaddition reactions (high atom economy), microwave-assisted organic synthesis (energy efficiency), and copper nanoparticle catalysis, thus featuring Earth-abundant metals. The molecules obtained showed FXa inhibition, with IC50-values as low as 17.2 μM and no associated cytotoxicity in HEK293 and HeLa cells. These results showcase the environmental potential and chemical implications of the applied methodologies for the development of new molecules with pharmacological potential.

Revisiting the covalent nature of halogen bonding: a polarized three-center four-electron bond

As an important intermolecular interaction, halogen bonding has been studied extensively, but its nature still suffers from controversy without one uniform essence. Electrostatics, charge transfer, polarization and dispersion are emphasized, but the covalent nature is usually overlooked except for the strong halogen bonding species I₃⁻, which is widely accepted as a result of a three-center four-electron (3c-4e) interaction. In our study, the potential energy surface of I₃⁻ has been evaluated to explore the dissociation from I₃⁻ to I₂⋯I⁻. We found that different from an equivalent 3c-4e bond in I₃⁻, I₂⋯I⁻ can be rationalized by a polarized one. In addition, when the orbitals are polarized, it is exactly what traditional charge transfer or the popular σ-hole picture describes. I₃⁻ can be described by the Lewis theory model with the middle I⁺ cation serving as the Lewis acid and two terminal I⁻ anions acting as Lewis base. Therefore, we further extended this model to a series of I-containing species with chemical composition of L–I⁺–L, F⁻–I⁺–L and H₃P–I⁺–L (L = OH⁻, F⁻, Cl⁻, Br⁻, I⁻, PH₃, NH₃, H₂S, HI, H₂O, HBr and HCl) to explore the nature of halogen bonding. When the forces of two bases around I⁺ are the same, it corresponds to an equivalent 3c-4e bond, such as I₃⁻. Otherwise, it is a polarized multicenter bond, such as I₂⋯I⁻. This work gives a new insight into the nature of halogen bonding compounds: besides the well-known I₃⁻, the nature of the other species is also a multicenter bond, existing as equivalent and polarized 3c-4e bonds, respectively.

The halogen bond could be described with a polarized 3c-4e bond.

Controlling the molecular arrangement of racemates through weak interactions: the synergy between π-interactions and halogen bonds

POX and NX halogen bonds in combination with π-stacking interactions lead to the sorting of π-extended R- and S-isomers. Theoretical calculations point to a positive synergistic effect between the π-interactions and the halogen bonds to be the origin of such phenomena. As a result, enantiomeric building blocks form homoleptically connected quadrangular structures.

Unraveling Halogen Effects in Supramolecular Polymerization

Halogens play a crucial role in numerous natural processes and synthetic materials due to their unique physicochemical properties and the diverse interactions they can engage in. In the field of supramolecular polymerization, however, halogen effects remain poorly understood, and investigations have been restricted to halogen bonding or the inclusion of polyfluorinated side groups. Recent contributions from our group have revealed that chlorine ligands greatly influence molecular packing and pathway complexity phenomena of various metal complexes. These results prompted us to explore the role of the halogen nature on supramolecular polymerization, a phenomenon that has remained unexplored to date. To address this issue, we have designed a series of archetypal bispyridyldihalogen Pt^II complexes bearing chlorine (1), bromine (2), or iodine (3) and systematically compared their supramolecular polymerization in nonpolar media using various experimental methods and theory. Our studies reveal a remarkably different supramolecular polymerization for the three compounds, which can undergo two competing pathways with either slipped (kinetic) or parallel (thermodynamic) molecular packing. The halogen exerts an inverse effect on the energetic levels of the two self-assembled states, resulting in a single thermodynamic pathway for 3, a transient kinetic species for 2, and a hidden thermodynamic state for 1. This seesaw-like bias of the energy landscape can be traced back to the involvement of the halogens in weak N-H···X hydrogen-bonding interactions in the kinetic pathway, whereas in the thermodynamic pathway the halogens are not engaged in the stabilizing interaction motif but rather amplify solvophobic effects.

Quantum-mechanical property prediction of solvated drug molecules: what have we learned from a decade of SAMPL blind prediction challenges?

Joint academic-industrial projects supporting drug discovery are frequently pursued to deploy and benchmark cutting-edge methodical developments from academia in a real-world industrial environment at different scales. The dimensionality of tasks ranges from small molecule physicochemical property assessment over protein-ligand interaction up to statistical analyses of biological data. This way, method development and usability both benefit from insights gained at both ends, when predictiveness and readiness of novel approaches are confirmed, but the pharmaceutical drug makers get early access to novel tools for the quality of drug products and benefit of patients. Quantum-mechanical and simulation methods particularly fall into this group of methods, as they require skills and expense in their development but also significant resources in their application, thus are comparatively slowly dripping into the realm of industrial use. Nevertheless, these physics-based methods are becoming more and more useful. Starting with a general overview of these and in particular quantum-mechanical methods for drug discovery we review a decade-long and ongoing collaboration between Sanofi and the Kast group focused on the application of the embedded cluster reference interaction site model (EC-RISM), a solvation model for quantum chemistry, to study small molecule chemistry in the context of joint participation in several SAMPL (Statistical Assessment of Modeling of Proteins and Ligands) blind prediction challenges. Starting with early application to tautomer equilibria in water (SAMPL2) the methodology was further developed to allow for challenge contributions related to predictions of distribution coefficients (SAMPL5) and acidity constants (SAMPL6) over the years. Particular emphasis is put on a frequently overlooked aspect of measuring the quality of models, namely the retrospective analysis of earlier datasets and predictions in light of more recent and advanced developments. We therefore demonstrate the performance of the current methodical state of the art as developed and optimized for the SAMPL6 pKa and octanol-water log P challenges when re-applied to the earlier SAMPL5 cyclohexane-water log D and SAMPL2 tautomer equilibria datasets. Systematic improvement is not consistently found throughout despite the similarity of the problem class, i.e. protonation reactions and phase distribution. Hence, it is possible to learn about hidden bias in model assessment, as results derived from more elaborate methods do not necessarily improve quantitative agreement. This indicates the role of chance or coincidence for model development on the one hand which allows for the identification of systematic error and opportunities toward improvement and reveals possible sources of experimental uncertainty on the other. These insights are particularly useful for further academia-industry collaborations, as both partners are then enabled to optimize both the computational and experimental settings for data generation.

Theory of Long-Lived Room-Temperature Phosphorescence in Organic Aggregates

ConspectusRoom-temperature phosphorescence (RTP) with a long afterglow from purely organic molecular aggregates has recently attracted many investigations because traditionally only inorganic and transition-metal complexes can emit phosphorescence at room temperature. Purely organic molecules can exhibit phosphorescence only at cryogenic temperatures and under inert conditions in solution. However, recently, a number of organic compounds have been found to demonstrate bright RTP upon aggregation, sometimes with a remarkable morphology dependence. We intended to rationalize such aggregation-induced organic RTP through theoretical investigation and quantum chemistry calculations by invoking intermolecular interaction effects. And we have identified the molecular descriptors for the molecular design of RTP materials.In this Account, we started with the proposition of the mechanism of intermolecular electrostatic-interaction-induced RTP at the molecular level by using molecular dynamics simulations, hybrid quantum mechanics, and molecular mechanics (QM/MM) coupled with the thermal vibration correlation function (TVCF) formalism we developed earlier. The effective intermolecular electrostatic interactions could stem from a variety of interactions in different organic RTP crystals, such as hydrogen bonding, π-halogen bonding, anion-π⁺ interaction, and d-pπ bonds and so forth. We find that these interactions can change the molecular orbital compositions involved in the lowest-lying singlet and triplet excited states that are responsible for phosphorescence, either through facilitating intersystem crossing from the excited-state singlet to the triplet and/or suppressing the nonradiative decay process from the lowest triplet to the ground state. This underlying RTP mechanism is believed to be very helpful in systematically and comprehensively understanding the aggregation/crystal-induced persistent organic RTP, which has been applied to explain a number of experiments.We then propose the molecular descriptors to characterize the phosphorescence efficiency and lifetime, respectively, derived from fundamental photophysical processes and requirements to obey the El-Sayed rule and generate phosphorescence. For a prototypical RTP system consisting of a carbonyl group and π-conjugated segments, the excited states can be regarded as an admixture of n → π* (with portion α) and π → π* (with portion β). The intersystem crossing (ISC) rate of S₁ → T_n is mostly governed by the modification of the product of α and β, and the nonradiative rate of T₁ → S₀ is determined by the β value of T₁. Thus, we employ γ = α × β and β to describe the phosphorescence efficiency and lifetime, respectively, which have been successfully applied in the molecular design of efficient and long-lived RTP systems in experiments. The molecular descriptors outlined in this Account, which are easily obtained from simple quantum chemistry calculations, are expected to play important roles in the machine-learning-based molecular screening in the future.

Supramolecular architectures sustained by delocalised C–I⋯π(arene) interactions in molecular crystals and the propensity of their formation

A survey of delocalised C–I⋯π(chelate ring) interactions is presented.

Experimental and computational evidence for a stabilising C–Cl(lone-pair)⋯π(chelate-ring) interaction

Stabilising C–Cl(lone-pair)⋯π(chelate ring) interactions are described.

Probing Halogen-π versus CH-π Interactions in Molecular Balance

Molecular balances based on the dibenzobicyclo[3.2.2]nonane template enable probing of the competition between halogen-π and CH-π interactions. Structural, NMR spectroscopic, and computational analyses revealed that the π system can favorably interact both with C-X or C-H functionalities, depending on the size of the functional group.

Rising trend on the halogen and non-halogen derivatives (Br, Cl, CF3, F, CH3 and NH2) in ruminal β-d-Xylopyranose - a quantum chemical perspective

The utmost aim of the current study is to find significance of the binding affinity in the halogen and non-halogen derivatives: Br, Cl, CF₃, F, CH₃ and NH₂ of β-d-Xylopyranose with the hinge region amino acids of ruminant-β-glycosidase. The interaction energy analysis was carried out in detail through various density functional studies as M062X/def2-QZVP, M062X/LANL2DZ, B3LYP/LANL2DZ and M06HF/LANL2DZ level of theories. The total interaction energy of halogen derivatives: Br, Cl, F and CF₃ are -618.21, -599.00, -720.45 and -553.08 kcal/mol respectively, and non-halogen derivative: amine group (NH₂) is -87.96 kcal/mol at M062X/def2-QZVP level of theory, which exist with strong binding affinity. Ligand properties: dipole moment, polarizability, volume, molecular mass, electrostatic potential map was evaluated to understand its electrostatic and structural behavior. The nature of the bonds was inferred from the electrostatic potential map for all the halogen and non-halogen derivatives ligand. The stabilization energy from NBO analysis reveals the stability of single hydrogen and halogen bonds (N-H…Br, C-Br…O, N-H…Cl, C-Cl…O, O-H…F, C-H…F, N-H…F, C-F…O, N-H…O, O-H…O, N-H…N, O-H…N) in β-d-Xylopyranose and its derivatives. Overall, this study paves way for scientist and medicinal chemist in modelling new drugs. Further, it suggests mutations that increase the binding and may enhance the catalytic action and strengthen the complex diet in animals and hence recommended for experimental synthesis.Communicated by Ramaswamy H. Sarma.

The Intriguing Effects of Substituents in the N-Phenethyl Moiety of Norhydromorphone: A Bifunctional Opioid from a Set of "Tail Wags Dog" Experiments

(−)-N-Phenethyl analogs of optically pure N-norhydromorphone were synthesized and pharmacologically evaluated in several in vitro assays (opioid receptor binding, stimulation of [³⁵S]GTPγS binding, forskolin-induced cAMP accumulation assay, and MOR-mediated β-arrestin recruitment assays). “Body” and “tail” interactions with opioid receptors (a subset of Portoghese’s message-address theory) were used for molecular modeling and simulations, where the “address” can be considered the “body” of the hydromorphone molecule and the “message” delivered by the substituent (tail) on the aromatic ring of the N-phenethyl moiety. One compound, N-p-chloro-phenethynorhydromorphone ((7aR,12bS)-3-(4-chlorophenethyl)-9-hydroxy-2,3,4,4a,5,6-hexahydro-1H-4,12-methanobenzofuro[3,2-e]isoquinolin-7(7aH)-one, 2i), was found to have nanomolar binding affinity at MOR and DOR. It was a potent partial agonist at MOR and a full potent agonist at DOR with a δ/μ potency ratio of 1.2 in the ([³⁵S]GTPγS) assay. Bifunctional opioids that interact with MOR and DOR, the latter as agonists or antagonists, have been reported to have fewer side-effects than MOR agonists. The p-chlorophenethyl compound 2i was evaluated for its effect on respiration in both mice and squirrel monkeys. Compound 2i did not depress respiration (using normal air) in mice or squirrel monkeys. However, under conditions of hypercapnia (using air mixed with 5% CO₂), respiration was depressed in squirrel monkeys.

Synthesis of furocoumarin-stilbene hybrids as potential multifunctional drugs against multiple biochemical targets associated with Alzheimer's disease

A series of furocoumarin-stilbene hybrids has been synthesized and evaluated in vitro for inhibitory effect against acetylcholinesterase (AChE), butyrylcholinestarase (BChE), β-secretase, cyclooxygenase-2 (COX-2), and lipoxygenase-5 (LOX-5) activities including free radical-scavenging properties. Among these hybrids, 8-(3,5-dimethoxyphenyl)-4-(3,5-dimethoxystyryl)furochromen-2-one 4h exhibited significant anticholinesterase activity and inhibitory effect against β-secretase, COX-2 and LOX-5 activities. 2,2-Diphenyl-1-picrylhydrazyl (DPPH) radical scavenging activity and an in vitro cell-based antioxidant activity assay involving lipopolysaccharide induced reactive oxygen species production revealed that 4h has capability of scavenging free radicals. Molecular docking into AChE, BChE, β-secretase, COX-2 and LOX-5 active sites has also been performed.

Chlorpromazine and Amitriptyline Are Substrates and Inhibitors of the AcrB Multidrug Efflux Pump

Efflux is an important mechanism in Gram-negative bacteria conferring multidrug resistance. Inhibition of efflux is an encouraging strategy to restore the antibacterial activity of antibiotics. Chlorpromazine and amitriptyline have been shown to behave as efflux inhibitors. However, their mode of action is poorly understood. Exposure of Salmonella enterica serovar Typhimurium and Escherichia coli to chlorpromazine selected for mutations within genes encoding RamR and MarR, regulators of the multidrug tripartite efflux pump AcrAB-TolC. Further experiments with S. Typhimurium containing AcrB D408A (a nonfunctional efflux pump) and chlorpromazine or amitriptyline resulted in the reversion of the mutant acrB allele to the wild type. Together, this suggests these drugs are AcrB efflux substrates. Subsequent docking studies with AcrB from S. Typhimurium and E. coli, followed by molecular dynamics simulations and free energy calculations showed that chlorpromazine and amitriptyline bind at the hydrophobic trap, a preferred binding site for substrates and inhibitors within the distal binding pocket of AcrB. Based on these simulations, we suggest that chlorpromazine and amitriptyline inhibit AcrB-mediated efflux by interfering with substrate binding. Our findings provide evidence that these drugs are substrates and inhibitors of AcrB, yielding molecular details of their mechanism of action and informing drug discovery of new efflux inhibitors.IMPORTANCE Efflux pumps of the resistance nodulation-cell division (RND) superfamily are major contributors to multidrug resistance for most of the Gram-negative ESKAPE (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter species) pathogens. The development of inhibitors of these pumps would be highly desirable; however, several issues have thus far hindered all efforts at designing new efflux inhibitory compounds devoid of adverse effects. An alternative route to de novo design relies on the use of marketed drugs, for which side effects on human health have been already assessed. In this work, we provide experimental evidence that the antipsychotic drugs chlorpromazine and amitriptyline are inhibitors of the AcrB transporter, the engine of the major RND efflux pumps in Escherichia coli and Salmonella enterica serovar Typhimurium. Furthermore, in silico calculations have provided a molecular-level picture of the inhibition mechanism, allowing rationalization of experimental data and paving the way for similar studies with other classes of marketed compounds.

A theoretical insight into non-covalent supramolecular interactions in the solid state structures of two octahedral iron(iii) complexes

The nature and characteristics of the C–H⋯π interactions that play an important role in crystal packing of two iron(iii) complexes have been discussed. The DFT calculations have been conducted to determine the interaction energies in these complexes.

Halogen-sensitive solvatochromism based on a phenolic polymer of tetraphenylethene

Herein, we describe the successful preparation of a methylene-bonded tetraphenylethene polymer using a phenolic-resin synthesis protocol. Our novel phenolic polymer showed solvatochromism in response to halogenated organic solvents. Solvatochromism is induced by halogen/π interactions between the polymer and the organic halide.

Herein, we describe novel phenolic polymer showed solvatochromism in response to halogenated organic solvents through halogen/π interactions.

Separation of halogenated benzenes enabled by investigation of halogen-π interactions with carbon materials

We report the existence of bimodal interactions, the π–π and halogen–π interactions, between the halogenated benzenes and aromatic materials.

The halogen–π (X–π) interaction is an intermolecular interaction between the electron-poor region of bonded halogen atoms and aromatic rings. We report an experimental evaluation of the halogen–π (X–π) interaction using liquid chromatography with carbon-material coated columns providing strong π interactions in the normal phase mode. A C₇₀-fullerene (C70)-coated column showed higher retentions for halogenated benzenes as the number of halogen substitutions increased as a result of X–π interactions. In addition, the strength of the X–π interaction increased in the order of F < Cl < Br < I. Changes to the UV absorption of C70 and the brominated benzenes suggested that the intermolecular interaction changed from the π–π interaction to X–π interaction as the number of bromo substitutions increased. Computer simulations also showed that the difference in dipole moments among structural isomers affected the strength of the π–π interaction. Furthermore, we concluded from small peak shifts in ¹H NMR and from computer simulations that the orbital interaction contributes to the X–π interactions. Finally, we succeeded in the one-pot separation of all isomers of brominated benzenes using the C70-coated column by optimizing the mobile phase conditions.

How Significant Are Unusual Protein-Ligand Interactions? Insights from Database Mining

We present a new approach to derive interaction propensities of protein-ligand atom pairs from mining of the Protein Data Bank. To ensure solid statistics, we use a line-of-sight contact filter and normalize the observed frequency of hits by a statistical null model based on exposed surface areas of atom types in the protein-ligand binding site. This allows us to investigate which intermolecular interactions and geometries are found more often than expected by chance in protein-ligand complexes. We focus our study on some of the unusual interactions that were postulated to be favorable, including σ-hole bonding of halogen and sulfur atoms, weak hydrogen bonding with fluorine as acceptor, and different types of dipolar interactions. Our results confirm some and challenge other common assumptions on these interactions and highlight other contact types that are yet underexplored in structure-based drug design.

Strategies for Late-Stage Optimization: Profiling Thermodynamics by Preorganization and Salt Bridge Shielding

Structural fixation of a ligand in its bioactive conformation may, due to entropic reasons, improve affinity. We present a congeneric series of thrombin ligands with a variety of functional groups triggering preorganization prior to binding. Fixation in solution and complex formation have been characterized by crystallography, isothermal titration calorimetry (ITC), and molecular dynamics (MD) simulations. First, we show why these preorganizing modifications do not affect the overall binding mode and how key interactions are preserved. Next, we demonstrate how preorganization thermodynamics can be largely dominated by enthalpy rather than entropy because of the significant population of low-energy conformations. Furthermore, a salt bridge is shielded by actively reducing its surface exposure, thus leading to an enhanced enthalpic binding profile. Our results suggest that the consideration of the ligand solution ensemble by MD simulation is necessary to predict preorganizing modifications that enhance the binding behavior of already promising binders.

Quantum Monte Carlo Study of the Water Dimer Binding Energy and Halogen-π Interactions

Halogen-π systems are involved with competition between halogen bonding and π-interaction. Using the diffusion quantum Monte Carlo (DMC) method, we compare the interaction of benzene with diatomic halogens (X₂: Cl₂/Br₂) with the typical hydrogen bonding in the water dimer, taking into account explicit correlations of up to three bodies. The benzene-Cl₂/Br₂ binding energies (13.07 ± 0.42/16.62 ± 0.02 kJ/mol) attributed to both halogen bonding and dispersion are smaller than but comparable to the typical hydrogen bonding in the water dimer binding energy (20.88 ± 0.27 kJ/mol). All of the above values are in good agreement with those from the coupled-cluster with single, double, and noniterative triple excitations (CCSD(T)) results at the complete basis set limit (benzene-Cl₂/Br₂: 12.78/16.17 kJ/mol; water dimer: 21.0 kJ/mol).

Relationships between Interaction Energy and Electron Density Properties for Homo Halogen Bonds of the [(A)nY-X···X-Z(B)m] Type (X = Cl, Br, I)

Relationships between interaction energy (Eint) and electron density properties at the X···X bond critical point or the d(X···X) distance were established for the large set of structures [(A)nY-X···X-Z(B)m] bearing the halogen bonds Cl···Cl, Br···Br, and I···I (640 structures in total). The best estimator of Eint is the kinetic energy density (Gb), which reasonably approximates the whole set of the structures as -Eint = 0.128Gb2 - 0.82Gb + 1.66 (R2 = 0.91, mean absolute deviation 0.39 kcal/mol) and demonstrates low dispersion. The potential and kinetic energy densities, electron density, and the d(X···X) distance behave similarly as estimators of Eint for the individual series Cl···Cl, Br···Br, and I···I. A number of the Eint(property) correlations are recommended for the practical application in the express estimates of the strength of the homo-halogen bonds.