Phosphorylation Regulates the Bound Structure of an Intrinsically Disordered Protein: The p53-TAZ2 Case

Peak Scores Significantly Depend on the Relationships between Contextual Signals in ChIP-Seq Peaks

Chromatin immunoprecipitation followed by massively parallel DNA sequencing (ChIP-seq) is a central genome-wide method for in vivo analyses of DNA-protein interactions in various cellular conditions. Numerous studies have demonstrated the complex contextual organization of ChIP-seq peak sequences and the presence of binding sites for transcription factors in them. We assessed the dependence of the ChIP-seq peak score on the presence of different contextual signals in the peak sequences by analyzing these sequences from several ChIP-seq experiments using our fully enumerative GPU-based de novo motif discovery method, Argo_CUDA. Analysis revealed sets of significant IUPAC motifs corresponding to the binding sites of the target and partner transcription factors. For these ChIP-seq experiments, multiple regression models were constructed, demonstrating a significant dependence of the peak scores on the presence in the peak sequences of not only highly significant target motifs but also less significant motifs corresponding to the binding sites of the partner transcription factors. A significant correlation was shown between the presence of the target motifs FOXA2 and the partner motifs HNF4G, which found experimental confirmation in the scientific literature, demonstrating the important contribution of the partner transcription factors to the binding of the target transcription factor to DNA and, consequently, their important contribution to the peak score.

Structural impairment of p53 C-terminal due to the effect of phosphorylation and acetylation: a study on the interdependence of PTM

The C-terminal of tumor suppressor protein p53 is intrinsically disordered while unbound. This particular segment often shows structural plasticity when bound to other binding partners. The disordered component undergoes a disordered to ordered transition upon recognition. Post-translational modifications (PTMs), namely phosphorylation and acetylation, significantly alter the structural motifs of the segment. Among the various types of PTMs, phosphorylation, and acetylation of p53 at both N- and C- terminals lead to stabilization and activation. It has been noted experimentally that phosphorylation often regulates (enhances or reduces) the acetylation at specific sites. The phosphorylation of Thr377 and Ser378 reduces the acetylation of Lys373 and Lys382. Mutations of Thr377 and Ser378 to neutral Ala enhance and phospho mimic Asp reduce the acetylation of Lys373 and Lys382. Simulations of several single-point and pair-wise mutated systems have been generated to compare how the presence or absence of phosphorylation favors or disfavors the acetylation by thermodynamic and conformational analysis. We are using implicit solvent replica exchange molecular dynamics simulations to get 200 ns well-converged conformational ensembles of each system. Different sets of systems having both single and double PTMs are simulated. The results admit the appreciable change in the secondary structural level upon specific PTM. Also, the residual structure of the unbound p53 with single-point PTM varies significantly with pair-wise modifications. These observations further shed light on the relationship between the interdependencies of the specific PTM sites and the secondary structural levels.Communicated by Ramaswamy H. Sarma.

New simulation insights on the structural transition mechanism of bovine rhodopsin activation

Inactive rhodopsin can absorb photons, which induces different structural transitions that finally activate rhodopsin. We have examined the change in spatial configurations and physicochemical factors that result during the transition mechanism from the inactive to the active rhodopsin state via intermediates. During the activation process, many existing atomic contacts are disrupted, and new ones are formed. This is related to the movement of Helix 5, which tilts away from Helix 3 in the intermediate state in lumirhodopsin and moves closer to Helix 3 again in the active state. Similar patterns of changing atomic contacts are observed between Helices 3 and 5 of the adenosine and neurotensin receptors. In addition, residues 220-238 of rhodopsin, which are disordered in the inactive state, fold in the active state before binding to the Gα, where it catalyzes GDP/GTP exchange on the Gα subunit. Finally, molecular dynamics simulations in the membrane environment revealed that the arrestin binding region adopts a more flexible extended conformation upon phosphorylation, likely promoting arrestin binding and inactivation. In summary, our results provide additional structural understanding of specific rhodopsin activation which might be relevant to other Class A G protein-coupled receptor proteins.

