Small molecule conformational preferences derived from crystal structure data. A medicinal chemistry focused analysis

Binding sites and design strategies for small molecule GLP-1R agonists

Glucagon-like peptide-1 receptor (GLP-1R) is a pivotal receptor involved in blood glucose regulation and influencing feeding behavior. It has received significant attention in the treatment of obesity and diabetes due to its potent incretin effect. Peptide GLP-1 receptor agonists (GLP-1RAs) have achieved tremendous success in the market, driving the vigorous development of small molecule GLP-1RAs. Currently, several small molecules have entered the clinical research stage. Additionally, recent discoveries of GLP-1R positive allosteric modulators (PAMs) are also unveiling new regulatory patterns and treatment methods. This article reviews the structure and functional mechanisms of GLP-1R, recent reports on small molecule GLP-1RAs and PAMs, as well as the optimization process. Furthermore, it combines computer simulations to analyze structure-activity relationships (SAR) studies, providing a foundation for exploring new strategies for designing small molecule GLP-1RAs.

Understanding the Cis-Trans Amide Bond Isomerization of N,N'-Diacylhydrazines to Develop Guidelines for A Priori Prediction of Their Most Stable Solution Conformers

N,N'-diacylhydrazines (R1CO-NR3-NR4-COR2) are a class of small molecules with a wide range of applications in chemistry and biology. They are structurally unique in the sense that their two amide groups are connected via a N-N single bond, and as a result, these molecules can exist in eight different isomeric forms. Four of these are amide isomers [trans-trans (t-t), trans-cis (t-c), cis-trans (c-t), and cis-cis (c-c)] arising from C-N bond restricted rotation. In addition, each of these amide isomers can exist in two different isomeric forms due to N-N bond restricted rotation, especially when R3 and R4 groups are relatively bigger. Herein, we have systematically investigated the conformations of 55 N,N'-diacylhydrazines using a combination of solution NMR spectroscopy, X-ray crystallography, and density functional theory calculations. Our data suggest that when the substituents R3 and R4 on the nitrogen atoms are both hydrogens. These molecules prefer twisted trans-trans (t-t) (>90%) geometries (H-N-C═O ∼ 180°), whereas the N-alkylated and N,N'-dialkylated molecules prefer twisted trans-cis (t-c) geometries. Herein, we have analyzed the stabilization of the various isomers of these molecules in light of steric and stereoelectronic effects. We provide a guideline to a priori predict the most stable conformers of the N,N'-diacylhydrazines just by examining their substituents (R1-R4).

Synthetic techniques for thermodynamically disfavoured substituted six-membered rings

Six-membered rings are ubiquitous structural motifs in bioactive compounds and multifunctional materials. Notably, their thermodynamically disfavoured isomers, like disubstituted cyclohexanes featuring one substituent in an equatorial position and the other in an axial position, often exhibit enhanced physical and biological activities in comparison with their opposite isomers. However, the synthesis of thermodynamically disfavoured isomers is, by its nature, challenging, with only a limited number of possible approaches. In this Review, we summarize and compare synthetic methodologies that produce substituted six-membered rings with thermodynamically disfavoured substitution patterns. We place particular emphasis on elucidating the crucial stereoinduction factors within each transformation. Our aim is to stimulate interest in the synthesis of these unique structures, while simultaneously providing synthetic chemists with a guide to approaching this synthetic challenge.

Conformational isomerism in trans-3-methoxycinnamic acid: From solid to gas phase

The rotational spectrum of laser ablated trans-3-methoxycinnamic acid has been observed in the 2-8 GHz range using chirped-pulse Fourier transform microwave spectroscopy coupled to a supersonic jet and adapted to support a laser ablation vaporization system (LA-CP-FTMW). Eight stable conformers were theoretically predicted to exist at B3LYP-D3BJ/6-311++(2d,p) level, all of which were experimentally detected. The experimental rotational parameters data evidence the essentially planar structures for all the conformers. The relative population distribution of conformers in the supersonic jet was investigated from relative intensity measurements. Cooling in the jet brings rotational temperatures close to 1 K for all the conformers. The theoretical predictions for the rotational constants and electric dipole moments show good agreement with the experimental constants and selection rules observed. The population distribution of conformers in the supersonic jet was found to be close to the equilibrium distribution calculated at temperatures lower than the stagnation temperature. Finally, the correlation of the observed conformers structures with those found in condensed phases was investigated.

The Science and Art of Structure-Based Virtual Screening

Structure-based virtual screening has gained momentum again as the high attrition rate at every stage of drug discovery drives the need to explore a greater chemical space. From the Bayesian perspective, its shortcomings as a viable strategy for sustainable hit discovery are discussed, with regard to the prior hit rates of screening libraries and the performance of computational methods. Lessons are shared in selecting virtual hits for experimental validation learned from a series of eight successful campaigns, one of which impacted the discovery of a drug candidate currently in clinical trials.

Optimization Efforts for Identification of Novel Highly Potent Keap1-Nrf2 Protein-Protein Interaction Inhibitors

In research focused on protein-protein interaction (PPI) inhibitors, the optimization process to achieve both high inhibitory activity and favorable physicochemical properties remains challenging. Our previous study reported the discovery of novel and bioavailable Keap1-Nrf2 PPI inhibitor 8 which exhibited moderate in vivo activity in rats. In this work, we present our subsequent efforts to optimize this compound. Two distinct approaches were employed, targeting high energy water molecules and Ser602 as "hot spots" from the anchor with good aqueous solubility, metabolic stability, and membrane permeability. Through ligand efficiency (LE)-guided exploration, we identified two novel inhibitors 22 and 33 with good pharmacokinetics (PK) profiles and more potent in vivo activities, which appear to be promising chemical probes among the existing inhibitors.

