Fernanda Duarte

Force Field Independent Metal Parameters Using a Nonbonded Dummy Model

The cationic dummy atom approach provides a powerful nonbonded description for a range of alkaline-earth and transition-metal centers, capturing both structural and electrostatic effects. In this work we refine existing literature parameters for octahedrally coordinated Mn2+, Zn2+, Mg2+, and Ca2+, as well as providing new parameters for Ni2+, Co2+, and Fe2+. In all the cases, we are able to reproduce both M2+–O distances and experimental solvation free energies, which has not been achieved to date for transition metals using any other model. The parameters have also been tested using two different water models and show consistent performance. Therefore, our parameters are easily transferable to any force field that describes nonbonded interactions using Coulomb and Lennard-Jones potentials. Finally, we demonstrate the stability of our parameters in both the human and Escherichia coli variants of the enzyme glyoxalase I as showcase systems, as both enzymes are active with a range of transition metals. The parameters presented in this work provide a valuable resource for the molecular simulation community, as they extend the range of metal ions that can be studied using classical approaches, while also providing a starting point for subsequent parametrization of new metal centers.

Resolving apparent conflicts between theoretical and experimental models of phosphate monoester hydrolysis.

Understanding phosphoryl and sulfuryl transfer is central to many biochemical processes. However, despite decades of experimental and computational studies, a consensus concerning the precise mechanistic details of these reactions has yet to be reached. In this work we perform a detailed comparative theoretical study of the hydrolysis of p-nitrophenyl phosphate, methyl phosphate and p-nitrophenyl sulfate, all of which have served as key model systems for understanding phosphoryl and sulfuryl transfer reactions, respectively. We demonstrate the existence of energetically similar but mechanistically distinct possibilities for phosphate monoester hydrolysis. The calculated kinetic isotope effects for p-nitrophenyl phosphate provide a means to discriminate between substrate- and solvent-assisted pathways of phosphate monoester hydrolysis, and show that the solvent-assisted pathway dominates in solution. This preferred mechanism for p-nitrophenyl phosphate hydrolysis is difficult to find computationally due to the limitations of compressing multiple bonding changes onto a 2-dimensional energy surface. This problem is compounded by the need to include implicit solvation to at least microsolvate the system and stabilize the highly charged species. In contrast, methyl phosphate hydrolysis shows a preference for a substrate-assisted mechanism. For p-nitrophenyl sulfate hydrolysis there is only one viable reaction pathway, which is similar to the solvent-assisted pathway for phosphate hydrolysis, and the substrate-assisted pathway is not accessible. Overall, our results provide a unifying mechanistic framework that is consistent with the experimentally measured kinetic isotope effects and reconciles the discrepancies between theoretical and experimental models for these biochemically ubiquitous classes of reaction.

Cooperative Electrostatic Interactions Drive Functional Evolution in the Alkaline Phosphatase Superfamily

It is becoming widely accepted that catalytic promiscuity, i.e., the ability of a single enzyme to catalyze the turnover of multiple, chemically distinct substrates, plays a key role in the evolution of new enzyme functions. In this context, the members of the alkaline phosphatase superfamily have been extensively studied as model systems in order to understand the phenomenon of enzyme multifunctionality. In the present work, we model the selectivity of two multiply promiscuous members of this superfamily, namely the phosphonate monoester hydrolases from Burkholderia caryophylli and Rhizobium leguminosarum. We have performed extensive simulations of the enzymatic reaction of both wild-type enzymes and several experimentally characterized mutants. Our computational models are in agreement with key experimental observables, such as the observed activities of the wild-type enzymes, qualitative interpretations of experimental pH-rate profiles, and activity trends among several active site mutants. In all cases the substrates of interest bind to the enzyme in similar conformations, with largely unperturbed transition states from their corresponding analogues in aqueous solution. Examination of transition-state geometries and the contribution of individual residues to the calculated activation barriers suggest that the broad promiscuity of these enzymes arises from cooperative electrostatic interactions in the active site, allowing each enzyme to adapt to the electrostatic needs of different substrates. By comparing the structural and electrostatic features of several alkaline phosphatases, we suggest that this phenomenon is a generalized feature driving selectivity and promiscuity within this superfamily and can be in turn used for artificial enzyme design.

