Attack Conformation Research Articles

Applications of the Near Attack Conformation (NAC) approach in the search for Acetylcholinesterase reactivators

The Near Attack Conformation (NAC) approach states that the efficiency of an enzyme-catalyzed reaction depends on the prior attainment of optimal conditions for substrate atom organization and positioning for bond formation. These conditions are prerequisites for the transition state (TS) in which the involved atoms are within the van der Waals range of contact and positioned at an angle similar to that achieved after bond formation. The successful application of this approach to investigate the reactivation mechanism of acetylcholinesterase inhibited by nerve agents has contributed to a better understanding of this mechanism and demonstrated consistent corroboration with experimental data. In this article, we summarize the accomplishments achieved thus far and outline future perspectives.

Virtual screening of carboxylic acid reductases for biocatalytic synthesis of 6-aminocaproic acid and 1,6-hexamethylenediamine.

The key precursors for nylon synthesis, that is, 6-aminocaproic acid (6-ACA) and 1,6-hexamethylenediamine (HMD), are produced from petroleum-based feedstocks. A sustainable biocatalytic alternative method from bio-based adipic acid has been demonstrated recently. However, the low efficiency and specificity of carboxylic acid reductases (CARs) used in the process hampers its further application. Herein, we describe a highly accurate protein structure prediction-based virtual screening method for the discovery of new CARs, which relies on near attack conformation frequency and the Rosetta Energy Score. Through virtual screening and functional detection, five new CARs were selected, each with a broad substrate scope and the highest activities toward various di- and ω-aminated carboxylic acids. Compared with the reported CARs, KiCAR was highly specific with regard to adipic acid without detectable activity to 6-ACA, indicating a potential for 6-ACA biosynthesis. In addition, MabCAR3 had a lower Km with regard to 6-ACA than the previously validated CAR MAB4714, resulting in twice conversion in the enzymatic cascade synthesis of HMD. The present work highlights the use of structure-based virtual screening for the rapid discovery of pertinent new biocatalysts.

Structural and Energetic Origin of Different Product Specificities and Activities for SETD3 and Its Mutants on the Methylation of the β-Actin H73K Peptide: Insights from a QM/MM Study.

The methylation of the lysine residue can affect some fundamental biological processes, and specific biological effects of the methylations are often related to product specificity of methyltransferases. The question remains concerning how active-site structural features and dynamics control the activity as well as the number (1, 2, or 3) of methyl groups on methyl lysine products. SET domain containing protein 3 (SETD3) has been identified recently as the β-actin histidine73-N3 methyltransferase, and also, it has a weak methylation activity on the H73K β-actin peptide for which the target H73 residue is mutated into K73. Interestingly, the K73 methylation activity of SETD3 increases significantly as a result of the N255 → A or N255 → F/W273 → A mutation, and the N255A product specificity also differs from that of wild-type. Here, we performed QM/MM molecular dynamics and potential of mean force (PMF) simulations for SETD3 and its mutants (N255A and N255F/W273A) to study how SETD3 and its mutants could have different product specificities and activities for the K73 methylation. The PMF simulations show that the barrier for the first methylation of K73 is higher compared to the barrier of the H73 methylation in SETD3. Moreover, the second methylation of K73 has been found to have a barrier from the free energy simulation that is higher by 2.2 kcal/mol compared to the barrier of the first methyl transfer to K73, agreeing with the suggestion that SETD3 is a monomethylase. For the first, second, and third methylations of K73 in the N255A mutant, the barriers obtained from the PMF simulations for transferring the second and third methyl groups are found to be lower relative to the barrier for the first methyl transfer. Thus, N255A can be considered as a trimethyl lysine methyltransferase. In addition, for the first K73 methylation, the activities from the PMF simulations follow the order of N255F/W273A > N255A > WT, in agreement with experiments. The examination of the structural and dynamic results at the active sites provides better understanding of different product specificities and activities for the K73 methylations in SETD3 and its mutants. It is demonstrated that the existence of well-balanced interactions at the active site leading to the near attack conformation is of crucial importance for the efficient methyl transfers. Moreover, the presence of potential interactions (e.g., the C-H···O and cation-π interactions) that are strengthening at the transition state can also be important. Furthermore, the activity as well as product specificity of the K73 methylation also seems to be controlled by certain active-site water molecules which may be released to provide extra space for the addition of more methyl groups on K73.

