Methyl effects on protein-ligand binding

In Silico Positional Analogue Scanning with Amber GPU-TI

Positional analogue scanning (PAS) is an accepted strategy for multiparameter lead optimization (MPO) in drug discovery. Small structural changes as introduced by PAS can lead to 10-fold changes in binding potency in ∼10-20% of cases, a significant parameter shift irrespective of other MPO objectives. Sometimes performing a complete PAS is challenging due to resource and time constraints, building block availability, or difficulty in synthesis. Calculating relative binding free energies (RBFEs) for all positions can contribute to prioritizing the most promising analogues for synthesis. We tested a well-established RBFE calculation method, Amber GPU-TI, for 20 positional analogue scans in 14 test systems (cyclin-dependent kinase 8 (CDK8), hepatitis C virus nonstructural protein 5B (HCV NS5B), tankyrase, RAC-α serine/threonine-protein kinase (Akt), phosphodiesterase 1B (PDE1B), orexin/hypocretin receptor type 1 (OX1R), orexin/hypocretin receptor type 2 (OX2R), histone acetyltransferase K (lysine) acetyltransferase 6A (KAT6A), peroxisome proliferator-activated receptor γ (PPARγ), extracellular signal-regulated kinases (ERK1/2), coactivator-associated arginine methyltransferase 1 (PRMT4), α_vβ₆, bromodomain 1 (BD1), human immunodeficiency virus-1 (HIV-1) entry) involving nitrogen, methyl, halogen, methoxy, and hydroxyl scans with at least four analogues per set. Among the 66 analogue positions explored, we found that in 18 cases Amber GPU-TI calculations predicted a more than 10-fold change in potency. In all of these cases, the experimentally observed direction of potency changes agreed with the predictions. In 16 cases, more than 10-fold changes in experimental potency were observed. Again, in all of these cases, Amber GPU-TI predicted the direction of the potency changes correctly. In none of these cases would a decision made for or against synthesis based on a 10-fold change in potency have resulted in missing an important analogue. Therefore, in silico RBFE calculations using Amber GPU-TI can meaningfully contribute to the prioritization of positional analogues before synthesis.

Implications of Additivity and Nonadditivity for Machine Learning and Deep Learning Models in Drug Design

Matched molecular pairs (MMPs) are nowadays a commonly applied concept in drug design. They are used in many computational tools for structure-activity relationship analysis, biological activity prediction, or optimization of physicochemical properties. However, until now it has not been shown in a rigorous way that MMPs, that is, changing only one substituent between two molecules, can be predicted with higher accuracy and precision in contrast to any other chemical compound pair. It is expected that any model should be able to predict such a defined change with high accuracy and reasonable precision. In this study, we examine the predictability of four classical properties relevant for drug design ranging from simple physicochemical parameters (log D and solubility) to more complex cell-based ones (permeability and clearance), using different data sets and machine learning algorithms. Our study confirms that additive data are the easiest to predict, which highlights the importance of recognition of nonadditivity events and the challenging complexity of predicting properties in case of scaffold hopping. Despite deep learning being well suited to model nonlinear events, these methods do not seem to be an exception of this observation. Though they are in general performing better than classical machine learning methods, this leaves the field with a still standing challenge.

Surprising lipophilicity observations identify unexpected conformational effects

Contrary to expectation N-aryl pyrrolidinones (and isosteric imidazolinones and oxazolinones) are more lipophilic and less soluble than the corresponding piperidinones (tetrahydropyrimidinones and oxazinones). Exploration of the basis for these results uncovered a subtle interplay of steric and electronic effects that result in different conformations for the two classes of compounds which drive the observed effects.

Discovery of Potent and Orally Bioavailable Pyridine N-Oxide-Based Factor XIa Inhibitors through Exploiting Nonclassical Interactions

Activated factor XI (FXIa) inhibitors are promising novel anticoagulants with low bleeding risk compared with current anticoagulants. The discovery of potent FXIa inhibitors with good oral bioavailability has been challenging. Herein, we describe our discovery effort, utilizing nonclassical interactions to improve potency, cellular permeability, and oral bioavailability by enhancing the binding while reducing polar atoms. Beginning with literature-inspired pyridine N-oxide-based FXIa inhibitor 1, the imidazole linker was first replaced with a pyrazole moiety to establish a polar C-H···water hydrogen-bonding interaction. Then, structure-based drug design was employed to modify lead molecule 2d in the P1' and P2' regions, with substituents interacting with key residues through various nonclassical interactions. As a result, a potent FXIa inhibitor 3f (K_i = 0.17 nM) was discovered. This compound demonstrated oral bioavailability in preclinical species (rat 36.4%, dog 80.5%, and monkey 43.0%) and displayed a dose-dependent antithrombotic effect in a rabbit arteriovenous shunt model of thrombosis.

