Experimental NOE, Chemical Shift, and Proline Isomerization Data Provide Detailed Insights into Amelotin Oligomerization.

Amelotin is an intrinsically disordered protein (IDP) rich in Pro residues and is involved in hydroxyapatite mineralization. It rapidly oligomerizes under physiological conditions of pH and pressure but reverts to its monomeric IDP state at elevated pressure. We identified a 105-residue segment of the protein that becomes ordered upon oligomerization, and we used pressure-jump NMR spectroscopy to measure long-range NOE contacts that exist exclusively in the oligomeric NMR-invisible state. The kinetics of oligomerization and dissociation were probed at the residue-specific level, revealing that the oligomerization process is initiated in the C-terminal half of the segment. Using pressure-jump NMR, the degree of order in the oligomer at the sites of Pro residues was probed by monitoring changes in cis/trans equilibria relative to the IDP state after long-term equilibration under oligomerizing conditions. Whereas most Pro residues revert to trans in the oligomeric state, Pro-49 favors a cis configuration and three Pro residues retain an unchanged cis fraction, pointing to their local lack of order in the oligomeric state. NOE contacts and secondary 13C chemical shifts in the oligomeric state indicate the presence of an 11-residue α-helix, preceded by a small intramolecular antiparallel ß-sheet, with slower formation of long-range intermolecular interactions to N-terminal residues. Although none of the models generated by AlphaFold2 for the amelotin monomer was consistent with experimental data, subunits of a hexamer generated by AlphaFold-Multimer satisfied intramolecular NOE and chemical shift data and may provide a starting point for developing atomic models for the oligomeric state.

NMR Observation of Sulfhydryl Signals in SARS-CoV-2 Main Protease Aids Structural Studies.

The 68-kDa homodimeric 3C-like protease of SARS-CoV-2, Mpro (3CLpro /Nsp5), is a key antiviral drug target. NMR spectroscopy of this large system proved challenging and resonance assignments have remained incomplete. Here we present the near-complete (>97 %) backbone assignments of a C145A variant of Mpro (Mpro C145A ) both with, and without, the N-terminal auto-cleavage substrate sequence, in its native homodimeric state. We also present SILLY (Selective Inversion of thioL and Ligand for NOESY), a simple yet effective pseudo-3D NMR experiment that utilizes NOEs to identify interactions between Cys-thiol or aliphatic protons, and their spatially proximate backbone amides in a perdeuterated protein background. High protection against hydrogen exchange is observed for 10 of the 11 thiol groups in Mpro C145A , even those that are partially accessible to solvent. A combination of SILLY methods and high-resolution triple-resonance NMR experiments reveals site-specific interactions between Mpro , its substrate peptides, and other ligands, which present opportunities for competitive binding studies in future drug design efforts.

Concordance of X-ray and AlphaFold2 Models of SARS-CoV-2 Main Protease with Residual Dipolar Couplings Measured in Solution.

The 68-kDa homodimeric 3C-like protease of SARS-CoV-2, Mpro (3CLpro/Nsp5), is a promising antiviral drug target. We evaluate the concordance of models generated by the newly introduced AlphaFold2 structure prediction program with residual dipolar couplings (RDCs) measured in solution for 15N-1HN and 13C'-1HN atom pairs. The latter were measured using a new, highly precise TROSY-AntiTROSY Encoded RDC (TATER) experiment. Three sets of AlphaFold2 models were evaluated: (1) MproAF, generated using the standard AlphaFold2 input structural database; (2) MproAFD, where the AlphaFold2 implementation was modified to exclude all candidate template X-ray structures deposited after Jan 1, 2020; and (3) MproAFS, which excluded all structures homologous to coronaviral Mpro. Close agreement between all three sets of AlphaFold models and experimental RDC data is found for most of the protein. For residues in well-defined secondary structure, the agreement decreases somewhat upon Amber relaxation. For these regions, MproAF agreement exceeds that of most high-resolution X-ray structures. Residues from domain 2 that comprise elements of both the active site and the homo-dimerization interface fit less well across all structures. These results indicate novel opportunities for combining experimentation with molecular dynamics simulations, where solution RDCs provide highly precise input for QM/MM simulations of substrate binding/reaction trajectories.

