Stacking Interactions between 9-Methyladenine and Heterocycles Commonly Found in Pharmaceuticals

The role of the methoxy group in approved drugs

The methoxy substituent is prevalent in natural products and, consequently, is present in many natural product-derived drugs. It has also been installed in modern drug molecules with no remnant of natural product features because medicinal chemists have been taking advantage of the benefits that this small functional group can bestow on ligand-target binding, physicochemical properties, and ADME parameters. Herein, over 230 methoxy-containing small-molecule drugs, as well as several fluoromethoxy-containing drugs, are presented from the vantage point of the methoxy group. Biochemical mechanisms of action, medicinal chemistry SAR studies, and numerous X-ray cocrystal structures are analyzed to identify the precise role of the methoxy group for many of the drugs and drug classes. Although the methoxy substituent can be considered as the hybridization of a hydroxy and a methyl group, the combination of these functionalities often results in unique effects that can amount to more than the sum of the individual parts.

Triazoles in Medicinal Chemistry: Physicochemical Properties, Bioisosterism, and Application

Triazole demonstrates distinctive physicochemical properties, characterized by weak basicity, various dipole moments, and significant dual hydrogen bond acceptor and donor capabilities. These features are poised to play a pivotal role in drug-target interactions. The inherent polarity of triazole contributes to its lower logP, suggesting the potential improvement in water solubility. The metabolic stability of triazole adds additional value to drug discovery. Moreover, the metal-binding capacity of the nitrogen atom lone pair electrons of triazole has broad applications in the development of metal chelators and antifungal agents. This Perspective aims to underscore the unique physicochemical attributes of triazole and its application. A comparative analysis involving triazole isomers and other heterocycles provides guiding insights for the subsequent design of triazoles, with the hope of offering valuable considerations for designing other heterocycles in medicinal chemistry.

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.

"Magic Chloro": Profound Effects of the Chlorine Atom in Drug Discovery

Chlorine is one of the most common atoms present in small-molecule drugs beyond carbon, hydrogen, nitrogen, and oxygen. There are currently more than 250 FDA-approved chlorine-containing drugs, yet the beneficial effect of the chloro substituent has not yet been reviewed. The seemingly simple substitution of a hydrogen atom (R = H) with a chlorine atom (R = Cl) can result in remarkable improvements in potency of up to 100,000-fold and can lead to profound effects on pharmacokinetic parameters including clearance, half-life, and drug exposure in vivo. Following the literature terminology of the "magic methyl effect" in drugs, the term "magic chloro effect" has been coined herein. Although reports of 500-fold or 1000-fold potency improvements are often serendipitous discoveries that can be considered "magical" rather than planned, hypotheses made to explain the magic chloro effect can lead to lessons that accelerate the cycle of drug discovery.

The pyridazine heterocycle in molecular recognition and drug discovery

The pyridazine ring is endowed with unique physicochemical properties, characterized by weak basicity, a high dipole moment that subtends π-π stacking interactions and robust, dual hydrogen-bonding capacity that can be of importance in drug-target interactions. These properties contribute to unique applications in molecular recognition while the inherent polarity, low cytochrome P450 inhibitory effects and potential to reduce interaction of a molecule with the cardiac hERG potassium channel add additional value in drug discovery and development. The recent approvals of the gonadotropin-releasing hormone receptor antagonist relugolix (24) and the allosteric tyrosine kinase 2 inhibitor deucravacitinib (25) represent the first examples of FDA-approved drugs that incorporate a pyridazine ring. In this review, the properties of the pyridazine ring are summarized in comparison to the other azines and its potential in drug discovery is illustrated through vignettes that explore applications that take advantage of the inherent physicochemical properties as an approach to solving challenges associated with candidate optimization.

Graphical Abstract

Synthesis, Biological Evaluation, and Computational Analysis of Biaryl Side-Chain Analogs of Solithromycin

There is an urgent need for new antibiotics to mitigate the existential threat posed by antibiotic resistance. Within the ketolide class, solithromycin has emerged as one of the most promising candidates for further development. Crystallographic studies of bacterial ribosomes and ribosomal subunits complexed with solithromycin have shed light on the nature of molecular interactions (π-stacking and H-bonding) between from the biaryl side-chain of the drug and key residues in the 50S ribosomal subunit. We have designed and synthesized a library of solithromycin analogs to study their structure-activity relationships (SAR) in tandem with new computational studies. The biological activity of each analog was evaluated in terms of ribosomal affinity (K_d determined by fluorescence polarization), as well as minimum inhibitory concentration assays (MICs). Density functional theory (DFT) studies of a simple binding site model identify key H-bonding interactions that modulate the potency of solithromycin analogs.

