Cooperative fluctuations point to the dimerization interface of p53 core domain

Structural and Drug Targeting Insights on Mutant p53

p53 is a transcription factor with a pivotal role in cell homeostasis and fate. Its impairment is a major event in tumor onset and development. In fact, about half of human cancers bear TP53 mutations that not only halt the normal function of p53, but also may acquire oncogenic gain of functions that favor tumorigenesis. Although considered undruggable for a long time, evidence has proven the capability of many compounds to restore a wild-type (wt)-like function to mutant p53 (mutp53). However, they have not reached the clinic to date. Structural studies have strongly contributed to the knowledge about p53 structure, stability, dynamics, function, and regulation. Importantly, they have afforded relevant insights into wt and mutp53 pharmacology at molecular levels, fostering the design and development of p53-targeted anticancer therapies. Herein, we provide an integrated view of mutp53 regulation, particularly focusing on mutp53 structural traits and on targeting agents capable of its reactivation, including their biological, biochemical and biophysical features. With this, we expect to pave the way for the development of improved small molecules that may advance precision cancer therapy by targeting p53.

Modulation of Toroidal Proteins Dynamics in Favor of Functional Mechanisms upon Ligand Binding

Toroidal proteins serve as molecular machines and play crucial roles in biological processes such as DNA replication and RNA transcription. Despite progress in the structural characterization of several toroidal proteins, we still lack a mechanistic understanding of the significance of their architecture, oligomerization states, and intermolecular interactions in defining their biological function. In this work, we analyze the collective dynamics of toroidal proteins with different oligomerization states, namely, dimeric and trimeric DNA sliding clamps, nucleocapsid proteins (4-, 5-, and 6-mers) and Trp RNA-binding attenuation proteins (11- and 12-mers). We observe common global modes, among which cooperative rolling stands out as a mechanism enabling DNA processivity, and clamshell motions as those underlying the opening/closure of the sliding clamps. Alterations in global dynamics due to complexation with DNA or the clamp loader are shown to assist in enhancing motions to enable robust function. The analysis provides new insights into the differentiation and enhancement of functional motions upon intersubunit and intermolecular interactions.

Comparative Study of Elastic Network Model and Protein Contact Network for Protein Complexes: The Hemoglobin Case

The overall topology and interfacial interactions play key roles in understanding structural and functional principles of protein complexes. Elastic Network Model (ENM) and Protein Contact Network (PCN) are two widely used methods for high throughput investigation of structures and interactions within protein complexes. In this work, the comparative analysis of ENM and PCN relative to hemoglobin (Hb) was taken as case study. We examine four types of structural and dynamical paradigms, namely, conformational change between different states of Hbs, modular analysis, allosteric mechanisms studies, and interface characterization of an Hb. The comparative study shows that ENM has an advantage in studying dynamical properties and protein-protein interfaces, while PCN is better for describing protein structures quantitatively both from local and from global levels. We suggest that the integration of ENM and PCN would give a potential but powerful tool in structural systems biology.

Comparing the intrinsic dynamics of multiple protein structures using elastic network models

BACKGROUND

Elastic network models (ENMs) are based on the simple idea that a protein can be described as a set of particles connected by springs, which can then be used to describe its intrinsic flexibility using, for example, normal mode analysis. Since the introduction of the first ENM by Monique Tirion in 1996, several variants using coarser protein models have been proposed and their reliability for the description of protein intrinsic dynamics has been widely demonstrated. Lately an increasing number of studies have focused on the meaning of slow dynamics for protein function and its potential conservation through evolution. This leads naturally to comparisons of the intrinsic dynamics of multiple protein structures with varying levels of similarity.

SCOPE OF REVIEW

We describe computational strategies for calculating and comparing intrinsic dynamics of multiple proteins using elastic network models, as well as a selection of examples from the recent literature.

MAJOR CONCLUSIONS

The increasing interest for comparing dynamics across protein structures with various levels of similarity, has led to the establishment and validation of reliable computational strategies using ENMs. Comparing dynamics has been shown to be a viable way for gaining greater understanding for the mechanisms employed by proteins for their function. Choices of ENM parameters, structure alignment or similarity measures will likely influence the interpretation of the comparative analysis of protein motion.