Discovery of Potent, Orally Bioavailable, Tricyclic NLRP3 Inhibitors

NLRP3 is a molecular sensor recognizing a wide range of danger signals. Its activation leads to the assembly of an inflammasome that allows for activation of caspase-1 and subsequent maturation of IL-1β and IL-18, as well as cleavage of Gasdermin-d and pyroptotic cell death. The NLRP3 inflammasome has been implicated in a plethora of diseases including gout, type 2 diabetes, atherosclerosis, Alzheimer's disease, and cancer. In this publication, we describe the discovery of a novel, tricyclic, NLRP3-binding scaffold by high-throughput screening. The hit (1) could be optimized into an advanced compound NP3-562 demonstrating excellent potency in human whole blood and full inhibition of IL-1β release in a mouse acute peritonitis model at 30 mg/kg po dose. An X-ray structure of NP3-562 bound to the NLRP3 NACHT domain revealed a unique binding mode as compared to the known sulfonylurea-based inhibitors. In addition, NP3-562 shows also a good overall development profile.

Crystal Clear: Decoding Isocyanide Intermolecular Interactions through Crystallography

The isocyanide group is the chameleon among the functional groups in organic chemistry. Unlike other multiatom functional groups, where the electrophilic and nucleophilic moieties are typically separated, isocyanides combine both functionalities in the terminal carbon. This unique feature can be rationalized using the frontier orbital concept and has significant implications for its intermolecular interactions and the reactivity of the functional group. In this study, we perform a Cambridge Crystallographic Database-supported analysis of isocyanide intramolecular interactions to investigate the intramolecular interactions of isocyanides in the solid state, excluding isocyanide-metal complexes. We discuss examples of different interaction classes, including the isocyanide as a hydrogen bond acceptor (RNC···HX), halogen bonding (RNC···X), and interactions involving the isocyanide and carbon atoms (RNC···C). The latter interaction serves as an intriguing illustration of a Bürgi-Dunitz trajectory and represents a crucial experimental detail in the well-known multicomponent reactions such as the Ugi- and Passerini-type mechanisms. Understanding the spectrum of intramolecular interactions that isocyanides can undergo holds significant implications in fields such as medicinal chemistry, materials science, and asymmetric catalysis.

Conformational energies of reference organic molecules: benchmarking of common efficient computational methods against coupled cluster theory

We selected 145 reference organic molecules that include model fragments used in computer-aided drug design. We calculated 158 conformational energies and barriers using force fields, with wide applicability in commercial and free softwares and extensive application on the calculation of conformational energies of organic molecules, e.g. the UFF and DREIDING force fields, the Allinger's force fields MM3-96, MM3-00, MM4-8, the MM2-91 clones MMX and MM+, the MMFF94 force field, MM4, ab initio Hartree-Fock (HF) theory with different basis sets, the standard density functional theory B3LYP, the second-order post-HF MP2 theory and the Domain-based Local Pair Natural Orbital Coupled Cluster DLPNO-CCSD(T) theory, with the latter used for accurate reference values. The data set of the organic molecules includes hydrocarbons, haloalkanes, conjugated compounds, and oxygen-, nitrogen-, phosphorus- and sulphur-containing compounds. We reviewed in detail the conformational aspects of these model organic molecules providing the current understanding of the steric and electronic factors that determine the stability of low energy conformers and the literature including previous experimental observations and calculated findings. While progress on the computer hardware allows the calculations of thousands of conformations for later use in drug design projects, this study is an update from previous classical studies that used, as reference values, experimental ones using a variety of methods and different environments. The lowest mean error against the DLPNO-CCSD(T) reference was calculated for MP2 (0.35 kcal mol-1), followed by B3LYP (0.69 kcal mol-1) and the HF theories (0.81-1.0 kcal mol-1). As regards the force fields, the lowest errors were observed for the Allinger's force fields MM3-00 (1.28 kcal mol-1), ΜΜ3-96 (1.40 kcal mol-1) and the Halgren's MMFF94 force field (1.30 kcal mol-1) and then for the MM2-91 clones MMX (1.77 kcal mol-1) and MM+ (2.01 kcal mol-1) and MM4 (2.05 kcal mol-1). The DREIDING (3.63 kcal mol-1) and UFF (3.77 kcal mol-1) force fields have the lowest performance. These model organic molecules we used are often present as fragments in drug-like molecules. The values calculated using DLPNO-CCSD(T) make up a valuable data set for further comparisons and for improved force field parameterization.

Moving out of CF3 -Land: Synthesis, Receptor Affinity, and in silico Studies of NK1 Receptor Ligands Containing a Pentafluorosulfanyl (SF5 ) Group

The NK1 receptor (NK1R) is a molecular target for both approved and experimental drugs intended for a variety of conditions, including emesis, pain, and cancers. While contemplating modifications to the typical NK1R pharmacophore, we wondered whether the CF3 groups common for many NK1R ligands, could be replaced with some other moiety. Our attention was drawn by the SF5 group, and so we designed, synthesized, and tested ten novel SF5 -containing compounds for NK1R affinity. All analogues exhibit detectable NK1R binding, with the best of them, compound 5 a, (3-bromo-5-(pentafluoro-λ6 -sulfanyl)benzyl acetyl-L-tryptophanate) binding only slightly worse (IC50 =34.3 nM) than the approved NK1R-targeting drug, aprepitant (IC50 =27.7 nM). Molecular docking provided structural explanation of SAR. According to our analysis, the SF5 group in our compounds occupies a position similar to that of one of the CF3 groups of aprepitant as found in the crystal structure. Additionally, we checked whether the docking scoring function or energies derived from Fragment Molecular Orbital quantum chemical calculations may be helpful in explaining and predicting the experimental receptor affinities for our analogues. Both these methods produce moderately good results. Overall, this is the first demonstration of the utility of the SF5 group in the design of NK1R ligands.

