Deconvoluting Protein (Un)folding Structural Ensembles Using X-Ray Scattering, Nuclear Magnetic Resonance Spectroscopy and Molecular Dynamics Simulation

The folding and unfolding of protein domains is an apparently cooperative process, but transient intermediates have been detected in some cases. Such (un)folding intermediates are challenging to investigate structurally as they are typically not long-lived and their role in the (un)folding reaction has often been questioned. One of the most well studied (un)folding pathways is that of Drosophila melanogaster Engrailed homeodomain (EnHD): this 61-residue protein forms a three helix bundle in the native state and folds via a helical intermediate. Here we used molecular dynamics simulations to derive sample conformations of EnHD in the native, intermediate, and unfolded states and selected the relevant structural clusters by comparing to small/wide angle X-ray scattering data at four different temperatures. The results are corroborated using residual dipolar couplings determined by NMR spectroscopy. Our results agree well with the previously proposed (un)folding pathway. However, they also suggest that the fully unfolded state is present at a low fraction throughout the investigated temperature interval, and that the (un)folding intermediate is highly populated at the thermal midpoint in line with the view that this intermediate can be regarded to be the denatured state under physiological conditions. Further, the combination of ensemble structural techniques with MD allows for determination of structures and populations of multiple interconverting structures in solution.

Design of a PDZbody, a bivalent binder of the E6 protein from human papillomavirus

Chronic infection by high risk human papillomavirus (HPV) strains may lead to cancer. Expression of the two viral oncoproteins E6 and E7 is largely responsible for immortalization of infected cells. The HPV E6 is a small (approximately 150 residues) two domain protein that interacts with a number of cellular proteins including the ubiquitin ligase E6-associated protein (E6AP) and several PDZ-domain containing proteins. Our aim was to design a high-affinity binder for HPV E6 by linking two of its cellular targets. First, we improved the affinity of the second PDZ domain from SAP97 for the C-terminus of HPV E6 from the high-risk strain HPV18 using phage display. Second, we added a helix from E6AP to the N-terminus of the optimized PDZ variant, creating a chimeric bivalent binder, denoted PDZbody. Full-length HPV E6 proteins are difficult to express and purify. Nevertheless, we could measure the affinity of the PDZbody for E6 from another high-risk strain, HPV16 (Kd = 65 nM). Finally, the PDZbody was used to co-immunoprecipitate E6 protein from HPV18-immortalized HeLa cells, confirming the interaction between PDZbody and HPV18 E6 in a cellular context.

Thermodynamic aspects of cAMP dependent protein kinase catalytic subunit allostery

Kinetics of thermal inactivation of acrylodan-labeled cAMP dependent protein kinase catalytic subunit, its binary complexes with ATP and peptide inhibitor PKI[5-24], respectively, and the ternary complex involving both of these ligands were studied at different temperatures (5-50 °C). The thermodynamic parameters ΔH and ΔS for ligand binding equilibria as well as for the allosteric interaction between the binding sites of these ligands were obtained by using the Van't Hoff analysis. The results indicated that more inter- and intra-molecular non-covalent bonds were involved in ATP binding with the protein when compared to the peptide binding. Similarly, nucleotide and peptide binding steps were accompanied with different entropy effects, while almost no entropy change accompanied PKI[5-24] binding, suggesting that the protein flexibility was not affected in this case. Differently from the binary complex formation the ternary complex formation was accompanied by a significant entropy change and with intensive formation of new non-covalent interactions (ΔH). At the same time both ligand binding steps as well as the allosteric interaction between ligand binding sites could be described by a common entropy-enthalpy compensation plot, pointing to a similar mechanism of these phenomena. It was concluded that numerous weak interactions govern the allostery of cAMP dependent protein kinase catalytic subunit.

Targeting protein-protein interactions with trimeric ligands: high affinity inhibitors of the MAGUK protein family

