Search

“Plant-type” ferredoxins (Fds) in the thylakoid membranes of plants, algae, and cyanobacteria possess a single [2Fe-2S] cluster in active sites and mediate light-induced electron transfer from Photosystem I reaction centers to various Fd-dependent enzymes. Structural knowledge of plant-type Fds is relatively limited to static structures, and the detailed behavior of oxidized and reduced Fds has not been fully elucidated. It is important that the investigations of the effects of active-center reduction on the structures and dynamics for elucidating electron-transfer mechanisms. In this study, model systems of oxidized and reduced Fds were constructed from the high-resolution crystal structure of Chlamydomonas reinhardtii Fd1, and three 200 ns molecular dynamics simulations were performed for each system. The force field parameters of the oxidized and reduced active centers were independently obtained using quantum chemical calculations. There were no substantial differences in the global conformations of the oxidized and reduced forms. In contrast, active-center reduction affected the hydrogen-bond network and compactness of the surrounding residues, leading to the increased flexibility of the side chain of Phe61, which is essential for the interaction between Fd and the target protein. These computational results will provide insight into the electron-transfer mechanisms in the Fds.

The coiled-coil domains of the putative yeast cargo receptors Emp46p and Emp47p (Emp46pcc and Emp47pcc) assemble into heterocomplexes at neutral pH. Upon lowering the pH, the complex dissociates and reassembles into homo-oligomers. A glutamate residue (E303) located on the hydrophobic surface of Emp46pcc serves as the pH-sensing switch for assembly and segregation, and we have suggested that its side chains are protonated in the heterocomplex, even at neutral pH. To examine this hypothesis, we constructed two structural models in which the side chains of E303 were negatively charged or protonated and analyzed the effects of these charged states on the structure of the heterocomplex using molecular dynamics (MD) simulations. The calculated structures suggested the side chains of E303 to be protonated in the heterocomplex, even at neutral pH. Based on these computational results, the pH dependence of Emp47pcc homo-oligomer assembly was experimentally modified by a glutamate mutation on its hydrophobic surface. The Q306E mutant of Emp47pcc underwent a structural transition at physiological pH. Our results suggest a method for modifying pH-dependent protein-protein interactions.

CYP3A4 is among the most abundant liver and intestinal drug-metabolizing cytochrome P450 enzymes, contributing to the metabolism of more than 30% of clinically used drugs. Therefore, interindividual variability in CYP3A4 activity is a frequent cause of reduced drug efficacy and adverse effects. In this study, we characterized wild-type CYP3A4 and 40 CYP3A4 variants, including 11 new variants, detected among 4773 Japanese individuals by assessing CYP3A4 enzymatic activities for two representative substrates (midazolam and testosterone). The reduced carbon monoxide–difference spectra of wild-type CYP3A4 and 31 CYP3A4 variants produced with our established mammalian cell expression system were determined by measuring the increase in maximum absorption at 450 nm after carbon monoxide treatment. The kinetic parameters of midazolam and testosterone hydroxylation by wild-type CYP3A4 and 29 CYP3A4 variants (Km, kcat, and catalytic efficiency) were determined, and the causes of their kinetic differences were evaluated by three-dimensional structural modeling. Our findings offer insight into the mechanism underlying interindividual differences in CYP3A4-dependent drug metabolism. Moreover, our results provide guidance for improving drug administration protocols by considering the information on CYP3A4 genetic polymorphisms. SIGNIFICANCE STATEMENT CYP3A4 metabolizes more than 30% of clinically used drugs. Interindividual differences in drug efficacy and adverse-effect rates have been linked to ethnicity-specific differences in CYP3A4 gene variants in Asian populations, including Japanese individuals, indicating the presence of CYP3A4 polymorphisms resulting in the increased expression of loss-of-function variants. This study detected alterations in CYP3A4 activity due to amino acid substitutions by assessing the enzymatic activities of coding variants for two representative CYP3A4 substrates.

