Kinetic Recognition of the Retinoblastoma Tumor Suppressor by a Specific Protein Target

Encounter Complexes Between the N-terminal of Neurotensin with the Extracellular Loop 2 of the Neurotensin Receptor 1 Steer Neurotensin to the Orthosteric Binding Pocket

Neurotensin (NT) is a linear disordered peptide that activates two different class A GPCRs, neurotensin receptor 1 (NTS1) and NTS2. Resolved structures of the complex of the C-terminal fragment of NT, NT8-13, with NTS1 shows the peptide takes a well-defined structure in the bound state. However, the mechanisms underlying NT recognition of NTS1, and the conformational transition of NT upon binding NTS1 is an open question that if answered may aid discovery of highly selective drugs and reveal potential secondary binding sites on the surface of the receptor. Herein we investigated the interactions guiding NT to the orthosteric binding pocket of NTS1 by combining NMR experiments with kinetic analysis of the binding pathway using stopped-flow fluorescence and mutagenesis on both NT and NTS1. We show the presence of transient structures in the middle part of NT that kinetically regulate the binding of NT to NTS1. Moreover, our results indicate that the binding pathway of NT onto NTS1 is mediated via electrostatic interactions between the N-terminal region of NT with the extracellular loop 2 of NTS1. These interactions induce backbone conformational changes in neurotensin similar to the bound-state neurotensin, suggesting that the N-terminal region of NT and these interactions should be considered for development of selective drugs against NTS1.

Vital for Viruses: Intrinsically Disordered Proteins

Viruses infect all kingdoms of life; their genomes vary from DNA to RNA and in size from 2kB to 1 MB or more. Viruses frequently employ disordered proteins, that is, protein products of virus genes that do not themselves fold into independent three-dimensional structures, but rather, constitute a versatile molecular toolkit to accomplish a range of functions necessary for viral infection, assembly, and proliferation. Interestingly, disordered proteins have been discovered in almost all viruses so far studied, whether the viral genome consists of DNA or RNA, and whatever the configuration of the viral capsid or other outer covering. In this review, I present a wide-ranging set of stories illustrating the range of functions of IDPs in viruses. The field is rapidly expanding, and I have not tried to include everything. What is included is meant to be a survey of the variety of tasks that viruses accomplish using disordered proteins.

Evolution of SLiM-mediated hijack functions in intrinsically disordered viral proteins

Viruses and their hosts are involved in an ‘arms race’ where they continually evolve mechanisms to overcome each other. It has long been proposed that intrinsic disorder provides a substrate for the evolution of viral hijack functions and that short linear motifs (SLiMs) are important players in this process. Here, we review evidence in support of this tenet from two model systems: the papillomavirus E7 protein and the adenovirus E1A protein. Phylogenetic reconstructions reveal that SLiMs appear and disappear multiple times across evolution, providing evidence of convergent evolution within individual viral phylogenies. Multiple functionally related SLiMs show strong coevolution signals that persist across long distances in the primary sequence and occur in unrelated viral proteins. Moreover, changes in SLiMs are associated with changes in phenotypic traits such as host range and tropism. Tracking viral evolutionary events reveals that host switch events are associated with the loss of several SLiMs, suggesting that SLiMs are under functional selection and that changes in SLiMs support viral adaptation. Fine-tuning of viral SLiM sequences can improve affinity, allowing them to outcompete host counterparts. However, viral SLiMs are not always competitive by themselves, and tethering of two suboptimal SLiMs by a disordered linker may instead enable viral hijack. Coevolution between the SLiMs and the linker indicates that the evolution of disordered regions may be more constrained than previously thought. In summary, experimental and computational studies support a role for SLiMs and intrinsic disorder in viral hijack functions and in viral adaptive evolution.

Structural basis for tunable affinity and specificity of LxCxE-dependent protein interactions with the retinoblastoma protein family

The retinoblastoma protein (Rb) and its homologs p107 and p130 are critical regulators of gene expression during the cell cycle and are commonly inactivated in cancer. Rb proteins use their "pocket domain" to bind an LxCxE sequence motif in other proteins, many of which function with Rb proteins to co-regulate transcription. Here, we present binding data and crystal structures of the p107 pocket domain in complex with LxCxE peptides from the transcriptional co-repressor proteins HDAC1, ARID4A, and EID1. Our results explain why Rb and p107 have weaker affinity for cellular LxCxE proteins compared with the E7 protein from human papillomavirus, which has been used as the primary model for understanding LxCxE motif interactions. Our structural and mutagenesis data also identify and explain differences in Rb and p107 affinities for some LxCxE-containing sequences. Our study provides new insights into how Rb proteins bind their cell partners with varying affinity and specificity.