Enabling Role of Ligand-Driven Conformational Changes in Enzyme Evolution

Many enzymes that show a large specificity in binding the enzymatic transition state with a higher affinity than the substrate utilize substrate binding energy to drive protein conformational changes to form caged substrate complexes. These protein cages provide strong stabilization of enzymatic transition states. Using part of the substrate binding energy to drive the protein conformational change avoids a similar strong stabilization of the Michaelis complex and irreversible ligand binding. A seminal step in the development of modern enzyme catalysts was the evolution of enzymes that couple substrate binding to a conformational change. These include enzymes that function in glycolysis (triosephosphate isomerase), the biosynthesis of lipids (glycerol phosphate dehydrogenase), the hexose monophosphate shunt (6-phosphogluconate dehydrogenase), and the mevalonate pathway (isopentenyl diphosphate isomerase), catalyze the final step in the biosynthesis of pyrimidine nucleotides (orotidine monophosphate decarboxylase), and regulate the cellular levels of adenine nucleotides (adenylate kinase). The evolution of enzymes that undergo ligand-driven conformational changes to form active protein-substrate cages is proposed to proceed by selection of variants, in which the selected side chain substitutions destabilize a second protein conformer that shows compensating enhanced binding interactions with the substrate. The advantages inherent to enzymes that incorporate a conformational change into the catalytic cycle provide a strong driving force for the evolution of flexible protein folds such as the TIM barrel. The appearance of these folds represented a watershed event in enzyme evolution that enabled the rapid propagation of enzyme activities within enzyme superfamilies.

Covalent Allosteric Inhibitors of Akt Generated Using a Click Fragment Approach

Akt is a protein kinase that has been implicated in the progression of cancerous tumours. A number of covalent allosteric Akt inhibitors are known, and based on these scaffolds, a small library of novel potential covalent allosteric imidazopyridine-based inhibitors was designed. The envisaged compounds were synthesised, with click chemistry enabling a modular approach to a number of the target compounds. The binding modes, potencies and antiproliferative activities of these synthesised compounds were explored, thereby furthering the structure activity relationship knowledge of this class of Akt inhibitors. Three novel covalent inhibitors were identified, exhibiting moderate activity against Akt1 and various cancer cell lines, potentially paving the way for future covalent allosteric inhibitors with improved properties.

The quest for magic: recent advances in C(sp3)–H methylation

Frequently referred to as the “magic methyl” effect, the introduction of a methyl group into a biologically active molecule has the potential to drastically alter its physical and biological properties and significantly increase potency. This effect is most pronounced when the methyl group is added at the α-position of an aliphatic heterocycle or ortho to a large rotatable group on an aromatic ring. Although seminal developments in C–H activation strategies offered solutions to the latter, until recent years there had been no selective and functional-group-tolerant method for C(sp3)–H methylation at late stages of synthesis. For many years, the lack of a generally applicable methylation strategy necessitated arduous de novo synthesis approaches to access methylated drug candidates, and discouraged further investigation and understandings of the magic methyl effect. This review will provide a summary of the most recent advances that enabled non-directed late-stage C(sp3)–H methylation, including through hydride transfer, chemical or anodic oxidation, and photocatalytic hydrogen atom transfer.

Facilitating Rh-Catalyzed C-H Alkylation of (Hetero)arenes and 6-Arylpurine Nucleosides (Nucleotides) with Electrochemistry

An electrochemical approach to promote the ortho-C-H alkylation of (hetero)arenes via rhodium catalysis under mild conditions is described. This approach features mild conditions with high levels of regio- and monoselectivity that tolerate a variety of aromatic and heteroaromatic groups and offers a widely applicable method for late-stage diversification of complex molecular architectures including tryptophan, estrone, diazepam, nucleosides, and nucleotides. Alkyl boronic acids and esters and alkyl trifluoroborates are demonstrated as suitable coupling partners. The isolation of key rhodium intermediates and mechanistic studies provided strong support for a rhodium(III/IV or V) regime.

AutoDesigner, a De Novo Design Algorithm for Rapidly Exploring Large Chemical Space for Lead Optimization: Application to the Design and Synthesis of d-Amino Acid Oxidase Inhibitors

