An improved method for distinguishing between anisotropic tumbling and chemical exchange in analysis of 15N relaxation parameters

Alternative low-populated conformations prompt phase transitions in polyalanine repeat expansions

Abnormal trinucleotide repeat expansions alter protein conformation causing malfunction and contribute to a significant number of incurable human diseases. Scarce structural insights available on disease-related homorepeat expansions hinder the design of effective therapeutics. Here, we present the dynamic structure of human PHOX2B C-terminal fragment, which contains the longest polyalanine segment known in mammals. The major α-helical conformation of the polyalanine tract is solely extended by polyalanine expansions in PHOX2B, which are responsible for most congenital central hypoventilation syndrome cases. However, polyalanine expansions in PHOX2B additionally promote nascent homorepeat conformations that trigger length-dependent phase transitions into solid condensates that capture wild-type PHOX2B. Remarkably, HSP70 and HSP90 chaperones specifically seize PHOX2B alternative conformations preventing phase transitions. The precise observation of emerging polymorphs in expanded PHOX2B postulates unbalanced phase transitions as distinct pathophysiological mechanisms in homorepeat expansion diseases, paving the way towards the search of therapeutics modulating biomolecular condensates in central hypoventilation syndrome.

Metamorphism in TDP-43 prion-like domain determines chaperone recognition

The RNA binding protein TDP-43 forms cytoplasmic inclusions via its C-terminal prion-like domain in several neurodegenerative diseases. Aberrant TDP-43 aggregation arises upon phase de-mixing and transitions from liquid to solid states, following still unknown structural conversions which are primed by oxidative stress and chaperone inhibition. Despite the well-established protective roles for molecular chaperones against protein aggregation pathologies, knowledge on the determinants of chaperone recognition in disease-related prions is scarce. Here we show that chaperones and co-chaperones primarily recognize the structured elements in TDP-43´s prion-like domain. Significantly, while HSP70 and HSP90 chaperones promote TDP-43 phase separation, co-chaperones from the three classes of the large human HSP40 family (namely DNAJA2, DNAJB1, DNAJB4 and DNAJC7) show strikingly different effects on TDP-43 de-mixing. Dismantling of the second helical element in TDP-43 prion-like domain by methionine sulfoxidation impacts phase separation and amyloid formation, abrogates chaperone recognition and alters phosphorylation by casein kinase-1δ. Our results show that metamorphism in the post-translationally modified TDP-43 prion-like domain encodes determinants that command mechanisms with major relevance in disease.

Conformational Priming of RepA-WH1 for Functional Amyloid Conversion Detected by NMR Spectroscopy

How proteins with a stable globular fold acquire the amyloid state is still largely unknown. RepA, a versatile plasmidic DNA binding protein from Pseudomonas savastanoi, is functional as a transcriptional repressor or as an initiator or inhibitor of DNA replication, the latter via assembly of an amyloidogenic oligomer. Its N-terminal domain (WH1) is responsible for discrimination between these functional abilities by undergoing insufficiently understood structural changes. RepA-WH1 is a stable dimer whose conformational dynamics had not been explored. Here, we have studied it through NMR {¹H}-¹⁵N relaxation and H/D exchange kinetics measurements. The N- and the C-terminal α-helices, and the internal amyloidogenic loop, are partially unfolded in solution. S4-indigo, a small inhibitor of RepA-WH1 amyloidogenesis, binds to and tethers the N-terminal α-helix to a β-hairpin that is involved in dimerization, thus providing evidence for a priming role of fraying ends and dimerization switches in the amyloidogenesis of folded proteins.

Solution structure of the autophagy-related protein LC3C reveals a polyproline II motif on a mobile tether with phosphorylation site

(Macro-)autophagy is a compartmental degradation pathway conserved from yeast to mammals. The yeast protein Atg8 mediates membrane tethering/hemifusion and cargo recruitment and is essential for autophagy. The human MAP1LC3/GABARAP family proteins show high sequence identity with Atg8, but MAP1LC3C is distinguished by a conspicuous amino-terminal extension with unknown functional significance. We have determined the high-resolution three-dimensional structure and measured the backbone dynamics of MAP1LC3C by NMR spectroscopy. From Ser18 to Ala120, MAP1LC3C forms an α-helix followed by the ubiquitin-like tertiary fold with two hydrophobic binding pockets used by MAP1LC3/GABARAP proteins to recognize targets presenting LC3-interacting regions (LIRs). Unlike other MAP1LC3/GABARAP proteins, the amino-terminal region of MAP1LC3C does not form a stable helix α₁ but a "sticky arm" consisting of a polyproline II motif on a flexible linker. Ser18 at the interface between this linker and the structural core can be phosphorylated in vitro by protein kinase A, which causes additional conformational heterogeneity as monitored by NMR spectroscopy and molecular dynamics simulations, including changes in the LIR-binding interface. Based on these results we propose that the amino-terminal polyproline II motif mediates specific interactions with the microtubule cytoskeleton and that Ser18 phosphorylation modulates the interplay of MAP1LC3C with its various target proteins.

Measuring Entropy in Molecular Recognition by Proteins

Molecular recognition by proteins is fundamental to the molecular basis of biology. Dissection of the thermodynamic landscape governing protein-ligand interactions has proven difficult because determination of various entropic contributions is quite challenging. Nuclear magnetic resonance relaxation measurements, theory, and simulations suggest that conformational entropy can be accessed through a dynamical proxy. Here, we review the relationship between measures of fast side-chain motion and the underlying conformational entropy. The dynamical proxy reveals that the contribution of conformational entropy can range from highly favorable to highly unfavorable and demonstrates the potential of this key thermodynamic variable to modulate protein-ligand interactions. The dynamical so-called entropy meter also refines the role of solvent entropy and directly determines the loss in rotational-translational entropy that occurs upon formation of high-affinity complexes. The ability to quantify the roles of entropy through an entropy meter based on measurable dynamical properties promises to highlight its role in protein function.

Methods of probing the interactions between small molecules and disordered proteins

It is generally recognized that a large fraction of the human proteome is made up of proteins that remain disordered in their native states. Despite the fact that such proteins play key biological roles and are involved in many major human diseases, they still represent challenging targets for drug discovery. A major bottleneck for the identification of compounds capable of interacting with these proteins and modulating their disease-promoting behaviour is the development of effective techniques to probe such interactions. The difficulties in carrying out binding measurements have resulted in a poor understanding of the mechanisms underlying these interactions. In order to facilitate further methodological advances, here we review the most commonly used techniques to probe three types of interactions involving small molecules: (1) those that disrupt functional interactions between disordered proteins; (2) those that inhibit the aberrant aggregation of disordered proteins, and (3) those that lead to binding disordered proteins in their monomeric states. In discussing these techniques, we also point out directions for future developments.

Point mutations in the N-terminal domain of transactive response DNA-binding protein 43 kDa (TDP-43) compromise its stability, dimerization, and functions

