p38γ is essential for cell cycle progression and liver tumorigenesis.

The cell cycle is a tightly regulated process that is controlled by the conserved cyclin-dependent kinase (CDK)-cyclin protein complex1. However, control of the G0-to-G1 transition is not completely understood. Here we demonstrate that p38 MAPK gamma (p38γ) acts as a CDK-like kinase and thus cooperates with CDKs, regulating entry into the cell cycle. p38γ shares high sequence homology, inhibition sensitivity and substrate specificity with CDK family members. In mouse hepatocytes, p38γ induces proliferation after partial hepatectomy by promoting the phosphorylation of retinoblastoma tumour suppressor protein at known CDK target residues. Lack of p38γ or treatment with the p38γ inhibitor pirfenidone protects against the chemically induced formation of liver tumours. Furthermore, biopsies of human hepatocellular carcinoma show high expression of p38γ, suggesting that p38γ could be a therapeutic target in the treatment of this disease.

Decrypting Allostery in Membrane-Bound K-Ras4B Using Complementary In Silico Approaches Based on Unbiased Molecular Dynamics Simulations.

Protein functions are dynamically regulated by allostery, which enables conformational communication even between faraway residues, and expresses itself in many forms, akin to different "languages": allosteric control pathways predominating in an unperturbed protein are often unintuitively reshaped whenever biochemical perturbations arise (e.g., mutations). To accurately model allostery, unbiased molecular dynamics (MD) simulations require integration with a reliable method able to, e.g., detect incipient allosteric changes or likely perturbation pathways; this is because allostery can operate at longer time scales than those accessible by plain MD. Such methods are typically applied singularly, but we here argue their joint application─as a "multilingual" approach─could work significantly better. We successfully prove this through unbiased MD simulations (∼100 µs) of the widely studied, allosterically active oncotarget K-Ras4B, solvated and embedded in a phospholipid membrane, from which we decrypt allostery using four showcase "languages": Distance Fluctuation analysis and the Shortest Path Map capture allosteric hotspots at equilibrium; Anisotropic Thermal Diffusion and Dynamical Non-Equilibrium MD simulations assess perturbations upon, respectively, either superheating or hydrolyzing the GTP that oncogenically activates K-Ras4B. Chosen "languages" work synergistically, providing an articulate, mutually coherent, experimentally consistent picture of K-Ras4B allostery, whereby distinct traits emerge at equilibrium and upon GTP cleavage. At equilibrium, combined evidence confirms prominent allosteric communication from the membrane-embedded hypervariable region, through a hub comprising helix α5 and sheet ß5, and up to the active site, encompassing allosteric "switches" I and II (marginally), and two proposed pockets. Upon GTP cleavage, allosteric perturbations mostly accumulate on the switches and documented interfaces.

Harnessing conformational dynamics in enzyme catalysis to achieve nature-like catalytic efficiencies: the shortest path map tool for computational enzyme redesign.

Enzymes exhibit diverse conformations, as represented in the free energy landscape (FEL). Such conformational diversity provides enzymes with the ability to evolve towards novel functions. The challenge lies in identifying mutations that enhance specific conformational changes, especially if located in distal sites from the active site cavity. The shortest path map (SPM) method, which we developed to address this challenge, constructs a graph based on the distances and correlated motions of residues observed in nanosecond timescale molecular dynamics (MD) simulations. We recently introduced a template based AlphaFold2 (tAF2) approach coupled with 10 nanosecond MD simulations to quickly estimate the conformational landscape of enzymes and assess how the FEL is shifted after mutation. In this study, we evaluate the potential of SPM when coupled with tAF2-MD in estimating conformational heterogeneity and identifying key conformationally-relevant positions. The selected model system is the beta subunit of tryptophan synthase (TrpB). We compare how the SPM pathways differ when integrating tAF2 with different MD simulation lengths from as short as 10 ns until 50 ns and considering two distinct Amber forcefield and water models (ff14SB/TIP3P versus ff19SB/OPC). The new methodology can more effectively capture the distal mutations found in laboratory evolution, thus showcasing the efficacy of tAF2-MD-SPM in rapidly estimating enzyme dynamics and identifying the key conformationally relevant hotspots for computational enzyme engineering.