Comparative Analysis of Structural Features in SLiMs from Eukaryotes, Bacteria, and Viruses with Importance for Host-Pathogen Interactions

Protein-protein interactions drive functions in eukaryotes that can be described by short linear motifs (SLiMs). Conservation of SLiMs help illuminate functional SLiMs in eukaryotic protein families. However, the simplicity of eukaryotic SLiMs makes them appear by chance due to mutational processes not only in eukaryotes but also in pathogenic bacteria and viruses. Further, functional eukaryotic SLiMs are often found in disordered regions. Although proteomes from pathogenic bacteria and viruses have less disorder than eukaryotic proteomes, their proteins can successfully mimic eukaryotic SLiMs and disrupt host cellular function. Identifying important SLiMs in pathogens is difficult but essential for understanding potential host-pathogen interactions. We performed a comparative analysis of structural features for experimentally verified SLiMs from the Eukaryotic Linear Motif (ELM) database across viruses, bacteria, and eukaryotes. Our results revealed that many viral SLiMs and specific motifs found across viruses and eukaryotes, such as some glycosylation motifs, have less disorder. Analyzing the disorder and coil properties of equivalent SLiMs from pathogens and eukaryotes revealed that some motifs are more structured in pathogens than their eukaryotic counterparts and vice versa. These results support a varying mechanism of interaction between pathogens and their eukaryotic hosts for some of the same motifs.

The Intrinsically Disordered Proteins MLLT3 (AF9) and MLLT1 (ENL) - Multimodal Transcriptional Switches With Roles in Normal Hematopoiesis, MLL Fusion Leukemia, and Kidney Cancer

AF9 (MLLT3) and ENL (MLLT1) are members of the YEATS family (named after the five proteins first shown to contain this domain: Yaf9, ENL, AF9, Taf14, Sas5) defined by the presence of a YEATS domain. The YEATS domain is an epigenetic reader that binds to acetylated and crotonylated lysines, unlike the bromodomain which can only bind to acetylated lysines. All members of this family have been shown to be components of various complexes with roles in chromatin remodeling, histone modification, histone variant deposition, and transcriptional regulation. MLLT3 is a critical regulator of hematopoiesis with a role in maintaining the hematopoietic stem or progenitor cell (HSPC) population. Approximately 10% of acute myeloid leukemia (AML) and acute lymphocytic leukemia (ALL) patients harbor a translocation involving MLL (mixed lineage leukemia). In the context of MLL fusion patients with AML and ALL, MLL-AF9 and MLL-ENL fusions are observed in 34 and 31% of the patients, respectively. The intrinsically disordered C-terminal domain of MLLT3 (AHD, ANC1 homology domain) undergoes coupled binding and folding upon interaction with partner proteins AF4, DOT1L, BCOR, and CBX8. Backbone dynamics studies of the complexes suggest a role for dynamics in function. Inhibitors of the interaction of the intrinsically disordered AHD with partner proteins have been described, highlighting the feasibility of targeting intrinsically disordered regions. MLLT1 undergoes phase separation to enhance recruitment of the super elongation complex (SEC) and drive transcription. Mutations in MLLT1 observed in Wilms tumor patients enhance phase separation and transcription to drive an aberrant gene expression program.

Unveiling induced folding of intrinsically disordered proteins - Protein engineering, frustration and emerging themes

Intrinsically disordered proteins (IDPs) can be generally described as a class of proteins that lack a well-defined ordered structure in isolation at physiological conditions. Upon binding to their physiological ligands, IDPs typically undergo a disorder-to-order transition, which may or may not lead to the complete folding of the IDP. In this short review, we focus on some of the key findings pertaining to the mechanisms of such induced folding. In particular, first we describe the general features of the reaction; then, we discuss some of the most remarkable findings obtained from applying protein engineering in synergy with kinetic studies to induced folding; and finally, we offer a critical view on some of the emerging themes when considering the structural heterogeneity of IDPs vis-à-vis to their inherent frustration.

