Impact of Ser17 Phosphorylation on the Conformational Dynamics of the Oncoprotein MDM2

Deciphering the Multi-state Conformational Equilibrium of HDM2 in the Regulation of p53 Binding: Perspectives from Molecular Dynamics Simulation and NMR Analysis

HDM2 negatively regulates the activity of the tumor suppressor p53. Previous NMR studies have shown that apo-HDM2 interconverts between an "open" state in which the N-terminal "lid" is disordered and a "closed" state in which the lid covers the p53-binding site in the core region. Molecular dynamics (MD) simulation studies have been performed to elucidate the conformational dynamics of HDM2, but the direct relevance of the experimental and computational analyses is unclear. In addition, how the phosphorylation of S17 in the lid contributes to the inhibition of p53 binding remains controversial. Here, we used both NMR and MD simulations to investigate the conformational dynamics of apo-HDM2. The NMR analysis revealed that apo-HDM2 exists in a fast-exchanging equilibrium within two closed states, closed 1 and closed 2, in addition to a previously demonstrated slow-exchanging "open-closed" equilibrium. MD simulations visualized two characteristic closed states, where the spatial orientation of the key residues corresponds well to the chemical shift changes of the NMR spectra. Furthermore, the phosphorylation of S17 induced an equilibrium shift toward closed 1, thereby suppressing the binding of p53 to HDM2. This study reveals a multi-state equilibrium of apo-HDM2 and provides new insights into the regulation mechanism of HDM2-p53 interactions.

Energetics of a protein disorder-order transition in small molecule recognition

Many proteins recognise other proteins via mechanisms that involve the folding of intrinsically disordered regions upon complex formation. Here we investigate how the selectivity of a drug-like small molecule arises from its modulation of a protein disorder-to-order transition. Binding of the compound AM-7209 has been reported to confer order upon an intrinsically disordered ‘lid’ region of the oncoprotein MDM2. Calorimetric measurements revealed that truncation of the lid region of MDM2 increases the apparent dissociation constant of AM-7209 250-fold. By contrast, lid truncation has little effect on the binding of the ligand Nutlin-3a. Insights into these differential binding energetics were obtained via a complete thermodynamic analysis that featured adaptive absolute alchemical free energy of binding calculations with enhanced-sampling molecular dynamics simulations. The simulations reveal that in apo MDM2 the ordered lid state is energetically disfavoured. AM-7209, but not Nutlin-3a, shows a significant energetic preference for ordered lid conformations, thus shifting the balance towards ordering of the lid in the AM-7209/MDM2 complex. The methodology reported herein should facilitate broader targeting of intrinsically disordered regions in medicinal chemistry.

Molecular simulations and biophysical measurements elucidate why the ligand AM-7209 orders a disordered region of the protein MDM2 on binding. This work expands strategies available to medicinal chemists for targeting disordered proteins.

Uncovering the anti-angiogenic effect of semisynthetic triterpenoid CDDO-Im on HUVECs by an integrated network pharmacology approach

AIM

To reveal the molecular mechanism of anti-angiogenic activity of semisynthetic triterpenoid CDDO-Im.

MATERIALS AND METHODS

Using re-analysis of cDNA microarray data of CDDO-Im-treated human vascular endothelial cells (HUVECs) (GSE71622), functional annotation of revealed differentially expressed genes (DEGs) and analysis of their co-expression, the key processes induced by CDDO-Im in HUVECs were identified. Venn diagram analysis was further performed to reveal the common DEGs, i.e. genes both susceptible to CDDO-Im and involved in the regulation of angiogenesis. A list of probable protein targets of CDDO-Im was prepared based on Connectivity Map/cheminformatics analysis and chemical proteomics data, among which the proteins that were most associated with the angiogenesis-related regulome were identified. Finally, identified targets were validated by molecular docking and text mining approaches.

KEY FINDINGS

