A deconstruction-reconstruction strategy for pyrimidine diversification.

Structure-Activity Relationship (SAR) studies are fundamental to drug and agrochemical development, yet only a few synthetic strategies apply to the nitrogen heteroaromatics frequently encountered in small molecule candidates.1-3 Here, we present an alternative approach where we convert pyrimidine-containing compounds various other nitrogen heteroaromatics. Transforming pyrmidines into their corresponding N-arylpyrimidinium salts enables cleavage into a three-carbon iminoenamine building block, used for various heterocycle-forming reactions. This deconstruction-reconstruction sequence diversifies the initial pyrimidine core and enables access to various heterocycles, such as azoles.4 In effect, this approach allows heterocycle formation on complex molecules, resulting in analogs that would be challenging to obtain by other methods. We anticipate this deconstruction-reconstruction strategy will extend to other heterocycle classes.

Predicting Lewis Acidity: Machine Learning the Fluoride Ion Affinity of p-Block-Atom-Based Molecules.

"How strong is this Lewis acid?" is a question researchers often approach by calculating its fluoride ion affinity (FIA) with quantum chemistry. Here, we present FIA49k, an extensive FIA dataset with 48,986 data points calculated at the RI-DSD-BLYP-D3(BJ)/def2-QZVPP//PBEh-3c level of theory, including 13 different p-block atoms as the fluoride accepting site. The FIA49k dataset was used to train FIA-GNN, two message-passing graph neural networks, which predict gas and solution phase FIA values of molecules excluded from training with a mean absolute error of 14 kJ mol-1 (r2=0.93) from the SMILES string of the Lewis acid as the only input. The level of accuracy is notable, given the wide energetic range of 750 kJ mol-1 spanned by FIA49k. The model's value was demonstrated with four case studies, including predictions for molecules extracted from the Cambridge Structural Database and by reproducing results from catalysis research available in the literature. Weaknesses of the model are evaluated and interpreted chemically. FIA-GNN and the FIA49k dataset can be reached via a free web app (www.grebgroup.de/fia-gnn).

Understanding the infection severity and epidemiological characteristics of mpox in the UK.

In May 2022, individuals infected with the monkeypox virus were detected in the UK without clear travel links to endemic areas. Understanding the clinical characteristics and infection severity of mpox is necessary for effective public health policy. The study period of this paper, from the 1st June 2022 to 30th September 2022, included 3,375 individuals that tested positive for the monkeypox virus. The posterior mean times from infection to hospital admission and length of hospital stay were 14.89 days (95% Credible Intervals (CrI): 13.60, 16.32) and 7.07 days (95% CrI: 6.07, 8.23), respectively. We estimated the modelled Infection Hospitalisation Risk to be 4.13% (95% CrI: 3.04, 5.02), compared to the overall sample Case Hospitalisation Risk (CHR) of 5.10% (95% CrI: 4.38, 5.86). The overall sample CHR was estimated to be 17.86% (95% CrI: 6.06, 33.11) for females and 4.99% (95% CrI: 4.27, 5.75) for males. A notable difference was observed between the CHRs that were estimated for each sex, which may be indicative of increased infection severity in females or a considerably lower infection ascertainment rate. It was estimated that 74.65% (95% CrI: 55.78, 86.85) of infections with the monkeypox virus in the UK were captured over the outbreak.

Designing solvent systems using self-evolving solubility databases and graph neural networks.

