DszA Catalyzes C-S Bond Cleavage through N5-Hydroperoxyl Formation.

Due to its detrimental impact on human health and the environment, regulations demand ultralow sulfur levels on fossil fuels, in particular in diesel. However, current desulfurization techniques are expensive and cannot efficiently remove heteroaromatic sulfur compounds, which are abundant in crude oil and concentrate in the diesel fraction after distillation. Biodesulfurization via the four enzymes of the metabolic 4S pathway of the bacterium Rhodococcus erythropolis (DszA-D) is a possible solution. However, the 4S pathway needs to operate at least 500 times faster for industrial applicability, a goal currently pursued through enzyme engineering. In this work, we unveil the catalytic mechanism of the flavin monooxygenase DszA. Surprisingly, we found that this enzyme follows a recently proposed atypical mechanism that passes through the formation of an N5OOH intermediate at the re side of the cofactor, aided by a well-defined, predominantly hydrophobic O2 pocket. Besides clarifying the unusual chemical mechanism of the complex DszA enzyme, with obvious implications for understanding the puzzling chemistry of flavin-mediated catalysis, the result is crucial for the rational engineering of DszA, contributing to making biodesulfurization attractive for the oil refining industry.

Engineering DszC Mutants from Transition State Macrodipole Considerations and Evolutionary Sequence Analysis.

We describe an approach to identify enzyme mutants with increased turnover using the enzyme DszC as a case study. Our approach is based on recalculating the barriers of alanine mutants through single-point energy calculations at the hybrid QM/MM level in the wild-type reactant and transition state geometries. We analyze the difference in the electron density between the reactant and transition state to identify sites/residues where electrostatic interactions stabilize the transition state over the reactants. We also assess the insertion of a unit probe charge to identify positions in which the introduction of charged residues lowers the barrier.

QM/MM Study of the Reaction Mechanism of Thermophilic Glucuronoyl Esterase for Biomass Treatment.

Hydrolysis of lignocellulosic biomass, composed of a lignin-carbohydrate-complex (LCC) matrix, is critical for producing bioethanol from glucose. However, current methods for LCC processing require costly and polluting processes. The fungal Thermothelomyces thermophila glucuronoyl esterase (TtGE) is a promising thermophilic enzyme that hydrolyses LCC ester bonds. This study describes the TtGE catalytic mechanism using QM/MM methods. Two nearly-degenerate rate-determining transition states were found, with barriers of 16 and 17 kcal ⋅ mol-1 , both with a zwitterionic nature that results from a proton interplay from His346 to either the Ser213-hydroxyl or the lignin leaving group and the rehybridisation of the ester moiety of the substrate to an alkoxide. An oxyanion hole, characteristic of esterases, was provided by the conserved Arg214 through its backbone and sidechain. Our work further suggests that a mutation of Glu267 to a non-negative residue will decrease the energetic barrier in ca. -5 kcal ⋅ mol-1 , improving the catalytic rate of TtGE.

Role of Enzyme and Active Site Conformational Dynamics in the Catalysis by α-Amylase Explored with QM/MM Molecular Dynamics.

We assessed enzyme:substrate conformational dynamics and the rate-limiting glycosylation step of a human pancreatic α-amylase:maltopentose complex. Microsecond molecular dynamics simulations suggested that the distance of the catalytic Asp197 nucleophile to the anomeric carbon of the buried glucoside is responsible for most of the enzyme active site fluctuations and that both Asp197 and Asp300 interact the most with the buried glucoside unit. The buried glucoside binds either in a 4C1 chair or 2SO skew conformations, both of which can change to TS-like conformations characteristic of retaining glucosidases. Starting from four distinct enzyme:substrate complexes, umbrella sampling quantum mechanics/molecular mechanics simulations (converged within less than 1 kcal·mol-1 within a total simulation time of 1.6 ns) indicated that the reaction occurrs with a Gibbs barrier of 13.9 kcal·mol -1, in one asynchronous concerted step encompassing an acid-base reaction with Glu233 followed by a loose SN2-like nucleophilic substitution by the Asp197. The transition state is characterized by a 2H3 half-chair conformation of the buried glucoside that quickly changes to the E3 envelope conformation preceding the attack of the anomeric carbon by the Asp197 nucleophile. Thermodynamic analysis of the reaction supported that a water molecule tightly hydrogen bonded to the glycosidic oxygen of the substrate at the reactant state (∼1.6 Å) forms a short hydrogen bond with Glu233 at the transition state (∼1.7 Å) and lowers the Gibbs barrier in over 5 kcal·mol-1. The resulting Asp197-glycosyl was mostly found in the 4C1 conformation, although the more endergonic B3,O conformation was also observed. Altogether, the combination of short distances for the acid-base reaction with the Glu233 and for the nucleophilic attack by the Asp197 nucleophile and the availability of water within hydrogen bonding distance of the glycosidic oxygen provides a reliable criteria to identify reactive conformations of α-amylase complexes.