The lead optimization stage of a drug discovery program generally involves the design, synthesis, and assaying of hundreds to thousands of compounds. The design phase is usually carried out via traditional medicinal chemistry approaches and/or structure-based drug design (SBDD) when suitable structural information is available. Two of the major limitations of this approach are (1) difficulty in rapidly designing potent molecules that adhere to myriad project criteria, or the multiparameter optimization (MPO) problem, and (2) the relatively small number of molecules explored compared to the vast size of chemical space. To address these limitations, we have developed AutoDesigner, a de novo design algorithm. AutoDesigner employs a cloud-native, multistage search algorithm to carry out successive rounds of chemical space exploration and filtering. Millions to billions of virtual molecules are explored and optimized while adhering to a customizable set of project criteria such as physicochemical properties and potency. Additionally, the algorithm only requires a single ligand with measurable affinity and a putative binding model as a starting point, making it amenable to the early stages of an SBDD project where limited data are available. To assess the effectiveness of AutoDesigner, we applied it to the design of novel inhibitors of d-amino acid oxidase (DAO), a target for the treatment of schizophrenia. AutoDesigner was able to generate and efficiently explore over 1 billion molecules to successfully address a variety of project goals. The compounds generated by AutoDesigner that were synthesized and assayed (1) simultaneously met not only physicochemical criteria, clearance, and central nervous system (CNS) penetration (Kp,uu) cutoffs but also potency thresholds and (2) fully utilize structural data to discover and explore novel interactions and a previously unexplored subpocket in the DAO active site. The reported data demonstrate that AutoDesigner can play a key role in accelerating the discovery of novel, potent chemical matter within the constraints of a given drug discovery lead optimization campaign.

A Study on the Effect of the Substituent against PAK4 Inhibition Using In Silico Methods

The intrinsic inductive properties of atoms or functional groups depend on the chemical properties of either electron-withdrawing groups (EWGs) or electron-donating groups (EDGs). This study aimed to evaluate in silico methods to determine whether changes in chemical properties of the compound by single atomic substitution affect the biological activity of target proteins and whether the results depend on the properties of the functional groups. We found an imidazo[4,5-b]pyridine-based PAK4 inhibitor, compound 1, as an initial hit compound with the well-defined binding mode for PAK4. In this study, we used both experimental and in silico methods to investigate the effect of atomic substitution on biological activity to optimize the initial hit compound. In biological assays, in the case of EWG, as the size of the halogen atom became smaller and the electronegativity increased, the biological activity IC₅₀ value ranged from 5150 nM to inactive; in the case of EDG, biological activity was inactive. Furthermore, we analyzed the interactions of PAK4 with compounds, focusing on the hinge region residues, L398 and E399, and gatekeeper residues, M395 and K350, of the PAK4 protein using molecular docking studies and fragment molecular orbital (FMO) methods to determine the differences between the effect of EWG and EDG on the activity of target proteins. These results of the docking score and binding energy did not explain the differences in biological activity. However, the pair-interaction energy obtained from the results of the FMO method indicated that there was a difference in the interaction energy between the EWG and EDG in the hinge region residues, L398 and E399, as well as in M395 and K350. The two groups with different properties exhibited opposite electrostatic energy and charge transfer energy between L398 and E399. Additionally, we investigated the electron distribution of the parts interacting with the hinge region by visualizing the molecular electrostatic potential (MEP) surface of the compounds. In conclusion, we described the properties of functional groups that affect biological activity using an in silico method, FMO.

Facile synthesis of 5-Isopropyl-2,3-dimethylbenzene-1,4-diol by Friedel-Crafts and Determination of Partition Coefficient in n- Octanol/Water

The wide therapeutic effect of quinone-based drugs has received considerable interest for a long time. In this research, Friedel-Crafts performed a facile synthesis of quinone derivatives using the mixture of Brønsted acid. Reflux of 2,3-dimethylhydroquinone (1), isopropanol, glacial acetic acid, and H2SO4 for 15 minutes gave yellow oil product of 5-isopropyl-2,3-dimethylbenzene-1,4-diol (2) as a major product. Characterization using Nuclear Magnetic Resonance (NMR) revealed the methine proton splitting for isopropyl at δ 3.13 ppm, which has a cross-coupling with aromatic carbon at δ 119.6 ppm suggested the substitution of a proton on quinone ring with isopropyl group. Analysis Fourier Transform Infra-Red (FT-IR) showed the broad spectrum of –OH, the vibration of CH sp3, and isopropyl groups. The minor products identified as 5-isopropyl-2,3-dimethylcyclohexa-2,5-diene-1,4-dione (3), 5-isopropyl-2,3-dimethyl-1,4 phenylene diacetate (4), and 2,3-dimethylbenzene-5,6-isopropyl-1,4-diol (5) confirmed from 2D HETCOR and MS analysis. The partition coefficient (log P) of compound 2 showed a higher solubility by 1.9-fold compared to hydroquinone 1. It is suggested that an additional methyl group increased the partition into the organic phase.