Designing solvent systems is key to achieving the facile synthesis and separation of desired products from chemical processes, so many machine learning models have been developed to predict solubilities. However, breakthroughs are needed to address deficiencies in the model's predictive accuracy and generalizability; this can be addressed by expanding and integrating experimental and computational solubility databases. To maximize predictive accuracy, these two databases should not be trained separately, and they should not be simply combined without reconciling the discrepancies from different magnitudes of errors and uncertainties. Here, we introduce self-evolving solubility databases and graph neural networks developed through semi-supervised self-training approaches. Solubilities from quantum-mechanical calculations are referred to during semi-supervised learning, but they are not directly added to the experimental database. Dataset augmentation is performed from 11 637 experimental solubilities to >900 000 data points in the integrated database, while correcting for the discrepancies between experiment and computation. Our model was successfully applied to study solvent selection in organic reactions and separation processes. The accuracy (mean absolute error around 0.2 kcal mol-1 for the test set) is quantitatively useful in exploring Linear Free Energy Relationships between reaction rates and solvation free energies for 11 organic reactions. Our model also accurately predicted the partition coefficients of lignin-derived monomers and drug-like molecules. While there is room for expanding solubility predictions to transition states, radicals, charged species, and organometallic complexes, this approach will be attractive to predictive chemistry areas where experimental, computational, and other heterogeneous data should be combined.

Catalytic Effects of Active Site Conformational Change in the Allosteric Activation of Imidazole Glycerol Phosphate Synthase.