Energy, orbital and structural stacking landscape of a purine homodimer system

The multidimensional study, combining the extensive calculations of potential energy surfaces for the parallel-displaced configurations and methods such as energy decomposition and natural bond orbital analysis, has been carried out. The resulted data give an energy, orbital and structural landscapes of this biologically essential system. The balance of the two energy sources, electrostatic and dispersion, is clearly visible. The obtained results, taken as a whole, provide an insight into the hierarchy of intermolecular interactions in the purine system, together with their sources.

Converting SMILES to Stacking Interaction Energies

Predicting the strength of stacking interactions involving heterocycles is vital for several fields, including structure-based drug design. While quantum chemical computations can provide accurate stacking interaction energies, these come at a steep computational cost. To address this challenge, we recently developed quantitative predictive models of stacking interactions between druglike heterocycles and the aromatic amino acids Phe, Tyr, and Trp (DOI: 10.1021/jacs.9b00936 ). These models depend on heterocycle descriptors derived from electrostatic potentials (ESPs) computed using density functional theory and provide accurate stacking interactions without the need for expensive computations on stacked dimers. Herein, we show that these ESP-based descriptors can be reliably evaluated directly from the atom connectivity of the heterocycle, providing a means of predicting both the descriptors and the potential for a given heterocycle to engage in stacking interactions without resorting to any quantum chemical computations. This enables the rapid conversion of simple molecular representations (e.g., SMILES) directly into accurate stacking interaction energies using a freely available online tool, thereby providing a way to rank the stacking abilities of large sets of heterocycles.

Tuning Stacking Interactions between Asp-Arg Salt Bridges and Heterocyclic Drug Fragments

Stacking interactions can play an integral role in the strength and selectivity of protein-drug binding and are of particular interest given the ubiquity and variety of heterocyclic fragments in drugs. In addition to traditional stacking interactions between aromatic rings, stacking interactions involving heterocyclic drug fragments and protein salt bridges have also been observed. These "salt-bridge stacking interactions" can be quite strong but are not well understood. We studied stacked dimers of the acetate···guanidinium ion pair with a diverse set of 63 heterocycles using robust ab initio methods. The computed interaction energies span more than 10 kcal mol^-1, indicating the sensitivity of these salt-bridge stacking interactions to heterocycle features. Trends in both the strength and preferred geometry of these interactions can be understood through analyses of the electrostatic potentials and electric fields above the heterocycles. We have developed new heterocycle descriptors that quantify these effects and used them to create robust predictors of the strength of salt-bridge stacking interactions both in the gas phase and a protein-like dielectric environment. These predictive tools, combined with a set of qualitative guidelines, should facilitate the choice of heterocycles that maximize salt-bridge stacking interactions in drug binding sites.

Development and Validation of Quantum Mechanically Derived Force-Fields: Thermodynamic, Structural, and Vibrational Properties of Aromatic Heterocycles

A selection of several aromatic molecules, representative of the important class of heterocyclic compounds, has been considered for testing and validating an automated Force Field (FF) parametrization protocol, based only on Quantum Mechanical data. The parametrization is carried out separately for the intra- and intermolecular contributions, employing respectively the Joyce and Picky software packages, previously implemented and refined in our research group. The whole approach is here automated and integrated with a computationally effective yet accurate method, devised very recently ( J. Chem.

THEORY

Comput., 2018, 14, 543-556) to evaluate a large number of dimer interaction energies. The resulting quantum mechanically derived FFs are then used in extensive molecular dynamics simulations, in order to evaluate a number of thermodynamic, structural, and dynamic properties of the heterocycle's gas and liquid phases. The comparison with the available experimental data is good and furnishes a validation of the presented approach, which can be confidently exploited for the design of novel and more complex materials.

Stacking Interactions of Heterocyclic Drug Fragments with Protein Amide Backbones

Stacking interactions can be important enthalpic contributors to drug binding. Among the less well-studied stacking interactions are those occurring between an arene and the π-face of an amide group. Given the ubiquity of heterocycles in drugs, combined with the abundance of amides in the protein backbone, optimizing these noncovalent interactions can provide a potential route to enhanced drug binding. Previously, Diederich et al. (ChemMedChem 2013, 8, 397-404) studied stacked dimers of a model amide with a set of 18 heterocycles, showing that computed interaction energies correlate with the dipole moments of the heterocycles and providing guidelines for the optimization of these interactions. We considered stacked dimers of the same model amide with a larger set of 28 heterocycles common in pharmaceuticals, by using more robust ab initio methods. While the overall trends in these new data corroborate many of the results of Diederich et al., these data provide a more refined view of the nature of amide stacking interactions. We present a robust scoring function for amide stacking interaction energies based on the molecular dipole moment and strength of the electric field above the arene.