Transactive response DNA-binding protein 43 (TDP-43) performs multiple tasks in mRNA processing, transport, and translational regulation, but it also forms aggregates implicated in amyotrophic lateral sclerosis. TDP-43's N-terminal domain (NTD) is important for these activities and dysfunctions; however, there is an open debate about whether or not it adopts a specifically folded, stable structure. Here, we studied NTD mutations designed to destabilize its structure utilizing NMR and fluorescence spectroscopies, analytical ultracentrifugation, splicing assays, and cell microscopy. The substitutions V31R and T32R abolished TDP-43 activity in splicing and aggregation processes, and even the rather mild L28A mutation severely destabilized the NTD, drastically reducing TDP-43's in vitro splicing activity and inducing aberrant localization and aggregation in cells. These findings strongly support the idea that a stably folded NTD is essential for correct TDP-43 function. The stably folded NTD also promotes dimerization, which is pertinent to the protein's activities and pathological aggregation, and we present an atomic-level structural model for the TDP-43 dimer based on NMR data. Leu-27 is evolutionarily well conserved even though it is exposed in the monomeric NTD. We found here that Leu-27 is buried in the dimer and that the L27A mutation promotes monomerization. In conclusion, our study sheds light on the structural and biological properties of the TDP-43 NTD, indicating that the NTD must be stably folded for TDP-43's physiological functions, and has implications for understanding the mechanisms promoting the pathological aggregation of this protein.

Spider wrapping silk fibre architecture arising from its modular soluble protein precursor

Spiders store spidroins in their silk glands as high concentration aqueous solutions, spinning these dopes into fibres with outstanding mechanical properties. Aciniform (or wrapping) silk is the toughest spider silk and is devoid of the short amino acid sequence motifs characteristic of the other spidroins. Using solution-state NMR spectroscopy, we demonstrate that the 200 amino acid Argiope trifasciata AcSp1 repeat unit contrasts with previously characterized spidroins, adopting a globular 5-helix bundle flanked by intrinsically disordered N- and C-terminal tails. Split-intein-mediated segmental NMR-active isotope-enrichment allowed unambiguous demonstration of modular and malleable "beads-on-a-string" concatemeric behaviour. Concatemers form fibres upon manual drawing with silk-like morphology and mechanical properties, alongside secondary structuring and orientation consistent with native AcSp1 fibres. AcSp1 structural stability varies locally, with the fifth helix denaturing most readily. The structural transition of aciniform spidroin from a mostly α-helical dope to a mixed α-helix/β-sheet-containing fibre can be directly related to spidroin architecture and stability.

Back to the Future: Nuclear Magnetic Resonance and Bioinformatics Studies on Intrinsically Disordered Proteins

From the 1970s to the present, regions of missing electron density in protein structures determined by X-ray diffraction and the characterization of the functions of these regions have suggested that not all protein regions depend on prior 3D structure to carry out function. Motivated by these observations, in early 1996 we began to use bioinformatics approaches to study these intrinsically disordered proteins (IDPs) and IDP regions. At just about the same time, several laboratory groups began to study a collection of IDPs and IDP regions using nuclear magnetic resonance. The temporal overlap of the bioinformatics and NMR studies played a significant role in the development of our understanding of IDPs. Here the goal is to recount some of this history and to project from this experience possible directions for future work.

Structural and functional studies of Stf76 from the Sulfolobus islandicus plasmid-virus pSSVx: a novel peculiar member of the winged helix-turn-helix transcription factor family

The hybrid plasmid–virus pSSVx from Sulfolobus islandicus presents an open reading frame encoding a 76 amino acid protein, namely Stf76, that does not show significant sequence homology with any protein with known 3D structure. The recombinant protein recognizes specifically two DNA-binding sites located in its own promoter, thus suggesting an auto-regulated role of its expression. Circular dichroism, spectrofluorimetric, light scattering and isothermal titration calorimetry experiments indicated a 2:1 molar ratio (protein:DNA) upon binding to the DNA target containing a single site. Furthermore, the solution structure of Stf76, determined by nuclear magnetic resonance (NMR) using chemical shift Rosetta software, has shown that the protein assumes a winged helix–turn–helix fold. NMR chemical shift perturbation analysis has been performed for the identification of the residues responsible for DNA interaction. In addition, a model of the Stf76–DNA complex has been built using as template a structurally related homolog.

A surprising role for conformational entropy in protein function

Formation of high-affinity complexes is critical for the majority of enzymatic reactions involving proteins. The creation of the family of Michaelis and other intermediate complexes during catalysis clearly involves a complicated manifold of interactions that are diverse and complex. Indeed, computing the energetics of interactions between proteins and small molecule ligands using molecular structure alone remains a great challenge. One of the most difficult contributions to the free energy of protein-ligand complexes to access experimentally is that due to changes in protein conformational entropy. Fortunately, recent advances in solution nuclear magnetic resonance (NMR) relaxation methods have enabled the use of measures-of-motion between conformational states of a protein as a proxy for conformational entropy. This review briefly summarizes the experimental approaches currently employed to characterize fast internal motion in proteins, how this information is used to gain insight into conformational entropy, what has been learned, and what the future may hold for this emerging view of protein function.

An assignment of intrinsically disordered regions of proteins based on NMR structures

Intrinsically disordered proteins (IDPs) do not adopt stable three-dimensional structures in physiological conditions, yet these proteins play crucial roles in biological phenomena. In most cases, intrinsic disorder manifests itself in segments or domains of an IDP, called intrinsically disordered regions (IDRs), but fully disordered IDPs also exist. Although IDRs can be detected as missing residues in protein structures determined by X-ray crystallography, no protocol has been developed to identify IDRs from structures obtained by Nuclear Magnetic Resonance (NMR). Here, we propose a computational method to assign IDRs based on NMR structures. We compared missing residues of X-ray structures with residue-wise deviations of NMR structures for identical proteins, and derived a threshold deviation that gives the best correlation of ordered and disordered regions of both structures. The obtained threshold of 3.2Å was applied to proteins whose structures were only determined by NMR, and the resulting IDRs were analyzed and compared to those of X-ray structures with no NMR counterpart in terms of sequence length, IDR fraction, protein function, cellular location, and amino acid composition, all of which suggest distinct characteristics. The structural knowledge of IDPs is still inadequate compared with that of structured proteins. Our method can collect and utilize IDRs from structures determined by NMR, potentially enhancing the understanding of IDPs.

Guanine nucleotides differentially modulate backbone dynamics of the STAS domain of the SulP/SLC26 transport protein Rv1739c of Mycobacterium tuberculosis