Imidazole glycerol phosphate synthase (IGPS) is a class-I glutamine amidotransferase (GAT) that hydrolyzes glutamine. Ammonia is produced and transferred to a second active site, where it reacts with N1-(5'-phosphoribosyl)-formimino-5-aminoimidazole-4-carboxamide ribonucleotide (PrFAR) to form precursors to purine and histidine biosynthesis. Binding of PrFAR over 25 Šaway from the active site increases glutaminase efficiency by ∼4500-fold, primarily altering the glutamine turnover number. IGPS has been the focus of many studies on allosteric communication; however, atomic details for how the glutamine hydrolysis rate increases in the presence of PrFAR are lacking. We present a density functional theory study on 237-atom active site cluster models of IGPS based on crystallized structures representing the inactive and allosterically active conformations and investigate the multistep reaction leading to thioester formation and ammonia production. The proposed mechanism is supported by similar, well-studied enzyme mechanisms, and the corresponding energy profile is consistent with steady-state kinetic studies of PrFAR + IGPS. Additional active site models are constructed to examine the relationship between active site structural change and transition-state stabilization via energy decomposition schemes. The results reveal that the inactive IGPS conformation does not provide an adequately formed oxyanion hole structure and that repositioning of the oxyanion strand relative to the substrate is vital for a catalysis-competent oxyanion hole, with or without the hVal51 dihedral flip. These findings are valuable for future endeavors in modeling the IGPS allosteric mechanism by providing insight into the atomistic changes required for rate enhancement that can inform suitable reaction coordinates for subsequent investigations.

Forecasting influenza hospital admissions within English sub-regions using hierarchical generalised additive models.

BACKGROUND: Seasonal influenza places a substantial burden annually on healthcare services. Policies during the COVID-19 pandemic limited the transmission of seasonal influenza, making the timing and magnitude of a potential resurgence difficult to ascertain and its impact important to forecast. METHODS: We have developed a hierarchical generalised additive model (GAM) for the short-term forecasting of hospital admissions with a positive test for the influenza virus sub-regionally across England. The model incorporates a multi-level structure of spatio-temporal splines, weekly cycles in admissions, and spatial correlation. Using multiple performance metrics including interval score, coverage, bias, and median absolute error, the predictive performance is evaluated for the 2022-2023 seasonal wave. Performance is measured against autoregressive integrated moving average (ARIMA) and Prophet time series models. RESULTS: Across the epidemic phases the hierarchical GAM shows improved performance, at all geographic scales relative to the ARIMA and Prophet models. Temporally, the hierarchical GAM has overall an improved performance at 7 and 14 day time horizons. The performance of the GAM is most sensitive to the flexibility of the smoothing function that measures the national epidemic trend. CONCLUSIONS: This study introduces an approach to short-term forecasting of hospital admissions for the influenza virus using hierarchical, spatial, and temporal components. The methodology was designed for the real time forecasting of epidemics. This modelling framework was used across the 2022-2023 winter for healthcare operational planning by the UK Health Security Agency and the National Health Service in England.

Seasonal influenza causes a burden for hospitals and therefore it is useful to be able to accurately predict how many patients might be admitted with the disease. We attempted to predict influenza admissions up to 14 days in the future by creating a computational model that incorporates how the disease is reported and how it spreads. We evaluated our optimised model on data acquired during the winter of 2022-2023 data in England and compared it with previously developed models. Our model was better at modelling how influenza spreads and predicting future hospital admissions than the models we compared it to. Improving how influenza admissions are forecast can enable hospitals to prepare better for increased admissions, enabling improved treatment and reduced death for all patients in hospital over winter.

Electrochemical Modification of Polypeptides at Selenocysteine.

Mild strategies for the selective modification of peptides and proteins are in demand for applications in therapeutic peptide and protein discovery, and in the study of fundamental biomolecular processes. Herein, we describe the development of an electrochemical selenoetherification (e-SE) platform for the efficient site-selective functionalization of polypeptides. This methodology utilizes the unique reactivity of the 21st amino acid, selenocysteine, to effect formation of valuable bioconjugates through stable selenoether linkages under mild electrochemical conditions. The power of e-SE is highlighted through late-stage C-terminal modification of the FDA-approved cancer drug leuprolide and assembly of a library of anti-HER2 affibody conjugates bearing complex cargoes. Following assembly by e-SE, the utility of functionalized affibodies for in vitro imaging and targeting of HER2 positive breast and lung cancer cell lines is also demonstrated.

Exploring Cuneanes as Potential Benzene Isosteres and Energetic Materials: Scope and Mechanistic Investigations into Regioselective Rearrangements from Cubanes.

Cuneane is a strained hydrocarbon that can be accessed via metal-catalyzed isomerization of cubane. The carbon atoms of cuneane define a polyhedron of the C2v point group with six faces─two triangular, two quadrilateral, and two pentagonal. The rigidity, strain, and unique exit vectors of the cuneane skeleton make it a potential scaffold of interest for the synthesis of functional small molecules and materials. However, the limited previous synthetic efforts toward cuneanes have focused on monosubstituted or redundantly substituted systems such as permethylated, perfluorinated, and bis(hydroxymethylated) cuneanes. Such compounds, particularly rotationally symmetric redundantly substituted cuneanes, have limited potential as building blocks for the synthesis of complex molecules. Reliable, predictable, and selective syntheses of polysubstituted cuneanes bearing more complex substitution patterns would facilitate the study of this ring system in myriad applications. Herein, we report the regioselective, AgI-catalyzed isomerization of asymmetrically 1,4-disubstituted cubanes to cuneanes. In-depth DFT calculations provide a charge-controlled regioselectivity model, and direct dynamics simulations indicate that the nonclassical carbocation invoked is short-lived and dynamic effects augment the charge model.

Iridium-Catalyzed Asymmetric Difunctionalization of C-C σ-Bonds Enabled by Ring-Strained Boronate Complexes.

Enantioenriched organoboron intermediates are important building blocks in organic synthesis and drug discovery. Recently, transition metal-catalyzed enantioselective 1,2-metalate rearrangements of alkenylboronates have emerged as an attractive protocol to access these valuable reagents by installing two different carbon fragments across C═C π-bonds. Herein, we report the development of an iridium-catalyzed asymmetric allylation-induced 1,2-metalate rearrangement of bicyclo[1.1.0]butyl (BCB) boronate complexes enabled by strain release, which allows asymmetric difunctionalization of C-C σ-bonds, including dicarbonation and carboboration. This protocol provides a variety of enantioenriched three-dimensional 1,1,3-trisubstituted cyclobutane products bearing a boronic ester that can be readily derivatized. Notably, the reaction gives trans diastereoisomers that result from an anti-addition across the C-C σ-bond, which is in contrast to the syn-additions observed for reactions promoted by PdII-aryl complexes and other electrophiles in our previous works. The diastereoselectivity has been rationalized based on a combination of experimental data and density functional theory calculations, which suggest that the BCB boronate complexes are highly nucleophilic and react via early transition states with low activation barriers.

Regiodivergent Nucleophilic Fluorination under Hydrogen Bonding Catalysis: A Computational and Experimental Study.

The controlled programming of regiochemical outcomes in nucleophilic fluorination reactions with alkali metal fluoride is a problem yet to be solved. Herein, two synergistic approaches exploiting hydrogen bonding catalysis are presented. First, we demonstrate that modulating the charge density of fluoride with a hydrogen-bond donor urea catalyst directly influences the kinetic regioselectivity in the fluorination of dissymmetric aziridinium salts with aryl and ester substituents. Moreover, we report a urea-catalyzed formal dyotropic rearrangement, a thermodynamically controlled regiochemical editing process consisting of C-F bond scission followed by fluoride rebound. These findings offer a route to access enantioenriched fluoroamine regioisomers from a single chloroamine precursor, and more generally, new opportunities in regiodivergent asymmetric (bis)urea-based organocatalysis.

Metal-Free Arylation of Benzothiophenes at C4 by Activation as their Benzothiophene S-Oxides.

Benzothiophenes, activated by oxidation to the corresponding S-oxides, undergo C-H/C-H-type coupling with phenols to give C4 arylation products. While an electron-withdrawing group at C3 of the benzothiophene is important, the process operates without a directing group and a metal catalyst, thus rendering it compatible with sensitive functionalities-e.g. halides and formyl groups. Quantum chemical calculations suggest a formal stepwise mechanism involving heterolytic cleavage of an aryloxysulfur species to give a π-complex of the corresponding benzothiophene and a phenoxonium cation. Subsequent addition of the phenoxonium cation to the C4 position of the benzothiophene is favored over the addition to C3; Fukui functions predict that the major regioisomer is formed at the more electron-rich position between C3 and C4. Varied selective manipulation of the benzothiophene products showcase the synthetic utility of the metal-free arylation process.

Control of stereogenic oxygen in a helically chiral oxonium ion.

The control of tetrahedral carbon stereocentres remains a focus of modern synthetic chemistry and is enabled by their configurational stability. By contrast, trisubstituted nitrogen1, phosphorus2 and sulfur compounds3 undergo pyramidal inversion, a fundamental and well-recognized stereochemical phenomenon that is widely exploited4. However, the stereochemistry of oxonium ions-compounds bearing three substituents on a positively charged oxygen atom-is poorly developed and there are few applications of oxonium ions in synthesis beyond their existence as reactive intermediates5,6. There are no examples of configurationally stable oxonium ions in which the oxygen atom is the sole stereogenic centre, probably owing to the low barrier to oxygen pyramidal inversion7 and the perception that all oxonium ions are highly reactive. Here we describe the design, synthesis and characterization of a helically chiral triaryloxonium ion in which inversion of the oxygen lone pair is prevented through geometric restriction to enable it to function as a determinant of configuration. A combined synthesis and quantum calculation approach delineates design principles that enable configurationally stable and room-temperature isolable salts to be generated. We show that the barrier to inversion is greater than 110 kJ mol-1 and outline processes for resolution. This constitutes, to our knowledge, the only example of a chiral non-racemic and configurationally stable molecule in which the oxygen atom is the sole stereogenic centre.

Transmission dynamics of monkeypox in the United Kingdom: contact tracing study.

OBJECTIVE: To analyse the transmission dynamics of the monkeypox outbreak in the UK, declared a Public Health Emergency of International Concern in July 2022. DESIGN: Contact tracing study, linking data on case-contact pairs and on probable exposure dates. SETTING: Case questionnaires from the UK Health Security Agency (UKHSA), United Kingdom. PARTICIPANTS: 2746 people with polymerase chain reaction confirmed monkeypox virus in the UK between 6 May and 1 August 2022. MAIN OUTCOME MEASURES: The incubation period and serial interval of a monkeypox infection using two bayesian time delay models-one corrected for interval censoring (ICC-interval censoring corrected) and one corrected for interval censoring, right truncation, and epidemic phase bias (ICRTC-interval censoring right truncation corrected). Growth rates of cases by reporting date, when monkeypox virus was confirmed and reported to UKHSA, were estimated using generalised additive models. RESULTS: The mean age of participants was 37.8 years and 95% reported being gay, bisexual, and other men who have sex with men (1160 out of 1213 reporting). The mean incubation period was estimated to be 7.6 days (95% credible interval 6.5 to 9.9) using the ICC model and 7.8 days (6.6 to 9.2) using the ICRTC model. The estimated mean serial interval was 8.0 days (95% credible interval 6.5 to 9.8) using the ICC model and 9.5 days (7.4 to 12.3) using the ICRTC model. Although the mean serial interval was longer than the incubation period for both models, short serial intervals were more common than short incubation periods, with the 25th centile and the median of the serial interval shorter than the incubation period. For the ICC and ICRTC models, the corresponding estimates ranged from 1.8 days (95% credible interval 1.5 to 1.8) to 1.6 days (1.4 to 1.6) shorter at the 25th centile and 1.6 days (1.5 to 1.7) to 0.8 days (0.3 to 1.2) shorter at the median. 10 out of 13 linked patients had documented pre-symptomatic transmission. Doubling times of cases declined from 9.07 days (95% confidence interval 12.63 to 7.08) on the 6 May, when the first case of monkeypox was reported in the UK, to a halving time of 29 days (95% confidence interval 38.02 to 23.44) on 1 August. CONCLUSIONS: Analysis of the instantaneous growth rate of monkeypox incidence indicates that the epidemic peaked in the UK as of 9 July and then started to decline. Short serial intervals were more common than short incubation periods suggesting considerable pre-symptomatic transmission, which was validated through linked patient level records. For patients who could be linked through personally identifiable data, four days was the maximum time that transmission was detected before symptoms manifested. An isolation period of 16 to 23 days would be required to detect 95% of people with a potential infection. The 95th centile of the serial interval was between 23 and 41 days, suggesting long infectious periods.

Halogenation of the 3-position of pyridines through Zincke imine intermediates.

Pyridine halogenation reactions are crucial for obtaining the vast array of derivatives required for drug and agrochemical development. However, despite more than a century of synthetic endeavors, halogenation processes that selectively functionalize the carbon-hydrogen bond in the 3-position of a broad range of pyridine precursors remain largely elusive. We report a reaction sequence of pyridyl ring opening, halogenation, and ring closing whereby the acyclic Zincke imine intermediates undergo highly regioselective halogenation reactions under mild conditions. Experimental and computational mechanistic studies indicate that the nature of the halogen electrophile can modify the selectivity-determining step. Using this method, we produced a diverse set of 3-halopyridines and demonstrated late-stage halogenation of complex pharmaceuticals and agrochemicals.

Site-selective photocatalytic functionalization of peptides and proteins at selenocysteine.

The importance of modified peptides and proteins for applications in drug discovery, and for illuminating biological processes at the molecular level, is fueling a demand for efficient methods that facilitate the precise modification of these biomolecules. Herein, we describe the development of a photocatalytic method for the rapid and efficient dimerization and site-specific functionalization of peptide and protein diselenides. This methodology, dubbed the photocatalytic diselenide contraction, involves irradiation at 450 nm in the presence of an iridium photocatalyst and a phosphine and results in rapid and clean conversion of diselenides to reductively stable selenoethers. A mechanism for this photocatalytic transformation is proposed, which is supported by photoluminescence spectroscopy and density functional theory calculations. The utility of the photocatalytic diselenide contraction transformation is highlighted through the dimerization of selenopeptides, and by the generation of two families of protein conjugates via the site-selective modification of calmodulin containing the 21st amino acid selenocysteine, and the C-terminal modification of a ubiquitin diselenide.

Modeling Catalysis in Allosteric Enzymes: Capturing Conformational Consequences.

Greater understanding of enzymatic mechanisms aids the discovery of new targets for biologics, the development of biocatalytic transformations, and de novo enzyme design. Methods using quantum mechanical (QM) potentials, such as Density Functional Theory (DFT), have enabled complex multistep enzymatic mechanisms to be studied, often in quantitative detail. Nevertheless, the dynamic interconversion of enzyme conformations between active and inactive catalytic forms, involving length- and timescales inaccessible to QM treatments, presents a formidable challenge for the development of computational models for allosterically modulated enzymes. We present an overview of the key concepts underlying multistate models of enzyme catalysis, enzyme allostery, and the challenge that large-scale conformational changes pose for methods using QM, QM/MM, and MM potentials. Structural clustering is highlighted as a valuable approach to bridge molecular dynamics conformational sampling of MM potentials and quantum chemical cluster models of catalysis. Particularly relevant to this discussion is structural allostery, which serves as the exemplar of conformational consequences. Here, a well-characterized allosteric enzyme, Imidazole Glycerol Phosphate Synthase (IGPS), is used to showcase the importance of multiple conformations and guide a new direction for qualitative understanding and quantitative modeling in enzyme catalysis.

Umpolung Synthesis of Pyridyl Ethers by BiV -Mediated O-Arylation of Pyridones.

We report that O-selective arylation of 2- and 4-pyridones with arylboronic acids is affected by a modular, bismacycle-based system. The utility of this umpolung approach to pyridyl ethers, which is complementary to conventional methods based on SN Ar or cross-coupling, is demonstrated through the concise synthesis of Ki6783 and picolinafen, and the formal synthesis of cabozantib and golvatinib. Computational investigations reveal that arylation proceeds in a concerted fashion via a 5-membered transition state. The kinetically-controlled regioselectivity for O-arylation-which is reversed relative to previous BiV -mediated pyridone arylations-is attributed primarily to the geometric constraints imposed by the bismacyclic scaffold.

Reading and erasing of the phosphonium analogue of trimethyllysine by epigenetic proteins.

N ε-Methylation of lysine residues in histones plays an essential role in the regulation of eukaryotic transcription. The 'highest' methylation mark, N ε-trimethyllysine, is specifically recognised by N ε-trimethyllysine binding 'reader' domains, and undergoes demethylation, as catalysed by 2-oxoglutarate dependent JmjC oxygenases. We report studies on the recognition of the closest positively charged N ε-trimethyllysine analogue, i.e. its trimethylphosphonium derivative (KPme3), by N ε-trimethyllysine histone binding proteins and Nε-trimethyllysine demethylases. Calorimetric and computational studies with histone binding proteins reveal that H3KP4me3 binds more tightly than the natural H3K4me3 substrate, though the relative differences in binding affinity vary. Studies with JmjC demethylases show that some, but not all, of them can accept the phosphonium analogue of their natural substrates and that the methylation state selectivity can be changed by substitution of nitrogen for phosphorus. The combined results reveal that very subtle changes, e.g. substitution of nitrogen for phosphorus, can substantially affect interactions between ligand and reader domains / demethylases, knowledge that we hope will inspire the development of highly selective small molecules modulating their activity.

Mechanistic Aspects on [3+2] Cycloaddition (32CA) Reactions of Azides to Nitroolefins: A Computational and Kinetic Study.

[3+2] cycloadditions of nitroolefins have emerged as a selective and catalyst-free alternative for the synthesis of 1,2,3-triazoles from azides. We describe mechanistic studies into the cycloaddition/rearomatization reaction sequence. DFT calculations revealed a rate-limiting cycloaddition step proceeding via an asynchronous TS with high kinetic selectivity for the 1,5-triazole. Kinetic studies reveal a second-order rate law, and 13 C kinetic isotopic effects at natural abundance were measured with a significant normal effect at the conjugated olefinic centers of 1.0158 and 1.0216 at the α and ß-carbons of ß-nitrostyrene. Distortion/interaction-activation strain and energy decomposition analyses revealed that the major regioisomeric pathway benefits from an earlier and less-distorted TS, while intermolecular interaction terms dominate the preference for 1,5- over 1,4-cycloadducts. In addition, the major regioisomer also has more favorable electrostatic and dispersion terms. Additionally, while static DFT calculations suggest a concerted but highly asynchronous Ei-type HNO2 elimination mechanism, quasiclassical direct-dynamics calculations reveal the existence of a dynamic intermediate.