Binding modes of phosphotriesterase-like lactonase complexed with δ-nonanoic lactone and paraoxon using molecular dynamics simulations

3,5,6-trichloro-2-pyridinol intensifies the effect of chlorpyrifos on the paracrine function of Sertoli cells by preventing binding of testosterone and the androgen receptor

3,5,6-Trichloro-2-pyridinol (TCP) is an important biomarker and one of the final metabolites of chlorpyrifos (CPF). TCP inhibits secretion of sex hormones. Similar to CPF, TCP can bind to sex steroid hormone receptors and decrease the secretion of sex hormones. However, little attention has been paid to the ability of TCP and CPF to interfere with androgen receptor (AR) in Sertoli cells. This study aimed to explain how TCP promotes the inhibitory effect of CPF on the paracrine function of Sertoli cells. Western blotting indicated that after 20 weeks of exposure, expression of AR in testes was significantly reduced by CPF. An in vitro assay measured the cytotoxicity of CPF, TCP and diethylphosphate (DEP) on viability of Sertoli cells by Cell Counting Kit-8. CPF cytotoxicity was greater than that of TCP, and TCP cytotoxicity was greater than that of DEP at concentrations of 1000 μmol/L. Western blotting indicated that TCP and CPF both decreased expression of AR and cAMP-response element binding protein phosphorylation, while DEP had no effect in Sertoli cells, which are important in regulating paracrine function of Sertoli cells. The fluorescence measurements and docking studies revealed that testosterone, CPF and TCP showed four types of intermolecular interactions with AR, highlighting alkyl bonds with some of the same amino acids. Compared with testosterone, CPF and TCP also showed significant synergistic interaction with AR. CPF interacted with more amino acids and interaction energy than TCP did. This research elucidates TCP in the antiandrogenic effect of CPF on the paracrine function and suggests that TCP or chemicals with a trichloropyridine structure must be considered during reproductive toxicity assessment of potential environmental pollutants.

How oncogenic mutations activate human MAP kinase 1 (MEK1): a molecular dynamics simulation study

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

Computational Study on the Effect of Inactivating/Activating Mutations on the Inhibition of MEK1 by Trametinib

Activation of the mitogen-activated protein kinase (MAPK) signaling pathway regulated by human MAP kinase 1 (MEK1) is associated with the carcinogenesis and progression of numerous cancers. In addition, two active mutations (P124S and E203K) have been reported to enhance the activity of MEK1, thereby eventually leading to the tumorigenesis of cancer. Trametinib is an MEK1 inhibitor for treating EML4-ALK-positive, EGFR-activated, and KRAS-mutant lung cancers. Therefore, in this study, molecular docking and molecular dynamic (MD) simulations were performed to explore the effects of inactive/active mutations (A52V/P124S and E203K) on the conformational changes of MEK1 and the changes in the interaction of MEK1 with trametinib. Moreover, steered molecular dynamic (SMD) simulations were further utilized to compare the dissociation processes of trametinib from the wild-type (WT) MEK1 and two active mutants (P124S and E203K). As a result, trametinib had stronger interactions with the non-active MEK1 (WT and A52V mutant) than the two active mutants (P124S and E203K). Moreover, two active mutants may make the allosteric channel of MEK1 wider and shorter than that of the non-active types (WT and A52V mutant). Hence, trametinib could dissociate from the active mutants (P124S and E203K) more easily compared with the WT MEK1. In summary, our theoretical results demonstrated that the active mutations may attenuate the inhibitory effects of MEK inhibitor (trametinib) on MEK1, which could be crucial clues for future anti-cancer treatment.

Construction of Novel Aspartokinase Mutant A380I and Its Characterization by Molecular Dynamics Simulation

In this study, a novel monomer aspartokinase (AK) from Corynebacterium pekinense was identified, and its monomer model was constructed. Site 380 was identified by homologous sequencing and monomer model comparison as the key site which was conserved and located around the binding site of the inhibitor Lys. Furthermore, the mutant A380I with enzyme activity 11.32-fold higher than wild type AK (WT-AK), was obtained by site-directed mutagenesis and high throughput screening. In the mutant A380I, the optimal temperature was raised from 26 °C (WT-AK) to 28 °C, the optimal pH remained unchanged at 8.0, and the half-life was prolonged from 4.5 h (WT-AK) to 6.0 h, indicating enhanced thermal stability. The inhibition of A380I was weakened at various inhibitor concentrations and even activated at certain inhibitor concentrations (10 mM of Lys, 5 mM or 10 mM of Lys + Thr, 10 mM of Lys + Met, 5 mM of Lys + Thr + Met). Molecular dynamics simulation results indicated that the occupancy rate of hydrogen bond between A380I and ATP was enhanced, the effect of Lys (inhibitor) on the protein was weakened, and the angle between Ser281-Tyre358 and Asp359-Gly427 was increased after mutation, leading to an open conformation (R-state) that favored the binding of substrate.

Effects of Tyr555 and Trp678 on the processivity of cellobiohydrolase A from Ruminiclostridium thermocellum: A simulation study

