Structure of the p53 transactivation domain in complex with the nuclear receptor coactivator binding domain of CREB binding protein

The activity and stability of the tumor suppressor p53 are regulated by interactions with key cellular proteins such as MDM2 and CBP/p300. The transactivation domain (TAD) of p53 contains two subdomains (AD1 and AD2) and interacts directly with the N-terminal domain of MDM2 and with several domains of CBP/p300. Here we report the NMR structure of the full-length p53 TAD in complex with the nuclear coactivator binding domain (NCBD) of CBP. Both the p53 TAD and NCBD are intrinsically disordered and fold synergistically upon binding, as evidenced by the observed increase in helicity and increased level of dispersion of the amide proton resonances. The p53 TAD folds to form a pair of helices (denoted Pα1 and Pα2), which extend from Phe19 to Leu25 and from Pro47 to Trp53, respectively. In the complex, the NCBD forms a bundle of three helices (Cα1, residues 2066-2075; Cα2, residues 2081-2092; and Cα3, residues 2095-2105) with a hydrophobic groove into which p53 helices Pα1 and Pα2 dock. The polypeptide chain between the p53 helices remains flexible and makes no detectable intermolecular contacts with the NCBD. Complex formation is driven largely by hydrophobic contacts that form a stable intermolecular hydrophobic core. A salt bridge between D49 of p53 and R2105 of NCBD may contribute to the binding specificity. The structure provides the first insights into simultaneous binding of the AD1 and AD2 motifs to a target protein.     

Residual Structure Accelerates Binding of Intrinsically Disordered ACTR by Promoting Efficient Folding upon Encounter

Intrinsically disordered proteins (IDPs) often fold into stable structures upon specific binding. The roles of residual structure of unbound IDPs in coupling binding and folding have been under much debate. While many studies emphasize the importance of conformational flexibility for IDP recognition, it was recently demonstrated that stabilization the N-terminal helix of intrinsically disordered ACTR accelerated its binding to another IDP, NCBD of the CREB-binding protein. To understand how enhancing ACTR helicity accelerates binding, we derived a series of topology-based coarse-grained models that mimicked various ACTR mutants with increasing helical contents and reproduced their NCBD binding affinities. Molecular dynamics simulations were then performed to sample hundreds of reversible coupled binding and folding transitions. The results show that increasing ACTR helicity does not alter the baseline mechanism of synergistic folding, which continues to follow “extended conformational selection” with multiple stages of selection and induced folding. Importantly, these coarse-grained models, while only calibrated based on binding thermodynamics, recapitulate the observed kinetic acceleration with increasing ACTR helicity. However, the residual helices do not enhance the association kinetics via more efficient seeding of productive collisions. Instead, they allow the nonspecific collision complexes to evolve more efficiently into the final bound and folded state, which is the primary source of accelerated association kinetics. Meanwhile, reduced dissociation kinetics with increasing ACTR helicity can be directly attributed to smaller entropic cost of forming the bound state. Altogether, this study provides important mechanistic insights into how residual structure may modulate thermodynamics and kinetics of IDP interactions.     

THz time scale structural rearrangements and binding modes in lysozyme-ligand interactions

Predicting the conformational changes in proteins that are relevant for substrate binding is an ongoing challenge in the aim of elucidating the functional states of proteins. The motions that are induced by protein-ligand interactions are governed by the protein global modes. Our measurements indicate that the detected changes in the global backbone motion of the enzyme upon binding reflect a shift from the large-scale collective dominant mode in the unbound state towards a functional twisting deformation that assists in closing the binding cleft. Correlated motion in lysozyme has been implicated in enzyme function in previous studies, but detailed characterization of the internal fluctuations that enable the protein to explore the ensemble of conformations that ultimately foster large-scale conformational change is yet unknown. For this reason, we use THz spectroscopy to investigate the picosecond time scale binding modes and collective structural rearrangements that take place in hen egg white lysozyme (HEWL) when bound by the inhibitor (NAG) ₃. These protein thermal motions correspond to fluctuations that have a role in both selecting and sampling from the available protein intrinsic conformations that communicate function. Hence, investigation of these fast, collective modes may provide knowledge about the mechanism leading to the preferred binding process in HEWL-(NAG) ₃. Specifically, in this work we find that the picosecond time scale hydrogen-bonding rearrangements taking place in the protein hydration shell with binding modify the packing density within the hydrophobic core on a local level. These localized, intramolecular contact variations within the protein core appear to facilitate the large cooperative movements within the interfacial region separating the α- and β- domain that mediate binding. The THz time-scale fluctuations identified in the protein-ligand system may also reveal a molecular mechanism for substrate recognition.     

