Solution structure of the 30 kDa N-terminal domain of enzyme I of the Escherichia coli phosphoenolpyruvate:sugar phosphotransferase system by multidimensional NMR

Enabling site-specific NMR investigations of therapeutic Fab using a cell-free based isotopic labeling approach: application to anti-LAMP1 Fab

Monoclonal antibodies (mAbs) are biotherapeutics that have achieved outstanding success in treating many life-threatening and chronic diseases. The recognition of an antigen is mediated by the fragment antigen binding (Fab) regions composed by four different disulfide bridge-linked immunoglobulin domains. NMR is a powerful method to assess the integrity, the structure and interaction of Fabs, but site specific analysis has been so far hampered by the size of the Fabs and the lack of approaches to produce isotopically labeled samples. We proposed here an efficient in vitro method to produce [15N, 13C, 2H]-labeled Fabs enabling high resolution NMR investigations of these powerful therapeutics. As an open system, the cell-free expression mode enables fine-tuned control of the redox potential in presence of disulfide bond isomerase to enhance the formation of native disulfide bonds. Moreover, inhibition of transaminases in the S30 cell-free extract offers the opportunity to produce perdeuterated Fab samples directly in 1H2O medium, without the need for a time-consuming and inefficient refolding process. This specific protocol was applied to produce an optimally labeled sample of a therapeutic Fab, enabling the sequential assignment of 1HN, 15N, 13C', 13Cα, 13Cβ resonances of a full-length Fab. 90% of the backbone resonances of a Fab domain directed against the human LAMP1 glycoprotein were assigned successfully, opening new opportunities to study, at atomic resolution, Fabs' higher order structures, dynamics and interactions, using solution-state NMR.

Solution Structure Ensembles of the Open and Closed Forms of the ∼130 kDa Enzyme I via AlphaFold Modeling, Coarse Grained Simulations, and NMR

