Structural basis for p53 binding-induced DNA bending

Molecular dynamics simulation reveals DNA-specific recognition mechanism via c-Myb in pseudo-palindromic consensus of mim-1 promoter

This study aims to gain insight into the DNA-specific recognition mechanism of c-Myb transcription factor during the regulation of cell early differentiation and proliferation. Therefore, we chose the chicken myeloid gene, mitochondrial import protein 1 (mim-1), as a target to study the binding specificity between potential dual-Myb-binding sites. The c-Myb-binding site in mim-1 is a pseudo-palindromic sequence AACGGTT, which contains two AACNG consensuses. Simulation studies in different biological scenarios revealed that c-Myb binding with mim-1 in the forward strand (complex F) ismore stable than that inthereverse strand (complex R). The principal component analysis (PCA) dynamics trajectory analyses suggested an opening motion of the recognition helices of R2 and R3 (R2R3), resulting in the dissociation of DNA from c-Myb in complex R at 330 K, triggered by the reduced electrostatic potential on the surface of R2R3. Furthermore, the DNA confirmation and hydrogen-bond interaction analyses indicated that the major groove width of DNA increased in complex R, which affected on the hydrogen-bond formation ability between R2R3 and DNA, and directly resulted in the dissociation of DNA from R2R3. The steered molecular dynamics (SMD) simulation studies also suggested that the electrostatic potential, major groove width, and hydrogen bonds made major contribution to the DNA‍-specific recognition. In vitro trials confirmed the simulation results that c-Myb specifically bound to mim-1 in the forward strand. This study indicates that the three-dimensional (3D) structure features play an important role in the DNA-specific recognition mechanism by c-Myb besides the AACNG consensuses, which is beneficial to understanding the cell early differentiation and proliferation regulated by c-Myb, as well as the prediction of novel c-Myb-binding motifs in tumorigenesis.

Single molecule studies reveal that p53 tetramers dynamically bind response elements containing one or two half sites

The tumor suppressor protein p53 is critical for cell fate decisions, including apoptosis, senescence, and cell cycle arrest. p53 is a tetrameric transcription factor that binds DNA response elements to regulate transcription of target genes. p53 response elements consist of two decameric half-sites, and data suggest one p53 dimer in the tetramer binds to each half-site. Despite a broad literature describing p53 binding DNA, unanswered questions remain, due partly to the need for more quantitative and structural studies with full length protein. Here we describe a single molecule fluorescence system to visualize full length p53 tetramers binding DNA in real time. The data revealed a dynamic interaction in which tetrameric p53/DNA complexes assembled and disassembled without a dimer/DNA intermediate. On a wild type DNA containing two half sites, p53/DNA complexes existed in two kinetically distinct populations. p53 tetramers bound response elements containing only one half site to form a single population of complexes with reduced kinetic stability. Altering the spacing and helical phasing between two half sites affected both the population distribution of p53/DNA complexes and their kinetic stability. Our real time single molecule measurements of full length p53 tetramers binding DNA reveal the parameters that define the stability of p53/DNA complexes, and provide insight into the pathways by which those complexes assemble.

Roles of computational modelling in understanding p53 structure, biology, and its therapeutic targeting

The transcription factor p53 plays pivotal roles in numerous biological processes, including the suppression of tumours. The rich availability of biophysical data aimed at understanding its structure–function relationships since the 1990s has enabled the application of a variety of computational modelling techniques towards the establishment of mechanistic models. Together they have provided deep insights into the structure, mechanics, energetics, and dynamics of p53. In parallel, the observation that mutations in p53 or changes in its associated pathways characterize several human cancers has resulted in a race to develop therapeutic modulators of p53, some of which have entered clinical trials. This review describes how computational modelling has played key roles in understanding structural-dynamic aspects of p53, formulating hypotheses about domains that are beyond current experimental investigations, and the development of therapeutic molecules that target the p53 pathway.

p53 dynamics upon response element recognition explored by molecular simulations