Different Enzyme Conformations Induce Different Mechanistic Traits in HIV-1 Protease.

The influence of the dynamical flexibility of enzymes on reaction mechanisms is a cornerstone in biological sciences. In this study, we aim to 1) study the convergence of the activation free energy by using the first step of the reaction catalysed by HIV-1 protease as a case study, and 2) provide further evidence for a mechanistic divergence in this enzyme, as two different reaction pathways were seen to contribute to this step. We used quantum mechanics/molecular mechanics molecular dynamics simulations, on four different initial conformations that led to different barriers in a previous study. Despite the sampling, the four activation free energies still spanned a range of 5.0 kcal ⋅ mol-1 . Furthermore, the new simulations did confirm the occurrence of an unusual mechanistic divergence, with two different mechanistic pathways displaying equivalent barriers. An active-site water molecule is proposed to influence the mechanistic pathway.

Towards the Accurate Thermodynamic Characterization of Enzyme Reaction Mechanisms.

We employed QM/MM molecular dynamics (MD) simulations to characterize the rate-limiting step of the glycosylation reaction of pancreatic α-amylase with combined DFT/molecular dynamics methods (PBE/def2-SVP : AMBER). Upon careful choice of four starting active site conformations based on thorough reactivity criteria, Gibbs energy profiles were calculated with umbrella sampling simulations within a statistical convergence of 1-2 kcal ⋅ mol-1 . Nevertheless, Gibbs activation barriers and reaction energies still varied from 11.0 to 16.8 kcal ⋅ mol-1 and -6.3 to +3.8 kcal ⋅ mol-1 depending on the starting conformations, showing that despite significant state-of-the-art QM/MM MD sampling (0.5 ns/profile) the result still depends on the starting structure. The results supported the one step dissociative mechanism of Asp197 glycosylation preceded by an acid-base reaction by the Glu233, which are qualitatively similar to those from multi-PES QM/MM studies, and thus support the use of the latter to determine enzyme reaction mechanisms.

Animal Fatty Acid Synthase: A Chemical Nanofactory.

Fatty acids are crucial molecules for most living beings, very well spread and conserved across species. These molecules play a role in energy storage, cell membrane architecture, and cell signaling, the latter through their derivative metabolites. De novo synthesis of fatty acids is a complex chemical process that can be achieved either by a metabolic pathway built by a sequence of individual enzymes, such as in most bacteria, or by a single, large multi-enzyme, which incorporates all the chemical capabilities of the metabolic pathway, such as in animals and fungi, and in some bacteria. Here we focus on the multi-enzymes, specifically in the animal fatty acid synthase (FAS). We start by providing a historical overview of this vast field of research. We follow by describing the extraordinary architecture of animal FAS, a homodimeric multi-enzyme with seven different active sites per dimer, including a carrier protein that carries the intermediates from one active site to the next. We then delve into this multi-enzyme's detailed chemistry and critically discuss the current knowledge on the chemical mechanism of each of the steps necessary to synthesize a single fatty acid molecule with atomic detail. In line with this, we discuss the potential and achieved FAS applications in biotechnology, as biosynthetic machines, and compare them with their homologous polyketide synthases, which are also finding wide applications in the same field. Finally, we discuss some open questions on the architecture of FAS, such as their peculiar substrate-shuttling arm, and describe possible reasons for the emergence of large megasynthases during evolution, questions that have fascinated biochemists from long ago but are still far from answered and understood.

Combined in silico and in vitro studies to identify novel antidiabetic flavonoids targeting glycogen phosphorylase.

