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 the Mechanisms of Nonenzymatic Intramolecular Cyclization of the Glutamine Residues Located at N-Termini Catalyzed by Inorganic Phosphate Species

Glutamine (Gln) residues located at N-termini undergo spontaneous intramolecular cyclization, causing the formation of pyroglutamic acid (pGlu) residues. pGlu residues have been detected at the N-termini in various peptides and proteins. The formation of pGlu residues during the fermentation and purification processes of antibody drugs is one of the concerns in the design and formulation of these drugs and has been reported to proceed rapidly in a phosphate buffer. In this study, we have examined the phosphate-catalyzed mechanisms of the pGlu residue formation from N-terminal Gln residues via quantum chemical calculations using B3LYP density functional methods. Single-point energies were calculated using the second-order Møller-Plesset perturbation theory. We performed the calculations for the model compound in which an uncharged N-terminal Gln residue is capped with a methyl amino group on the C-terminal. The activation energy of the formation of pGlu residues was calculated as 83.8 kJ mol^-1, which was lower than that of the typical nonenzymatic reaction of amino acid residues. In addition, the computational results indicate that the flexibility of the main and side chains in N-terminal Gln residues was necessary for the formation of pGlu residues to proceed. In the obtained pathway, inorganic phosphate species act as the catalyst by mediating the proton transfer.

Computational Studies on Water-Catalyzed Mechanisms for Stereoinversion of Glutarimide Intermediates Formed from Glutamic Acid Residues in Aqueous Phase

Aspartic acid (Asp) residues are prone to non-enzymatic stereoinversion, and Asp-residue stereoinversion is believed to be mediated via a succinimide (SI) intermediate. The stereoinverted Asp residues are believed to cause several age-related diseases. However, in peptides and proteins, few studies have reported the stereoinversion of glutamic acid (Glu) residues whose structures are similar to that of Asp. We previously presumed that Glu-residue stereoinversion proceeds via a glutarimide (GI) intermediate and showed that the calculated activation barriers of SI- and GI-intermediate stereoinversion are almost equivalent in the gas phase. In this study, we investigated the stereoinversion pathways of the l-GI intermediate in the aqueous phase using B3LYP density functional methods. The calculated activation barrier of l-GI-intermediate stereoinversion in the aqueous phase was approximately 36 kcal·mol-1, which was much higher than that in the gas phase. Additionally, as this activation barrier exceeded that of Asp-residue stereoinversion, it is presumed that Glu-residue stereoinversion has a lower probability of proceeding under physiological conditions than Asp-residue stereoinversion.

Phosphate-Catalyzed Succinimide Formation from an NGR-Containing Cyclic Peptide: A Novel Mechanism for Deammoniation of the Tetrahedral Intermediate

Spontaneous deamidation in the Asn-Gly-Arg (NGR) motif that yields an isoAsp-Gly-Arg (isoDGR) sequence has recently attracted considerable attention because of the possibility of application to dual tumor targeting. It is well known that Asn deamidation reactions in peptide chains occur via the five-membered ring succinimide intermediate. Recently, we computationally showed by the B3LYP density functional theory method, that inorganic phosphate and the Arg side chain can catalyze the NGR deamidation using a cyclic peptide, c[CH₂CO⁻NGRC]⁻NH₂. In this previous study, the tetrahedral intermediate of the succinimide formation was assumed to be readily protonated at the nitrogen originating from the Asn side chain by the solvent water before the release of an NH₃ molecule. In the present study, we found a new mechanism for the decomposition of the tetrahedral intermediate that does not require the protonation by an external proton source. The computational method is the same as in the previous study. In the new mechanism, the release of an NH₃ molecule occurs after a proton exchange between the peptide and the phosphate and conformational changes. The rate-determining step of the overall reaction course is the previously reported first step, i.e., the cyclization to form the tetrahedral intermediate.

Influences of conformations of peptides on stereoinversions and/or isomerizations of aspartic acid residues

Recently, non-enzymatic stereoinversions of aspartic acid (Asp) residues in proteins and peptides have been reported. Here, we performed replica exchange molecular dynamics (REMD) simulations of model peptides (exon 6, 26A-1, and 26A-2) extracted from elastin to investigate their structural features, thereby revealing the factor that influences stereoinversions. For REMD trajectories, we calculated distances between carboxyl carbon in Asp and amide nitrogen in the (n + 1) residue (CN distances). Because bond formation between carbon and nitrogen is indispensable to the formation of a succinimide intermediate the distance between them seems to play an important role in stereoinversion. Moreover, we calculated polar surface areas (PSAs) for the trajectories, finding that CN distances and PSA were different for each peptide, with the longest CN distance and smallest PSA observed for exon 6 peptide, where stereoinversion of Asp is the slowest. Although the average CN distance was shorter for exon 26A-1 peptide than for exon 26A-2 peptide, the number of conformations with CN distances <3.0 Å was greater for exon 26A-2 peptide than for exon 26A-1 peptide. Furthermore, PSA for amide nitrogen of the (n + 1) residue was larger for exon 26A-2 peptide than for exon 26A-1 peptide. These results indicated that the flexibility of Asp and (n + 1) residues and hydrophilicity of peptides, especially in the (n + 1) residue, play important roles in the stereoinversion of Asp. This article is part of a Special Issue entitled: D-Amino acids: biology in the mirror, edited by Dr. Loredano Pollegioni, Dr. Jean-Pierre Mothet and Dr. Molla Gianluca.

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 Serine Residues Catalyzed by Dihydrogen Phosphate Ion: A Computational Study