Synthesis and Structure–Activity Relationships of Novel Non-Steroidal CYP17A1 Inhibitors as Potential Prostate Cancer Agents

Twenty new compounds, targeting CYP17A1, were synthesized, based on our previous work on a benzimidazole scaffold, and their biological activity evaluated. Inhibition of CYP17A1 is an important modality in the treatment of prostate cancer, which remains the most abundant cancer type in men. The biological assessment included CYP17A1 hydroxylase and lyase inhibition, CYP3A4 and P450 oxidoreductase (POR) inhibition, as well as antiproliferative activity in PC3 prostate cancer cells. The most potent compounds were selected for further analyses including in silico modeling. This combined effort resulted in a compound (comp 2, IC₅₀ 1.2 µM, in CYP17A1) with a potency comparable to abiraterone and selectivity towards the other targets tested. In addition, the data provided an understanding of the structure–activity relationship of this novel non-steroidal compound class.

Palladium-catalyzed nucleomethylation of alkynes for synthesis of methylated heteroaromatic compounds

Herein, we disclosed a novel and efficient palladium-catalyzed nucleomethylation of alkynes for the simultaneous construction of the heteroaromatic ring and methyl group. The 3-methylindoles, 3-methylbenzofurans and 4-methylisoquinolines were obtained in moderate to excellent yields. Notably, this methodology was employed as a key step for synthesis of a pregnane X receptor antagonist, zindoxifene, bazedoxifene and AFN-1252. The kinetic studies revealed that reductive elimination might be the rate-determining step.

A novel palladium-catalyzed nucleomethylation of alkynes is developed, affording 3-methylindoles, 3-methylbenzofurans and 4-methylisoquinolines in moderate to excellent yields.

Discovery of "Molecular Switches" within a Series of mGlu5 Allosteric Ligands Driven by a "Magic Methyl" Effect Affording Both PAMs and NAMs with In Vivo Activity, Derived from an M1 PAM Chemotype

In the course of optimizing an M₁ PAM chemotype, introduction of an ether moiety unexpectedly abolished M₁ PAM activity while engendering a "molecular switch" to afford a weak, pure mGlu₅ PAM. Further optimization was able to deliver a potent (mGlu₅ EC₅₀ = 520 nM, 63% Glu Max), centrally penetrant (K_p = 0.83), MPEP-site binding mGlu₅ PAM 17a (VU6036486) that reversed amphetamine-induced hyperlocomotion. A pronounced "magic methyl" effect was noted with a regioisomeric methyl congener, leading to a change in pharmacology to afford a potent (mGlu₅ IC₅₀ = 110 nM, 3% Glu Min), centrally penetrant (K_p = 0.94), MPEP-site binding NAM 28d (VU6044766) that displayed anxiolytic activity in a mouse marble burying assay. These data further support the growing body of literature concerning the existence of G protein-coupled receptor (GPCR) allosteric privileged structures, and the value and impact of subtle methyl group walks, as well as the highly productive fluorine walk, around allosteric ligand cores to stabilize unique GPCR conformations.

Systematic identification of activity cliffs with dual-atom replacements and their rationalization on the basis of single-atom replacement analogs and X-ray structures

Very small chemical changes in active compounds causing large potency effects are of particular interest in medicinal chemistry and drug design. We have systematically searched active compounds with available high-confidence activity data for pairs of structural analogs with dual-atom replacements and additional analogs with corresponding single-atom replacements. From ~287,000 unique qualifying compounds with activity against nearly 1900 unique targets, ~3500 target-based analog pairs with dual-atom replacements were identified. These included 852 pairs with significant differences in compound potency, representing a set of previously unobserved activity cliffs. Comparing these pairs with corresponding single-atom replacement analogs, which were frequently identified, made it possible to systematically analyze how potency changes propagated from single- to dual-atom replacements. The analysis uncovered different potency effects and revealed that individual atom replacements were often decisive for activity cliff formation. For a limited number of activity cliffs, X-ray structures of targets in complex with cliff compounds were available, which aided in rationalizing potency alterations among analogs with single- or dual-atom replacements. The analog pairs identified herein provide a rich resource of structure-activity relationship information and attractive test cases for calibrating computational methods.