Empirical valence bond simulations of the hydride transfer step in the monoamine oxidase B catalyzed metabolism of dopamine: EVB Simulations of Hydride Transfer in MAO B

Monoamine oxidases (MAOs) A and B are flavoenzymes responsible for the metabolism of biogenic amines such as dopamine, serotonin and noradrenaline. In this work, we present a comprehensive study of the rate‐limiting step of dopamine degradation by MAO B, which consists in the hydride transfer from the methylene group of the substrate to the flavin moiety of the FAD prosthetic group. This article builds on our previous quantum chemical study of the same reaction using a cluster model (Vianello et al., Eur J Org Chem 2012; 7057), but now considering the full dimensionality of the hydrated enzyme with extensive configurational sampling. We show that MAO B is specifically tuned to catalyze the hydride transfer step from the substrate to the flavin moiety of the FAD prosthetic group and that it lowers the activation barrier by 12.3 kcal mol−1 compared to the same reaction in aqueous solution, a rate enhancement of more than nine orders of magnitude. Taking into account the deprotonation of the substrate prior to the hydride transfer reaction, the activation barrier in the enzyme is calculated to be 16.1 kcal mol−1, in excellent agreement with the experimental value of 16.5 kcal mol−1. Additionally, we demonstrate that the protonation state of the active site residue Lys296 does not have an influence on the hydride transfer reaction. Proteins 2014; 82:3347–3355.

How does Pin1 catalyze the cis-trans prolyl peptide bond isomerization? A QM/MM and mean reaction force study.

Pin1 represents an enzyme that specifically catalyzes the isomerization of peptide bonds between phosphorylated threonine or serine residues and proline. Despite its relevance as molecular timer in a number of biological processes related to cancer and Alzheimer disease, a detailed understanding of the factors contributing to the catalysis is still missing. In this study, we employ extensive QM/MM molecular dynamics simulations in combination with the mean reaction force (MRF) to discern the influence of the enzyme on the reaction mechanism and the origin of the catalysis. As a recently introduced method, the MRF separates the activation free energy barrier to reach the transition state into structural and electronic contributions providing a more detailed description of the enzymes function. As a reference, we first study the isomerization starting from the cis form in solution and obtain a free energy barrier and a reaction free energy, which are in agreement with previous studies and experiment. With the new mean reaction force method, intramolecular hydrogen bonds in the peptide were identified that stabilize the transition state and reduce the electronic contribution to the free energy barrier. To elucidate the mechanism of catalysis of Pin1, the reaction in solution and in the catalytic cavity of the enzyme were compared. Both yield the same free energy barrier for the isomerization of the cis form, but with different decomposition in structural and electronic contributions by the mean reaction force. The enzyme reduces the energy required for structural rearrangements to reach the transition state, pointing to a destabilization of the reactant, but increases the electronic contribution to the barrier through specific enzyme-peptide hydrogen bonds. In the reverse reaction, the isomerization of the trans form, the enzyme alters the energetics and the mechanism of the reaction considerably. Unfavorable enzyme-peptide interactions in the catalytic cavity during the isomerization change the reaction coordinate, resulting in two minima with small energy differences to the transition state. These small free energy barriers should in principle make the reaction feasible at room temperature once the conformer is bound in the right conformation.

The Alkaline Hydrolysis of Sulfonate Esters: Challenges in Interpreting Experimental and Theoretical Data

Sulfonate ester hydrolysis has been the subject of recent debate, with experimental evidence interpreted in terms of both stepwise and concerted mechanisms. In particular, a recent study of the alkaline hydrolysis of a series of benzene arylsulfonates (Babtie et al., Org. Biomol. Chem.10, 2012, 8095) presented a nonlinear Brønsted plot, which was explained in terms of a change from a stepwise mechanism involving a pentavalent intermediate for poorer leaving groups to a fully concerted mechanism for good leaving groups and supported by a theoretical study. In the present work, we have performed a detailed computational study of the hydrolysis of these compounds and find no computational evidence for a thermodynamically stable intermediate for any of these compounds. Additionally, we have extended the experimental data to include pyridine-3-yl benzene sulfonate and its N-oxide and N-methylpyridinium derivatives. Inclusion of these compounds converts the Brønsted plot to a moderately scattered but linear correlation and gives a very good Hammett correlation. These data suggest a concerted pathway for this reaction that proceeds via an early transition state with little bond cleavage to the leaving group, highlighting the care that needs to be taken with the interpretation of experimental and especially theoretical data.

