Search

AbstractChinese hazelnut was chosen to become a probable precursor of biological active peptides via computer simulations in this article. There were a large number of bioactive peptides in Chinese hazelnut sequences according to analytical results from the BIOPEP database. The most prominent of these was the inhibitory peptide for dipeptidyl peptidase-IV (DPP-IV; EC 3.4.14.5), which can be used to treat type 2 diabetes, so the theoretical method to obtain DPP-IV inhibitory peptides by hydrolysis with a single or combination of enzymes was studied. Cytotoxicity analysis performed by ToxinPred showed that all of the DPP-IV inhibitory peptides generated from protein hydrolysis were not cytotoxic. Structural interaction fingerprint analysis revealed that Asp663 and Phe357 may be important residues for ligand binding. In order to further understand the inhibitory mechanism of peptide, VR with lowest half maximum inhibitory concentration (IC50) and IPI (inhibitors have been reported) were selected as ligand of DPP-IV to perform steered molecular dynamics simulations and PMF calculations. The results showed that P1 is the preferred (un)binding tunnel for the inhibitors obtained. Our findings help in the development of new DPP-IV inhibitors which were derived from common food.Communicated by Ramaswamy H. Sarma

Communicated by Ramaswamy H. Sarma

AbstractApproximately 30% of all types of human cancers possess a constitutively activated the mitogen-activated protein kinase (MAPK) signaling pathway while MAP kinase 1 (MEK1) is a critical component of this pathway. It has been reported mutations could improve the activity of MEK1 to result in cell proliferation and transformation, which is a known oncogenic event in various cancer types. In this study, eight molecular dynamics simulations, molecular mechanics Poisson–Boltzmann surface area (MM-PBSA), combined with protein structure network were performed to explore the mechanism that mutations activate MEK1. Protein structure networks and hydrogen bonds analysis demonstrated that active mutations broke the interaction between activation segments (residues 216–222) and C-helix (residues 105–121) in MEK1, leading to it transform inactive form to active form. Moreover, hydrogen bond analysis and MM-PBSA calculation indicated that activating mutations decrease the binding affinity between MEK1 and inhibitor to reduce the inhibitory effect of inhibitors. In addition, some active mutations cause structural changes in the Pro-rich loop (residues 261–268) of MEK1. These changes may stabilize the interaction between the MEK1 mutants and the ligands by increasing the number of exposed hydrophobic residues in the active site of MEK1. Our results may provide useful theoretical evidences for the mechanism underlying the role of human MEK1 in human cancers.Communicated by Ramaswamy H. Sarma

AbstractHypoxia-inducible factor-1 (HIF-1) is a transcription factor that plays an important role in the expression of genes, whose function is exerted through protein–protein interactions (PPIs), such as the transcriptional co-activator (CREB)-binding protein (CBP) and p300. Under hypoxic conditions, HIF-1is stabilized and translocated to CBP or p300, leading to the hypoxic response cascade. Furthermore, the PPI between HIF and p300/CBP is a potential cancer target for their role in the hypoxic response. In this study, molecular dynamics (MD) simulation was used to explore the conformational change for the p300 binding to one subunit of HIF-1, namely HIF-1α. Results indicated that HIF-1α–p300 complex was stable during MD simulation. New H-bonds were made in the intra-chain of p300 with HIF-1α binding. Inhibiting the HIF-1α–p300 interaction modulated the HIF-1α identification of selective molecules, which may indicate the target metabolic and cellular processes that enable the survival and growth of tumors in cancer chemotherapy. CAVER 3.0 results suggested that three main tunnels were present, according to helices 1, 2 and 3 of p300. To explore the unbinding pathway for HIF-1α via p300, we selected helices 1, 2 and 3 on the HIF-1α as a new ligand to explore the unbinding pathway via its own tunnel. For helix 1, R368 in p300 formed a H-bond with E816 in HIF1-α. A345 and D346 in p300 formed H-bonds with N803 in HIF-1α. A H-bond existed between K351(p300) and E789 (Hif1-α). These molecules may be the key residues in the unbinding pathway via its tunnel.Communicated by Ramaswamy H. Sarma