Structural modulation of p53TAD1-TAZ2 complex upon mutations and post-translational modification

The tumour suppressing p53 is a target for genetic alterations in human cancer. Native p53, found in latent state in cells, gets activated following various intracellular or extracellular responses. It plays imperative role in cell-cycle control, via growth-arrest, DNA repair and apoptosis, mainly regulated by post-translational modifications (PTM). However, the influence of PTMs on the activity of p53 is still under extensive experimental and computational study. There are numerous PTM sites in p53, which are reported to regulate its binding affinities with other proteins. Of the many, Thr18 at transactivational domain (TAD) of p53 is reported to amplify p53 activity upon phosphorylation. To understand the molecular basis of p53 recognition by its binding partner upon mutations and PTMs, we have exploited all atom molecular dynamic (MD) simulation of p53TAD1 bound to TAZ2 domain of p300. The MD simulation inferred that phosphorylated and mutated Thr18, as a phospho-mimic, bound with TAZ2, redistributed the charge environment of the interface, thereby modulating the stronger interactions with TAZ2 to enhance the binding efficiency. The electrostatic interactions due to different charge environment together with H-bonding and hydrophobic interaction dictate diverse binding approach between the two. The results of this computational study further explain the importance of the Thr18 as a PTM site in atomistic detail, hence shedding further light to the understanding of how PTMs are imperative for p53 activity to protect the cellular world.Communicated by Ramaswamy H. Sarma.

Biophysical insights into glucose-dependent transcriptional regulation by PDX1

The pancreatic and duodenal homeobox 1 (PDX1) is a central regulator of glucose-dependent transcription of insulin in pancreatic β cells. PDX1 transcription factor activity is integral to the development and sustained health of the pancreas; accordingly, deciphering the complex network of cellular cues that lead to PDX1 activation or inactivation is an important step toward understanding the etiopathologies of pancreatic diseases and the development of novel therapeutics. Despite nearly 3 decades of research into PDX1 control of Insulin expression, the molecular mechanisms that dictate the function of PDX1 in response to glucose are still elusive. The transcriptional activation functions of PDX1 are regulated, in part, by its two intrinsically disordered regions, which pose a barrier to its structural and biophysical characterization. Indeed, many studies of PDX1 interactions, clinical mutations, and posttranslational modifications lack molecular level detail. Emerging methods for the quantitative study of intrinsically disordered regions and refined models for transactivation now enable us to validate and interrogate the biochemical and biophysical features of PDX1 that dictate its function. The goal of this review is to summarize existing PDX1 studies and, further, to generate a comprehensive resource for future studies of transcriptional control via PDX1.

Recognizing the Binding Pattern and Dissociation Pathways of the p300 Taz2-p53 TAD2 Complex

The dynamic association and dissociation between proteins are the basis of cellular signal transduction. This process becomes much more complicated if one or both interaction partners are intrinsically disordered because intrinsically disordered proteins can undergo disorder-to-order transitions upon binding to their partners. p53, a transcription factor with disordered regions, plays significant roles in many cellular signaling pathways. It is critical to understand the binding/unbinding mechanism involving these disordered regions of p53 at the residue level to reveal how p53 performs its biological functions. Here, we studied the dissociation process of the intrinsically disordered N-terminal transactivation domain 2 (TAD2) of p53 and the transcriptional adaptor zinc-binding 2 (Taz2) domain of transcriptional coactivator p300 using a combination of classical molecular dynamics, steered molecular dynamics, self-organizing maps, and time-resolved force distribution analysis (TRFDA). We observed two different dissociation pathways with different probabilities. One dissociation pathway starts from the TAD2 N-terminus and propagates to the α-helix and finally the C-terminus. The other dissociation pathway is in the opposite order. Subsequent TRFDA results reveal that key residues in TAD2 play critical roles. Besides the residues in agreement with previous experimental results, we also highlighted some other residues that play important roles in the disassociation process. In the dissociation process, non-native interactions were formed to partially compensate for the energy loss due to the breaking of surrounding native interactions. Moreover, our statistical analysis results of other experimentally determined complex structures involving either Taz2 or TAD2 suggest that the binding of the Taz2-TAD2 complex is mainly governed by the binding site of Taz2, which includes three main binding regions. Therefore, the complexes involving Taz2 may follow similar binding/unbinding behaviors, which could be studied together to generate common principles.

