Acetic acid can catalyze succinimide formation from aspartic acid residues by a concerted bond reorganization mechanism: a computational study

Predicting Reaction Mechanisms for the Threonine-Residue Stereoinversion Catalyzed by a Dihydrogen Phosphate Ion

The stereoinversion of amino acid residues in proteins is considered to trigger various age-related diseases. Serine (Ser) residues are relatively prone to stereoinversion. It is assumed that threonine (Thr) residues also undergo stereoinversion, which results in the formation of the d-allo-Thr residue, by the same mechanisms as those for Ser-residue stereoinversion; however, d-allo-Thr residues have not been detected in vivo. To date, although Ser-residue stereoinversion has been suggested to progress via enolization, plausible reaction mechanisms for Thr-residue stereoinversion have not been proposed. In this study, we investigated the pathway of Thr-residue enolization and successfully identified the three types of plausible reaction pathways of Thr-residue stereoinversion catalyzed by a dihydrogen phosphate ion. The geometries of reactant complexes, transition states, and enolized product complexes were optimized using B3LYP density functional methods, and single-point calculations were performed for all optimized geometries using Møller-Plesset perturbation theory to obtain reliable energies. As a result, the calculated activation energies of Thr-residue stereoinversion were 105-106 kJ mol^-1, which were comparable with those of Ser-residue stereoinversion reported previously. The infrequency of Thr-residue stereoinversion may be due to other factors, such as the hydrophobicity and/or the steric hindrance of the γ-methyl group, rather than the high activation energies.

Molecular Mechanisms of Succinimide Formation from Aspartic Acid Residues Catalyzed by Two Water Molecules in the Aqueous Phase

Aspartic acid (Asp) residues are prone to nonenzymatic isomerization via a succinimide (Suc) intermediate. The formation of isomerized Asp residues is considered to be associated with various age-related diseases, such as cataracts and Alzheimer’s disease. In the present paper, we describe the reaction pathway of Suc residue formation from Asp residues catalyzed by two water molecules using the B3LYP/6-31+G(d,p) level of theory. Single-point energies were calculated using the MP2/6-311+G(d,p) level of theory. For these calculations, we used a model compound in which an Asp residue was capped with acetyl and methylamino groups on the N- and C-termini, respectively. In the aqueous phase, Suc residue formation from an Asp residue was roughly divided into three steps, namely, iminolization, cyclization, and dehydration, with the activation energy estimated to be 109 kJ mol⁻¹. Some optimized geometries and reaction modes in the aqueous phase were observed that differed from those in the gas phase.

Influence of the conformations of αA-crystallin peptides on the isomerization rates of aspartic acid residues

The isomerization rate of aspartic acid (Asp) residue is known to be affected by the three-dimensional structures of peptides and proteins. Although the isomerized Asp residues were experimentally observed, structural features which affect the isomerization cannot be elucidated sufficiently because of protein denaturation and aggregation. In this study, molecular dynamics (MD) simulations were conducted on three αA-crystallin peptides (T6, T10, and T18), each containing a single Asp residue with different isomerization rate (T18 > T6 > T10) to clarify the structural factors of Asp isomerization tendency. For MD trajectories, distances between side-chain carboxyl carbon of Asp and main-chain amide nitrogen of (n + 1) residue (C_γ-N distances), root mean square fluctuations (RMSFs), and polar surface areas for main-chain amide nitrogen of (n + 1) residues (PSA_N) were calculated, because these structural features are considered to relate to the formations of cyclic imide intermediates. RMSFs and PSA_N are indexes of peptide backbone flexibilities and solvent exposure of the amide nitrogen, respectively. The average C_γ-N distances of T10 was longer than those of the other two peptides. In addition, the peptide containing Asp residue with a higher isomerization rate showed higher flexibility of the peptide backbone around the Asp residue. PSA_N for amide nitrogen in T18 were much larger than those of other two peptides. The computational results suggest that Asp-residue isomerization rates are affected by these factors.

Computational studies on nonenzymatic succinimide-formation mechanisms of the aspartic acid residues catalyzed by two water molecules

In the biological proteins, aspartic acid (Asp) residues are prone to nonenzymatic isomerization via a succinimide (Suc) intermediate. Asp-residue isomerization causes the aggregation and the insolubilization of proteins, and is considered to be involved in various age-related diseases. Although Suc intermediate was considered to be formed by nucleophilic attack of the main-chain amide nitrogen of N-terminal side adjacent residue to the side-chain carboxyl carbon of Asp residue, previous studies have shown that the nucleophilic attack is more likely to proceed via iminol tautomer when the water molecules act as catalysts. However, the full pathway to Suc-intermediate formation has not been investigated, and the experimental analyses for the Asp-residue isomerization mechanism at atomic and molecular levels, such as the analysis of the transition state geometry, are difficult. In the present study, we computationally explored the full pathways for Suc-intermediate formation from Asp residues. The calculations were performed two types of reactant complexes, and all energy minima and TS geometries were optimized using B3LYP density functional methods. As a result, the SI-intermediate formation was divided into three processes, i.e., iminolization, cyclization, and dehydration processes, and the activation energies were calculated to be 26.1 or 28.4 kcal mol^-1. These values reproduce the experimental data. The computational results show that abundant water molecules in living organisms are effective catalysts for the Asp-residue isomerization.

