Search

β-N-Acetyl-d-hexosaminidase from Ostrinia furnacalis (OfHex1) is a new target for the design of insecticides. Although some of its inhibitors have been found, there is still no commercial drug available at present. The residence time of the ligand may be important for its pharmacodynamic effect. However, the unbinding routes of ligands from OfHex1 still remain largely unexplored. In the present study, we first simulated the six dissociation routes of N,N,N-trimethyl-d-glucosamine-chitotriomycin (TMG-chitotriomycin, a highly selective inhibitor of OfHex1) from the active pocket of OfHex1 by steered molecular dynamics simulations. By comparing the potential of mean forces (PMFs) of six routes, Route 1 was considered as the most possible route with the lowest energy barrier. Furthermore, the structures of six different states for Route 1 were snapshotted, and the key amino acid residues affecting the dissociated time were analyzed in the unbinding pathway. Moreover, we also analyzed the “open–close„ mechanism of Glu368 and Trp448 and found that their conformational changes directly affected the dissociation of TMG-chitotriomycin. Our findings would be helpful to understanding and identifying novel inhibitors against OfHex1 from virtual screening or lead-optimization.

ToxIN is a triangular structure formed by three protein toxins (ToxNs) and three specific noncoding RNA antitoxins (ToxIs). To respond to stimuli, ToxI is preferentially degraded, releasing the ToxN. Thus, the dynamic character is essential in the normal function interactions between ToxN and ToxI. Here, equilibrated molecular dynamics (MD) simulations were performed to study the stability of ToxN and ToxI. The results indicate that ToxI adjusts the conformation of 3′ and 5′ termini to bind to ToxN. Steered molecular dynamics (SMD) simulations combined with the recently developed thermodynamic integration in 3nD (TI3nD) method were carried out to investigate ToxN unbinding from the ToxIN complex. The potentials of mean force (PMFs) and atomistic pictures suggest the unbinding mechanism as follows: (1) dissociation of the 5′ terminus from ToxN, (2) missing the interactions involved in the 3′ terminus of ToxI without three nucleotides (G31, A32, and A33), (3) starting to unfold for ToxI, (4) leaving the binding package of ToxN for three nucleotides of ToxI, (5) unfolding of ToxI. This work provides information on the structure-function relationship at the atomistic level, which is helpful for designing new potent antibacterial drugs in the future.

SCH772984, VTX-11e, FR180204 and 5-iTU are four promising inhibitors targeting ERK2 kinase with high bioactivity. These four inhibitors also have different residence times and binding modes with the ERK2 kinase. Revealing the molecular mechanism of different residence times of ERK2 inhibitors is helpful for designing more efficient inhibitors. The molecular mechanics/generalized Born surface area (MM/GBSA) method was used to calculate the binding free energy and identify the key residues for the ERK2 protein binding to the four inhibitors. Steered molecular dynamics (SMD) and adaptive biasing force (ABF) simulations were employed to investigate the molecular mechanism behind this difference in residence time and binding mode. The binding free energy decomposition by the MM/GBSA method reveals the residues Y27, K45, I47, P49, Y55, R58, T59, Q96 and G160 located around the allosteric binding pocket play an important role in determining the longer residence time of SCH772984. The results from the SMD and the ABF simulations show SCH772984 has different unbinding mechanism compared with the other three inhibitors. SCH772984 needs to overcome two energy barriers: one is the π-π stacking interaction formed by the piperazine-phenyl-pyrimidine of SCH772984 and the residue Y55 of the ERK2 kinase; the other is the hydrophobic interaction at the ATP active site. VTX-11e, FR180204 and 5-iTU just need to overcome the hydrophobic interaction at the ATP active site. Our simulation results are useful to understand the interaction mechanism between four inhibitors and ERK2 kinase and are helpful for designing more potent ERK2 inhibitors. [ABSTRACT FROM AUTHOR]