Is a Malleable Protein Necessarily Highly Dynamic? The Hydrophobic Core of the Nuclear Coactivator Binding Domain Is Well Ordered

The nuclear coactivator binding domain of CREB binding protein folds into remarkably different structures in complex with different ligands. To understand the mechanism of the structural adaptability in the nuclear coactivator binding domain (NCBD), we have compared the dynamics of the hydrophobic core of NCBD in the ligand-free state and in a well-folded complex with the ligand activator for thyroid hormone and retinoid receptors using multiple NMR methods including methyl chemical shifts, coupling constants, and methyl order parameters. From all NMR measures, the aliphatic side chains in the hydrophobic core are slightly more dynamic in the free protein than in the complex, but have mobility comparable to the hydrophobic cores of average folded proteins. Urea titration monitored by NMR reveals that all parts of the protein, including the side-chain packing in the hydrophobic core, denatures in a single cooperative process. The molten globule characteristics of NCBD are thus restricted to a slowly fluctuating tertiary structure. Consequently, the conformational plasticity of the protein is most likely related to its low overall stability rather than an intrinsically flexible protein structure. The well-defined structure supports a model of molecular recognition dominated by conformational selection, whereas only minor structural adjustments are necessary after the association.     

A preformed binding interface in the unbound ensemble of an intrinsically disordered protein: evidence from molecular simulations

Intrinsically disordered proteins play an important role in cellular signalling, mediated by their interactions with other biomolecules. A key question concerns the nature of their binding mechanism, and whether the bound structure is induced only by proximity to the binding partner. This is difficult to answer through experiment alone because of the very heterogeneous nature of the unbound ensemble, and the probable rapid interconversion of the various unbound structures. Here we report the most extensive set of simulations on NCBD to date: we use large-scale replica exchange molecular dynamics to explore the unbound state. An important feature of the study is the use of an atomistic force field that has been parametrised against experimental data for weakly structured peptides, together with an accurate explicit water model. Neither the force field nor the starting conformations are biased towards a particular structure. The regions of NCBD that have high helical propensity in the simulations correspond closely to helices in the 'core' unbound conformation determined by NMR, although no single member of the simulated unbound ensemble closely resembles the core conformation, or either of the two known bound conformations. We have validated the results against NMR spectroscopy and SAXS measurements, obtaining reasonable agreement. The two helices which most stabilise the binding of NCBD with ACTR are formed readily; the third helix, which is less important for binding but is involved in most of the intraprotein contacts of NCBD in the bound conformation, is formed more rarely, and tends not to coexist with the other helices. These results support a mechanism by which NCBD gains the advantages of disorder, while forming binding-competent structures in the unbound state. We obtain support for this mechanism from coarse-grained simulations of NCBD with, and without, its binding partner.     

Applying conformational selection theory to improve crossdocking efficiency in 3‐phosphoinositide dependent protein kinase‐1

The emerging picture of biomolecular recognition is that of conformational selection followed by induced‐fit. Conformational selection theory states that binding partners exist in various conformations in solution, with binding involving a “selection” between complementary conformers. In this study, we devise a docking protocol that mimics conformational selection in protein–ligand binding and demonstrate that it significantly enhances crossdocking accuracy over Glide's flexible docking protocol, which is widely used in the pharmaceutical industry. Our protocol uses a pregenerated conformational ensemble to simulate ligand flexibility. The ensemble was generated by thorough conformational sampling coupled with conformer minimization. The generated conformers were then rigidly docked in the active site of the protein along with a postdocking minimization step that allows limited induced fit effects to be modeled for the ligand. We illustrate the improved performance of our protocol through crossdocking of 31 ligands to cocomplexed proteins of the kinase 3‐phosphoinositide dependent protein kinase‐1 extracted from the crystal structures 1H1W (ATP bound), 1OKY (staurosporine bound) and 3QD0 (bound to a potent inhibitor). Consistent with conformational selection theory, the performance of our protocol was the best for crossdocking to the cognate protein bound to the natural ligand, ATP. Proteins 2014; 82:436–451. © 2013 Wiley Periodicals, Inc.     