Controlling Monoterpene Isomerization by Guiding Challenging Carbocation Rearrangement Reactions in Engineered Squalene-Hopene Cyclases.

The interconversion of monoterpenes is facilitated by a complex network of carbocation rearrangement pathways. Controlling these isomerization pathways is challenging when using common Brønsted and Lewis acid catalysts, which often produce product mixtures that are difficult to separate. In contrast, natural monoterpene cyclases exhibit high control over the carbocation rearrangement reactions but are reliant on phosphorylated substrates. In this study, we present engineered squalene-hopene cyclases from Alicyclobacillus acidocaldarius (AacSHC) that catalyze the challenging isomerization of monoterpenes with unprecedented precision. Starting from a promiscuous isomerization of (+)-ß-pinene, we first demonstrate noticeable shifts in the product distribution solely by introducing single point mutations. Furthermore, we showcase the tuneable cation steering by enhancing (+)-borneol selectivity from 1 % to >90 % (>99 % de) aided by iterative saturation mutagenesis. Our combined experimental and computational data suggest that the reorganization of key aromatic residues leads to the restructuring of the water network that facilitates the selective termination of the secondary isobornyl cation. This work expands our mechanistic understanding of carbocation rearrangements and sets the stage for target-oriented skeletal reorganization of broadly abundant terpenes.

C-H Bonds as Functional Groups: Simultaneous Generation of Multiple Stereocenters by Enantioselective Hydroxylation at Unactivated Tertiary C-H Bonds.

Enantioselective C-H oxidation is a standing chemical challenge foreseen as a powerful tool to transform readily available organic molecules into precious oxygenated building blocks. Here, we describe a catalytic enantioselective hydroxylation of tertiary C-H bonds in cyclohexane scaffolds with H2O2, an evolved manganese catalyst that provides structural complementary to the substrate similarly to the lock-and-key recognition operating in enzymatic active sites. Theoretical calculations unveil that enantioselectivity is governed by the precise fitting of the substrate scaffold into the catalytic site, through a network of complementary weak non-covalent interactions. Stereoretentive C(sp3)-H hydroxylation results in a single-step generation of multiple stereogenic centers (up to 4) that can be orthogonally manipulated by conventional methods providing rapid access, from a single precursor to a variety of chiral scaffolds.

Brave new surfactant world revisited by thermoalkalophilic lipases: computational insights into the role of SDS as a substrate analog.

Non-ionic surfactants were shown to stabilize the active conformation of thermoalkalophilic lipases by mimicking the lipid substrate while the catalytic interactions formed by anionic surfactants have not been well characterized. In this study, we combined µs-scale molecular dynamics (MD) simulations and lipase activity assays to analyze the effect of ionic surfactant, sodium dodecyl sulfate (SDS), on the structure and activity of thermoalkalophilic lipases. Both the open and closed lipase conformations that differ in geometry were recruited to the MD analysis to provide a broader understanding of the molecular effect of SDS on the lipase structure. Simulations at 298 K showed the potential of SDS for maintaining the active lipase through binding to the sn-1 acyl-chain binding pocket in the open conformation or transforming the closed conformation to an open-like state. Consistent with MD findings, experimental analysis showed increased lipase activity upon SDS incubation at ambient temperature. Notably, the lipase cores stayed intact throughout 2 µs regardless of an increase in the simulation temperature or SDS concentration. However, the surface structures were unfolded in the presence of SDS and at elevated temperature for both conformations. Simulations of the dimeric lipase were also carried out and showed reduced flexibility of the surface structures which were unfolded in the monomer, indicating the insulating role of dimer interactions against SDS. Taken together, this study provides insights into the possible substrate mimicry by the ionic surfactant SDS for the thermoalkalophilic lipases without temperature elevation, underscoring SDS's potential for interfacial activation at ambient temperatures.

Harnessing the Structure and Dynamics of the Squalene-Hopene Cyclase for (-)-Ambroxide Production.