Enzymatic catalysis and protein signaling are dynamic processes that involve local and/or global conformational changes occurring across a broad range of time scales. (1) H-(15) N relaxation NMR provides a comprehensive understanding of protein backbone dynamics both in the apo (unliganded) and ligand-bound conformations, enabling both fast and slow internal motions of individual amino acid residues to be observed. We recently reported the structure and nucleotide binding properties of the sulfate transporter and anti-sigma factor antagonist (STAS) domain of Rv1739c, a SulP anion transporter protein of Mycobacterium tuberculosis. In the present study, we report (1) H-(15) N NMR backbone dynamics measurements [longitudinal (T(1) ), transverse (T(2) ) and steady-state ({(1) H}-(15) N) heteronuclear NOE] of the Rv1739c STAS domain, in the absence and presence of saturating concentrations of GTP and GDP. Analysis of measured relaxation data and estimated dynamic parameters indicated distinct features differentiating the binding of GTP and GDP to Rv1739c STAS. The 9.55 ns overall rotational correlation time of Rv1739c STAS increased to 10.48 ns in the presence of GTP, and to 13.25 ns in the presence of GDP, indicating significant nucleotide-induced conformational changes. These conformational changes were accompanied by slow time scale (μs to ms) motions in discrete regions of the protein, as reflected by guanine nucleotide-induced changes in relaxation parameters. The observed nucleotide-specific alterations in the relaxation properties of individual STAS residues reflect an increased molecular anisotropy and/or the emergence of conformational equilibria governing functional properties of the STAS domain.

Internal dynamics of the tryptophan repressor (TrpR) and two functionally distinct TrpR variants, L75F-TrpR and A77V-TrpR, in their l-Trp-bound forms

Backbone amide dynamics of the Escherichia coli tryptophan repressor protein (WT-TrpR) and two functionally distinct variants, L75F-TrpR and A77V-TrpR, in their holo (l-tryptophan corepressor-bound) form have been characterized using (15)N nuclear magnetic resonance (NMR) relaxation. The three proteins possess very similar structures, ruling out major conformational differences as the source of their functional differences, and suggest that changes in protein flexibility are at the origin of their distinct functional properties. Comparison of site specific (15)N-T(1), (15)N-T(2), (15)N-{(1)H} nuclear Overhauser effect, reduced spectral density, and generalized order (S(2)) parameters indicates that backbone dynamics in the three holo-repressors are overall very similar with a few notable and significant exceptions for backbone atoms residing within the proteins' DNA-binding domain. We find that flexibility is highly restricted for amides in core α-helices (i.e., helices A-C and F), and a comparable "stiffening" is observed for residues in the DNA recognition helix (helix E) of the helix D-turn-helix E (HTH) DNA-binding domain of the three holo-repressors. Unexpectedly, amides located in helix D and in adjacent turn regions remain flexible. These data support the concept that residual flexibility in TrpR is essential for repressor function, DNA binding, and molecular recognition of target operators. Comparison of the (15)N NMR relaxation parameters of the holo-TrpRs with those of the apo-TrpRs indicates that the single-point amino acid substitutions, L75F and A77V, perturb the flexibility of backbone amides of TrpR in very different ways and are most pronounced in the apo forms of the three repressors. Finally, we present these findings in the context of other DNA-binding proteins and the role of protein flexibility in molecular recognition.

Apparent tradeoff of higher activity in MMP-12 for enhanced stability and flexibility in MMP-3

The greater activity of MMP-12 than MMP-3 toward substrates from protein fibrils has been quantified. Why is MMP-12 the more active protease? We looked for behaviors associated with the higher activity of MMP-12 than MMP-3, using nuclear magnetic resonance to monitor backbone dynamics and residue-specific stabilities of their catalytic domain. The proteolytic activities are likely to play important roles in inflammatory diseases of arteries, lungs, joints, and intestines. Nuclear magnetic resonance line broadening indicates that regions surrounding the active sites of both proteases sample conformational substates within milliseconds. The more extensive line broadening in MMP-3 suggests greater sampling of conformational substates, affecting the full length of helix B and beta-strand IV forming the active site, and more remote sites. This could suggest more excursions to functionally incompetent substates. MMP-3 also has enhanced subnanosecond fluctuations in helix A, in the beta-hairpin of strands IV and V, and before and including helix C. Hydrogen exchange protection in the EX2 regime suggests that MMP-3 possesses 2.8 kcal/mol higher folding stability than MMP-12(E219A). The beta-sheet of MMP-3 appears to be stabilized still more. The higher stability of MMP-3 relative to MMP-12 coincides with the former's considerably lower proteolytic activity. This relationship is consistent with the hypothesis that enzymes often trade stability for higher activity.

Backbone amide dynamics studies of Apo-L75F-TrpR, a temperature-sensitive mutant of the tryptophan repressor protein (TrpR): comparison with the (15)N NMR relaxation profiles of wild-type and A77V mutant Apo-TrpR repressors

Backbone amide dynamics studies were conducted on a temperature-sensitive mutant (L75F-TrpR) of the tryptophan repressor protein (TrpR) of Escherichia coli in its apo (i.e., no l-tryptophan corepressor-bound) form. The (15)N NMR relaxation profiles of apo-L75F-TrpR were analyzed and compared to those of wild-type (WT) and super-repressor mutant (A77V) TrpR proteins, also in their apo forms. The (15)N NMR relaxation data ((15)N-T(1), (15)N-T(2), and heteronuclear (15)N-{(1)H}-nOe) recorded on all three aporepressors at a magnetic field strength of 600 MHz ((1)H Larmor frequency) were analyzed to extract dynamics parameters, including diffusion tensor ratios (D(∥)/D(⊥)), correlation times (τ(m)) for overall reorientations of the proteins in solution, reduced spectral density terms [J(eff)(0), J(0.87ω(H)), J(ω(N))], and generalized order parameters (S(2)), which report on protein internal motions on the picosecond to nanosecond and slower microsecond to millisecond chemical exchange time scales. Our results indicate that all three aporepressors exhibit comparable D(∥)/D(⊥) ratios and characteristic time constants, τ(m), for overall global reorientation, indicating that in solution, all three apoproteins display very similar overall shape, structure, and rotational diffusion properties. Comparison of (15)N NMR relaxation data, reduced spectral density profiles, and generalized S(2) order parameters indicated that these parameters are quite uniform for backbone amides positioned within the four (A-C and F) core α-helices of all three aporepressors. In contrast, small but noticeable differences in internal dynamics were observed for backbone amides located within the helix D-turn-helix E DNA-binding domain of the apo-TrpR proteins. The significance of these dynamics differences in terms of the biophysical characteristics and ligand binding properties of the three apo-TrpR proteins is discussed.

NMR backbone dynamics studies of human PED/PEA-15 outline protein functional sites

PED/PEA-15 (phosphoprotein enriched in diabetes/phosphoprotein enriched in astrocytes) is a ubiquitously expressed protein and a key regulator of cell growth and glucose metabolism. PED/PEA-15 mediates both homotypic and heterotypic interactions and is constituted by an N-terminal canonical death effector domain and a C-terminal tail. In the present study, the backbone dynamics of PED/PEA-15 via (15)N R(1) and R(2) and steady-state [(1)H]-(15)N NOE measurements is reported. The dynamic parameters were analyzed using both Lipari-Szabo model-free formalism and a reduced spectral density mapping approach. The results obtained define a polar and charged surface of the death effector domain characterized by internal motions in the micro- to millisecond timescale, which is crucial for the multiple heterotypic functional protein-protein interactions in which PED/PEA-15 is involved. The present study contributes to a better understanding of the molecular basis of the PED/PEA-15 functional interactions and provides a more detailed surface for the design and development of PED/PEA-15 binders.