Evidence of Conformational Selection Driving the Formation of Ligand Binding Sites in Protein-Protein Interfaces

Many protein-protein interactions (PPIs) are compelling targets for drug discovery, and in a number of cases can be disrupted by small molecules. The main goal of this study is to examine the mechanism of binding site formation in the interface region of proteins that are PPI targets by comparing ligand-free and ligand-bound structures. To avoid any potential bias, we focus on ensembles of ligand-free protein conformations obtained by nuclear magnetic resonance (NMR) techniques and deposited in the Protein Data Bank, rather than on ensembles specifically generated for this study. The measures used for structure comparison are based on detecting binding hot spots, i.e., protein regions that are major contributors to the binding free energy. The main tool of the analysis is computational solvent mapping, which explores the surface of proteins by docking a large number of small “probe” molecules. Although we consider conformational ensembles obtained by NMR techniques, the analysis is independent of the method used for generating the structures. Finding the energetically most important regions, mapping can identify binding site residues using ligand-free models based on NMR data. In addition, the method selects conformations that are similar to some peptide-bound or ligand-bound structure in terms of the properties of the binding site. This agrees with the conformational selection model of molecular recognition, which assumes such pre-existing conformations. The analysis also shows the maximum level of similarity between unbound and bound states that is achieved without any influence from a ligand. Further shift toward the bound structure assumes protein-peptide or protein-ligand interactions, either selecting higher energy conformations that are not part of the NMR ensemble, or leading to induced fit. Thus, forming the sites in protein-protein interfaces that bind peptides and can be targeted by small ligands always includes conformational selection, although other recognition mechanisms may also be involved.     

How maltose influences structural changes to bind to maltose‐binding protein: Results from umbrella sampling simulation

A well‐studied periplasmic‐binding protein involved in the abstraction of maltose is maltose‐binding protein (MBP), which undergoes a ligand‐induced conformational transition from an open (ligand‐free) to a closed (ligand‐bound) state. Umbrella sampling simulations have been us to estimate the free energy of binding of maltose to MBP and to trace the potential of mean force of the unbinding event using the center‐of‐mass distance between the protein and ligand as the reaction coordinate. The free energy thus obtained compares nicely with the experimentally measured value justifying our theoretical basis. Measurement of the domain angle (N‐terminal‐domain – hinge – C‐terminal‐domain) along the unbinding pathway established the existence of three different states. Starting from a closed state, the protein shifts to an open conformation during the initial unbinding event of the ligand then resides in a semi‐open conformation and later resides predominantly in an open‐state. These transitions along the ligand unbinding pathway have been captured in greater depth using principal component analysis. It is proposed that in mixed‐model, both conformational selection and an induced‐fit mechanism combine to the ligand recognition process in MBP. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Modulation of p53 N-terminal transactivation domain 2 conformation ensemble and kinetics by phosphorylation

Abstract

Phosphorylation of protein is critical for various cell processes, which preferentially happens in intrinsically disordered proteins (IDPs). How phosphorylation modulates structural ensemble of disordered peptide remains largely unexplored. Here, using replica exchange molecular dynamics (REMD) and Markov state model (MSM), the conformational distribution and kinetics of p53 N-terminal transactivation domain (TAD) 2 as well as its dual-site phosphorylated form (pSer46, pThr55) were simulated. It reveals that the dual phosphorylation does not change overall size and secondary structure element fraction, while a change in the distribution of hydrogen bonds induces slightly more pre-existing bound helical conformations. MSM analysis indicates that the dual phosphorylation accelerates conformation exchange between disordered and order-like states in target-binding region. It suggests that p53 TAD2 after phosphorylation would be more apt to bind to both the human p62 pleckstrin homology (PH) domain and the yeast tfb1?PH domain through different binding mechanism, where experimentally it exhibits an extended and α-helix conformation, respectively, with increased binding strength in both complexes. Our study implies except binding interface, both conformation ensemble and kinetics should be considered for the effects of phosphorylation on IDPs. Abbreviations IDPs intrinsically disordered proteins

