Search

Recently, inhibitors of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease (Mpro) have been proposed as potential therapeutic agents for COVID-19. Studying effects of amino acid mutations in the conformation of drug targets is necessary for anticipating drug resistance. In this study, with the structure of the SARS-CoV-2 Mpro complexed with a non-covalent inhibitor, we performed molecular dynamics (MD) simulations to determine the conformation of the complex when single amino acid residue in the active site is mutated. As a model of amino acid mutation, we constructed mutant proteins with one residue in the active site mutated to alanine. This method is called virtual alanine scan. The results of the MD simulations showed that the conformation and configuration of the ligand was changed for mutants H163A and E166A, although the structure of the whole protein and of the catalytic dyad did not change significantly, suggesting that mutations in His163 and Glu166 may be linked to drug resistance.

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&rsquo, 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&minus, 1. Some optimized geometries and reaction modes in the aqueous phase were observed that differed from those in the gas phase.

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 &gamma, 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.

Cytochrome P450 (CYP) 2A6 is a monooxygenase involved in the metabolism of various endogenous and exogenous chemicals, such as nicotine and therapeutic drugs. The genetic polymorphisms in CYP2A6 are a cause of individual variation in smoking behavior and drug toxicities. The enzymatic activities of the allelic variants of CYP2A6 were analyzed in previous studies. However, the three-dimensional structures of the mutants were not investigated, and the mechanisms underlying activity reduction remain unknown. In this study, to investigate the structural changes involved in the reduction in enzymatic activities, we performed molecular dynamics simulations for ten allelic mutants of CYP2A6. For the calculated wild type structure, no significant structural changes were observed in comparison with the experimental structure. On the other hand, the mutations affected the interaction with heme, substrates, and the redox partner. In CYP2A6.44, a structural change in the substrate access channel was also observed. Those structural effects could explain the alteration of enzymatic activity caused by the mutations. The results of simulations provide useful information regarding the relationship between genotype and phenotype.

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&middot, mol&minus, 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.