Androgen receptor (AR) has a key role in the development and progression of prostate cancer, and AR antagonists are used for its remedy. Recently, carborane derivatives, which are carbon-containing boron clusters have attracted attention as new AR ligands. Here we determined the force field parameters of 10-vertex and 12-vertex p-carborane to facilitate in silico drug design of boron clusters. Then, molecular dynamics (MD) simulations of complexes of AR-carborane derivatives were performed to evaluate the parameters and investigate the influences of carborane derivatives on the three-dimensional structure of AR. Energy profiles were obtained using quantum chemical calculations, and the force-field parameters were determined by curve fitting of the energy profiles. The results of MD simulations indicated that binding of the antagonist-BA341 changed some hydrogen-bond formations involved in the structure and location of helix 12. Those results were consistent with previously reported data. The determined parameters are also useful for refining the structure of the carborane-receptor complex obtained by docking simulations and development of new ligands with carborane cages not only for AR but also for various receptors.

Proteins of modern terrestrial organisms are composed of nearly 20 amino acids; however, the amino acid sets of primitive organisms may have contained fewer than 20 amino acids. Furthermore, the full set of 20 amino acids is not required by some proteins to encode their function. Indeed, simplified variants of Escherichia coli (E. coli) orotate phosphoribosyltransferase (OPRTase) constructed by Akanuma et al. and composed of a limited amino acid set exhibit significant catalytic activity for the growth of E. coli. However, its structural details are currently unclear. Here, we predict the structures of simplified variants of OPRTase using molecular dynamics (MD) simulations and evaluate the accuracy of the MD simulations for simplified proteins. The three-dimensional structure of the wild-type was largely maintained in the simplified variants, but differences in the catalyst loop and C-terminal helix were observed. These results are considered sufficient to elucidate the differences in catalytic activity between the wild-type and simplified OPRTase variants. Thus, using MD simulations to make structural predictions appears to be a useful strategy when investigating non-wild-type proteins composed of reduced amino acid sets.

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

The deamidation of glutamine (Gln) residues, which occurs non-enzymatically under physiological conditions, triggers protein denaturation and aggregation. Gln residues are deamidated via the cyclic glutarimide intermediates to l-α-, d-α-, l-β-, and d-β-glutamate residues. The production of these biologically uncommon amino acid residues is implicated in the pathogenesis of autoimmune diseases. The reaction rate of Gln deamidation is influenced by the C-terminal adjacent (N +1) residue and is highest in the Gln-glycine (Gly) sequence. Here, we investigated the effect of the (N + 1) Gly on the mechanism of Gln deamidation and the activation barrier using quantum chemical calculations. Energy-minima and transition-state geometries were optimized by the B3LYP density functional theory, and MP2 calculations were used to obtain the single-point energy. The calculated activation barrier (85.4 kJ mol−1) was sufficiently low for the reactions occurring under physiological conditions. Furthermore, the hydrogen bond formation between the catalytic ion and the main chain of Gly on the C-terminal side was suggested to accelerate Gln deamidation by stabilizing the transition state.

Dihydropyrimidinase (DHP) is an enzyme that catabolizes the degradation of pyrimidine and fluoropyrimidine drugs such as 5-fluorouracil. DHP deficiency triggers various clinical symptoms and increases the risk of fluoropyrimidine drug toxicity. Various pathogenic variants of DHP cause DHP deficiency, and their catalytic activities have been well studied. However, the three-dimensional structures of DHP variants have not been clarified. In this study, we investigated the effects of mutations on DHP structures using the molecular dynamics simulations. Simulations of the wild type and 10 variants were performed and compared. In the T68R, D81G, G278D, and L337P variants, the flexibilities of structures related to the interaction for oligomer formation increased in comparison with those of the wild type. W117R, T343A, and R412 M mutations affected the structures of stereochemistry gate loops or the substrate-binding pocket. The three-dimensional structures of W360R and G435R variants were suggested to collapse. On the other hand, only slight structural changes were observed in the R181W variant, whose experimentally observed activity was similar to that of the wild type. The computational results are expected to clarify the relationship between clinical mutations and structural effects of drug-metabolizing enzymes.

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.