Possible Mechanisms of Nonenzymatic Formation of Dehydroalanine Residue Catalyzed by Dihydrogen Phosphate Ion

Uncommon crosslinked amino acids have been identified in several aging tissues and are suspected to trigger various age-related diseases. Several uncommon residues are formed when the dehydroalanine (Dha) residue undergoes a nucleophilic attack by surrounding residues. Dha residues are considered to be formed by posttranslational modification of serine (Ser) and cysteine residues. In the present study, we investigated the Dha residue formation mechanism catalyzed by dihydrogen phosphate ion (H₂PO₄^-) using quantum chemical calculations. We obtained optimized geometries using the B3LYP density functional method and carried out single-point energy calculations using the second-order Møller-Plesset perturbation method. All calculations were performed using Ace-Ser-Nme (Ace = acetyl, Nme = methylamino) as a model compound. Results of the computational analysis suggest that the mechanism underlying the Dha residue formation from Ser consists of two steps: enolization and 1,3-elimination. The H₂PO₄^- catalyzed both reactions as a proton-relay mediator. The calculated activation barrier for Dha residue formation was estimated as 30.4 kcal mol^-1. In this pathway, the catalytic H₂PO₄^- interacts with the Ser residue α-proton, carbonyl oxygen of Ser, and C-terminal side adjacent residues, and the calculated activation energy produced was the same as the experimentally reported value for nonenzymatic modifications of amino acid residues. Therefore, our calculation suggests that H₂PO₄^--catalyzed Ser residue dehydration can proceed nonenzymatically.

Phosphate-Catalyzed Succinimide Formation from Asp Residues: A Computational Study of the Mechanism

Aspartic acid (Asp) residues in proteins and peptides are prone to the non-enzymatic reactions that give biologically uncommon l-β-Asp, d-Asp, and d-β-Asp residues via the cyclic succinimide intermediate (aminosuccinyl residue, Suc). These abnormal Asp residues are known to have relevance to aging and pathologies. Despite being non-enzymatic, the Suc formation is thought to require a catalyst under physiological conditions. In this study, we computationally investigated the mechanism of the Suc formation from Asp residues that were catalyzed by the dihydrogen phosphate ion, H₂PO₄⁻. We used Ac–l-Asp–NHMe (Ac = acetyl, NHMe = methylamino) as a model compound. The H₂PO₄⁻ ion (as a catalyst) and two explicit water molecules (as solvent molecules stabilizing the negative charge) were included in the calculations. All of the calculations were performed by density functional theory with the B3LYP functional. We revealed a phosphate-catalyzed two-step mechanism (cyclization–dehydration) of the Suc formation, where the first step is predicted to be rate-determining. In both steps, the reaction involved a proton relay mediated by the H₂PO₄⁻ ion. The calculated activation barrier for this mechanism (100.3 kJ mol⁻¹) is in reasonable agreement with an experimental activation energy (107 kJ mol⁻¹) for the Suc formation from an Asp-containing peptide in a phosphate buffer, supporting the catalytic mechanism of the H₂PO₄⁻ ion that is revealed in this study.

Comparison of the activation energy barrier for succinimide formation from α- and β-aspartic acid residues obtained from density functional theory calculations

The l-α-Asp residues in peptides or proteins are prone to undergo nonenzymatic reactions to form l-β-Asp, d-α-Asp, and d-β-Asp residues via a succinimide five-membered ring intermediate. From these three types of isomerized aspartic acid residues, particularly d-β-Asp has been widely detected in aging tissue. In this study, we computationally investigated the cyclization of α- and β-Asp residues to form succinimide with dihydrogen phosphate ion as a catalyst (H₂PO₄^-). We performed the study using B3LYP/6-31+G(d,p) density functional theory calculations. The comparison of the activation barriers of both residues is discussed. All the calculations were performed using model compounds in which an α/β-Asp-Gly sequence is capped with acetyl and methylamino groups on the N- and C-termini, respectively. Moreover, H₂PO₄^- catalyzes all the steps of the succinimide formation (cyclization-dehydration) acting as a proton-relay mediator. The calculated activation energy barriers for succinimide formation of α- and β-Asp residues are 26.9 and 26.0kcalmol^-1, respectively. Although it was experimentally confirmed that β-Asp has higher stability than α-Asp, there was no clear difference between the activation barriers. Therefore, the higher stability of β-Asp residue than α-Asp residue may be caused by an entropic effect associated with the succinimide formation.

Racemization of the Succinimide Intermediate Formed in Proteins and Peptides: A Computational Study of the Mechanism Catalyzed by Dihydrogen Phosphate Ion