Parallel and Sequential Pathways of Molecular Recognition of a Tandem-Repeat Protein and Its Intrinsically Disordered Binding Partner

The Wnt signalling pathway plays an important role in cell proliferation, differentiation, and fate decisions in embryonic development and the maintenance of adult tissues. The twelve armadillo (ARM) repeat-containing protein β-catenin acts as the signal transducer in this pathway. Here, we investigated the interaction between β-catenin and the intrinsically disordered transcription factor TCF7L2, comprising a very long nanomolar-affinity interface of approximately 4800 Å2 that spans ten of the twelve ARM repeats of β-catenin. First, a fluorescence reporter system for the interaction was engineered and used to determine the kinetic rate constants for the association and dissociation. The association kinetics of TCF7L2 and β-catenin were monophasic and rapid (7.3 ± 0.1 × 107 M-1·s-1), whereas dissociation was biphasic and slow (5.7 ± 0.4 × 10-4 s-1, 15.2 ± 2.8 × 10-4 s-1). This reporter system was then combined with site-directed mutagenesis to investigate the striking variability in the conformation adopted by TCF7L2 in the three different crystal structures of the TCF7L2-β-catenin complex. We found that the mutation had very little effect on the association kinetics, indicating that most interactions form after the rate-limiting barrier for association. Mutations of the N- and C-terminal subdomains of TCF7L2 that adopt relatively fixed conformations in the crystal structures had large effects on the dissociation kinetics, whereas the mutation of the labile sub-domain connecting them had negligible effect. These results point to a two-site avidity mechanism of binding with the linker region forming a "fuzzy" complex involving transient contacts that are not site-specific. Strikingly, the two mutations in the N-terminal subdomain that had the largest effects on the dissociation kinetics showed two additional phases, indicating partial flux through an alternative dissociation pathway that is inaccessible to the wild type. The results presented here provide insights into the kinetics of the molecular recognition of a long intrinsically disordered region with an elongated repeat-protein surface, a process found to involve parallel routes with sequential steps in each.

Kinetic Methods of Deducing Binding Mechanisms Involving Intrinsically Disordered Proteins

There are multiple examples of protein-protein interactions involving one intrinsically disordered protein region binding to an ordered protein domain in a coupled binding and folding reaction. Similarly to protein folding studies, much effort has been devoted to understanding the mechanisms of such coupled binding and folding reactions. In this chapter, we describe how kinetics can be used to assess binding mechanisms with focus on fluorescence-monitored stopped-flow experiments. The approach can be applied more generally to any protein interaction with or without a coupled conformational change and to other kinetic techniques. Determining binding mechanisms is a great challenge and while "proving" a mechanism may be futile, it is possible to deduce the simplest scenarios, which are consistent with experimental data.

A disordered encounter complex is central to the yeast Abp1p SH3 domain binding pathway

Protein-protein interactions are involved in a wide range of cellular processes. These interactions often involve intrinsically disordered proteins (IDPs) and protein binding domains. However, the details of IDP binding pathways are hard to characterize using experimental approaches, which can rarely capture intermediate states present at low populations. SH3 domains are common protein interaction domains that typically bind proline-rich disordered segments and are involved in cell signaling, regulation, and assembly. We hypothesized, given the flexibility of SH3 binding peptides, that their binding pathways include multiple steps important for function. Molecular dynamics simulations were used to characterize the steps of binding between the yeast Abp1p SH3 domain (AbpSH3) and a proline-rich IDP, ArkA. Before binding, the N-terminal segment 1 of ArkA is pre-structured and adopts a polyproline II helix, while segment 2 of ArkA (C-terminal) adopts a 310 helix, but is far less structured than segment 1. As segment 2 interacts with AbpSH3, it becomes more structured, but retains flexibility even in the fully engaged state. Binding simulations reveal that ArkA enters a flexible encounter complex before forming the fully engaged bound complex. In the encounter complex, transient nonspecific hydrophobic and long-range electrostatic contacts form between ArkA and the binding surface of SH3. The encounter complex ensemble includes conformations with segment 1 in both the forward and reverse orientation, suggesting that segment 2 may play a role in stabilizing the correct binding orientation. While the encounter complex forms quickly, the slow step of binding is the transition from the disordered encounter ensemble to the fully engaged state. In this transition, ArkA makes specific contacts with AbpSH3 and buries more hydrophobic surface. Simulating the binding between ApbSH3 and ArkA provides insight into the role of encounter complex intermediates and nonnative hydrophobic interactions for other SH3 domains and IDPs in general.