Human RegIV protein adopts a typical C-type lectin fold but binds mannan with two calcium-independent sites

Human RegIV protein, which contains a sequence motif homologous to calcium-dependent (C-type) lectin-like domain, is highly expressed in mucosa cells of the gastrointestinal tract during pathogen infection and carcinogenesis and may be applied in both diagnosis and treatment of gastric and colon cancers. Here, we provide evidence that, unlike other C-type lectins, human RegIV binds to polysaccharides, mannan, and heparin in the absence of calcium. To elucidate the structural basis for carbohydrate recognition by NMR, we generated the mutant with Pro91 replaced by Ser (hRegIV-P91S) and showed that the structural property and carbohydrate binding ability of hRegIV-P91S are almost identical with those of wild-type protein. The solution structure of hRegIV-P91S was determined, showing that it adopts a typical fold of C-type lectin. Based on the chemical shift perturbations of amide resonances, two calcium-independent mannan-binding sites were proposed. One site is similar to the calcium-independent sugar-binding site on human RegIII and Langerin. Interestingly, the other site is adjacent to the conserved calcium-dependent site at position Ca-2 of typical C-type lectins. Moreover, model-free analysis of (15)N relaxation parameters and simplified Carr-Purcell-Meiboom-Gill relaxation dispersion experiments showed that a slow microsecond-to-millisecond time-scale backbone motion is involved in mannan binding by this site, suggesting a potential role for specific carbohydrate recognition. Our findings shed light on the sugar-binding mode of Reg family proteins, and we postulate that Reg family proteins evolved to bind sugar without calcium to keep the carbohydrate recognition activity under low-pH environments in the gastrointestinal tract.

Conformational dynamics and structural plasticity play critical roles in the ubiquitin recognition of a UIM domain

Ubiquitin-interacting motifs (UIMs) are an important class of protein domains that interact with ubiquitin or ubiquitin-like proteins. These approximately 20-residue-long domains are found in a variety of ubiquitin receptor proteins and serve as recognition modules towards intracellular targets, which may be individual ubiquitin subunits or polyubiquitin chains attached to a variety of proteins. Previous structural studies of interactions between UIMs and ubiquitin have shown that UIMs adopt an extended structure of a single alpha-helix, containing a hydrophobic surface with a conserved sequence pattern that interacts with key hydrophobic residues on ubiquitin. In light of this large body of structural studies, details regarding the presence and the roles of structural dynamics and plasticity are surprisingly lacking. In order to better understand the structural basis of ubiquitin-UIM recognition, we have characterized changes in the structure and dynamics of ubiquitin upon binding of a UIM domain from the yeast Vps27 protein. The solution structure of a ubiquitin-UIM fusion protein designed to study these interactions is reported here and found to consist of a well-defined ubiquitin core and a bipartite UIM helix. Moreover, we have studied the plasticity of the docking interface, as well as global changes in ubiquitin due to UIM binding at the picoseconds-to-nanoseconds and microseconds-to-milliseconds protein motions by nuclear magnetic resonance relaxation. Changes in generalized-order parameters of amide groups show a distinct trend towards increased structural rigidity at the UIM-ubiquitin interface relative to values determined in unbound ubiquitin. Analysis of (15)N Carr-Purcell-Meiboom-Gill relaxation dispersion measurements suggests the presence of two types of motions: one directly related to the UIM-binding interface and the other induced to distal parts of the protein. This study demonstrates a case where localized interactions among protein domains have global effects on protein motions at timescales ranging from picoseconds to milliseconds.

General theoretical/computational tool for interpreting NMR spin relaxation in proteins

We developed in recent years the slowly relaxing local structure (SRLS) approach for analyzing NMR spin relaxation in proteins. SRLS is a two-body coupled rotator model which accounts rigorously for mode-coupling between the global motion of the protein and the local motion of the spin-bearing probe and allows for general properties of the second rank tensors involved. We showed that a general tool of data analysis requires both capabilities. Several important functionalities were missing in our previous implementations of SRLS in data fitting schemes, and in some important cases, the calculations were tedious. Here we present a general implementation which allows for asymmetric local and global diffusion tensors, distinct local ordering and local diffusion frames, and features a rhombic local potential which includes Wigner matrix element terms of ranks 2 and 4. A recently developed hydrodynamics-based approach for calculating global diffusion tensors has been incorporated into the data-fitting scheme. The computational efficiency of the latter has been increased significantly through object-oriented programming within the scope of the C++ programming language, and code parallelization. A convenient graphical user interface is provided. Currently autocorrelated (15)N spin relaxation data can be analyzed effectively. Adaptation to any autocorrelated and cross-correlated relaxation analysis is straightforward. New physical insight is gleaned on largely preserved local structure in solution, even in chain segments which experience slow local motion. Prospects associated with improved dynamic models, and new applications made possible by the current implementation of SRLS, are delineated.

Evidence for domain motion in proteins affecting global diffusion properties: a nuclear magnetic resonance study

The rotational diffusion of proteins is an important hydrodynamic property. Compact protein structures were found previously to exhibit hydration layer viscosity, etaloc, higher than the viscosity of bulk water, eta. This implies an apparent activation energy for rotational diffusion higher than the activation energy of water viscosity, Eeta=15.4+/-0.3 kJ/mol. In this study we examine etaloc of internally mobile proteins using 15N spin relaxation methods. We also examine the activation enthalpy, DeltaH#, and activation entropy, DeltaS#, for rotational diffusion. Of particular relevance are internally mobile ligand-free forms and compact ligand-bound forms of multidomain proteins. Adenylate kinase (AKeco) and Ca2+-calmodulin (Ca2+-CaM) are typical examples. For AKeco (Ca2+-CaM) we find that DeltaH# is 14.5+/-0.5 (15.7+/-0.4) kJ/mol. For the complex of AKeco with the inhibitor AP5A (the complex of Ca2+-CaM with the peptide smMLCKp), we find that DeltaH# is 18.1+/-0.7 (18.2+/-0.5) kJ/mol. The internally mobile outer surface protein A has DeltaH#=12.6+/-0.8 kJ/mol, and the compact protein Staphylococcal nuclease has DeltaH#=18.8+/-0.6 kJ/mol. For the internally mobile and compact proteins studied, <|DeltaS(|> equals 62+/-7 J/(mol K) and 44+/-5 J/(mol K), respectively. The fact is that etaloc>eta (DeltaH#>Eeta) for compact proteins was ascribed previously to electrostatic interactions between surface sites and water rigidifying the hydration layer. We find herein that obliteration of these interactions by domain motion leads to etaloc approximately eta, DeltaH# approximately Eeta, and large activation entropy for internally mobile protein structures.