Terpene cyclases offer enormous synthetic potential, given their unique ability to forge complex hydrocarbon scaffolds from achiral precursors within a single cationic rearrangement cascade. Harnessing their synthetic power, however, has proved to be challenging owing to their generally low catalytic performance. In this study, we unveiled the catalytic potential of the squalene-hopene cyclase (SHC) by harnessing its structure and dynamics. First, we synergistically tailored the active site and entrance tunnel of the enzyme to generate a 397-fold improved (-)-ambroxide synthase. Our computational investigations explain how the introduced mutations work in concert to improve substrate acquisition, flow, and chaperoning. Kinetics, however, showed terpene-induced inactivation of the membrane-bound SHC to be the major turnover limitation in vivo. Merging this insight with the improved and stereoselective catalysis of the enzyme, we applied a feeding strategy to exceed 105 total turnovers. We believe that our results may bridge the gap for broader application of SHCs in synthetic chemistry.

Time Evolution of the Millisecond Allosteric Activation of Imidazole Glycerol Phosphate Synthase.

Deciphering the molecular mechanisms of enzymatic allosteric regulation requires the structural characterization of functional states and also their time evolution toward the formation of the allosterically activated ternary complex. The transient nature and usually slow millisecond time scale interconversion between these functional states hamper their experimental and computational characterization. Here, we combine extensive molecular dynamics simulations, enhanced sampling techniques, and dynamical networks to describe the allosteric activation of imidazole glycerol phosphate synthase (IGPS) from the substrate-free form to the active ternary complex. IGPS is a heterodimeric bienzyme complex whose HisH subunit is responsible for hydrolyzing glutamine and delivering ammonia for the cyclase activity in HisF. Despite significant advances in understanding the underlying allosteric mechanism, essential molecular details of the long-range millisecond allosteric activation of IGPS remain hidden. Without using a priori information of the active state, our simulations uncover how IGPS, with the allosteric effector bound in HisF, spontaneously captures glutamine in a catalytically inactive HisH conformation, subsequently attains a closed HisF:HisH interface, and finally forms the oxyanion hole in HisH for efficient glutamine hydrolysis. We show that the combined effector and substrate binding dramatically decreases the conformational barrier associated with oxyanion hole formation, in line with the experimentally observed 4500-fold activity increase in glutamine hydrolysis. The allosteric activation is controlled by correlated time-evolving dynamic networks connecting the effector and substrate binding sites. This computational strategy tailored to describe millisecond events can be used to rationalize the effect of mutations on the allosteric regulation and guide IGPS engineering efforts.

SBMOpenMM: A Builder of Structure-Based Models for OpenMM.

Molecular dynamics (MD) simulations have become a standard tool to correlate the structure and function of biomolecules and significant advances have been made in the study of proteins and their complexes. A major drawback of conventional MD simulations is the difficulty and cost of obtaining converged results, especially when exploring potential energy surfaces containing considerable energy barriers. This limits the wide use of MD calculations to determine the thermodynamic properties of biomolecular processes. Indeed, this is true when considering the conformational entropy of such processes, which is ultimately critical in assessing the simulations' convergence. Alternatively, a wide range of structure-based models (SBMs) has been used in the literature to unravel the basic mechanisms of biomolecular dynamics. These models introduce simplifications that focus on the relevant aspects of the physical process under study. Because of this, SBMs incorporate the need to modify the force field definition and parameters to target specific biophysical simulations. Here we introduce SBMOpenMM, a Python library to build force fields for SBMs, that uses the OpenMM framework to create and run SBM simulations. The code is flexible and user-friendly and profits from the high customizability and performance provided by the OpenMM platform.

Regioselective Synthesis and Characterization of Tris- and Tetra-Prato Adducts of M3N@C80 (M = Y, Gd).