Energetics and kinetics of substrate analog-coupled staphylococcal nuclease folding revealed by a statistical mechanical approach

Protein conformational changes associated with ligand binding, especially those involving intrinsically disordered proteins, are mediated by tightly coupled intra- and intermolecular events. Such reactions are often discussed in terms of two limiting kinetic mechanisms, conformational selection (CS), where folding precedes binding, and induced fit (IF), where binding precedes folding. It has been shown that coupled folding/binding reactions can proceed along both CS and IF pathways with the flux ratio depending on conditions such as ligand concentration. However, the structural and energetic basis of such complex reactions remains poorly understood. Therefore, we used experimental, theoretical, and computational approaches to explore structural and energetic aspects of the coupled-folding/binding reaction of staphylococcal nuclease in the presence of the substrate analog adenosine-3',5'-diphosphate. Optically monitored equilibrium and kinetic data, combined with a statistical mechanical model, gave deeper insight into the relative importance of specific and Coulombic protein-ligand interactions in governing the reaction mechanism. We also investigated structural aspects of the reaction at the residue level using NMR and all-atom replica-permutation molecular dynamics simulations. Both approaches yielded clear evidence for accumulation of a transient protein-ligand encounter complex early in the reaction under IF-dominant conditions. Quantitative analysis of the equilibrium/kinetic folding revealed that the ligand-dependent CS-to-IF shift resulted from stabilization of the compact transition state primarily by weakly ligand-dependent Coulombic interactions with smaller contributions from specific binding energies. At a more macroscopic level, the CS-to-IF shift was represented as a displacement of the reaction "route" on the free energy surface, which was consistent with a flux analysis.

Templated folding of intrinsically disordered proteins

Much of our current knowledge of biological chemistry is founded in the structure-function relationship, whereby sequence determines structure that determines function. Thus, the discovery that a large fraction of the proteome is intrinsically disordered, while being functional, has revolutionized our understanding of proteins and raised new and interesting questions. Many intrinsically disordered proteins (IDPs) have been determined to undergo a disorder-to-order transition when recognizing their physiological partners, suggesting that their mechanisms of folding are intrinsically different from those observed in globular proteins. However, IDPs also follow some of the classic paradigms established for globular proteins, pointing to important similarities in their behavior. In this review, we compare and contrast the folding mechanisms of globular proteins with the emerging features of binding-induced folding of intrinsically disordered proteins. Specifically, whereas disorder-to-order transitions of intrinsically disordered proteins appear to follow rules of globular protein folding, such as the cooperative nature of the reaction, their folding pathways are remarkably more malleable, due to the heterogeneous nature of their folding nuclei, as probed by analysis of linear free-energy relationship plots. These insights have led to a new model for the disorder-to-order transition in IDPs termed "templated folding," whereby the binding partner dictates distinct structural transitions en route to product, while ensuring a cooperative folding.

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.

A structurally heterogeneous transition state underlies coupled binding and folding of disordered proteins

Many intrinsically disordered proteins (IDPs) attain a well-defined structure in a coupled folding and binding reaction with another protein. Such reactions may involve early to late formation of different native structural regions along the reaction pathway. To obtain insights into the transition state for a coupled binding and folding reaction, we performed restrained molecular dynamics simulations using previously determined experimental binding Φ_b values of the interaction between two IDP domains: the activation domain from the p160 transcriptional co-activator for thyroid hormone and retinoid receptors (ACTR) and the nuclear co-activator binding domain (NCBD) of CREB-binding protein, each forming three well-defined α-helices upon binding. These simulations revealed that both proteins are largely disordered in the transition state for complex formation, except for two helices, one from each domain, that display a native-like structure. The overall transition state structure was extended and largely dynamic with many weakly populated contacts. To test the transition state model, we combined site-directed mutagenesis with kinetic experiments, yielding results consistent with overall diffuse interactions and formation of native intramolecular interactions in the third NCBD helix during the binding reaction. Our findings support the view that the transition state and, by inference, any encounter complex in coupled binding and folding reactions are structurally heterogeneous and largely independent of specific interactions. Furthermore, experimental Φ_b values and Brønsted plots suggested that the transition state is globally robust with respect to most mutations but can display more native-like features for some highly destabilizing mutations, possibly because of Hammond behavior or ground-state effects.