Picomolar FKBP inhibitors enabled by a single water-displacing methyl group in bicyclic [4.3.1] aza-amides

Methyl groups can have profound effects in drug discovery but the underlying mechanisms are diverse and incompletely understood. Here we report the stereospecific effect of a single, solvent-exposed methyl group in bicyclic [4.3.1] aza-amides, robustly leading to a 2 to 10-fold increase in binding affinity for FK506-binding proteins (FKBPs). This resulted in the most potent and efficient FKBP ligands known to date. By a combination of co-crystal structures, isothermal titration calorimetry (ITC), density-functional theory (DFT), and 3D reference interaction site model (3D-RISM) calculations we elucidated the origin of the observed affinity boost, which was purely entropically driven and relied on the displacement of a water molecule at the protein-ligand-bulk solvent interface. The best compounds potently occupied FKBPs in cells and enhanced bone morphogenic protein (BMP) signaling. Our results show how subtle manipulation of the solvent network can be used to design atom-efficient ligands for difficult, solvent-exposed binding pockets.

Structural Optimization and Biological Activity of Pyrazole Derivatives: Virtual Computational Analysis, Recovery Assay and 3D Culture Model as Potential Predictive Tools of Effectiveness against Trypanosoma cruzi

Chagas disease, a chronic and silent disease caused by Trypanosoma cruzi, is currently a global public health problem. The treatment of this neglected disease relies on benznidazole and nifurtimox, two nitroheterocyclic drugs that show limited efficacy and severe side effects. The failure of potential drug candidates in Chagas disease clinical trials highlighted the urgent need to identify new effective chemical entities and more predictive tools to improve translational success in the drug development pipeline. In this study, we designed a small library of pyrazole derivatives (44 analogs) based on a hit compound, previously identified as a T. cruzi cysteine protease inhibitor. The in vitro phenotypic screening revealed compounds 3g, 3j, and 3m as promising candidates, with IC₅₀ values of 6.09 ± 0.52, 2.75 ± 0.62, and 3.58 ± 0.25 µM, respectively, against intracellular amastigotes. All pyrazole derivatives have good oral bioavailability prediction. The structure–activity relationship (SAR) analysis revealed increased potency of 1-aryl-1H-pyrazole-imidazoline derivatives with the Br, Cl, and methyl substituents in the para-position. The 3m compound stands out for its trypanocidal efficacy in 3D microtissue, which mimics tissue microarchitecture and physiology, and abolishment of parasite recrudescence in vitro. Our findings encourage the progression of the promising candidate for preclinical in vivo studies.

Engineering Fluorophore Recycling in a Fluorogenic RNA Aptamer

Fluorogenic aptamers can potentially show minimal photobleaching during continuous irradiation since any photobleached fluorophore can exchange with fluorescent dyes in the media. However, fluorophores have not been designed to maximize "fluorophore recycling." Here we describe TBI, a novel fluorophore for the Broccoli fluorogenic aptamer. Previous fluorophores either fail to rapidly dissociate when they undergo photobleaching via cis-trans isomerization, or bind slowly, resulting in extended periods after dissociation of the photobleached fluorophore when no fluorophore is bound. By contrast, photobleached TBI dissociates rapidly from Broccoli, and TBI from the media rapidly replaces dissociated photobleached fluorophore. Using TBI, Broccoli exhibits markedly enhanced fluorescence in cells during continuous imaging. These data show that designing fluorophores to optimize fluorophore recycling can lead to enhanced fluorescence of fluorogenic aptamers.

Homogeneous catalytic C(sp3)-H functionalization of gaseous alkanes

The conversion of light alkanes into bulk chemicals is becoming an important challenge as it effectively avoids the use of prefunctionalized alkylating reagents. The implementation of such processes is, however, hampered by their gaseous nature and low solubility, as well as the low reactivity of the C-H bonds. Efforts have been made to enable both polar and radical processes to activate these inert compounds. In addition, these methodologies also benefit significantly from the development of a suitable reactor technology that intensifies gas-liquid mass transfer. In this review, we critically highlight these developments, both from a conceptual and a practical point of view. The recent expansion of these mechanistically-different methods have enabled the use of various gaseous alkanes for the development of different bond-forming reactions, including C-C, C-B, C-N, C-Si and C-S bonds.

Influence of ring size in conformationally restricted ring I analogs of paromomycin on antiribosomal and antibacterial activity