Molecular replacement for small-molecule crystal structure determination from X-ray and electron diffraction data with reduced resolution

The resolution of 3D electron diffraction (ED) data of small-molecule crystals is often relatively poor, due to either electron-beam radiation damage during data collection or poor crystallinity of the material. Direct methods, used as standard for crystal structure determination, are not applicable when the data resolution falls below the commonly accepted limit of 1.2 Å. Therefore an evaluation was carried out of the performance of molecular replacement (MR) procedures, regularly used for protein structure determination, for structure analysis of small-molecule crystal structures from 3D ED data. In the course of this study, two crystal structures of Bi-3812, a highly potent inhibitor of the oncogenic transcription factor BCL6, were determined: the structure of α-Bi-3812 was determined from single-crystal X-ray data, the structure of β-Bi-3812 from 3D ED data, using direct methods in both cases. These data were subsequently used for MR with different data types, varying the data resolution limit (1, 1.5 and 2 Å) and by using search models consisting of connected or disconnected fragments of BI-3812. MR was successful with 3D ED data at 2 Å resolution using a search model that represented 74% of the complete molecule.

A quantum chemical interaction energy dataset for accurately modeling protein-ligand interactions

Fast and accurate calculation of intermolecular interaction energies is desirable for understanding many chemical and biological processes, including the binding of small molecules to proteins. The Splinter ["Symmetry-adapted perturbation theory (SAPT0) protein-ligand interaction"] dataset has been created to facilitate the development and improvement of methods for performing such calculations. Molecular fragments representing commonly found substructures in proteins and small-molecule ligands were paired into >9000 unique dimers, assembled into numerous configurations using an approach designed to adequately cover the breadth of the dimers' potential energy surfaces while enhancing sampling in favorable regions. ~1.5 million configurations of these dimers were randomly generated, and a structurally diverse subset of these were minimized to obtain an additional ~80 thousand local and global minima. For all >1.6 million configurations, SAPT0 calculations were performed with two basis sets to complete the dataset. It is expected that Splinter will be a useful benchmark dataset for training and testing various methods for the calculation of intermolecular interaction energies.

The Role of Allylic Strain for Conformational Control in Medicinal Chemistry

It is axiomatic in medicinal chemistry that optimization of the potency of a small molecule at a macromolecular target requires complementarity between the ligand and target. In order to minimize the conformational penalty on binding, both enthalpically and entropically, it is therefore preferred to have the ligand preorganized in the bound conformation. In this Perspective, we highlight the role of allylic strain in controlling conformational preferences. Allylic strain was originally described for carbon-based allylic systems, but the same principles apply to other types of structure with sp² or pseudo-sp² arrangements. These systems include benzylic (including heteroaryl methyl) positions, amides, N-aryl groups, aryl ethers, and nucleotides. We have derived torsion profiles from small molecule X-ray structures for these systems. Through multiple examples, we show how these effects have been applied in drug discovery and how they can be used prospectively to influence conformation in the design process.

Synthesis of 2,3-Disubstituted Pyrroles by Lewis Acid Promoted Cyclization of N-Sulfonyl Vinylogous Carbamates and Amides

A three-step sequence to construct pyrroles from three components including 2,2-dimethoxyethylamine, aryl/alkyl sulfonyl chlorides and alkynes bearing electron-withdrawing groups is presented. The pyrroles are proposed to arise via a 5-exo-trig cyclization proceeding through both oxocarbenium and N-sulfonyliminium ions. This modular route allows for the variability at the N-sulfonyl group, the C2 and C3 substituents for rational vectoring of the pyrrole nucleus for downstream processes.

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

A rod- and tessellation-based comparative analysis of polymorphic and structurally-invariant molecular crystals: application to sulfathiazole and 2-benzyl-5-benzylidenecyclopentanones

A rationalization of the alternative crystal structures adopted by a given molecular compound or by a set of substitutionally related molecular compounds is provided by reference to the five known polymorphs of sulfathiazole and 16 substituted 2-benzyl-5-benzylidene cyclopentanones (BBCPs), respectively. Two-dimensional (2D) packing fractions (ϕ_2D) take space-group symmetry into account, with a clear demarcation of closed-packed zones (CPZ) and molecular junction zones (JZ) in all Z' = 1 structures. Representation of the molecules as two linked rods allows a concise treatment of conformation and rapid visualization of crystal packing. Combined with calculations of intermolecular potential energies, the rod method provides insight into the stabilization mechanisms of alternative polymorphs. In sulfathiazole, the primary factor is to obtain satisfactory hydrogen bonding, with close packing a secondary consideration. In BBCP derivatives, by comparison, close packing is the primary mechanism of stabilization. Whereas the 2D structures arising in CPZ can be analysed as tessellations of molecular-based cells, a method based on 2D Dirichlet cells is required for the JZ. These are calculated from the centroids of the molecular envelopes in high-symmetry planes. It is shown that these centroid coordinates, when combined with space-group symmetry and unit cell coordinates, provide a concise parameterization of all structures containing JZ. It is anticipated that this parameterization may be exploited to predict such crystal structures from powder diffraction data.

Graph-Convolutional Neural Net Model of the Statistical Torsion Profiles for Small Organic Molecules

