Modeling a new water channel that allows SET9 to dimethylate p53

New Mechanism for the Apoptosis of Human Neuroblastoma Cells by the Interaction between Fluorene-9-Bisphenol and the G Protein-Coupled Estrogen Receptor 1

Fluorene-9-bisphenol (BHPF) is an emerging contaminant. Presently, there is no report on its interaction with G protein-coupled estrogen receptor 1 (GPER). By using an integrated toxicity research scenario that combined theoretical study with experimental methods, BHPF was found to inhibit the GPER-mediated effect via direct receptor binding. Molecular dynamics simulations found that Trp2726.48 and Glu2756.51 be the key amino acids of BHPF binding with GPER. Moreover, the calculation indicated that BHPF was a suspected GPER inhibitor, which neither can activate GPER nor is able to form water channels of GPER. The role of two residues was successfully verified by following gene knockout and site-directed mutagenesis assays. Further in vitro assays showed that BHPF could attenuate the increase in intracellular concentration of free Ca2+ induced by G1-activated GPER. Besides, BHPF showed an enhanced cytotoxicity compared with G15, indicating that BHPF might be a more potent GPER inhibitor than G15. In addition, a statistically significant effect on the mRNA level of GPER was observed for BHPF. In brief, the present study proposes that BHPF be a GPER inhibitor, and its GPER molecular recognition mechanism has been revealed, which is of great significance for the health risk and assessment of BHPF.

Approaching the catalytic mechanism of protein lysine methyltransferases by biochemical and simulation techniques

Protein lysine methyltransferases (PKMTs) transfer up to three methyl groups to the side chains of lysine residues in proteins and fulfill important regulatory functions by controlling protein stability, localization and protein/protein interactions. The methylation reactions are highly regulated, and aberrant methylation of proteins is associated with several types of diseases including neurologic disorders, cardiovascular diseases, and various types of cancer. This review describes novel insights into the catalytic machinery of various PKMTs achieved by the combined application of biochemical experiments and simulation approaches during the last years, focusing on clinically relevant and well-studied enzymes of this group like DOT1L, SMYD1-3, SET7/9, G9a/GLP, SETD2, SUV420H2, NSD1/2, different MLLs and EZH2. Biochemical experiments have unraveled many mechanistic features of PKMTs concerning their substrate and product specificity, processivity and the effects of somatic mutations observed in PKMTs in cancer cells. Structural data additionally provided information about the substrate recognition, enzyme-substrate complex formation, and allowed for simulations of the substrate peptide interaction and mechanism of PKMTs with atomistic resolution by molecular dynamics and hybrid quantum mechanics/molecular mechanics methods. These simulation technologies uncovered important mechanistic details of the PKMT reaction mechanism including the processes responsible for the deprotonation of the target lysine residue, essential conformational changes of the PKMT upon substrate binding, but also rationalized regulatory principles like PKMT autoinhibition. Further developments are discussed that could bring us closer to a mechanistic understanding of catalysis of this important class of enzymes in the near future. The results described here illustrate the power of the investigation of enzyme mechanisms by the combined application of biochemical experiments and simulation technologies.

Histone Lysine Methylation Modification and Its Role in Vascular Calcification

Histone methylation is an epigenetic change mediated by histone methyltransferase, and has been connected to the beginning and progression of several diseases. The most common ailments that affect the elderly are cardiovascular and cerebrovascular disorders. They are the leading causes of death, and their incidence is linked to vascular calcification (VC). The key mechanism of VC is the transformation of vascular smooth muscle cells (VSMCs) into osteoblast-like phenotypes, which is a highly adjustable process involving a variety of complex pathophysiological processes, such as metabolic abnormalities, apoptosis, oxidative stress and signalling pathways. Many researchers have investigated the mechanism of VC and related targets for the prevention and treatment of cardiovascular and cerebrovascular diseases. Their findings revealed that histone lysine methylation modification may play a key role in the various stages of VC. As a result, a thorough examination of the role and mechanism of lysine methylation modification in physiological and pathological states is critical, not only for identifying specific molecular markers of VC and new therapeutic targets, but also for directing the development of new related drugs. Finally, we provide this review to discover the association between histone methylation modification and VC, as well as diverse approaches with which to investigate the pathophysiology of VC and prospective treatment possibilities.

