Ricardo A. Mata

The inhibition mechanism of human 20S proteasomes enables next-generation inhibitor design

Insights into proteasome inhibition Proteasomes are large protein complexes that degrade and remove proteins to maintain proper cellular physiology and growth. Proteasomes are a validated target for anticancer therapy, but drug design has been hampered by poor understanding of how inhibitors interact with the active site. Schrader et al. succeeded in crystallizing various proteasome-inhibitor complexes. They subsequently obtained crystal structures for the native human proteasome and eight different inhibitor complexes at resolutions between 1.9 and 2.1 Å. The inhibitors sampled include drugs that are approved or in trial for cancer treatment. Science, this issue p. 594 High-resolution structures of human 20S proteasomes reveal chemical principles for next-generation drug design. The proteasome is a validated target for anticancer therapy, and proteasome inhibition is employed in the clinic for the treatment of tumors and hematological malignancies. Here, we describe crystal structures of the native human 20S proteasome and its complexes with inhibitors, which either are drugs approved for cancer treatment or are in clinical trials. The structure of the native human 20S proteasome was determined at an unprecedented resolution of 1.8 angstroms. Additionally, six inhibitor-proteasome complex structures were elucidated at resolutions between 1.9 and 2.1 angstroms. Collectively, the high-resolution structures provide new insights into the catalytic mechanisms of inhibition and necessitate a revised description of the proteasome active site. Knowledge about inhibition mechanisms provides insights into peptide hydrolysis and can guide strategies for the development of next-generation proteasome-based cancer therapeutics.

Manganese(I)-Catalyzed Dispersion-Enabled C−H/C−C Activation

C-H/C-C Functionalizations were achieved with the aid of a versatile manganese(I) catalyst. Thus, an organometallic manganese-catalyzed C-H activation set the stage for silver-free C-H/C-C transformations with ample substrate scope and excellent levels of chemo-, site-, and diastereo-selectivities. The robust nature of the manganese(I) catalysis regime was reflected by the first C-H/C-C functionalization on amino acids under racemization-free reaction conditions. Detailed experimental and computational mechanistic studies provided strong evidence for a facile C-H activation and a rate-determining C-C cleavage, with considerable contribution from London dispersion interactions.

Benchmarking Quantum Chemical Methods: Are We Heading in the Right Direction?

Theoreticians and experimentalists should work together more closely to establish reliable rankings and benchmarks for quantum chemical methods. Comparison to carefully designed experimental benchmark data should be a priority. Guidelines to improve the situation for experiments and calculations are proposed.

Visualizing dispersion interactions through the use of local orbital spaces

The interpretation of chemical properties/phenomena can often be aided through the use of imagery. The mapping of molecular electrostatic potentials is a prime example, serving as a guideline in the design of novel compounds or understanding transition state stabilization effects. It is today a common tool for theoreticians and experimentalists alike. With the emergence of concepts such as dispersion energy donors, and the overall importance of dispersion in chemical systems, representations targeting such a class of interactions are warranted. In this work, we make use of local orbital analysis to extract dispersion interactions and represent them in a scalar quantity, the Dispersion Interaction Density (DID). A particular advantage of the method is the possibility to represent at the same footing intermolecular and intramolecular interactions in a straightforward fashion from wave function calculations. We present examples for the benzene dimer, several substituted benzenes and a coupled diamondoid molecule.

The furan microsolvation blind challenge for quantum chemical methods: First steps

Herein we present the results of a blind challenge to quantum chemical methods in the calculation of dimerization preferences in the low temperature gas phase. The target of study was the first step of the microsolvation of furan, 2-methylfuran and 2,5-dimethylfuran with methanol. The dimers were investigated through IR spectroscopy of a supersonic jet expansion. From the measured bands, it was possible to identify a persistent hydrogen bonding OH-O motif in the predominant species. From the presence of another band, which can be attributed to an OH-π interaction, we were able to assert that the energy gap between the two types of dimers should be less than or close to 1 kJ/mol across the series. These values served as a first evaluation ruler for the 12 entries featured in the challenge. A tentative stricter evaluation of the challenge results is also carried out, combining theoretical and experimental results in order to define a smaller error bar. The process was carried out in a double-blind fashion, with both theory and experimental groups unaware of the results on the other side, with the exception of the 2,5-dimethylfuran system which was featured in an earlier publication.

Efficient parallelization of perturbative Monte Carlo QM/MM simulations in heterogeneous platforms