In order to further investigate the importance of the conformation of the ring I side chain in aminoglycoside antibiotic binding to the ribosomal target several derivatives of paromomycin were designed with conformationally locked side chains. By changing the size of the appended ring between O-4' and C-6' used to restrict the motion of the side chain, the position of the C-6' hydroxy group was fine tuned to probe for the optimal conformation for inhibition of the ribosome. While the changes in orientation of the 6'-hydroxy group cannot be completely dissociated from the size and hydrophobicity of the conformation-restricting ring, overall, it is apparent that the preferred conformation of the ring I side chain for interaction with A1408 in the decoding A site of the bacterial ribosome is an ideal gt conformation, which results in the highest antimicrobial activity as well as increased selectivity for bacterial over eukaryotic ribosomes.

Deaminative Reductive Methylation of Alkylpyridinium Salts

Methyl groups can imbue valuable properties in organic molecules, often leading to enhanced bioactivity. To enable efficient installation of methyl groups on simple building blocks and in late-stage functionalization, a nickel-catalyzed reductive coupling of secondary Katritzky alkylpyridinium salts with methyl iodide was developed. When coupled with formation of the pyridinium salt from an alkyl amine, this method allows amino groups to be readily transformed to methyl groups with broad functional group and heterocycle tolerance.

The Deuterated "Magic Methyl" Group: A Guide to Site-Selective Trideuteromethyl Incorporation and Labeling by Using CD3 Reagents

In the field of medicinal chemistry, the precise installation of a trideuteromethyl group is gaining ever-increasing attention. Site-selective incorporation of the deuterated "magic methyl" group can provide profound pharmacological benefits and can be considered an important tool for drug optimization and development. This review provides a structured overview, according to trideuteromethylation reagent, of currently established methods for site-selective trideuteromethylation of carbon atoms. In addition to CD3 , the selective introduction of CD2 H and CDH2 groups is also considered. For all methods, the corresponding mechanism and scope are discussed whenever reported. As such, this review can be a starting point for synthetic chemists to further advance trideuteromethylation methodologies. At the same time, this review aims to be a guide for medicinal chemists, offering them the available C-CD3 formation strategies for the preparation of new or modified drugs.

Low rotational barriers for the most dynamically active methyl groups in the proposed antiviral drugs for treatment of SARS-CoV-2, apilimod and tetrandrine

A recent screening study highlighted a molecular compound, apilimod, for its efficacy against the SARS-CoV-2 virus, while another compound, tetrandrine, demonstrated a remarkable synergy with the benchmark antiviral drug, remdesivir. Here, we find that because of significantly reduced potential energy barriers, which also give rise to pronounced quantum effects, the rotational dynamics of the most dynamically active methyl groups in apilimod and tetrandrine are much faster than those in remdesivir. Because dynamics of methyl groups are essential for biochemical activity, screening studies based on the computed potential energy profiles may help identify promising candidates within a given class of drugs.

Sustainable manganese catalysis for late-stage C-H functionalization of bioactive structural motifs

The late-stage C–H functionalization of bioactive structural motifs is a powerful synthetic strategy for accessing advanced agrochemicals, bioimaging materials, and drug candidates, among other complex molecules. While traditional late-stage diversification relies on the use of precious transition metals, the utilization of 3d transition metals is an emerging approach in organic synthesis. Among the 3d metals, manganese catalysts have gained increasing attention for late-stage diversification due to the sustainability, cost-effectiveness, ease of operation, and reduced toxicity. Herein, we summarize recent manganese-catalyzed late-stage C–H functionalization reactions of biologically active small molecules and complex peptides.

N,N,N',N'-Tetramethylethylenediamine-Enabled Photoredox-Catalyzed C-H Methylation of N-Heteroarenes

Aiming at the valuable methylation process, readily available and inexpensive N,N,N',N'-tetramethylethylenediamine (TMEDA) was first identified as a new methyl source in photoredox-catalyzed transformation in this work. By virtue of this simple methylating reagent, a facile and practical protocol for the direct C-H methylation of N-heteroarenes was developed, featuring mild reaction conditions, broad substrate scope, and scalability. Mechanistic studies disclosed that a sequential photoredox, base-assisted proton shift, fragmentation, and tautomerization process was essentially involved.

Copper catalyzed late-stage C(sp3)-H functionalization of nitrogen heterocycles

Nitrogen heterocycle represents a ubiquitous skeleton in natural products and drugs. Late-stage C(sp3)-H bond functionalization of N-heterocycles with broad substrate scope remains a challenge and of particular significance to modern chemical synthesis and pharmaceutical chemistry. Here, we demonstrate copper-catalysed late-stage C(sp3)-H functionalizaion of N-heterocycles using commercially available catalysts under mild reaction conditions. We have investigated 8 types of N-heterocycles which are usually found as medicinally important skeletons. The scope and utility of this approach are demonstrated by late-stage C(sp3)-H modification of these heterocycles including a number of pharmaceuticals with a broad range of nucleophiles, e.g. methylation, arylation, azidination, mono-deuteration and glycoconjugation etc. Preliminary mechanistic studies reveal that the reaction undergoes a C-H fluorination process which is followed by a nucleophilic substitution.