Resveratrol induces p53 in colorectal cancer through SET7/9

Resveratrol is one of the most promising phytoalexins for use as an anti-cancer agent, which is present in the skin of red grapes and berries. Resveratrol has been demonstrated to modulate a number of signalling pathways that are involved in carcinogenesis. In the present study, the function of resveratrol as a pro-apoptotic agent in colorectal cancer cell lines, including HCT116, CO115 and SW48, was investigated. The results revealed that resveratrol supressed cell viability. Additionally, resveratrol enhanced the expression of tumour protein p53 (p53) and p53 target genes, including Bcl2 associated X, apoptosis regulator and Bcl2 binding component 3 that have a pivotal role in p53-dependent apoptosis. Furthermore, treating cells with resveratrol upregulated SET domain containing lysine methyltransferase 7/9 (SET7/9) expression, which positively regulates p53 through its mono-methylation at lysine 372, compared with untreated cells. Furthermore, treating cells with resveratrol induced the expression of apoptotic markers including cleaved caspase-3 and poly (ADP-ribose) polymerases (PARP) compared with untreated cells. However, the genetic knockdown of SET7/9 by short hairpin RNA attenuated the resveratrol-driven overexpression of p53, cleaved caspase-3 and PARP. Collectively, these results reveal the molecular mechanisms by which resveratrol induces p53 stability in colon cancer that results in the activation of p53-mediated apoptosis.

Chemical and Biochemical Perspectives of Protein Lysine Methylation

Protein lysine methylation is a distinct posttranslational modification that causes minimal changes in the size and electrostatic status of lysine residues. Lysine methylation plays essential roles in regulating fates and functions of target proteins in an epigenetic manner. As a result, substrates and degrees (free versus mono/di/tri) of protein lysine methylation are orchestrated within cells by balanced activities of protein lysine methyltransferases (PKMTs) and demethylases (KDMs). Their dysregulation is often associated with neurological disorders, developmental abnormalities, or cancer. Methyllysine-containing proteins can be recognized by downstream effector proteins, which contain methyllysine reader domains, to relay their biological functions. While numerous efforts have been made to annotate biological roles of protein lysine methylation, limited work has been done to uncover mechanisms associated with this modification at a molecular or atomic level. Given distinct biophysical and biochemical properties of methyllysine, this review will focus on chemical and biochemical aspects in addition, recognition, and removal of this posttranslational mark. Chemical and biophysical methods to profile PKMT substrates will be discussed along with classification of PKMT inhibitors for accurate perturbation of methyltransferase activities. Semisynthesis of methyllysine-containing proteins will also be covered given the critical need for these reagents to unambiguously define functional roles of protein lysine methylation.

Ligand induced change of β2 adrenergic receptor from active to inactive conformation and its implication for the closed/open state of the water channel: insight from molecular dynamics simulation, free energy calculation and Markov state model analysis

The reported crystal structures of β2 adrenergic receptor (β2AR) reveal that the open and closed states of the water channel are correlated with the inactive and active conformations of β2AR. However, more details about the process by which the water channel states are affected by the active to inactive conformational change of β2AR remain illusive. In this work, molecular dynamics simulations are performed to study the dynamical inactive and active conformational change of β2AR induced by inverse agonist ICI 118,551. Markov state model analysis and free energy calculation are employed to explore the open and close states of the water channel. The simulation results show that inverse agonist ICI 118,551 can induce water channel opening during the conformational transition of β2AR. Markov state model (MSM) analysis proves that the energy contour can be divided into seven states. States S1, S2 and S5, which represent the active conformation of β2AR, show that the water channel is in the closed state, while states S4 and S6, which correspond to the intermediate state conformation of β2AR, indicate the water channel opens gradually. State S7, which represents the inactive structure of β2AR, corresponds to the full open state of the water channel. The opening mechanism of the water channel is involved in the ligand-induced conformational change of β2AR. These results can provide useful information for understanding the opening mechanism of the water channel and will be useful for the rational design of potent inverse agonists of β2AR.

