Predicting multiple conformations via sequence clustering and AlphaFold2.

AlphaFold2 (ref. 1) has revolutionized structural biology by accurately predicting single structures of proteins. However, a protein's biological function often depends on multiple conformational substates2, and disease-causing point mutations often cause population changes within these substates3,4. We demonstrate that clustering a multiple-sequence alignment by sequence similarity enables AlphaFold2 to sample alternative states of known metamorphic proteins with high confidence. Using this method, named AF-Cluster, we investigated the evolutionary distribution of predicted structures for the metamorphic protein KaiB5 and found that predictions of both conformations were distributed in clusters across the KaiB family. We used nuclear magnetic resonance spectroscopy to confirm an AF-Cluster prediction: a cyanobacteria KaiB variant is stabilized in the opposite state compared with the more widely studied variant. To test AF-Cluster's sensitivity to point mutations, we designed and experimentally verified a set of three mutations predicted to flip KaiB from Rhodobacter sphaeroides from the ground to the fold-switched state. Finally, screening for alternative states in protein families without known fold switching identified a putative alternative state for the oxidoreductase Mpt53 in Mycobacterium tuberculosis. Further development of such bioinformatic methods in tandem with experiments will probably have a considerable impact on predicting protein energy landscapes, essential for illuminating biological function.

From primordial clocks to circadian oscillators.

Circadian rhythms play an essential part in many biological processes, and only three prokaryotic proteins are required to constitute a true post-translational circadian oscillator1. The evolutionary history of the three Kai proteins indicates that KaiC is the oldest member and a central component of the clock2. Subsequent additions of KaiB and KaiA regulate the phosphorylation state of KaiC for time synchronization. The canonical KaiABC system in cyanobacteria is well understood3-6, but little is known about more ancient systems that only possess KaiBC. However, there are reports that they might exhibit a basic, hourglass-like timekeeping mechanism7-9. Here we investigate the primordial circadian clock in Rhodobacter sphaeroides, which contains only KaiBC, to elucidate its inner workings despite missing KaiA. Using a combination of X-ray crystallography and cryogenic electron microscopy, we find a new dodecameric fold for KaiC, in which two hexamers are held together by a coiled-coil bundle of 12 helices. This interaction is formed by the carboxy-terminal extension of KaiC and serves as an ancient regulatory moiety that is later superseded by KaiA. A coiled-coil register shift between daytime and night-time conformations is connected to phosphorylation sites through a long-range allosteric network that spans over 140 Å. Our kinetic data identify the difference in the ATP-to-ADP ratio between day and night as the environmental cue that drives the clock. They also unravel mechanistic details that shed light on the evolution of self-sustained oscillators.

A prolyl-isomerase mediates dopamine-dependent plasticity and cocaine motor sensitization.

Synaptic plasticity induced by cocaine and other drugs underlies addiction. Here we elucidate molecular events at synapses that cause this plasticity and the resulting behavioral response to cocaine in mice. In response to D1-dopamine-receptor signaling that is induced by drug administration, the glutamate-receptor protein metabotropic glutamate receptor 5 (mGluR5) is phosphorylated by microtubule-associated protein kinase (MAPK), which we show potentiates Pin1-mediated prolyl-isomerization of mGluR5 in instances where the product of an activity-dependent gene, Homer1a, is present to enable Pin1-mGluR5 interaction. These biochemical events potentiate N-methyl-D-aspartate receptor (NMDAR)-mediated currents that underlie synaptic plasticity and cocaine-evoked motor sensitization as tested in mice with relevant mutations. The findings elucidate how a coincidence of signals from the nucleus and the synapse can render mGluR5 accessible to activation with consequences for drug-induced dopamine responses and point to depotentiation at corticostriatal synapses as a possible therapeutic target for treating addiction.

Structure determination of high-energy states in a dynamic protein ensemble.

Macromolecular function frequently requires that proteins change conformation into high-energy states1-4. However, methods for solving the structures of these functionally essential, lowly populated states are lacking. Here we develop a method for high-resolution structure determination of minorly populated states by coupling NMR spectroscopy-derived pseudocontact shifts5 (PCSs) with Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion6 (PCS-CPMG). Our approach additionally defines the corresponding kinetics and thermodynamics of high-energy excursions, thereby characterizing the entire free-energy landscape. Using a large set of simulated data for adenylate kinase (Adk), calmodulin and Src kinase, we find that high-energy PCSs accurately determine high-energy structures (with a root mean squared deviation of less than 3.5 angström). Applying our methodology to Adk during catalysis, we find that the high-energy excursion involves surprisingly small openings of the AMP and ATP lids. This previously unresolved high-energy structure solves a longstanding controversy about conformational interconversions that are rate-limiting for catalysis. Primed for either substrate binding or product release, the high-energy structure of Adk suggests a two-step mechanism combining conformational selection to this state, followed by an induced-fit step into a fully closed state for catalysis of the phosphoryl-transfer reaction. Unlike other methods for resolving high-energy states, such as cryo-electron microscopy and X-ray crystallography, our solution PCS-CPMG approach excels in cases involving domain rearrangements of smaller systems (less than 60 kDa) and populations as low as 0.5%, and enables the simultaneous determination of protein structure, kinetics and thermodynamics while proteins perform their function.