In proteins and peptides, d-aspartic acid (d-Asp) and d-β-Asp residues can be spontaneously formed via racemization of the succinimide intermediate formed from l-Asp and l-asparagine (l-Asn) residues. These biologically uncommon amino acid residues are known to have relevance to aging and pathologies. Although nonenzymatic, the succinimide racemization will not occur without a catalyst at room or biological temperature. In the present study, we computationally investigated the mechanism of succinimide racemization catalyzed by dihydrogen phosphate ion, H₂PO₄⁻, by B3LYP/6-31+G(d,p) density functional theory calculations, using a model compound in which an aminosuccinyl (Asu) residue is capped with acetyl (Ace) and NCH₃ (Nme) groups on the N- and C-termini, respectively (Ace–Asu–Nme). It was shown that an H₂PO₄⁻ ion can catalyze the enolization of the H_α–C_α–C=O portion of the Asu residue by acting as a proton-transfer mediator. The resulting complex between the enol form and H₂PO₄⁻ corresponds to a very flat intermediate region on the potential energy surface lying between the initial reactant complex and its mirror-image geometry. The calculated activation barrier (18.8 kcal·mol⁻¹ after corrections for the zero-point energy and the Gibbs energy of hydration) for the enolization was consistent with the experimental activation energies of Asp racemization.

A Computational Study of the Mechanism of Succinimide Formation in the Asn-His Sequence: Intramolecular Catalysis by the His Side Chain

The rates of deamidation reactions of asparagine (Asn) residues which occur spontaneously and nonenzymatically in peptides and proteins via the succinimide intermediate are known to be strongly dependent on the nature of the following residue on the carboxyl side (Xxx). The formation of the succinimide intermediate is by far the fastest when Xxx is glycine (Gly), the smallest amino acid residue, while extremely slow when Xxx is bulky such as isoleucine (Ile) and valine (Val). In this respect, it is very interesting to note that the succinimide formation is definitely accelerated when Xxx is histidine (His) despite its large size. In this paper, we computationally show that, in an Asn-His sequence, the His side-chain imidazole group (in the neutral Nε-protonated form) can specifically catalyze the formation of the tetrahedral intermediate in the succinimide formation by mediating a proton transfer. The calculations were performed for Ace-Asn-His-Nme (Ace = acetyl, Nme = methylamino) as a model compound by the density functional theory with the B3LYP functional and the 6-31+G(d,p) basis set. We also show that the tetrahedral intermediate, once protonated at the NH₂ group, easily releases an ammonia molecule to give the succinimide species.

Acetic acid-catalyzed formation of N-phenylphthalimide from phthalanilic acid: a computational study of the mechanism

In glacial acetic acid, phthalanilic acid and its monosubstituents are known to be converted to the corresponding phthalimides in relatively good yields. In this study, we computationally investigated the experimentally proposed two-step (addition-elimination or cyclization-dehydration) mechanism at the second-order Møller-Plesset perturbation (MP2) level of theory for the unsubstituted phthalanilic acid, with an explicit acetic acid molecule included in the calculations. In the first step, a gem-diol tetrahedral intermediate is formed by the nucleophilic attack of the amide nitrogen. The second step is dehydration of the intermediate to give N-phenylphthalimide. In agreement with experimental findings, the second step has been shown to be rate-determining. Most importantly, both of the steps are catalyzed by an acetic acid molecule, which acts both as proton donor and acceptor. The present findings, along with those from our previous studies, suggest that acetic acid and other carboxylic acids (in their undissociated forms) can catalyze intramolecular nucleophilic attacks by amide nitrogens and breakdown of the resulting tetrahedral intermediates, acting simultaneously as proton donor and acceptor. In other words, double proton transfers involving a carboxylic acid molecule can be part of an extensive bond reorganization process from cyclic hydrogen-bonded complexes.

Glycolic acid-catalyzed deamidation of asparagine residues in degrading PLGA matrices: a computational study

Poly(lactic-co-glycolic acid) (PLGA) is a strong candidate for being a drug carrier in drug delivery systems because of its biocompatibility and biodegradability. However, in degrading PLGA matrices, the encapsulated peptide and protein drugs can undergo various degradation reactions, including deamidation at asparagine (Asn) residues to give a succinimide species, which may affect their potency and/or safety. Here, we show computationally that glycolic acid (GA) in its undissociated form, which can exist in high concentration in degrading PLGA matrices, can catalyze the succinimide formation from Asn residues by acting as a proton-transfer mediator. A two-step mechanism was studied by quantum-chemical calculations using Ace-Asn-Nme (Ace = acetyl, Nme = NHCH₃) as a model compound. The first step is cyclization (intramolecular addition) to form a tetrahedral intermediate, and the second step is elimination of ammonia from the intermediate. Both steps involve an extensive bond reorganization mediated by a GA molecule, and the first step was predicted to be rate-determining. The present findings are expected to be useful in the design of more effective and safe PLGA devices.