Spontaneous, nonenzymatic reactions in proteins are known to have relevance to aging and age-related diseases, such as cataract and Alzheimer’s disease. Among such reactions is the racemization of Ser residues, but its mechanism in vivo remains to be clarified. The most likely intermediate is an enol. Although being nonenzymatic, the enolization would need to be catalyzed to occur at a biologically relevant rate. In the present study, we computationally found plausible reaction pathways for the enolization of a Ser residue where a dihydrogen phosphate ion, H2PO4−, acts as a catalyst. The H2PO4− ion mediates the proton transfer required for the enolization by acting simultaneously as both a general base and a general acid. Using the B3LYP density functional theory method, reaction pathways were located in the gas phase and hydration effects were evaluated by single-point calculations using the SM8 continuum model. The activation barriers calculated for the reaction pathways found were around 100 kJ mol−1, which is consistent with spontaneous reactions occurring at physiological temperature. Our results are also consistent with experimental observations that Ser residue racemization occurs more readily in flexible regions in proteins.

Identification of Sequence Similarities among Isomerization Hotspots in Crystallin Proteins

The eye lens crystallins represent an ideal target for studying the effects of aging on protein structure. Herein we examine separately the water-soluble (WS) and water-insoluble (WI) crystallin fractions and identify sites of isomerization and epimerization. Both collision-induced dissociation and radical-directed dissociation are needed for detection of these non-mass-shifting post-translational modifications. Isomerization levels differ significantly between the WS and the WI fractions from sheep, pig, and cow eye lenses. Residues that are most susceptible to isomerization are identified site-specifically and are found to reside in structurally disordered regions. However, isomerization in structured domains, although less common, often yields more dramatic effects on solubility. Numerous isomerization hotspots were also identified and occur in regions with aspartic acid and serine repeats. For example, ¹²⁸KMEIVDDDVPSLW¹⁴⁰ in βB3 crystallin contains three sequential aspartic acid residues and is isomerized heavily in the WI fractions, while it is not modified at all in the WS fractions. Potential causes for enhanced isomerization at sites with acidic residue repeats are presented. The importance of acidic residue repeats extends beyond the lens, as they are found in many other long-lived proteins associated with disease.

Succinimide Formation from an NGR-Containing Cyclic Peptide: Computational Evidence for Catalytic Roles of Phosphate Buffer and the Arginine Side Chain

The Asn-Gly-Arg (NGR) motif and its deamidation product isoAsp-Gly-Arg (isoDGR) have recently attracted considerable attention as tumor-targeting ligands. Because an NGR-containing peptide and the corresponding isoDGR-containing peptide target different receptors, the spontaneous NGR deamidation can be used in dual targeting strategies. It is well known that the Asn deamidation proceeds via a succinimide derivative. In the present study, we computationally investigated the mechanism of succinimide formation from a cyclic peptide, c[CH₂CO-NGRC]-NH₂, which has recently been shown to undergo rapid deamidation in a phosphate buffer. An H₂PO₄⁻ ion was explicitly included in the calculations. We employed the density functional theory using the B3LYP functional. While geometry optimizations were performed in the gas phase, hydration Gibbs energies were calculated by the SM8 (solvation model 8) continuum model. We have found a pathway leading to the five-membered ring tetrahedral intermediate in which both the H₂PO₄⁻ ion and the Arg side chain act as catalyst. This intermediate, once protonated at the NH₂ group on the five-membered ring, was shown to easily undergo NH₃ elimination leading to the succinimide formation. This study is the first to propose a possible catalytic role for the Arg side chain in the NGR deamidation.

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.

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.

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

Succinimide formation from aspartic acid (Asp) residues is a concern in the formulation of protein drugs. Based on density functional theory calculations using Ace-Asp-Nme (Ace = acetyl, Nme = NHMe) as a model compound, we propose the possibility that acetic acid (AA), which is often used in protein drug formulation for mildly acidic buffer solutions, catalyzes the succinimide formation from Asp residues by acting as a proton-transfer mediator. The proposed mechanism comprises two steps: cyclization (intramolecular addition) to form a gem-diol tetrahedral intermediate and dehydration of the intermediate. Both steps are catalyzed by an AA molecule, and the first step was predicted to be rate-determining. The cyclization results from a bond formation between the amide nitrogen on the C-terminal side and the side-chain carboxyl carbon, which is part of an extensive bond reorganization (formation and breaking of single bonds and the interchange of single and double bonds) occurring concertedly in a cyclic structure formed by the amide NH bond, the AA molecule and the side-chain C=O group and involving a double proton transfer. The second step also involves an AA-mediated bond reorganization. Carboxylic acids other than AA are also expected to catalyze the succinimide formation by a similar mechanism.

Roles of intramolecular and intermolecular hydrogen bonding in a three-water-assisted mechanism of succinimide formation from aspartic acid residues

Aspartic acid (Asp) residues in peptides and proteins are prone to isomerization to the β-form and racemization via a five-membered succinimide intermediate. These nonenzymatic reactions have relevance to aging and age-related diseases. In this paper, we report a three water molecule-assisted, six-step mechanism for the formation of succinimide from Asp residues found by density functional theory calculations. The first two steps constitute a stepwise iminolization of the C-terminal amide group. This iminolization involves a quintuple proton transfer along intramolecular and intermolecular hydrogen bonds formed by the C-terminal amide group, the side-chain carboxyl group, and the three water molecules. After a conformational change (which breaks the intramolecular hydrogen bond involving the iminol nitrogen) and a reorganization of water molecules, the iminol nitrogen nucleophilically attacks the carboxyl carbon of the Asp side chain to form a five-membered ring. This cyclization is accompanied by a triple proton transfer involving two water molecules, so that a gem-diol tetrahedral intermediate is formed. The last step is dehydration of the gem-diol group catalyzed by one water molecule, and this is the rate-determining step. The calculated overall activation barrier (26.7 kcal mol⁻¹) agrees well with an experimental activation energy.