Anomeric Effect in “High Energy” Phosphate Bonds. Selective Destabilization of the Scissile Bond and Modulation of the Exothermicity of Hydrolysis

The Complex Roles of Adenosine Triphosphate in Bioenergetics

ATP is generally defined as the "energy currency" of the cell. Its phosphoanhydride P-O bonds are often considered to be "high energy" linkages that release free energy when broken, and its hydrolysis is described as "strongly exergonic". However, breaking bonds cannot release energy and ATP hydrolysis in motor and active transport proteins is not "strongly exergonic". So, the relevance of ATP resides elsewhere. As important as the nucleotide are the proteins that undergo functionally relevant conformational changes upon both ATP binding and release of ADP and inorganic phosphate. ATP phosphorylates proteins for signaling, active transport, and substrates in condensation reactions. The ensuing dephosphorylation has different consequences in each case. In signaling and active transport the phosphate group is hydrolyzed whereas in condensation reactions the phosphoryl fragment acts as a dehydrating agent. As it will be discussed in this article, ATP does much more than simply contribute free energy to biological processes.

Anomeric effect, hyperconjugation and electrostatics: lessons from complexity in a classic stereoelectronic phenomenon

Understanding the interplay of multiple components (steric, electrostatic, stereoelectronic, dispersive, etc.) that define the overall energy, structure, and reactivity of organic molecules can be a daunting task. The task becomes even more difficult when multiple approaches based on different physical premises disagree in their analysis of a multicomponent molecular system. Herein, we will use a classic conformational "oddity", the anomeric effect, to discuss the value of identifying the key contributors to reactivity that can guide chemical predictions. After providing the background related to the relevant types of hyperconjugation and a brief historic outline of the origins of the anomeric effect, we outline variations of its patterns and provide illustrative examples for the role of the anomeric effect in structure, stability, and spectroscopic properties. We show that the complete hyperconjugative model remains superior in explaining the interplay between structure and reactivity. We will use recent controversies regarding the origin of the anomeric effect to start a deeper discussion relevant to any electronic effect. Why are such questions inherently controversial? How to describe a complex quantum system using a model that is "as simple as possible, but no simpler"? What is a fair test for such a model? Perhaps, instead of asking "who is right and who is wrong?" one should ask "why do we disagree?". Stereoelectronic thinking can reconcile quantum complexity with chemical intuition and build the conceptual bridge between structure and reactivity. Even when many factors contribute to the observed structural and conformational trends, electron delocalization is a dominating force when the electronic demand is high (i.e., bonds are breaking as molecules distort from their equilibrium geometries). In these situations, the role of orbital interactions increases to the extent where they can define reactivity. For example, negative hyperconjugation can unleash the "underutilized" stereoelectronic power of unshared electrons (i.e., the lone pairs) to stabilize a developing positive charge at an anomeric carbon. This analysis paves the way for the broader discussion of the omnipresent importance of negative hyperconjugation in oxygen-containing functional groups. From that point of view, the stereoelectronic component of the anomeric effect plays a unique role in guiding reaction design.

Stereoelectronic power of oxygen in control of chemical reactivity: the anomeric effect is not alone

Although carbon is the central element of organic chemistry, oxygen is the central element of stereoelectronic control in organic chemistry. Generally, a molecule with a C-O bond has both a strong donor (a lone pair) and a strong acceptor (e.g., a σ*_C-O orbital), a combination that provides opportunities to influence chemical transformations at both ends of the electron demand spectrum. Oxygen is a stereoelectronic chameleon that adapts to the varying situations in radical, cationic, anionic, and metal-mediated transformations. Arguably, the most historically important stereoelectronic effect is the anomeric effect (AE), i.e., the axial preference of acceptor groups at the anomeric position of sugars. Although AE is generally attributed to hyperconjugative interactions of σ-acceptors with a lone pair at oxygen (negative hyperconjugation), recent literature reports suggested alternative explanations. In this context, it is timely to evaluate the fundamental connections between the AE and a broad variety of O-functional groups. Such connections illustrate the general role of hyperconjugation with oxygen lone pairs in reactivity. Lessons from the AE can be used as the conceptual framework for organizing disjointed observations into a logical body of knowledge. In contrast, neglect of hyperconjugation can be deeply misleading as it removes the stereoelectronic cornerstone on which, as we show in this review, the chemistry of organic oxygen functionalities is largely based. As negative hyperconjugation releases the "underutilized" stereoelectronic power of unshared electrons (the lone pairs) for the stabilization of a developing positive charge, the role of orbital interactions increases when the electronic demand is high and molecules distort from their equilibrium geometries. From this perspective, hyperconjugative anomeric interactions play a unique role in guiding reaction design. In this manuscript, we discuss the reactivity of organic O-functionalities, outline variations in the possible hyperconjugative patterns, and showcase the vast implications of AE for the structure and reactivity. On our journey through a variety of O-containing organic functional groups, from textbook to exotic, we will illustrate how this knowledge can predict chemical reactivity and unlock new useful synthetic transformations.