We present a graph-convolutional neural network (GCNN)-based method for learning and prediction of statistical torsional profiles (STP) in small organic molecules based on the experimental X-ray structure data. A specialized GCNN torsion profile model is trained using the structures in the Crystallography Open Database (COD). The GCNN-STP model captures torsional preferences over a wide range of torsion rotor chemotypes and correctly predicts a variety of effects from the vicinal atoms and moieties. GCNN-STP statistical profiles also show good agreement with quantum chemically (DFT) calculated torsion energy profiles. Furthermore, we demonstrate the application of the GCNN-STP statistical profiles for conformer generation. A web server that allows interactive profile prediction and viewing is made freely available at https://www.molsoft.com/tortool.html.

Auto3D: Automatic Generation of the Low-Energy 3D Structures with ANI Neural Network Potentials

Computational programs accelerate the chemical discovery processes but often need proper three-dimensional molecular information as part of the input. Getting optimal molecular structures is challenging because it requires enumerating and optimizing a huge space of stereoisomers and conformers. We developed the Python-based Auto3D package for generating the low-energy 3D structures using SMILES as the input. Auto3D is based on state-of-the-art algorithms and can automatize the isomer enumeration and duplicate filtering process, 3D building process, geometry optimization, and ranking process. Tested on 50 molecules with multiple unspecified stereocenters, Auto3D is guaranteed to find the stereoconfiguration that yields the lowest-energy conformer. With Auto3D, we provide an extension of the ANI model. The new model, dubbed ANI-2xt, is trained on a tautomer-rich data set. ANI-2xt is benchmarked with DFT methods on geometry optimization and electronic and Gibbs free energy calculations. Compared with ANI-2x, ANI-2xt provides a 42% error reduction for tautomeric reaction energy calculations when using the gold-standard coupled-cluster calculation as the reference. ANI-2xt can accurately predict the energies and is several orders of magnitude faster than DFT methods.

Investigation of tadalafil molecular arrangement in solid dispersions using inverse gas chromatography and Raman mapping

The aim of this study was to investigate the molecular structures of tadalafil solid dispersions prepared by different techniques and further to relate them to surface free energy information indicating the final amorphousness of the product. Thus, we tried to complement the existing knowledge of solid dispersion formation. Poorly water-soluble tadalafil was combined with different polymers, i.e. Kollidon® 12 PF, Kollidon® VA 64 and Soluplus®, to form model systems. To assess the extent of drug-polymer miscibility, we studied model solid dispersion surface energy using inverse gas chromatography and phase micro-structure using confocal Raman microscopy. The selection of the preparation method was found to play a crucial role in the molecular arrangement of the incorporated drug and the polymer in resulting solid dispersion. Our results showed that a lower surface free energy indicated the formation of a more homogeneous solid dispersion. Conversely, a higher surface free energy corresponded to the heterogeneous systems containing tadalafil amorphous clusters that were captured by Raman mapping. Thus, we successfully introduced a novel evaluation approach of the drug molecular arrangement in solid dispersions that is especially useful for examining the miscibility of the components when the conventional characterizing techniques are inconclusive or yield variable results.

A Small-Molecule Oral Agonist of the Human Glucagon-like Peptide-1 Receptor

Peptide agonists of the glucagon-like peptide-1 receptor (GLP-1R) have revolutionized diabetes therapy, but their use has been limited because they require injection. Herein, we describe the discovery of the orally bioavailable, small-molecule, GLP-1R agonist PF-06882961 (danuglipron). A sensitized high-throughput screen was used to identify 5-fluoropyrimidine-based GLP-1R agonists that were optimized to promote endogenous GLP-1R signaling with nanomolar potency. Incorporation of a carboxylic acid moiety provided considerable GLP-1R potency gains with improved off-target pharmacology and reduced metabolic clearance, ultimately resulting in the identification of danuglipron. Danuglipron increased insulin levels in primates but not rodents, which was explained by receptor mutagensis studies and a cryogenic electron microscope structure that revealed a binding pocket requiring a primate-specific tryptophan 33 residue. Oral administration of danuglipron to healthy humans produced dose-proportional increases in systemic exposure (NCT03309241). This opens an opportunity for oral small-molecule therapies that target the well-validated GLP-1R for metabolic health.

The Shapes of Sulfonamides: A Rotational Spectroscopy Study

Benzenesulfonamides are a class of molecules of extreme interest in the biochemical field because many of them are active against a variety of diseases. In this work, the pharmacophoric group benzensulfonamide, its derivatives para-toluensulfonamide and ortho-toluensulfonamide, and the bioactive molecule sulfanilamide, were investigated using rotational spectroscopy to determine their conformations and the influence of different substituents on their structures. For all species, the hyperfine structure due to the ¹⁴N atom was analyzed, and this provided crucial information for the unambiguous identification of the observed conformation of all molecules. In addition, for ortho-toluensulfonamide, the vibration–rotation hyperfine structure related to the methyl torsion was analyzed, and the methyl group rotation barrier was determined. For benzensulfonamide, partial r_S and r₀ structures were established from the experimental rotational constants of the parent and two deuterated isotopic species. In all compounds except ortho-toluensulfonamide, the amino group of the sulfonamide group lies perpendicular to the benzene plane with the aminic hydrogens eclipsing the oxygen atoms. In ortho-toluensulfonamide, where weak attractive interactions occur between the nitrogen lone pair and the methyl hydrogen atoms, the amino group lies in a gauche orientation, retaining the eclipsed configuration with respect to the SO₂ frame. A comparison of the geometrical arrangements found in the PDB database allowed us to understand that the bioactive conformations are different from those found in isolated conditions. The conformations within the receptor are reached with an energy cost, which is balanced by the interactions established in the receptor.

Computationally driven discovery of SARS-CoV-2 Mpro inhibitors: from design to experimental validation