p53 is a representative transcription factor that activates multiple target genes. To realize stimulus-dependent specificities, p53 has to recognize targets with structural variety, of which molecular mechanisms are largely unknown. Here, we conducted a series of long-time scale (totally more than 100-ms) coarse-grained molecular dynamics simulations, uncovering structure and dynamics of full-length p53 tetramer that recognizes its response element (RE). We obtained structures of a full-length p53 tetramer that binds to the RE, which is strikingly different from the structure of p53 at search. These structures are not only consistent with previous low-resolution or partial structural information, but also give access to previously unreachable detail, such as the preferential distribution of intrinsically disordered regions, the contacts between core domains, the DNA bending, and the connectivity of linker regions. We also explored how the RE variation affects the structure of the p53-RE complex. Further analysis of simulation trajectories revealed how p53 finds out the RE and how post-translational modifications affect the search mechanism.

Structural and sequential context of p53: A review of experimental and theoretical evidence

Approximately 27 million people are suffering from cancer that contains either an inactivating missense mutation of TP53 gene or partially abrogated p53 signaling pathway. Concerted action of folded and intrinsically disordered domains accounts for multi-faceted role of p53. The intricacy of dynamic p53 structure is believed to shed light on its cellular activity for developing new cancer therapies. In this review, insights into structural details of p53, diverse single point mutations affecting its core domain, thermodynamic understanding and therapeutic strategies for pharmacological rescue of p53 function has been illustrated. An effort has been made here to bridge the structural and sequential evidence of p53 from experimental to computational studies. First, we focused on the individual domains and the crucial protein-protein or DNA-protein contacts that determine conformation and dynamic behavior of p53. Next, the oncogenic mutations associated with cancer and its contribution to thermodynamic fluctuation has been discussed. Thus the emerging anti-cancer strategies include targeting of destabilized cancer mutants with selective inhibition of its negative regulators. Recent advances in development of small molecule inhibitors and peptides exploiting p53-MDM2 interaction has been included. In a nutshell, this review attempts to describe structural biology of p53 which provide new openings for structure-guided rescue.

Understanding the recognition mechanisms of Zα domain of human editing enzyme ADAR1 (hZα(ADAR1)) and various Z-DNAs from molecular dynamics simulation

The Z-DNA-binding domain of human double-stranded RNA adenosine deaminase I (hZαADAR1) can specifically recognize the left-handed Z-DNA which preferentially occurs at alternating purine-pyrimidine repeats, especially the CG-repeats. The interactions of hZαADAR1 and Z-DNAs in different sequence contexts can affect many important biological functions including gene regulation and chromatin remodeling. Therefore it is of great necessity to fully understand their recognition mechanisms. However, most existing studies are aimed at the standard CG-repeat Z-DNA rather than the non-CG-repeats, and whether the molecular basis of hZαADAR1 binding to various Z-DNAs are identical or not is still unclear on the atomic level. Here, based on the recently determined crystal structures of three representative non-CG-repeat Z-DNAs (d(CACGTG)2, d(CGTACG)2 and d(CGGCCG)2) in complex with hZαADAR1, 40 ns molecular dynamics simulation together with binding free energy calculation were performed for each system. For comparison, the standard CG-repeat Z-DNA (d(CGCGCG)2) complexed with hZαADAR1 was also simulated. The consistent results demonstrate that nonpolar interaction is the driving force during the protein-DNA binding process, and that polar interaction mainly from helix α3 also provides important contributions. Five common hot-spot residues were identified, namely Lys169, Lys170, Asn173, Arg174 and Tyr177. Hydrogen bond analysis coupled with surface charge distribution further reveal the interfacial information between hZαADAR1 and Z-DNA in detail. All of the analysis illustrate that four complexes share the common key features and the similar binding modes irrespective of Z-DNA sequences, suggesting that Z-DNA recognition by hZαADAR1 is conformation-specific rather than sequence-specific. Additionally, by analyzing the conformational changes of hZαADAR1, we found that the binding of Z-DNA could effectively stabilize hZαADAR1 protein. Our study can provide some valuable information for better understanding the binding mechanism between hZαADAR1 or even other Z-DNA-binding protein and Z-DNA.

Mapping the structural and dynamical features of multiple p53 DNA binding domains: insights into loop 1 intrinsic dynamics