REMD replica exchange molecular dynamics

MSM Markov state model

TAD transactivation domain

PH pleckstrin homology

PRR proline-rich region

DBD DNA-binding domain

TET Tetramerization domain

REG regulatory domain

MD molecular dynamics

PME particle-mesh Ewald

TICA time-lagged independent component analysis

CK Chapman–Kolmogorov

GMRQ generalized matrix Rayleigh quotient

SARW self-avoiding random walk

KID kinase-inducible domain

MFPT mean first passage time

DSSP definition of secondary structure of proteins

RMSD root mean square deviation

R_g radius of gyration

R_ee end to end distance

Communicated by Ramaswamy H. Sarma     

Streptomycin interferes with conformational coupling between codon recognition and GTPase activation on the ribosome

Aminoacyl-tRNAs (aa-tRNAs) are selected by the ribosome through a kinetically controlled induced fit mechanism. Cognate codon recognition induces a conformational change in the decoding center and a domain closure of the 30S subunit. We studied how these global structural rearrangements are related to tRNA discrimination by using streptomycin to restrict the conformational flexibility of the 30S subunit. The antibiotic stabilized aa-tRNA on the ribosome both with a cognate and with a near-cognate codon in the A site. Streptomycin altered the rates of GTP hydrolysis by elongation factor Tu (EF-Tu) on cognate and near-cognate codons, resulting in almost identical rates of GTP hydrolysis and virtually complete loss of selectivity. These results indicate that movements within the 30S subunit at the streptomycin-binding site are essential for the coupling between base pair recognition and GTP hydrolysis, thus modulating the fidelity of aa-tRNA selection.     

A role for both conformational selection and induced fit in ligand binding by the LAO protein

Molecular recognition is determined by the structure and dynamics of both a protein and its ligand, but it is difficult to directly assess the role of each of these players. In this study, we use Markov State Models (MSMs) built from atomistic simulations to elucidate the mechanism by which the Lysine-, Arginine-, Ornithine-binding (LAO) protein binds to its ligand. We show that our model can predict the bound state, binding free energy, and association rate with reasonable accuracy and then use the model to dissect the binding mechanism. In the past, this binding event has often been assumed to occur via an induced fit mechanism because the protein's binding site is completely closed in the bound state, making it impossible for the ligand to enter the binding site after the protein has adopted the closed conformation. More complex mechanisms have also been hypothesized, but these have remained controversial. Here, we are able to directly observe roles for both the conformational selection and induced fit mechanisms in LAO binding. First, the LAO protein tends to form a partially closed encounter complex via conformational selection (that is, the apo protein can sample this state), though the induced fit mechanism can also play a role here. Then, interactions with the ligand can induce a transition to the bound state. Based on these results, we propose that MSMs built from atomistic simulations may be a powerful way of dissecting ligand-binding mechanisms and may eventually facilitate a deeper understanding of allostery as well as the prediction of new protein-ligand interactions, an important step in drug discovery.     

Molecular Mechanisms of Tight Binding through Fuzzy Interactions

Many intrinsically disordered proteins (IDPs) form fuzzy complexes upon binding to their targets. Although many IDPs are weakly bound in fuzzy complexes, some IDPs form high-affinity complexes. One example is the nonstructural protein 1 (NS1) of the 1918 Spanish influenza A virus, which hijacks cellular CRKII through the strong binding affinity (K_d ～10 nM) of its proline-rich motif (PRM^NS1) to the N-terminal Src-homology 3 domain of CRKII. However, its molecular mechanism remains elusive. Here, we examine the interplay between structural disorder of a bound PRM^NS1 and its long-range electrostatic interactions. Using x-ray crystallography and NMR spectroscopy, we found that PRM^NS1 retains substantial conformational flexibility in the bound state. Moreover, molecular dynamics simulations showed that structural disorder of the bound PRM^NS1 increases the number of electrostatic interactions and decreases the mean distances between the positively charged residues in PRM^NS1 and the acidic residues in the N-terminal Src-homology 3 domain. These results are analyzed using a polyelectrostatic model. Our results provide an insight into the molecular recognition mechanism for a high-affinity fuzzy complex.     