Redox-dependent dynamics of a dual thioredoxin fold protein: evolution of specialized folds

An enzyme system protecting bacteria from oxidative stress includes the flavoprotein AhpF and the peroxiredoxin AhpC. The N-terminal domain of AhpF (NTD), with two fused thioredoxin (Trx) folds, belongs to the hyperthermophilic protein disulfide oxidoreductase family. The NTD is distinct in that it contains a redox active a fold with a CxxC sequence and a redox inactive b fold that has lost the CxxC motif. Here we characterize the stability, the (15)N backbone relaxation, and the hydrogen-deuterium exchange properties of reduced [NTD-(SH)(2)] and oxidized (NTD-S(2)) NTD from Salmonella typhimurium. While both NTD-(SH)(2) and NTD-S(2) exhibit similar equilibrium unfolding transitions and order parameters, R(ex) relaxation terms are quite distinct with considerably more intermediate time scale motions in NTD-S(2). Hydrogen exchange protection factors show that the slowly exchanging core corresponds to residues in the b fold in both NTD-(SH)(2) and NTD-S(2). Interestingly, folded-state dynamic fluctuations in the catalytic a fold are significantly increased for residues in NTD-S(2) compared to NTD-(SH)(2). Taken together, these data demonstrate that oxidation of the active site disulfide does not significantly increase stability but results in a dramatic increase in conformational heterogeneity in residues primarily in the redox active a fold. Differences in dynamics between the two folds of the NTD suggest that each evolved a specialized function which, in the a fold, couples redox state to internal motions which may enhance catalysis and specificity and, in the b fold, provides a redox insensitive stable core.

Phosphorylation-induced changes in backbone dynamics of the dematin headpiece C-terminal domain

Dematin is an actin-binding protein abundant in red blood cells and other tissues. It contains a villin-type 'headpiece' F-actin-binding domain at its extreme C-terminus. The isolated dematin headpiece domain (DHP) undergoes a significant conformational change upon phosphorylation. The mutation of Ser74 to Glu closely mimics the phosphorylation of DHP. We investigated motions in the backbone of DHP and its mutant DHPS74E using several complementary NMR relaxation techniques: laboratory frame (15)N NMR relaxation, which is sensitive primarily to the ps-ns time scale, cross-correlated chemical shift modulation NMR relaxation detecting correlated mus-ms time scale motions of neighboring (13)C' and (15)N nuclei, and cross-correlated relaxation of two (15)N-(1)H dipole-dipole interactions detecting slow motions of backbone NH vectors in successive amino acid residues. The results indicate a reduction in mobility upon the mutation in several regions of the protein. The additional salt bridge formed in DHPS74E that links the N- and C-terminal subdomains is likely to be responsible for these changes.

A stabilizing alpha/beta-hydrophobic core greatly contributes to hyperthermostability of archaeal [P62A]Ssh10b

The hyperthermophilic Ssh10b from Sulfolobus shibatae is a member of the Sac10b family, which has been postulated to play a role in chromosomal organization in Archaea. Ssh10b is capable of significantly constraining negative DNA supercoils at elevated temperatures. In this study, the solution structure of the dimeric P62A mutant Ssh10b ([P62A]Ssh10b) was determined by multidimensional NMR spectroscopy. The backbone 15N dynamics, H/D exchange with and without the denaturant GdmSCN, and chemical and thermal denaturation experiments were performed to investigate the molecular basis of high thermostability of [P62A]Ssh10b. Data analysis has revealed an alpha/beta-hydrophobic core consisting of two alpha-helices and one beta-sheet which are stabilized by cooperative hydrophobic and hydrogen-bonding interactions. This stabilizing alpha/beta-hydrophobic core of [P62A]Ssh10b exhibiting highly restricted internal motions is composed of residues having highly protected amide protons which exchange with solvent mostly by means of a global unfolding process. The K40N mutation greatly destabilizes the mutant [P62A]Ssh10b because this mutation disturbs the packing of alpha-helix against the beta-sheet reducing the stability of the alpha/beta-hydrophobic core in the mutant protein. In comparison with homologous mesophilic and thermophilic proteins, it can be presumed that the stabilizing alpha/beta-hydrophobic core in the [P62A]Ssh10b structure greatly contributes to the high thermostability of the protein.

Folding kinetics and thermodynamics of Pseudomonas syringae effector protein AvrPto provide insight into translocation via the type III secretion system

In order to infect their hosts, many Gram-negative bacteria translocate agents of infection, called effector proteins, through the type III secretion system (TTSS) into the host cytoplasm. This process is thought to require at least partial unfolding of these agents, raising the question of how an effector protein might unfold to enable its translocation and then refold once it reaches the host cytoplasm. AvrPto is a well-studied effector protein of Pseudomonas syringae pv tomato. The presence of a readily observed unfolded population of AvrPto in aqueous solution and the lack of a known secretion chaperone make it ideal for studying the kinetic and thermodynamic characteristics that facilitate translocation. Application of Nzz exchange spectroscopy revealed a global, two-state folding equilibrium with 16% unfolded population, a folding rate of 1.8 s(-1), and an unfolding rate of 0.33 s(-1) at pH 6.1. TrAvrPto stability increases with increasing pH, with only 2% unfolded population observed at pH 7.0. The R(1) relaxation of TrAvrPto, which is sensitive to both the global anisotropy of folded TrAvrPto and slow exchange between folded and unfolded conformations, provided independent verification of the global kinetic rate constants. Given the acidic apoplast in which the pathogen resides and the more basic host cytoplasm, these results offer an intriguing mechanism by which the pH dependence of stability and slow folding kinetics of AvrPto would allow efficient translocation of the unfolded form through the TTSS and refolding into its functional folded form once inside the host.

Molecular basis of histone H3K4me3 recognition by ING4

The inhibitors of growth (ING) family of tumor suppressors consists of five homologous proteins involved in chromatin remodeling. They form part of different acetylation and deacetylation complexes and are thought to direct them to specific regions of the chromatin, through the recognition of H3K4me3 (trimethylated K4 in the histone 3 tail) by their conserved plant homeodomain (PHD). We have determined the crystal structure of ING4-PHD bound to H3K4me3, which reveals a tight complex stabilized by numerous interactions. NMR shows that there is a reduction in the backbone mobility on the regions of the PHD that participate in the peptide binding, and binding affinities differ depending on histone tail lengths Thermodynamic analysis reveals that the discrimination in favor of methylated lysine is entropy-driven, contrary to what has been described for chromodomains. The molecular basis of H3K4me3 recognition by ING4 differs from that of ING2, which is consistent with their different affinities for methylated histone tails. These differences suggest a distinct role in transcriptional regulation for these two ING family members because of the antagonistic effect of the complexes that they recruit onto chromatin. Our results illustrate the versatility of PHD fingers as readers of the histone code.