A novel perturbative Monte Carlo mixed quantum mechanics (QM)/molecular mechanics (MM) approach has been recently developed to simulate molecular systems in complex environments. However, the required accuracy to efficiently simulate such complex molecular systems is usually granted at the cost of long executing times. To alleviate this problem, a new parallelization strategy of multi-level Monte Carlo molecular simulations is herein proposed for heterogeneous systems. It simultaneously exploits fine-grained (at the data level), coarse-grained (at the Markov chain level) and task-grained (pure QM, pure MM and QM/MM procedures) parallelism to ensure an efficient execution in heterogeneous systems composed of central processing units and multiple and possibly different graphical processing units. This is achieved by making use of the OpenCL library, together with appropriate dynamic load balancing schemes. From the conducted evaluation with real benchmarking data, a speed-up of 56x in the computational bottleneck part was observed, which results in a global speed-up of 38x for the whole simulation, reducing the time of a typical simulation from 80 hours to only 2 hours.

QM/MM study of the reaction mechanism of sulfite oxidase

Sulfite oxidase is a mononuclear molybdenum enzyme that oxidises sulfite to sulfate in many organisms, including man. Three different reaction mechanisms have been suggested, based on experimental and computational studies. Here, we study all three with combined quantum mechanical (QM) and molecular mechanical (QM/MM) methods, including calculations with large basis sets, very large QM regions (803 atoms) and QM/MM free-energy perturbations. Our results show that the enzyme is set up to follow a mechanism in which the sulfur atom of the sulfite substrate reacts directly with the equatorial oxo ligand of the Mo ion, forming a Mo-bound sulfate product, which dissociates in the second step. The first step is rate limiting, with a barrier of 39–49 kJ/mol. The low barrier is obtained by an intricate hydrogen-bond network around the substrate, which is preserved during the reaction. This network favours the deprotonated substrate and disfavours the other two reaction mechanisms. We have studied the reaction with both an oxidised and a reduced form of the molybdopterin ligand and quantum-refinement calculations indicate that it is in the normal reduced tetrahydro form in this protein.

The phenyl vinyl ether–methanol complex: a model system for quantum chemistry benchmarking

The structure of the isolated aggregate of phenyl vinyl ether and methanol is studied by combining a multi-spectroscopic approach and quantum-chemical calculations in order to investigate the delicate interplay of noncovalent interactions. The complementary results of vibrational and rotational spectroscopy applied in molecular beam experiments reveal the preference of a hydrogen bond of the methanol towards the ether oxygen (OH∙∙∙O) over the π-docking motifs via the phenyl and vinyl moieties, with an additional less populated OH∙∙∙P(phenyl)-bound isomer detected only by microwave spectroscopy. The correct prediction of the energetic order of the isomers using quantum-chemical calculations turns out to be challenging and succeeds with a sophisticated local coupled cluster method. The latter also yields a quantification as well as a visualization of London dispersion, which prove to be valuable tools for understanding the role of dispersion on the docking preferences. Beyond the structural analysis of the electronic ground state (S0), the electronically excited (S1) state is analyzed, in which a destabilization of the OH∙∙∙O structure compared to the S0 state is observed experimentally and theoretically.

Lipoxygenase 2 from Cyanothece sp. controls dioxygen insertion by steric shielding and substrate fixation

The biological function of lipoxygenases depends on the regio and stereo specific formation of fatty acid-derived hydroperoxides and different concepts exist to explain the mechanism that directs dioxygen to a specific carbon atom within the substrate. Here, we report the 1.8 Å resolution crystal structure of a cyanobacterial lipoxygenase that produces bis-allylic hydroperoxides (CspLOX2). Site directed mutagenesis experiments combined with computational approaches reveal that residues around the active site direct dioxygen to a preferred carbon atom and stereo configuration in the substrate fatty acid. Modulating the cavity volume around the pentadiene system of linoleic acid shifted the product formation towards 9S-, 9R-, 13S- or 13R-hydroperoxides in correlation with the site of mutation, thus decreasing the amount of the bis-allylic 11R-hydroperoxide. Decreasing the channel size of a 9R-lipoxygenase (CspLOX1) on the other hand could in turn induce formation of the bis-allylic 11R-hydroperoxide. Together this study suggests that an active site clamp fixing the pentadiene system of the substrate together with steric shielding controls the stereo and regio specific positioning of dioxygen at all positions of the reacting pentadiene system of substrate fatty acids.

Optimization and benchmarking of a perturbative Metropolis Monte Carlo quantum mechanics/molecular mechanics program

In this work, we present an optimized perturbative quantum mechanics/molecular mechanics (QM/MM) method for use in Metropolis Monte Carlo simulations. The model adopted is particularly tailored for the simulation of molecular systems in solution but can be readily extended to other applications, such as catalysis in enzymatic environments. The electrostatic coupling between the QM and MM systems is simplified by applying perturbation theory to estimate the energy changes caused by a movement in the MM system. This approximation, together with the effective use of GPU acceleration, leads to a negligible added computational cost for the sampling of the environment. Benchmark calculations are carried out to evaluate the impact of the approximations applied and the overall computational performance.