Nonadditivity in public and inhouse data: implications for drug design

Numerous ligand-based drug discovery projects are based on structure-activity relationship (SAR) analysis, such as Free-Wilson (FW) or matched molecular pair (MMP) analysis. Intrinsically they assume linearity and additivity of substituent contributions. These techniques are challenged by nonadditivity (NA) in protein-ligand binding where the change of two functional groups in one molecule results in much higher or lower activity than expected from the respective single changes. Identifying nonlinear cases and possible underlying explanations is crucial for a drug design project since it might influence which lead to follow. By systematically analyzing all AstraZeneca (AZ) inhouse compound data and publicly available ChEMBL25 bioactivity data, we show significant NA events in almost every second assay among the inhouse and once in every third assay in public data sets. Furthermore, 9.4% of all compounds of the AZ database and 5.1% from public sources display significant additivity shifts indicating important SAR features or fundamental measurement errors. Using NA data in combination with machine learning showed that nonadditive data is challenging to predict and even the addition of nonadditive data into training did not result in an increase in predictivity. Overall, NA analysis should be applied on a regular basis in many areas of computational chemistry and can further improve rational drug design.

Synthesis and biological activity of cyclopropyl Δ7-dafachronic acids as DAF-12 receptor ligands

The four cyclopropyl stereoisomers of Δ7-dafachronic acids were prepared from the bile acid hyodeoxycholic acid and employed as chemical tools to exploit the importance of the orientation and spatial disposition of the carboxyl tail and the C25-methyl group for the binding at the DAF-12 receptor. The synthesis route was based on (a) Walden inversion and stereoselective PtO₂-hydrogenation to convert the L-shaped 5β-cholanoid scaffold into the planar 5α-sterol intermediate; (b) two-carbon homologation of the side chain by Wittig and cyclopropanation reaction; and (c) formation of the 3-keto group and Δ7 double bond. The synthesized isomers were isolated and tested for their activity as DAF-12 ligands by AlphaScreen assays. Results showed a significant loss of potency and efficacy for all the four stereoisomers when compared to the parent endogenous ligand. Computational analysis has evidenced the configurational and conformational arrangement of both the carboxylic and the C25-methyl group of dafachronic acids as key structural determinants for DAF-12 binding and activation.

Discovery of a Biased Allosteric Modulator for Cannabinoid 1 Receptor: Preclinical Anti-Glaucoma Efficacy

We apply the magic methyl effect to improve the potency/efficacy of GAT211, the prototypic 2-phenylindole-based cannabinoid type-1 receptor (CB1R) agonist-positive allosteric modulator (ago-PAM). Introducing a methyl group at the α-position of nitro group generated two diastereomers, the greater potency and efficacy of erythro, (±)-9 vs threo, (±)-10 constitutes the first demonstration of diastereoselective CB1R-allosteric modulator interaction. Of the (±)-9 enantiomers, (-)-(S,R)-13 evidenced improved potency over GAT211 as a CB1R ago-PAM, whereas (+)-(R,S)-14 was a CB1R allosteric agonist biased toward G protein- vs β-arrestin1/2-dependent signaling. (-)-(S,R)-13 and (+)-(R,S)-14 were devoid of undesirable side effects (triad test), and (+)-(R,S)-14 reduced intraocular pressure with an unprecedentedly long duration of action in a murine glaucoma model. (-)-(S,R)-13 docked into both a CB1R extracellular PAM and intracellular allosteric-agonist site(s), whereas (+)-(R,S)-14 preferentially engaged only the latter. Exploiting G-protein biased CB1R-allosteric modulation can offer safer therapeutic candidates for glaucoma and, potentially, other diseases.

Metal-, Photocatalyst-, and Light-Free Minisci C-H Acetylation of N-Heteroarenes with Vinyl Ethers

Herein, we report a mild, operationally simple method for Minisci C-H acetylation of N-heteroarenes using vinyl ethers as robust, inexpensive acetyl sources. The reactions do not require a conventional photocatalysis, electrocatalysis, metal catalysis, light activation, or high temperature. This method is thus significantly more sustainable than previously reported methods in terms of cost, reagent toxicity, and waste generation. This protocol can be expected to obtain medically relevant molecules from abundant feedstock materials.