Backbone dynamics in an intramolecular prolylpeptide-SH3 complex from the diphtheria toxin repressor, DtxR

The diphtheria toxin repressor contains an SH3-like domain that forms an intramolecular complex with a proline-rich (Pr) peptide segment and stabilizes the inactive state of the repressor. Upon activation of diphtheria toxin repressor (DtxR) by transition metals, this intramolecular complex must dissociate as the SH3 domain and Pr segment form different interactions in the active repressor. Here we investigate the dynamics of this intramolecular complex using backbone amide nuclear spin relaxation rates determined using NMR spectroscopy and molecular dynamics trajectories. The SH3 domain in the unbound and bound states showed typical dynamics in that the secondary structures were fairly ordered with high generalized order parameters and low effective correlation times, while residues in the loops connecting beta-strands exhibited reduced generalized order parameters and required additional motional terms to adequately model the relaxation rates. Residues forming the Pr segment exhibited low-order parameters with internal rotational correlation times on the order of 0.6 ns-1 ns. Further analysis showed that the SH3 domain was rich in millisecond time scale motions while the Pr segment exhibited motions on the 100 mus time scale. Molecular dynamics simulations indicated structural rearrangements that may contribute to the observed relaxation rates and, together with the observed relaxation rate data, suggested that the Pr segment exhibits a binding<-->unbinding equilibrium. The results here provide new insights into the nature of the intramolecular complex and provide a better understanding of the biological role of the SH3 domain in regulating DtxR activity.

Protein dynamics from NMR: the slowly relaxing local structure analysis compared with model-free analysis

(15)N-(1)H spin relaxation is a powerful method for deriving information on protein dynamics. The traditional method of data analysis is model-free (MF), where the global and local N-H motions are independent and the local geometry is simplified. The common MF analysis consists of fitting single-field data. The results are typically field-dependent, and multifield data cannot be fit with standard fitting schemes. Cases where known functional dynamics has not been detected by MF were identified by us and others. Recently we applied to spin relaxation in proteins the slowly relaxing local structure (SRLS) approach, which accounts rigorously for mode mixing and general features of local geometry. SRLS was shown to yield MF in appropriate asymptotic limits. We found that the experimental spectral density corresponds quite well to the SRLS spectral density. The MF formulas are often used outside of their validity ranges, allowing small data sets to be force-fitted with good statistics but inaccurate best-fit parameters. This paper focuses on the mechanism of force-fitting and its implications. It is shown that MF analysis force-fits the experimental data because mode mixing, the rhombic symmetry of the local ordering and general features of local geometry are not accounted for. Combined multifield multitemperature data analyzed with the MF approach may lead to the detection of incorrect phenomena, and conformational entropy derived from MF order parameters may be highly inaccurate. On the other hand, fitting to more appropriate models can yield consistent physically insightful information. This requires that the complexity of the theoretical spectral densities matches the integrity of the experimental data. As shown herein, the SRLS spectral densities comply with this requirement.

Model-free analysis for large proteins at high magnetic field strengths

Protein backbone dynamics is often characterized using model-free analysis of three sets of (15)N relaxation data: longitudinal relaxation rate (R1), transverse relaxation rate (R2), and (15)N-{H} NOE values. Since the experimental data is limited, a simplified model-free spectral density function is often used that contains one Lorentzian describing overall rotational correlation but not one describing internal motion. The simplified spectral density function may be also used in estimating the overall rotational correlation time, by making the R2/R1 largely insensitive to internal motions, as well as used as one of the choices in the model selection protocol. However, such approximation may not be valid for analysis of relaxation data of large proteins recorded at high magnetic field strengths since the contribution to longitudinal relaxation from the Lorentzian describing the overall rotational diffusion of the molecule is comparably small relative to that describing internal motion. Here, we quantitatively estimate the errors introduced by the use of the simplified spectral density in model-free analysis for large proteins at high magnetic field strength.

Conformational dynamics in loop swap mutants of homologous fibronectin type III domains

Fibronectin type III (FN-III) domains are autonomously folded modules found in a variety of multidomain proteins. The 10th FN-III domain from fibronectin (fnFN10) and the 3rd FN-III domain from tenascin-C (tnFN3) have 27% sequence identity and the same overall fold; however, the CC' loop has a different pattern of backbone hydrogen bonds and the FG loop is longer in fnFN10 compared to tnFN3. To examine the influence of length, sequence, and context in determining dynamical properties of loops, CC' and FG loops were swapped between fnFN10 and tnFN3 to generate four mutant proteins and backbone conformational dynamics on ps-ns and mus-ms timescales were characterized by solution (15)N-NMR spin relaxation spectroscopy. The grafted loops do not strongly perturb the properties of the protein scaffold; however, specific effects of the mutations are observed for amino acids that are proximal in space to the sites of mutation. The amino acid sequence primarily dictates conformational dynamics when the wild-type and grafted loop have the same length, but both sequence and context contribute to conformational dynamics when the loop lengths differ. The results suggest that changes in conformational dynamics of mutant proteins must be considered in both theoretical studies and protein design efforts.

NMR investigation of Tyr105 mutants in TEM-1 beta-lactamase: dynamics are correlated with function

The existence of coupled residue motions on various time scales in enzymes is now well accepted, and their detailed characterization has become an essential element in understanding the role of dynamics in catalysis. To this day, a handful of enzyme systems has been shown to rely on essential residue motions for catalysis, but the generality of such phenomena remains to be elucidated. Using NMR spectroscopy, we investigated the electronic and dynamic effects of several mutations at position 105 in TEM-1 beta-lactamase, an enzyme responsible for antibiotic resistance. Even in absence of substrate, our results show that the number and magnitude of short and long range effects on (1)H-(15)N chemical shifts are correlated with the catalytic efficiencies of the various Y105X mutants investigated. In addition, (15)N relaxation experiments on mutant Y105D show that several active-site residues of TEM-1 display significantly altered motions on both picosecond-nanosecond and microsecond-millisecond time scales despite many being far away from the site of mutation. The altered motions among various active-site residues in mutant Y105D may account for the observed decrease in catalytic efficiency, therefore suggesting that short and long range residue motions could play an important catalytic role in TEM-1 beta-lactamase. These results support previous observations suggesting that internal motions play a role in promoting protein function.

Set theory formulation of the model-free problem and the diffusion seeded model-free paradigm

Model-free analysis of NMR relaxation data, which describes the motion of individual atoms, is a problem intricately linked to the Brownian rotational diffusion of the macromolecule. The diffusion tensor parameters strongly influence the optimisation of the various model-free models and the subsequent model selection between them. Finding the optimal model of the dynamics of the system among the numerous diffusion and model-free models is hence quite complex. Using set theory, the entirety of this global problem has been encapsulated by the universal set Ll, and its resolution mathematically formulated as the universal solution Ll. Ever since the original Lipari and Szabo papers the model-free dynamics of a molecule has most often been solved by initially estimating the diffusion tensor. The model-free models which depend on the diffusion parameter values are then optimised and the best model is chosen to represent the dynamics of the residue. Finally, the global model of all diffusion and model-free parameters is optimised. These steps are repeated until convergence. For simplicity this approach to Ll will be labelled the diffusion seeded model-free paradigm. Although this technique suffers from a number of problems many have been solved. All aspects of the diffusion seeded paradigm and its consequences, together with a few alternatives to the paradigm, will be reviewed through the use of set notation.