The effect of CDDO-Im in HUVECs can be divided into two main phases: the short early phase (0.5-3 h) with an acute FOXD1/CEBPA/JUNB-regulated pro-angiogenic response induced by xenobiotic stress, and the second anti-angiogenic step (6-24 h) with massive suppression of various angiogenesis-related processes, accompanied by the activation of cytoprotective mechanisms. Our analysis showed that the anti-angiogenic activity of CDDO-Im is mediated by its inhibition of the expression of PLAT, ETS1, A2M, SPAG9, RASGRP3, FBXO32, GCNT1 and HDGFRP3 and its direct interactions with EGFR, mTOR, NOS2, HSP90AA1, MDM2, SYK, IRF3, ATR and KIF14.

SIGNIFICANCE

Our findings provide valuable insights into the understanding of the molecular mechanisms of the anti-angiogenic activity of cyano enone-bearing triterpenoids and revealed a range of novel promising therapeutic targets to control pathological neovascularization.

Combining Virtual Reality Visualization with Ensemble Molecular Dynamics to Study Complex Protein Conformational Changes

Molecular dynamics (MD) simulations are increasingly used to elucidate relationships between protein structure, dynamics, and their biological function. Currently, it is extremely challenging to perform MD simulations of large-scale structural rearrangements in proteins that occur on millisecond timescales or beyond, as this requires very significant computational resources, or the use of cumbersome "collective variable" enhanced sampling protocols. Here, we describe a framework that combines ensemble MD simulations and virtual reality visualization (eMD-VR) to enable users to interactively generate realistic descriptions of large amplitude, millisecond timescale protein conformational changes in proteins. Detailed tests demonstrate that eMD-VR substantially decreases the computational cost of folding simulations of a WW domain, without the need to define collective variables a priori. We further show that eMD-VR generated pathways can be combined with Markov state models to describe the thermodynamics and kinetics of large-scale loop motions in the enzyme cyclophilin A. Our results suggest eMD-VR is a powerful tool for exploring protein energy landscapes in bioengineering efforts.

Mechanism of allosteric activation of human mRNA cap methyltransferase (RNMT) by RAM: insights from accelerated molecular dynamics simulations

The RNA guanine-N7 methyltransferase (RNMT) in complex with RNMT-activating miniprotein (RAM) catalyses the formation of a N7-methylated guanosine cap structure on the 5' end of nascent RNA polymerase II transcripts. The mRNA cap protects the primary transcript from exonucleases and recruits cap-binding complexes that mediate RNA processing, export and translation. By using microsecond standard and accelerated molecular dynamics simulations, we provide for the first time a detailed molecular mechanism of allosteric regulation of RNMT by RAM. We show that RAM selects the RNMT active site conformations that are optimal for binding of substrates (AdoMet and the cap), thus enhancing their affinity. Furthermore, our results strongly suggest the likely scenario in which the cap binding promotes the subsequent AdoMet binding, consistent with the previously suggested cooperative binding model. By employing the network community analyses, we revealed the underlying long-range allosteric networks and paths that are crucial for allosteric regulation by RAM. Our findings complement and explain previous experimental data on RNMT activity. Moreover, this study provides the most complete description of the cap and AdoMet binding poses and interactions within the enzyme's active site. This information is critical for the drug discovery efforts that consider RNMT as a promising anti-cancer target.

Structural Basis for Graded Inhibition of CREB:DNA Interactions by Multisite Phosphorylation

Phosphorylation of the kinase inducible domain (KID) of the cyclic AMP response element binding transcription factor (CREB) regulates its function through several mechanisms. Transcriptional activation occurs following phosphorylation at serine 133, but multisite phosphorylation in a neighboring region termed the CK cassette, residues 108-117, results in inhibition of CREB-mediated transcription. A molecular-level understanding of the mechanism of these opposing reactions has been lacking, in part because of the difficulty of preparing multiply phosphorylated CREB in vitro. By substituting a single residue, we have generated an engineered mammalian CREB in which the CK cassette can be phosphorylated in vitro by casein kinases and have characterized its interactions with cyclic AMP response element DNA. Phosphorylation of the CK cassette promotes an intramolecular interaction between the KID domain and the site of DNA binding, the basic region of the C-terminal basic leucine zipper (bZip) domain. Competition between the phosphorylated KID domain and DNA for bZip binding results in a decreased affinity of CREB for DNA. The binding free energy calculated from the dissociation constant is directly proportional to the number of phosphate groups in the CK cassette, indicating that the DNA binding is regulated by a rheostat-like mechanism. The rheostat is modulated by variation of the concentration of cations such as Mg2+ and by alternative isoforms such as the natural CREB isoform that lacks residues 162-272. Multisite phosphorylation of CREB represents a versatile mechanism by which transcription can be tuned to meet the variable needs of the cell.

