Deamidation of asparagine (Asn) residues is a nonenzymatic post-translational modification of proteins. Asn deamidation is associated with pathogenesis of age-related diseases and hypofunction of monoclonal antibodies. Deamidation rate is known to be affected by the residue following Asn on the carboxyl side and by secondary structure. Information about main-chain conformation of Asn residues is necessary to accurately predict deamidation rate. In this study, the effect of main-chain conformation of Asn residues on deamidation rate was computationally investigated using molecular dynamics (MD) simulations and quantum chemical calculations. The results of MD simulations for γS-crystallin suggested that frequently deamidated Asn residues have common main-chain conformations on the N-terminal side. Based on the simulated structure, initial structures for the quantum chemical calculations were constructed and optimized geometries were obtained using the B3LYP density functional method. Structures that were frequently deamidated had a lower activation energy barrier than that of the little deamidated structure. We also showed that dihydrogen phosphate and bicarbonate ions are important catalysts for deamidation of Asn residues.

Although various higher-order protein structure prediction methods have been developed, almost all of them were developed based on the three-dimensional (3D) structure information of known proteins. Here we predicted the short protein structures by molecular dynamics (MD) simulations in which only Newton’s equations of motion were used and 3D structural information of known proteins was not required. To evaluate the ability of MD simulationto predict protein structures, we calculated seven short test protein (10–46 residues) in the denatured state and compared their predicted and experimental structures. The predicted structure for Trp-cage (20 residues) was close to the experimental structure by 200-ns MD simulation. For proteins shorter or longer than Trp-cage, root-mean square deviation values were larger than those for Trp-cage. However, secondary structures could be reproduced by MD simulations for proteins with 10–34 residues. Simulations by replica exchange MD were performed, but the results were similar to those from normal MD simulations. These results suggest that normal MD simulations can roughly predict short protein structures and 200-ns simulations are frequently sufficient for estimating the secondary structures of protein (approximately 20 residues). Structural prediction method using only fundamental physical laws are useful for investigating non-natural proteins, such as primitive proteins and artificial proteins for peptide-based drug delivery systems.

The nonenzymatic deamidation reactions of asparagine (Asn) and glutamine (Gln) residues in proteins are associated with protein turnover and age-related diseases. The reactions are also believed to provide a molecular clock for biological processes. Although Gln deamidation is assumed to occur through the glutarimide intermediate, the mechanisms for this are unclear because under normal physiological conditions, Gln deamidation occurs relatively less frequently and at a lower rate than Asn deamidation. We investigate the mechanisms underlying glutarimide formation from Gln residues, which proceeds in two steps (cyclization and deammoniation) catalyzed by phosphate and carbonate. We also compare these reactions with noncatalytic mechanisms and water-catalyzed mechanisms. The calculations were performed on the model compound Ace–Gln–Nme (Ace = acetyl, Nme = methylamino) using the density functional theory with the B3LYP/6-31+G(d,p) level of theory. Our results suggest that all the catalysts used in our study can mediate the proton relays required for glutarimide formation. We further determined that the calculated activation barriers of the reactions catalyzed by phosphate ions (115 kJ mol^–1) and carbonate ions (112 kJ mol^–1) are sufficiently low for the reactions to occur under normal physiological conditions. We also show that nucleophilic enhancement of Nme nitrogen is essential for the cyclization of Gln residues.

Objective The blood–brain barrier (BBB) limits the entry of some therapeutics into the brain, resulting in reduced efficacy. BBB-opening techniques have been developed to enhance the entry into the brain. However, a noninvasive, highly sensitive and quantitative method for evaluating the changes in BBB permeability induced by such techniques is needed to optimize treatment protocols. We evaluated 2-amino-[3-¹¹C]isobutyric acid ([3-¹¹C]AIB) as a PET probe to quantify BBB permeability in model rats. Methods BBB opening was induced by a lipopolysaccharide injection or focused ultrasound (FUS) sonication. [3-¹¹C]AIB distribution in the brain was evaluated by autoradiography and PET and compared with that of Evans blue, a traditional BBB permeability marker. Kinetics of [3-¹¹C]AIB was compared with that of gadolinium-diethylenetriamine pentaacetic acid (Gd-DTPA)-enhanced MRI. The unidirectional blood–brain transfer constant (K_i) of [3-¹¹C]AIB was estimated using the Patlak plot. Results [3-¹¹C]AIB uptake in the lesion area was significantly higher than that in the control area and radioactivity colocalized with Evans blue in both models. [3-¹¹C]AIB uptake in the FUS-sonicated region decreased over time after sonication. The ratio of [3-¹¹C]AIB accumulation in the FUS-treated to the contralateral side increased during the experimental period, whereas that of the Gd-DTPA intensity reached a maximum at 10 min after injection and decreased thereafter. The [3-¹¹C]AIB K_i values were significantly higher in the lesion area than the control area. Conclusion [3-¹¹C]AIB PET is a promising, highly sensitive and quantitative imaging method for assessment of BBB permeability.

The nonenzymatic peptide bond cleavage at the C-terminal side of Asn residues is a protein post-translational modification that occurs under physiological conditions. This reaction proceeds much slower than the deamidation of the Asn side chain and causes denaturation and hypofunction of proteins. The peptide bond cleavage of Asn is detected primarily in crystallins and aquaporin 0 in the eye lens. Therefore, cleavage is thought to be involved in age-related cataracts. In this study, to clarify the mechanism underlying succinimide formation for the peptide bond cleavage of the Asn residue, we performed quantum chemical calculations on the model compound Ace-Asn-Gly-Nme (Ace = acetyl and Nme = methylamino). The density functional theory with the B3LYP/6-31+G(d,p) level of theory was used to obtain optimized geometries. The results suggested that the reaction proceeds through two steps, cyclization and C-terminal fragment release, and the required proton transfers can be mediated by H₂PO₄^– and HCO₃^– ions. The conformational change of the main chain on the N-terminal side of Asn was needed for the C-terminal fragmentation step, and a separate conformational change at the C-terminal side was required for the cyclization step. Furthermore, the calculated activation barriers of the reactions catalyzed by the H₂PO₄^– ion (130 kJ mol^–1) and the HCO₃^– ion (123 kJ mol^–1) were sufficiently low for the reactions to occur under normal physiological conditions.

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.