We report a fast-track computationally driven discovery of new SARS-CoV-2 main protease (M^pro) inhibitors whose potency ranges from mM for the initial non-covalent ligands to sub-μM for the final covalent compound (IC₅₀ = 830 ± 50 nM). The project extensively relied on high-resolution all-atom molecular dynamics simulations and absolute binding free energy calculations performed using the polarizable AMOEBA force field. The study is complemented by extensive adaptive sampling simulations that are used to rationalize the different ligand binding poses through the explicit reconstruction of the ligand–protein conformation space. Machine learning predictions are also performed to predict selected compound properties. While simulations extensively use high performance computing to strongly reduce the time-to-solution, they were systematically coupled to nuclear magnetic resonance experiments to drive synthesis and for in vitro characterization of compounds. Such a study highlights the power of in silico strategies that rely on structure-based approaches for drug design and allows the protein conformational multiplicity problem to be addressed. The proposed fluorinated tetrahydroquinolines open routes for further optimization of M^pro inhibitors towards low nM affinities.

The dominant binding mode of the QUB-00006-Int-07 main protease inhibitor during absolute binding free energy simulations.

Modulating Conformational Preferences by Allylic Strain toward Improved Physical Properties and Binding Interactions

The preference of the axial over the equatorial orientation of 2-substitutent for both phenyl-1-piperidines and N-acylpiperidines is studied at the M06-2X level of theory. For phenyl-1-piperidines, the axial 2-substituent is modestly favored over the equatorial one. In contrast, the pseudoallylic strain in N-acylpiperidines dictates the axial orientation of 2-substituent with a ΔG up to -3.2 kcal/mol. The calculations agree well with the statistics from both the Cambridge Structural Database of small-molecule organic crystal structures and the Protein Data Bank. The equilibrium between the twist-boat and chair conformations for N-acylpiperidines with a 2-substituent was further investigated. The twist-boat conformation is found to be around 1.5 kcal/mol less favorable. Finally, the three-dimensionality in shape resulting from minimization of the pseudoallylic strain is characterized, and its implication in protein-ligand interactions is briefly reviewed.

Discovery of a Partial Glucokinase Activator Clinical Candidate: Diethyl ((3-(3-((5-(Azetidine-1-carbonyl)pyrazin-2-yl)oxy)-5-isopropoxybenzamido)-1H-pyrazol-1-yl)methyl)phosphonate (BMS-820132)

Glucokinase (GK) is a key regulator of glucose homeostasis, and its small-molecule activators represent a promising opportunity for the treatment of type 2 diabetes. Several GK activators have been advanced into clinical trials and have demonstrated promising efficacy; however, hypoglycemia represents a key risk for this mechanism. In an effort to mitigate this hypoglycemia risk while maintaining the efficacy of the GK mechanism, we have investigated a series of amino heteroaryl phosphonate benzamides as ''partial" GK activators. The structure-activity relationship studies starting from a "full GK activator" 11, which culminated in the discovery of the "partial GK activator" 31 (BMS-820132), are discussed. The synthesis and in vitro and in vivo preclinical pharmacology profiles of 31 and its pharmacokinetics (PK) are described. Based on its promising in vivo efficacy and preclinical ADME and safety profiles, 31 was advanced into human clinical trials.

TorsionNet: A Deep Neural Network to Rapidly Predict Small-Molecule Torsional Energy Profiles with the Accuracy of Quantum Mechanics

Fast and accurate assessment of small-molecule dihedral energetics is crucial for molecular design and optimization in medicinal chemistry. Yet, accurate prediction of torsion energy profiles remains challenging as the current molecular mechanics (MM) methods are limited by insufficient coverage of drug-like chemical space and accurate quantum mechanical (QM) methods are too expensive. To address this limitation, we introduce TorsionNet, a deep neural network (DNN) model specifically developed to predict small-molecule torsion energy profiles with QM-level accuracy. We applied active learning to identify nearly 50k fragments (with elements H, C, N, O, F, S, and Cl) that maximized the coverage of our corporate compound library and leveraged massively parallel cloud computing resources for density functional theory (DFT) torsion scans of these fragments, generating a training data set of 1.2 million DFT energies. After training TorsionNet on this data set, we obtain a model that can rapidly predict the torsion energy profile of typical drug-like fragments with DFT-level accuracy. Importantly, our method also provides an uncertainty estimate for the predicted profiles without any additional calculations. In this report, we show that TorsionNet can accurately identify the preferred dihedral geometries observed in crystal structures. Our TorsionNet-based analysis of a diverse set of protein-ligand complexes with measured binding affinity shows a strong association between high ligand strain and low potency. We also present practical applications of TorsionNet that demonstrate how consideration of DNN-based strain energy leads to substantial improvement in existing lead discovery and design workflows. TorsionNet500, a benchmark data set comprising 500 chemically diverse fragments with DFT torsion profiles (12k MM- and DFT-optimized geometries and energies), has been created and is made publicly available.

Recent applications of seven-membered rings in drug design

This short review aims at highlighting recent design strategies hinged on using seven-membered rings. Analyses of the different selected examples coupled with torsion profiles derived from the CCDC suggest some of these strategies could have broad applications.

Control of Molecular Conformation and Crystal Packing of Biphenyl Derivatives

This Review presents a discussion of the conformation of biphenyl derivatives in different chemical environments. The interplay between aromatic stabilization and steric repulsion, normally considered to explain the conformation of the molecule, is contrasted with the interpretation provided by models not based on molecular orbitals. The electronic control of conformation by means of appropriate hydrogen substitution is discussed by examples taken from chemistry and molecular electronics. Supramolecular synthons involving biphenyl are critically analyzed in terms of the molecular conformation, crystal packing and intermolecular forces. Some directions for future research on the control of the conformation of biphenyls are also presented.