Analysis of Density Functional Tight Binding with Natural Bonding Orbitals

The description of chemical bonding by the density functional tight binding (DFTB) model is analyzed using natural bonding orbitals (NBOs) and compared to results from density functional theory (B3LYP/aug-cc-pVTZ) calculations. Several molecular systems have been chosen to represent fairly diverse bonding scenarios that include standard covalent bonds, hypervalent interactions, multicenter bonds, metal-ligand interactions (with and without the pseudo-Jahn-Teller effect), and through-space donor-acceptor interactions. Overall, the results suggest that DFTB3/3OB provides physically sound descriptions for the different bonding scenarios analyzed here, as reflected by the general agreement between DFTB3 and B3LYP NBO properties, such as the nature of the NBOs, the magnitudes of natural charges and bond orders, and the dominant donor-acceptor interactions. The degree of ligand-to-metal charge transfer and the ionic nature of pentavalent phosphate are overestimated, likely reflecting the minimal-basis nature of DFTB3/3OB. Moreover, certain orbital interactions, such as geminal interactions, are observed to be grossly overestimated by DFTB3 for hypervalent phosphate and several transition metal compounds that involve copper and nickel. The study indicates that results from NBO analysis can be instructive for identifying electronic structure descriptions at the approximate quantum-mechanical level that require improvement and thus for guiding the systematic improvement of these methods.

Drastic Compensation of Electronic and Solvation Effects on ATP Hydrolysis Revealed through Large-Scale QM/MM Simulations Combined with a Theory of Solutions

Hydrolysis of adenosine triphosphate (ATP) is the "energy source" for a variety of biochemical processes. In the present work, we address key features of ATP hydrolysis: the relatively moderate value (about -10 kcal/mol) of the standard free energy, ΔG_hyd, of reaction and the insensitivity of ΔG_hyd to the number of excess electrons on ATP. We conducted quantum mechanical/molecular mechanical simulation combined with the energy-representation theory of solutions to analyze the electronic-state and solvation contributions to ΔG_hyd. It was revealed that the electronic-state contribution in ΔG_hyd is largely negative (favorable) upon hydrolysis, due to the reduction of electrostatic repulsion accompanying the breakage of the P-O bond. In contrast, the solvation effect was found to be strongly more favorable on the reactant side. Thus, we showed that a drastic compensation of the two opposite effects takes place, leading to the modest value of ΔG_hyd at each number of excess electrons examined. The computational analyses were also conducted for pyrophosphate ions (PPi), and the parallelism between the ATP and PPi hydrolyses was confirmed. Classical molecular dynamics simulation was further carried out to discuss the effect of the solvent environment; the insensitivity of ΔG_hyd to the number of excess electrons was seen to hold in solvent water and ethanol.

Use of restrained molecular dynamics to predict the conformations of phosphorylated receiver domains in two-component signaling systems