Molecular Mechanism of Protein Arginine Deiminase 2: A Study Involving Multiple Microsecond Long Molecular Dynamics Simulations.

Peptidylarginine deiminase 2 (PAD2) is a Ca2+-dependent enzyme that catalyzes the conversion of protein arginine residues to citrulline. This kind of structural modification in histone molecules may affect gene regulation, leading to effects that may trigger several diseases, including breast cancer, which makes PAD2 an attractive target for anticancer drug development. To design new effective inhibitors to control activation of PAD2, improving our understanding of the molecular mechanisms of PAD2 using up-to-date computational techniques is essential. We have designed five different PAD2–substrate complex systems based on varying protonation states of the active site residues. To search the conformational space broadly, multiple independent molecular dynamics simulations of the complexes have been performed. In total, 50 replica simulations have been performed, each of 1 μs, yielding a total simulation time of 50 μs. Our findings identify that the protonation states of Cys647, Asp473, and His471 are critical for the binding and localization of the N-α-benzoyl-l-arginine ethyl ester substrate within the active site. A novel mechanism for enzyme activation is proposed according to near attack conformers. This represents an important step in understanding the mechanism of citrullination and developing PAD2-inhibiting drugs for the treatment of breast cancer.

Extensive modelling and quantum chemical study of sterol C-22 desaturase mechanism: A commercially important cytochrome P450 family

CYP450s, first wave biocatalysts and nature’s tool used for catalysing complex reactions like C-H bond hydroxylation, N-dealkylation, N-hydroxylation, O-dealkylation, etc. Engineered CYP450s catalyse new-to-nature carbene and nitrene transfer some of which have no precedent using chemical methods. Due to broad substrate range and immense biocatalytic scope, CYP450s are extensively used in synthesising of active pharmaceutical ingredients (APIs) and drugs. CYP450 desaturase adds double bonds to form mono and poly desaturated sterols. The addition of a double bond at a specific position is conserved within a particular structural fold type of CYP450. Minimal information is currently available on the mechanism of enzyme action of CYP450 desaturase, especially no information is available concerning experimental details on C-22 desaturase. Here, for the first time, we report the structure-function of desaturase and the detailed mechanism of the enzyme using quantum chemical (QC) calculations. The study provides atomistic details from precise near attack conformations of different sterols in the active site of CYP450 desaturase to the active site residues that define substrate specificity. A novel finding on second hydrogen abstraction; QC studies revealed transition states involving hydrogen abstractions. The second hydrogen abstraction takes place in a relay mode moving from carbon to carbon and then to electronegative oxygen attached to haem. Sterol diffusion path was predicted, and its plasticity was studied using extended molecular dynamic simulation experiments. This report is also a source of information for N-terminal modification and its effect on heterologous expression and the role of reductase binding for electron transfer. The results presented here will be valuable for guiding the experimental work to further improve the biocatalytic applicability of CYP450 desaturase enzymes.

Effect of the acyl-group length on the chemoselectivity of the lipase-catalyzed acylation of propranolol-a computational study.

The selective N-acylation of 1,2-amino alcohols has been proposed to occur through the proton shuttle mechanism. However, the O-acetylation of propranolol catalyzed by Candida antarctica lipase B is an exception. We investigated the relation between the chemoselectivity of this reaction and the acyl group length. For this purpose, we compared the acyl groups: ethanoyl, butanoyl, octanoyl, and hexadecanoyl. We studied the Michaelis complexes between serine-acylated Candida antarctica lipase B and propranolol, employing a computational approach that involved sampling Michaelis complex conformations through ensemble docking plus consensus scoring and molecular dynamics simulations. The conformations were then classified as near attack conformations for acylation of the amino or hydroxy group. The relative populations of these two classes of conformations were found to be consistent with the experimentally observed chemoselective O-acetylation. We predict that increasing the length of the hydrocarbon chain of the acyl group will cause O-acylation to be unfavorable with respect to N-acylation. The nucleophilic attack of propranolol to the acylated lipase was found to be more favorable through the classical mechanism when compared with the proton shuttle mechanism.

Strict conformational demands of RNA cleavage in bulge-loops created by peptidyl-oligonucleotide conjugates.