The transcription factor p53 regulates cellular integrity in response to stress. p53 is mutated in more than half of cancerous cells, with a majority of the mutations localized to the DNA binding domain (DBD). In order to map the structural and dynamical features of the DBD, we carried out multiple copy molecular dynamics simulations (totaling 0.8 μs). Simulations show the loop 1 to be the most dynamic element among the DNA-contacting loops (loops 1-3). Loop 1 occupies two major conformational states: extended and recessed; the former but not the latter displays correlations in atomic fluctuations with those of loop 2 (~24 Å apart). Since loop 1 binds to the major groove whereas loop 2 binds to the minor groove of DNA, our results begin to provide some insight into the possible mechanism underpinning the cooperative nature of DBD binding to DNA. We propose (1) a novel mechanism underlying the dynamics of loop 1 and the possible tread-milling of p53 on DNA and (2) possible mutations on loop 1 residues to restore the transcriptional activity of an oncogenic mutation at a distant site.

The role of response elements organization in transcription factor selectivity: the IFN-β enhanceosome example

What is the mechanism through which transcription factors (TFs) assemble specifically along the enhancer DNA? The IFN-β enhanceosome provides a good model system: it is small; its components' crystal structures are available; and there are biochemical and cellular data. In the IFN-β enhanceosome, there are few protein-protein interactions even though consecutive DNA response elements (REs) overlap. Our molecular dynamics (MD) simulations on different motif combinations from the enhanceosome illustrate that cooperativity is achieved via unique organization of the REs: specific binding of one TF can enhance the binding of another TF to a neighboring RE and restrict others, through overlap of REs; the order of the REs can determine which complexes will form; and the alternation of consensus and non-consensus REs can regulate binding specificity by optimizing the interactions among partners. Our observations offer an explanation of how specificity and cooperativity can be attained despite the limited interactions between neighboring TFs on the enhancer DNA. To date, when addressing selective TF binding, attention has largely focused on RE sequences. Yet, the order of the REs on the DNA and the length of the spacers between them can be a key factor in specific combinatorial assembly of the TFs on the enhancer and thus in function. Our results emphasize cooperativity via RE binding sites organization.

Predicting the coordination geometry for Mg2+ in the p53 DNA-binding domain: insights from computational studies

Zn(2+) in the tumor-suppressor protein p53 DNA-binding domain (DBD) is essential for its structural stability and DNA-binding specificity. Mg(2+) has also been recently reported to bind to the p53DBD and influence its DNA-binding activity. In this contribution, the binding geometry of Mg(2+) in the p53DBD and the mechanism of how Mg(2+) affects its DNA-binding activity were investigated using density functional theory (DFT) calculations and molecular dynamics (MD) simulations. Various possible coordination geometries of Mg(2+) binding to histidines (His), cysteines (Cys), and water molecules were studied at the B3LYP/6-311+g** level of theory. The protonation state of Cys and the environment were taken into account to explore the factors governing the coordination geometry. The free energy of the reaction to form the Mg(2+) complexes was estimated, suggesting that the favorable binding mode changes from a four- to six-coordinated geometry as the number of the protonated Cys increases. Furthermore, MD simulations were employed to explore the binding modes of Mg(2+) in the active site of the p53DBD. The simulation results of the Mg(2+) system and the native Zn(2+) system show that the binding affinity of Mg(2+)to the p53DBD is weaker than that of Zn(2+), in agreement with the DFT calculation results and experiments. In addition, the two metal ions are found to make a significant contribution to maintain a favorable orientation for Arg248 to interact with putative DNA, which is critically important to the sequence-specific DNA-binding activity of the p53DBD. However, the effect of Mg(2+) is less marked. Additionally, analysis of the natural bond orbital (NBO) charge transfer reveals that Mg(2+) has a higher net positive charge than Zn(2+), leading to a stronger electrostatic attractive interaction between Mg(2+) and putative DNA. This may partly explain the higher sequence-independent DNA-binding affinity of p53DBD-Mg(2+) compared to p53DBD-Zn(2+) observed in experiment.

Lysine120 interactions with p53 response elements can allosterically direct p53 organization

