Theoretical Studies of the Acid-Base Equilibria in a Model Active Site of the Human 20S Proteasome

Thermodynamics-informed neural networks and extensive data sets: key factors to accurate blind predictions of apparent pKa values in the euroSAMPL challenge

Microscopic and macroscopic pKa values for 35 compounds selected by the organizers of euroSAMPL 1 challenge were blindly predicted with our thermodynamics-informed empirical S + pKa model (ranked submission 0x4cb7101f). Our results have received the first overall rank from the challenge organizers. We describe our methodology and discuss evaluation methods.

Dynamic Protonation States Underlie Carbene Formation in ThDP-Dependent Enzymes: A Theoretical Study

The activation mechanism of thiamine diphosphate (ThDP) in enzymes has long been the subject of intense research and controversial discussion. Particularly contentious is the formation of a carbene intermediate, the first one observed in an enzyme. For the formation of the carbene to take place, both intramolecular and intermolecular proton transfer pathways have been proposed. However, the physiologically relevant pH of ThDP-dependent enzymes around neutrality does not seem to be suitable for the formation of such reactive chemical species. Herein, we investigate the general mechanism of activation of the ThDP cofactor in human transketolase (TKT), by means of electronic structure methods. We show that in the case of the human TKT, the carbene species is accessible through a pKa shift induced by the electrostatics of a neighboring histidine residue (H110), whose protonation state change modulates the pKa of ThDP and suppresses the latter by more than 6 pH units. Our findings highlight that ThDP enzymes activate the cofactor beyond simple geometric constraints and the canonical glutamate. Such observations in nature can pave the way for the design of biomimetic carbene catalysts and the engineering of tailored enzymatic carbenes.

The Catalytic Mechanism of Acetoacetate Decarboxylase: A Detailed Study of Schiff Base Formation, Protonation States, and Their Impact on Catalysis

The enzyme acetoacetate decarboxylase (AAD) has a crucial function in the process of decarboxylating the substrate acetoacetate (AA). It has been extensively studied over the years, but its exact catalytic mechanism has remained partly unsolved due to the difficulty in assessing reaction intermediates. In this study, we combine molecular dynamics and electronic structure calculations to rediscover its catalytic mechanism. Our results show that the presence of the substrate, the acetoacetate, significantly influences the electrostatic potential of the active site. Furthermore, our simulations show that the decarboxylation reaction can take place by means of a direct proton transfer instead of via an enamine intermediate, which is thought to be strictly necessary. This work provides new insights into the role of the electrostatic interactions on the catalytic activity of AAD and for the first time connects it to the catalytic mechanism of other decarboxylases.

Uncertainty Quantification of Reactivity Scales

According to Mayr, polar organic synthesis can be rationalized by a simple empirical relationship linking bimolecular rate constants to as few as three reactivity parameters. Here, we propose an extension to Mayr's reactivity method that is rooted in uncertainty quantification and transforms the reactivity parameters into probability distributions. Through uncertainty propagation, these distributions can be transformed into uncertainty estimates for bimolecular rate constants. Chemists can exploit these virtual error bars to enhance synthesis planning and to decrease the ambiguity of conclusions drawn from experimental data. We demonstrate the above at the example of the reference data set released by Mayr and co-workers [J. Am. Chem. Soc. 2001, 123, 9500; J. Am. Chem. Soc. 2012, 134, 13902]. As by-product of the new approach, we obtain revised reactivity parameters for 36 π-nucleophiles and 32 benzhydrylium ions.

Mechanistic and thermodynamic characterization of oxathiazolones as potent and selective covalent immunoproteasome inhibitors

The ubiquitin–proteasome system is responsible for the degradation of proteins and plays a critical role in key cellular processes. While the constitutive proteasome (cPS) is expressed in all eukaryotic cells, the immunoproteasome (iPS) is primarily induced during disease processes, and its inhibition is beneficial in the treatment of cancer, autoimmune disorders and neurodegenerative diseases. Oxathiazolones were reported to selectively inhibit iPS over cPS, and the inhibitory activity of several oxathiazolones against iPS was experimentally determined. However, the detailed mechanism of the chemical reaction leading to irreversible iPS inhibition and the key selectivity drivers are unknown, and separate characterization of the noncovalent and covalent inhibition steps is not available for several compounds. Here, we investigate the chemical reaction between oxathiazolones and the Thr1 residue of iPS by quantum mechanics/molecular mechanics (QM/MM) simulations to establish a plausible reaction mechanism and to determine the rate-determining step of covalent complex formation. The modelled binding mode and reaction mechanism are in line with the selective inhibition of iPS versus cPS by oxathiazolones. The k_inact value of several ligands was estimated by constructing the potential of mean force of the rate-determining step by QM/MM simulations coupled with umbrella sampling. The equilibrium constant K_i of the noncovalent complex formation was evaluated by classical force field-based thermodynamic integration. The calculated K_i and k_inact values made it possible to analyse the contribution of the noncovalent and covalent steps to the overall inhibitory activity. Compounds with similar intrinsic reactivities exhibit varying selectivities for iPS versus cPS owing to subtle differences in the binding modes that slightly affect K_i, the noncovalent affinity, and importantly alter k_inact, the covalent reactivity of the bound compounds. A detailed understanding of the inhibitory mechanism of oxathiazolones is useful in designing iPS selective inhibitors with improved drug-like properties.