Acylaminoacylpeptidase (AAP) belongs to peptidase protein family, which can degrade amyloid β-peptide forms in the brains of patients, and hence leads to Alzheimer’s disease. And so, AAP is considered to be a novel target in the design of drugs against Alzheimer’s disease. In this investigation, six molecular dynamics simulations were used to find that the interaction between the wild-type and R526V AAP with two different substrates (p-nitrophenylcaprylate and Ac-Leu-p-nitroanilide). Our results were as follows: firstly, Ac-Leu-p-nitroanilide bound to R526V AAP to form a more disordered loop (residues 552–562) in the α/β-hydrolase fold like of AAP, which caused an open and inactive AAP domain form, secondly, binding p-nitrophenylcaprylate and Ac-Leu-p-nitroanilide to AAP can decrease the flexibility of residues 225–250, 260–270, and 425–450, in which the ordered secondary structures may contain the suitable geometrical structure and so it is useful to serine attack. Our theoretical results showed that the binding of the two substrates can induce specific conformational changes responsible for the diverse AAP catalytic specificity. These theoretical substrate-induced structural diversities can help explain the abilities of AAPs to recognize and hydrolyze extremely different substrates.

Wild-type phosphotriesterase (PTE) prefers the SP-enantiomers over the corresponding RP-enantiomers by factors ranging from 10 to 90. To satisfy the binding modes of the PTE of SP- and RP-enantiomers, all-atom molecular dynamics simulations were carried out on two paraoxon SPand RPderivatives, namely, Sp-1 and Rp-1. Molecular mechanics Poisson–Boltzmann surface area and molecular mechanics generalized Born surface area (MM-PBSA and MM-GBSA) calculations indicated that His230 in Sp-1-PTE had a closer interaction with the substrate than that in Rp-1-PTE and that such interaction increased the catalytic efficiency of PTE for Sp-1. The steered molecular dynamics simulation indicated that, compared with Sp-1, Rp-1 in the unbinding (binding) may hinder some residue displacement, thus requiring more effort to escape the binding pocket of PTE. In addition, Trp131, Phe306, and Tyr309 are deemed important residues for the Sp-1 unbinding pathway via PTE, whereas Tyr309 alone is considered an important residue for the Rp-1 unbinding pathway. These results demonstrate the possibility of dramatically altering the stereoselectivity and overall reactivity of the native enzyme toward chiral substrates by modifying specific residues located within the active site of PTE.

Phosphotriesterase-like lactonases (PLLs) have received much attention because of their physical and chemical properties. They may have widespread applications in various fields. For example, they show potential for quorum-sensing signaling pathways and organophosphorus (OP) detoxification in agricultural science. However, the mechanism by which PLLs hydrolyze, which involves OP compounds and lactones and a variety of distinct catalytic efficiencies, has only rarely been explored. In the present study, molecular dynamics (MD) simulations were performed to characterize and contrast the structural dynamics of DrPLL, a member of the PLL superfamily in Deinococcus radiodurans, bound to two substrates, δ-nonanoic lactone and paraoxon. It has been observed that there is a 16-fold increase in the catalytic efficiency of the two mutant strains of DrPLL (F26G/C72I) vs. the wild-type enzyme toward the hydrolysis of paraoxon, but an explanation for this behavior is currently lacking. The analysis of the molecular trajectories of DrPLL bound to δ-nonanoic lactone indicated that lactone-induced conformational changes take place in loop 8, which is near the active site. Binding to paraoxon may lead to conformational displacement of loop 1 residues, which could lead to the deformation of the active site and so trigger the entry of the paraoxon into the active site. The efficiency of the F26G/C72I mutant was increased by decreasing the displacement of loop 1 residues and increasing the flexibility of loop 8 residues. These results provide a molecular-level explanation for the experimental behavior.

The bacterial enzyme phosphotriesterase (PTE) exhibits stereoselectivity toward hydrolysis of chiral substrates with a preference for the Spenantiomer. In this work, docking analysis and two explicit-solvent molecular dynamics (MD) simulations were performed to characterize and differentiate the structural dynamics of PTE bound to the Spand Rpparaoxon derivative enantiomers (Rp-1 and Sp-1) hydrolyzed with distinct catalytic efficiencies. Comparative analysis of the molecular trajectories for PTE bound to Rp-1 and Sp-1 suggested that substrate binding induced conformational changes in the loops near the active site. After 100 ns of MD simulation, the Zn β2+metal ion formed hexacoordinated- and tetracoordinated geometries in the Sp-1-PTE and Rp-1-PTE ensembles, respectively. Simulation results further showed that the hydrogen bond between Asp301 and His254 occurred with a higher probability after Sp-1 binding to PTE (47.5%) than that after Rp-1 binding (22.2%). These results provide a qualitative and molecular-level explanation for the 10 orders of magnitude increase in the catalytic efficiency of PTE toward the Spenantiomer of paraoxon.