Changes in structure and flexibility of p53 TAD2 upon binding to p300 Taz2

p53 is a transcription factor with intrinsically disordered regions that plays an essential role in many cellular processes. As a tumor suppressor, the dysfunction of p53 causes various cancers. p53 can be activated by binding with cofactors in the cell due to stresses or DNA damages. The N-terminal transactivation domain (TAD) of p53 can regulate cell apoptosis by interacting with its binding partners, such as the transcriptional adaptor zinc-binding 2 (Taz2) domain of p300, and cofactors, i.e. cyclic-AMP response element-binding protein (CREB)-binding protein (CBP). The experimentally solved structure of p300 Taz2 and p53 TAD2 has provided insights into the interactions that potentially lead to the structural changes of p53 TAD2 and stabilize the complex. To explore the structural changes as well as the residues that lead to such changes from an isolated state to a bound state in p53 TAD2, we used all-atom molecular dynamics (MD) to simulate two different systems: (1) the p300 Taz2–p53 TAD2 complex and (2) isolated p53 TAD2 (residues 35–59). Although still largely unstructured, residues across the p53 TAD2 contribute significantly to stabilizing the binding between p300 Taz2 and p53 TAD2. Our results suggest that the binding affinity of the p300 Taz2–p53 TAD2 complex originates from hydrophobic and electrostatic interactions. The results are in agreement with previous reports and experimental data. By comparing the two simulated systems, our results not only demonstrate the structural changes of p53 TAD2 after binding with Taz2 but also identify the key residues leading to such changes. We also identify the critical residues that can provide insight into the interaction network between p300 Taz2 and p53 TAD2.

Phosphorylation of a Disordered Peptide-Structural Effects and Force Field Inconsistencies

Phosphorylation is one of the most abundant types of post-translational modifications of intrinsically disordered proteins (IDPs). This study examines the conformational changes in the 15-residue-long N-terminal fragment of the IDP statherin upon phosphorylation, using computer simulations with two different force fields: AMBER ff99SB-ILDN and CHARMM36m. The results from the simulations are compared with experimental small-angle X-ray scattering (SAXS) and circular dichroism data. In the unphosphorylated state, the two force fields are in excellent agreement regarding global structural properties such as size and shape. However, they exhibit some differences in the extent and type of the secondary structure. In the phosphorylated state, neither of the force fields performs well compared to the experimental data. Both force fields show a compaction of the peptide upon phosphorylation, greater than what is seen in SAXS experiments, although they differ in the local structure. While the CHARMM force field increases the fraction of bends in the peptide as a response to strong interactions between the phosphorylated residues and arginines, the AMBER force field shows an increase of the helical content in the N-terminal part of the peptide, where the phosphorylated residues reside, in better agreement with circular dichroism results.

Generation of the configurational ensemble of an intrinsically disordered protein from unbiased molecular dynamics simulation