The tris- and tetra-adducts of M3N@Ih-C80 metallofullerenes were synthesized and characterized for the first time. The 1,3-dipolar cycloaddition (Prato reaction) of Y3N@Ih-C80 and Gd3N@Ih-C80 with an excess of N-ethylglycine and formaldehyde provided tris- and tetra-fulleropyrrolidine adducts in a regioselective manner. Purification by HPLC and analyses of the isolated peaks by NMR, MS, and vis-NIR spectra revealed that the major products were four tris- and one tetra-isomers for both Y3N@Ih-C80 and Gd3N@Ih-C80. Considering the large number of possible isomers (e.g., at least 1140 isomers for the tris-adduct), the limited number of isomers obtained indicated that the reactions proceeded with high regioselectivity. NMR analyses of the Y3N@Ih-C80 adducts found that the tris-adducts were all-[6,6]- or [6,6][6,6][5,6]-isomers and that some showed mutual isomerization or remained intact at room temperature. The tetra-adduct obtained as a major product was all-[6,6] and stable. For the structural elucidation of Gd3N@Ih-C80 tris- and tetra-adducts, density functional theory (DFT) calculations were performed to estimate the relative stabilities of tris- and tetra-adducts formed upon Prato functionalization of the most pyramidalized regions of the fullerene structure. The most stable structures corresponded to additions on the most pyramidalized (i.e., strained) bonds. Taking together the experimental vis-NIR spectra, NMR assignments, and the computed relative DFT stabilities of the potential tris- and tetra-adducts, the structures of the isolated adducts were elucidated. Electron resonance (ESR) measurements measurements of pristine, bis-, and tris-adducts of Gd3N@C80 suggested that the rotation of the endohedral metal cluster slowed upon increase of the addition numbers to C80 cage, which is favored for accommodating the Gd atoms of the relatively large Gd3N cluster inner space at the sp3 addition sites. This is presumably related to the high regioselectivity in the Prato addition reaction driven by the strain release of the Gd3N@C80 fullerene structure.

Computational NMR Spectra of o-Benzyne and Stable Guests and Their Hemicarceplexes.

The incarceration of o-benzyne and 27 other guest molecules within hemicarcerand 1, as reported experimentally by Warmuth, and Cram and co-workers, has been studied by density functional theory (DFT). The 1 H NMR chemical shifts, rotational mobility, and conformational preference of the guests within the supramolecular cage were determined, which showed intriguing correlations of the chemical shifts with structural parameters of the host-guest system. Furthermore, based on the computed chemical shifts reassignments of some NMR signals are proposed. This affects, in particular, the putative characterization of the volatile benzyne molecule inside a hemicarcerand, for which our CCSD(T) and KT2 results indicate that the experimentally observed signals are most likely not resulting from an isolated o-benzyne within the supramolecular host. Instead, it is shown that the guest reacted with an aromatic ring of the host, and this adduct is responsible for the experimentally observed signals.

Enzyme Conformation Influences the Performance of Lipase-powered Nanomotors.

Enzyme-powered micro/nanomotors have myriads of potential applications in various areas. To efficiently reach those applications, it is necessary and critical to understand the fundamental aspects affecting the motion dynamics. Herein, we explored the impact of enzyme orientation on the performance of lipase-powered nanomotors by tuning the lipase immobilization strategies. The influence of the lipase orientation and lid conformation on substrate binding and catalysis was analyzed using molecular dynamics simulations. Besides, the motion performance indicates that the hydrophobic binding (via OTES) represents the best orienting strategy, providing 48.4 % and 95.4 % increase in diffusion coefficient compared to hydrophilic binding (via APTES) and Brownian motion (no fuel), respectively (with C[triacetin] of 100 mm). This work provides vital evidence for the importance of immobilization strategy and corresponding enzyme orientation for the catalytic activity and in turn, the motion performance of nanomotors, and is thus helpful to future applications.

Regio- and Stereoselective Steroid Hydroxylation at C7 by Cytochrome†P450 Monooxygenase Mutants.