The π–π stacking of tanshinone I and isotanshinone I with phenylalanine: The effects of isomerization, complexation and environment

The misfolding and aggregation of amyloid-[Formula: see text] (A[Formula: see text] peptides into amyloid fibrils is regarded as one of the possible causes of Alzheimer’s disease (AD). Aromatic interactions (i.e. [Formula: see text]–[Formula: see text] interactions) between tanshinone drugs extracted from Chinese herb Salvia Miltiorrhiza (SM) and aromatic residues of A[Formula: see text] peptides have been shown to prevent further growth of amyloid aggregates. In this work, the effects of isomerization, complexation and polarity of environment on the strength of [Formula: see text]–[Formula: see text] stacking interactions of tanshinone I (TS1) and isotanshinone I (IS1), as the two diterpenoid quinones in the SM herb, and their complexes with Mg[Formula: see text] cation, IS1-Mg[Formula: see text] and TS1-Mg[Formula: see text], with phenylalanine (PHE), as an aromatic amino acid in A[Formula: see text] structure have been investigated using the quantum mechanical calculations in the gas-phase, ether, acetone, DMSO, and water solvent. Molecular electrostatic potential (MEP), which are used to predict the nucleophile active sites, electron densities calculated at the bond critical points ([Formula: see text] and ring critical points ([Formula: see text] by the atoms in molecules (AIM) method and the donor–acceptor interaction energies ([Formula: see text] calculated using the natural bond orbital (NBO) method were used to investigate the interplay between [Formula: see text]–[Formula: see text] stacking and complexation. The results show that the IS1/TS1-Mg[Formula: see text] compounds have a more protective role compared to TS1/IS1-Mg[Formula: see text] compounds due to the stronger interaction with PHE of A[Formula: see text] in antiaggregation for further development of A[Formula: see text] inhibitors to prevent and disaggregate amyloid formation. Complexation with Mg[Formula: see text] increases the interaction diterpenoid drugs with PHE and makes notable changes in structural and electronic properties of diterpenoids. Also, the interactions between diterpenoid and PHE in less polar environment are more than other environments. Low-polarity environments have the best mimics of the A[Formula: see text] binding site.

Discovering the stacking landscape of a pyridine-pyridine system

Extremely extensive calculations of potential energy surfaces for the parallel-displaced configuration of pyridine dimer systems have been carried out using a dispersion-corrected density functional. Instead of focusing on stationary geometries these calculations provide much deeper insight into the "landscape" of the interaction energies of the particular systems-one can learn how the pyridine dimer stability changes along with various geometrical parameters. Other calculations such as natural bond orbital and energy decomposition have also been applied. The interplay of two significant factors, electrostatic forces and electron correlation effects, have been evaluated. The role of π···π interactions in the stacked pyridine systems has also been confirmed, and surprisingly, this happened to be true even for the geometries where the formation of C-H···π interactions might be proposed instead. The combination of many different methods has revealed the complexity of the stacking interactions. Apart from providing a "literal new look" into pyridine interaction patterns another picture has emerged. A stacking interaction in a pyridine dimer system is perceived as a combination of many different sources of the interaction energy, including orbital ones, and this is true for many different geometries.

Noncovalent Interactions in the Catechol Dimer

Noncovalent interactions play a significant role in a wide variety of biological processes and bio-inspired species. It is, therefore, important to have at hand suitable computational methods for their investigation. In this paper, we report on the contribution of dispersion and hydrogen bonds in both stacked and T-shaped catechol dimers, with the aim of delineating the respective role of these classes of interactions in determining the most stable structure. By using second-order Møller⁻Plesset (MP2) calculations with a small basis set, specifically optimized for these species, we have explored a number of significant sections of the interaction potential energy surface and found the most stable structures for the dimer, in good agreement with the highly accurate, but computationally more expensive coupled cluster single and double excitation and the perturbative triples (CCSD(T))/CBS) method.

Cation–π interactions in CREBBP bromodomain inhibition: an electrostatic model for small-molecule binding affinity and selectivity

A new model is presented to explain and predict binding affinity of aromatic and heteroaromatic ligands for the CREBBP bromodomain based on cation–π interaction strength.