Intrinsically disordered proteins (IDPs) are abundant in eukaryotic proteomes, play a major role in cell signaling, and are associated with human diseases. To understand IDP function it is critical to determine their configurational ensemble, i.e., the collection of 3-dimensional structures they adopt, and this remains an immense challenge in structural biology. Attempts to determine this ensemble computationally have been hitherto hampered by the necessity of reweighting molecular dynamics (MD) results or biasing simulation in order to match ensemble-averaged experimental observables, operations that reduce the precision of the generated model because different structural ensembles may yield the same experimental observable. Here, by employing enhanced sampling MD we reproduce the experimental small-angle neutron and X-ray scattering profiles and the NMR chemical shifts of the disordered N terminal (SH4UD) of c-Src kinase without reweighting or constraining the simulations. The unbiased simulation results reveal a weakly funneled and rugged free energy landscape of SH4UD, which gives rise to a heterogeneous ensemble of structures that cannot be described by simple polymer theory. SH4UD adopts transient helices, which are found away from known phosphorylation sites and could play a key role in the stabilization of structural regions necessary for phosphorylation. Our findings indicate that adequately sampled molecular simulations can be performed to provide accurate physical models of flexible biosystems, thus rationalizing their biological function.

AMBER and CHARMM Force Fields Inconsistently Portray the Microscopic Details of Phosphorylation

Phosphorylation of serine, threonine, and tyrosine is one of the most frequently occurring and crucial post-translational modifications of proteins often associated with important structural and functional changes. We investigated the direct effect of phosphorylation on the intrinsic conformational preferences of amino acids as a potential trigger of larger structural events. We conducted a comparative study of force fields on terminally capped amino acids (dipeptides) as the simplest model for phosphorylation. Our bias-exchange metadynamics simulations revealed that all model dipeptides sampled a great heterogeneity of ensembles affected by introduction of mono- and dianionic phosphate groups. However, the detected changes in populations of backbone conformers and side-chain rotamers did not reveal a strong discriminatory shift in preferences, as could be anticipated for the bulky, charged phosphate group. Furthermore, the AMBER and CHARMM force fields provided inconsistent populations of individual conformers as well as net structural trends upon phosphorylation. Detailed analysis of ensembles revealed competition between hydration and formation of internal hydrogen bonds involving amide hydrogens and the phosphate group. The observed difference in hydration free energy and potential for hydrogen bonding in individual force fields could be attributed to the different partial atomic charges used in each force field and, hence, the different parametrization strategies. Nevertheless, conformational propensities and net structural changes upon phosphorylation are difficult to extract from experimental measurements, and existing experimental data provide limited guidance for force field assessment and further development.

Eukaryotic transcription factors: paradigms of protein intrinsic disorder

Gene-specific transcription factors (TFs) are key regulatory components of signaling pathways, controlling, for example, cell growth, development, and stress responses. Their biological functions are determined by their molecular structures, as exemplified by their structured DNA-binding domains targeting specific cis-acting elements in genes, and by the significant lack of fixed tertiary structure in their extensive intrinsically disordered regions. Recent research in protein intrinsic disorder (ID) has changed our understanding of transcriptional activation domains from ‘negative noodles’ to ID regions with function-related, short sequence motifs and molecular recognition features with structural propensities. This review focuses on molecular aspects of TFs, which represent paradigms of ID-related features. Through specific examples, we review how the ID-associated flexibility of TFs enables them to participate in large interactomes, how they use only a few hydrophobic residues, short sequence motifs, prestructured motifs, and coupled folding and binding for their interactions with co-activators, and how their accessibility to post-translational modification affects their interactions. It is furthermore emphasized how classic biochemical concepts like allostery, conformational selection, induced fit, and feedback regulation are undergoing a revival with the appreciation of ID. The review also describes the most recent advances based on computational simulations of ID-based interaction mechanisms and structural analysis of ID in the context of full-length TFs and suggests future directions for research in TF ID.

Stimulatory effects of curcumin and quercetin on posttranslational modifications of p53 during lung carcinogenesis