Two‐component signaling (TCS) is the primary means by which bacteria, as well as certain plants and fungi, respond to external stimuli. Signal transduction involves stimulus‐dependent autophosphorylation of a sensor histidine kinase and phosphoryl transfer to the receiver domain of a downstream response regulator. Phosphorylation acts as an allosteric switch, inducing structural and functional changes in the pathway's components. Due to their transient nature, phosphorylated receiver domains are challenging to characterize structurally. In this work, we provide a methodology for simulating receiver domain phosphorylation to predict conformations that are nearly identical to experimental structures. Using restrained molecular dynamics, phosphorylated conformations of receiver domains can be reliably sampled on nanosecond timescales. These simulations also provide data on conformational dynamics that can be used to identify regions of functional significance related to phosphorylation. We first validated this approach on several well‐characterized receiver domains and then used it to compare the upstream and downstream components of the fungal Sln1 phosphorelay. Our results demonstrate that this technique provides structural insight, obtained in the absence of crystallographic or NMR information, regarding phosphorylation‐induced conformational changes in receiver domains that regulate the output of their associated signaling pathway. To our knowledge, this is the first time such a protocol has been described that can be broadly applied to TCS proteins for predictive purposes. Proteins 2016; 85:155–176. © 2016 Wiley Periodicals, Inc.

Supramolecular phosphate transfer catalysis by pillar[5]arene

Phosphate ester destabilization inside pillar[5]arene due to disrupting the H-bonding to water speed its hydrolysis reaction.

Restricting the ψ Torsion Angle Has Stereoelectronic Consequences on a Scissile Bond: An Electronic Structure Analysis

Protein motion is intimately linked to enzymatic catalysis, yet the stereoelectronic changes that accompany different conformational states of a substrate are poorly defined. Here we investigate the relationship between conformation and stereoelectronic effects of a scissile amide bond. Structural studies have revealed that the C-terminal glycine of ubiquitin and ubiquitin-like proteins adopts a syn (ψ ∼ 0°) or gauche (ψ ∼ ±60°) conformation upon interacting with deubiquitinases/ubiquitin-like proteases. We used hybrid density functional theory and natural bond orbital analysis to understand how the stereoelectronic effects of the scissile bond change as a function of φ and ψ torsion angles. This led to the discovery that when ψ is between 30° and -30° the scissile bond becomes geometrically and electronically deformed. Geometric distortion occurs through pyramidalization of the carbonyl carbon and amide nitrogen. Electronic distortion is manifested by a decrease in the strength of the donor-acceptor interaction between the amide nitrogen and antibonding orbital (π*) of the carbonyl. Concomitant with the reduction in nN → π* delocalization energy, the sp(2) hybrid orbital of the carbonyl carbon becomes richer in p-character, suggesting the syn configuration causes the carbonyl carbon hybrid orbitals to adopt a geometry reminiscent of a tetrahedral-like intermediate. Our work reveals important insights into the role of substrate conformation in activating the reactive carbonyl of a scissile bond. These findings have implications for designing potent active site inhibitors based on the concept of transition state analogues.

Fundamentals of phosphate transfer

Historically, the chemistry of phosphate transfer-a class of reactions fundamental to the chemistry of Life-has been discussed almost exclusively in terms of the nucleophile and the leaving group. Reactivity always depends significantly on both factors; but recent results for reactions of phosphate triesters have shown that it can also depend strongly on the nature of the nonleaving or "spectator" groups. The extreme stabilities of fully ionised mono- and dialkyl phosphate esters can be seen as extensions of the same effect, with one or two triester OR groups replaced by O(-). Our chosen lead reaction is hydrolysis-phosphate transfer to water: because water is the medium in which biological chemistry takes place; because the half-life of a system in water is an accepted basic index of stability; and because the typical mechanisms of hydrolysis, with solvent H2O providing specific molecules to act as nucleophiles and as general acids or bases, are models for reactions involving better nucleophiles and stronger general species catalysts. Not least those available in enzyme active sites. Alkyl monoester dianions compete with alkyl diester monoanions for the slowest estimated rates of spontaneous hydrolysis. High stability at physiological pH is a vital factor in the biological roles of organic phosphates, but a significant limitation for experimental investigations. Almost all kinetic measurements of phosphate transfer reactions involving mono- and diesters have been followed by UV-visible spectroscopy using activated systems, conveniently compounds with good leaving groups. (A "good leaving group" OR* is electron-withdrawing, and can be displaced to generate an anion R*O(-) in water near pH 7.) Reactivities at normal temperatures of P-O-alkyl derivatives-better models for typical biological substrates-have typically had to be estimated: by extended extrapolation from linear free energy relationships, or from rate measurements at high temperatures. Calculation is free from these limitations, able to handle very slow reactions as readily as very fast ones, and capable of predicting rate constants with levels of accuracy acceptable to the experimentalist. We present an updated overview of phosphate transfer, with particular reference to the mechanisms of the reactions of alkyl derivatives and triesters. The intention is to present a holistic (not comprehensive!) overview of the reactivity of typical phosphate esters, in terms familiar to the working chemist, at a level sufficient to support informed predictions of reactivity for structures of interest.