The Competing Mechanisms of Phosphate Monoester Dianion Hydrolysis

Despite the numerous experimental and theoretical studies on phosphate monoester hydrolysis, significant questions remain concerning the mechanistic details of these biologically critical reactions. In the present work we construct a linear free energy relationship for phosphate monoester hydrolysis to explore the effect of modulating leaving group pKa on the competition between solvent- and substrate-assisted pathways for the hydrolysis of these compounds. Through detailed comparative electronic-structure studies of methyl phosphate and a series of substituted aryl phosphate monoesters, we demonstrate that the preferred mechanism is dependent on the nature of the leaving group. For good leaving groups, a strong preference is observed for a more dissociative solvent-assisted pathway. However, the energy difference between the two pathways gradually reduces as the leaving group pKa increases and creates mechanistic ambiguity for reactions involving relatively poor alkoxy leaving groups. Our calculations show that the transition-state structures vary smoothly across the range of pKas studied and that the pathways remain discrete mechanistic alternatives. Therefore, while not impossible, a biological catalyst would have to surmount a significantly higher activation barrier to facilitate a substrate-assisted pathway than for the solvent-assisted pathway when phosphate is bonded to good leaving groups. For poor leaving groups, this intrinsic preference disappears.

Concerted or Stepwise: How Much Do Free-Energy Landscapes Tell Us about the Mechanisms of Elimination Reactions?

The base-catalyzed dehydration of benzene cis-1,2-dihydrodiols is driven by formation of an aromatic product as well as intermediates potentially stabilized by hyperaromaticity. Experiments exhibit surprising shifts in isotope effects, indicating an unusual mechanistic balance on the E2-E1cB continuum. In this study, both 1- and 2-dimensional free energy surfaces are generated for these compounds with various substituents, using density functional theory and a mixed implicit/explicit solvation model. The computational data help unravel hidden intermediates along the reaction coordinate and provide a novel conceptual framework for distinguishing between competing pathways in this and any other system with borderline reaction mechanisms.

Promiscuity in the Enzymatic Catalysis of Phosphate and Sulfate Transfer

The enzymes that facilitate phosphate and sulfate hydrolysis are among the most proficient natural catalysts known to date. Interestingly, a large number of these enzymes are promiscuous catalysts that exhibit both phosphatase and sulfatase activities in the same active site and, on top of that, have also been demonstrated to efficiently catalyze the hydrolysis of other additional substrates with varying degrees of efficiency. Understanding the factors that underlie such multifunctionality is crucial both for understanding functional evolution in enzyme superfamilies and for the development of artificial enzymes. In this Current Topic, we have primarily focused on the structural and mechanistic basis for catalytic promiscuity among enzymes that facilitate both phosphoryl and sulfuryl transfer in the same active site, while comparing this to how catalytic promiscuity manifests in other promiscuous phosphatases. We have also drawn on the large number of experimental and computational studies of selected model systems in the literature to explore the different features driving the catalytic promiscuity of such enzymes. Finally, on the basis of this comparative analysis, we probe the plausible origins and determinants of catalytic promiscuity in enzymes that catalyze phosphoryl and sulfuryl transfer.

Molecular Recognition in Asymmetric Counteranion Catalysis: Understanding Chiral Phosphate-Mediated Desymmetrization

We describe the first theoretical study of a landmark example of chiral anion phase-transfer catalysis: the enantioselective ring-opening of meso-aziridinium and episulfonium cations promoted by asymmetric counteranion-directed catalysis (ACDC). The mechanism of ion-pairing, ring-opening, and catalyst deactivation have been studied in the condensed phase with both classical and quantum methods using explicitly and implicitly solvated models. We find that the stability of chiral ion-pairs, a prerequisite for asymmetric catalysis, is dominated by electrostatic interactions at long range and by CH···O interactions at short range. The decisive role of solvent upon ion-pair formation and of nonbonding interactions upon enantioselectivity are quantified by complementary computational approaches. The major enantiomer is favored by a smaller distortion of the substrate, demonstrated by a distortion/interaction analysis. Our computational results rationalize the stereoselectivity for several experimental results and demonstrate a combined classical/quantum approach to perform realistic modeling of chiral counterion catalysis in solution.