Synthesis and applications of RNAs with position-selective labelling and mosaic composition.

Knowledge of the structure and dynamics of RNA molecules is critical to understanding their many biological functions. Furthermore, synthetic RNAs have applications as therapeutics and molecular sensors. Both research and technological applications of RNA would be dramatically enhanced by methods that enable incorporation of modified or labelled nucleotides into specifically designated positions or regions of RNA. However, the synthesis of tens of milligrams of such RNAs using existing methods has been impossible. Here we develop a hybrid solid-liquid phase transcription method and automated robotic platform for the synthesis of RNAs with position-selective labelling. We demonstrate its use by successfully preparing various isotope- or fluorescently labelled versions of the 71-nucleotide aptamer domain of an adenine riboswitch for nuclear magnetic resonance spectroscopy or single-molecule Förster resonance energy transfer, respectively. Those RNAs include molecules that were selectively isotope-labelled in specific loops, linkers, a helix, several discrete positions, or a single internal position, as well as RNA molecules that were fluorescently labelled in and near kissing loops. These selectively labelled RNAs have the same fold as those transcribed using conventional methods, but they greatly simplify the interpretation of NMR spectra. The single-position isotope- and fluorescently labelled RNA samples reveal multiple conformational states of the adenine riboswitch. Lastly, we describe a robotic platform and the operation that automates this technology. Our selective labelling method may be useful for studying RNA structure and dynamics and for making RNA sensors for a variety of applications including cell-biological studies, substance detection, and disease diagnostics.

Study of protein folding under native conditions by rapidly switching the hydrostatic pressure inside an NMR sample cell.

In general, small proteins rapidly fold on the timescale of milliseconds or less. For proteins with a substantial volume difference between the folded and unfolded states, their thermodynamic equilibrium can be altered by varying the hydrostatic pressure. Using a pressure-sensitized mutant of ubiquitin, we demonstrate that rapidly switching the pressure within an NMR sample cell enables study of the unfolded protein under native conditions and, vice versa, study of the native protein under denaturing conditions. This approach makes it possible to record 2D and 3D NMR spectra of the unfolded protein at atmospheric pressure, providing residue-specific information on the folding process. 15N and 13C chemical shifts measured immediately after dropping the pressure from 2.5 kbar (favoring unfolding) to 1 bar (native) are close to the random-coil chemical shifts observed for a large, disordered peptide fragment of the protein. However, 15N relaxation data show evidence for rapid exchange, on a ∼100-µs timescale, between the unfolded state and unstable, structured states that can be considered as failed folding events. The NMR data also provide direct evidence for parallel folding pathways, with approximately one-half of the protein molecules efficiently folding through an on-pathway kinetic intermediate, whereas the other half fold in a single step. At protein concentrations above ∼300 µM, oligomeric off-pathway intermediates compete with folding of the native state.

Modulating the Stiffness of the Myosin VI Single α-Helical Domain.

Highly charged, single α-helical (SAH) domains contain a high percentage of Arg, Lys, and Glu residues. Their dynamic salt bridge pairing creates the exceptional stiffness of these helical rods, with a persistence length of more than 200 Å for the myosin VI SAH domain. With the aim of modulating the stiffness of the helical structure, we investigated the effect, using NMR spectroscopy, of substituting key charged Arg, Lys, Glu, and Asp residues by Gly or His. Results indicate that such mutations result in the transient breaking of the helix at the site of mutation but with noticeable impact on amide hydrogen exchange rates extending as far as ±2 helical turns, pointing to a substantial degree of cooperativity in SAH stability. Whereas a single Gly substitution caused transient breaks ∼20% of the time, two consecutive Gly substitutions break the helix ∼65% of the time. NMR relaxation measurements indicate that the exchange rate between an intact and a broken helix is fast (>300,000 s-1) and that for the wild-type sequence, the finite persistence length is dominated by thermal fluctuations of backbone torsion angles and H-bond lengths, not by transient helix breaking. The double mutation D27H/E28H causes a pH-dependent fraction of helix disruption, in which the helix breakage increases from 26% at pH 7.5 to 53% at pH 5.5. The ability to modulate helical integrity by pH may enable incorporation of externally tunable dynamic components in the design of molecular machines.