Cytochrome P450 1A2 (CYP1A2), which accounts for approximately 13% of the total hepatic cytochrome content, catalyzes the metabolic reactions of approximately 9% of frequently used drugs, including theophylline and olanzapine. Substantial inter-individual differences in enzymatic activity have been observed among patients, which could be caused by genetic polymorphisms. Therefore, we functionally characterized 21 novel CYP1A2 variants identified in 4773 Japanese individuals by determining the kinetic parameters of phenacetin O-deethylation. Our results showed that most of the evaluated variants exhibited decreased or no enzymatic activity, which may be attributed to potential structural alterations. Notably, the Leu98Gln, Gly233Arg, Ser380del Gly454Asp, and Arg457Trp variants did not exhibit quantifiable enzymatic activity. Additionally, three-dimensional (3D) docking analyses were performed to further understand the underlying mechanisms behind variant pharmacokinetics. Our data further suggest that despite mutations occurring on the protein surface, accumulating interactions could result in the impairment of protein function through the destabilization of binding regions and changes in protein folding. Therefore, our findings provide additional information regarding rare CYP1A2 genetic variants and how their underlying effects could clarify discrepancies noted in previous phenotypical studies. This would allow the improvement of personalized therapeutics and highlight the importance of identifying and characterizing rare variants.

Isomerized aspartic acid (Asp) residues have previously been identified in various aging tissues, and are suspected to contribute to age-related diseases. Asp-residue isomerization occurs nonenzymatically under physiological conditions, resulting in the formation of three types of isomerized Asp (i.e., L-isoAsp, D-Asp, and D-isoAsp) residues. Asp-residue isomerization often accelerates protein aggregation and insolubilization, making structural biology analyses difficult. Recently, Sakaue et al. reported the synthesis of a ribonuclease A (RNase A) in which Asp121 was artificially replaced with different isomerized Asp residues, and experimentally demonstrated that the enzymatic activities of these artificial mutants were completely lost. However, their structural features have not yet been elucidated. In the present study, the three-dimensional (3D) structures of these artificial-mutant RNases A were predicted using molecular dynamics (MD) simulations. The 3D structures of wild-type and artificial-mutant RNases A were converged by 3000-ns MD simulations. Our computational data show that the structures of the active site and the formation frequencies of the appropriate catalytic dyad structures in the artificial-mutant RNases A were quite different from wild-type RNase A. These computational findings may provide an explanation for the experimental data which show that artificial-mutant RNases A lack enzymatic activity. Herein, MD simulations have been used to evaluate the influences of isomerized Asp residues on the 3D structures of proteins.

A new biological scaffold was produced by replacing the 6π-electron phenyl ring-B of a natural flavone skeleton with a 10π-electron benzothiophene (BT). Since aromatic rings are important for ligand protein interactions, this expansion of the π-electron system of ring-B might change the bioactivity profile. One of the resulting novel natural product-inspired compounds, 2-(benzo[b]thiophen-3-yl)-5-hydroxy-7-isopropoxy-6-methoxyflavone (6), effectively arrested the cell cycle at the G2/M phase and displayed significant antiproliferative effects with IC50 values of 0.05–0.08 μM against multiple human tumor cell lines, including a multidrug resistant line. A structure-activity relationship study revealed that a 10π-electron system with high aromaticity, juxtaposed 4-oxo and 5-hydroxy groups, and 7-alkoxy groups were important for potent antimitotic activity. Interestingly, two BT-flavonols (3-hydroxyflavone), 16 and 20, with 3-hydroxy and 5-alkoxy groups, induced distinct biological profiles affecting the cell cycle at the G1/S phase by inhibition of DNA replication through an interaction with topoisomerase I.

Cytochrome P450 2C9 (CYP2C9) is an important drug-metabolizing enzyme that contributes to the metabolism of approximately 15% of clinically used drugs, including warfarin, which is known for its narrow therapeutic window. Interindividual differences in CYP2C9 enzymatic activity caused by CYP2C9 genetic polymorphisms lead to inconsistent treatment responses in patients. Thus, in this study, we characterized the functional differences in CYP2C9 wild-type (CYP2C9.1), CYP2C9.2, CYP2C9.3, and 12 rare novel variants identified in 4773 Japanese individuals. These CYP2C9 variants were heterologously expressed in 293FT cells, and the kinetic parameters (Km, kcat, Vmax, catalytic efficiency, and CLint) of (S)-warfarin 7-hydroxylation and tolbutamide 4-hydroxylation were estimated. From this analysis, almost all novel CYP2C9 variants showed significantly reduced or null enzymatic activity compared with that of the CYP2C9 wild-type. A strong correlation was found in catalytic efficiencies between (S)-warfarin 7-hydroxylation and tolbutamide 4-hydroxylation among all studied CYP2C9 variants. The causes of the observed perturbation in enzyme activity were evaluated by three-dimensional structural modeling. Our findings could clarify a part of discrepancies among genotype–phenotype associations based on the novel CYP2C9 rare allelic variants and could, therefore, improve personalized medicine, including the selection of the appropriate warfarin dose.

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.