A biophysical framework for double-drugging kinases.

Selective orthosteric inhibition of kinases has been challenging due to the conserved active site architecture of kinases and emergence of resistance mutants. Simultaneous inhibition of distant orthosteric and allosteric sites, which we refer to as "double-drugging", has recently been shown to be effective in overcoming drug resistance. However, detailed biophysical characterization of the cooperative nature between orthosteric and allosteric modulators has not been undertaken. Here, we provide a quantitative framework for double-drugging of kinases employing isothermal titration calorimetry, Förster resonance energy transfer, coupled-enzyme assays, and X-ray crystallography. We discern positive and negative cooperativity for Aurora A kinase (AurA) and Abelson kinase (Abl) with different combinations of orthosteric and allosteric modulators. We find that a conformational equilibrium shift is the main principle governing cooperativity. Notably, for both kinases, we find a synergistic decrease of the required orthosteric and allosteric drug dosages when used in combination to inhibit kinase activities to clinically relevant inhibition levels. X-ray crystal structures of the double-drugged kinase complexes reveal the molecular principles underlying the cooperative nature of double-drugging AurA and Abl with orthosteric and allosteric inhibitors. Finally, we observe a fully closed conformation of Abl when bound to a pair of positively cooperative orthosteric and allosteric modulators, shedding light on the puzzling abnormality of previously solved closed Abl structures. Collectively, our data provide mechanistic and structural insights into rational design and evaluation of double-drugging strategies.

Transient non-native hydrogen bonds promote activation of a signaling protein.

Phosphorylation is a common mechanism for activating proteins within signaling pathways. Yet, the molecular transitions between the inactive and active conformational states are poorly understood. Here we quantitatively characterize the free-energy landscape of activation of a signaling protein, nitrogen regulatory protein C (NtrC), by connecting functional protein dynamics of phosphorylation-dependent activation to protein folding and show that only a rarely populated, pre-existing active conformation is energetically stabilized by phosphorylation. Using nuclear magnetic resonance (NMR) dynamics, we test an atomic scale pathway for the complex conformational transition, inferred from molecular dynamics simulations (Lei et al., 2009). The data show that the loss of native stabilizing contacts during activation is compensated by non-native transient atomic interactions during the transition. The results unravel atomistic details of native-state protein energy landscapes by expanding the knowledge about ground states to transition landscapes.

Cumulative mechanism of several major imatinib-resistant mutations in Abl kinase.

Despite the outstanding success of the cancer drug imatinib, one obstacle in prolonged treatment is the emergence of resistance mutations within the kinase domain of its target, Abl. We noticed that many patient-resistance mutations occur in the dynamic hot spots recently identified to be responsible for imatinib's high selectivity toward Abl. In this study, we provide an experimental analysis of the mechanism underlying drug resistance for three major resistance mutations (G250E, Y253F, and F317L). Our data settle controversies, revealing unexpected resistance mechanisms. The mutations alter the energy landscape of Abl in complex ways: increased kinase activity, altered affinity, and cooperativity for the substrates, and, surprisingly, only a modestly decreased imatinib affinity. Only under cellular adenosine triphosphate (ATP) concentrations, these changes cumulate in an order of magnitude increase in imatinib's half-maximal inhibitory concentration (IC50). These results highlight the importance of characterizing energy landscapes of targets and its changes by drug binding and by resistance mutations developed by patients.

Allosteric modulation of a human protein kinase with monobodies.

Despite being the subject of intense effort and scrutiny, kinases have proven to be consistently challenging targets in inhibitor drug design. A key obstacle has been promiscuity and consequent adverse effects of drugs targeting the ATP binding site. Here we introduce an approach to controlling kinase activity by using monobodies that bind to the highly specific regulatory allosteric pocket of the oncoprotein Aurora A (AurA) kinase, thereby offering the potential for more specific kinase modulators. Strikingly, we identify a series of highly specific monobodies acting either as strong kinase inhibitors or activators via differential recognition of structural motifs in the allosteric pocket. X-ray crystal structures comparing AurA bound to activating vs inhibiting monobodies reveal the atomistic mechanism underlying allosteric modulation. The results reveal 3 major advantages of targeting allosteric vs orthosteric sites: extreme selectivity, ability to inhibit as well as activate, and avoidance of competing with ATP that is present at high concentrations in the cells. We envision that exploiting allosteric networks for inhibition or activation will provide a general, powerful pathway toward rational drug design.

Antiparallel EmrE exports drugs by exchanging between asymmetric structures.

Small multidrug resistance transporters provide an ideal system to study the minimal requirements for active transport. EmrE is one such transporter in Escherichia coli. It exports a broad class of polyaromatic cation substrates, thus conferring resistance to drug compounds matching this chemical description. However, a great deal of controversy has surrounded the topology of the EmrE homodimer. Here we show that asymmetric antiparallel EmrE exchanges between inward- and outward-facing states that are identical except that they have opposite orientation in the membrane. We quantitatively measure the global conformational exchange between these two states for substrate-bound EmrE in bicelles using solution NMR dynamics experiments. Förster resonance energy transfer reveals that the monomers within each dimer are antiparallel, and paramagnetic relaxation enhancement NMR experiments demonstrate differential water accessibility of the two monomers within each dimer. Our experiments reveal a 'dynamic symmetry' that reconciles the asymmetric EmrE structure with the functional symmetry of residues in the active site.

A minor conformation of a lanthanide tag on adenylate kinase characterized by paramagnetic relaxation dispersion NMR spectroscopy.

NMR relaxation dispersion techniques provide a powerful method to study protein dynamics by characterizing lowly populated conformations that are in dynamic exchange with the major state. Paramagnetic NMR is a versatile tool for investigating the structures and dynamics of proteins. These two techniques were combined here to measure accurate and precise pseudocontact shifts of a lowly populated conformation. This method delivers valuable long-range structural restraints for higher energy conformations of macromolecules in solution. Another advantage of combining pseudocontact shifts with relaxation dispersion is the increase in the amplitude of dispersion profiles. Lowly populated states are often involved in functional processes, such as enzyme catalysis, signaling, and protein/protein interactions. The presented results also unveil a critical problem with the lanthanide tag used to generate paramagnetic relaxation dispersion effects in proteins, namely that the motions of the tag can interfere severely with the observation of protein dynamics. The two-point attached CLaNP-5 lanthanide tag was linked to adenylate kinase. From the paramagnetic relaxation dispersion only motion of the tag is observed. The data can be described accurately by a two-state model in which the protein-attached tag undergoes a 23° tilting motion on a timescale of milliseconds. The work demonstrates the large potential of paramagnetic relaxation dispersion and the challenge to improve current tags to minimize relaxation dispersion from tag movements.

Hidden alternative structures of proline isomerase essential for catalysis.

A long-standing challenge is to understand at the atomic level how protein dynamics contribute to enzyme catalysis. X-ray crystallography can provide snapshots of conformational substates sampled during enzymatic reactions, while NMR relaxation methods reveal the rates of interconversion between substates and the corresponding relative populations. However, these current methods cannot simultaneously reveal the detailed atomic structures of the rare states and rationalize the finding that intrinsic motions in the free enzyme occur on a timescale similar to the catalytic turnover rate. Here we introduce dual strategies of ambient-temperature X-ray crystallographic data collection and automated electron-density sampling to structurally unravel interconverting substates of the human proline isomerase, cyclophilin A (CYPA, also known as PPIA). A conservative mutation outside the active site was designed to stabilize features of the previously hidden minor conformation. This mutation not only inverts the equilibrium between the substates, but also causes large, parallel reductions in the conformational interconversion rates and the catalytic rate. These studies introduce crystallographic approaches to define functional minor protein conformations and, in combination with NMR analysis of the enzyme dynamics in solution, show how collective motions directly contribute to the catalytic power of an enzyme.

A highly conserved cysteine of neuronal calcium-sensing proteins controls cooperative binding of Ca2+ to recoverin.

Recoverin, a 23-kDa Ca(2+)-binding protein of the neuronal calcium sensing (NCS) family, inhibits rhodopsin kinase, a Ser/Thr kinase responsible for termination of photoactivated rhodopsin in rod photoreceptor cells. Recoverin has two functional EF hands and a myristoylated N terminus. The myristoyl chain imparts cooperativity to the Ca(2+)-binding sites through an allosteric mechanism involving a conformational equilibrium between R and T states of the protein. Ca(2+) binds preferentially to the R state; the myristoyl chain binds preferentially to the T state. In the absence of myristoylation, the R state predominates, and consequently, binding of Ca(2+) to the non-myristoylated protein is not cooperative. We show here that a mutation, C39A, of a highly conserved Cys residue among NCS proteins, increases the apparent cooperativity for binding of Ca(2+) to non-myristoylated recoverin. The binding data can be explained by an effect on the T/R equilibrium to favor the T state without affecting the intrinsic binding constants for the two Ca(2+) sites.