Potent knockdown of pathogenic RNA in vivo is an urgent health need unmet by both small-molecule and biologic drugs. ‘Smart’ supramolecular assembly of catalysts offers precise recognition and potent destruction of targeted RNA, hitherto not found in nature. Peptidyl-oligonucleotide ribonucleases are here chemically engineered to create and attack bulge-loop regions upon hybridization to target RNA. Catalytic peptide was incorporated either via a centrally modified nucleotide (Type 1) or through an abasic sugar residue (Type 2) within the RNA-recognition motif to reveal striking differences in biological performance and strict structural demands of ribonuclease activity. None of the Type 1 conjugates were catalytically active, whereas all Type 2 conjugates cleaved RNA target in a sequence-specific manner, with up to 90% cleavage from 5-nt bulge-loops (BC5-α and BC5L-β anomers) through multiple cuts, including in folds nearby. Molecular dynamics simulations provided structural explanation of accessibility of the RNA cleavage sites to the peptide with adoption of an ‘in-line’ attack conformation for catalysis. Hybridization assays and enzymatic probing with RNases illuminated how RNA binding specificity and dissociation after cleavage can be balanced to permit turnover of the catalytic reaction. This is an essential requirement for inactivation of multiple copies of disease-associated RNA and therapeutic efficacy.

Near Attack Conformation as Strategy for ESIPT Modulation for White-Light Generation

Two derivatives salicylidene-4,6-(dimethylamino)pyridine (KM1) and salicylidene-5-bromo-4,6-(dimethylamino)pyridine (KM2) reactive to excited-state intramolecular proton transfer (ESIPT) were synth...

Catalytic Mechanism and Product Specificity of Protein Arginine Methyltransferase PRMT7: A Study from QM/MM Molecular Dynamics and Free Energy Simulations.

QM/MM molecular dynamics and potential of mean force (PMF) free-energy simulations are performed for wild-type PRMT7 and E172Q, E181Q, and Q329A mutants in this work, and the catalytic mechanism, product specificity, and the role of key residues for the PRMT7 activity are investigated. The main strategies of PRMT7 in reducing the activation barrier for methyl transfer that are found in this study include (1) formation of reactive (near attack) conformations for the substrate Arg, (2) strengthening the active-site interactions at the transition state, and (3) generation of more effective nucleophiles by changing charge distributions on the target Arg through active-site interactions. More importantly, it is shown that it is a combination of these different factors that determines the PRMT7 methylation activity and substrate/product specificity. By taking these factors into consideration, it is possible to provide explanations for the observed effects of some mutations. For E172Q, E181Q, and Q329A, the simulation results suggest that E172Q has the least activity among the three mutants. The free energy barrier increases by 7 and 3 kcal/mol, respectively, as a result of the E181 → Q and Q329 → A mutations. The results showed that PRMT7 has a preference of adding a methyl group to the ω-guanidino nitrogen Nη2 atom of the substrate Arg and that the second methylation reactions cannot occur, which are consistent with previous investigations.

Fine-Tuning of Sequence Specificity by Near Attack Conformations in Enzyme-Catalyzed Peptide Hydrolysis

The catalytic role of near attack conformations (NACs), molecular states that lie on the pathway between the ground state (GS) and transition state (TS) of a chemical reaction, is not understood completely. Using a computational approach that combines Bürgi–Dunitz theory with all-atom molecular dynamics simulations, the role of NACs in catalyzing the first stages of HIV-1 protease peptide hydrolysis was previously investigated using a substrate that represents the recognized SP1-NC cleavage site of the HIV-1 Gag polyprotein. NACs were found to confer no catalytic effect over the uncatalyzed reaction there ( Δ Δ G N ‡ ∼ 0 kcal/mol). Here, using the same approach, the role of NACs across multiple substrates that each represent a further recognized cleavage site is investigated. Overall rate enhancement varies by | Δ Δ G ‡ | ∼ 12–15 kcal/mol across this set, and although NACs contribute a small and approximately constant barrier to the uncatalyzed reaction (< Δ G N ‡ u > = 4.3 ± 0.3 kcal/mol), they are found to contribute little significant catalytic effect ( | Δ Δ G N ‡ | ∼ 0–2 kcal/mol). Furthermore, no correlation is exhibited between NAC contributions and the overall energy barrier ( R 2 = 0.01). However, these small differences in catalyzed NAC contributions enable rates to match those required for the kinetic order of processing. Therefore, NACs may offer an alternative and subtle mode compared to non-NAC contributions for fine-tuning reaction rates during complex evolutionary sequence selection processes—in this case across cleavable polyproteins whose constituents exhibit multiple functions during the virus life-cycle.