p53 can serve as a paradigm in studies aiming to figure out how allosteric perturbations in transcription factors (TFs) triggered by small changes in DNA response element (RE) sequences, can spell selectivity in co-factor recruitment. p53-REs are 20-base pair (bp) DNA segments specifying diverse functions. They may be located near the transcription start sites or thousands of bps away in the genome. Their number has been estimated to be in the thousands, and they all share a common motif. A key question is then how does the p53 protein recognize a particular p53-RE sequence among all the similar ones? Here, representative p53-REs regulating diverse functions including cell cycle arrest, DNA repair, and apoptosis were simulated in explicit solvent. Among the major interactions between p53 and its REs involving Lys120, Arg280 and Arg248, the bps interacting with Lys120 vary while the interacting partners of other residues are less so. We observe that each p53-RE quarter site sequence has a unique pattern of interactions with p53 Lys120. The allosteric, DNA sequence-induced conformational and dynamic changes of the altered Lys120 interactions are amplified by the perturbation of other p53-DNA interactions. The combined subtle RE sequence-specific allosteric effects propagate in the p53 and in the DNA. The resulting amplified allosteric effects far away are reflected in changes in the overall p53 organization and in the p53 surface topology and residue fluctuations which play key roles in selective co-factor recruitment. As such, these observations suggest how similar p53-RE sequences can spell the preferred co-factor binding, which is the key to the selective gene transactivation and consequently different functional effects.

Preferred drifting along the DNA major groove and cooperative anchoring of the p53 core domain: mechanisms and scenarios

While the importance of specific p53-DNA binding is broadly accepted, the recognition process is still not fully understood. Figuring out the initial tetrameric p53-DNA association and the swift and cooperative search for specific binding sites is crucial for understanding the transactivation mechanism and selectivity. To gain insight into the p53-DNA binding process, here we have carried out explicit solvent molecular dynamic (MD) simulations of several p53 core domain-DNA conformations with the p53 and the DNA separated by varying distances. p53 approached the DNA, bound non-specifically, and quickly drifted along the DNA surface to find the major groove, cooperatively anchoring in a way similar to the specific binding observed in the crystal structure. Electrostatics was the major driving force behind the p53 movement. Mechanistically, this is a cooperative process: key residues, particularly Lys120 and Arg280 acted as sensors; upon finding their hydrogen-bonding partners, they lock in, anchoring p53 into the major groove. Concomitantly, the DNA adopted a conformation that facilitated p53 easy access. The initial non-specific core domain-DNA contacts assist in shifting the DNA and the p53 substrates toward conformations "ready" for specific major groove binding, with subsequent optimization of the interactions. This work is an invited contribution for the special issue of the Journal of Molecular Recognition dedicated to Professor Martin Karplus.

Crystal structure of the p53 core domain bound to a full consensus site as a self-assembled tetramer

Recent studies suggest that p53 binds predominantly to consensus sites composed of two decameric half-sites with zero spacing in vivo. Here we report the crystal structure of the p53 core domain bound to a full consensus site as a tetramer at 2.13A resolution. Comparison with previously reported structures of p53 dimer:DNA complexes and a chemically trapped p53 tetramer:DNA complex reveals that DNA binding by the p53 core domain is a cooperative self-assembling process accompanied by structural changes of the p53 dimer and DNA. Each p53 monomer interacts with its two neighboring subunits through two different protein-protein interfaces. The DNA is largely B-form and shows no discernible bend, but the central base-pairs between the two half-sites display a significant slide. The extensive protein-protein and protein-DNA interactions explain the high cooperativity and kinetic stability of p53 binding to contiguous decameric sites and the conservation of such binding-site configuration in vivo.

Molecular mechanisms of functional rescue mediated by P53 tumor suppressor mutations

We have utilized both molecular dynamics simulations and solution biophysical measurements to investigate the rescue mechanism of mutation N235K, which plays a key role in the recently identified global suppressor motif of K235/Y239/R240 in the human p53 DNA-binding domain (DBD). Previous genetic analysis indicates that N235K alone rescues five out of six destabilized cancer mutants. However, the solution biophysical measurement shows that N235K generates only a slight increase to the stability of DBD, implying a rescue mechanism that is not a simple additive contribution to thermodynamic stability. Our molecular simulations show that the N235K substitution generates two non-native salt bridges with residues D186 and E198. We find that the nonnative salt bridges, D186-K235 and E198-K235, and a native salt bridge, E171-R249, are mutually exclusive, thus resulting in only a marginal increase in stability as compared to the wild type protein. When a destabilized V157F is paired with N235K, the native salt bridge E171-R249 is retained. In this context, the non-native salt bridges, D186-K235 and E198-K235, produce a net increase in stability as compared to V157F alone. A similar rescue mechanism may explain how N235K stabilize other highly unstable beta-sandwich cancer mutants.