The conformational landscape of fold-switcher KaiB is tuned to the circadian rhythm timescale.

How can a single protein domain encode a conformational landscape with multiple stably-folded states, and how do those states interconvert? Here, we use real-time and relaxation-dispersion NMR to characterize the conformational landscape of the circadian rhythm protein KaiB from Rhodobacter sphaeroides. Unique among known natural metamorphic proteins, this KaiB variant spontaneously interconverts between two monomeric states: the "Ground" and "Fold-switched" (FS) state. KaiB in its FS state interacts with multiple binding partners, including the central KaiC protein, to regulate circadian rhythms. We find that KaiB itself takes hours to interconvert between the Ground and FS state, underscoring the ability of a single sequence to encode the slow process needed for function. We reveal the rate-limiting step between the Ground and FS state is the cis-trans isomerization of three prolines in the fold-switching region by demonstrating interconversion acceleration by the prolyl isomerase CypA. The interconversion proceeds through a "partially disordered" (PD) state, where the C-terminal half becomes disordered while the N-terminal half remains stably folded. We discovered two additional properties of KaiB's landscape. Firstly, the Ground state experiences cold denaturation: at 4°C, the PD state becomes the majorly populated state. Secondly, the Ground state exchanges with a fourth state, the "Enigma" state, on the millisecond timescale. We combine AlphaFold2-based predictions and NMR chemical shift predictions to predict this "Enigma" state is a beta-strand register shift that eases buried charged residues, and support this structure experimentally. These results provide mechanistic insight in how evolution can design a single sequence that achieves specific timing needed for its function.

Dynamic personalities of proteins.

Because proteins are central to cellular function, researchers have sought to uncover the secrets of how these complex macromolecules execute such a fascinating variety of functions. Although static structures are known for many proteins, the functions of proteins are governed ultimately by their dynamic character (or 'personality'). The dream is to 'watch' proteins in action in real time at atomic resolution. This requires addition of a fourth dimension, time, to structural biology so that the positions in space and time of all atoms in a protein can be described in detail.

A hierarchy of timescales in protein dynamics is linked to enzyme catalysis.

The synergy between structure and dynamics is essential to the function of biological macromolecules. Thermally driven dynamics on different timescales have been experimentally observed or simulated, and a direct link between micro- to milli-second domain motions and enzymatic function has been established. However, very little is understood about the connection of these functionally relevant, collective movements with local atomic fluctuations, which are much faster. Here we show that pico- to nano-second timescale atomic fluctuations in hinge regions of adenylate kinase facilitate the large-scale, slower lid motions that produce a catalytically competent state. The fast, local mobilities differ between a mesophilic and hyperthermophilic adenylate kinase, but are strikingly similar at temperatures at which enzymatic activity and free energy of folding are matched. The connection between different timescales and the corresponding amplitudes of motions in adenylate kinase and their linkage to catalytic function is likely to be a general characteristic of protein energy landscapes.

Intrinsic motions along an enzymatic reaction trajectory.

The mechanisms by which enzymes achieve extraordinary rate acceleration and specificity have long been of key interest in biochemistry. It is generally recognized that substrate binding coupled to conformational changes of the substrate-enzyme complex aligns the reactive groups in an optimal environment for efficient chemistry. Although chemical mechanisms have been elucidated for many enzymes, the question of how enzymes achieve the catalytically competent state has only recently become approachable by experiment and computation. Here we show crystallographic evidence for conformational substates along the trajectory towards the catalytically competent 'closed' state in the ligand-free form of the enzyme adenylate kinase. Molecular dynamics simulations indicate that these partially closed conformations are sampled in nanoseconds, whereas nuclear magnetic resonance and single-molecule fluorescence resonance energy transfer reveal rare sampling of a fully closed conformation occurring on the microsecond-to-millisecond timescale. Thus, the larger-scale motions in substrate-free adenylate kinase are not random, but preferentially follow the pathways that create the configuration capable of proficient chemistry. Such preferred directionality, encoded in the fold, may contribute to catalysis in many enzymes.

Intrinsic dynamics of an enzyme underlies catalysis.