Novel pharmacological strategies for the treatment of diabetic patients are now focusing on inhibiting glycogenolysis steps. In this regard, glycogen phosphorylase (GP) is a validated target for the discovery of innovative antihyperglycemic molecules. Natural products, and in particular flavonoids, have been reported as potent inhibitors of GP at the cellular level. Herein, free-energy calculations and microscale thermophoresis approaches were performed to get an in-depth assessment of the binding affinities and elucidate intermolecular interactions of several flavonoids at the inhibitor site of GP. To our knowledge, this is the first study indicating genistein, 8-prenylgenistein, apigenin, 8-prenylapigenin, 8-prenylnaringenin, galangin and valoneic acid dilactone as natural molecules with high inhibitory potency toward GP. We identified: i) the residues Phe285, Tyr613, Glu382 and/or Arg770 as the most relevant for the binding of the best flavonoids to the inhibitor site of GP, and ii) the 5-OH, 7-OH, 8-prenyl substitutions in ring A and the 4'-OH insertion in ring B to favor flavonoid binding at this site. Our results are invaluable to plan further structural modifications through organic synthesis approaches and develop more effective pharmaceuticals for Type 2 Diabetes treatment, and serve as the starting point for the exploration of food products for therapeutic usage, as well as for the development of novel bio-functional food and dietary supplements/herbal medicines.

Structural and mechanistic aspects of S-S bonds in the thioredoxin-like family of proteins.

Disulfide bonds play a critical role in a variety of structural and mechanistic processes associated with proteins inside the cells and in the extracellular environment. The thioredoxin family of proteins like thioredoxin (Trx), glutaredoxin (Grx) and protein disulfide isomerase, are involved in the formation, transfer or isomerization of disulfide bonds through a characteristic thiol-disulfide exchange reaction. Here, we review the structural and mechanistic determinants behind the thiol-disulfide exchange reactions for the different enzyme types within this family, rationalizing the known experimental data in light of the results from computational studies. The analysis sheds new atomic-level insight into the structural and mechanistic variations that characterize the different enzymes in the family, helping to explain the associated functional diversity. Furthermore, we review here a pattern of stabilization/destabilization of the conserved active-site cysteine residues presented beforehand, which is fully consistent with the observed roles played by the thioredoxin family of enzymes.

Properties that rank protein:protein docking poses with high accuracy.

The development of docking algorithms to predict near-native structures of protein:protein complexes from the structure of the isolated monomers is of paramount importance for molecular biology and drug discovery. In this study, we assessed the capacity of the interfacial area of protein:protein complexes and of Molecular Mechanics-Poisson Boltzmann Surface Area (MM-PBSA)-derived properties, to rank docking poses. We used a set of 48 protein:protein complexes, and a total of 67 docking experiments distributed among bound:bound, bound:unbound, and unbound:unbound test cases. The MM-PBSA binding free energy of protein monomers has been shown to be very convenient to predict high-quality structures with a high success rate. In fact, considering solely the top-ranked pose of more than 200 docking solutions of each of 39 protein:protein complexes, the success rate was 77% in the prediction of high-quality poses, or 90% if considering high- or medium-quality poses. If considering high- or medium-quality poses as the top-one prediction, a success rate of 87% was obtained for a scoring scheme based on computational alanine scanning mutagenesis data. Such ranking accuracy highlights the ability of these properties to predict near-native poses in protein:protein docking.

Modeling of Human Fatty Acid Synthase and in Silico Docking of Acyl Carrier Protein Domain and Its Partner Catalytic Domains.