Cellobiohydrolase A from Ruminiclostridium thermocellum (Cbh9A) is a processive exoglucanase from family 9 and is an important cellobiohydrolase that hydrolyzes cello-oligosaccharide into cellobiose. Residues Tyr555 and Trp678 considerably affect catalytic activity, but their mechanisms are still unknown. To investigate how the Tyr555 and Trp678 affect the processivity of Cbh9A, conventional molecular dynamics, steered molecular dynamics, and free energy calculation were performed to simulate the processive process of wild type (WT)-Cbh9A, Y555S mutant, and W678G mutant. Analysis of simulation results suggests that the binding free energies between the substrate and WT-Cbh9A are lower than those of Y555S and W678G mutants. The pull forces and energy barrier in Y555S and W678G mutants also reduced significantly during the steered molecular dynamics (SMD) simulation compared with that of the WT-Cbh9A. And the potential mean force calculations showed that the pulling energy barrier of Y555S and W678G mutants is much lower than that of WT-Cbh9A.

Theoretical Study on Zearalenol Compounds Binding with Wild Type Zearalenone Hydrolase and V153H Mutant

Zearalenone hydrolase (ZHD) is the only reported α/β-hydrolase that can detoxify zearalenone (ZEN). ZHD has demonstrated its potential as a treatment for ZEN contamination that will not result in damage to cereal crops. Recent researches have shown that the V153H mutant ZHD increased the specific activity against α-ZOL, but decreased its specific activity to β-ZOL. To understand whyV153H mutation showed catalytic specificity for α-ZOL, four molecular dynamics simulations combining with protein network analysis for wild type ZHD α-ZOL, ZHD β-ZOL, V153H α-ZOL, and V153H β-ZOL complexes were performed using Gromacs software. Our theoretical results indicated that the V153H mutant could cause a conformational switch at the cap domain (residues Gly161–Thr190) to affect the relative position catalytic residue (H242). Protein network analysis illustrated that the V153H mutation enhanced the communication with the whole protein and residues with high betweenness in the four complexes, which were primarily assembled in the cap domain and residues Met241 to Tyr245 regions. In addition, the existence of α-ZOL binding to V153H mutation enlarged the distance from the O_AE atom in α-ZOL to the NE2 atom in His242, which prompted the side chain of H242 to the position with catalytic activity, thereby increasing the activity of V153H on the α-ZOL. Furthermore, α-ZOL could easily form a right attack angle and attack distance in the ZHD and α-ZOL complex to guarantee catalytic reaction. The alanine scanning results indicated that modifications of the residues in the cap domain produced significant changes in the binding affinity for α-ZOL and β-ZOL. Our results may provide useful theoretical evidence for the mechanism underlying the catalytic specificity of ZHD.

Exploration of binding and inhibition mechanism of a small molecule inhibitor of influenza virus H1N1 hemagglutinin by molecular dynamics simulation

Influenza viruses are a major public health threat worldwide. The influenza hemagglutinin (HA) plays an essential role in the virus life cycle. Due to the high conservation of the HA stem region, it has become an especially attractive target for inhibitors for therapeutics. In this study, molecular simulation was applied to study the mechanism of a small molecule inhibitor (MBX2329) of influenza HA. Behaviors of the small molecule under neutral and acidic conditions were investigated, and an interesting dynamic binding mechanism was found. The results suggested that the binding of the inhibitor with HA under neutral conditions facilitates only its intake, while it interacts with HA under acidic conditions using a different mechanism at a new binding site. After a series of experiments, we believe that binding of the inhibitor can prevent the release of HA1 from HA2, further maintaining the rigidity of the HA2 loop and stabilizing the distance between the long helix and short helices. The investigated residues in the new binding site show high conservation, implying that the new binding pocket has the potential to be an effective drug target. The results of this study will provide a theoretical basis for the mechanism of new influenza virus inhibitors.

Structural basis of pesticide detection by enzymatic biosensing: a molecular docking and MD simulation study

Designing of rapid, facile, selective, and cost-effective biosensor technology is a growing area for the detection of various classes of pesticides. The biosensor with these features can be achieved only through the various bio-components using different transducers. This study, therefore, focuses on the usage of molecular docking, specificity tendencies, and capabilities of proteins for the detection of pesticides. Accordingly, the four transducers, acetylcholinesterase (ACH), cytochromes P450 (CYP), glutathione S-transferase (GST), and protein kinase C (PKC) were selected based on their applications including neurotransmitter, metabolism, detoxification enzyme, and protein phosphorylation. Then after molecular docking of the pesticides, fenobucarb, dichlorodiphenyltrichloroethane (DDT), and parathion onto each enzyme, the conformational behavior of the most stable complexes was further analyzed using 50 ns Molecular Dynamics (MD) simulations carried out under explicit water conditions. In the case of protein kinase C (PKC) and cytochrome P450 3A4 enzyme (CYP), the fenobucarb complex showed the most suitable combination of free energy of binding and inhibition constant -4.42 kcal/mol (573.73 μM) and -5.1 kcal/mol (183.49 μM), respectively. Parathion dominated for acetylcholinesterase (ACH) with -4.57 kcal/mol (448.09 μM) and lastly dichlorodiphenyltrichloroethane for glutathione S-transferase (GST), -5.43 kcal/mol (103.88 μM). The RMSD variations were critical for understanding the impact of pesticides as they distinctively influence the energetic attributes of the proteins. Overall, the outcomes from the extensive analysis provide an insight into the structural features of the proteins studied, thereby highlighting their potential use as a substrate in biorecognition sensing of pesticide compounds.