Cooperativity dominates the genomic organization of p53-response elements: a mechanistic view

p53-response elements (p53-REs) are organized as two repeats of a palindromic DNA segment spaced by 0 to 20 base pairs (bp). Several experiments indicate that in the vast majority of the human p53-REs there are no spacers between the two repeats; those with spacers, particularly with sizes beyond two nucleotides, are rare. This raises the question of what it indicates about the factors determining the p53-RE genomic organization. Clearly, given the double helical DNA conformation, the orientation of two p53 core domain dimers with respect to each other will vary depending on the spacer size: a small spacer of 0 to 2 bps will lead to the closest p53 dimer-dimer orientation; a 10-bp spacer will locate the p53 dimers on the same DNA face but necessitate DNA looping; while a 5-bp spacer will position the p53 dimers on opposite DNA faces. Here, via conformational analysis we show that when there are 0-2 bp spacers, p53-DNA binding is cooperative; however, cooperativity is greatly diminished when there are spacers with sizes beyond 2 bp. Cooperative binding is broadly recognized to be crucial for biological processes, including transcriptional regulation. Our results clearly indicate that cooperativity of the p53-DNA association dominates the genomic organization of the p53-REs, raising questions of the structural organization and functional roles of p53-REs with larger spacers. We further propose that a dynamic landscape scenario of p53 and p53-REs can better explain the selectivity of the degenerate p53-REs. Our conclusions bear on the evolutionary preference of the p53-RE organization and as such, are expected to have broad implications to other multimeric transcription factor response element organization.

Nuance in the double-helix and its role in protein-DNA recognition

It has been known for some time that the double-helix is not a uniform structure but rather exhibits sequence-specific variations that, combined with base-specific intermolecular interactions, offer the possibility of numerous modes of protein-DNA recognition. All-atom simulations have revealed mechanistic insights into the structural and energetic basis of various recognition mechanisms for a number of protein-DNA complexes while coarser grained simulations have begun to provide an understanding of the function of larger assemblies. Molecular simulations have also been applied to the prediction of transcription factor binding sites, while empirical approaches have been developed to predict nucleosome positioning. Studies that combine and integrate experimental, statistical and computational data offer the promise of rapid advances in our understanding of protein-DNA recognition mechanisms.

LIM and SH3 protein 1 (Lasp1) is a novel p53 transcriptional target involved in hepatocellular carcinoma

BACKGROUND/AIMS

Hepatocellular carcinoma (HCC) is one of the leading causes of cancer-related death worldwide with poor prognosis associated with tumor invasion and metastasis. The tumor suppressor p53 plays critical roles in tumor development, but there is increasing evidence for its involvement in tumor metastasis with the underlying mechanisms largely unexplored.

METHODS

Using combinatorial analysis of a p53 binding database with HCC microarray expression profile, we identified a novel metastasis-related gene Lasp1 as a potential p53 target.

RESULTS

In this study, we demonstrate that Lasp1 is indeed a bona fide p53 target by validating the functional repression effect of p53 on Lasp1 via a p53 response element. Transient transfection of wild-type p53 but not the mutant form suppressed Lasp1 in Hep3B (p53-/-) cells, while p53 siRNA up-regulated its expression in HepG2 (p53+/+) cells. p53 mutations at key residues involved in DNA binding abrogates the p53-mediated suppression of Lasp1 expression. In addition, Lasp1 regulates HCC cell growth as well as cell migration and invasion ability.

CONCLUSIONS

p53 transcriptionally represses Lasp1, which is a partner protein in affecting HCC cell motility. This suggests that p53 may play a role in influencing tumor metastasis through Lasp1.

Crystal structure of a p53 core tetramer bound to DNA

The tumor suppressor p53 regulates downstream genes in response to many cellular stresses and is frequently mutated in human cancers. Here, we report the use of a crosslinking strategy to trap a tetrameric p53 DNA binding domain (p53DBD) bound to DNA and the X-ray crystal structure of the protein/DNA complex. The structure reveals that two p53DBD dimers bind to B form DNA with no relative twist and that a p53 tetramer can bind to DNA without introducing significant DNA bending. The numerous dimer-dimer interactions involve several strictly conserved residues thus suggesting a molecular basis for p53DBD-DNA binding cooperativity. Surface residue conservation of the p53DBD tetramer bound to DNA highlights possible regions of other p53 domain or p53 cofactor interactions.