Interplay between sequence, structure and linear motifs in the adenovirus E1A hub protein

E1A is the main transforming protein in mastadenoviruses. This work uses bioinformatics to extrapolate experimental knowledge from Human adenovirus serotype 5 and 12 E1A proteins to all known serotypes. A conserved domain architecture with a high degree of intrinsic disorder acts as a scaffold for multiple linear motifs with variable occurrence mediating the interaction with over fifty host proteins. While linear motifs contribute strongly to sequence conservation within intrinsically disordered E1A regions, motif repertoires can deviate significantly from those found in prototypical serotypes. Close to one hundred predicted residue-residue contacts suggest the presence of stable structure in the CR3 domain and of specific conformational ensembles involving both short- and long-range intramolecular interactions. Our computational results suggest that E1A sequence conservation and co-evolution reflect the evolutionary pressure to maintain a mainly disordered, yet non-random conformation harboring a high number of binding motifs that mediate viral hijacking of the cell machinery.

A Mathematical Model of Cell Cycle Dysregulation Due to Human Papillomavirus Infection

Human papillomaviruses (HPVs) that infect mucosal epithelium can be classified as high risk or low risk based on their propensity to cause lesions that can undergo malignant progression. HPVs produce the E7 protein that binds to cell cycle regulatory proteins including the retinoblastoma tumor suppressor protein (RB) to modulate cell cycle control. Generally, high-risk HPV E7 proteins bind to RB with a higher affinity than low-risk HPV E7s, but both are able to deactivate RB and trigger S phase progression. In uninfected cells, RB inactivation is a tightly controlled process that must coincide with growth factor stimulation to commit cells to division. High-risk HPV E7 proteins short-circuit this control by decreasing growth factor requirement for cell division. We develop a mathematical model to examine the role that RB binding affinity, growth factor concentration, and E7 concentration have on cell cycle progression. Our model predicts that high RB binding affinity and E7 concentration accelerate the [Formula: see text] to S phase transition and weaken the dependence on growth factor. This model thus captures a key step in high-risk HPV oncogenesis.

APC/C and retinoblastoma interaction: cross-talk of retinoblastoma protein with the ubiquitin proteasome pathway

The ubiquitin (Ub) ligase anaphase promoting complex/cyclosome (APC/C) and the tumour suppressor retinoblastoma protein (pRB) play key roles in cell cycle regulation. APC/C is a critical regulator of mitosis and G₁-phase of the cell cycle whereas pRB keeps a check on proliferation by inhibiting transition to the S-phase. APC/C and pRB interact with each other via the co-activator of APC/C, FZR1, providing an alternative pathway of regulation of G₁ to S transition by pRB using a post-translational mechanism. Both pRB and FZR1 have complex roles and are implicated not only in regulation of cell proliferation but also in differentiation, quiescence, apoptosis, maintenance of chromosomal integrity and metabolism. Both are also targeted by transforming viruses. We discuss recent advances in our understanding of the involvement of APC/C and pRB in cell cycle based decisions and how these insights will be useful for development of anti-cancer and anti-viral drugs.

Fuzzy regions in an intrinsically disordered protein impair protein-protein interactions

Despite the partial disorder-to-order transition that intrinsically disordered proteins often undergo upon binding to their partners, a considerable amount of residual disorder may be retained in the bound form, resulting in a fuzzy complex. Fuzzy regions flanking molecular recognition elements may enable partner fishing through non-specific, transient contacts, thereby facilitating binding, but may also disfavor binding through various mechanisms. So far, few computational or experimental studies have addressed the effect of fuzzy appendages on partner recognition by intrinsically disordered proteins. In order to shed light onto this issue, we used the interaction between the intrinsically disordered C-terminal domain of the measles virus (MeV) nucleoprotein (NTAIL ) and the X domain (XD) of the viral phosphoprotein as model system. After binding to XD, the N-terminal region of NTAIL remains conspicuously disordered, with α-helical folding taking place only within a short molecular recognition element. To study the effect of the N-terminal fuzzy region on NTAIL /XD binding, we generated N-terminal truncation variants of NTAIL , and assessed their binding abilities towards XD. The results revealed that binding increases with shortening of the N-terminal fuzzy region, with this also being observed with hsp70 (another MeV NTAIL binding partner), and for the homologous NTAIL /XD pairs from the Nipah and Hendra viruses. Finally, similar results were obtained when the MeV NTAIL fuzzy region was replaced with a highly dissimilar artificial disordered sequence, supporting a sequence-independent inhibitory effect of the fuzzy region.

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.

Cysteine-rich positions outside the structural zinc motif of human papillomavirus E7 provide conformational modulation and suggest functional redox roles

The E7 protein from high-risk human papillomavirus is essential for cell transformation in cervical, oropharyngeal, and other HPV-related cancers, mainly through the inactivation of the retinoblastoma (Rb) tumor suppressor. Its high cysteine content (~7%) and the observation that HPV-transformed cells are under oxidative stress prompted us to investigate the redox properties of the HPV16 E7 protein under biologically compatible oxidative conditions. The seven cysteines in HPV16 E7 remain reduced in conditions resembling the basal reduced state of a cell. However, under oxidative stress, a stable disulfide bridge forms between cysteines 59 and 68. Residue 59 has a protective effect on the other cysteines, and its mutation leads to an overall increase in the oxidation propensity of E7, including cysteine 24 central to the Rb binding motif. Gluthationylation of Cys 24 abolishes Rb binding, which is reversibly recovered upon reduction. Cysteines 59 and 68 are located 18.6 Å apart, and the formation of the disulfide bridge leads to a large structural rearrangement while retaining strong Zn association. These conformational and covalent changes are fully reversible upon restoration of the reductive environment. In addition, this is the first evidence of an interaction between the N-terminal intrinsically disordered and the C-terminal globular domains, known to be highly and separately conserved among human papillomaviruses. The significant conservation of such noncanonical cysteines in HPV E7 proteins leads us to propose a functional redox activity. Such an activity adds to the previously discovered chaperone activity of E7 and supports the picture of a moonlighting pathological role of this paradigmatic viral oncoprotein.

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).

Conformational dissection of a viral intrinsically disordered domain involved in cellular transformation

Intrinsic disorder is abundant in viral genomes and provides conformational plasticity to its protein products. In order to gain insight into its structure-function relationships, we carried out a comprehensive analysis of structural propensities within the intrinsically disordered N-terminal domain from the human papillomavirus type-16 E7 oncoprotein (E7N). Two E7N segments located within the conserved CR1 and CR2 regions present transient α-helix structure. The helix in the CR1 region spans residues L8 to L13 and overlaps with the E2F mimic linear motif. The second helix, located within the highly acidic CR2 region, presents a pH-dependent structural transition. At neutral pH the helix spans residues P17 to N29, which include the retinoblastoma tumor suppressor LxCxE binding motif (residues 21-29), while the acidic CKII-PEST region spanning residues E33 to I38 populates polyproline type II (PII) structure. At pH 5.0, the CR2 helix propagates up to residue I38 at the expense of loss of PII due to charge neutralization of acidic residues. Using truncated forms of HPV-16 E7, we confirmed that pH-induced changes in α-helix content are governed by the intrinsically disordered E7N domain. Interestingly, while at both pH the region encompassing the LxCxE motif adopts α-helical structure, the isolated 21-29 fragment including this stretch is unable to populate an α-helix even at high TFE concentrations. Thus, the E7N domain can populate dynamic but discrete structural ensembles by sampling α-helix-coil-PII-ß-sheet structures. This high plasticity may modulate the exposure of linear binding motifs responsible for its multi-target binding properties, leading to interference with key cell signaling pathways and eventually to cellular transformation by the virus.

Folding of a cyclin box: linking multitarget binding to marginal stability, oligomerization, and aggregation of the retinoblastoma tumor suppressor AB pocket domain

The retinoblastoma tumor suppressor (Rb) controls the proliferation, differentiation, and survival of cells in most eukaryotes with a role in the fate of stem cells. Its inactivation by mutation or oncogenic viruses is required for cellular transformation and eventually carcinogenesis. The high conservation of the Rb cyclin fold prompted us to investigate the link between conformational stability and ligand binding properties of the RbAB pocket domain. RbAB unfolding presents a three-state transition involving cooperative secondary and tertiary structure changes and a partially folded intermediate that can oligomerize. The first transition corresponds to unfolding of the metastable B subdomain containing the binding site for the LXCXE motif present in cellular and viral targets, and the second transition corresponds to the stable A subdomain. The low thermodynamic stability of RbAB translates into a propensity to rapidly oligomerize and aggregate at 37 °C (T50 = 28 min) that is suppressed by human papillomavirus E7 and E2F peptide ligands, suggesting that Rb is likely stabilized in vivo through binding to target proteins. We propose that marginal stability and associated oligomerization may be conserved for function as a "hub" protein, allowing the formation of multiprotein complexes, which could constitute a robust mechanism to retain its cell cycle regulatory role throughout evolution. Decreased stability and oligomerization are shared with the p53 tumor suppressor, suggesting a link between folding and function in these two essential cell regulators that are inactivated in most cancers and operate within multitarget signaling pathways.

Minute time scale prolyl isomerization governs antibody recognition of an intrinsically disordered immunodominant epitope

Conformational rearrangements in antibody·antigen recognition are essential events where kinetic discrimination of isomers expands the universe of combinations. We investigated the interaction mechanism of a monoclonal antibody, M1, raised against E7 from human papillomavirus, a prototypic viral oncoprotein and a model intrinsically disordered protein. The mapped 12-amino acid immunodominant epitope lies within a "hinge" region between the N-terminal intrinsically disordered and the C-terminal globular domains. Kinetic experiments show that despite being within an intrinsically disordered region, the hinge E7 epitope has at least two populations separated by a high energy barrier. Nuclear magnetic resonance traced the origin of this barrier to a very slow (t(1/2)∼4 min) trans-cis prolyl isomerization event involving changes in secondary structure. The less populated (10%) cis isomer is the binding-competent species, thus requiring the 90% of molecules in the trans configuration to isomerize before binding. The association rate for the cis isomer approaches 6 × 10(7) M(-1) s(-1), a ceiling for antigen-antibody interactions. Mutagenesis experiments showed that Pro-41 in E7Ep was required for both binding and isomerization. After a slow postbinding unimolecular rearrangement, a consolidated complex with K(D) = 1.2 × 10(-7) M is reached. Our results suggest that presentation of this viral epitope by the antigen-presenting cells would have to be "locked" in the cis conformation, in opposition to the most populated trans isomer, in order to select the specific antibody clone that goes through affinity and kinetic maturation.

Slow, reversible, coupled folding and binding of the spectrin tetramerization domain

Many intrinsically disordered proteins (IDPs) are significantly unstructured under physiological conditions. A number of these IDPs have been shown to undergo coupled folding and binding reactions whereby they can gain structure upon association with an appropriate partner protein. In general, these systems display weaker binding affinities than do systems with association between completely structured domains, with micromolar K(d) values appearing typical. One such system is the association between α- and β-spectrin, where two partially structured, incomplete domains associate to form a fully structured, three-helix bundle, the spectrin tetramerization domain. Here, we use this model system to demonstrate a method for fitting association and dissociation kinetic traces where, using typical biophysical concentrations, the association reactions are expected to be highly reversible. We elucidate the unusually slow, two-state kinetics of spectrin assembly in solution. The advantages of studying kinetics in this regime include the potential for gaining equilibrium constants as well as rate constants, and for performing experiments with low protein concentrations. We suggest that this approach would be particularly appropriate for high-throughput mutational analysis of two-state reversible binding processes.

Sequence evolution of the intrinsically disordered and globular domains of a model viral oncoprotein

In the present work, we have used the papillomavirus E7 oncoprotein to pursue structure-function and evolutionary studies that take into account intrinsic disorder and the conformational diversity of globular domains. The intrinsically disordered (E7N) and globular (E7C) domains of E7 show similar degrees of conservation and co-evolution. We found that E7N can be described in terms of conserved and coevolving linear motifs separated by variable linkers, while sequence evolution of E7C is compatible with the known homodimeric structure yet suggests other activities for the domain. Within E7N, inter-residue relationships such as residue co-evolution and restricted intermotif distances map functional coupling and co-occurrence of linear motifs that evolve in a coordinate manner. Within E7C, additional cysteine residues proximal to the zinc-binding site may allow redox regulation of E7 function. Moreover, we describe a conserved binding site for disordered domains on the surface of E7C and suggest a putative target linear motif. Both homodimerization and peptide binding activities of E7C are also present in the distantly related host PHD domains, showing that these two proteins share not only structural homology but also functional similarities, and strengthening the view that they evolved from a common ancestor. Finally, we integrate the multiple activities and conformations of E7 into a hierarchy of structure-function relationships.

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.

Why is F19Ap53 unable to bind MDM2? Simulations suggest crack propagation modulates binding

Why doesn't the F19A mutant of p53 bind to MDM2? Binding thermodynamics have suggested that the loss of packing interactions upon mutating Phe into Ala sidechain results in destabilizing the binding free energy between p53 and MDM2. Does this mutation also modulate the initial recognition between p53 and MDM2? We look at atomistic computer simulations of the process of the initial encounter between wild type p53 peptide and its F19A mutant with the N-terminal domain of MDM2. These simulations show that binding is characterized by a complex multistep process. It starts with the capture of F19 of wild type p53 by certain residues in the MDM2 binding pocket. This initial step anchors the peptide onto the surface of MDM2, and with the consequent reduction in the search space of the peptide, the peptide docks into the partially occluded surface of MDM2. This is similar to a crack forming in an otherwise occluded hydrophobic cavity in MDM2, and the peptide, docked through F19, modulates the propagation of this crack, which subsequently results in the stepwise docking of the rest of the peptide through insertions of W23 and L26. The lack of the bulky sidechain of F in the F19A mutant results in the absence of the initial "grasp" complex, and hence the mutant peptide diffuses randomly on the surface of MDM2 without binding. This is the first such demonstration of the possibility that a "kinetic" effect may partly underlie the destabilized thermodynamics of binding of F19A and is a feature that appears to be conserved in evolution. The observations by Wallace et al. (Mol Cell 2006; 23:251-63) that despite the inability of F19A to bind at the N-terminal domain of MDM2, it gets ubiquitinated, can now be partly understood based on a mechanism whereby the occupation of the binding pocket by ligands/peptides induces, via crack propagation and the dynamics of gatekeeper Y100, the ubiquitination signal for interactions between the acidic domain of MDM2 and the DNA binding domain of p53.

Evolution of linear motifs within the papillomavirus E7 oncoprotein

Many protein functions can be traced to linear sequence motifs of less than five residues, which are often found within intrinsically disordered domains. In spite of their prevalence, their role in protein evolution is only beginning to be understood. The study of papillomaviruses has provided many insights on the evolution of protein structure and function. We have chosen the papillomavirus E7 oncoprotein as a model system for the evolution of functional linear motifs. The multiple functions of E7 proteins from paradigmatic papillomavirus types can be explained to a large extent in terms of five linear motifs within the intrinsically disordered N-terminal domain and two linear motifs within the globular homodimeric C-terminal domain. We examined the motif inventory of E7 proteins from over 200 known papillomavirus types and found that the motifs reported for paradigmatic papillomavirus types are absent from many uncharacterized E7 proteins. Several motif pairs occur more often than expected, suggesting that linear motifs may evolve and function in a cooperative manner. The E7 linear motifs have appeared or disappeared multiple times during papillomavirus evolution, confirming the evolutionary plasticity of short functional sequences. Four of the motifs appeared several times during papillomavirus evolution, providing direct evidence for convergent evolution. Interestingly, the evolution pattern of a motif is independent of its location in a globular or disordered domain. The correlation between the presence of some motifs and virus host specificity and tissue tropism suggests that linear motifs play a role in the adaptive evolution of papillomaviruses.

The transition state of coupled folding and binding for a flexible β-finger

Flexible and fully disordered protein regions that fold upon binding mediate numerous protein-protein interactions. However, little is known about their mechanism of interaction. One such coupled folding and binding occurs when a flexible region of neuronal nitric oxide synthase adopts a β-finger structure upon binding to its protein ligand, a PDZ [PSD-95 (postsynaptic density protein-95)/Discs large/ZO-1] domain from PSD-95. We have analyzed this binding reaction by protein engineering combined with kinetic experiments. Mutational destabilization of the β-finger changed mainly the dissociation rate constant of the proteins and, to a lesser extent, the association rate constant. Thus, mutation affected late events in the coupled folding and binding reaction. Our results therefore suggest that the native binding interactions of the β-finger are not present in the rate-limiting transition state for binding but form on the downhill side in a cooperative manner. However, by mutation, we could destabilize the β-finger further and change the rate-limiting step such that an initial conformational change becomes rate limiting. This switch in rate-limiting step shows that multistep binding mechanisms are likely to be found among flexible and intrinsically disordered regions of proteins.