A unique feature of chemical catalysis mediated by enzymes is that the catalytically reactive atoms are embedded within a folded protein. Although current understanding of enzyme function has been focused on the chemical reactions and static three-dimensional structures, the dynamic nature of proteins has been proposed to have a function in catalysis. The concept of conformational substates has been described; however, the challenge is to unravel the intimate linkage between protein flexibility and enzymatic function. Here we show that the intrinsic plasticity of the protein is a key characteristic of catalysis. The dynamics of the prolyl cis-trans isomerase cyclophilin A (CypA) in its substrate-free state and during catalysis were characterized with NMR relaxation experiments. The characteristic enzyme motions detected during catalysis are already present in the free enzyme with frequencies corresponding to the catalytic turnover rates. This correlation suggests that the protein motions necessary for catalysis are an intrinsic property of the enzyme and may even limit the overall turnover rate. Motion is localized not only to the active site but also to a wider dynamic network. Whereas coupled networks in proteins have been proposed previously, we experimentally measured the collective nature of motions with the use of mutant forms of CypA. We propose that the pre-existence of collective dynamics in enzymes before catalysis is a common feature of biocatalysts and that proteins have evolved under synergistic pressure between structure and dynamics.

Probing microsecond time scale dynamics in proteins by methyl (1)H Carr-Purcell-Meiboom-Gill relaxation dispersion NMR measurements. Application to activation of the signaling protein NtrC(r).

To study microsecond processes by relaxation dispersion NMR spectroscopy, low power deposition and short pulses are crucial and encourage the development of experiments that employ (1)H Carr-Purcell-Meiboom-Gill (CPMG) pulse trains. Herein, a method is described for the comprehensive study of microsecond to millisecond time scale dynamics of methyl groups in proteins, exploiting their high abundance and favorable relaxation properties. In our approach, protein samples are produced using [(1)H, (13)C]-d-glucose in ∼100% D(2)O, which yields CHD(2) methyl groups for alanine, valine, threonine, isoleucine, leucine, and methionine residues with high abundance, in an otherwise largely deuterated background. Methyl groups in such samples can be sequence-specifically assigned to near completion, using (13)C TOCSY NMR spectroscopy, as was recently demonstrated (Otten, R.; et al. J. Am. Chem. Soc. 2010, 132, 2952-2960). In this Article, NMR pulse schemes are presented to measure (1)H CPMG relaxation dispersion profiles for CHD(2) methyl groups, in a vein similar to that of backbone relaxation experiments. Because of the high deuteration level of methyl-bearing side chains, artifacts arising from proton scalar coupling during the CPMG pulse train are negligible, with the exception of Ile-δ1 and Thr-γ2 methyl groups, and a pulse scheme is described to remove the artifacts for those residues. Strong (13)C scalar coupling effects, observed for several leucine residues, are removed by alternative biochemical and NMR approaches. The methodology is applied to the transcriptional activator NtrC(r), for which an inactive/active state transition was previously measured and the motions in the microsecond time range were estimated through a combination of backbone (15)N CPMG dispersion NMR spectroscopy and a collection of experiments to determine the exchange-free component to the transverse relaxation rate. Exchange contributions to the (1)H line width were detected for 21 methyl groups, and these probes were found to collectively report on a local structural rearrangement around the phosphorylation site, with a rate constant of (15.5 ± 0.5) × 10(3) per second (i.e., τ(ex) = 64.7 ± 1.9 µs). The affected methyl groups indicate that, already before phosphorylation, a substantial, transient rearrangement takes place between helices 3 and 4 and strands 4 and 5. This conformational equilibrium allows the protein to gain access to the active, signaling state in the absence of covalent modification through a shift in a pre-existing dynamic equilibrium. Moreover, the conformational switching maps exactly to the regions that differ between the solution NMR structures of the fully inactive and active states. These results demonstrate that a cost-effective and quantitative study of protein methyl group dynamics by (1)H CPMG relaxation dispersion NMR spectroscopy is possible and can be applied to study functional motions on the microsecond time scale that cannot be accessed by backbone (15)N relaxation dispersion NMR. The use of methyl groups as dynamics probes extends such applications also to larger proteins.

Response to Comment on "Ancient origins of allosteric activation in a Ser-Thr kinase".

Park et al question one out of seven findings from Hadzipasic et al: whether TPX2 allosterically regulates the oldest Aurora. We had already addressed the two concerns raised-sparse sequence sampling and not forcing the gene to the species tree-before publication. Moreover, we believe their ancestral sequence reconstruction would be consistent with a nonallosteric common ancestor, and we show large sequence differences caused by species tree-enforced gene trees.