beta-Hairpins,alpha-helices,and the intermediates among the secondary structures in the energy landscape of a peptide from a distal beta-hairpin of SH3 domain

Energy landscape of a peptide, extracted from a distal beta-hairpin of src SH3 domain, in explicit water was obtained with the multicanonical molecular dynamics. A variety of beta-hairpins with various strand-strand hydrogen bonds were found in the energy landscape at 300 K. There was no energy barrier between random-coil and hairpins. Thus, the peptide conformation can easily change from the random-coil to the hairpins in the thermal fluctuations at 300 K. The landscape also included two clusters of alpha-helices, among which an energy barrier existed, and besides, these helix clusters were separated from the other conformations. Thus, the free-energy barrier exists among the helices and the other conformations. Intermediate clusters were found between the helix and the hairpin clusters. The current study showed that the isolated state of this peptide in water fluctuates among random-coil, beta-hairpin, and alpha-helix. In SH3 domain, which has a topology of mainly beta-protein, the whole-protein folding may proceed when the segment is folded in the beta-hairpin and the other parts of the protein are coupled with the beta-hairpin in an energetically or kinetically favorite way.     

Unfolding of heptameric co-chaperonin protein follows "fly casting" mechanism: observation of transient nonnative heptamer

Protein recognition and binding play a fundamental role in living systems but sometimes also result in pathological aggregates. To probe the coupling between folding and binding in a homoheptameric system, we have characterized the time-resolved unfolding/disassembly mechanism of human co-chaperonin protein 10 (cpn10) by a combination of experimental and computational methods. The results from both approaches are in excellent agreement and make obvious that the kinetic process is three-state: an initial polypeptide-unfolding step, resulting in a non-native heptamer, is followed by a slower heptamer-dissociation step. This demonstrates that the barriers on the kinetic free-energy landscape are defined by thermodynamic stability. cpn10 is one of few, and the only heptameric, experimentally characterized system that follows the "fly casting scenario" of molecular recognition.     

非天然相互作用对均聚多肽链螺旋形成动力学过程的影响的计算机模拟研究

以螺旋结构的形成过程为研究对象,基于粗粒化的格子模型和动态蒙特卡罗模拟方法,初步探讨了非天然氢键相互作用对均聚多肽链螺旋折叠动力学过程的影响.研究发现,非天然氢键的引入虽然延缓了其热力学转变的发生,但也从整体上降低了折叠动力学过程的能垒,在某一特定温度之下,反而可以提高折叠速率.对其折叠路径分布的分析表明,非天然氢键可以减少慢速折叠路径的发生,而后者是导致折叠时间增加的主要因素.另一方面,比较特定温度下多肽链链构象及螺旋片断随时间的演化进程,发现非天然氢键在一定程度上影响了天然氢键的形成以及天然态构象的稳定存在,同时也加快了其部分解折叠过程.这说明,非天然相互作用的存在有利于蛋白质构象的快速动态调整,从而行使其相应的生物功能.     

Mapping Multivalency and Differential Affinities within Large Intrinsically Disordered Protein Complexes with Segmental Motion Analysis

Intrinsically disordered proteins (IDPs) can bind to multiple interaction partners. Numerous binding regions in the IDP that act in concert through complex cooperative effects facilitate such interactions, but complicate studying IDP complexes. To address this challenge we developed a combined fluorescence correlation and time‐resolved polarization spectroscopy approach to study the binding properties of the IDP nucleoporin153 (Nup153) to nuclear transport receptors (NTRs). The detection of segmental backbone mobility of Nup153 within the unperturbed complex provided a readout of local, region‐specific binding properties that are usually masked in measurements of the whole IDP. The binding affinities of functionally and structurally diverse NTRs to distinct regions of Nup153 can differ by orders of magnitudes—a result with implications for the diversity of transport routes in nucleocytoplasmic transport.     

Rate constants and mechanisms of intrinsically disordered proteins binding to structured targets

The binding of intrinsically disordered proteins (IDPs) to structured targets is gaining increasing attention. Here we review experimental and computational studies on the binding kinetics of IDPs. Experiments have yielded both the binding rate constants and the binding mechanisms, the latter via mutation and deletion studies and NMR techniques. Most computational studies have aimed at qualitative understanding of the binding rate constants or at mapping the free energy surfaces after the IDPs are engaged with their targets. The experiments and computation show that IDPs generally gain structures after they are engaged with their targets; that is, interactions with the targets facilitate the IDPs' folding. It also seems clear that the initial contact of an IDP with the target is formed by just a segment, not the entire IDP. The docking of one segment to its sub-site followed by coalescing of other segments around the corresponding sub-sites emerges as a recurring feature in the binding of IDPs. Such a dock-and-coalesce model forms the basis for quantitative calculation of binding rate constants. For both disordered and ordered proteins, strong electrostatic attraction with their targets can enhance the binding rate constants by several orders of magnitude. There are now tremendous opportunities in narrowing the gap in our understanding of IDPs relative to ordered proteins with regard to binding kinetics.     