Long-Timescale Simulations Revealed Critical Non-Conserved Residues of Phosphodiesterases Affecting Selectivity of BAY60-7550

A major obstacle of the selective inhibitor design for specific human phosphodiesterase (PDE) is that highly conserved catalytic pockets are difficult to be distinguished by inhibitor molecules. To overcome this, a feasible path is to understand the molecular determinants underlying the selectivity of current inhibitors. BAY60-7550 (BAY for short; IC₅₀ = 4.7 nM) is a highly selective inhibitor targeting PDE2A which is a dual-specificity PDE and an attractive target for therapeutic intervention of the central nervous system (CNS) disorders. Recent studies suggest that molecular determinants may be in binding processes of BAY. However, a detailed understanding of these processes are still lacking. To explore these processes, High-Throughput Molecular Dynamics (HTMD) simulations were performed to reproduce the spontaneous association of BAY with catalytic pockets of 4 PDE isoforms; Ligand Gaussian Accelerated Molecular Dynamics (LiGaMD) simulations were performed to reproduce the unbinding-rebinding processes of FKG and 10.13039/100016266MC2, two pyrazolopyrimidinone PDE2A selective inhibitors, in the PDE2A system. The produced molecular trajectories were analyzed by the Markov state model (MSM) and the molecular mechanics/generalized Born surface area (MM/GBSA). The results showed that the non-covalent interactions between the non-conserved residues and BAY, especially the hydrogen bonds, determined the unique binding pathways of BAY on the surface of PDE2A. These pathways were different from those of BAY on the surface of the other three PDE isoforms and the binding pathways of the other two PDE2A inhibitors in PDE2A systems. These differences were ultimately reflected in the high selectivity of this inhibitor for PDE2A. As a result, this study demonstrates the critical role of the binding processes in the selectivity of BAY, and also identifies the key non-conserved residues affecting the binding processes of BAY. Thus, this study provides a new perspective and data support for the further development of BAY-derived inhibitors targeting PDE2A.

Unravelling conformational and crystal packing preferences of cyclohexane-5-spirohydantoin derivatives incorporating a halogenated benzoyl group

The conformational and supramolecular diversity of spirohydantoins have been investigated.

Navigating complex peptide structures using macrocycle conformational maps

Identification of turn motifs that are stabilized by intramolecular hydrogen bonds can be useful in describing the conformation of peptide systems. However, this approach is somewhat insufficient for cyclic peptides because peptide regions that are not positioned within a hydrogen bond can be left with no description. Furthermore, non-regular secondary structures and other rarely-observed conformations can be left without detailed evaluation. Herein, we describe “higher-order” ϕ/ψ plots termed macrocycle conformational maps (MCMs) as a tool for evaluating and comparing the conformations of a series of structurally related macrocyclic peptides.

Identification of turn motifs that are stabilized by hydrogen bonds can be useful in describing the conformation of peptides. Herein, we describe “higher-order” ϕ/ψ plots termed macrocycle conformational maps (MCMs) as a tool to evaluate and compare the conformations of related macrocycles.

The reorganization energy of compounds upon binding to proteins, from dynamic and solvated bound and unbound states

The intramolecular reorganization energy (ΔE_Reorg) of compounds upon binding to proteins is a component of the binding free energy, which has long received particular attention, for fundamental and practical reasons. Understanding ΔE_Reorg would benefit the science of molecular recognition and drug design. For instance, the tolerable strain energy of compounds upon binding has been elusive. Prior studies found some large ΔE_Reorg values (e.g. > 10 kcal/mol), received with skepticism since they imply excessive opposition to binding. Indeed, estimating ΔE_Reorg is technically difficult. Typically, ΔE_Reorg has been approached by taking two energy-minimized conformers representing the bound and unbound states, and subtracting their conformational energy. This is a drastic oversimplification, liable to conformational collapse of the unbound conformer. Instead, the present work applies extensive molecular dynamics (MD) and the modern OPLS3 force-field to simulate compounds bound and unbound states, in explicit solvent under physically relevant conditions. The thermalized unbound compounds populate multiple conformations, not reducible to one or a few energy-minimized conformers. The intramolecular energies in the bound and unbound states were averaged over pertinent conformational ensembles, and the reorganization enthalpy upon binding (ΔH_Reorg) deduced by subtraction. This was applied to 76 systems, including 43 approved drugs, carefully selected for i) the quality of the bioactive X-ray structures and ii) the diversity of the chemotypes, their properties and protein targets. It yielded comparatively low ΔH_Reorg values (median = 1.4 kcal/mol, mean = 3.0 kcal/mol). A new finding is the observation of negative ΔH_Reorg values. Indeed, reorganization energies do not have to oppose binding, e.g. when intramolecular interactions stabilize preferentially the bound state. Conversely, even with competing water molecules, intramolecular interactions can occur predominantly in the unbound compound, and be replaced by intermolecular counterparts upon protein binding. Such disruption of intramolecular interactions upon binding gives rise to occasional larger ΔH_Reorg values. Such counterintuitive larger ΔH_Reorg values may be rationalized as a redistribution of interactions upon binding, qualitatively compatible with binding.

Comprehensive Strategies to Bicyclic Prolines: Applications in the Synthesis of Potent Arginase Inhibitors

Comprehensive synthetic strategies afforded a diverse set of structurally unique bicyclic proline-containing arginase inhibitors with a high degree of three-dimensionality. The analogs that favored the Cγ-exo conformation of the proline improved the arginase potency over the initial lead. The novel synthetic strategies reported here not only enable access to previously unknown stereochemically complex proline derivatives but also provide a foundation for the future synthesis of bicyclic proline analogs, which incorporate inherent three-dimensional character into building blocks, medicine, and catalysts and could have a profound impact on the conformation of proline-containing peptides and macrocycles.