Remarkable Rigidity of the Single α-Helical Domain of Myosin-VI As Revealed by NMR Spectroscopy.

Although the α-helix has long been recognized as an all-important element of secondary structure, it generally requires stabilization by tertiary interactions with other parts of a protein's structure. Highly charged single α-helical (SAH) domains, consisting of a high percentage (>75%) of Arg, Lys, and Glu residues, are exceptions to this rule but have been difficult to characterize structurally. Our study focuses on the 68-residue medial tail domain of myosin-VI, which is found to contain a highly ordered α-helical structure extending from Glu-6 to Lys-63. High hydrogen exchange protection factors (15-150), small (ca. 4 Hz) 3 JHNHα couplings, and a near-perfect fit to an ideal model α-helix for its residual dipolar couplings (RDCs), measured in a filamentous phage medium, support the high regularity of this helix. Remarkably, the hydrogen exchange rates are far more homogeneous than the protection factors derived from them, suggesting that for these transiently broken helices the intrinsic exchange rates derived from the amino acid sequence are not appropriate reference values. 15N relaxation data indicate a very high degree of rotational diffusion anisotropy ( D∥/ D⊥ ≈ 7.6), consistent with the hydrodynamic behavior predicted for such a long, nearly straight α-helix. Alignment of the helix by a paramagnetic lanthanide ion attached to its N-terminal region shows a decrease in alignment as the distance from the tagging site increases. This decrease yields a precise measure for the persistence length of 224 ± 10 Å at 20 °C, supporting the idea that the role of the SAH helix is to act as an extension of the myosin-VI lever arm.

Importance of time-ordered non-uniform sampling of multi-dimensional NMR spectra of Aß1-42 peptide under aggregating conditions.

Although the order of the time steps in which the non-uniform sampling (NUS) schedule is implemented when acquiring multi-dimensional NMR spectra is of limited importance when sample conditions remain unchanged over the course of the experiment, it is shown to have major impact when samples are unstable. In the latter case, time-ordering of the NUS data points by the normalized radial length yields a reduction of sampling artifacts, regardless of the spectral reconstruction algorithm. The disadvantage of time-ordered NUS sampling is that halting the experiment prior to its completion will result in lower spectral resolution, rather than a sparser data matrix. Alternatively, digitally correcting for sample decay prior to reconstruction of randomly ordered NUS data points can mitigate reconstruction artifacts, at the cost of somewhat lower sensitivity. Application of these sampling schemes to the Alzheimer's amyloid beta (Aß1-42) peptide at an elevated concentration, low temperature, and 3 kbar of pressure, where approximately 75% of the peptide reverts to an NMR-invisible state during the collection of a 3D 15N-separated NOESY spectrum, highlights the improvement in artifact suppression and reveals weak medium-range NOE contacts in several regions, including the C-terminal region of the peptide.

Propensity for cis-Proline Formation in Unfolded Proteins.

In unfolded proteins, peptide bonds involving Pro residues exist in equilibrium between the minor cis and major trans conformations. Folded proteins predominantly contain trans-Pro bonds, and slow cis-trans Pro isomerization in the unfolded state is often found to be a rate-limiting step in protein folding. Moreover, kinases and phosphatases that act upon Ser/Thr-Pro motifs exhibit preferential recognition of either the cis- or trans-Pro conformer. Here, NMR spectra obtained at both atmospheric and high pressures indicate that the population of cis-Pro falls well below previous estimates, an effect attributed to the use of short peptides with charged termini in most prior model studies. For the intrinsically disordered protein α-synuclein, cis-Pro populations at all of its five X-Pro bonds are less than 5 %, with only modest ionic strength dependence and no detectable effect of the previously demonstrated interaction between the N- and C-terminal halves of the protein. Comparison to small peptides with the same amino-acid sequence indicates that peptides, particularly those with unblocked, oppositely charged amino and carboxyl end groups, strongly overestimate the amount of cis-Pro.