GENERAL SIGNIFICANCE

Understanding the relation between protein function and dynamics is relevant to the fundamental understanding of protein structure-dynamics-function relationship. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Changes in dynamics upon oligomerization regulate substrate binding and allostery in amino acid kinase family members

Oligomerization is a functional requirement for many proteins. The interfacial interactions and the overall packing geometry of the individual monomers are viewed as important determinants of the thermodynamic stability and allosteric regulation of oligomers. The present study focuses on the role of the interfacial interactions and overall contact topology in the dynamic features acquired in the oligomeric state. To this aim, the collective dynamics of enzymes belonging to the amino acid kinase family both in dimeric and hexameric forms are examined by means of an elastic network model, and the softest collective motions (i.e., lowest frequency or global modes of motions) favored by the overall architecture are analyzed. Notably, the lowest-frequency modes accessible to the individual subunits in the absence of multimerization are conserved to a large extent in the oligomer, suggesting that the oligomer takes advantage of the intrinsic dynamics of the individual monomers. At the same time, oligomerization stiffens the interfacial regions of the monomers and confers new cooperative modes that exploit the rigid-body translational and rotational degrees of freedom of the intact monomers. The present study sheds light on the mechanism of cooperative inhibition of hexameric N-acetyl-L-glutamate kinase by arginine and on the allosteric regulation of UMP kinases. It also highlights the significance of the particular quaternary design in selectively determining the oligomer dynamics congruent with required ligand-binding and allosteric activities.

Protein function requires a three-dimensional structure with specific dynamic features for catalytic and binding events, and, in many cases, the structure results from the assembly of more than one polypeptide chain (also called monomer or subunit) to form an oligomer or multimer. Proteins such as hemoglobin or chaperonin GroEL are oligomers formed by 2 and 14 subunits, respectively, whereas virus capsids are multimers composed of hundreds of monomers. In these cases, the architecture of the interface between the subunits and the overall assembly geometry are essential in determining the functional motions that these sophisticated structures are able to perform under physiological conditions. Here we present results from our computational study of the large-amplitude motions of dimeric and hexameric proteins that belong to the Amino Acid Kinase family. Our study reveals that the monomers in these oligomeric proteins are arranged in such a way that the oligomer inherits the intrinsic dynamic features of its components. The packing geometry additionally confers the ability to perform highly cooperative conformational changes that involve all monomers and enable the biological activity of the multimer. The study highlights the significance of the quaternary design in favoring the oligomer dynamics that enables ligand-binding and allosteric regulation functions.

The elastic network model reveals a consistent picture on intrinsic functional dynamics of type II restriction endonucleases

The vibrational dynamics of various type II restriction endonucleases, in complex with cognate/non-cognate DNA and in the apo form, are investigated with the elastic network model in order to reveal common functional mechanisms in this enzyme family. Scissor-like and tong-like motions observed in the slowest modes of all enzymes and their complexes point to common DNA recognition and cleavage mechanisms. Normal mode analysis further points out that the scissor-like motion has an important role in differentiating between cognate and non-cognate sequences at the recognition site, thus implying its catalytic relevance. Flexible regions observed around the DNA-binding site of the enzyme usually concentrate on the highly conserved β-strands, especially after DNA binding. These β-strands may have a structurally stabilizing role in functional dynamics for target site recognition and cleavage. In addition, hot spot residues based on high-frequency modes reveal possible communication pathways between the two distant cleavage sites in the enzyme family. Some of these hot spots also exist on the shortest path between the catalytic sites and are highly conserved.

Folding of tetrameric p53: oligomerization and tumorigenic mutations induce misfolding and loss of function