Charge density distribution and theoretical analysis of low and high energy phosphate esters

There is a strict relation between the energy of hydrolysis of phosphate esters and the extent of interactions between the p ester oxygen lone pair and the antibonding orbitals of the rest of the molecule. Its impact on experimental charge density distribution is analyzed.

Common hydrogen bond interactions in diverse phosphoryl transfer active sites

Phosphoryl transfer reactions figure prominently in energy metabolism, signaling, transport and motility. Prior detailed studies of selected systems have highlighted mechanistic features that distinguish different phosphoryl transfer enzymes. Here, a top-down approach is developed for comparing statistically the active site configurations between populations of diverse structures in the Protein Data Bank, and it reveals patterns of hydrogen bonding that transcend enzyme families. Through analysis of large samples of structures, insights are drawn at a level of detail exceeding the experimental precision of an individual structure. In phosphagen kinases, for example, hydrogen bonds with the O3β of the nucleotide substrate are revealed as analogous to those in unrelated G proteins. In G proteins and other enzymes, interactions with O3β have been understood in terms of electrostatic favoring of the transition state. Ground state quantum mechanical calculations on model compounds show that the active site interactions highlighted in our database analysis can affect substrate phosphate charge and bond length, in ways that are consistent with prior experimental observations, by modulating hyperconjugative orbital interactions that weaken the scissile bond. Testing experimentally the inference about the importance of O3β interactions in phosphagen kinases, mutation of arginine kinase Arg280 decreases kcat, as predicted, with little impact upon KM.

Hyperconjugation-mediated solvent effects in phosphoanhydride bonds

Density functional theory and natural bond orbital analysis are used to explore the impact of solvent on hyperconjugation in methyl triphosphate, a model for "energy rich" phosphoanhydride bonds, such as found in ATP. As expected, dihedral rotation of a hydroxyl group vicinal to the phosphoanhydride bond reveals that the conformational dependence of the anomeric effect involves modulation of the orbital overlap between the donor and acceptor orbitals. However, a conformational independence was observed in the rotation of a solvent hydrogen bond. As one lone pair orbital rotates away from an optimal antiperiplanar orientation, the overall magnitude of the anomeric effect is compensated approximately by the other lone pair as it becomes more antiperiplanar. Furthermore, solvent modulation of the anomeric effect is not restricted to the antiperiplanar lone pair; hydrogen bonds involving gauche lone pairs also affect the anomeric interaction and the strength of the phosphoanhydride bond. Both gauche and anti solvent hydrogen bonds lengthen nonbridging O-P bonds, increasing the distance between donor and acceptor orbitals and decreasing orbital overlap, which leads to a reduction of the anomeric effect. Solvent effects are additive with greater reduction in the anomeric effect upon increasing water coordination. By controlling the coordination environment of substrates in an active site, kinases, phosphatases, and other enzymes important in metabolism and signaling may have the potential to modulate the stability of individual phosphoanhydride bonds through stereoelectronic effects.

Mechanistic insights into the hydrolysis of a nucleoside triphosphate model in neutral and acidic solution

Nucleoside triphosphate hydrolysis is an essential component of all living systems. Despite extensive research, the exact modus and mechanism of this ubiquitous reaction still remain elusive. In this work, we examined the detailed hydrolysis mechanisms of a model nucleoside triphosphate in acidic and neutral solution by means of ab initio simulations. The timescale of the reaction was accessed through use of an accelerated sampling method, metadynamics. Both hydrolyses were found to proceed via different mechanisms; the acidic system reacted by means of concerted general acid catalysis (found to be a so-called D(N)A(N)A(H)D(xh) mechanism), whereas the neutral system reacted by way of a different mechanism (namely, D(N)*A(N)D(xh)A(H)). A neighboring water molecule took on the role of a general base in both systems, which has not been seen before but is a highly plausible reaction path, meaning that substrate-assisted catalysis was not observed in the bulk water environment.