Photomediated ring contraction of saturated heterocycles

Saturated heterocycles are found in numerous therapeutics and bioactive natural products and are abundant in many medicinal and agrochemical compound libraries. To access new chemical space and function, many methods for functionalization on the periphery of these structures have been developed. Comparatively fewer methods are known for restructuring their core framework. Herein, we describe a visible light-mediated ring contraction of α-acylated saturated heterocycles. This unconventional transformation is orthogonal to traditional ring contractions, challenging the paradigm for diversification of heterocycles including piperidine, morpholine, thiane, tetrahydropyran, and tetrahydroisoquinoline derivatives. The success of this Norrish type II variant rests on reactivity differences between photoreactive ketone groups in specific chemical environments. This strategy was applied to late-stage remodeling of pharmaceutical derivatives, peptides, and sugars.

Enantioselective Cyclobutenylation of Olefins Using N-Sulfonyl-1,2,3-Triazoles as Vicinal Dicarbene Equivalents

Cyclobutenes are highly useful synthetic intermediates as well as important motifs in bioactive small molecules. Herein, we report a regio-, chemo-, and enantioselective synthesis of cyclobutenes from olefins using N-sulfonyl-1,2,3-triazoles as vicinal dicarbene equivalents or alkyne [2 + 2] cycloaddition surrogates. Terminal and cis-olefins can be transformed into enantioenriched cyclopropanes via rhodium catalysis. Then, in one pot, treatment of these intermediates with tosyl hydrazide and base effects diazo formation followed by rhodium-catalyzed ring expansion to yield enantioenriched cyclobutenes. These cyclobutenes can be transformed into highly substituted, enantioenriched cyclobutanes, including structures relevant to natural product scaffolds.

What We Learned in 25 Years of Interactive Molecular Design Sessions

We retrace Prof. François Diederich's consultancy work for Roche and its impact over the years he worked with us. François Diederich uniquely shaped our approach to molecular design, and interactions with him and his research group at ETH Zurich have created deep insights into molecular recognition. Herein we share how his style and approach continue to inspire us.

Discovery of a novel series of GPR119 agonists: Design, synthesis, and biological evaluation of N-(Piperidin-4-yl)-N-(trifluoromethyl)pyrimidin-4-amine derivatives

We undertook an optimization effort involving propan-2-yl 4-({6-[5-(methanesulfonyl)-2,3-dihydro-1H-indol-1-yl]pyrimidin-4-yl}oxy)piperidine-1-carboxylate 1, which we had previously discovered as a novel G protein-coupled receptor 119 (GPR119) agonist. To occupy a presumed hydrophobic space between the pyrimidine and piperidine rings in interaction with GPR119, we replaced the linker oxygen with nitrogen. Subsequently, the introduction of a substituent at the bridging nitrogen atom was explored. We found that the installation of N-trifluoromethyl group 10 not only enhanced GPR119 agonist activity but also considerably improved the human ether-à-go-go-related gene (hERG) inhibition profile. These improvements were not observed for non-fluorinated substituents, such as ethyl analog 8b. The next optimization effort focused on the exploration of a new surrogate structure for the indoline ring and the isosteric replacements of the piperidine N-Boc group to improve solubility, metabolic stability, and oral bioavailability. As a result, N-{1-[3-(2-fluoropropan-2-yl)-1,2,4-oxadiazol-5-yl]piperidin-4-yl}-6-{[1-(methanesulfonyl)piperidin-4-yl]oxy}-N-(trifluoromethyl)pyrimidin-4-amine (27) was identified as a potent and orally bioavailable GPR119 agonist. This compound augmented insulin secretion and effectively lowered plasma glucose excursion in a diabetic animal model after oral administration. In this study, we discuss the designs, syntheses, and biological activities of a novel series of N-(piperidin-4-yl)-N-(trifluoromethyl)pyrimidin-4-amine derivatives as GPR119 agonists, and to determine the distinctive effect of the N-trifluoromethyl group on hERG inhibition, we also discuss the conformational preference of representative compounds.

An Exploration of Chemical Properties Required for Cooperative Stabilization of the 14-3-3 Interaction with NF-κB-Utilizing a Reversible Covalent Tethering Approach

Protein–protein modulation has emerged as a proven approach to drug discovery. While significant progress has been gained in developing protein–protein interaction (PPI) inhibitors, the orthogonal approach of PPI stabilization lacks established methodologies for drug design. Here, we report the systematic ″bottom-up″ development of a reversible covalent PPI stabilizer. An imine bond was employed to anchor the stabilizer at the interface of the 14-3-3/p65 complex, leading to a molecular glue that elicited an 81-fold increase in complex stabilization. Utilizing protein crystallography and biophysical assays, we deconvoluted how chemical properties of a stabilizer translate to structural changes in the ternary 14-3-3/p65/molecular glue complex. Furthermore, we explore how this leads to high cooperativity and increased stability of the complex.

The Good, the Bad, and the Twisted Revisited: An Analysis of Ligand Geometry in Highly Resolved Protein-Ligand X-ray Structures

An analysis of the rotatable bond geometry of drug-like ligand models is reported for high-resolution (<1.1 Å) crystallographic protein-ligand complexes. In cases where the ligand fit to the electron density is very good, unusual torsional geometry is rare and, most often, though not exclusively, associated with strong polar, metal, or covalent ligand-protein interactions. It is rarely associated with a torsional strain of greater than 2 kcal mol^-1 by calculation. An unusual torsional geometry is more prevalent where the fit to electron density is not perfect. Multiple low-strain conformer bindings were observed in 21% of the set and, it is suggested, may also lie behind many of the 35% of single-occupancy cases, where a poor fit to the e-density was found. It is concluded that multiple conformer ligand binding is an under-recognized phenomenon in structure-based drug design and that there is a need for more robust crystallographic refinement methods to better handle such cases.

