Shuichi Fukuyoshi

Molecular Dynamics Simulations to Investigate the Influences of Amino Acid Mutations on Protein Three-Dimensional Structures of Cytochrome P450 2D6.1, 2, 10, 14A, 51, and 62

Many natural mutants of the drug metabolizing enzyme cytochrome P450 (CYP) 2D6 have been reported. Because the enzymatic activities of many mutants are different from that of the wild type, the genetic polymorphism of CYP2D6 plays an important role in drug metabolism. In this study, the molecular dynamics simulations of the wild type and mutants of CYP2D6, CYP2D6.1, 2, 10, 14A, 51, and 62 were performed, and the predictions of static and dynamic structures within them were conducted. In the mutant CYP2D6.10, 14A, and 61, dynamic properties of the F-G loop, which is one of the components of the active site access channel of CYP2D6, were different from that of the wild type. The F-G loop acted as the “hatch” of the channel, which was closed in those mutants. The structure of CYP2D6.51 was not converged by the simulation, which indicated that the three-dimensional structure of CYP2D6.51 was largely different from that of the wild type. In addition, the intramolecular interaction network of CYP2D6.10, 14A, and 61 was different from that of the wild type, and it is considered that these structural changes are the reason for the decrease or loss of enzymatic activities. On the other hand, the static and dynamic properties of CYP2D6.2, whose activity was normal, were not considerably different from those of the wild type.

Prediction of three-dimensional structures and structural flexibilities of wild-type and mutant cytochrome P450 1A2 using molecular dynamics simulations.

In this study, we investigated the effect of genetic polymorphism on the three-dimensional (3D) conformation of cytochrome P450 1A2 (CYP1A2) using molecular dynamics (MD) simulations. CYP1A2, a major drug-metabolizing enzyme among cytochrome P450 enzymes (CYPs), is known to have many variant alleles. The genetic polymorphism of CYP1A2 may cause individual differences in the pharmacokinetics of medicines. By performing 100ns or longer MD simulations, we investigated the influence of amino acid mutation on the 3D structures and the dynamic properties of proteins. The results show that the static structures were changed by the mutations of amino acid residues, not only near the mutated residues but also in distant portions of the proteins. Moreover, the mutation of only one amino acid was shown to change the structural flexibility of proteins, which may influence the substrate recognition and enzymatic activity. Our results clearly suggest that it is necessary to investigate the dynamic property as well as the static 3D structure for understanding the change of the enzymatic activity of mutant CYP1A2.

Predicting Three-Dimensional Conformations of Peptides Constructed of Only Glycine, Alanine, Aspartic Acid, and Valine

The GADV hypothesis is a form of the protein world hypothesis, which suggests that life originated from proteins (Lacey et al. 1999; Ikehara 2002; Andras 2006). In the GADV hypothesis, life is thought to have originated from primitive proteins constructed of only glycine, alanine, aspartic acid, and valine ([GADV]-proteins). In this study, the three-dimensional (3D) conformations of randomly generated short [GADV]-peptides were computationally investigated using replica-exchange molecular dynamics (REMD) simulations (Sugita and Okamoto 1999). Because the peptides used in this study consisted of only 20 residues each, they could not form certain 3D structures. However, the conformational tendencies of the peptides were elucidated by analyzing the conformational ensembles generated by REMD simulations. The results indicate that secondary structures can be formed in several randomly generated [GADV]-peptides. A long helical structure was found in one of the hydrophobic peptides, supporting the conjecture of the GADV hypothesis that many peptides aggregated to form peptide multimers with enzymatic activity in the primordial soup. In addition, these results indicate that REMD simulations can be used for the structural investigation of short peptides.

Effects of initial settings on computational protein–ligand docking accuracies for several docking programs

In this study, the influences of initial settings, i.e. initial conformations, configurations and docking parameters, on docking results were investigated. The conformations used in the study were generated by the CAMDAS program. After the conformational search calculations, five structures were selected from the conformer groups according to their conformation energies and root mean square deviations against crystal structures; for example, the lowest energy conformer, as well as the closest and farthest conformers to the crystal structure, was retrieved. Several docking parameter settings were used (default, high speed, generating 50 poses). In this study, docking calculations were conducted using the GOLD, eHiTS, AutoDock, AutoDock vina, FRED and DOCK programs. The success rates of GOLD, eHiTS and FRED were better than those of AutoDock, AutoDock vina and DOCK. The docking results using the farthest conformations were worse than those obtained using other conformations, indicating that some conformation search for the ligand molecule should be performed before the docking calculations.