From Induced Fit to Conformational Selection: A Continuum of Binding Mechanism Controlled by the Timescale of Conformational Transitions

In receptor-ligand binding, a question that generated considerable interest is whether the mechanism is induced fit or conformational selection. This question is addressed here by a solvable model, in which a receptor undergoes transitions between active and inactive forms. The inactive form is favored while unbound but the active form is favored while a ligand is loosely bound. As the active-inactive transition rates increase, the binding mechanism gradually shifts from conformational selection to induced fit. The timescale of conformational transitions thus plays a crucial role in controlling binding mechanisms.     

Features of molecular recognition of intrinsically disordered proteins via coupled folding and binding

The sequence–structure–function paradigm of proteins has been revolutionized by the discovery of intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs). In contrast to traditional ordered proteins, IDPs/IDRs are unstructured under physiological conditions. The absence of well‐defined three‐dimensional structures in the free state of IDPs/IDRs is fundamental to their function. Folding upon binding is an important mode of molecular recognition for IDPs/IDRs. While great efforts have been devoted to investigating the complex structures and binding kinetics and affinities, our knowledge on the binding mechanisms of IDPs/IDRs remains very limited. Here, we review recent advances on the binding mechanisms of IDPs/IDRs. The structures and kinetic parameters of IDPs/IDRs can vary greatly, and the binding mechanisms can be highly dependent on the structural properties of IDPs/IDRs. IDPs/IDRs can employ various combinations of conformational selection and induced fit in a binding process, which can be templated by the target and/or encoded by the IDP/IDR. Further studies should provide deeper insights into the molecular recognition of IDPs/IDRs and enable the rational design of IDP/IDR binding mechanisms in the future.     

Conformational selection or induced fit? A critical appraisal of the kinetic mechanism

For almost five decades, two competing mechanisms of ligand recognition, conformational selection and induced fit, have dominated our interpretation of ligand binding in biological macromolecules. When binding-dissociation events are fast compared to conformational transitions, the rate of approach to equilibrium, k(obs), becomes diagnostic of conformational selection or induced fit based on whether it decreases or increases, respectively, with the ligand concentration, [L]. However, this simple conclusion based on the rapid equilibrium approximation is not valid in general. Here we show that conformational selection is associated with a rich repertoire of kinetic properties, with k(obs) decreasing or increasing with [L] depending on the relative magnitude of the rate of ligand dissociation, k(off), and the rate of conformational isomerization, k(r). We prove that, even for the simplest two-step mechanism of ligand binding, a decrease in k(obs) with [L] is unequivocal evidence of conformational selection, but an increase in k(obs) with [L] is not unequivocal evidence of induced fit. Ligand binding to glucokinase, thrombin, and its precursor prethrombin-2 are used as relevant examples. We conclude that conformational selection as a mechanism for a ligand binding to its target may be far more common than currently believed.     

A systematic study of the energetics involved in structural changes upon association and connectivity in protein interaction networks

The study of protein binding mechanisms is a major topic of research in structural biology. Here, we implement a combination of metrics to systematically assess the cost of backbone conformational changes that protein domains undergo upon association. Through the analyses of 2090 unique unbound → bound transitions, from over 12,000 structures, we show that two-thirds of these proteins do not suffer significant structural changes upon binding, and could thus fit the lock-and-key model well. Among the remaining proteins, one-third explores the bound conformation in the unbound state (conformational selection model) and, while most transitions are possible from an energetic perspective, a few do require external help to break the thermodynamic barrier (induced fit model). We also analyze the relationship between conformational transitions and protein connectivity, finding that, in general, domains interacting with many partners undergo smaller changes upon association, and are less likely to freely explore larger conformational changes.