Dataset Information

Efficient Atomistic Simulation of Pathways and Calculation of Rate Constants for a Protein-Peptide Binding Process: Application to the MDM2 Protein and an Intrinsically Disordered p53 Peptide.

ABSTRACT: The characterization of protein binding processes - with all of the key conformational changes - has been a grand challenge in the field of biophysics. Here, we have used the weighted ensemble path sampling strategy to orchestrate molecular dynamics simulations, yielding atomistic views of protein-peptide binding pathways involving the MDM2 oncoprotein and an intrinsically disordered p53 peptide. A total of 182 independent, continuous binding pathways were generated, yielding a kon that is in good agreement with experiment. These pathways were generated in 15 days using 3500 cores of a supercomputer, substantially faster than would be possible with "brute force" simulations. Many of these pathways involve the anchoring of p53 residue F19 into the MDM2 binding cleft when forming the metastable encounter complex, indicating that F19 may be a kinetically important residue. Our study demonstrates that it is now practical to generate pathways and calculate rate constants for protein binding processes using atomistic simulation on typical computing resources.

SUBMITTER: Zwier MC

PROVIDER: S-EPMC5008990 | biostudies-literature | 2016 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Efficient Atomistic Simulation of Pathways and Calculation of Rate Constants for a Protein-Peptide Binding Process: Application to the MDM2 Protein and an Intrinsically Disordered p53 Peptide.

Zwier Matthew C MC Pratt Adam J AJ Adelman Joshua L JL Kaus Joseph W JW Zuckerman Daniel M DM Chong Lillian T LT

The journal of physical chemistry letters 20160822 17

The characterization of protein binding processes - with all of the key conformational changes - has been a grand challenge in the field of biophysics. Here, we have used the weighted ensemble path sampling strategy to orchestrate molecular dynamics simulations, yielding atomistic views of protein-peptide binding pathways involving the MDM2 oncoprotein and an intrinsically disordered p53 peptide. A total of 182 independent, continuous binding pathways were generated, yielding a kon that is in go ...[more]

PMID: 27532687

Similar Datasets

Project description:Disordered regions and Intrinsically Disordered Proteins (IDPs) are involved in critical cellular processes and may acquire a stable three-dimensional structure only upon binding to their partners. IDPs may follow a folding-after-binding process, known as induced folding, or a folding-before-binding process, known as conformational selection. The transcription factor p53 is involved in the regulation of cellular events that arise upon stress or DNA damage. The p53 domain structure is composed of an N-terminal transactivation domain (p53TAD), a DNA Binding Domain and a tetramerization domain. The activity of TAD is tightly regulated by interactions with cofactors, inhibitors and phosphorylation. To initiate transcription, p53TAD binds to the TAZ2 domain of CBP, a co-transcription factor, and undergoes a folding and binding process, as revealed by the recent NMR structure of the complex. The activity of p53 is regulated by phosphorylation at multiple sites on the TAD domain and recent studies have shown that modifications at three residues affect the binding towards TAZ2. However, we still do not know how these phosphorylations affect the structure of the bound state and, therefore, how they regulate the p53 function. In this work, we have used computational simulations to understand how phosphorylation affects the structure of the p53TAD:TAZ2 complex and regulates the recognition mechanism. Phosphorylation has been proposed to enhance binding by direct interaction with the folded protein or by changing the unbound conformation of IDPs, for example by pre-folding the protein favoring the recognition mechanism. Here, we show an interesting turn in the p53 case: phosphorylation mainly affects the bound structure of p53TAD, highlighting the complexity of IDP protein-protein interactions. Our results are in agreement with previous experimental studies, allowing a clear picture of how p53 is regulated by phosphorylation and giving new insights into how post-translational modifications can regulate the function of IDPs.

Project description:Chemical denaturants are the most commonly used agents for unfolding proteins and are thought to act by better solvating the unfolded state. Improved solvation is expected to lead to an expansion of unfolded chains with increasing denaturant concentration, providing a sensitive probe of the denaturant action. However, experiments have so far yielded qualitatively different results concerning the effects of chemical denaturation. Studies using Förster resonance energy transfer (FRET) and other methods found an increase in radius of gyration with denaturant concentration, but most small-angle X-ray scattering (SAXS) studies found no change. This discrepancy therefore challenges our understanding of denaturation mechanism and more generally the accuracy of these experiments as applied to unfolded or disordered proteins. Here, we use all-atom molecular simulations to investigate the effect of urea and guanidinium chloride on the structure of the intrinsically disordered protein ACTR, which can be studied by experiment over a wide range of denaturant concentration. Using unbiased molecular simulations with a carefully calibrated denaturant model, we find that the protein chain indeed swells with increasing denaturant concentration. This is due to the favorable association of urea or guanidinium chloride with the backbone of all residues and with the side-chains of almost all residues, with denaturant-water transfer free energies inferred from this association in reasonable accord with experimental estimates. Interactions of the denaturants with the backbone are dominated by hydrogen bonding, while interactions with side-chains include other contributions. By computing FRET efficiencies and SAXS intensities at each denaturant concentration, we show that the simulation trajectories are in accord with both experiments on this protein, demonstrating that there is no fundamental inconsistency between the two types of experiment. Agreement with experiment also supports the picture of chemical denaturation described in our simulations, driven by weak association of denaturant with the protein. Our simulations support some assumptions needed for each experiment to accurately reflect changes in protein size, namely, that the commonly used FRET chromophores do not qualitatively alter the results and that possible effects such as preferential solvent partitioning into the interior of the chain do not interfere with the determination of radius of gyration from the SAXS experiments.

Dataset Information

Efficient Atomistic Simulation of Pathways and Calculation of Rate Constants for a Protein-Peptide Binding Process: Application to the MDM2 Protein and an Intrinsically Disordered p53 Peptide.

Publications

Efficient Atomistic Simulation of Pathways and Calculation of Rate Constants for a Protein-Peptide Binding Process: Application to the MDM2 Protein and an Intrinsically Disordered p53 Peptide.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure