Intrinsically disordered proteins (IDPs) are characterized by a lack of persistent structure. Since their identification more than a decade ago, many questions regarding their functional relevance and interaction mechanisms remain unanswered. Although most experiments have taken equilibrium and structural perspectives, fewer studies have investigated the kinetics of their interactions. Here we review and highlight the type of information that can be gained from kinetic studies. In particular, we show how kinetic studies of coupled folding and binding reactions, an important class of signaling event, are needed to determine mechanisms.

Understanding and control of structures and rates involved in protein ligand binding are essential for drug design. Unfortunately, atomistic molecular dynamics (MD) simulations cannot directly sample the excessively long residence and rearrangement times of tightly binding complexes. Here we exploit the recently developed multi-ensemble Markov model framework to compute full protein-peptide kinetics of the oncoprotein fragment 25–109Mdm2 and the nano-molar inhibitor peptide PMI. Using this system, we report, for the first time, direct estimates of kinetics beyond the seconds timescale using simulations of an all-atom MD model, with high accuracy and precision. These results only require explicit simulations on the sub-milliseconds timescale and are tested against existing mutagenesis data and our own experimental measurements of the dissociation and association rates. The full kinetic model reveals an overall downhill but rugged binding funnel with multiple pathways. The overall strong binding arises from a variety of conformations with different hydrophobic contact surfaces that interconvert on the milliseconds timescale.

Appropriate integration of cellular signals requires a delicate balance of ligand–target binding affinities. Increasing the level of residual structure in intrinsically disordered proteins (IDPs), which are overrepresented in these cellular processes, has been shown previously to enhance binding affinities and alter cellular function. Conserved proline residues are commonly found flanking regions of IDPs that become helical upon interacting with a partner protein. Here, we mutate these helix-flanking prolines in p53 and MLL and find opposite effects on binding affinity upon an increase in free IDP helicity. In both cases, changes in affinity were due to alterations in dissociation, not association, rate constants, which is inconsistent with conformational selection mechanisms. We conclude that, contrary to previous suggestions, helix-flanking prolines do not regulate affinity by modulating the rate of complex formation. Instead, they influence binding affinities by controlling the lifetime of the bound complex.

Intrinsically disordered regions are present in one-third of eukaryotic proteins and are overrepresented in cellular processes such as signaling, suggesting that intrinsically disordered proteins (IDPs) may have a functional advantage over folded proteins. Upon interacting with a partner macromolecule, a subset of IDPs can fold and bind to form a well-defined three-dimensional conformation. For example, disordered BH3-only proteins bind promiscuously to a large number of homologous BCL-2 family proteins, where they fold to a helical structure in a groove on the BCL-2–like protein surface. As two protein chains are involved in the folding reaction, and the structure is only formed in the presence of the partner macromolecule, this raises the question of where the folding information is encoded. Here, we examine these coupled folding and binding reactions to determine which component determines the folding and binding pathway. Using Φ value analysis to compare transition state interactions between the disordered BH3-only proteins PUMA and BID and the folded BCL-2–like proteins A1 and MCL-1, we found that, even though the BH3-only protein is disordered in isolation and requires a stabilizing partner to fold, its folding and binding pathway is encoded in the IDP itself; the reaction is not templated by the folded partner. We suggest that, by encoding both its transition state and level of residual structure, an IDP can evolve a specific kinetic profile, which could be a crucial functional advantage of disorder.

Many intrinsically disordered proteins (IDPs) participate in coupled folding and binding reactions and form alpha helical structures in their bound complexes. Alanine, glycine, or proline scanning mutagenesis approaches are often used to dissect the contributions of intrinsic helicities to coupled folding and binding. These experiments can yield confounding results because the mutagenesis strategy changes the amino acid compositions of IDPs. Therefore, an important next step in mutagenesis-based approaches to mechanistic studies of coupled folding and binding is the design of sequences that satisfy three major constraints. These are (i) achieving a target intrinsic alpha helicity profile; (ii) fixing the positions of residues corresponding to the binding interface; and (iii) maintaining the native amino acid composition. Here, we report the development of a G: enetic A: lgorithm for D: esign of I: ntrinsic secondary S: tructure (GADIS) for designing sequences that satisfy the specified constraints. We describe the algorithm and present results to demonstrate the applicability of GADIS by designing sequence variants of the intrinsically disordered PUMA system that undergoes coupled folding and binding to Mcl-1. Our sequence designs span a range of intrinsic helicity profiles. The predicted variations in sequence-encoded mean helicities are tested against experimental measurements.

The realization of the therapeutic potential of targeting the M1 muscarinic acetylcholine receptor (mAChR) for the treatment of cognitive decline in Alzheimer’s disease has prompted the discovery of M1 mAChR ligands showing efficacy in alleviating cognitive dysfunction in both rodents and humans. Among these is GSK1034702 (7-fluoro-5-methyl-3-[1-(oxan-4-yl)piperidin-4-yl]-1H-benzimidazol-2-one), described previously as a potent M1 receptor allosteric agonist, which showed procognitive effects in rodents and improved immediate memory in a clinical nicotine withdrawal test but induced significant side effects. Here we provide evidence using ligand binding, chemical biology and functional assays to establish that rather than the allosteric mechanism claimed, GSK1034702 interacts in a bitopic manner at the M1 mAChR such that it can concomitantly span both the orthosteric and an allosteric binding site. The bitopic nature of GSK1034702, together with the intrinsic agonist activity and a lack of muscarinic receptor subtype selectivity reported here, all likely contribute to the adverse effects of this molecule in clinical trials. Although they impart beneficial effects on learning and memory, we conclude that these properties are undesirable in a clinical candidate due to the likelihood of adverse side effects. Rather, our data support the notion that “pure” positive allosteric modulators showing selectivity for the M1 mAChR with low levels of intrinsic activity would be preferable to provide clinical efficacy with low adverse responses.

RalA and RalB are members of the Ras family of small G proteins and are activated downstream of Ras via RalGEFs. The RalGEF-Ral axis represents one of the major effector pathways controlled by Ras and as such is an important pharmacological target. RalA and RalB are approximately 80% identical at the amino acid level; despite this, they have distinct roles both in normal cells and in the disease state. We have used our structure of RalB-RLIP76 to guide an analysis of Ral-effector interaction interfaces, creating panels of mutant proteins to probe the energetics of these interactions. The data provide a physical mechanism that underpins the effector selective mutations commonly employed to dissect Ral G protein function. Comparing the energetic landscape of the RalB-RLIP76 and RalB-Sec5 complexes reveals mutations in RalB that lead to differential binding of the two effector proteins. A panel of RLIP76 mutants was used to probe the interaction between RLIP76 and RalA and -B. Despite 100% sequence identity in the RalA and -B contact residues with RLIP76, differences still exist in the energetic profiles of the two complexes. Therefore, we have revealed properties that may account for some of the functional separation observed with RalA and RalB at the cellular level. Our mutations, in both the Ral isoforms and RLIP76, provide new tools that can be employed to parse the complex biology of Ral G protein signaling networks. The combination of these thermodynamic and structural data can also guide efforts to ablate RalA and -B activity with small molecules and peptides.

In the search for improved symptomatic treatment options for neurodegenerative and neuropsychiatric diseases, muscarinic acetylcholine M1 receptors (M1 mAChRs) have received significant attention. Drug development efforts have identified a number of novel ligands, some of which have advanced to the clinic. However, a significant issue for progressing these therapeutics is the lack of robust, translatable, and validated biomarkers. One valuable approach to assessing target engagement is to use positron emission tomography (PET) tracers. In this study we describe the pharmacological characterization of a selective M1 agonist amenable for in vivo tracer studies. We used a novel direct binding assay to identify nonradiolabeled ligands, including LSN3172176, with the favorable characteristics required for a PET tracer. In vitro functional and radioligand binding experiments revealed that LSN3172176 was a potent partial agonist (EC 2.4-7.0 nM, E 43%-73%), displaying binding selectivity for M1 mAChRs (K = 1.5 nM) that was conserved across species (native tissue K = 1.02, 2.66, 8, and 1.03 at mouse, rat, monkey, and human, respectively). Overall selectivity of LSN3172176 appeared to be a product of potency and stabilization of the high-affinity state of the M1 receptor, relative to other mAChR subtypes (M1 > M2, M4, M5 > M3). In vivo, use of wild-type and mAChR knockout mice further supported the M1-preferring selectivity profile of LSN3172176 for the M1 receptor (78% reduction in cortical occupancy in M1 KO mice). These findings support the development of LSN3172176 as a potential PET tracer for assessment of M1 mAChR target engagement in the clinic and to further elucidate the function of M1 mAChRs in health and disease.