PDZ domains in general, and those of PSD-95 in particular, are emerging as promising drug targets for diseases such as ischemic stroke. We have previously shown that dimeric ligands that simultaneously target PDZ1 and PDZ2 of PSD-95 are highly potent inhibitors of PSD-95. However, PSD-95 and the related MAGUK proteins contain three consecutive PDZ domains, hence we envisioned that targeting all three PDZ domains simultaneously would lead to more potent and potentially more specific interactions with the MAGUK proteins. Here we describe the design, synthesis and characterization of a series of trimeric ligands targeting all three PDZ domains of PSD-95 and the related MAGUK proteins, PSD-93, SAP-97 and SAP-102. Using our dimeric ligands targeting the PDZ1-2 tandem as starting point, we designed novel trimeric ligands by introducing a PDZ3-binding peptide moiety via a cysteine-derivatized NPEG linker. The trimeric ligands generally displayed increased affinities compared to the dimeric ligands in fluorescence polarization binding experiments and optimized trimeric ligands showed low nanomolar inhibition towards the four MAGUK proteins, thus being the most potent inhibitors described. Kinetic experiments using stopped-flow spectrometry showed that the increase in affinity is caused by a decrease in the dissociation rate of the trimeric ligand as compared to the dimeric ligands, likely reflecting the lower probability of simultaneous dissociation of all three PDZ ligands. Thus, we have provided novel inhibitors of the MAGUK proteins with exceptionally high affinity, which can be used to further elucidate the therapeutic potential of these proteins.

Nonequilibrium capture rates induce protein accumulation and enhanced adsorption to solid-state nanopores

Single molecule capturing of analytes using an electrically biased nanopore is the fundamental mechanism in which nearly all nanopore experiments are conducted. With pore dimensions being on the order of a single molecule, the spatial zone of sensing only contains approximately a zeptoliter of volume. As a result, nanopores offer high precision sensing within the pore but provide little to no information about the analytes outside the pore. In this study, we use capture frequency and rate balance theory to predict and study the accumulation of proteins at the entrance to the pore. Protein accumulation is found to have positive attributes such as capture rate enhancement over time but can additionally lead to negative effects such as long-term blockages typically attributed to protein adsorption on the surface of the pore. Working with the folded and unfolded states of the protein domain PDZ2 from SAP97, we show that applying short (e.g., 3-25 s in duration) positive voltage pulses, rather than a constant voltage, can prevent long-term current blockades (i.e., adsorption events). By showing that the concentration of proteins around the pore can be controlled in real time using modified voltage protocols, new experiments can be explored which study the role of concentration on single molecular kinetics including protein aggregation, folding, and protein binding.

Deciphering the mechanisms of binding induced folding at nearly atomic resolution: The Φ value analysis applied to IDPs

The Φ value analysis is a method to analyze the structure of metastable states in reaction pathways. Such a methodology is based on the quantitative analysis of the effect of point mutations on the kinetics and thermodynamics of the probed reaction. The Φ value analysis is routinely used in protein folding studies and is potentially an extremely powerful tool to analyze the mechanism of binding induced folding of intrinsically disordered proteins. In this review we recapitulate the key equations and experimental advices to perform the Φ value analysis in the perspective of the possible caveats arising in intrinsically disordered systems. Finally, we briefly discuss some few examples already available in the literature.

Kinetics of acrylodan-labelled cAMP-dependent protein kinase catalytic subunit denaturation

Fluorescence spectroscopy was used to study denaturation of cAMP-dependent protein kinase catalytic subunit labeled with an acrylodan moiety. The dye was covalently bound to a cystein residue introduced into the enzyme by replacement of arginine in position 326 in the native sequence, located near the enzyme active center. This labeling had no effect on catalytic activity of the enzyme, but provided possibility to monitor changes in protein structure through measuring the fluorescence spectrum of the dye, which is sensitive to changes in its environment. This method was used to monitor denaturation of the protein kinase catalytic subunit and study the kinetics of this process as well as influence of specific ligands on stability of the protein. Stabilization of the enzyme structure was observed in the presence of adenosine triphosphate, peptide substrate RRYSV and inhibitor peptide PKI[5-24].

The role of backbone hydrogen bonds in the transition state for protein folding of a PDZ domain

Backbone hydrogen bonds are important for the structure and stability of proteins. However, since conventional site-directed mutagenesis cannot be applied to perturb the backbone, the contribution of these hydrogen bonds in protein folding and stability has been assessed only for a very limited set of small proteins. We have here investigated effects of five amide-to-ester mutations in the backbone of a PDZ domain, a 90-residue globular protein domain, to probe the influence of hydrogen bonds in a β-sheet for folding and stability. The amide-to-ester mutation removes NH-mediated hydrogen bonds and destabilizes hydrogen bonds formed by the carbonyl oxygen. The overall stability of the PDZ domain generally decreased for all amide-to-ester mutants due to an increase in the unfolding rate constant. For this particular region of the PDZ domain, it is therefore clear that native hydrogen bonds are formed after crossing of the rate-limiting barrier for folding. Moreover, three of the five amide-to-ester mutants displayed an increase in the folding rate constant suggesting that the hydrogen bonds are involved in non-native interactions in the transition state for folding.