Steroidal C7ß alcohols and their respective esters have shown significant promise as neuroprotective and anti-inflammatory agents to treat chronic neuronal damage like stroke, brain trauma, and cerebral ischemia. Since C7 is spatially far away from any functional groups that could direct C-H activation, these transformations are not readily accessible using modern synthetic organic techniques. Reported here are P450-BM3 mutants that catalyze the oxidative hydroxylation of six different steroids with pronounced C7 regioselectivities and ß stereoselectivities, as well as high activities. These challenging transformations were achieved by a focused mutagenesis strategy and application of a novel technology for protein library construction based on DNA assembly and USER (Uracil-Specific Excision Reagent) cloning. Upscaling reactions enabled the purification of the respective steroidal alcohols in moderate to excellent yields. The high-resolution X-ray structure and molecular dynamics simulations of the best mutant unveil the origin of regio- and stereoselectivity.

Deciphering the Allosterically Driven Conformational Ensemble in Tryptophan Synthase Evolution.

Multimeric enzyme complexes are ubiquitous in nature and catalyze a broad range of useful biological transformations. They are often characterized by a tight allosteric coupling between subunits, making them highly inefficient when isolated. A good example is Tryptophan synthase (TrpS), an allosteric heterodimeric enzyme in the form of an αßßα complex that catalyzes the biosynthesis of L-tryptophan. In this study, we decipher the allosteric regulation existing in TrpS from Pyrococcus furiosus (PfTrpS), and how the allosteric conformational ensemble is recovered in laboratory-evolved stand-alone ß-subunit variants. We find that recovering the conformational ensemble of a subdomain of TrpS affecting the relative stabilities of open, partially closed, and closed conformations is a prerequisite for enhancing the catalytic efficiency of the ß-subunit in the absence of its binding partner. The distal mutations resuscitate the allosterically driven conformational regulation and alter the populations and rates of exchange between these multiple conformational states, which are essential for the multistep reaction pathway of the enzyme. Interestingly, these distal mutations can be a priori predicted by careful analysis of the conformational ensemble of the TrpS enzyme through computational methods. Our study provides the enzyme design field with a rational approach for evolving allosteric enzymes toward improved stand-alone function for biosynthetic applications.

Structures of Gd3N@C80 Prato Bis-Adducts: Crystal Structure, Thermal Isomerization, and Computational Study.

The structures of two bis-ethylpyrrolidinoadducts of Gd3N@Ih-C80, obtained by regioselective 1,3-dipolar cycloadditions, were elucidated by single crystal X-ray, visible-near infrared (vis-NIR) spectra, studies on their thermal isomerization, and theoretical calculations. The structure of the minor-bis-adduct reveals a C2-symmetric carbon cage with [6,6][6,6]-addition sites and with an endohedral Gd3N cluster that is completely flattened. This is the first example of a crystal structure of Gd3N@Ih-C80 derivatives. The structure of the major-bis-adduct was inferred by the vis-NIR spectrum being corresponded to the structure of a previously reported major-bis-adduct of Y3N@Ih-C80 known to have an asymmetric [6,6][6,6]-structure. Based on experimental results showing that the minor-bis-adduct of Gd3N@Ih-C80 isomerized to the major-adduct, a possible second addition site was elucidated with support from density functional theory calculations.

Hidden Conformations in Aspergillus niger Monoamine Oxidase are Key for Catalytic Efficiency.

Enzymes exist as an ensemble of conformational states, whose populations can be shifted by substrate binding, allosteric interactions, but also by introducing mutations to their sequence. Tuning the populations of the enzyme conformational states through mutation enables evolution towards novel activity. Herein, Markov state models are used to unveil hidden conformational states of monoamine oxidase from Aspergillus niger (MAO-N). These hidden conformations, not previously observed by any other technique, play a crucial role in substrate binding and enzyme activity. This reveals how distal mutations regulate MAO-N activity by stabilizing these hidden, catalytically important conformational states, but also by modulating the communication pathway between both MAO-N subunits.

Epoxide Hydrolase Conformational Heterogeneity for the Resolution of Bulky Pharmacologically Relevant Epoxide Substrates.