Sparse multidimensional iterative lineshape-enhanced (SMILE) reconstruction of both non-uniformly sampled and conventional NMR data.

Implementation of a new algorithm, SMILE, is described for reconstruction of non-uniformly sampled two-, three- and four-dimensional NMR data, which takes advantage of the known phases of the NMR spectrum and the exponential decay of underlying time domain signals. The method is very robust with respect to the chosen sampling protocol and, in its default mode, also extends the truncated time domain signals by a modest amount of non-sampled zeros. SMILE can likewise be used to extend conventional uniformly sampled data, as an effective multidimensional alternative to linear prediction. The program is provided as a plug-in to the widely used NMRPipe software suite, and can be used with default parameters for mainstream application, or with user control over the iterative process to possibly further improve reconstruction quality and to lower the demand on computational resources. For large data sets, the method is robust and demonstrated for sparsities down to ca 1%, and final all-real spectral sizes as large as 300 Gb. Comparison between fully sampled, conventionally processed spectra and randomly selected NUS subsets of this data shows that the reconstruction quality approaches the theoretical limit in terms of peak position fidelity and intensity. SMILE essentially removes the noise-like appearance associated with the point-spread function of signals that are a default of five-fold above the noise level, but impacts the actual thermal noise in the NMR spectra only minimally. Therefore, the appearance and interpretation of SMILE-reconstructed spectra is very similar to that of fully sampled spectra generated by Fourier transformation.

Atomic-resolution dynamics on the surface of amyloid-ß protofibrils probed by solution NMR.

Exchange dynamics between molecules free in solution and bound to the surface of a large supramolecular structure, a polymer, a membrane or solid support are important in many phenomena in biology and materials science. Here we present a novel and generally applicable solution NMR technique, known as dark-state exchange saturation transfer (DEST), to probe such exchange phenomena with atomic resolution. This is illustrated by the exchange reaction between amyloid-ß (Aß) monomers and polydisperse, NMR-invisible ('dark') protofibrils, a process of significant interest because the accumulation of toxic, aggregated forms of Aß, from small oligomers to very large assemblies, has been implicated in the aetiology of Alzheimer's disease. The (15)N-DEST experiment imprints with single-residue-resolution dynamic information on the protofibril-bound species in the form of (15)N transverse relaxation rates ((15)N-R(2)) and exchange kinetics between monomers and protofibrils onto the easily observed two-dimensional (1)H-(15)N correlation spectrum of the monomer. The exchanging species on the protofibril surface comprise an ensemble of sparsely populated states where each residue is either tethered to (through other residues) or in direct contact with the surface. The first eight residues exist predominantly in a mobile tethered state, whereas the largely hydrophobic central region and part of the carboxy (C)-terminal hydrophobic region are in direct contact with the protofibril surface for a significant proportion of the time. The C-terminal residues of both Aß40 and Aß42 display lower affinity for the protofibril surface, indicating that they are likely to be surface exposed rather than buried as in structures of Aß fibrils, and might therefore comprise the critical nucleus for fibril formation. The values, however, are significantly larger for the C-terminal residues of Aß42 than Aß40, which might explain the former's higher propensity for rapid aggregation and fibril formation.

Dissociation of the trimeric gp41 ectodomain at the lipid-water interface suggests an active role in HIV-1 Env-mediated membrane fusion.

The envelope glycoprotein gp41 mediates the process of membrane fusion that enables entry of the HIV-1 virus into the host cell. The actual fusion process involves a switch from a homotrimeric prehairpin intermediate conformation, consisting of parallel coiled-coil helices, to a postfusion state where the ectodomains are arranged as a trimer of helical hairpins, adopting a six-helix bundle (6HB) state. Here, we show by solution NMR spectroscopy that a water-soluble 6HB gp41 ectodomain binds to zwitterionic detergents that contain phosphocholine or phosphatidylcholine head groups and phospholipid vesicles that mimic T-cell membrane composition. Binding results in the dissociation of the 6HB and the formation of a monomeric state, where its two α-helices, N-terminal heptad repeat (NHR) and C-terminal heptad repeat (CHR), become embedded in the lipid-water interface of the virus and host cell. The atomic structure of the gp41 ectodomain monomer, based on NOE distance restraints and residual dipolar couplings, shows that the NHR and CHR helices remain mostly intact, but they completely lose interhelical contacts. The high affinity of the ectodomain helices for phospholipid surfaces suggests that unzippering of the prehairpin intermediate leads to a state where the NHR and CHR helices become embedded in the host cell and viral membranes, respectively, thereby providing a physical force for bringing these membranes into close juxtaposition before actual fusion.

Nuclear Magnetic Resonance Observation of α-Synuclein Membrane Interaction by Monitoring the Acetylation Reactivity of Its Lysine Side Chains.

The interaction between α-synuclein (αS) protein and lipid membranes is key to its role in synaptic vesicle homeostasis and plays a role in initiating fibril formation, which is implicated in Parkinson's disease. The natural state of αS inside the cell is generally believed to be intrinsically disordered, but chemical cross-linking experiments provided evidence of a tetrameric arrangement, which was reported to be rich in α-helical secondary structure based on circular dichroism (CD). Cross-linking relies on chemical modification of the protein's Lys C(ε) amino groups, commonly by glutaraldehyde, or by disuccinimidyl glutarate (DSG), with the latter agent preferred for cellular assays. We used ultra-high-resolution homonuclear decoupled nuclear magnetic resonance experiments to probe the reactivity of the 15 αS Lys residues toward N-succinimidyl acetate, effectively half the DSG cross-linker, which results in acetylation of Lys. The intensities of both side chain and backbone amide signals of acetylated Lys residues provide direct information about the reactivity, showing a difference of a factor of 2.5 between the most reactive (K6) and the least reactive (K102) residue. The presence of phospholipid vesicles decreases reactivity of most Lys residues by up to an order of magnitude at high lipid:protein stoichiometries (500:1), but only weakly at low ratios. The decrease in Lys reactivity is found to be impacted by lipid composition, even for vesicles that yield similar αS CD signatures. Our data provide new insight into the αS-bilayer interaction, including the pivotal state in which the available lipid surface is limited. Protection of Lys C(ε) amino groups by αS-bilayer interaction will strongly impact quantitative interpretation of DSG cross-linking experiments.

Monomeric Aß(1-40) and Aß(1-42) Peptides in Solution Adopt Very Similar Ramachandran Map Distributions That Closely Resemble Random Coil.

The pathogenesis of Alzheimer's disease is characterized by the aggregation and fibrillation of amyloid peptides Aß(1-40) and Aß(1-42) into amyloid plaques. Despite strong potential therapeutic interest, the structural pathways associated with the conversion of monomeric Aß peptides into oligomeric species remain largely unknown. In particular, the higher aggregation propensity and associated toxicity of Aß(1-42) compared to that of Aß(1-40) are poorly understood. To explore in detail the structural propensity of the monomeric Aß(1-40) and Aß(1-42) peptides in solution, we recorded a large set of nuclear magnetic resonance (NMR) parameters, including chemical shifts, nuclear Overhauser effects (NOEs), and J couplings. Systematic comparisons show that at neutral pH the Aß(1-40) and Aß(1-42) peptides populate almost indistinguishable coil-like conformations. Nuclear Overhauser effect spectra collected at very high resolution remove assignment ambiguities and show no long-range NOE contacts. Six sets of backbone J couplings ((3)JHNHα, (3)JC'C', (3)JC'Hα, (1)JHαCα, (2)JNCα, and (1)JNCα) recorded for Aß(1-40) were used as input for the recently developed MERA Ramachandran map analysis, yielding residue-specific backbone ϕ/ψ torsion angle distributions that closely resemble random coil distributions, the absence of a significantly elevated propensity for ß-conformations in the C-terminal region of the peptide, and a small but distinct propensity for αL at K28. Our results suggest that the self-association of Aß peptides into toxic oligomers is not driven by elevated propensities of the monomeric species to adopt ß-strand-like conformations. Instead, the accelerated disappearance of Aß NMR signals in D2O over H2O, particularly pronounced for Aß(1-42), suggests that intermolecular interactions between the hydrophobic regions of the peptide dominate the aggregation process.