Inhibition of p53-MDM2 interaction by small molecules is considered to be a promising approach to re-activate wild-type p53 for tumor suppression. Several inhibitors of the MDM2-p53 interaction were designed and studied by the experimental methods and the molecular dynamics simulation. However, the unbinding mechanism was still unclear. The steered molecular dynamics simulations combined with Brownian dynamics fluctuation-dissipation theorem were employed to obtain the free-energy landscape of unbinding between MDM2 and their four ligands. It was shown that compounds 4 and 8 dissociate faster than compounds 5 and 7. The absolute binding free energies for these four ligands are in close agreement with experimental results. The open movement of helix II and helix IV in the MDM2 protein-binding pocket upon unbinding is also consistent with experimental MDM2-unbound conformation. We further found that different binding mechanisms among different ligands are associated with H-bond with Lys51 and Glu25. These mechanistic results may be useful for improving ligand design. [ABSTRACT FROM AUTHOR]

β-N-Acetyl-d-hexosaminidase from Ostrinia furnacalis (OfHex1) is a new target for the design of insecticides. Although some of its inhibitors have been found, there is still no commercial drug available at present. The residence time of the ligand may be important for its pharmacodynamic effect. However, the unbinding routes of ligands from OfHex1 still remain largely unexplored. In the present study, we first simulated the six dissociation routes of N,N,N-trimethyl-d-glucosamine-chitotriomycin (TMG-chitotriomycin, a highly selective inhibitor of OfHex1) from the active pocket of OfHex1 by steered molecular dynamics simulations. By comparing the potential of mean forces (PMFs) of six routes, Route 1 was considered as the most possible route with the lowest energy barrier. Furthermore, the structures of six different states for Route 1 were snapshotted, and the key amino acid residues affecting the dissociated time were analyzed in the unbinding pathway. Moreover, we also analyzed the "open–close" mechanism of Glu368 and Trp448 and found that their conformational changes directly affected the dissociation of TMG-chitotriomycin. Our findings would be helpful to understanding and identifying novel inhibitors against OfHex1 from virtual screening or lead-optimization. [ABSTRACT FROM AUTHOR]

ToxIN is a triangular structure formed by three protein toxins (ToxNs) and three specific noncoding RNA antitoxins (ToxIs). To respond to stimuli, ToxI is preferentially degraded, releasing the ToxN. Thus, the dynamic character is essential in the normal function interactions between ToxN and ToxI. Here, equilibrated molecular dynamics (MD) simulations were performed to study the stability of ToxN and ToxI. The results indicate that ToxI adjusts the conformation of 3′ and 5′ termini to bind to ToxN. Steered molecular dynamics (SMD) simulations combined with the recently developed thermodynamic integration in 3nD (TI3nD) method were carried out to investigate ToxN unbinding from the ToxIN complex. The potentials of mean force (PMFs) and atomistic pictures suggest the unbinding mechanism as follows: (1) dissociation of the 5′ terminus from ToxN, (2) missing the interactions involved in the 3′ terminus of ToxI without three nucleotides (G31, A32, and A33), (3) starting to unfold for ToxI, (4) leaving the binding package of ToxN for three nucleotides of ToxI, (5) unfolding of ToxI. This work provides information on the structure-function relationship at the atomistic level, which is helpful for designing new potent antibacterial drugs in the future. [ABSTRACT FROM AUTHOR]

Ecdysone receptor (EcR) is an important target for pesticide design. Ligand binding regulates EcR transcriptional activity similar to other nuclear receptors; however, the pathways by which ligands enter and leave the EcR remain poorly understood. Here, we performed computational studies to identify unbinding pathways of an ecdysone agonist [the selective ecdysone agonist, BYI06830] from the EcR ligand binding domain (EcR LBD). BYI06830 can dissociate from EcR LBD via four different pathways with little effect on receptor structure. By comparing the potential of mean force (PMF) of four pathways, path 2 was considered to be the most likely exit path for BYI06830, which was located in the cleft formed by the H3-H4 loop, H6-H7 loop, and the H11 C-terminus. Furthermore, structural features along path 2 were analyzed and the structural snapshots of the metastable and transition states were isolated to illustrate the unbinding mechanism of ecdysone agonist from EcR LBD.