Human fatty acid synthase (hFAS) is a megasynthase whose main function is de novo biosynthesis of saturated fatty acids. Interest has been drawn to this enzyme beyond its physiological role due to the association between high levels of hFAS and clinical conditions such as obesity, diabetes, and cancer. Thus, it has become an undeniably attractive pharmacological target. Until now, no crystal structure of the complete hFAS is available, hindering attempts to fully understand this protein. Using homology modeling, we built a model of the entire megasynthase, encompassing all of its domains, including the acyl carrier protein (ACP) and thioesterase (TE) mobile domains absent in the crystal structure of mammalian fatty acid synthase (FAS). On a second stage, we used data-driven protein-protein docking between the substrate shuttling domain ACP and every catalytic domain in the protein. We also propose sets of amino acids at the interface of each domain that we believe are important to favor the interaction between ACP and each domain of hFAS. After inspection, we validated each complex between ACP and MAT/KS/KR/DH/ER domains through classical molecular dynamics simulations and RMSd analysis. Additionally, we mapped the interactions between the residues at the active site of each catalytic domain and its intermediaries. In every docking, we ensured that the distance between catalytic residues and the intermediaries was maintained. Until now, there was not a complete 3D model of this megasynthase. This study is the first to present a homology model for the whole hFAS, including its two mobile domains and possible poses of ACP throughout the cycle of fatty acid biosynthesis, thus mapping obligatory checkpoints in its trajectory. Hence, we believe that these structural insights will allow for future studies of the catalytic mechanism of the overall hFAS.

Mechanistic insights on the reduction of glutathione disulfide by protein disulfide isomerase.

We explore the enzymatic mechanism of the reduction of glutathione disulfide (GSSG) by the reduced a domain of human protein disulfide isomerase (hPDI) with atomistic resolution. We use classical molecular dynamics and hybrid quantum mechanics/molecular mechanics calculations at the mPW1N/6-311+G(2d,2p):FF99SB//mPW1N/6-31G(d):FF99SB level. The reaction proceeds in two stages: (i) a thiol-disulfide exchange through nucleophilic attack of the Cys53-thiolate to the GSSG-disulfide followed by the deprotonation of Cys56-thiol by Glu47-carboxylate and (ii) a second thiol-disulfide exchange between the Cys56-thiolate and the mixed disulfide intermediate formed in the first step. The Gibbs activation energy for the first stage was 18.7 kcal·mol-1, and for the second stage, it was 7.2 kcal·mol-1, in excellent agreement with the experimental barrier (17.6 kcal·mol-1). Our results also suggest that the catalysis by protein disulfide isomerase (PDI) and thiol-disulfide exchange is mostly enthalpy-driven (entropy changes below 2 kcal·mol-1 at all stages of the reaction). Hydrogen bonds formed between the backbone of His55 and Cys56 and the Cys56-thiol result in an increase in the Gibbs energy barrier of the first thiol-disulfide exchange. The solvent plays a key role in stabilizing the leaving glutathione thiolate formed. This role is not exclusively electrostatic, because an explicit inclusion of several water molecules at the density-functional theory level is a requisite to form the mixed disulfide intermediate. In the intramolecular oxidation of PDI, a transition state is only observed if hydrogen bond donors are nearby the mixed disulfide intermediate, which emphasizes that the thermochemistry of thiol-disulfide exchange in PDI is influenced by the presence of hydrogen bond donors.

A QM/MM study of the reaction mechanism of human ß-ketoacyl reductase.

Human fatty acid synthase (hFAS) is a multifunctional enzyme involved in a wide diversity of biological functions. For instance, it is a precursor of phospholipids and other complex processes such as the de novo synthesis of long chain fatty acid. Human FAS is also a component of biological membranes and it is implicated in the overexpression of several types of cancers. In this work, we describe the catalytic mechanism of ß-ketoreductase (KR), which is a catalytic domain of the hFAS enzyme that catalyzes the reduction of ß-ketoacyl to ß-hydroxyacyl with the concomitant oxidation of the NADPH cofactor. The catalysis by KR is an intermediate step in the cycle of reactions that elongate the substrate's carbon chain until the final product is obtained. We study and propose the catalytic mechanism of the KR domain determined using the hybrid QM/MM methodology, at the ONIOM(B3LYP/6-311+G(2d,2p):AMBER) level of theory. The results indicate that the reaction mechanism occurs in two stages: (i) nucleophilic attack by a NADPH hydride to the ß-carbon of the substrate, together with an asynchronous deprotonation of the Tyr2034 by the oxygen of the ß-alkoxide to hold the final alcohol product; and (ii) an asynchronous deprotonation of the hydroxyl in the NADP+'s ribose by Tyr2034, and of the Lys1995 by the resulting alkoxide in the former ribose to restore the protonation state of Tyr2034. The reduction step occurs with a Gibbs energy barrier of 11.7 kcal mol-1 and a Gibbs reaction energy of -10.6 kcal mol-1. These results have provided an understanding of the catalytic mechanism of the KR hFAS domain, a piece of the heavy hFAS biosynthetic machinery.