Lipase mediated enzymatic kinetic resolution of phenylethyl halohydrins acetates: A case of study and rationalization

Racemic phenylethyl halohydrins acetates containing several groups attached to the aromatic ring were resolved via hydrolysis reaction in the presence of lipase B from Candida antarctica (Novozym® 435). In all cases, the kinetic resolution was highly selective (E > 200) leading to the corresponding (S)-β-halohydrin with ee > 99 %. However, the time required for an ideal 50 % conversion ranged from 15 min for 2,4-dichlorophenyl chlorohydrin acetate to 216 h for 2-chlorophenyl bromohydrin acetate. Six chlorohydrins and five bromohydrins were evaluated, the latter being less reactive. For the β-brominated substrates, steric hindrance on the aromatic ring played a crucial role, which was not observed for the β-chlorinated derivatives. To shed light on the different reaction rates, docking studies were carried out with all the substrates using MD simulations. The computational data obtained for the β-brominated substrates, based on the parameters analysed such as NAC (near attack conformation), distance between Ser-O and carbonyl-C and oxyanion site stabilization were in agreement with the experimental results. On the other hand, the data obtained for β-chlorinated substrates suggested that physical aspects such as high hydrophobicity or induced change in the conformation of the enzymatic active site are more relevant aspects when compared to steric hindrance effects.

Reactivation of VX-Inhibited Human Acetylcholinesterase by Deprotonated Pralidoxime. A Complementary Quantum Mechanical Study.

In the present work, we performed a complementary quantum mechanical (QM) study to describe the mechanism by which deprotonated pralidoxime (2-PAM) could reactivate human (Homo sapiens sapiens) acetylcholinesterase (HssAChE) inhibited by the nerve agent VX. Such a reaction is proposed to occur in subsequent addition–elimination steps, starting with a nucleophile bimolecular substitution (SN2) mechanism through the formation of a trigonal bipyramidal transition state (TS). A near attack conformation (NAC), obtained in a former study using molecular mechanics (MM) calculations, was taken as a starting point for this project, where we described the possible formation of the TS. Together, this combined QM/MM study on AChE reactivation shows the feasibility of the reactivation occurring via attack of the deprotonated form of 2-PAM against the Ser203-VX adduct of HssAChE.

Rational modulation of the enzymatic intermediates for tuning the phosphatase activity of histidine kinase HK853

Histidine kinase (HK) of two-component signal transduction system (TCS) is a potential drug target for treating bacterial infections, and most HKs are bifunctional. We have previously identified the HXXXT motif of HK in HisKA subfamily to perform the phosphatase activity, but the specific working mechanism of the threonine is not well understood. In this paper, we use the phosphate group analog BeF3− to capture the enzymatic intermediates between HK853 and RR468 from Thermotoga maritima during dephosphorylation, and demonstrate that the T264 site is essential for populating capable near attack conformers (NAC) between enzyme and substrate to facilitate catalysis. Importantly, mutations at this site can modulate the phosphatase activity of HK. Our results help to understand the TCS signal transduction mechanisms and provide a reference for drug design.

Molecular modelling studies on the interactions of 7-methoxytacrine-4-pyridinealdoxime with VX-inhibited human acetylcholinesterase. A near attack approach to assess different spacer-lengths

The novel prophylactic agent 7-methoxytacrine-4-pyridinealdoxime is a hybrid compound formerly designed to keep acetylcholinesterase resistant to organophosphates by reactivating it in case of intoxication by such inhibitors. In rational design, a 5-carbon length-spacer hybrid compound was synthesized to evaluate its inhibitory and reactivation capabilities. In this work, theoretical results were achieved through molecular modelling techniques, taking for granted the enzymatic reactivation reaction through nucleophilic substitution. Based on the near attack conformation approach, docking studies were performed to assess the spacer-length from 1 to 10 carbons long of a series of analogues of 7-methoxytacrine-4-pyridinealdoxime. Consequently, the hybrids with length-spacer of 4 and 5 carbons long were the best assessed and subsequently subjected to further molecular dynamics simulations, complemented by Poisson-Boltzmann surface area calculations. As a result, intermolecular interactions with the different binding sites inside human acetylcholinesterase were elucidated. Besides, thermodynamics and kinetics concepts pointed to the 4-carbon linker as optimum for enzymatic reactivation. Further studies, based on quantum mechanics in conjunction with molecular mechanics, were recommended to the presented near attack conformations to achieve more thermodynamics results between the hybrids with 4- and 5-carbon linkers, like values of activation energy for the reactivation reaction. All of those in silico evaluations could be considered as a set of tools for theoretically investigate novel enzymatic reactivators with different shape of spacers.