Large-scale interdomain rearrangements are essential to protein function, governing the activity of large enzymes and molecular machineries. Yet, obtaining an atomic-resolution understanding of how the relative domain positioning is affected by external stimuli is a hard task in modern structural biology. Here, we show that combining structural modeling by AlphaFold2 with coarse-grained molecular dynamics simulations and NMR residual dipolar coupling data is sufficient to characterize the spatial domain organization of bacterial enzyme I (EI), a ∼130 kDa multidomain oligomeric protein that undergoes large-scale conformational changes during its catalytic cycle. In particular, we solve conformational ensembles for EI at two different experimental temperatures and demonstrate that a lower temperature favors sampling of the catalytically competent closed state of the enzyme. These results suggest a role for conformational entropy in the activation of EI and demonstrate the ability of our protocol to detect and characterize the effect of external stimuli (such as mutations, ligand binding, and post-translational modifications) on the interdomain organization of multidomain proteins. We expect the ensemble refinement protocol described here to be easily transferrable to the investigation of the structure and dynamics of other uncharted multidomain systems and have assembled a Google Colab page (https://potoyangroup.github.io/Seq2Ensemble/) to facilitate implementation of the presented methodology elsewhere.

Conformation and Affinity Modulations by Multiple Phosphorylation Occurring in the BIN1 SH3 Domain Binding Site of the Tau Protein Proline-Rich Region

An increase in phosphorylation of the Tau protein is associated with Alzheimer's disease (AD) progression through unclear molecular mechanisms. In general, phosphorylation modifies the interaction of intrinsically disordered proteins, such as Tau, with other proteins; however, elucidating the structural basis of this regulation mechanism remains challenging. The bridging integrator-1 gene is an AD genetic determinant whose gene product, BIN1, directly interacts with Tau. The proline-rich motif recognized within a Tau(210-240) peptide by the SH3 domain of BIN1 (BIN1 SH3) is defined as ₂₁₆PTPP₂₁₉, and this interaction is modulated by phosphorylation. Phosphorylation of T217 within the Tau(210-240) peptide led to a 6-fold reduction in the affinity, while single phosphorylation at either T212, T231, or S235 had no effect on the interaction. Nonetheless, combined phosphorylation of T231 and S235 led to a 3-fold reduction in the affinity, although these phosphorylations are not within the BIN1 SH3-bound region of the Tau peptide. Using nuclear magnetic resonance (NMR) spectroscopy, these phosphorylations were shown to affect the local secondary structure and dynamics of the Tau(210-240) peptide. Models of the (un)phosphorylated peptides were obtained from molecular dynamics (MD) simulation validated by experimental data and showed compaction of the phosphorylated peptide due to increased salt bridge formation. This dynamic folding might indirectly impact the BIN1 SH3 binding by a decreased accessibility of the binding site. Regulation of the binding might thus not only be due to local electrostatic or steric effects from phosphorylation but also to the modification of the conformational properties of Tau.

Progress toward automated methyl assignments for methyl-TROSY applications

Methyl-TROSY spectroscopy has extended the reach of solution-state NMR to supra-molecular machineries over 100 kDa in size. Methyl groups are ideal probes for studying structure, dynamics, and protein-protein interactions in quasi-physiological conditions with atomic resolution. Successful implementation of the methodology requires accurate methyl chemical shift assignment, and the task still poses a significant challenge in the field. In this work, we outline the current state of technology for methyl labeling, data collection, data analysis, and nuclear Overhauser effect (NOE)-based automated methyl assignment approaches. We present MAGIC-Act and MAGIC-View, two Python extensions developed as part of the popular NMRFAM-Sparky package, and MAGIC-Net a standalone structure-based network analysis program. MAGIC-Act conducts statistically driven amino acid typing, Leu/Val pairing guided by 3D HMBC-HMQC, and NOESY cross-peak symmetry checking. MAGIC-Net provides model-based NOE statistics to aid in selection of a methyl labeling scheme. The programs provide a versatile, semi-automated framework for rapid methyl assignment.

Automatic structure-based NMR methyl resonance assignment in large proteins

Isotopically labeled methyl groups provide NMR probes in large, otherwise deuterated proteins. However, the resonance assignment constitutes a bottleneck for broader applicability of methyl-based NMR. Here, we present the automated MethylFLYA method for the assignment of methyl groups that is based on methyl-methyl nuclear Overhauser effect spectroscopy (NOESY) peak lists. MethylFLYA is applied to five proteins (28–358 kDa) comprising a total of 708 isotope-labeled methyl groups, of which 612 contribute NOESY cross peaks. MethylFLYA confidently assigns 488 methyl groups, i.e. 80% of those with NOESY data. Of these, 459 agree with the reference, 6 were different, and 23 were without reference assignment. MethylFLYA assigns significantly more methyl groups than alternative algorithms, has an average error rate of 1%, modest runtimes of 0.4–1.2 h, and can handle arbitrary isotope labeling patterns and data from other types of NMR spectra.

The structures and dynamics of large proteins can be studied with methyl-based NMR but peak assignment is still challenging. Here the authors present MethylFLYA that allows automated assignment of methyl groups and apply it to five proteins with molecular weights in the range from 28 to 358 kDa.

Impact of lipid binding on the tertiary structure and allergenic potential of Jug r 3, the non-specific lipid transfer protein from walnut

Plant non-specific lipid transfer proteins type 1 (nsLTP1) are small basic proteins with a hydrophobic cavity able to host a number of different ligands: i.e. fatty acids, fatty acyl-CoA, phospholipids, glycolipids, and hydroxylated fatty acids. However, ligand binding specificity differs among nsLTPs. Within this protein family, Jug r 3 from walnut has been identified as a major allergen. So far, data on the structural characterization of Jug r 3 and its lipid binding capacity are lacking. We report the results from a fluorescence-based ligand-binding assay and ligand-based NMR experiments, to study the binding interactions between Jug r 3 and the 18-carbon monounsaturated oleic acid. Furthermore, protein-based NMR experiments were employed to detect the oleate binding site of Jug r 3. The NMR data were used to dock the oleate molecule into the structural model of Jug r 3. Finally, the impact of the interaction on the allergenic potential of Jug r 3 was investigated by IgE ELISA with 6 sera from walnut allergic patients. Our data corroborate the hypothesis of direct impact of food-derived matrix on the IgE reactivity of nsLTPs.

Selective Methyl Labeling of Proteins: Enabling Structural and Mechanistic Studies As Well As Drug Discovery Applications by Solution-State NMR

Escherichia coli expression protocols for selective labeling of methyl groups in proteins have been essential in expanding the size range of targets that can be studied by biomolecular NMR. Based on the initial work achieving selective labeling of isoleucine, leucine, and valine residues, additional methods were developed over the past years which enabled the individual and/or simultaneous combinatorial labeling of all methyl containing amino acids. Together with the introduction of new methyl-optimized NMR experiments, this now allows the detailed characterization of protein-ligand interactions as well as mechanistic and dynamic processes of protein-protein complexes up to 1MDa in size. In this chapter, we provide a general introduction to selective labeling of proteins using E. coli-based expression systems, describe the considerations taken into account prior to the selective labeling of a protein, and include the protocols used to produce such proteins. An overview of applications using selectively labeled proteins with an emphasis on examples relevant to the drug discovery process is then presented.

Structural Basis of Tau Interaction With BIN1 and Regulation by Tau Phosphorylation

Bridging integrator-1 (BIN1) gene is associated with an increased risk to develop Alzheimer's disease, a tauopathy characterized by intra-neuronal accumulation of phosphorylated Tau protein as paired helical filaments. Direct interaction of BIN1 and Tau proteins was demonstrated to be mediated through BIN1 SH3 C-terminal domain and Tau (210-240) peptide within Tau proline-rich domain. We previously showed that BIN1 SH3 interaction with Tau is decreased by phosphorylation within Tau proline-rich domain, of at least T231. In addition, the BIN1/Tau interaction is characterized by a dynamic equilibrium between a closed and open conformations of BIN1 isoform 1, involving an intramolecular interaction with its C-terminal BIN1 SH3 domain. However, the role of the BIN1/Tau interaction, and its potential dysregulation in Alzheimer's disease, is not yet fully understood. Here we showed that within Tau (210-240) peptide, among the two proline-rich motifs potentially recognized by SH3 domains, only motif P²¹⁶TPPTR²²¹ is bound by BIN1 SH3. A structural model of the complex between BIN1 SH3 and Tau peptide (213-229), based on nuclear magnetic resonance spectroscopy data, revealed the molecular detail of the interaction. P216 and P219 within the proline-rich motif were in direct contact with the aromatic F588 and W562 of the BIN1 SH3 domain. The contact surface is extended through electrostatic interactions between the positively charged R221 and K224 residues of Tau peptide and those negatively charged of BIN1 SH3, corresponding to E556 and E557. We next investigated the impact of multiple Tau phosphorylations within Tau (210-240) on its interaction with BIN1 isoform 1. Tau (210-240) phosphorylated at four different sites (T212, T217, T231, and S235), contrary to unphosphorylated Tau, was unable to compete with the intramolecular interaction of BIN1 SH3 domain with its CLAP domain. In accordance, the affinity of BIN1 SH3 for phosphorylated Tau (210-240) peptide was reduced, with a five-fold increase in the dissociation constant, from a Kd of 44 to 256 μM. This study highlights the complexity of the regulation of BIN1 isoform 1 with Tau. As abnormal phosphorylation of Tau is linked to the pathology development, this regulation by phosphorylation might have important functional consequences.

Active Site Gate Dynamics Modulate the Catalytic Activity of the Ubiquitination Enzyme E2-25K

The ubiquitin proteasome system (UPS) signals for degradation of proteins through attachment of K48-linked polyubiquitin chains, or alterations in protein-protein recognition through attachment of K63-linked chains. Target proteins are ubiquitinated in three sequential chemical steps by a three-component enzyme system. Ubiquitination, or E2 enzymes, catalyze the central step by facilitating reaction of a target protein lysine with the C-terminus of Ub that is attached to the active site cysteine of the E2 through a thioester bond. E2 reactivity is modulated by dynamics of an active site gate, whose central residue packs against the active site cysteine in a closed conformation. Interestingly, for the E2 Ubc13, which specifically catalyzes K63-linked ubiquitination, the central gate residue adopts an open conformation. We set out to determine if active site gate dynamics play a role in catalysis for E2-25K, which adopts the canonical, closed gate conformation, and which selectively synthesizes K48-linked ubiquitin chains. Gate dynamics were characterized using mutagenesis of key residues, combined with enzyme kinetics measurements, and main chain NMR relaxation. The experimental data were interpreted with all atom MD simulations. The data indicate that active site gate opening and closing rates for E2-25K are precisely balanced.

High quality NMR structures: a new force field with implicit water and membrane solvation for Xplor-NIH

Structure determination of proteins by NMR is unique in its ability to measure restraints, very accurately, in environments and under conditions that closely mimic those encountered in vivo. For example, advances in solid-state NMR methods enable structure determination of membrane proteins in detergent-free lipid bilayers, and of large soluble proteins prepared by sedimentation, while parallel advances in solution NMR methods and optimization of detergent-free lipid nanodiscs are rapidly pushing the envelope of the size limit for both soluble and membrane proteins. These experimental advantages, however, are partially squandered during structure calculation, because the commonly used force fields are purely repulsive and neglect solvation, Van der Waals forces and electrostatic energy. Here we describe a new force field, and updated energy functions, for protein structure calculations with EEFx implicit solvation, electrostatics, and Van der Waals Lennard-Jones forces, in the widely used program Xplor-NIH. The new force field is based primarily on CHARMM22, facilitating calculations with a wider range of biomolecules. The new EEFx energy function has been rewritten to enable OpenMP parallelism, and optimized to enhance computation efficiency. It implements solvation, electrostatics, and Van der Waals energy terms together, thus ensuring more consistent and efficient computation of the complete nonbonded energy lists. Updates in the related python module allow detailed analysis of the interaction energies and associated parameters. The new force field and energy function work with both soluble proteins and membrane proteins, including those with cofactors or engineered tags, and are very effective in situations where there are sparse experimental restraints. Results obtained for NMR-restrained calculations with a set of five soluble proteins and five membrane proteins show that structures calculated with EEFx have significant improvements in accuracy, precision, and conformation, and that structure refinement can be obtained by short relaxation with EEFx to obtain improvements in these key metrics. These developments broaden the range of biomolecular structures that can be calculated with high fidelity from NMR restraints.

Structure of the NPr:EINNtr Complex: Mechanism for Specificity in Paralogous Phosphotransferase Systems

Paralogous enzymes arise from gene duplication events that confer a novel function, although it is unclear how cross-reaction between the original and duplicate protein interaction network is minimized. We investigated HPr:EIsugar and NPr:EINtr, the initial complexes of paralogous phosphorylation cascades involved in sugar import and nitrogen regulation in bacteria, respectively. Although the HPr:EIsugar interaction has been well characterized, involving multiple complexes and transient interactions, the exact nature of the NPr:EINtr complex was unknown. We set out to identify the key features of the interaction by performing binding assays and elucidating the structure of NPr in complex with the phosphorylation domain of EINtr (EINNtr), using a hybrid approach involving X-ray, homology, and sparse nuclear magnetic resonance. We found that the overall fold and active-site structure of the two complexes are conserved in order to maintain productive phosphorylation, however, the interface surface potential differs between the two complexes, which prevents cross-reaction.

Automated assignment of NMR chemical shifts based on a known structure and 4D spectra

Apart from their central role during 3D structure determination of proteins the backbone chemical shift assignment is the basis for a number of applications, like chemical shift perturbation mapping and studies on the dynamics of proteins. This assignment is not a trivial task even if a 3D protein structure is known and needs almost as much effort as the assignment for structure prediction if performed manually. We present here a new algorithm based solely on 4D [(1)H,(15)N]-HSQC-NOESY-[(1)H,(15)N]-HSQC spectra which is able to assign a large percentage of chemical shifts (73-82 %) unambiguously, demonstrated with proteins up to a size of 250 residues. For the remaining residues, a small number of possible assignments is filtered out. This is done by comparing distances in the 3D structure to restraints obtained from the peak volumes in the 4D spectrum. Using dead-end elimination, assignments are removed in which at least one of the restraints is violated. Including additional information from chemical shift predictions, a complete unambiguous assignment was obtained for Ubiquitin and 95 % of the residues were correctly assigned in the 251 residue-long N-terminal domain of enzyme I. The program including source code is available at https://github.com/thomasexner/4Dassign .

Hybrid Approaches to Structural Characterization of Conformational Ensembles of Complex Macromolecular Systems Combining NMR Residual Dipolar Couplings and Solution X-ray Scattering

Solving structures or structural ensembles of large macromolecular systems in solution poses a challenging problem. While NMR provides structural information at atomic resolution, increased spectral complexity, chemical shift overlap, and short transverse relaxation times (associated with slow tumbling) render application of the usual techniques that have been so successful for medium sized systems (<50 kDa) difficult. Solution X-ray scattering, on the other hand, is not limited by molecular weight but only provides low resolution structural information related to the overall shape and size of the system under investigation. Here we review how combining atomic resolution structures of smaller domains with sparse experimental data afforded by NMR residual dipolar couplings (which yield both orientational and shape information) and solution X-ray scattering data in rigid-body simulated annealing calculations provides a powerful approach for investigating the structural aspects of conformational dynamics in large multidomain proteins. The application of this hybrid methodology is illustrated for the 128 kDa dimer of bacterial Enzyme I which exists in a variety of open and closed states that are sampled at various points in the catalytic cycles, and for the capsid protein of the human immunodeficiency virus.

Dynamic equilibrium between closed and partially closed states of the bacterial Enzyme I unveiled by solution NMR and X-ray scattering

Enzyme I (EI) is the first component in the bacterial phosphotransferase system, a signal transduction pathway in which phosphoryl transfer through a series of bimolecular protein-protein interactions is coupled to sugar transport across the membrane. EI is a multidomain, 128-kDa homodimer that has been shown to exist in two conformational states related to one another by two large (50-90°) rigid body domain reorientations. The open conformation of apo EI allows phosphoryl transfer from His189 located in the N-terminal domain α/β (EIN(α/β)) subdomain to the downstream protein partner bound to the EIN(α) subdomain. The closed conformation, observed in a trapped phosphoryl transfer intermediate, brings the EIN(α/β) subdomain into close proximity to the C-terminal dimerization domain (EIC), thereby permitting in-line phosphoryl transfer from phosphoenolpyruvate (PEP) bound to EIC to His189. Here, we investigate the solution conformation of a complex of an active site mutant of EI (H189A) with PEP. Simulated annealing refinement driven simultaneously by solution small angle X-ray scattering and NMR residual dipolar coupling data demonstrates unambiguously that the EI(H189A)-PEP complex exists in a dynamic equilibrium between two approximately equally populated conformational states, one corresponding to the closed structure and the other to a partially closed species. The latter likely represents an intermediate in the open-to-closed transition.

Large interdomain rearrangement triggered by suppression of micro- to millisecond dynamics in bacterial Enzyme I

Enzyme I (EI), the first component of the bacterial phosphotransfer signal transduction system, undergoes one of the largest substrate-induced interdomain rearrangements documented to date. Here, we characterize the perturbations generated by two small molecules, the natural substrate phosphoenolpyruvate (PEP) and the inhibitor α-ketoglutarate (αKG), on the structure and dynamics of EI using NMR, small-angle X-ray scattering (SAXS) and biochemical techniques. The results indicate unambiguously that the open-to-closed conformational switch of EI is triggered by complete suppression of micro- to millisecond dynamics within the C-terminal domain of EI. Indeed, we show that a ligand-induced transition from a dynamic to a more rigid conformational state of the C-terminal domain stabilizes the interface between the N- and C-terminal domains observed in the structure of the closed state, thereby promoting the resulting conformational switch and autophosphorylation of EI. The mechanisms described here may be common to several other multidomain proteins and allosteric systems.

Using small angle solution scattering data in Xplor-NIH structure calculations

This contribution describes the use of small and wide angle X-ray and small angle neutron scattering for biomolecular structure calculation using the program Xplor-NIH, both with and without NMR data. The current algorithms used for calculating scattering curves are described, and the use of scattering data as a structural restraint is given concrete form as a fragment of an Xplor-NIH structure calculation script. We review five examples of the use of scattering data in structure calculation, including the treatment of single domain proteins, nucleic acids, structure determination of large proteins, and the use of ensemble representations to characterize small and large amplitude motions.

A practical implicit solvent potential for NMR structure calculation

The benefits of protein structure refinement in water are well documented. However, performing structure refinement with explicit atomic representation of the solvent molecules is computationally expensive and impractical for NMR-restrained structure calculations that start from completely extended polypeptide templates. Here we describe a new implicit solvation potential, EEFx (Effective Energy Function for XPLOR-NIH), for NMR-restrained structure calculations of proteins in XPLOR-NIH. The key components of EEFx are an energy term for solvation energy that works together with other nonbonded energy functions, and a dedicated force field for conformational and nonbonded protein interaction parameters. The initial results obtained with EEFx show that significant improvements in structural quality can be obtained. EEFx is computationally efficient and can be used both to fold and refine structures. Overall, EEFx improves the quality of protein conformation and nonbonded atomic interactions. Moreover, such benefits are accompanied by enhanced structural precision and enhanced structural accuracy, reflected in improved agreement with the cross-validated dipolar coupling data. Finally, implementation of EEFx calculations is straightforward and computationally efficient. Overall, EEFx provides a useful method for the practical calculation of experimental protein structures in a physically realistic environment.

A NMR experiment for simultaneous correlations of valine and leucine/isoleucine methyls with carbonyl chemical shifts in proteins

A methyl-detected 'out-and-back' NMR experiment for obtaining simultaneous correlations of methyl resonances of valine and isoleucine/leucine residues with backbone carbonyl chemical shifts, SIM-HMCM(CGCBCA)CO, is described. The developed pulse-scheme serves the purpose of convenience in recording a single data set for all Ile(δ1), Leu(δ) and Val(γ) (ILV) methyl positions instead of acquiring two separate spectra selective for valine or leucine/isoleucine residues. The SIM-HMCM(CGCBCA)CO experiment can be used for ILV methyl assignments in moderately sized protein systems (up to ~100 kDa) where the backbone chemical shifts of (13)C(α), (13)Cβ and (13)CO are known from prior NMR studies and where some losses in sensitivity can be tolerated for the sake of an overall reduction in NMR acquisition time.

Four-bond deuterium isotope effects on the chemical shifts of amide nitrogens in proteins

An approach towards precision NMR measurements of four-bond deuterium isotope effects on the chemical shifts of backbone amide nitrogen nuclei in proteins is described. Three types of four-bond (15) N deuterium isotope effects are distinguished depending on the site of proton-to-deuterium substitution: (4)ΔN(N(i-1)D), (4)ΔN(N(i+1)D) and (4)ΔN(Cβ,(i-1)D). All the three types of isotope shifts are quantified in the (partially) deuterated protein ubiquitin. The (4)ΔN(N(i+1)D) and (4)ΔN(C(β,i-1)D) effects are by far the largest in magnitude and vary between 16 and 75 ppb and -18 and 46 ppb, respectively. A semi-quantitative correlation between experimental (4)ΔN(N(i+1)D) and (4)ΔN(C(β,i-1)D) values and the distances between nitrogen nuclei and the sites of (1)H-to-D substitution is noted. The largest isotope shifts in both cases correspond to the shortest inter-nuclear distances.

Structure, dynamics and biophysics of the cytoplasmic protein-protein complexes of the bacterial phosphoenolpyruvate: sugar phosphotransferase system

The bacterial phosphotransferase system (PTS) couples phosphoryl transfer, via a series of bimolecular protein-protein interactions, to sugar transport across the membrane. The multitude of complexes in the PTS provides a paradigm for studying protein interactions, and for understanding how the same binding surface can specifically recognize a diverse array of targets. Fifteen years of work aimed at solving the solution structures of all soluble protein-protein complexes of the PTS has served as a test bed for developing NMR and integrated hybrid approaches to study larger complexes in solution and to probe transient, spectroscopically invisible states, including encounter complexes. We review these approaches, highlighting the problems that can be tackled with these methods, and summarize the current findings on protein interactions.

Structural basis for enzyme I inhibition by α-ketoglutarate

Creating new bacterial strains in which carbon and nitrogen metabolism are uncoupled is potentially very useful for optimizing yields of microbial produced chemicals from renewable carbon sources. However, the mechanisms that balance carbon and nitrogen consumption in bacteria are poorly understood. Recently, α-ketoglutarate (αKG), the carbon substrate for ammonia assimilation, has been observed to inhibit Escherichia coli enzyme I (EI), the first component of the bacterial phosphotransferase system (PTS), thereby providing a direct biochemical link between central carbon and nitrogen metabolism. Here we investigate the EI-αKG interaction by NMR and enzymatic assays. We show that αKG binds with a KD of ∼2.2 mM at the active site of EI, acting as a competitive inhibitor. In addition, we use molecular docking simulations to derive a structural model of the enzyme-inhibitor complex that is fully consistent with NMR and analytical ultracentrifugation data. We expect that the EI-αKG structure presented here will provide a starting point for structure-based design of EI mutants resistant to αKG.

Smooth statistical torsion angle potential derived from a large conformational database via adaptive kernel density estimation improves the quality of NMR protein structures

Statistical potentials that embody torsion angle probability densities in databases of high-quality X-ray protein structures supplement the incomplete structural information of experimental nuclear magnetic resonance (NMR) datasets. By biasing the conformational search during the course of structure calculation toward highly populated regions in the database, the resulting protein structures display better validation criteria and accuracy. Here, a new statistical torsion angle potential is developed using adaptive kernel density estimation to extract probability densities from a large database of more than 10⁶ quality-filtered amino acid residues. Incorporated into the Xplor-NIH software package, the new implementation clearly outperforms an older potential, widely used in NMR structure elucidation, in that it exhibits simultaneously smoother and sharper energy surfaces, and results in protein structures with improved conformation, nonbonded atomic interactions, and accuracy.

Stability and binding of the phosphorylated species of the N-terminal domain of enzyme I and the histidine phosphocarrier protein from the Streptomyces coelicolor phosphoenolpyruvate:sugar phosphotransferase system

The phosphotransferase system (PTS) is involved in the use of carbon sources in bacteria. It is formed by two general proteins: enzyme I (EI) and the histidine phosphocarrier (HPr), and various sugar-specific permeases. EI is formed by two domains, with the N-terminal domain (EIN) being responsible for the binding to HPr. In low-G+C Gram-positive bacteria, HPr becomes phosphorylated not only by phosphoenolpyruvate (PEP) at the active-site histidine, but also by ATP at a serine. In this work, we have characterized: (i) the stability and binding affinities between the active-site-histidine phosphorylated species of HPr and the EIN from Streptomyces coelicolor; and (ii) the stability and binding affinities of the species involving the phosphorylation at the regulatory serine of HPr(sc). Our results show that the phosphorylated active-site species of both proteins are less stable than the unphosphorylated counterparts. Conversely, the Hpr-S47D, which mimics phosphorylation at the regulatory serine, is more stable than wild-type HPr(sc) due to helical N-capping effects, as suggested by the modeled structure of the protein. Binding among the phosphorylated and unphosphorylated species is always entropically driven, but the affinity and the enthalpy vary widely.

Deuterium isotope shifts for backbone ¹H, ¹⁵N and ¹³C nuclei in intrinsically disordered protein α-synuclein

Intrinsically disordered proteins (IDPs) are abundant in nature and characterization of their potential structural propensities remains a widely pursued but challenging task. Analysis of NMR secondary chemical shifts plays an important role in such studies, but the output of such analyses depends on the accuracy of reference random coil chemical shifts. Although uniform perdeuteration of IDPs can dramatically increase spectral resolution, a feature particularly important for the poorly dispersed IDP spectra, the impact of deuterium isotope shifts on random coil values has not yet been fully characterized. Very precise (2)H isotope shift measurements for (13)C(α), (13)C(β), (13)C', (15)N, and (1)H(N) have been obtained by using a mixed sample of protonated and uniformly perdeuterated α-synuclein, a protein with chemical shifts exceptionally close to random coil values. Decomposition of these isotope shifts into one-bond, two-bond and three-bond effects as well as intra- and sequential residue contributions shows that such an analysis, which ignores conformational dependence, is meaningful but does not fully describe the total isotope shift to within the precision of the measurements. Random coil (2)H isotope shifts provide an important starting point for analysis of such shifts in structural terms in folded proteins, where they are known to depend strongly on local geometry.

Calorimetric and spectroscopic investigation of the interaction between the C-terminal domain of Enzyme I and its ligands

Enzyme I initiates a series of phosphotransfer reactions during sugar uptake in the bacterial phosphotransferase system. Here, we have isolated a stable recombinant C-terminal domain of Enzyme I (EIC) of Escherichia coli and characterized its interaction with the N-terminal domain of Enzyme I (EIN) and also with various ligands. EIC can phosphorylate EIN, but their binding is transient regardless of the presence of phosphoenolpyruvate (PEP). Circular dichroism and NMR indicate that ligand binding to EIC induces changes near aromatic groups but not in the secondary structure of EIC. Binding of PEP to EIC is an endothermic reaction with the equilibrium dissociation constant (K(D) ) of 0.28 mM, whereas binding of the inhibitor oxalate is an exothermic reaction with K(D) of 0.66 mM from calorimetry. The binding thermodynamics of EIC and PEP compared to that of Enzyme I (EI) and PEP reveals that domain-domain motion in EI can contribute as large as ∼-3.2 kcal/mol toward PEP binding.

Conformational selection and substrate binding regulate the monomer/dimer equilibrium of the C-terminal domain of Escherichia coli enzyme I

The bacterial phosphotransferase system (PTS) is a signal transduction pathway that couples phosphoryl transfer to active sugar transport across the cell membrane. The PTS is initiated by the binding of phosphoenolpyruvate (PEP) to the C-terminal domain (EIC) of enzyme I (EI), a highly conserved protein that is common to all sugar branches of the PTS. EIC exists in a dynamic monomer/dimer equilibrium that is modulated by ligand binding and is thought to regulate the overall PTS. Isolation of EIC has proven challenging, and conformational dynamics within the EIC domain during the catalytic cycle are still largely unknown. Here, we present a robust protocol for expression and purification of recombinant EIC from Escherichia coli and show that isolated EIC is capable of hydrolyzing PEP. NMR analysis and residual dipolar coupling measurements indicate that the isolated EIC domain in solution adopts a stable tertiary fold and quaternary structure that is consistent with previously reported crystallographic data. NMR relaxation dispersion measurements indicate that residues around the PEP binding site and in the β3α3 turn (residues 333-366), which is located at the dimer interface, undergo a rapid transition on the sub-millisecond time scale (with an exchange rate constant of ∼1500 s(-1)) between major open (∼97%) and minor closed (∼3%) conformations. Upon PEP binding, the β3α3 turn is effectively locked in the closed state by the formation of salt bridges between the phosphate group of PEP and the side chains of Lys(340) and Arg(358), thereby stabilizing the dimer.

Active site phosphoryl groups in the biphosphorylated phosphotransferase complex reveal dynamics in a millisecond time scale

The N-terminal domain of Enzyme I (EIN) and phosphocarrier HPr can form a biphosphorylated complex when they are both phosphorylated by excess cellular phosphoenolpyruvate. Here we show that the electrostatic repulsion between the phosphoryl groups in the biphosphorylated complex results in characteristic dynamics at the active site in a millisecond time scale. The dynamics is localized to phospho-His15 and the stabilizing backbone amide groups of HPr, and does not impact on the phospho-His189 of EIN. The dynamics occurs with the k(ex) of ~500 s(-1) which compares to the phosphoryl transfer rate of ~850 s(-1) between EIN and HPr. The conformational dynamics in HPr may be important for its phosphotransfer reactions with multiple partner proteins.

Solution structure of the IIAChitobiose-HPr complex of the N,N'-diacetylchitobiose branch of the Escherichia coli phosphotransferase system

The solution structure of the complex of enzyme IIA of the N,N'-diacetylchitobiose (Chb) transporter with the histidine phosphocarrier protein HPr has been solved by NMR. The IIA(Chb)-HPr complex completes the structure elucidation of representative cytoplasmic complexes for all four sugar branches of the bacterial phosphoryl transfer system (PTS). The active site His-89 of IIA(Chb) was mutated to Glu to mimic the phosphorylated state. IIA(Chb)(H89E) and HPr form a weak complex with a K(D) of ~0.7 mM. The interacting binding surfaces, concave for IIA(Chb) and convex for HPr, complement each other in terms of shape, residue type, and charge distribution, with predominantly hydrophobic residues, interspersed by some uncharged polar residues, located centrally, and polar and charged residues at the periphery. The active site histidine of HPr, His-15, is buried within the active site cleft of IIA(Chb) formed at the interface of two adjacent subunits of the IIA(Chb) trimer, thereby coming into close proximity with the active site residue, H89E, of IIA(Chb). A His89-P-His-15 pentacoordinate phosphoryl transition state can readily be modeled without necessitating any significant conformational changes, thereby facilitating rapid phosphoryl transfer. Comparison of the IIA(Chb)-HPr complex with the IIA(Chb)-IIB(Chb) complex, as well as with other cytoplasmic complexes of the PTS, highlights a unifying mechanism for recognition of structurally diverse partners. This involves generating similar binding surfaces from entirely different underlying structural elements, large interaction surfaces coupled with extensive redundancy, and side chain conformational plasticity to optimize diverse sets of intermolecular interactions.

Spatial arrangement of an RNA zipcode identifies mRNAs under post-transcriptional control

How RNA-binding proteins recognize specific sets of target mRNAs remains poorly understood because current approaches depend primarily on sequence information. In this study, we demonstrate that specific recognition of messenger RNAs (mRNAs) by RNA-binding proteins requires the correct spatial positioning of these sequences. We characterized both the cis-acting sequence elements and the spatial restraints that define the mode of RNA binding of the zipcode-binding protein 1 (ZBP1/IMP1/IGF2BP1) to the β-actin zipcode. The third and fourth KH (hnRNP K homology) domains of ZBP1 specifically recognize a bipartite RNA element comprised of a 5' element (CGGAC) followed by a variable 3' element (C/A-CA-C/U) that must be appropriately spaced. Remarkably, the orientation of these elements is interchangeable within target transcripts bound by ZBP1. The spatial relationship of this consensus binding site identified conserved transcripts that were verified to associate with ZBP1 in vivo. The dendritic localization of one of these transcripts, spinophilin, was found to be dependent on both ZBP1 and the RNA elements recognized by ZBP1 KH34.

Assignment strategies for aliphatic protons in the solid-state in randomly protonated proteins

Biological solid-state nuclear magnetic resonance spectroscopy developed rapidly in the past two decades and emerged as an important tool for structural biology. Resonance assignment is an essential prerequisite for structure determination and the characterization of motional properties of a molecule. Experiments, which rely on carbon or nitrogen detection, suffer, however, from low sensitivity. Recently, we introduced the RAP (Reduced Adjoining Protonation) labeling scheme, which allows to detect backbone and sidechain protons with high sensitivity and resolution. We present here a (1)H-detected 3D (H)CCH experiment for assignment of backbone and sidechain proton resonances. Resolution is significantly improved by employing simultaneous (13)CO and (13)Cβ J-decoupling during evolution of the (13)Cα chemical shift. In total, ~90% of the (1)Hα-(13)Cα backbone resonances of chicken α-spectrin SH3 could be assigned.

Automated sequence- and stereo-specific assignment of methyl-labeled proteins by paramagnetic relaxation and methyl-methyl nuclear Overhauser enhancement spectroscopy

Methyl-transverse relaxation optimized spectroscopy is rapidly becoming the preferred NMR technique for probing structure and dynamics of very large proteins up to ~1 MDa in molecular size. Data interpretation, however, necessitates assignment of methyl groups which still presents a very challenging and time-consuming process. Here we demonstrate that, in combination with a known 3D structure, paramagnetic relaxation enhancement (PRE), induced by nitroxide spin-labels incorporated at only a few surface-exposed engineered cysteines, provides fast, straightforward and robust access to methyl group resonance assignments, including stereoassignments for the methyl groups of leucine and valine. Neither prior assignments, including backbone assignments, for the protein, nor experiments that transfer magnetization between methyl groups and the protein backbone, are required. PRE-derived assignments are refined by 4D methyl-methyl nuclear Overhauser enhancement data, eliminating ambiguities and errors that may arise due to the high sensitivity of PREs to the potential presence of sparsely-populated transient states.

Impact of 15N R2/R1 relaxation restraints on molecular size, shape, and bond vector orientation for NMR protein structure determination with sparse distance restraints

(15)N R(2)/R(1) relaxation data contain information on molecular shape and size as well as on bond vector orientations relative to the diffusion tensor. Since the diffusion tensor can be directly calculated from the molecular coordinates, direct inclusion of (15)N R(2)/R(1) restraints in NMR structure calculations without any a priori assumptions is possible. Here we show that (15)N R(2)/R(1) restraints are particularly valuable when only sparse distance restraints are available. Using three examples of proteins of varying size, namely, GB3 (56 residues), ubiquitin (76 residues), and the N-terminal domain of enzyme I (EIN, 249 residues), we show that incorporation of (15)N R(2)/R(1) restraints results in large and significant increases in coordinate accuracy that can make the difference between being able or unable to determine an approximate global fold. For GB3 and ubiquitin, good coordinate accuracy was obtained using only backbone hydrogen-bond restraints supplemented by (15)N R(2)/R(1) relaxation restraints. For EIN, the global fold could be determined using sparse nuclear Overhauser enhancement (NOE) distance restraints involving only NH and methyl groups in conjunction with (15)N R(2)/R(1) restraints. These results are of practical significance in the study of larger and more complex systems, where the increasing spectral complexity and number of chemical shift degeneracies reduce the number of unambiguous NOE assignments that can be readily obtained, resulting in progressively reduced NOE coverage as the size of the protein increases.

Combined use of residual dipolar couplings and solution X-ray scattering to rapidly probe rigid-body conformational transitions in a non-phosphorylatable active-site mutant of the 128 kDa enzyme I dimer

The first component of the bacterial phosphotransferase system, enzyme I (EI), is a multidomain 128 kDa dimer that undergoes large rigid-body conformational transitions during the course of its catalytic cycle. Here we investigate the solution structure of a non-phosphorylatable active-site mutant in which the active-site histidine is substituted by glutamine. We show that perturbations in the relative orientations and positions of the domains and subdomains can be rapidly and reliably determined by conjoined rigid-body/torsion angle/Cartesian simulated annealing calculations driven by orientational restraints from residual dipolar couplings and shape and translation information afforded by small- and wide-angle X-ray scattering. Although histidine and glutamine are isosteric, the conformational space available to a Gln side chain is larger than that for the imidazole ring of His. An additional hydrogen bond between the side chain of Gln189 located on the EIN(α/β) subdomain and an aspartate (Asp129) on the EIN(α) subdomain results in a small (∼9°) reorientation of the EIN(α) and EIN(α/β) subdomains that is in turn propagated to a larger reorientation (∼26°) of the EIN domain relative to the EIC dimerization domain, illustrating the positional sensitivity of the EIN domain and its constituent subdomains to small structural perturbations.

Solution structure of the 128 kDa enzyme I dimer from Escherichia coli and its 146 kDa complex with HPr using residual dipolar couplings and small- and wide-angle X-ray scattering

The solution structures of free Enzyme I (EI, ∼128 kDa, 575 × 2 residues), the first enzyme in the bacterial phosphotransferase system, and its complex with HPr (∼146 kDa) have been solved using novel methodology that makes use of prior structural knowledge (namely, the structures of the dimeric EIC domain and the isolated EIN domain both free and complexed to HPr), combined with residual dipolar coupling (RDC), small- (SAXS) and wide- (WAXS) angle X-ray scattering and small-angle neutron scattering (SANS) data. The calculational strategy employs conjoined rigid body/torsion/Cartesian simulated annealing, and incorporates improvements in calculating and refining against SAXS/WAXS data that take into account complex molecular shapes in the description of the solvent layer resulting in a better representation of the SAXS/WAXS data. The RDC data orient the symmetrically related EIN domains relative to the C(2) symmetry axis of the EIC dimer, while translational, shape, and size information is provided by SAXS/WAXS. The resulting structures are independently validated by SANS. Comparison of the structures of the free EI and the EI-HPr complex with that of the crystal structure of a trapped phosphorylated EI intermediate reveals large (∼70-90°) hinge body rotations of the two subdomains comprising the EIN domain, as well as of the EIN domain relative to the dimeric EIC domain. These large-scale interdomain motions shed light on the structural transitions that accompany the catalytic cycle of EI.

The N-terminal domain of the enzyme I is a monomeric well-folded protein with a low conformational stability and residual structure in the unfolded state

The bacterial phosphoenolpyruvate-dependent sugar phosphotransferase system is a multiprotein complex that phosphorylates and, concomitantly, transports carbohydrates across the membrane into the cell. The first protein of the cascade is a multidomain protein so-called enzyme I (EI). The N-terminal domain of EI from Streptomyces coelicolor, EIN(sc), responsible for the binding to the second protein in the cascade (the histidine phosphocarrier, HPr), was cloned and successfully expressed and purified. We have previously shown that EI(sc) binds to HPr(sc) with smaller affinity than other members of the EI and HPr families [Hurtado-Gómez et al. (2008) Biophys. J., 95, 1336-1348]. We think that the study of the isolated binding HPr(sc) domain, that is EIN(sc), could shed light on the small affinity value measured. Therefore, in this work we present a detailed description of the structural features of the EIN domain, as a first step towards a complete characterization of the molecular recognition process between the two proteins. We show that EIN(sc) is a folded protein, with alpha-helix and beta-sheet structures and also random-coil conformations, as shown by circular dichroism (CD), FTIR and NMR spectroscopies. The acquisition of secondary and tertiary structures, and the burial of hydrophobic regions, occurred concomitantly at acidic pHs, but at very low pH, the domain acquired a molten-globule conformation. The EIN(sc) protein was not very stable, with an apparent conformational free energy change upon unfolding, DeltaG, of 4.1 +/- 0.4 kcal mol(-1), which was pH independent in the range explored (from pH 6.0 to 8.5). The thermal denaturation midpoint, which was also pH invariant, was similar to that measured in the isolated intact EI(sc). Although EIN(sc) shows thermal- and chemical denaturations that seems to follow a two-state mechanism, there is evidence of residual structure in the chemical and thermally unfolded states, as indicated by differential scanning calorimetry and CD measurements.

Molecular basis for barbed end uncapping by CARMIL homology domain 3 of mouse CARMIL-1

Capping protein (CP) is a ubiquitously expressed, 62-kDa heterodimer that binds the barbed end of the actin filament with approximately 0.1 nm affinity to prevent further monomer addition. CARMIL is a multidomain protein, present from protozoa to mammals, that binds CP and is important for normal actin dynamics in vivo. The CARMIL CP binding site resides in its CAH3 domain (CARMIL homology domain 3) located at or near the protein's C terminus. CAH3 binds CP with approximately 1 nm affinity, resulting in a complex with weak capping activity (30-200 nm). Solution assays and single-molecule imaging show that CAH3 binds CP already present on the barbed end, causing a 300-fold increase in the dissociation rate of CP from the end (i.e. uncapping). Here we used nuclear magnetic resonance (NMR) to define the molecular interaction between the minimal CAH3 domain (CAH3a/b) of mouse CARMIL-1 and CP. Specifically, we show that the highly basic CAH3a subdomain is required for the high affinity interaction of CAH3 with a complementary "acidic groove" on CP opposite its actin-binding surface. This CAH3a-CP interaction orients the CAH3b subdomain, which we show is also required for potent anti-CP activity, directly adjacent to the basic patch of CP, shown previously to be required for CP association to and high affinity interaction with the barbed end. The importance of specific residue interactions between CP and CAH3a/b was confirmed by site-directed mutagenesis of both proteins. Together, these results offer a mechanistic explanation for the barbed end uncapping activity of CARMIL, and they identify the basic patch on CP as a crucial regulatory site.

SAIL – stereo-array isotope labeling

Optimal stereospecific and regiospecific labeling of proteins with stable isotopes enhances the nuclear magnetic resonance (NMR) method for the determination of the three-dimensional protein structures in solution. Stereo-array isotope labeling (SAIL) offers sharpened lines, spectral simplification without loss of information and the ability to rapidly collect and automatically evaluate the structural restraints required to solve a high-quality solution structure for proteins up to twice as large as before. This review gives an overview of stable isotope labeling methods for NMR spectroscopy with proteins and provides an in-depth treatment of the SAIL technology.

A new model for mapping the peptide backbone: predicting proton chemical shifts in proteins

This paper describes a methodology that correlates experimental chemical shifts (at the alpha proton) of proteins with their geometrical data (both dihedral angles and distances) obtained from 13 representative proteins, which are taken from the Protein Data Bank (PDB) and the BioMagRes Data Bank (BMRB). To this end, the experimentally measured proton chemical shifts of simple amides have been correlated with DFT-based calculated structures, at the B3PW91/6-31G* level. This results in a series of mathematical relationships, which are extrapolated to the above-mentioned proteins giving rise to a modified equation for such skeleta. It is relevant to note that the equation is also supported by a clear comparison with NMR data of a protein beyond the chosen set, such as insulin, even with lower errors. The model also relates the dependence of chemical shifts on hydrophobic and anisotropic effects at the amino acid residues.

Solution structure of the IIAChitobiose-IIBChitobiose complex of the N,N'-diacetylchitobiose branch of the Escherichia coli phosphotransferase system

The solution structure of the IIA-IIB complex of the N,N'-diacetylchitobiose (Chb) transporter of the Escherichia coli phosphotransferase system has been solved by NMR. The active site His-89 of IIA(Chb) was mutated to Glu to mimic the phosphorylated state and the active site Cys-10 of IIB(Chb) was substituted by serine to prevent intermolecular disulfide bond formation. Binding is weak with a K(D) of approximately 1.3 mm. The two complementary interaction surfaces are largely hydrophobic, with the protruding active site loop (residues 9-16) of IIB(Chb) buried deep within the active site cleft formed at the interface of two adjacent subunits of the IIA(Chb) trimer. The central hydrophobic portion of the interface is surrounded by a ring of polar and charged residues that provide a relatively small number of electrostatic intermolecular interactions that serve to correctly align the two proteins. The conformation of the active site loop in unphosphorylated IIB(Chb) is inconsistent with the formation of a phosphoryl transition state intermediate because of steric hindrance, especially from the methyl group of Ala-12 of IIB(Chb). Phosphorylation of IIB(Chb) is accompanied by a conformational change within the active site loop such that its path from residues 11-13 follows a mirror-like image relative to that in the unphosphorylated state. This involves a transition of the phi/psi angles of Gly-13 from the right to left alpha-helical region, as well as smaller changes in the backbone torsion angles of Ala-12 and Met-14. The resulting active site conformation is fully compatible with the formation of the His-89-P-Cys-10 phosphoryl transition state without necessitating any change in relative translation or orientation of the two proteins within the complex.

Robust structure-based resonance assignment for functional protein studies by NMR

High-throughput functional protein NMR studies, like protein interactions or dynamics, require an automated approach for the assignment of the protein backbone. With the availability of a growing number of protein 3D structures, a new class of automated approaches, called structure-based assignment, has been developed quite recently. Structure-based approaches use primarily NMR input data that are not based on J-coupling and for which connections between residues are not limited by through bonds magnetization transfer efficiency. We present here a robust structure-based assignment approach using mainly H(N)-H(N) NOEs networks, as well as (1)H-(15) N residual dipolar couplings and chemical shifts. The NOEnet complete search algorithm is robust against assignment errors, even for sparse input data. Instead of a unique and partly erroneous assignment solution, an optimal assignment ensemble with an accuracy equal or near to 100% is given by NOEnet. We show that even low precision assignment ensembles give enough information for functional studies, like modeling of protein-complexes. Finally, the combination of NOEnet with a low number of ambiguous J-coupling sequential connectivities yields a high precision assignment ensemble. NOEnet will be available under: http://www.icsn.cnrs-gif.fr/download/nmr.

NMR reveals pathway for ferric mineral precursors to the central cavity of ferritin

Ferritin is a multimeric nanocage protein that directs the reversible biomineralization of iron. At the catalytic ferroxidase site two iron(II) ions react with dioxygen to form diferric species. In order to study the pathway of iron(III) from the ferroxidase site to the central cavity a new NMR strategy was developed to manage the investigation of a system composed of 24 monomers of 20 kDa each. The strategy is based on (13)C-(13)C solution NOESY experiments combined with solid-state proton-driven (13)C-(13)C spin diffusion and 3D coherence transfer experiments. In this way, 75% of amino acids were recognized and 35% sequence-specific assigned. Paramagnetic broadening, induced by iron(III) species in solution (13)C-(13)C NOESY spectra, localized the iron within each subunit and traced the progression to the central cavity. Eight iron ions fill the 20-A-long iron channel from the ferrous/dioxygen oxidoreductase site to the exit into the cavity, inside the four-helix bundle of each subunit, contrasting with short paths in models. Magnetic susceptibility data support the formation of ferric multimers in the iron channels. Multiple iron channel exits are near enough to facilitate high concentration of iron that can mineralize in the ferritin cavity, illustrating advantages of the multisubunit cage structure.

Using the experimentally determined components of the overall rotational diffusion tensor to restrain molecular shape and size in NMR structure determination of globular proteins and protein-protein complexes

This paper describes an approach for making use of the components of the experimentally determined rotational diffusion tensor derived from NMR relaxation measurements in macromolecular structure determination. The parameters of the rotational diffusion tensor describe the shape and size of the macromolecule or macromolecular complex, and are therefore complementary to traditional NMR restraints. The structural information contained in the rotational diffusion tensor is not dissimilar to that present in the small-angle region of solution X-ray scattering profiles. We demonstrate the utility of rotational diffusion tensor restraints for protein structure refinement using the N-terminal domain of enzyme I (EIN) as an example and validate the results by solution small-angle X-ray scattering. We also show how rotational diffusion tensor restraints can be used for docking complexes using the dimeric HIV-1 protease and the EIN-HPr complexes as examples. In the former case, the rotational diffusion tensor restraints are sufficient in their own right to determine the position of one subunit relative to another. In the latter case, rotational diffusion tensor restraints complemented by highly ambiguous distance restraints derived from chemical shift perturbation mapping and a hydrophobic contact potential are sufficient to correctly dock EIN to HPr. In each case, the cluster containing the lowest-energy structure corresponds to the correct solution.