Cytochrome P450 2C9 (CYP2C9) is an important drug-metabolizing enzyme that contributes to the metabolism of approximately 15% of clinically used drugs, including warfarin, which is known for its narrow therapeutic window. Interindividual differences in CYP2C9 enzymatic activity caused by CYP2C9 genetic polymorphisms lead to inconsistent treatment responses in patients. Thus, in this study, we characterized the functional differences in CYP2C9 wild-type (CYP2C9.1), CYP2C9.2, CYP2C9.3, and 12 rare novel variants identified in 4773 Japanese individuals. These CYP2C9 variants were heterologously expressed in 293FT cells, and the kinetic parameters (Km, kcat, Vmax, catalytic efficiency, and CLint) of (S)-warfarin 7-hydroxylation and tolbutamide 4-hydroxylation were estimated. From this analysis, almost all novel CYP2C9 variants showed significantly reduced or null enzymatic activity compared with that of the CYP2C9 wild-type. A strong correlation was found in catalytic efficiencies between (S)-warfarin 7-hydroxylation and tolbutamide 4-hydroxylation among all studied CYP2C9 variants. The causes of the observed perturbation in enzyme activity were evaluated by three-dimensional structural modeling. Our findings could clarify a part of discrepancies among genotype–phenotype associations based on the novel CYP2C9 rare allelic variants and could, therefore, improve personalized medicine, including the selection of the appropriate warfarin dose.