Improving the Biodesulfurization of Crude Oil and Derivatives: A QM/MM Investigation of the Catalytic Mechanism of NADH-FMN Oxidoreductase (DszD).

The development of biocatalytic desulfurization strategies of petroleum and its derivatives could result in more economic alternatives than the widely used chemical desulfurization. The organism Rhodococcus erythropolis IGTS8 has been shown to metabolize organic sulfur compounds through a mechanism known as 4S pathway, which involves four enzymes (DszA, DszB, DszC, and DszD) and has been explored in biodesulfurization. Here we have applied QM/MM methods to study the catalytic mechanism of the enzyme DszD, a NADH-FMN oxidoreductase that occupies a central place on the 4S pathway by catalyzing the formation of the FMNH2 that is used by the two monooxynases in the cycle: DszA and DszC. In addition, to clarify the catalytic mechanism of this enzyme, this study analyzed in detail the role played by the active site Thr residue and of Asn and Ala enzyme mutants. The results help to explain previous experimental evidence and suggest new strategies for improving biodesulfurization through an increase in the activity of DszD.

Benchmarking of Density Functionals for the Accurate Description of Thiol-Disulfide Exchange.

A set of 92 density functionals was employed to accurately characterize thiol-disulfide exchange. The properties we have benchmarked throughout the study include the geometry of a 15 atoms model system, the potential energy surface, the activation barrier, and the energy of reaction for thiol-disulfide exchange. Reference energies were determined at the CCSD(T)/CBS//MP2/aug-cc-pVDZ level of theory, and reference geometries were calculated at the MP2/aug-cc-pVTZ level. M11-L, M06-2X, M06-HF, N12-SX, PBE1PBE, PBEh1PBE, and OHSE2PBE described better the geometry of the model system, with average deviations of 0.06 Å in bond lengths (0.06 Å in bond-breaking lengths) and 1.9° in bond angles. On the other hand, the potential energy surface and its gradient were more accurately described by the hybrid density functional BHandH, closely followed by mPW1N, mPW1K, and mPWB1K. The barrier height and energy of reaction were better reproduced by the BMK and M06-2X functionals (deviations of 0.17 and 0.07 kcal·mol(-1), respectively) for a set of 10 Pople's basis sets. MN12-SX and M11-L showed very good results for the widely used 6-311++G(2d,2p) basis set, with deviations of 0.02 and 0.05 kcal·mol(-1), respectively. We studied the effect of the split-valence, diffuse, and polarized functions in the activation barrier of thiol-disulfide exchange, for a set of 10 Pople's basis sets. While increasing the splitting and polarization may increase the activation barrier in approximately 1 kcal·mol(-1), diffuse functions generally contribute to decreasing it no more than 0.10 kcal·mol(-1). In general, 13 functionals provided energies within 1 kcal·mol(-1) of the reference value. The BB1K density functional is one of the best density functionals to characterize thiol-disulfide exchange reactions; however, several density functionals with modified Perdew-Wang exchange and about 40% Hartree-Fock exchange, such as mPW1K, mPW1N, and mPWB1K, show a good performance, too.

Parameters for Molecular Dynamics Simulations of Manganese-Containing Metalloproteins.

A set of geometrical parameters has been determined for single manganese metalloproteins for the AMBER force field, and ultimately to other force fields with a similar philosophy. Twelve (12) models from 9 different single-cluster manganese proteins were optimized and parametrized, using a bonded model approach. Mn-ligand bonds, Mn-ligand angles, and Restrained Electrostatic Potential charges for all the 74 residues in the first metal coordination sphere of each Mn metalloprotein were parametrized. The determined parameters were validated with molecular dynamics simulations and several statistics strategies were used to analyze the results. In addition, to validate the parametrized models, frequency and normal mode calculations were performed and comparisons were obtained for the overall structures both with quantum mechanics and molecular mechanics calculations. Linear and polynomial fittings were performed to estimate Mn-ligand bond force constants for generic manganese centers. Furthermore, averages are proposed for the main Mn-ligand angle interactions of typical manganese coordination centers: axial, square and triangular equatorial planes, and tetrahedral positions, for the different combinations of donor atoms from waters and hard ligands.