Interplay of secondary structures and side-chain contacts in the denatured state of BBA1

The denatured state of a miniprotein BBA1 is studied under the native condition with the AMBER/Poisson-Boltzmann energy model and with the self-guided enhanced sampling technique. Forty independent trajectories are collected to sample the highly diversified denatured structures. Our simulation data show that the denatured BBA1 contains high percentage of native helix and native turn, but low percentage of native hairpin. Conditional population analysis indicates that the native helix formation and the native hairpin formation are not cooperative in the denatured state. Side-chain analysis shows that the native hydrophobic contacts are more preferred than the non-native hydrophobic contacts in the denatured BBA1. In contrast, the salt-bridge contacts are more or less nonspecific even if their populations are higher than those of hydrophobic contacts. Analysis of the trajectories shows that the native helix mostly initiates near the N terminus and propagates to the C terminus, and mostly forms from 3(10)-helix/turn to alpha helix. The same analysis shows that the native turn is important but not necessary in its formation in the denatured BBA1. In addition, the formations of the two strands in the native hairpin are rather asymmetric, demonstrating the likely influence of the protein environment. Energetic analysis shows that the native helix formation is largely driven by electrostatic interactions in denatured BBA1. Further, the native helix formation is associated with the breakup of non-native salt-bridge contacts and the accumulation of native salt-bridge contacts. However, the native hydrophobic contacts only show a small increase upon the native helix formation while the non-native hydrophobic contacts stay essentially the same, different from the evolution of hydrophobic contacts observed in an isolated helix folding.     

Side-chain interactions form late and cooperatively in the binding reaction between disordered peptides and PDZ domains

Intrinsically disordered proteins are very common and mediate numerous protein-protein and protein-DNA interactions. While it is clear that these interactions are instrumental for the life of the mammalian cell, there is a paucity of data regarding their molecular binding mechanisms. Here we have used short peptides as a model system for intrinsically disordered proteins. Linear free energy relationships based on rate and equilibrium constants for the binding of these peptides to ordered target proteins, PDZ domains, demonstrate that native side-chain interactions form mainly after the rate-limiting barrier for binding and in a cooperative fashion. This finding suggests that these disordered peptides first form a weak encounter complex with non-native interactions. The data do not support the recent notion that the affinities of intrinsically disordered proteins toward their targets are generally governed by their association rate constants. Instead, we observed the opposite for peptide-PDZ interactions, namely, that changes in K(d) correlate with changes in k(off).     

Helical Propensity in an Intrinsically Disordered Protein Accelerates Ligand Binding

Many intrinsically disordered proteins fold upon binding to other macromolecules. The secondary structure present in the well‐ordered complex is often formed transiently in the unbound state. The consequence of such transient structure for the binding process is, however, not clear. The activation domain of the activator for thyroid hormone and retinoid receptors (ACTR) is intrinsically disordered and folds upon binding to the nuclear coactivator binding domain (NCBD) of the CREB binding protein. A number of mutants was designed that selectively perturbs the amount of secondary structure in unbound ACTR without interfering with the intermolecular interactions between ACTR and NCBD. Using NMR spectroscopy and fluorescence‐monitored stopped‐flow kinetic measurements we show that the secondary structure content in helix 1 of ACTR indeed influences the binding kinetics. The results thus support the notion of preformed secondary structure as an important determinant for molecular recognition in intrinsically disordered proteins.     

Computational strategy for intrinsically disordered protein ligand design leads to the discovery of p53 transactivation domain I binding compounds that activate the p53 pathway

Intrinsically disordered proteins or intrinsically disordered regions (IDPs) have gained much attention in recent years due to their vital roles in biology and prevalence in various human diseases. Although IDPs are perceived as attractive therapeutic targets, rational drug design targeting IDPs remains challenging because of their conformational heterogeneity. Here, we propose a hierarchical computational strategy for IDP drug virtual screening (IDPDVS) and applied it in the discovery of p53 transactivation domain I (TAD1) binding compounds. IDPDVS starts from conformation sampling of the IDP target, then it combines stepwise conformational clustering with druggability evaluation to identify potential ligand binding pockets, followed by multiple docking screening runs and selection of compounds that can bind multi-conformations. p53 is an important tumor suppressor and restoration of its function provides an opportunity to inhibit cancer cell growth. TAD1 locates at the N-terminus of p53 and plays key roles in regulating p53 function. No compounds that directly bind to TAD1 have been reported due to its highly disordered structure. We successfully used IDPDVS to identify two compounds that bind p53 TAD1 and restore wild-type p53 function in cancer cells. Our study demonstrates that IDPDVS is an efficient strategy for IDP drug discovery and p53 TAD1 can be directly targeted by small molecules.