Inhibition of CMP-sialic acid transport by endogenous 5-methyl CMP

Nucleotide-sugar transporters (NSTs) transport nucleotide-sugar conjugates into the Golgi lumen where they are then used in the synthesis of glycans. We previously reported crystal structures of a mammalian NST, the CMP-sialic acid transporter (CST) (Ahuja and Whorton 2019). These structures elucidated many aspects of substrate recognition, selectivity, and transport; however, one fundamental unaddressed question is how the transport activity of NSTs might be physiologically regulated as a means to produce the vast diversity of observed glycan structures. Here, we describe the discovery that an endogenous methylated form of cytidine monophosphate (m⁵CMP) binds and inhibits CST. The presence of m⁵CMP in cells results from the degradation of RNA that has had its cytosine bases post-transcriptionally methylated through epigenetic processes. Therefore, this work not only demonstrates that m⁵CMP represents a novel physiological regulator of CST, but it also establishes a link between epigenetic control of gene expression and regulation of glycosylation.

Discovery of the S1P2 Antagonist GLPG2938 (1-[2-Ethoxy-6-(trifluoromethyl)-4-pyridyl]-3-[[5-methyl-6-[1-methyl-3-(trifluoromethyl)pyrazol-4-yl]pyridazin-3-yl]methyl]urea), a Preclinical Candidate for the Treatment of Idiopathic Pulmonary Fibrosis

Mounting evidence from the literature suggests that blocking S1P2 receptor (S1PR2) signaling could be effective for the treatment of idiopathic pulmonary fibrosis (IPF). However, only a few antagonists have been so far disclosed. A chemical enablement strategy led to the discovery of a pyridine series with good antagonist activity. A pyridazine series with improved lipophilic efficiency and with no CYP inhibition liability was identified by scaffold hopping. Further optimization led to the discovery of 40 (GLPG2938), a compound with exquisite potency on a phenotypic IL8 release assay, good pharmacokinetics, and good activity in a bleomycin-induced model of pulmonary fibrosis.

Peroxide-mediated synthesis of benzimidazo[2,1-a]isoquinoline-6(5H)-ones via cascade methylation/ethylation and intramolecular cyclization

A metal-free oxidative radical methylation/arylation of 2-arylbenzoimidazoles with DTBP as the oxidant and methyl radical source was developed. The reaction proceeds through a sequential methyl radical addition/cyclization pathway and affords a series of methyl functionalized benzimidazo[2,1-a]isoquinoline-6(5H)-ones in moderate to good yields. Besides, the ethylation/arylation of 2-arylbenzoimidazoles was also achieved with DTAP.

Mechanochemical Solvent-Free Catalytic C-H Methylation

The mechanochemical, solvent-free, highly regioselective, rhodium-catalyzed C-H methylation of (hetero)arenes is reported. The reaction shows excellent functional-group compatibility and is demonstrated to work for the late-stage C-H methylation of biologically active compounds. The method requires no external heating and benefits from considerably shorter reaction times than previous solution-based C-H methylation protocols. Additionally, the mechanochemical approach is shown to enable the efficient synthesis of organometallic complexes that are difficult to generate conventionally.

Iodine-promoted ring-opening methylation of benzothiazoles with dimethyl sulfite

A halogen-bond promoted ring-opening methylation of benzothiazoles has been developed using dimethyl sulphite as a methylating reagent in the presence of a base. This approach represents a simple and efficient synthesis of N-methyl-N-(o-methylthio)phenyl amides, and features direct construction of both N-Me and S-Me bonds in a one-pot reaction through the decomposition of easily prepared benzothiazoles.

Selected Thoughts on Hydrophobicity in Drug Design

The fundamental aim of drug design in research and development is to invent molecules with selective affinity towards desired disease-associated targets. At the atomic loci of binding surfaces, systematic structural variations can define affinities between drug candidates and biomolecules, and thereby guide the optimization of safety, efficacy and pharmacologic properties. Hydrophobic interaction between biomolecules and drugs is integral to binding affinity and specificity. Examples of antiviral drug discovery are discussed.

Regioselective C(sp3)-H fluorination of ketones: from methyl to the monofluoromethyl group

Herein, we report a novel strategy to access CH2F-containing ketones through Pd-catalysed β-selective methyl C(sp3)-H fluorination. The reaction features high regioselectivity and a broad substrate scope, constituting a modular method for the late-stage transformation of the native methyl (CH3) into the monofluoromethyl (CH2F) group.