Dynamics of the SPRY domain-containing SOCS box protein 2: flexibility of key functional loops

The SPRY domain was identified originally as a sequence repeat in the dual-specificity kinase splA and ryanodine receptors and subsequently found in many other distinct proteins, including more than 70 encoded in the human genome. It is a subdomain of the B30.2/SPRY domain and is believed to function as a protein-protein interaction module. Three-dimensional structures of several B30.2/SPRY domain-containing proteins have been reported recently: murine SSB-2 in solution by NMR spectroscopy, a Drosophila SSB (GUSTAVUS), and human PRYSPRY protein by X-ray crystallography. The three structures share a core of two antiparallel beta-sheets for the B30.2/SPRY domain but show differences located mainly at one end of the beta-sandwich. Analysis of SSB-2 residues required for interactions with its intracellular ligands has provided insights into B30.2/SPRY binding specificity and identified loop residues critical for the function of this domain. We have investigated the backbone dynamics of SSB-2 by means of Modelfree analysis of its backbone (15)N relaxation parameters and carried out coarse-grained dynamics simulation of B30.2/SPRY domain-containing proteins using normal mode analysis. Translational self-diffusion coefficients of SSB-2 measured using pulsed field gradient NMR were used to confirm the monomeric state of SSB-2 in solution. These results, together with previously reported amide exchange data, highlight the underlying flexibility of the loop regions of B30.2/SPRY domain-containing proteins that have been shown to be important for protein-protein interactions. The underlying flexibility of certain regions of the B30.2/SPRY domain-containing proteins may also contribute to some apparent structural differences observed between GUSTAVUS or PRYSPRY and SSB-2.

The mechanism of discrimination between cognate and non-specific DNA by dimeric b/HLH/LZ transcription factors

The Myc/Max/Mad proteins are basic region-helix-loop-helix-leucine zipper (b/HLH/LZ) transcription factors that regulate the transcription of numerous genes involved in cell growth and proliferation. The Max protein is the obligate heterodimeric partner of the Myc and Mad proteins. Heterodimerization and DNA binding to target gene promoters are mediated by the b/HLH/LZ domains. Max can also form a homodimeric b/HLH/LZ. The enhanced expression of Myc and binding to promoters of target genes contribute to almost every aspect of tumor biology. However, the detailed mechanism by which dimeric and heterodimeric b/HLH/LZs discriminate cognate DNA (E-Box: CACGTG) from non-specific sequences in the target gene promoters is still unknown. Here, we use the Max b/HLH/LZ homodimer as a model for this class of transcription factors in the characterization and understanding of the mechanism of discrimination between the E-Box and non-specific DNA sequences. We report the characterization of a cognate and a non-specific Max b/HLH/LZ/DNA complex by EMSA, CD and NMR. Our results support a detailed mechanism by which dimeric b/HLH/LZ transcription factors can discriminate E-Box sequences from non-specific DNA. The mechanism proceeds via the conformational selection of fitting b/HLH/LZ homodimers with the basic region only partially helical. Next, the basic region undergoes a DNA-assisted folding or induced-fit. It is this step that provides the discrimination by stabilizing and destabilizing the alpha-helical conformation of the basic region in the cognate and non-specific complexes, respectively. This leads to a low affinity complex with a higher probability of being dissociated and hence to discrimination. A description of the side-chains and nucleotides proposed to be involved in the discrimination process is provided.

Methyl Dynamics in Proteins from NMR Slowly Relaxing Local Structure Spin Relaxation Analysis: A New Perspective

NMR spin relaxation of (2)H nuclei in (13)CH(2)D groups is a powerful method for studying side-chain motion in proteins. The analysis is typically carried out with the original model-free (MF) approach adapted to methyl dynamics. The latter is described in terms of axial local motions around, and of, the methyl averaging axis, mutually decoupled and independent of the global motion of the protein. Methyl motion is characterized primarily by the axial squared order parameter, S(axis)2, associated with fluctuations of the methyl averaging axis. This view is shown to be oversimplified by applying to typical experimental data the slowly relaxing local structure (SRLS) approach of Polimeno and Freed (Adv. Chem. Phys. 1993, 83, 89) which can be considered the generalization of the MF approach. Neglecting mode coupling and the asymmetry of the local ordering and treating approximately features of local geometry imply inaccurate values of S(axis)2, hence of the residual configurational entropy derived from it. S(axis)2, interpreted as amplitude of motion, was found to range from near disorder to almost complete order. Contrary to this picture, we find with the SRLS approach a moderate distribution in the magnitude of asymmetric local ordering and significant variation in its symmetry. The latter important property can be associated implicitly with the contribution of side-chain rotamer jumps. This is consistent with experimental residual dipolar coupling studies and theoretical work based on molecular dynamics simulations and molecular mechanics considerations. Configurational entropy is obtained in the SRLS approach directly from experimentally determined asymmetric potentials. Inconsistency between order parameters from 2H relaxation and from eta(HC-HH) cross-correlation and increase in order parameters with increasing temperature were observed with the MF approach. These discrepancies are reconciled, and physically tenable temperature dependence is obtained with the SRLS approach.

Backbone Dynamics of TEM-1 Determined by NMR: Evidence for a Highly Ordered Protein†

Backbone dynamics of TEM-1 beta-lactamase (263 amino acids, 28.9 kDa) were studied by 15N nuclear magnetic resonance relaxation at 11.7, 14.1, and 18.8 T. The high quality of the spectra allowed us to measure the longitudinal relaxation rate (R1), the transverse relaxation rate (R2), and the {1H}-15N NOE for up to 227 of the 250 potentially observable backbone amide groups. The model-free formalism was used to determine internal motional parameters using an axially anisotropic model. TEM-1 exhibits a small prolate axial anisotropy (D(parallel)/D(perpendicular) = 1.23 +/- 0.01) and a global correlation time (tau(m)) of 12.41 +/- 0.01 ns. The unusually high average generalized order parameter (S2) of 0.90 +/- 0.02 indicates that TEM-1 is one of the most ordered proteins studied by liquid-state NMR to date. Although the omega-loop has a high degree of order in the picosecond-to-nanosecond time scale (mean S2 value of 0.90 +/- 0.02), we observed the presence of microsecond-to-millisecond time scale motions for this loop, as for the vicinity of the active site. These motions could be relevant for the catalytic function of TEM-1. Amide exchange experiments were also performed, and several amide groups were not exchanged after 12 days, an indication that global motions in TEM-1 are also very limited. Although detailed dynamics characterization by NMR cannot be readily applied to TEM-1 in the presence of relevant substrates, the unusual picosecond-to-nanosecond dynamics behavior of TEM-1 presented here will be essential to the validation and improvement of future molecular dynamics simulations of TEM-1 in the presence of functionally relevant substrates.