p53-Induced DNA bending: the interplay between p53-DNA and p53-p53 interactions

Specific p53 binding-induced DNA bending and its underlying responsible forces are crucial for the understanding of selective transcription activation. Diverse p53-response elements exist in the genome; however, it is not known what determines the DNA bending and to what extent. In order to gain knowledge of the forces that govern the DNA bending, molecular dynamics simulations were performed on a series of p53 core domain tetramer-DNA complexes in which each p53 core domain was bound to a DNA quarter site specifically. By varying the sequence of the central 4-base pairs of each half-site, different DNA bending extents were observed. The analysis showed that the dimer-dimer interactions in p53 were similar for the complexes; on the other hand, the specific interactions between the p53 and DNA, including the interactions of Arg280, Lys120, and Arg248 with the DNA, varied more significantly. In particular, the Arg280 interactions were better maintained in the complex with the CATG-containing DNA sequence and were mostly lost in the complex with the CTAG-containing DNA sequence. Structural analysis shows that the base pairings for the CATG sequence were stable throughout the simulation trajectory, whereas those for the CTAG sequence were partially dissociated in part of the trajectory, which affected the stability of the nearby Arg280-Gua base interactions. Thus, DNA bending depends on the balance between the p53 dimer-dimer interactions and p53-DNA interactions, which is in turn related to the DNA sequence and DNA flexibility.

Molecular dynamics simulations of nucleic acid-protein complexes

Molecular dynamics simulation studies of protein-nucleic acid complexes are more complicated than studies of either component alone-the force field has to be properly balanced, the systems tend to become very large, and a careful treatment of solvent and of electrostatic interactions is necessary. Recent investigations into several protein-DNA and protein-RNA systems have shown the feasibility of the simulation approach, yielding results of biological interest not readily accessible to experimental methods.

Probing potential binding modes of the p53 tetramer to DNA based on the symmetries encoded in p53 response elements

Symmetries in the p53 response-element (p53RE) encode binding modes for p53 tetramer to recognize DNA. We investigated the molecular mechanisms and biological implications of the possible binding modes. The probabilities evaluated with molecular dynamics simulations and DNA sequence analyses were found to be correlated, indicating that p53 tetramer models studied here are able to read DNA sequence information. The traditionally believed mode with four p53 monomers binding at all four DNA quarter-sites does not cause linear DNA to bend. Alternatively, p53 tetramer can use only two monomers to recognize DNA sequence and induce DNA bending. With an arrangement of dimer of AB dimer observed in p53 trimer-DNA complex crystal, p53 can recognize supercoiled DNA sequence-specifically by binding to quarter-sites one and four (H14 mode) and recognize Holliday junction geometry-specifically. Examining R273H mutation and p53-DNA interactions, we found that at least three R273H monomers are needed to disable the p53 tetramer, consistent with experiments. But just one R273H monomer may greatly shift the binding mode probabilities. Our work suggests that p53 needs balanced binding modes to maintain genome stability. Inverse repeat p53REs favor the H14 mode and direct repeat p53REs may have high possibilities of other modes.

The zinc/thiolate redox biochemistry of metallothionein and the control of zinc ion fluctuations in cell signaling

Free zinc ions are potent effectors of proteins. Their tightly controlled fluctuations ("zinc signals") in the picomolar range of concentrations modulate cellular signaling pathways. Sulfur (cysteine) donors generate redox-active coordination environments in proteins for the redox-inert zinc ion and make it possible for redox signals to induce zinc signals. Amplitudes of zinc signals are determined by the cellular zinc buffering capacity, which itself is redox-sensitive. In part by interfering with zinc and redox buffering, reactive species, drugs, toxins, and metal ions can elicit zinc signals that initiate physiological and pathobiochemical changes or lead to cellular injury when free zinc ions are sustained at higher concentrations. These interactions establish redox-inert zinc as an important factor in redox signaling. At the center of zinc/redox signaling are the zinc/thiolate clusters of metallothionein. They can transduce zinc and redox signals and thereby attenuate or amplify these signals.