Influence of the molecular environment on phosphorylated amino acid models: a density functional theory study

A protein environment can affect the structure and charge distribution of substrate molecules. Here, the structure and partial charges were studied for different phosphorylated amino acid models in varying environments using density functional theory. The three systems investigated, acetyl phosphate, methyl phosphate, and p-tolyl phosphate are representative models for aspartyl phosphate, serine or threonine phosphate, and tyrosine phosphate, respectively. Combined with the CPCM continuum model, explicit HF and H(2)O molecules were added in order to model environmental effects and interactions that may occur in a protein matrix. We show how the different interactions affect the scissile P-O(R) bond and that the elongation can be explained by an anomeric effect. An increasing scissile bond length will result in transfer of negative charge to the leaving group and in a widening of the angle between the terminal oxygens of the phosphate molecule, features that can expose the phosphate group to attacking nucleophiles. Lastly, calculations were performed on the active site of the Ca(2+)-ATPase E2P intermediate, which provide an example of how a protein environment facilitates the formation of a destabilized ground state.

The hydrolysis activity of adenosine triphosphate in myosin: a theoretical analysis of anomeric effects and the nature of the transition state

Combined quantum mechanical/molecular mechanical (QM/MM) calculations with density functional theory are employed to analyze two issues related to the hydrolysis activity of adenosine triphosphate (ATP) in myosin. First, we compare the geometrical properties and electronic structure of ATP in the open (post-rigor) and closed (pre-powerstroke) active sites of the myosin motor domain. Compared to both solution and the open active site cases, the scissile P(gamma)-O(3beta) bond of ATP in the closed active site is shown to be substantially elongated. Natural bond orbital (NBO) analysis clearly shows that this structural feature is correlated with the stronger anomeric effects in the closed active site, which involve charge transfers from the lone pairs in the nonbridging oxygen in the gamma-phosphate to the antibonding orbital of the scissile bond. However, an energetic analysis finds that the ATP molecule is not significantly destabilized by the P(gamma)-O(3beta) bond elongation. Therefore, despite the notable perturbations in the geometry and electronic structure of ATP as its environment changes from solution to the hydrolysis-competent active site, ground-state destabilization is unlikely to play a major role in enhancing the hydrolysis activity in myosin. Second, two-dimensional potential energy maps are used to better characterize the energetic landscape near the hydrolysis transition state. The results indicate that the transition-state region is energetically flat and a range of structures representative of different mechanisms according to the classical nomenclature (e.g., "associative", "dissociative", and "concerted") are very close in energy. Therefore, at least in the case of ATP hydrolysis in myosin, the energetic distinction between different reaction mechanisms following the conventional nomenclature is likely small. This study highlights the importance of (i) explicitly evaluating the relevant energetic properties for determining whether a factor is essential to catalysis and (ii) broader explorations of the energy landscape beyond saddle points (even on free-energy surface) for characterizing the molecular mechanism of catalysis.

Hydrolysis of aspartic acid phosphoramidate nucleotides: a comparative quantum chemical study

L-Aspartic acid has recently been found to be a good leaving group during HIV reverse transcriptase catalyzed incorporation of deoxyadenosine monophosphate (dAMP) in DNA. This showed that L-Asp is a good mimic for the pyrophosphate moiety of deoxyadenosine triphosphate. The present work explores the thermochemistry and mechanism for hydrolysis of several models for L-aspartic-dAMP using B3LYP/DGDZVP, MP2/6-311++G** and G3MP2 level of theory. The effect of the new compound is gradually investigated: starting from a simple methyl amine leaving group up to the aspartic acid leaving group. The enzymatic environment was mimicked by involving two Mg(2+) ions and some important active site residues in the reaction. All reactions are compared to the corresponding O-coupled leaving group, which is methanol for methyl amine and malic acid for aspartic acid. With methyl amine as a leaving group a tautomeric associative or tautomeric dissociative mechanism is preferred and the barrier is lower than the comparable mechanism with methanol as a leaving group. The calculations on the aspartic acid in the enzymatic environment show that qualitatively the mechanism is the same as for triphosphate but the barrier for hydrolysis by the associative mechanism is higher for L-aspartic-dAMP than for L-malic-dAMP and pyrophosphate.