Search for β2 adrenergic receptor ligands by virtual screening via grid computing and investigation of binding modes by docking and molecular dynamics simulations

We designed a program called MolGridCal that can be used to screen small molecule database in grid computing on basis of JPPF grid environment. Based on MolGridCal program, we proposed an integrated strategy for virtual screening and binding mode investigation by combining molecular docking, molecular dynamics (MD) simulations and free energy calculations. To test the effectiveness of MolGridCal, we screened potential ligands for β2 adrenergic receptor (β2AR) from a database containing 50,000 small molecules. MolGridCal can not only send tasks to the grid server automatically, but also can distribute tasks using the screensaver function. As for the results of virtual screening, the known agonist BI-167107 of β2AR is ranked among the top 2% of the screened candidates, indicating MolGridCal program can give reasonable results. To further study the binding mode and refine the results of MolGridCal, more accurate docking and scoring methods are used to estimate the binding affinity for the top three molecules (agonist BI-167107, neutral antagonist alprenolol and inverse agonist ICI 118,551). The results indicate agonist BI-167107 has the best binding affinity. MD simulation and free energy calculation are employed to investigate the dynamic interaction mechanism between the ligands and β2AR. The results show that the agonist BI-167107 also has the lowest binding free energy. This study can provide a new way to perform virtual screening effectively through integrating molecular docking based on grid computing, MD simulations and free energy calculations. The source codes of MolGridCal are freely available at http://molgridcal.codeplex.com.

Molecular modeling study on the dynamical structural features of human smoothened receptor and binding mechanism of antagonist LY2940680 by metadynamics simulation and free energy calculation

BACKGROUND

The smoothened (SMO) receptor, one of the Class F G protein coupled receptors (GPCRs), is an essential component of the canonical hedgehog signaling pathway which plays a key role in the regulation of embryonic development in animals. The function of the SMO receptor can be modulated by small-molecule agonists and antagonists, some of which are potential antitumour agents. Understanding the binding mode of an antagonist in the SMO receptor is crucial for the rational design of new antitumour agents.

METHODS

Molecular dynamics (MD) simulation and dynamical network analysis are used to study the dynamical structural features of SMO receptor. Metadynamics simulation and free energy calculation are employed to explore the binding mechanism between the antagonist and SMO receptor.

RESULTS

The MD simulation results and dynamical network analysis show that the conserved KTXXXW motif in helix VIII has strong interaction with helix I. The α-helical extension of transmembrane 6 (TM6) is detected as part of the ligand-binding pocket and dissociation pathway of the antagonist. The metadynamics simulation results illustrate the binding mechanism of the antagonist in the pocket of SMO receptor, and free energy calculation shows the antagonist needs to overcome about 38kcal/mol of energy barrier to leave the binding pocket of SMO receptor.

CONCLUSIONS

The unusually long TM6 plays an important role on the binding behavior of the antagonist in the pocket of SMO receptor.

GENERAL SIGNIFICANCE

The results can not only profile the binding mechanism between the antagonist and Class F GPCRs, but also supply the useful information for the rational design of a more potential small molecule antagonist bound to SMO receptor.

Theoretical insights into catalytic mechanism of protein arginine methyltransferase 1