CYP3A4 is among the most abundant liver and intestinal drug-metabolizing cytochrome P450 enzymes, contributing to the metabolism of more than 30% of clinically used drugs. Therefore, interindividual variability in CYP3A4 activity is a frequent cause of reduced drug efficacy and adverse effects. In this study, we characterized wild-type CYP3A4 and 40 CYP3A4 variants, including 11 new variants, detected among 4773 Japanese individuals by assessing CYP3A4 enzymatic activities for two representative substrates (midazolam and testosterone). The reduced carbon monoxide-difference spectra of wild-type CYP3A4 and 31 CYP3A4 variants produced with our established mammalian cell expression system were determined by measuring the increase in maximum absorption at 450 nm after carbon monoxide treatment. The kinetic parameters of midazolam and testosterone hydroxylation by wild-type CYP3A4 and 29 CYP3A4 variants (

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.

In contrast with the common belief that all the amino acid residues in higher organisms are l-forms, d-amino acid residues have been recently detected in various aging tissues. Aspartic acid (Asp) residues are known to be the most prone to stereoinvert via cyclic imide intermediate. Although the glutamic acid (Glu) is similar in chemical structure to Asp, little has been reported to detect d-Glu residues in human proteins. In this study, we investigated the mechanism of the Glu-residue stereoinversion catalyzed by water molecules using B3LYP/6-31+G(d,p) density functional theory calculations. We propose that the Glu-residue stereoinversion proceeds via a cyclic imide intermediate, i.e., glutarimide (GI). All calculations were performed by using a model compound in which a Glu residue was capped with acetyl and methylamino groups on the N- and C-termini, respectively. We found that two water molecules catalyze the three steps involved in the GI formation: iminolization, cyclization, and dehydration. The activation energy required for the Glu residue to form a GI intermediate was estimated to be 32.3 kcal mol-1 , which was higher than that of the experimental Asp-residue stereoinversion. This calculation result suggests that the Glu-residue stereoinversion is not favored under the physiological condition.

Recently, stereoinversions and isomerizations of amino acid residues in the proteins of living beings have been observed. Because isomerized amino acids cause structural changes and denaturation of proteins, isomerizations of amino acid residues are suspected to cause age-related diseases. In this study, AMBER molecular force field parameters were tested by using computationally generated nonapeptides and tripeptides including stereoinverted and/or isomerized amino acid residues. Energy calculations by using density functional theory were also performed for comparison. Although the force field parameters were developed by parameter fitting for l-α-amino acids, the accuracy of the computational results for d-amino acids and β-amino acids was comparable to those for l-α-amino acids. The conformational energies for tripeptides calculated by using density functional theory were reproduced more accurately than those for nonapeptides calculated by using the molecular mechanical force field. The evaluations were performed for the ff99SB, ff03, ff12SB, and the latest ff14SB force field parameters.

Computational docking programs generally predict many candidates for ligand poses, and only one or few poses should be selected for later drug design. In this study, we performed molecular dynamics (MD) simulations on 81 test sets of protein–ligand complexes and determined how well MD simulations evaluated docking poses. Our results suggest that appropriate candidate poses can be chosen from the many docking poses by using MD simulations with explicit solvents. To distinguish the appropriate pose from other candidates, 5-ns and 10-ns MD simulations were effective.

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.

Although nearly all of the amino acids that constitute proteins are l-amino acids, d-amino acid residues in human proteins have been recently reported. d-amino acid residues cause a change in the three-dimensional structure of proteins, and d-aspartic acid (Asp) residues are considered to be one of the causes of age-related diseases. The stereoinversion of Asp residues in peptides and proteins is thought to proceed via a succinimide intermediate; however, it has been reported that stereoinversion can occur even under conditions where a succinimide intermediate cannot be formed. In order to elucidate the non-succinimide-mediated stereoinversion pathway, we investigated the stereoinversion of l-Asp to d-Asp catalysed by phosphate and estimated the activation barrier using B3LYP/6−31+G(d,p) density functional theory (DFT) calculations. For the DFT calculations, a model compound in which the Asp residue is capped with acetyl and methyl-amino groups on the N- and C-termini, respectively, was used. The ...

Spontaneous cyclisation of glutamic acid (Glu) residues located at N-termini in peptides and proteins is called ‘pyroglutamylation’ and is assumed to be involved in several neurodegenerative diseases. Although it has long been believed that N-terminal Glu residues undergo pyroglutamylation enzymatically, it has recently been experimentally confirmed that nonenzymatic pyroglutamylation can proceed in some types of aqueous buffer. However, the detailed mechanism has not been proposed or investigated, and even whether some small-molecule catalysts are required for pyroglutamylation has not been clarified. Therefore, we investigated three types of pyroglutamylation mechanism of N-terminal Glu residues using quantum chemical calculations: in the absence of any catalysts, catalysed by one water molecule, and catalysed by two water molecules. All calculations were performed using N-terminal Glu residues capped with a methylamino group on the C-terminal as a model compound. Optimised energy minima and transition state geometries were obtained using the B3LYP density functional method. The pyroglutamylation mechanism is roughly divided into two steps: cyclisation and dehydration, and the calculated activation barrier was 108 and 107 kJ mol−1 in the two- and three-water-assisted pathways, respectively. The results of computational analysis suggest that water molecules can act as catalysts for pyroglutamylation. The calculated activation barrier of two-water-assisted pyroglutamylation was 108 kJ mol−1, and the results of computational analysis indicate that water molecules can act as catalysts for pyroglutamylation.

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 Moller-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.

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.

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 (H2PO4-) 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 Moller-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 H2PO4- 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 H2PO4- 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 H2PO4--catalyzed Ser residue dehydration can proceed nonenzymatically.

While CYP2D6 allele and phenotype frequencies have been extensively studied, currently, very little ethnically specific data is available regarding the East African and South Pacific region, including Kenya and Vanuatu. The absence of information regarding gene polymorphisms and their resulting clinical effects in these populations may hinder treatment strategies and patient outcome. Given the scarceness of CYP2D6 related data in these populations, the purpose of this study was to perform a pharmacogenomic analysis of the Kenyan and Ni-Vanuatu population and ultimately characterize the enzymatic properties of eight novel CYP2D6 variant proteins expressed in 293FT cells in vitro using dextromethorphan as a substrate. Our study revealed a prevalence of functional alleles in both populations a low frequency for decreased function defining genotypes in the Ni-Vanuatu population, with approximately 36% of our Kenyan subjects presenting substrate-dependent decreased function alleles. Additionally, 6 variants (P171L, G306R, V402L, K1, K2, and K3) showed significantly reduced intrinsic clearance compared to wild-type CYP2D6.1. Our findings aid in efforts to bridge the gap between pharmacogenomic analysis and clinical application, by providing useful information in the development of ethnic-specific strategies as well as stressing the importance of population-specific genotyping when conducting multi-regional clinical trials and designing therapeutic strategies.

Erypoegin K, an isoflavone isolated from the stem bark of Erythrina poeppigiana, has a single chiral carbon in its structure and exists naturally as a racemic mixture. Our previous study showed (S)-erypoegin K selectively exhibits potent anti-proliferative and apoptosis-inducing activity against human leukemia HL-60 cells. To identify the target molecule of (S)-erypoegin K, we employed the human cancer cell panel analysis (termed JFCR39) coupled with a drug sensitivity database of pharmacologically well-characterized drugs for comparison using the COMPARE algorithm. (S)-erypoegin K exhibited a similar profile to that of etoposide, suggesting the molecular target for erypoegin K may be topoisomerase II (Topo II). Subsequent experiments using purified human Topo IIα established that the (S)-isomer selectively stabilizes the cleavage complex composed of double-stranded plasmid DNA and the enzyme. Moreover, (S)-erypoegin K inhibited decatenation of kinetoplast DNA. Molecular docking studies clearly indicated specific binding of the (S)-isomer to the active site of Topo IIα involving hydrogen bonds that help stabilize the cleavage complex. (S)-erypoegin K displayed potent cytotoxic activity against two human gastric cancer cells GCIY and MKN-1 with IC50 values of 0.270 and 0.327 μM, respectively, and induced enzyme activities of caspase 3 and 9. Cell cycle analysis showed marked cell cycle arrest at G2 phase in both cell lines. (S)-erypoegin K also displayed significant antitumor activity toward GCIY xenografted mice. The present study suggests (S)-erypoegin K acts as a Topo II inhibitor to block the G2/M transition of cancer cells.

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 (PSAN) were calculated, because these structural features are considered to relate to the formations of cyclic imide intermediates. RMSFs and PSAN 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. PSAN 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.

In peptides and proteins, nonenzymatic deamidation of asparagine (Asn) residues can trigger some age-related diseases by disrupting the conformation of the biological proteins. In addition, Asn-residue deamidation, which has been observed in various protein preparations, is an important determinant of the quality of protein preparations. In the present study, we investigated the molecular mechanisms of Asn-residue deamidation catalysed by acetic acid, which is frequently used as a buffer in protein preparations. The calculations were conducted using an Asn residue capped with acetyl and methylamino groups on the N- and C-termini, respectively. Energy minima and transition-state geometries were optimised using B3LYP density functional theory (DFT) calculations. The relative energies of all optimised geometries obtained by the MP2 single-point energy calculations were corrected for zero-point energies calculated using the B3LYP DFT method. Asn-residue deamidation was divided roughly into two processes (cyclisation and deammoniation). Computational results indicate that cyclisation is rate-determining. A catalytic acetic acid molecule acted as a proton-transfer mediator in both processes. These results provide useful information for improving formulations of protein preparation.

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.

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 (H2PO4-). 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, H2PO4- 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.

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.

In order to elucidate the reason why glutamic acid residues have lesser racemisation reactivity than asparaginic acid, we investigated the racemisation energy barrier of piperidinedione, which is the presumed intermediate of the isomerisation reaction of L-Glu to D-Glu, by density functional theory calculations. In two-water-molecule-assisted racemisation, the activation barrier for keto–enol isomerisation was 28.1 kcal/mol. The result showed that the activation barrier for the racemisation of glutamic acid residues was not different from that for the racemisation of aspartic acid residues. Thus, glutamic acid residues can possibly cause the racemisation reaction if the cyclic intermediate stably exists.