Interaction of the intrinsically unstructured phage lambda N Protein with Escherichia coli NusA

N protein of the Escherichia coli phage lambda (lambdaN) is involved in antitermination, a transcription regulatory process that is essential for the expression of delayed early genes during phage lytic development. lambdaN is an intrinsically unstructured protein that possesses three distinct binding sites interacting with the carboxy terminus of the E. coli host factor protein NusA, the viral nutBoxB-RNA, and RNA polymerase, respectively. Heteronuclear NMR experiments with lambdaN(1-53) in complex with NusA(339-495) revealed that upon complex formation the lambdaN-binding interface, lambdaN(34-47), adopts a rigid structure. NMR data also indicate the induction of a weak helical structure in the nutboxB RNA-binding region lambdaN(1-22) upon binding to NusA(339-495) even in the absence of RNA. Titration experiments of the lambdaN(1-53)-nutBoxB RNA complex with NusA(339-495) revealed that the ternary complex can be described in terms of two structurally independent binary interactions. Furthermore, chemical-shift perturbation experiments with different NusA constructs and different lambdaN peptides showed that only NusA(353-416) is involved in lambdaN binding. We found that only one molecule of NusA(339-426) binds to one molecule of lambdaN(1-53). We also clarified the role of the lambdaN-binding region and could show that N41-R47 also binds to NusA(339-495). Furthermore, we observe that lambdaN(1-22) adopts a helical fold upon binding to NusA(339-495), in agreement with one of the theoretical models of lambdaN action.

Activation Energy of Catalysis-Related Domain Motion in E. coli Adenylate Kinase

Adenylate kinase from E. coli (AKeco), folded into domains CORE, AMPbd, and LID, catalyzes the reaction AMP + ATP <--> 2ADP. Previous X-ray crystallography and optical solution methods showed that the domains AMPbd and LID, and the conserved P-loop, execute large-amplitude catalysis-related motions. We used (15)N NMR spin relaxation methods to find that the simplified model-free (MF) analysis does not, whereas our general Slowly Relaxing Local Structure analysis does, detect catalytic domain motion. SRLS set for the first time the correlation time for domain motion at tau(L)perpendicular = 8.2 ns, to be compared with tau(m) = 15.1 ns for global tumbling. These results were obtained at 303 K. Herein we conduct a temperature-dependent investigation of tau(L)perpendicular and tau(m) in the range of 288-310 K. We found that the activation energy for global tumbling is Ea = 16.9 +/- 0.5 kJ/mol, the hydrodynamic volume of hydrated AKeco is 65.6 +/- 2.1 nm3, its radius is 2.50 +/- 0.03 nm, and the number of hydration layers is 1.77. The average tau(L)perpendicular value decreases from 11 ns at 288 K to 4 ns at 310 K, with activation energies of 29.7 +/- 3.3, 32.1 +/- 4.3, and 30.4 +/- 4.3 kJ/mol for the domains AMPbd and LID, and the catalytic P-loop, respectively. These values are two-to-three times smaller than typical activation energies of enzymatic reactions. Hence kinase catalysis appears not to be controlled by domain motion in the ligand-free enzyme. However, the latter process clearly facilitates important mechanical aspects such as steric recognition and capturing of the AMP and ATP substrates, their proper positioning for phosphorylation, and the release of the ADP product.

NMR dynamic studies suggest that allosteric activation regulates ligand binding in chicken liver bile acid-binding protein

Apo chicken liver bile acid-binding protein has been structurally characterized by NMR. The dynamic behavior of the protein in its apo- and holo-forms, complexed with chenodeoxycholate, has been determined via (15)N relaxation and steady state heteronuclear (15)N((1)H) nuclear Overhauser effect measurements. The dynamic parameters were obtained at two pH values (5.6 and 7.0) for the apoprotein and at pH 7.0 for the holoprotein, using the model free approach. Relaxation studies, performed at three different magnetic fields, revealed a substantial conformational flexibility on the microsecond to millisecond time scales, mainly localized in the C-terminal face of the beta-barrel. The observed dynamics are primarily caused by the protonation/deprotonation of a buried histidine residue, His(98), located on this flexible face. A network of polar buried side chains, defining a spine going from the E to J strand, is likely to provide the long range connectivity needed to communicate motion from His(98) to the EF loop region. NMR data are accompanied by molecular dynamics simulations, suggesting that His(98) protonation equilibrium is the triggering event for the modulation of a functionally important motion, i.e. the opening/closing at the protein open end, whereas ligand binding stabilizes one of the preexisting conformations (the open form). The results presented here, complemented with an analysis of proteins belonging to the intracellular lipid-binding protein family, are consistent with a model of allosteric activation governing the binding mechanism. The functional role of this mechanism is thoroughly discussed within the framework of the mechanism for the enterohepatic circulation of bile acids.

PhosphoThr peptide binding globally rigidifies much of the FHA domain from Arabidopsis receptor kinase-associated protein phosphatase

A net increase in the backbone rigidity of the kinase-interacting FHA domain (KI-FHA) from the Arabidopsis receptor kinase-associated protein phosphatase (KAPP) accompanies the binding of a phosphoThr peptide from its CLV1 receptor-like kinase partner, according to (15)N NMR relaxation at 11.7 and 14.1 T. All of the loops of free KI-FHA display evidence of nanosecond-scale motions. Many of these same residues have residual dipolar couplings that deviate from structural predictions. Binding of the CLV1 pT868 peptide seems to reduce nanosecond-scale fluctuations of all loops, including half of the residues of recognition loops. Residues important for affinity are found to be rigid, i.e., conserved residues and residues of the subsite for the key pT+3 peptide position. This behavior parallels SH2 and PTB domain recognition of pTyr peptides. PhosphoThr peptide binding increases KI-FHA backbone rigidity (S(2)) of three recognition loops, a loop nearby, seven strands from the beta-sandwich, and a distal loop. Compensating the trend of increased rigidity, binding enhances fast mobility at a few sites in four loops on the periphery of the recognition surface and in two loops on the far side of the beta-sandwich. Line broadening evidence of microsecond- to millisecond-scale fluctuations occurs across the six-stranded beta-sheet and nearby edges of the beta-sandwich; this forms a network connected by packing of interior side chains and H-bonding. A patch of the slowly fluctuating residues coincides with the site of segment-swapped dimerization in crystals of the FHA domain of human Chfr. Phosphopeptide binding introduces microsecond- to millisecond-scale fluctuations to more residues of the long 8/9 recognition loop of KI-FHA. The rigidity of this FHA domain appears to couple as a whole to pThr peptide binding.