A Novel MHC-I Surface Targeted for Binding by the MCMV m06 Immunoevasin Revealed by Solution NMR.

As part of its strategy to evade detection by the host immune system, murine cytomegalovirus (MCMV) encodes three proteins that modulate cell surface expression of major histocompatibility complex class I (MHC-I) molecules: the MHC-I homolog m152/gp40 as well as the m02-m16 family members m04/gp34 and m06/gp48. Previous studies of the m04 protein revealed a divergent Ig-like fold that is unique to immunoevasins of the m02-m16 family. Here, we engineer and characterize recombinant m06 and investigate its interactions with full-length and truncated forms of the MHC-I molecule H2-L(d) by several techniques. Furthermore, we employ solution NMR to map the interaction footprint of the m06 protein on MHC-I, taking advantage of a truncated H2-L(d), "mini-H2-L(d)," consisting of only the α1α2 platform domain. Mini-H2-L(d) refolded in vitro with a high affinity peptide yields a molecule that shows outstanding NMR spectral features, permitting complete backbone assignments. These NMR-based studies reveal that m06 binds tightly to a discrete site located under the peptide-binding platform that partially overlaps with the ß2-microglobulin interface on the MHC-I heavy chain, consistent with in vitro binding experiments showing significantly reduced complex formation between m06 and ß2-microglobulin-associated MHC-I. Moreover, we carry out NMR relaxation experiments to characterize the picosecond-nanosecond dynamics of the free mini-H2-L(d) MHC-I molecule, revealing that the site of interaction is highly ordered. This study provides insight into the mechanism of the interaction of m06 with MHC-I, suggesting a structural manipulation of the target MHC-I molecule at an early stage of the peptide-loading pathway.

Accurate measurement of (3)J(HNHα) couplings in small or disordered proteins from WATERGATE-optimized TROSY spectra.

Provided that care is taken in adjusting the WATERGATE element of a (1)H-(15)N TROSY-HSQC experiment, such that neither the water magnetization nor the (1)H(α) protons are inverted by its final 180° pulse, (3)JHNHα couplings can be measured directly from splittings in the (1)H dimension of the spectrum. With band-selective (1)H decoupling, very high (15)N resolution can be achieved. A complete set of (3)JHNHα values, ranging from 3.4 to 10.1 Hz was measured for the 56-residue third domain of IgG-binding protein G (GB3). Using the H-N-C(α)-H(α) dihedral angles extracted from a RDC-refined structure of GB3, (3)JHNHα values predicted by a previously parameterized Karplus equation agree to within a root-mean-square deviation (rmsd) of 0.37 Hz with the experimental data. Values measured for the Alzheimer's implicated Aß(1-40) peptide fit to within an rmsd of 0.45 Hz to random coil (3)JHNHα values.

Quantitative evaluation of positive ϕ angle propensity in flexible regions of proteins from three-bond J couplings.