How phosphorylation influences E1 subunit pyruvate dehydrogenase: A computational study

Pyruvate (PYR) dehydrogenase complex (PDC) is an enzymatic system that plays a crucial role in cellular metabolism as it controls the entry of carbon into the Krebs cycle. From a structural point of view, PDC is formed by three different subunits (E1, E2 and E3) capable of catalyzing the three reaction steps necessary for the full conversion of pyruvate to acetyl-CoA. Recent investigations pointed out the crucial role of this enzyme in the replication and survival of specific cancer cell lines, renewing the interest of the scientific community. Here, we report the results of our molecular dynamics studies on the mechanism by which posttranslational modifications, in particular the phosphorylation of three serine residues (Ser-264-α, Ser-271-α, and Ser-203-α), influence the enzymatic function of the protein. Our results support the hypothesis that the phosphorylation of Ser-264-α and Ser-271-α leads to (1) a perturbation of the catalytic site structure and dynamics and, especially in the case of Ser-264-α, to (2) a reduction in the affinity of E1 for the substrate. Additionally, an analysis of the channels connecting the external environment with the catalytic site indicates that the inhibitory effect should not be due to the occlusion of the access/egress pathways to/from the active site.

Phosphorylation promotes binding affinity of Rap-Raf complex by allosteric modulation of switch loop dynamics

The effects of phosphorylation of a serine residue on the structural and dynamic properties of Ras-like protein, Rap, and its interactions with effector protein Ras binding domain (RBD) of Raf kinase, in the presence of GTP, are investigated via molecular dynamics simulations. The simulations show that phosphorylation significantly effects the dynamics of functional loops of Rap which participate in the stability of the complex with effector proteins. The effects of phosphorylation on Rap are significant and detailed conformational analysis suggest that the Rap protein, when phosphorylated and with GTP ligand, samples different conformational space as compared to non-phosphorylated protein. In addition, phosphorylation of SER11 opens up a new cavity in the Rap protein which can be further explored for possible drug interactions. Residue network analysis shows that the phosphorylation of Rap results in a community spanning both Rap and RBD and strongly suggests transmission of allosteric effects of local alterations in Rap to distal regions of RBD, potentially affecting the downstream signalling. Binding free energy calculations suggest that phosphorylation of SER11 residue increases the binding between Rap and Raf corroborating the network analysis results. The increased binding of the Rap-Raf complex can have cascading effects along the signalling pathways where availability of Raf can influence the oncogenic effects of Ras proteins. These simulations underscore the importance of post translational modifications like phosphorylation on the functional dynamics in proteins and can be an alternative to drug-targeting, especially in notoriously undruggable oncoproteins belonging to Ras-like GTPase family.

Extended Impact of Pin1 Catalytic Loop Phosphorylation Revealed by S71E Phosphomimetic

