A computational analysis of interaction mechanisms of peptide and non-peptide inhibitors with MDMX based on molecular dynamics simulation

Probing molecular mechanism of inhibitor bindings to bromodomain-containing protein 4 based on molecular dynamics simulations and principal component analysis

It is well known that bromodomain-containing protein 4 (BRD4) has been thought as a promising target utilized for treating various human diseases, such as inflammatory disorders, malignant tumours, acute myelogenous leukaemia (AML), bone diseases, etc. For this study, molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA) were integrated together to uncover binding modes of inhibitors 8P9, 8PU, and 8PX to BRD4(1). The results obtained from binding free energy calculations show that van der Waals interactions act as the main regulator in bindings of inhibitors to BRD4(1). The information stemming from PCA reveals that inhibitor associations extremely affect conformational changes, internal dynamics, and movement patterns of BRD4(1). Residue-based free energy decomposition method was wielded to unveil contributions of independent residues to inhibitor bindings and the data signify that hydrogen bonding interactions and hydrophobic interactions are decisive factors affecting bindings of inhibitors to BRD4(1). Meanwhile, eight residues Trp81, Pro82, Val87, Leu92, Leu94, Cys136, Asn140, and Ile146 are recognized as the common hot interaction spots of three inhibitors with BRD4(1). The results from this work are expected to provide a meaningfully theoretical guidance for design and development of effective inhibitors inhibiting of the activity of BRD4.

Electronic polarization stabilizes tertiary structure prediction of HP-36

Molecular dynamic (MD) simulations with both implicit and explicit solvent models have been carried out to study the folding dynamics of HP-36 protein. Starting from the extended conformation, the secondary structure of all three helices in HP-36 was formed in about 50 ns and remained stable in the remaining simulation. However, the formation of the tertiary structure was difficult. Although some intermediates were close to the native structure, the overall conformation was not stable. Further analysis revealed that the large structure fluctuation of loop and hydrophobic core regions was devoted mostly to the instability of the structure during MD simulation. The backbone root-mean-square deviation (RMSD) of the loop and hydrophobic core regions showed strong correlation with the backbone RMSD of the whole protein. The free energy landscape indicated that the distribution of main chain torsions in loop and turn regions was far away from the native state. Starting from an intermediate structure extracted from the initial AMBER simulation, HP-36 was found to generally fold to the native state under the dynamically adjusted polarized protein-specific charge (DPPC) simulation, while the peptide did not fold into the native structure when AMBER force filed was used. The two best folded structures were extracted and taken into further simulations in water employing AMBER03 charge and DPPC for 25 ns. Result showed that introducing polarization effect into interacting potential could stabilize the near-native protein structure.

Molecular modeling and molecular dynamics simulation studies on the interactions of hydroxylated polychlorinated biphenyls with estrogen receptor-β

Endocrine-disrupting chemicals have attracted great concern. As major metabolites of polychlorinated biphenyls (PCBs), hydroxylated polychlorinated biphenyls (HO-PCBs) may disrupt estrogen hormone status because of their structural similarity to estrogen endogenous compounds. However, interactions between HO-PCBs and estrogen receptors (ERs) are not fully understood. In the present work, a molecular modeling study combining molecular docking, molecular dynamics simulations, and binding free energy calculations was performed to characterize the interactions of three HO-PCBs (4'-HO-PCB50, 2'-HO-PCB65, and 4'-HO-PCB69) having much different estrogenic activities with ERβ. Docking results showed that binding between ligands and ERβ was stabilized by hydrogen bond and hydrophobic interactions. The binding free energies of three ligands with ERβ were calculated, and further binding free energy decomposition analysis indicated that the dominating driving force of the binding between the ligands and ERβ was the van der Waals interaction. Some key residues, such as Leu298, Phe356, Gly472, His475, and Leu476, played important roles in ligand-receptor interactions by forming hydrophobic and hydrogen bond interactions with ligands. The results may be beneficial to increase understanding of the interactions between HO-PCBs and ERβ.

Molecular dynamic simulation insights into the normal state and restoration of p53 function

As a tumor suppressor protein, p53 plays a crucial role in the cell cycle and in cancer prevention. Almost 50 percent of all human malignant tumors are closely related to a deletion or mutation in p53. The activity of p53 is inhibited by over-active celluar antagonists, especially by the over-expression of the negative regulators MDM2 and MDMX. Protein-protein interactions, or post-translational modifications of the C-terminal negative regulatory domain of p53, also regulate its tumor suppressor activity. Restoration of p53 function through peptide and small molecular inhibitors has become a promising strategy for novel anti-cancer drug design and development. Molecular dynamics simulations have been extensively applied to investigate the conformation changes of p53 induced by protein-protein interactions and protein-ligand interactions, including peptide and small molecular inhibitors. This review focuses on the latest MD simulation research, to provide an overview of the current understanding of interactions between p53 and its partners at an atomic level.