Predicting the structures of complexes between phosphoinositide 3-kinase (PI3K) and romidepsin-related compounds for the drug design of PI3K/histone deacetylase dual inhibitors using computational docking and the ligand-based drug design approach.

Predictions of the three-dimensional (3D) structures of the complexes between phosphoinositide 3-kinase (PI3K) and two inhibitors were conducted using computational docking and the ligand-based drug design approach. The obtained structures were refined by structural optimizations and molecular dynamics (MD) simulations. The ligands were located deep inside the ligand binding pocket of the p110α subunit of PI3K, and the hydrogen bond formations and hydrophobic effects of the surrounding amino acids were predicted. Although rough structures were obtained for the PI3K-inhibitor complexes before the MD simulations, the refinement of the structures by these simulations clarified the hydrogen bonding patterns of the complexes.

Evaluations of the conformational search accuracy of CAMDAS using experimental three-dimensional structures of protein-ligand complexes

CAMDAS is a conformational search program, through which high temperature molecular dynamics (MD) calculations are carried out. In this study, the conformational search ability of CAMDAS was evaluated using structurally known 281 protein-ligand complexes as a test set. For the test, the influences of initial settings and initial conformations on search results were validated. By using the CAMDAS program, reasonable conformations whose root mean square deviations (RMSDs) in comparison with crystal structures were less than 2.0 A could be obtained from 96% of the test set even though the worst initial settings were used. The success rate was comparable to those of OMEGA, and the errors of CAMDAS were less than those of OMEGA. Based on the results obtained using CAMDAS, the worst RMSD was around 2.5 A, although the worst value obtained was around 4.0 A using OMEGA. The results indicated that CAMDAS is a robust and versatile conformational search method and that it can be used for a wide variety of small molecules. In addition, the accuracy of a conformational search in relation to this study was improved by longer MD calculations and multiple MD simulations.

Theoretical study on keto–enol tautomerisation of glutarimide for exploration of the isomerisation reaction pathway of glutamic acid in proteins using density functional theory

ABSTRACT 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. GRAPHICAL ABSTRACT

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

Prediction of binding modes between protein L-isoaspartyl (D-aspartyl) O-methyltransferase and peptide substrates including isomerized aspartic acid residues using in silico analytic methods for the substrate screening.

Because the aspartic acid (Asp) residues in proteins are occasionally isomerized in the human body, not only l-α-Asp but also l-β-Asp, D-α-Asp and D-β-Asp are found in human proteins. In these isomerized aspartic acids, the proportion of D-β-Asp is the largest and the proportions of l-β-Asp and D-α-Asp found in human proteins are comparatively small. To explain the proportions of aspartic acid isomers, the possibility of an enzyme able to repair l-β-Asp and D-α-Asp is frequently considered. The protein L-isoaspartyl (D-aspartyl) O-methyltransferase (PIMT) is considered one of the possible repair enzymes for l-β-Asp and D-α-Asp. Human PIMT is an enzyme that recognizes both l-β-Asp and D-α-Asp, and catalyzes the methylation of their side chains. In this study, the binding modes between PIMT and peptide substrates containing l-β-Asp or D-α-Asp residues were investigated using computational protein-ligand docking and molecular dynamics simulations. The results indicate that carboxyl groups of both l-β-Asp and D-α-Asp were recognized in similar modes by PIMT and that the C-terminal regions of substrate peptides were located in similar positions on PIMT for both the l-β-Asp and D-α-Asp peptides. In contrast, for peptides containing l-α-Asp or D-β-Asp residues, which are not substrates of PIMT, the computationally constructed binding modes between PIMT and peptides greatly differed from those between PIMT and substrates. In the nonsubstrate peptides, not inter- but intra-molecular hydrogen bonds were observed, and the conformations of peptides were more rigid than those of substrates. Thus, the in silico analytical methods were able to distinguish substrates from nonsubstrates and the computational methods are expected to complement experimental analytical methods.

Computational studies on non-succinimide-mediated stereoinversion mechanism of aspartic acid residues assisted by phosphate

ABSTRACT 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 calculated activation barrier was not excessively high for the stereoinversion to occur in vivo. Therefore, this stereoinversion mechanism may compete with the succinimide-mediated mechanism. GRAPHICAL ABSTRACT