The physiologically active form of p53 consists of a tetramer of four identical 393-amino-acid subunits associated via their tetramerization domains (TDs; residues 325-355). One in two human tumors contains a point mutation in the DNA binding domain (DBD) of p53 (residues 94-312). Most existing studies on the effects of these mutations on p53 structure and function have been carried out on the isolated DBD fragment, which is monomeric. Recent structural evidence, however, suggests that DBDs may interact with each other in full-length tetrameric forms of p53. Here, we investigate the effects of tumorigenic DBD mutations on the folding of p53 in its tetrameric form. We employ the construct consisting of DBD and TD (amino acids 94-360). We characterize the stability and conformational state of the tumorigenic DBD mutants R248Q, R249S, and R282Q using equilibrium denaturation and functional assays. Destabilizing mutations cause DBD to misfold when it is part of the p53 tetramer, but not when it is monomeric. This conformation is populated under moderately destabilizing conditions (10 degrees C in 2 M urea, and at physiological temperature in the absence of denaturant). Under those same conditions, it is not present in the isolated DBD fragment or in the presence of the TD mutation L344P, which abolishes tetramerization. Misfolding appears to involve intramolecular DBD-DBD association within a single tetrameric molecule. This association is promoted by destabilization of DBD (caused by mutation or elevated temperature) and by the high local DBD concentration enforced by tetramerization of TD. Disrupting the nonnative DBD-DBD interaction or transiently inhibiting tetramerization and allowing p53 to fold as a monomer may be potential strategies for pharmacological intervention in cancer.

Influence of oligomerization on the dynamics of G-protein coupled receptors as assessed by normal mode analysis

The recently discovered impact of oligomerization on G-protein coupled receptor (GPCR) function further complicates the already challenging goal of unraveling the molecular and dynamic mechanisms of these receptors. To help understand the effect of oligomerization on the dynamics of GPCRs, we have compared the motion of monomeric, dimeric, and tetrameric arrangements of the prototypic GPCR rhodopsin, using an approximate-yet powerful-normal mode analysis (NMA) technique termed elastic network model (ENM). Moreover, we have used ENM to discriminate between putative dynamic mechanisms likely to account for the recently observed conformational rearrangement of the TM4,5-TM4,5 dimerization interface of GPCRs that occurs upon activation. Our results indicate: (1) significant perturbation of the normal modes (NMs) of the rhodopsin monomer upon oligomerization, which is mainly manifested at interfacial regions; (2) increased positive correlation among the transmembrane domains (TMs) and between the extracellular loop (EL) and TM regions of the rhodopsin protomer; (3) highest interresidue positive correlation at the interfaces between protomers; and (4) experimentally testable hypotheses of differential motional changes within different putative oligomeric arrangements.

Conformational changes in protein function

Conformational changes are the hallmarks of protein dynamics and are often intimately related to protein functions. Molecular dynamics (MD) simulation is a powerful tool to study the time-resolved properties of protein structure in atomic details. In this chapter, we discuss the various applications of MD simulation to the study of protein conformational changes, and introduce several selected advanced techniques that may significantly increase the sampling efficiencies, including locally enhanced sampling (LES), and grow-to-fit molecular dynamics (G2FMD).

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.

Dynamic property is a key determinant for protein-protein interactions

Dynamic property is highly correlated with the biological functions of macromolecules, such as the activity and specificity of enzymes and the allosteric regulation in the signal transduction process. Applications of the dynamic property to protein function researches have been discussed and encouraging progresses have been achieved, for example, in enzyme activity and protein-protein docking studies. However, how the global dynamic property contributes to protein-protein interaction was still unclear. We have studied the dynamic property in protein-protein interactions based on Gaussian Network Model and applied it to classify biological and nonbiological protein-protein complexes in crystal structures. The global motion correlation between residues from the two protomers was found to be remarkably different for biological and nonbiological complexes. This correlation has been used to discriminate biological and nonbiological complexes in crystal and gave a classification rate of 86.9% in the cross-validation test. The innovation of this feature is that it is a global dynamic property which does not rely directly on the interfacial properties of the complex. In addition, the correlation of the global motions was found to be weakly correlated with the dissociation rate constant of protein complexes. We suggest that the dynamic property is a key determinant for protein-protein interaction, which can be used to discriminate native and crystal complexes and potentially be applied in protein-protein dynamic rate constants estimations.