Energetic pathway sampling in a protein interaction domain

The affinity and specificity of protein-ligand interactions are influenced by energetic crosstalk within the protein domain. However, the molecular details of such intradomain allostery are still unclear. Here, we have experimentally detected and computationally predicted interaction pathways in the postsynaptic density 95/discs large/zonula occludens 1 (PDZ)-peptide ligand model system using wild-type and circularly permuted PDZ proteins. The circular permutant introduced small perturbations in the tertiary structure and a concomitant rewiring of allosteric pathways, allowing us to describe how subtle changes may reshape energetic signaling. The results were analyzed in the context of other members of the PDZ family, which were found to contain distinct interaction pathways for different peptide ligands. The data reveal a fascinating scenario whereby several energetic pathways are sampled within one single domain and distinct pathways are activated by specific protein ligands.

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.

A frustrated binding interface for intrinsically disordered proteins

Intrinsically disordered proteins are very common in the eukaryotic proteome, and many of them are associated with diseases. Disordered proteins usually undergo a coupled binding and folding reaction and often interact with many different binding partners. Using double mutant cycles, we mapped the energy landscape of the binding interface for two interacting disordered domains and found it to be largely suboptimal in terms of interaction free energies, despite relatively high affinity. These data depict a frustrated energy landscape for interactions involving intrinsically disordered proteins, which is likely a result of their functional promiscuity.

The binding mechanisms of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins are very common and instrumental for cellular signaling. Recently, a number of studies have investigated the kinetic binding mechanisms of IDPs and IDRs. These results allow us to draw conclusions about the energy landscape for the coupled binding and folding of disordered proteins. The association rate constants of IDPs cover a wide range (10(5)-10(9) M(-1) s(-1)) and are largely governed by long-range charge-charge interactions, similarly to interactions between well-folded proteins. Off-rate constants also differ significantly among IDPs (with half-lives of up to several minutes) but are usually around 0.1-1000 s(-1), allowing for rapid dissociation of complexes. Likewise, affinities span from pM to μM suggesting that the low-affinity high-specificity concept for IDPs is not straightforward. Overall, it appears that binding precedes global folding although secondary structure elements such as helices may form before the protein-protein interaction. Short IDPs bind in apparent two-state reactions whereas larger IDPs often display complex multi-step binding reactions. While the two extreme cases of two-step binding (conformational selection and induced fit) or their combination into a square mechanism is an attractive model in theory, it is too simplistic in practice. Experiment and simulation suggest a more complex energy landscape in which IDPs bind targets through a combination of conformational selection before binding (e.g., secondary structure formation) and induced fit after binding (global folding and formation of short-range intermolecular interactions).

The transition state structure for coupled binding and folding of disordered protein domains

Intrinsically disordered proteins are abundant in the eukaryotic proteome, and they are implicated in a range of different diseases. However, there is a paucity of experimental data on molecular details of the coupled binding and folding of such proteins. Two interacting and relatively well studied disordered protein domains are the activation domain from the p160 transcriptional co-activator ACTR and the nuclear co-activator binding domain (NCBD) of CREB binding protein. We have analyzed the transition state for their coupled binding and folding by protein engineering and kinetic experiments (Φ-value analysis) and found that it involves weak native interactions between the N-terminal helices of ACTR and NCBD, but is otherwise "disordered-like". Most native hydrophobic interactions in the interface between the two domains form later, after the rate-limiting barrier for association. Linear free energy relationships suggest a cooperative formation of native interactions, reminiscent of the nucleation-condensation mechanism in protein folding.

Probing the role of backbone hydrogen bonds in protein-peptide interactions by amide-to-ester mutations

