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.

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.

Comprehensive analysis of NMR data using advanced line shape fitting.

NMR spectroscopy is uniquely suited for atomic resolution studies of biomolecules such as proteins, nucleic acids and metabolites, since detailed information on structure and dynamics are encoded in positions and line shapes of peaks in NMR spectra. Unfortunately, accurate determination of these parameters is often complicated and time consuming, in part due to the need for different software at the various analysis steps and for validating the results. Here, we present an integrated, cross-platform and open-source software that is significantly more versatile than the typical line shape fitting application. The software is a completely redesigned version of PINT ( https://pint-nmr.github.io/PINT/ ). It features a graphical user interface and includes functionality for peak picking, editing of peak lists and line shape fitting. In addition, the obtained peak intensities can be used directly to extract, for instance, relaxation rates, heteronuclear NOE values and exchange parameters. In contrast to most available software the entire process from spectral visualization to preparation of publication-ready figures is done solely using PINT and often within minutes, thereby, increasing productivity for users of all experience levels. Unique to the software are also the outstanding tools for evaluating the quality of the fitting results and extensive, but easy-to-use, customization of the fitting protocol and graphical output. In this communication, we describe the features of the new version of PINT and benchmark its performance.

Subpolar addition of new cell wall is directed by DivIVA in mycobacteria.

Mycobacteria are surrounded by a complex multilayered envelope and elongate at the poles. The principles that organize the coordinated addition of chemically diverse cell wall layers during polar extension remain unclear. We show that enzymes mediating the terminal cytosolic steps of peptidoglycan, arabinogalactan, and mycolic acid synthesis colocalize at sites of cell growth or division. The tropomyosin-like protein, DivIVA, is targeted to the negative curvature of the pole, is enriched at the growing end, and determines cell shape from this site. In contrast, cell wall synthetic complexes are concentrated at a distinct subpolar location. When viewed at subdiffraction resolution, new peptidoglycan is deposited at this subpolar site, and inert cell wall covers the DivIVA-marked tip. The differentiation between polar tip and cell wall synthetic complexes is also apparent at the biochemical level. Enzymes that generate mycolate precursors interact with DivIVA, but the final condensation of mycolic acids occurs in a distinct protein complex at the site of nascent cell wall addition. We propose an ultrastructural model of mycobacterial polar growth where new cell wall is added in an annular zone below the cell tip. This model may be broadly applicable to other bacterial and fungal organisms that grow via polar extension.

Cavity as a source of conformational fluctuation and high-energy state: high-pressure NMR study of a cavity-enlarged mutant of T4 lysozyme.

Although the structure, function, conformational dynamics, and controlled thermodynamics of proteins are manifested by their corresponding amino acid sequences, the natural rules for molecular design and their corresponding interplay remain obscure. In this study, we focused on the role of internal cavities of proteins in conformational dynamics. We investigated the pressure-induced responses from the cavity-enlarged L99A mutant of T4 lysozyme, using high-pressure NMR spectroscopy. The signal intensities of the methyl groups in the (1)H/(13)C heteronuclear single quantum correlation spectra, particularly those around the enlarged cavity, decreased with the increasing pressure, and disappeared at 200 MPa, without the appearance of new resonances, thus indicating the presence of heterogeneous conformations around the cavity within the ground state ensemble. Above 200 MPa, the signal intensities of >20 methyl groups gradually decreased with the increasing pressure, without the appearance of new resonances. Interestingly, these residues closely matched those sensing a large conformational change between the ground- and high-energy states, at atmospheric pressure. (13)C and (1)H NMR line-shape simulations showed that the pressure-induced loss in the peak intensity could be explained by the increase in the high-energy state population. In this high-energy state, the aromatic side chain of F114 gets flipped into the enlarged cavity. The accommodation of the phenylalanine ring into the efficiently packed cavity may decrease the partial molar volume of the high-energy state, relative to the ground state. We suggest that the enlarged cavity is involved in the conformational transition to high-energy states and in the volume fluctuation of the ground state.

The solution structure, binding properties, and dynamics of the bacterial siderophore-binding protein FepB.

The periplasmic binding protein (PBP) FepB plays a key role in transporting the catecholate siderophore ferric enterobactin from the outer to the inner membrane in Gram-negative bacteria. The solution structures of the 34-kDa apo- and holo-FepB from Escherichia coli, solved by NMR, represent the first solution structures determined for the type III class of PBPs. Unlike type I and II PBPs, which undergo large "Venus flytrap" conformational changes upon ligand binding, both forms of FepB maintain similar overall folds; however, binding of the ligand is accompanied by significant loop movements. Reverse methyl cross-saturation experiments corroborated chemical shift perturbation results and uniquely defined the binding pocket for gallium enterobactin (GaEnt). NMR relaxation experiments indicated that a flexible loop (residues 225-250) adopted a more rigid and extended conformation upon ligand binding, which positioned residues for optimal interactions with the ligand and the cytoplasmic membrane ABC transporter (FepCD), respectively. In conclusion, this work highlights the pivotal role that structural dynamics plays in ligand binding and transporter interactions in type III PBPs.

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.

Optimized co-solute paramagnetic relaxation enhancement for the rapid NMR analysis of a highly fibrillogenic peptide.

Co-solute paramagnetic relaxation enhancement (PRE) is an attractive way to speed up data acquisition in NMR spectroscopy by shortening the T 1 relaxation time of the nucleus of interest and thus the necessary recycle delay. Here, we present the rationale to utilize high-spin iron(III) as the optimal transition metal for this purpose and characterize the properties of its neutral chelate form Fe(DO3A) as a suitable PRE agent. Fe(DO3A) effectively reduces the T 1 values across the entire sequence of the intrinsically disordered protein α-synuclein with negligible impact on line width. The agent is better suited than currently used alternatives, shows no specific interaction with the polypeptide chain and, due to its high relaxivity, is effective at low concentrations and in 'proton-less' NMR experiments. By using Fe(DO3A) we were able to complete the backbone resonance assignment of a highly fibrillogenic peptide from α1-antitrypsin by acquiring the necessary suite of multidimensional NMR datasets in 3 h.

Fluoride resistance and transport by riboswitch-controlled CLC antiporters.

A subclass of bacterial CLC anion-transporting proteins, phylogenetically distant from long-studied CLCs, was recently shown to be specifically up-regulated by F(-). We establish here that a set of randomly selected representatives from this "CLC(F)" clade protect Escherichia coli from F(-) toxicity, and that the purified proteins catalyze transport of F(-) in liposomes. Sequence alignments and membrane transport experiments using (19)F NMR, osmotic response assays, and planar lipid bilayer recordings reveal four mechanistic traits that set CLC(F) proteins apart from all other known CLCs. First, CLC(F)s lack conserved residues that form the anion binding site in canonical CLCs. Second, CLC(F)s exhibit high anion selectivity for F(-) over Cl(-). Third, at a residue thought to distinguish CLC channels and transporters, CLC(F)s bear a channel-like valine rather than a transporter-like glutamate, and yet are F(-)/H(+) antiporters. Finally, F(-)/H(+) exchange occurs with 1:1 stoichiometry, in contrast to the usual value of 2:1.

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.

Protein surface and core dynamics show concerted hydration-dependent activation.

By specifically labeling leucine/valine methyl groups and lysine side chains "inside" and "outside" dynamics of proteins on the nanosecond timescale are compared using neutron scattering. Surprisingly, both groups display similar dynamics as a function of temperature, and the buried hydrophobic core is sensitive to hydration and undergoes a dynamical transition.

Cystathionine ß-synthase (CBS) domains 1 and 2 fulfill different roles in ionic strength sensing of the ATP-binding cassette (ABC) transporter OpuA.

The cystathionine ß-synthase module of OpuA in conjunction with an anionic membrane surface acts as a sensor of internal ionic strength, which allows the protein to respond to osmotic stress. We now show by chemical modification and cross-linking studies that CBS2-CBS2 interface residues are critical for transport activity and/or ionic regulation of transport, whereas CBS1 serves no functional role. We establish that Cys residues in CBS1, CBS2, and the nucleotide-binding domain are more accessible for cross-linking at high than low ionic strength, indicating that these domains undergo conformational changes when transiting between the active and inactive state. Structural analyses suggest that the cystathionine ß-synthase module is largely unstructured. Moreover, we could substitute CBS1 by a linker and preserve ionic regulation of transport. These data suggest that CBS1 serves as a linker and the structured CBS2-CBS2 interface forms a hinge point for ionic strength-dependent rearrangements that are transmitted to the nucleotide-binding domain and thereby affect translocation activity.

Origin and removal of mixed-phase artifacts in gradient sensitivity enhanced heteronuclear single quantum correlation spectra.

Here we describe phasing anomalies observed in gradient sensitivity enhanced (15)N-(1)H HSQC spectra, and analyze their origin. It is shown that, as a result of (15)N off-resonance effects, dispersive contributions to the (1)H signal become detectable, and lead to (15)N-offset dependent phase errors. Strategies that effectively suppress these artifacts are presented.

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.

Comprehensive and cost-effective NMR spectroscopy of methyl groups in large proteins.

An NMR approach is described which yields the methyl resonance assignments of alanine, threonine, valine, leucine, and isoleucine residues in proteins with high sensitivity and excellent resolution. The method relies on protein samples produced by bacterial expression using [(1)H,(13)C]-D-glucose and approximately 100% D(2)O, which is cost-effective and ensures the isotopic enrichment of all possible methyl groups. Magnetization transfer throughout the methyl-containing side chains is possible with this labeling scheme due to the high level of deuteration along the amino acid side chain, coupled with the selection of the favorable CHD(2) methyl isotopomer for detection. In an application to the 34 kDa periplasmic binding protein FepB 164 out of 195 methyl groups (85%) were assigned sequence-specifically and stereospecifically. This percentage increases to 91% when taking into account that not all backbone assignments are available for this system. The remaining unassigned methyl groups belong to six leucine residues, caused by low cross-peak intensities, and four alanine residues due to degeneracy of the (13)C(alpha)/(13)C(beta) frequencies. Our results demonstrate that NMR spectroscopic investigations of protein structure, dynamics, and interactions can be extended to include all methyl-containing amino acids also for larger proteins.

How directed evolution reshapes the energy landscape in an enzyme to boost catalysis.

The advent of biocatalysts designed computationally and optimized by laboratory evolution provides an opportunity to explore molecular strategies for augmenting catalytic function. Applying a suite of nuclear magnetic resonance, crystallography, and stopped-flow techniques to an enzyme designed for an elementary proton transfer reaction, we show how directed evolution gradually altered the conformational ensemble of the protein scaffold to populate a narrow, highly active conformational ensemble and accelerate this transformation by nearly nine orders of magnitude. Mutations acquired during optimization enabled global conformational changes, including high-energy backbone rearrangements, that cooperatively organized the catalytic base and oxyanion stabilizer, thus perfecting transition-state stabilization. The development of protein catalysts for many chemical transformations could be facilitated by explicitly sampling conformational substates during design and specifically stabilizing productive substates over all unproductive conformations.

Comprehensive determination of (3)J (HNHalpha) for unfolded proteins using (13)C'-resolved spin-echo difference spectroscopy.

An experiment is presented to determine (3)J(HNHalpha) coupling constants, with significant advantages for applications to unfolded proteins. The determination of coupling constants for the peptide chain using 1D (1)H, or 2D and 3D (1)H-(15)N correlation spectroscopy is often hampered by extensive resonance overlap when dealing with flexible, disordered proteins. In the experiment detailed here, the overlap problem is largely circumvented by recording (1)H-(13)C' correlation spectra, which demonstrate superior resolution for unfolded proteins. J-coupling constants are extracted from the peak intensities in a pair of 2D spin-echo difference experiments, affording rapid acquisition of the coupling data. In an application to the cytoplasmic domain of human neuroligin-3 (hNlg3cyt) data were obtained for 78 residues, compared to 54 coupling constants obtained from a 3D HNHA experiment. The coupling constants suggest that hNlg3cyt is intrinsically disordered, with little propensity for structure.

Probing the Transition State in Enzyme Catalysis by High-Pressure NMR Dynamics.

Protein conformational changes are frequently essential for enzyme catalysis, and in several cases, shown to be the limiting factor for overall catalytic speed. However, a structural understanding of corresponding transition states, needed to rationalize the kinetics, remains obscure due to their fleeting nature. Here, we determine the transition-state ensemble of the rate-limiting conformational transition in the enzyme adenylate kinase, by a synergistic approach between experimental high-pressure NMR relaxation during catalysis and molecular dynamics simulations. By comparing homologous kinases evolved under ambient or high pressure in the deep-sea, we detail transition state ensembles that differ in solvation as directly measured by the pressure dependence of catalysis. Capturing transition-state ensembles begins to complete the catalytic energy landscape that is generally characterized by structures of all intermediates and frequencies of transitions among them.