Protein arginine methyltransferase 1 (PRMT1), the major arginine asymmetric dimethylation enzyme in mammals, is emerging as a potential drug target for cancer and cardiovascular disease. Understanding the catalytic mechanism of PRMT1 will facilitate inhibitor design. However, detailed mechanisms of the methyl transfer process and substrate deprotonation of PRMT1 remain unclear. In this study, we present a theoretical study on PRMT1 catalyzed arginine dimethylation by employing molecular dynamics (MD) simulation and quantum mechanics/molecular mechanics (QM/MM) calculation. Ternary complex models, composed of PRMT1, peptide substrate, and S-adenosyl-methionine (AdoMet) as cofactor, were constructed and verified by 30-ns MD simulation. The snapshots selected from the MD trajectory were applied for the QM/MM calculation. The typical SN2-favored transition states of the first and second methyl transfers were identified from the potential energy profile. Deprotonation of substrate arginine occurs immediately after methyl transfer, and the carboxylate group of E144 acts as proton acceptor. Furthermore, natural bond orbital analysis and electrostatic potential calculation showed that E144 facilitates the charge redistribution during the reaction and reduces the energy barrier. In this study, we propose the detailed mechanism of PRMT1-catalyzed asymmetric dimethylation, which increases insight on the small-molecule effectors design, and enables further investigations into the physiological function of this family.

Computational study on the different ligands induced conformation change of β2 adrenergic receptor-Gs protein complex

β₂ adrenergic receptor (β₂AR) regulated many key physiological processes by activation of a heterotrimeric GTP binding protein (Gs protein). This process could be modulated by different types of ligands. But the details about this modulation process were still not depicted. Here, we performed molecular dynamics (MD) simulations on the structures of β₂AR-Gs protein in complex with different types of ligands. The simulation results demonstrated that the agonist BI-167107 could form hydrogen bonds with Ser203^5.42, Ser207^5.46 and Asn293^6.55 more than the inverse agonist ICI 118,551. The different binding modes of ligands further affected the conformation of β₂AR. The energy landscape profiled the energy contour map of the stable and dissociated conformation of Gαs and Gβγ when different types of ligands bound to β₂AR. It also showed the minimum energy pathway about the conformational change of Gαs and Gβγ along the reaction coordinates. By using interactive essential dynamics analysis, we found that Gαs and Gβγ domain of Gs protein had the tendency to separate when the inverse agonist ICI 118,551 bound to β₂AR. The α5-helix had a relatively quick movement with respect to transmembrane segments of β₂AR when the inverse agonist ICI 118,551 bound to β₂AR. Besides, the analysis of the centroid distance of Gαs and Gβγ showed that the Gαs was separated from Gβγ during the MD simulations. Our results not only could provide details about the different types of ligands that induced conformational change of β₂AR and Gs protein, but also supplied more information for different efficacies of drug design of β₂AR.

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.

QM/MM MD and free energy simulations of G9a-like protein (GLP) and its mutants: understanding the factors that determine the product specificity

Certain lysine residues on histone tails could be methylated by protein lysine methyltransferases (PKMTs) using S-adenosyl-L-methionine (AdoMet) as the methyl donor. Since the methylation states of the target lysines play a fundamental role in the regulation of chromatin structure and gene expression, it is important to study the property of PKMTs that allows a specific number of methyl groups (one, two or three) to be added (termed as product specificity). It has been shown that the product specificity of PKMTs may be controlled in part by the existence of specific residues at the active site. One of the best examples is a Phe/Tyr switch found in many PKMTs. Here quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) and free energy simulations are performed on wild type G9a-like protein (GLP) and its F1209Y and Y1124F mutants for understanding the energetic origin of the product specificity and the reasons for the change of product specificity as a result of single-residue mutations at the Phe/Tyr switch as well as other positions. The free energy barriers of the methyl transfer processes calculated from our simulations are consistent with experimental data, supporting the suggestion that the relative free energy barriers may determine, at least in part, the product specificity of PKMTs. The changes of the free energy barriers as a result of the mutations are also discussed based on the structural information obtained from the simulations. The results suggest that the space and active-site interactions around the ε-amino group of the target lysine available for methyl addition appear to among the key structural factors in controlling the product specificity and activity of PKMTs.