One of the most frequent protein-protein interaction modules in mammalian cells is the postsynaptic density 95/discs large/zonula occludens 1 (PDZ) domain, involved in scaffolding and signaling and emerging as an important drug target for several diseases. Like many other protein-protein interactions, those of the PDZ domain family involve formation of intermolecular hydrogen bonds: C-termini or internal linear motifs of proteins bind as β-strands to form an extended antiparallel β-sheet with the PDZ domain. Whereas extensive work has focused on the importance of the amino acid side chains of the protein ligand, the role of the backbone hydrogen bonds in the binding reaction is not known. Using amide-to-ester substitutions to perturb the backbone hydrogen-bonding pattern, we have systematically probed putative backbone hydrogen bonds between four different PDZ domains and peptides corresponding to natural protein ligands. Amide-to-ester mutations of the three C-terminal amides of the peptide ligand severely affected the affinity with the PDZ domain, demonstrating that hydrogen bonds contribute significantly to ligand binding (apparent changes in binding energy, ΔΔG = 1.3 to >3.8 kcal mol(-1)). This decrease in affinity was mainly due to an increase in the dissociation rate constant, but a significant decrease in the association rate constant was found for some amide-to-ester mutations suggesting that native hydrogen bonds have begun to form in the transition state of the binding reaction. This study provides a general framework for studying the role of backbone hydrogen bonds in protein-peptide interactions and for the first time specifically addresses these for PDZ domain-peptide interactions.

Single-molecule studies of intrinsically disordered proteins using solid-state nanopores

Partially or fully disordered proteins are instrumental for signal-transduction pathways; however, many mechanistic aspects of these proteins are not well-understood. For example, the number and nature of intermediate states along the binding pathway is still a topic of intense debate. To shed light on the conformational heterogeneity of disordered protein domains and their complexes, we performed single-molecule experiments by translocating disordered proteins through a nanopore embedded within a thin dielectric membrane. This platform allows for single-molecule statistics to be generated without the need of fluorescent labels or other modification groups. These studies were performed on two different intrinsically disordered protein domains, a binding domain from activator of thyroid hormone and retinoid receptors (ACTR) and the nuclear coactivator binding domain of CREB-binding protein (NCBD), along with their bimolecular complex. Our results demonstrate that both ACTR and NCBD populate distinct conformations upon translocation through the nanopore. The folded complex of the two disordered domains, on the other hand, translocated as one conformation. Somewhat surprisingly, we found that NCBD undergoes a charge reversal under high salt concentrations. This was verified by both translocation statistics as well as by measuring the ζ-potential. Electrostatic interactions have been previously suggested to play a key role in the association of intrinsically disordered proteins, and the observed behavior adds further complexity to their binding reactions.

Single molecule unfolding and stretching of protein domains inside a solid-state nanopore by electric field

Single molecule methods have provided a significantly new look at the behavior of biomolecules in both equilibrium and non-equilibrium conditions. Most notable are the stretching experiments performed by atomic force microscopes and laser tweezers. Here we present an alternative single molecule method that can unfold a protein domain, observed at electric fields greater than 10(6) V/m, and is fully controllable by the application of increasing voltages across the membrane of the pore. Furthermore this unfolding mechanism is characterized by measuring both the residence time of the protein within the nanopore and the current blockade. The unfolding data supports a gradual unfolding mechanism rather than the cooperative transition observed by classical urea denaturation experiments. Lastly it is shown that the voltage-mediated unfolding is a function of the stability of the protein by comparing two mutationally destabilized variants of the protein.

Tolerance of protein folding to a circular permutation in a PDZ domain

Circular permutation is a common molecular mechanism for evolution of proteins. However, such re-arrangement of secondary structure connectivity may interfere with the folding mechanism causing accumulation of folding intermediates, which in turn can lead to misfolding. We solved the crystal structure and investigated the folding pathway of a circularly permuted variant of a PDZ domain, SAP97 PDZ2. Our data illustrate how well circular permutation may work as a mechanism for molecular evolution. The circular permutant retains the overall structure and function of the native protein domain. Further, unlike most examples in the literature, this circular permutant displays a folding mechanism that is virtually identical to that of the wild type. This observation contrasts with previous data on the circularly permuted PDZ2 domain from PTP-BL, for which the folding pathway was remarkably affected by the same mutation in sequence connectivity. The different effects of this circular permutation in two homologous proteins show the strong influence of sequence as compared to topology. Circular permutation, when peripheral to the major folding nucleus, may have little effect on folding pathways and could explain why, despite the dramatic change in primary structure, it is frequently tolerated by different protein folds.