(3)JHNHα and (3)JC'C' couplings can be readily measured in isotopically enriched proteins and were shown to contain precise information on the backbone torsion angles, ϕ, sampled in disordered regions of proteins. However, quantitative interpretation of these couplings required the population of conformers with positive ϕ angles to be very small. Here, we demonstrate that this restriction can be removed by measurement of (3)JC'Hα values. Even though the functional forms of the (3)JC'Hα and (3)JHNHα Karplus equations are the same, large differences in their coefficients enable accurate determination of the fraction of time that positive ϕ angles are sampled. A four-dimensional triple resonance HACANH[C'] E.COSY experiment is introduced to simultaneously measure (3)JC'Hα and (3)JHNC' in the typically very congested spectra of disordered proteins. High resolution in these spectra is obtained by non-uniform sampling (in the 0.1-0.5% range). Application to the intrinsically disordered protein α-synuclein shows that while most residues have close-to-zero positive ϕ angle populations, up to 16% positive ϕ population is observed for Asn residues. Positive ϕ angle populations determined with the new approach agree closely with consensus values from protein coil libraries and prior analysis of a large set of other NMR parameters. The combination of (3)JHNC' and (3)JC'C' provides information about the amplitude of ϕ angle dynamics.

Probing the transient dark state of substrate binding to GroEL by relaxation-based solution NMR.

The mechanism whereby the prototypical chaperonin GroEL performs work on substrate proteins has not yet been fully elucidated, hindered by lack of detailed structural and dynamic information on the bound substrate. Previous investigations have produced conflicting reports on the state of GroEL-bound polypeptides, largely due to the transient and dynamic nature of these complexes. Here, we present a unique approach, based on combined analysis of four complementary relaxation-based NMR experiments, to probe directly the "dark" NMR-invisible state of the model, intrinsically disordered, polypeptide amyloid ß (Aß40) bound to GroEL. The four NMR experiments, lifetime line-broadening, dark-state exchange saturation transfer, relaxation dispersion, and small exchange-induced chemical shifts, are dependent in different ways on the overall exchange rates and populations of the free and bound states of the substrate, as well as on residue-specific dynamics and structure within the bound state as reported by transverse magnetization relaxation rates and backbone chemical shifts, respectively. Global fitting of all the NMR data shows that the complex is transient with a lifetime of <1 ms, that binding involves two predominantly hydrophobic segments corresponding to predicted GroEL consensus binding sequences, and that the structure of the bound polypeptide remains intrinsically and dynamically disordered with minimal changes in secondary structure propensity relative to the free state. Our results establish a unique method to observe NMR-invisible dynamic states of GroEL-bound substrates and to describe at atomic resolution the events between substrate binding and encapsulation that are crucial for understanding the normal and stress-related metabolic function of chaperonins.

Side Chain Conformational Distributions of a Small Protein Derived from Model-Free Analysis of a Large Set of Residual Dipolar Couplings.

Accurate quantitative measurement of structural dispersion in proteins remains a prime challenge to both X-ray crystallography and NMR spectroscopy. Here we use a model-free approach based on measurement of many residual dipolar couplings (RDCs) in differentially orienting aqueous liquid crystalline solutions to obtain the side chain χ1 distribution sampled by each residue in solution. Applied to the small well-ordered model protein GB3, our approach reveals that the RDC data are compatible with a single narrow distribution of side chain χ1 angles for only about 40% of the residues. For more than half of the residues, populations greater than 10% for a second rotamer are observed, and four residues require sampling of three rotameric states to fit the RDC data. In virtually all cases, sampled χ1 values are found to center closely around ideal g(-), g(+) and t rotameric angles, even though no rotamer restraint is used when deriving the sampled angles. The root-mean-square difference between experimental (3)JHαHß couplings and those predicted by the Haasnoot-parametrized, motion-adjusted Karplus equation reduces from 2.05 to 0.75 Hz when using the new rotamer analysis instead of the 1.1-Å X-ray structure as input for the dihedral angles. A comparison between observed and predicted (3)JHαHß values suggests that the root-mean-square amplitude of χ1 angle fluctuations within a given rotamer well is ca. 20°. The quantitatively defined side chain rotamer equilibria obtained from our study set new benchmarks for evaluating improved molecular dynamics force fields, and also will enable further development of quantitative relations between side chain chemical shift and structure.

High accuracy of Karplus equations for relating three-bond J couplings to protein backbone torsion angles.

(3) JC'C' and (3) JHNHα couplings are related to the intervening backbone torsion angle ${\varphi }$ by standard Karplus equations. Although these couplings are known to be affected by parameters other than ${\varphi }$, including H-bonding, valence angles and residue type, experimental results and quantum calculations indicate that the impact of these latter parameters is typically very small. The solution NMR structure of protein GB3, newly refined by using extensive sets of residual dipolar couplings, yields 50-60 % better Karplus equation agreement between ${\varphi }$ angles and experimental (3) JC'C' and (3) JHNHα values than does the high-resolution X-ray structure. In intrinsically disordered proteins, (3) JC'C' and (3) JHNHα couplings can be measured at even higher accuracy, and the impact of factors other than the intervening torsion angle on (3) J will be smaller than in folded proteins, making these couplings exceptionally valuable reporters on the ensemble of ${\varphi }$ angles sampled by each residue.