Connecting the dots in the mechanism of action of Cucurbitacin E (CurE) – path analysis and steered molecular dynamics reveal the precise site of entry and the passage of CurE in filamentous actin

Cucurbitacin E (CurE) modulates actin cytoskeleton by forming an irreversible covalent bond with Cys257 of actin. The reported binding conformation of CurE is deeply buried in the subdomain 4 of actin and is closely situated to the ATP-binding site. The entry and the path taken by CurE to reach this buried site remain a mystery. In this study, steered molecular dynamics (SMD) simulations were conducted to delineate the diffusion of CurE to its binding site. SMD simulations reveal that the distinctive entry site of CurE found in subdomain 4 is by itself a closed and compact region of two loops lying beside each other like a closed door and CurE induces it to open. From this point, CurE moves toward its binding site through a path facilitated by Thr188, Leu261, Ile267, Ile309, Tyr306. This study is also an insight into how CurE distinctly differentiates its tunnel to the binding site from the ATP-binding site. The conformational changes of CurE along the path to the binding site are surprisingly very minimal and closely resemble the attack conformation at the end of simulation. The study reveals that the little energy spent by the molecule is compensated by the enthalpic contribution to binding-free energy barrier making it undoubtedly the most preferred path of CurE. This study is the first of its kind in which the SMD was used to derive the complete and continuous translocation of one of the most potent phytochemicals, CurE through the binding site gorge of actin. Communicated by Ramaswamy H. Sarma

Molecular modeling studies on the interactions of 7-methoxytacrine-4-pyridinealdoxime, 4-PA, 2-PAM, and obidoxime with VX-inhibited human acetylcholinesterase: a near attack conformation approach

7-methoxytacrine-4-pyridinealdoxime (7-MEOTA-4-PA, named hybrid 5C) is a compound formerly synthesized and evaluated in vitro, together with 4-pyridine aldoxime (4-PA) and commercial reactivators of acetylcholinesterase (AChE). This compound was designed with the purpose of being a prophylactic reactivator, capable of interacting with different subdomains of the active site of AChE. To investigate these interactions, theoretical results from docking were first compared with experimental data of hybrid 5C, 4-PA, and two commercial oximes, on the reactivation of human AChE (HssAChE) inhibited by VX. Then, further docking studies, molecular dynamics simulations, and molecular mechanics Poisson–Boltzmann surface area calculations, were carried out to investigate reactivation performances, considering the near attack conformation (NAC) approach, prior to the nucleophilic substitution mechanism. Our results helped to elucidate the interactions of such molecules with the different subdomains of the active site of HssAChE. Additionally, NAC poses of each oxime were suggested for further theoretical studies on the reactivation reaction.

Relative Binding Affinities of Monolignols to Horseradish Peroxidase.

Monolignol binding to the peroxidase active site is the first step in lignin polymerization in plant cell walls. Using molecular dynamics, docking, and free energy perturbation calculations, we investigate the binding of monolignols to horseradish peroxidase C. Our results suggest that p-coumaryl alcohol has the strongest binding affinity followed by sinapyl and coniferyl alcohol. Stacking interactions between the monolignol aromatic rings and nearby phenylalanine residues play an important role in determining the calculated relative binding affinities. p-Coumaryl and coniferyl alcohols bind in a pose productive for reaction in which a direct H-bond is formed between the phenolic -OH group and a water molecule (W2) that may facilitate proton transfer during oxidation. In contrast, in the case of sinapyl alcohol there is no such direct interaction, the phenolic -OH group instead interacting with Pro139. Since proton and electron transfer is the rate-limiting step in monolignol oxidation by peroxidase, the binding pose (and thus the formation of near attack conformation) appears to play a more important role than the overall binding affinity in determining the oxidation rate.

Structure-Based Prediction of Protease Multispecificity using Computational Protein Design