A hierarchical computational strategy for IDP drug virtual screening (IDPDVS) was proposed and successfully applied to identify compounds that bind p53 TAD1 and restore wild-type p53 function in cancer cells.     

Compensatory Adaptations of Structural Dynamics in an Intrinsically Disordered Protein Complex

Intrinsically disordered proteins (IDPs) play crucial roles in protein interaction networks and in this context frequently constitute important hubs and interfaces. Here we show by a combination of NMR and EPR spectroscopy that the binding of the cytokine osteopontin (OPN) to its natural ligand, heparin, is accompanied by thermodynamically compensating structural adaptations. The core segment of OPN expands upon binding. This “unfolding‐upon‐binding” is governed primarily through electrostatic interactions between heparin and charged patches along the protein backbone and compensates for entropic penalties due to heparin–OPN binding. It is shown how structural unfolding compensates for entropic losses through ligand binding in IDPs and elucidates the interplay between structure and thermodynamics of rapid substrate‐binding and ‐release events in IDP interaction networks.     

The Dynamic Multisite Interactions between Two Intrinsically Disordered Proteins

Protein interactions involving intrinsically disordered proteins (IDPs) comprise a variety of binding modes, from the well‐characterized folding upon binding to dynamic fuzzy complexes. To date, most studies concern the binding of an IDP to a structured protein, while the interaction between two IDPs is poorly understood. In this study, NMR, smFRET, and molecular dynamics (MD) simulation are combined to characterize the interaction between two IDPs, the C‐terminal domain (CTD) of protein 4.1G and the nuclear mitotic apparatus (NuMA) protein. It is revealed that CTD and NuMA form a fuzzy complex with remaining structural disorder. Multiple binding sites on both proteins were identified by molecular dynamics and mutagenesis studies. This study provides an atomic scenario in which two IDPs bearing multiple binding sites interact with each other in dynamic equilibrium. The combined approach employed here could be widely applicable for investigating IDPs and their dynamic interactions.     

Free-energy barriers in MbCO rebinding

The rebinding of CO to myoglobin (Mb) from locations around the active site is studied using a combination of molecular dynamics and stochastic simulations for native and L29F mutant Mb. The interaction between the dissociated ligand and the protein environment is described by the recently developed fluctuating three-point charge model for the CO molecule. Umbrella sampling along trajectories, previously found to sample the binding site (B) and the Xe4 pocket, is used to construct free-energy profiles for the ligand escape. On the basis of the Smoluchowski equation, the relaxation of different initial population distributions is followed in space and time. For native Mb at room temperature, the calculated rebinding times are in good agreement with experimental values and give an inner barrier of 4.3 kcal/mol between the docking site B (Mb...CO) and the A state (bound MbCO), compared to an effective barrier, Heff, of 4.5 kcal/mol and barriers into the majority conformation A1 and the minority conformation A3 of 2.4 and 4.3 kcal/mol, respectively. In the case of the L29F mutant, the free-energy surface is flatter and the dynamics is much more rapid. As was found in experiment, escape to the Xe4 pocket is facile for L29F whereas, for native Mb, the barriers to this site are larger. At lower temperatures, the rebinding dynamics is delayed by orders of magnitude also due to increased barriers between the docking sites.     

A computational analysis of binding modes and conformation changes of MDM2 induced by p53 and inhibitor bindings

Molecular dynamics (MD) simulations followed by principal component analysis were performed to study the conformational change of MDM2 induced by p53 and two inhibitor (P4 and MI63a) bindings. The results show that the hydrophobic cleft of MDM2 is very flexible and adaptive to different structural binding partners. The cleft tends to become wider and more stable as MDM2 binds to the three binding partners, while unbound MDM2 shows a narrower and pretty flexible cleft, which agrees with recent experimental data and theoretical studies. It was also found that the binding of P4 and p53 stabilizes the motion of the loop L2 linking the helix α2 and β strand (β3), but the presence of MI63a makes the motion of L2 disordered. In addition, the binding free energies of the three partners to MDM2 were calculated using molecular mechanics generalized Born surface area to explain the binding modes of these three partners to MDM2. This study will be helpful not only for better understanding the functional, concerted motion of MDM2, but also for the rational design of potent anticancer drugs targeting the p53–MDM2 interaction.     