Interactions outside the boundaries of the canonical binding groove of a PDZ domain influence ligand binding

The postsynaptic density protein-95/discs large/zonula occludens-1 (PDZ) domain is a protein-protein interaction module with a shallow binding groove where protein ligands bind. However, interactions that are not part of this canonical binding groove are likely to modulate peptide binding. We have investigated such interactions beyond the binding groove for PDZ3 from PSD-95 and a peptide derived from the C-terminus of the natural ligand CRIPT. We found via nuclear magnetic resonance experiments that up to eight residues of the peptide ligand interact with the PDZ domain, showing that the interaction surface extends far outside of the binding groove as defined by the crystal structure. PDZ3 contains an extra structural element, a C-terminal helix (α3), which is known to affect affinity. Deletion of this helix resulted in the loss of several intermolecular nuclear Overhauser enhancements from peptide residues outside of the binding pocket, suggesting that α3 forms part of the extra binding surface in wild-type PDZ3. Site-directed mutagenesis, isothermal titration calorimetry, and fluorescence intensity experiments confirmed the importance of both α3 and the N-terminal part of the peptide for the affinity. Our data suggest a general mechanism in which different binding surfaces outside of the PDZ binding groove could provide sites for specific interactions.

Ligand binding by PDZ domains

The postsynaptic density protein-95/disks large/zonula occludens-1 (PDZ) protein domain family is one of the most common protein-protein interaction modules in mammalian cells, with paralogs present in several hundred human proteins. PDZ domains are found in most cell types, but neuronal proteins, for example, are particularly rich in these domains. The general function of PDZ domains is to bring proteins together within the appropriate cellular compartment, thereby facilitating scaffolding, signaling, and trafficking events. The many functions of PDZ domains under normal physiological as well as pathological conditions have been reviewed recently. In this review, we focus on the molecular details of how PDZ domains bind their protein ligands and their potential as drug targets in this context.

Fast association and slow transitions in the interaction between two intrinsically disordered protein domains

Proteins that contain long disordered regions are prevalent in the proteome and frequently associated with diseases. However, the mechanisms by which such intrinsically disordered proteins (IDPs) recognize their targets are not well understood. Here, we report the first experimental investigation of the interaction kinetics of the nuclear co-activator binding domain of CREB-binding protein and the activation domain from the p160 transcriptional co-activator for thyroid hormone and retinoid receptors. Both protein domains are intrinsically disordered in the free state and synergistically fold upon binding each other. Using the stopped-flow technique, we found that the binding reaction is fast, with an association rate constant of 3 × 10(7) m(-1) s(-1) at 277 K. Mutation of a conserved buried intermolecular salt bridge showed that electrostatics govern the rapid association. Furthermore, upon mutation of the salt bridge or at high salt concentration, an additional kinetic phase was detected (∼20 and ∼40 s(-1), respectively, at 277 K), suggesting that the salt bridge may steer formation of the productive bimolecular complex in an intramolecular step. Finally, we directly measured slow kinetics for the IDP domains (∼1 s(-1) at 277 K) related to conformational transitions upon binding. Together, the experiments demonstrate that the interaction involves several steps and accumulation of intermediate states. Our data are consistent with an induced fit mechanism, in agreement with previous simulations. We propose that the slow transitions may be a consequence of the multipartner interactions of IDPs.

Characterization of the endopeptidase activity of tripeptidyl-peptidase II

Tripeptidyl-peptidase II (TPP II) is a giant cytosolic peptidase with a proposed role in cellular protein degradation and protection against apoptosis. Beside its well-characterised exopeptidase activity, TPP II also has an endopeptidase activity. Little is known about this activity, and since it could be important for the physiological role of TPP II, we have investigated it in more detail. Two peptides, Nef(69-87) and LL37, were incubated with wild-type murine TPP II and variants thereof as well as TPP II from human and Drosophila melanogaster. Two intrinsically disordered proteins were also included in the study. We conclude that the endopeptidase activity is more promiscuous than previously reported. It is also clear that TPP II can attack longer disordered peptides up to 75 amino acid residues. Using a novel FRET substrate, the catalytic efficiency of the endopeptidase activity could be determined to be 5 orders of magnitude lower than for the exopeptidase activity.