The substrate specificity of protease enzymes constitutes the molecular basis of their diverse and complex biological roles. Rapid and accurate prediction of extended substrate specificity would also enable the design of custom enzymes capable of selectively cleaving biotechnologically or therapeutically relevant proteins. However, current in silico approaches for protease specificity prediction rely on, and are therefore limited by, machine learning of sequence patterns in known experimental data. Here, we describe a general structure-based approach for predicting peptidase substrates. We develop a structure-guided sequence sampling technique using self-consistent mean field theory for rapid specificity profile prediction and implement it in the Rosetta software. We show that computationally calculated specificity profiles closely match experimentally determined ones. Further, we construct atomic resolution models of thousands of substrate-enzyme complexes for each of six model proteases belonging to the four major protease mechanistic classes, and develop a discriminatory scoring function using enzyme design modules from Rosetta and Amber-MMPBSA. We rank putative peptide substrates in order of their interaction energies with a modeled near attack conformation of the enzyme active site. We show that the energetic patterns obtained from these simulations can be used to robustly classify known cleaved and uncleaved peptides, and these patterns have greater discriminatory power compared to purely sequence-based statistical inference. Combining sequence and energetic patterns for a given protease using machine-learning algorithms further improves classification performance, and analysis of atomic resolution models affords physical insight into the structural basis for the observed specificities. We tested the predictive capability of the approach by designing and experimentally characterizing the cleavage of four previously unidentified, but less rapidly cleaved, novel substrates for the Hepatitis C virus NS3/4 protease using an in vivo assay.

Peptide Bond Cleavage through Asparagine Cyclization

Protein splicing is the post translational modification by which an internal peptide, or intein, separates itself from the surrounding residues, or exteins. In the third step of the mechanism, asparagine cyclization, coupled to peptide bond cleavage, cleaves the intein from the exteins. We have studied this cyclization both experimentally and computationally. In the lab, we use a model peptide assay in which the cyclization of Asn is coupled to fluorescence. Cleavage via Asn cyclization is much faster than via Gln cyclization. High temperature, low dielectric constant, and high pH promote cyclization. We also have observed cleavage after Asp, but only at low pH. Computationally, we used molecular dynamics simulations, which suggest that the Asn side chain more frequently samples optimal attack conformations than a Gln side chain. Density functional theory and nudged elastic band simulations confirm the side chain nitrogen, not the oxygen, is the nucleophile and that nitrogen deprotonation is likely necessary for the attack to occur. The activation barriers of both the deprotonation and cyclization are being calculated to determine the rate limiting step for the mechanism. Simulations of peptides containing Asp suggest that cyclization occurs via a different mechanism than that of the native Asn and Gln residues. This material is based upon work supported by the National Science Foundation under grant MCB-1244089 and by the Camille and Henry Dreyfus Foundation.

Why does Asn71 deamidate faster than Asn15 in the enzyme triosephosphate isomerase? Answers from microsecond molecular dynamics simulation and QM/MM free energy calculations.

Deamidation is the uncatalyzed process by which asparagine or glutamine can be transformed into aspartic acid or glutamic acid, respectively. In its active homodimeric form, mammalian triosephosphate isomerase (TPI) contains two deamidation sites per monomer. Experimental evidence shows that the primary deamidation site (Asn71-Gly72) deamidates faster than the secondary deamidation site (Asn15-Gly16). To evaluate the factors controlling the rates of these two deamidation sites in TPI, we have performed graphics processing unit-enabled microsecond long molecular dynamics simulations of rabbit TPI. The kinetics of asparagine dipeptide and two deamidation sites in mammalian TPI are also investigated using quantum mechanical/molecular mechanical tools with the umbrella sampling technique. Analysis of the simulations has been performed using independent global and local descriptors that can influence the deamidation rates: desolvation effects, backbone acidity, and side chain conformations. Our findings show that all the descriptors add up to favor the primary deamidation site over the secondary one in mammalian TPI: Asn71 deamidates faster because it is more solvent accessible, the adjacent glycine NH backbone acidity is enhanced, and the Asn side chain has a preferential near attack conformation. The crucial impact of the backbone amide acidity of the adjacent glycine on the deamidation rate is shown by kinetic analysis. Our findings also shed light on the effect of high-order structure on deamidation: the deamidation in a small peptide is favored first because of the higher reactivity of the asparagine residue and then because of the stronger stability of the tetrahedral intermediate.