Lung cancer is responsible for increase in mortality due to cancer-related deaths, and new approaches are being explored for the betterment of the situation. In the present study, chemopreventive efficacy of curcumin and quercetin was investigated against benzo(a)pyrene (BP)-induced lung carcinogenesis. The mice were segregated into five groups, which included normal control, BP-treated, BP + curcumin-treated, BP + quercetin-treated, and BP + curcumin + quercetin-treated groups. The morphological and histological analyses of tumor nodules confirmed lung carcinogenesis22 weeks after weeks single intraperitoneal injection of BP at a dose of 100 mg/kg body weight to mice. Curcumin and quercetin when administered individually as well as in combination significantly elevated the expression of acetylated-p53, which was otherwise depressed due to BP treatment. Also, both the phytochemicals significantly reduced the BP-inflicted increased levels of phosphorylated-p53. Furthermore, observed increase in the number of apoptotic cells by terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL), assay and increased activities of caspase 3 and 9 confirmed the induction of apoptosis by curcumin and quercetin. Moreover, the histological slides also showed noticeable improvement in the histoarchitecture of lungs by phytochemicals. The present study concludes that prophylactic treatment with curcumin and quercetin induces apoptosis in the lungs by modulation of p53 posttranslational modifications.

Mechanistic roles of protein disorder within transcription

Understanding the interactions of proteins involved in transcriptional regulation is critical to describing biological systems because they control the expression profile of the cell. Yet sadly they belong to a less well biophysically characterized subset of proteins; they frequently contain long disordered regions that are highly dynamic. A key question therefore is, why? What functional roles does protein disorder play in transcriptional regulation? Experimental data exemplifying these roles are starting to emerge, with common themes being enabling complexity within networks and quick responses. Most recently a role for disorder in mediating phase transitions of membrane-less organelles has been proposed.

R, C. S. Expediting dynamics approach to understand the influence of 14-3-3ζ causing metastatic cancer through the interaction of YAP1 and β-TRCP

The 14-3-3ζ protein acts as a molecular switch in regulating the TGF-β pathway, which alters from a tumor suppressor in the early stage of breast cancer to a promoter of metastasis in the late stage.

Biosimilars in rheumatic diseases: structural and functional variability that may impact clinical and regulatory decisions

Biologics as therapeutic interventions for human disease represent both a distinctly modern novelty and an echo of ancient, or at least old, medical practice. The similarity lies in the sense that in both the synthetic effort occurs in living organisms (an extract of a plant, animal tissue, or a cell culture) while the difference is apparent in the bioengineering required in modern methods and the corresponding flexibility to customize the therapeutic product. Although the concept of looking to living systems as a source of medically useful compounds either for research or for actual patient care has never vanished, the development of biochemistry and advances in medicinal chemistry made production by total synthesis the standard for a safe, reliable, and commercial drug production at sufficient scale. In this interval was where much of the modern apparatus for approving medical therapies came to be developed, and as such, the most proper extension of the regulatory regime to modern biologics is not entirely obvious. In particular, the notion of generics for off-patent conventional pharmaceuticals and their role in the marketplace with respect to increasing the accessibility of treatment is not congruent with the relationship between what are known as biosimilars and off-patent originating biologics. In this article, we review elements of the scientific basis for challenges in the production, use, and regulation of biosimilars. In light of these advances, we propose suggestions to modify constraints on biosimilar regulations in the interest of patient care and access to therapies.

Intrinsic protein disorder could be overlooked in cocrystallization conditions: An SRCD case study

X-ray diffractometry dominates protein studies, as it can provide 3D structures of these diverse macromolecules or their molecular complexes with interacting partners: substrates, inhibitors, and/or cofactors. Here, we show that under cocrystallization conditions the results could reflect induced protein folds instead of the (partially) disordered original structures. The analysis of synchrotron radiation circular dichroism spectra revealed that the Im7 immunity protein stabilizes the native-like solution structure of unfolded NColE7 nuclease mutants via complex formation. This is consistent with the fact that among the several available crystal structures with its inhibitor or substrate, all NColE7 structures are virtually the same. Our results draw attention to the possible structural consequence of protein modifications, which is often hidden by compensational effects of intermolecular interactions. The growing evidence on the importance of protein intrinsic disorder thus, demands more extensive complementary experiments in solution phase with the unligated form of the protein of interest.