Pin1 is a two-domain human protein that catalyzes the cis-trans isomerization of phospho-Ser/Thr-Pro (pS/T-P) motifs in numerous cell-cycle regulatory proteins. These pS/T-P motifs bind to Pin1's peptidyl-prolyl isomerase (PPIase) domain in a catalytic pocket, between an extended catalytic loop and the PPIase domain core. Previous studies showed that post-translational phosphorylation of S71 in the catalytic loop decreases substrate binding affinity and isomerase activity. To define the origins for these effects, we investigated a phosphomimetic Pin1 mutant, S71E-Pin1, using solution NMR. We find that S71E perturbs not only its host loop but also the nearby PPIase core. The perturbations identify a local network of hydrogen bonds and salt bridges that is more extended than previously thought, and includes interactions between the catalytic loop and the α2/α3 turn in the PPIase core. Explicit-solvent molecular dynamics simulations and phylogenetic analysis suggest that these interactions act as conserved "latches" between the loop and PPIase core that enhance binding of phosphorylated substrates, as they are absent in PPIases lacking pS/T-P specificity. Our results suggest that S71 is a hub residue within an electrostatic network primed for phosphorylation, and may illustrate a common mechanism of phosphorylation-mediated allostery.

Phosphorylation of the IDP KID Modulates Affinity for KIX by Increasing the Lifetime of the Complex

Intrinsically disordered proteins (IDPs) are known to undergo a range of posttranslational modifications, but by what mechanism do such modifications affect the binding of an IDP to its partner protein? We investigate this question using one such IDP, the kinase inducible domain (KID) of the transcription factor CREB, which interacts with the KIX domain of CREB-binding protein upon phosphorylation. As with many other IDPs, KID undergoes coupled folding and binding to form α-helical structure upon interacting with KIX. This single site phosphorylation plays an important role in the control of transcriptional activation in vivo. Here we show that, contrary to expectation, phosphorylation has no effect on association rates—unphosphorylated KID binds just as rapidly as pKID, the phosphorylated form—but rather, acts by increasing the lifetime of the complex. We propose that by controlling the lifetime of the bound complex of pKID:KIX via altering the dissociation rate, phosphorylation can facilitate effective control of transcription regulation.

Understanding the Phosphorylation Mechanism by Using Quantum Chemical Calculations and Molecular Dynamics Simulations

Phosphorylation is one of the most frequent post-translational modifications on proteins. It regulates many cellular processes by modulation of phosphorylation on protein structure and dynamics. However, the mechanism of phosphorylation-induced conformational changes of proteins is still poorly understood. Here, we report a computational study of three representative groups of tyrosine in ADP-ribosylhydrolase 1, serine in BTG2, and serine in Sp100C by using six molecular dynamics (MD) simulations and quantum chemical calculations. Added phosphorylation was found to disrupt hydrogen bond, and increase new weak interactions (hydrogen bond and hydrophobic interaction) during MD simulations, leading to conformational changes. Quantum chemical calculations further indicate that the phosphorylation on tyrosine, threonine, and serine could decrease the optical band gap energy (Egap), which can trigger electronic transitions to form or disrupt interactions easily. Our results provide an atomic and electronic description of how phosphorylation facilitates conformational and dynamic changes in proteins, which may be useful for studying protein function and protein design.

Markov models of the apo-MDM2 lid region reveal diffuse yet two-state binding dynamics and receptor poses for computational docking

MDM2 is a negative regulator of p53 activity and an important target for cancer therapeutics. The N-terminal lid region of MDM2 modulates interactions with p53 via competition for its binding cleft, exchanging slowly between docked and undocked conformations in the absence of p53. To better understand these dynamics, we constructed Markov State Models (MSMs) from large collections of unbiased simulation trajectories of apo-MDM2, and find strong evidence for diffuse, yet two-state folding and binding of the N-terminal region to the p53 receptor site. The MSM also identifies holo-like receptor conformations highly suitable for computational docking, despite initiating trajectories from closed-cleft receptor structures unsuitable for docking. Fixed-anchor docking studies using a test set of high-affinity small molecules and peptides show simulated receptor ensembles achieve docking successes comparable to cross-docking studies using crystal structures of receptors bound by alternative ligands. For p53, the best-scoring receptor structures have the N-terminal region lid region bound in a helical conformation mimicking the bound structure of p53, suggesting lid region association induces receptor conformations suitable for binding. These results suggest that MD + MSM approaches can sample binding-competent receptor conformations suitable for computational peptidomimetic design, and that inclusion of disordered regions may be essential to capturing the correct receptor dynamics.