Discovery of a novel series of indolinylpyrimidine-based GPR119 agonists: Elimination of ether-a-go-go-related gene liability using a hydrogen bond acceptor-focused approach

We previously identified a novel series of indolinylpyrimidine derivatives exemplified by 2 in Figure 1, which is an indoline based derivative, as potent GPR119 agonists. Despite the attractive potency of 2, this compound inhibited the human ether-a-go-go-related gene (hERG) K⁺ channel. We elucidated crucial roles of the methylsulfonyl group of 2 in its interaction with the hERG channel and the GPR119 receptor, presumably as a hydrogen bond acceptor (HBA). To remove the undesirable hERG inhibitory activity, a strategy was implemented to arrange an HBA on a less conformationally flexible framework at the indoline 5-position instead of the methylsulfonyl group. This successfully led to the discovery of a piperidinone ring as a desirable motif at the indoline 5-position, which could minimize hERG liability as shown by 24b. Further optimization focused on the reduction of lipophilicity in terms of more favorable drug-like properties. Consequently, the introduction of a hydroxy group at the 3-position of the piperidinone ring effectively reduced lipophilicity without compromising GPR119 potency, resulting in the identification of (3S)-3-hydroxy-1-{1-[6-({1-[3-(propan-2-yl)-1,2,4-oxadiazol-5-yl]piperidin-4-yl}oxy)pyrimidin-4-yl]- 2,3-dihydro-1H-indol-5-yl}piperidin-2-one ((S)-29) as a novel, potent, and orally bioavailable GPR119 agonist with a well-balanced profile. The pharmacological effects of this compound were also confirmed after single and chronic oral administration in diabetic animal models.

Force Field Optimization Guided by Small Molecule Crystal Lattice Data Enables Consistent Sub-Angstrom Protein-Ligand Docking

Accurate and rapid calculation of protein-small molecule interaction free energies is critical for computational drug discovery. Because of the large chemical space spanned by drug-like molecules, classical force fields contain thousands of parameters describing atom-pair distance and torsional preferences; each parameter is typically optimized independently on simple representative molecules. Here, we describe a new approach in which small molecule force field parameters are jointly optimized guided by the rich source of information contained within thousands of available small molecule crystal structures. We optimize parameters by requiring that the experimentally determined molecular lattice arrangements have lower energy than all alternative lattice arrangements. Thousands of independent crystal lattice-prediction simulations were run on each of 1386 small molecule crystal structures, and energy function parameters of an implicit solvent energy model were optimized, so native crystal lattice arrangements had the lowest energy. The resulting energy model was implemented in Rosetta, together with a rapid genetic algorithm docking method employing grid-based scoring and receptor flexibility. The success rate of bound structure recapitulation in cross-docking on 1112 complexes was improved by more than 10% over previously published methods, with solutions within <1 Å in over half of the cases. Our results demonstrate that small molecule crystal structures are a rich source of information for guiding molecular force field development, and the improved Rosetta energy function should increase accuracy in a wide range of small molecule structure prediction and design studies.

Illuminating the dark conformational space of macrocycles using dominant rotors

Three-dimensional conformation is the primary determinant of molecular properties. The thermal energy available at room temperature typically equilibrates the accessible conformational states. Here, we introduce a method for isolating unique and previously understudied conformations of macrocycles. The observation of unusual conformations of 16- to 22-membered rings has been made possible by controlling their interconversion using dominant rotors, which represent tunable atropisomeric constituents with relatively high rotational barriers. Density functional theory and in situ NMR measurements suggest that dominant rotor candidates for the amino-acid-based structures considered here should possess a rotational energy barrier of at least 25 kcal mol^-1. Notable differences in the geometries of the macrocycle conformations were identified by NMR spectroscopy and X-ray crystallography. There is evidence that amino acid residues can be forced into rare turn motifs not observed in the corresponding linear counterparts and homodetic rings. These findings should unlock new avenues for studying the conformation-activity relationships of bioactive molecules.

Biaryl sulfonamides as cisoid azosteres for photopharmacology

Biaryl sulfonamides are excellent candidates for the azologization approach that yields photoswitchable drugs more active in their metastable cis state, compared to the stable trans state.

Experimental and theoretical investigation of conformational states and noncovalent interactions in crystalline sulfonamides with a methoxyphenyl moiety

The conformational and noncovalent interaction properties of sulfonamides with a methoxyphenyl moiety were examined by both experimental and theoretical methods.

Discovery of Novel and Highly Selective Cyclopropane ALK Inhibitors through a Fragment-Assisted, Structure-Based Drug Design

Fragment screening is frequently used for hit identification. However, there was no report starting from a small fragment for the development of an anaplastic lymphoma kinase (ALK) inhibitor, despite the number of ALK inhibitors reported. We began our research with the fragment hit F-1 and our subsequent linker design, and its docking analysis yielded novel cis-1,2,2-trisubstituted cyclopropane 1. The fragment information was integrated with a structure-based approach to improve upon the selectivity over tropomyosin receptor kinase A, leading to the potent and highly selective ALK inhibitor, 4-trifluoromethylphenoxy-cis-1,2,2-trisubstituted cyclopropane 12. This work shows that fragments become a powerful tool for both lead generation and optimization, such as the improvement of selectivity, by combining them with a structure-based drug design approach, resulting in the fast and efficient development of a novel, potent, and highly selective compound.