The conformational landscape of Bacillus megaterium epoxide hydrolase (BmEH) and how it is altered by mutations that confer the enzyme the ability to accept bulky epoxide substrates has been investigated. Extensive molecular dynamics (MD) simulations coupled to active site volume calculations have unveiled relevant features of the enzyme conformational dynamics and function. Our long-timescale MD simulations identify key conformational states not previously observed by means of X-ray crystallography and short MD simulations that present the loop containing one of the catalytic residues, Asp239, in a wide-open conformation, which is likely involved in the binding of the epoxide substrate. Introduction of mutations M145S and F128A dramatically alters the conformational landscape of the enzyme. These singly mutated variants can accept bulky epoxide substrates due to the disorder induced by mutation in the α-helix containing the catalytic Tyr144 and some parts of the lid domain. These changes impact the enzyme active site, which is substantially wider and more complementary to the bulky pharmacologically relevant epoxide substrates.

Origins of stereoselectivity in evolved ketoreductases.

Mutants of Lactobacillus kefir short-chain alcohol dehydrogenase, used here as ketoreductases (KREDs), enantioselectively reduce the pharmaceutically relevant substrates 3-thiacyclopentanone and 3-oxacyclopentanone. These substrates differ by only the heteroatom (S or O) in the ring, but the KRED mutants reduce them with different enantioselectivities. Kinetic studies show that these enzymes are more efficient with 3-thiacyclopentanone than with 3-oxacyclopentanone. X-ray crystal structures of apo- and NADP(+)-bound selected mutants show that the substrate-binding loop conformational preferences are modified by these mutations. Quantum mechanical calculations and molecular dynamics (MD) simulations are used to investigate the mechanism of reduction by the enzyme. We have developed an MD-based method for studying the diastereomeric transition state complexes and rationalize different enantiomeric ratios. This method, which probes the stability of the catalytic arrangement within the theozyme, shows a correlation between the relative fractions of catalytically competent poses for the enantiomeric reductions and the experimental enantiomeric ratio. Some mutations, such as A94F and Y190F, induce conformational changes in the active site that enlarge the small binding pocket, facilitating accommodation of the larger S atom in this region and enhancing S-selectivity with 3-thiacyclopentanone. In contrast, in the E145S mutant and the final variant evolved for large-scale production of the intermediate for the antibiotic sulopenem, R-selectivity is promoted by shrinking the small binding pocket, thereby destabilizing the pro-S orientation.

Exploring the origins of selectivity in soluble epoxide hydrolase from Bacillus megaterium.

Epoxide hydrolase (EH) enzymes catalyze the hydration of racemic epoxides to yield their corresponding vicinal diols. These enzymes present different enantio- and regioselectivity depending upon either the substrate structure or the substitution pattern of the epoxide ring. In this study, we computationally investigate the Bacillus megaterium epoxide hydrolase (BmEH)-mediated hydrolysis of racemic styrene oxide (rac-SO) and its para-nitro styrene oxide (rac-p-NSO) derivative using density functional theory (DFT) and an active site cluster model consisting of 195 and 197 atoms, respectively. Full reaction mechanisms for epoxide ring opening were evaluated considering the attack at both oxirane carbons and considering two possible orientations of the substrate at the BmEH active site. Our results indicate that for both SO and p-NSO substrates the BmEH enantio- and regioselectivity is opposite to the inherent (R)-BmEH selectivity, the attack at the benzylic position (C1) of the (S)-enantiomer being the most favoured chemical outcome.

Exploring the reversal of enantioselectivity on a zinc-dependent alcohol dehydrogenase.

Alcohol Dehydrogenase (ADH) enzymes catalyse the reversible reduction of prochiral ketones to the corresponding alcohols. These enzymes present two differently shaped active site pockets, which dictate their substrate scope and selectivity. In this study, we computationally evaluate the effect of two commonly reported active site mutations (I86A, and W110T) on a secondary alcohol dehydrogenase from Thermoanaerobacter brockii (TbSADH) through Molecular Dynamics simulations. Our results indicate that the introduced mutations induce dramatic changes in the shape of the active site, but most importantly they impact the substrate-enzyme interactions. We demonstrate that the combination of Molecular Dynamics simulations with the tools POVME and NCIplot corresponds to a powerful strategy for rationalising and engineering the stereoselectivity of ADH variants.