The influence of variations of ligand protonation and tautomerism on protein-ligand recognition and binding energy landscape

We investigate the influence of variations of ligand protonation and tautomeric states on the protein-ligand binding energy landscape by applying the concept of structural consensus. In docking simulations, allowing full flexibility of the ligand, we explore whether the native binding mode could be successfully recovered using a non-native ligand protonation state. Here, we consider three proteins, dihydrofolate reductase, transketolase, and alpha-trichosanthin, complexed with ligands having multiple tautomeric forms. We find that for the majority of protonation and tautomeric states the native binding mode can be recovered without a great loss of accuracy.     

Time Window Expansion for HDX Analysis of an Intrinsically Disordered Protein

Application of typical HDX methods to examine intrinsically disordered proteins (IDP), proteins that are natively unstructured and highly dynamic at physiological pH, is limited because of the rapid exchange of unprotected amide hydrogens with solvent. The exchange rates of these fast exchanging amides are usually faster than the shortest time scale (10 s) employed in typical automated HDX-MS experiments. Considering the functional importance of IDPs and their association with many diseases, it is valuable to develop methods that allow the study of solution dynamics of these proteins as well as the ability to probe the interaction of IDPs with their wide range of binding partners. Here, we report the application of time window expansion to the millisecond range by altering the on-exchange pH of the HDX experiment to study a well-characterized IDP; the activation domain of the nuclear receptor coactivator, peroxisome proliferator-activated receptor gamma coactivator-1 alpha (PGC-1α). This method enabled mapping the regions of PGC-1α that are stabilized upon binding the ligand binding domain (LBD) of the nuclear receptor peroxisome proliferator-activated receptor gamma (PPARγ). We further demonstrate the method’s applicability to other binding partners of the IDP PGC-1α and pave the way for characterizing many other biologically important ID proteins.

Two‐Dimensional Crowding Uncovers a Hidden Conformation of α‐Synuclein

The intrinsically disordered protein (IDP), α‐synuclein (αS), is well‐known for phospholipid membrane binding‐coupled folding into tunable helical conformers. Here, using single‐molecule experiments in conjunction with ensemble assays and a theoretical model, we present a unique case demonstrating that the interaction–folding landscape of αS can be tuned by two‐dimensional (2D) crowding through simultaneous binding of a second protein on the bilayer surface. Unexpectedly, the experimental data show a clear deviation from a simple competitive inhibition model, but are consistent with a bimodal inhibition mechanism wherein membrane binding of a second protein (a membrane interacting chaperone, Hsp27, in this case) differentially inhibits two distinct modules of αS–membrane interaction. As a consequence, αS molecules are forced to access a hidden conformational state on the phospholipid bilayer in which only the higher‐affinity module remains membrane‐bound. Our results demonstrate that macromolecular crowding in two dimensions can play a significant role in shaping the conformational landscape of membrane‐binding IDPs with multiple binding modes.     

Folding transition-state and denatured-state ensembles of FSD-1 from folding and unfolding simulations

Characterization of the folding transition-state ensemble and the denatured-state ensemble is an important step toward a full elucidation of protein folding mechanisms. We report herein an investigation of the free-energy landscape of FSD-1 protein by a total of four sets of folding and unfolding molecular dynamics simulations with explicit solvent. The transition-state ensemble was initially identified from unfolding simulations at 500 K and was verified by simulations at 300 K starting from the ensemble structures. The denatured-state ensemble and the early-stage folding were studied by a combination of unfolding simulations at 500 K and folding simulations at 300 K starting from the extended conformation. A common feature of the transition-state ensemble was the substantial formation of the native secondary structures, including both the alpha-helix and beta-sheet, with partial exposure of the hydrophobic core in the solvent. Both the native and non-native secondary structures were observed in the denatured-state ensemble and early-stage folding, consistent with the smooth experimental melting curve. Interestingly, the contact orders of the transition-state ensemble structures were similar to that of the native structure and were notably lower than those of the compact structures found in early-stage folding, implying that chain and topological entropy might play significant roles in protein folding. Implications for FSD-1 folding mechanisms and the rate-limiting step are discussed. Analyses further revealed interesting non-native interactions in the denatured-state ensemble and early-stage folding and the possibility that destabilization of these interactions could help to enhance the stability and folding rate of the protein.