Lipari-Szabo mapping: A graphical approach to Lipari-Szabo analysis of NMR relaxation data using reduced spectral density mapping

Impact of the Hereditary P301L Mutation on the Correlated Conformational Dynamics of Human Tau Protein Revealed by the Paramagnetic Relaxation Enhancement NMR Experiments

Tau forms intracellular insoluble aggregates as a neuropathological hallmark of Alzheimer’s disease. Tau is largely unstructured, which complicates the characterization of the tau aggregation process. Recent studies have demonstrated that tau samples two distinct conformational ensembles, each of which contains the soluble and aggregation-prone states of tau. A shift to populate the aggregation-prone ensemble may promote tau fibrillization. However, the mechanism of this ensemble transition remains elusive. In this study, we explored the conformational dynamics of a tau fragment by using paramagnetic relaxation enhancement (PRE) and interference (PRI) NMR experiments. The PRE correlation map showed that tau is composed of segments consisting of residues in correlated motions. Intriguingly, residues forming the β-structures in the heparin-induced tau filament coincide with residues in these segments, suggesting that each segment behaves as a structural unit in fibrillization. PRI data demonstrated that the P301L mutation exclusively alters the transiently formed tau structures by changing the short- and long-range correlated motions among residues. The transient conformations of P301L tau expose the amyloid motif PHF6 to promote tau self-aggregation. We propose the correlated motions among residues within tau determine the population sizes of the conformational ensembles, and perturbing the correlated motions populates the aggregation-prone form.

Dynamics of the intrinsically disordered inhibitor IF7 of glutamine synthetase in isolation and in complex with its partner

Glutamine synthetase (GS) catalyzes the ATP-dependent formation of glutamine from glutamate and ammonia. The activity of Synechocystis sp. PCC 6803 GS is regulated, among other mechanisms, by protein-protein interactions with a 65-residue-long, intrinsically disordered protein (IDP), named IF7. IDPs explore diverse conformations in their free states and, in some cases, in their molecular complexes. We used both nuclear magnetic resonance (NMR) at 11.7 T and small angle X-ray scattering (SAXS) to study the size and the dynamics in the picoseconds-to-nanosecond (ps-ns) timescale of: (i) isolated IF7; and (ii) the IF7/GS complex. Our SAXS findings, together with MD results, show: (i) some of the possible IF7 structures in solution; and, (ii) that the presence of IF7 affected the structure of GS in solution. The joint use of SAXS and NMR shows that movements of each amino acid of IF7 were uncorrelated with those of its neighbors. Residues of IF7 with the largest values of the relaxation rates (R₁, R₂ and η_xy), in the free and bound species, were mainly clustered around: (i) the C terminus of the protein; and (ii) Ala30. These residues, together with Arg8 (which is a hot-spot residue in the interaction with GS), had a restricted mobility in the presence of GS. The C-terminal region, which appeared more compact in our MD simulations of isolated IF7, seemed to be involved in non-native contacts with GS that help in the binding between the two macromolecules.

Dynamics of the intrinsically disordered protein NUPR1 in isolation and in its fuzzy complexes with DNA and prothymosin α

Intrinsically disordered proteins (IDPs) explore diverse conformations in their free states and, a few of them, also in their molecular complexes. This functional plasticity is essential for the function of IDPs, although their dynamics in both free and bound states is poorly understood. NUPR1 is a protumoral multifunctional IDP, activated during the acute phases of pancreatitis. It interacts with DNA and other IDPs, such as prothymosin α (ProTα), with dissociation constants of ~0.5 μM, and a 1:1 stoichiometry. We studied the structure and picosecond-to-nanosecond (ps-ns) dynamics by using both NMR and SAXS in: (i) isolated NUPR1; (ii) the NUPR1/ProTα complex; and (iii) the NUPR1/double stranded (ds) GGGCGCGCCC complex. Our SAXS findings show that NUPR1 remained disordered when bound to either partner, adopting a worm-like conformation; the fuzziness of bound NUPR1 was also pinpointed by NMR. Residues with the largest values of the relaxation rates (R₁, R_1ρ, R₂ and η_xy), in the free and bound species, were mainly clustered around the 30s region of the sequence, which agree with one of the protein hot-spots already identified by site-directed mutagenesis. Not only residues in this region had larger relaxation rates, but they also moved slower than the rest of the molecule, as indicated by the reduced spectral density approach (RSDA). Upon binding, the energy landscape of NUPR1 was not funneled down to a specific, well-folded conformation, but rather its backbone flexibility was kept, with distinct motions occurring at the hot-spot region.

Continuous Interdomain Orientation Distributions Reveal Components of Binding Thermodynamics

The flexibility of biological macromolecules is an important structural determinant of function. Unfortunately, the correlations between different motional modes are poorly captured by discrete ensemble representations. Here, we present new ways to both represent and visualize correlated interdomain motions. Interdomain motions are determined directly from residual dipolar couplings, represented as a continuous conformational distribution, and visualized using the disk-on-sphere representation. Using the disk-on-sphere representation, features of interdomain motions, including correlations, are intuitively visualized. The representation works especially well for multidomain systems with broad conformational distributions.This analysis also can be extended to multiple probability density modes, using a Bingham mixture model. We use this new paradigm to study the interdomain motions of staphylococcal protein A, which is a key virulence factor contributing to the pathogenicity of Staphylococcus aureus. We capture the smooth transitions between important states and demonstrate the utility of continuous distribution functions for computing the reorientational components of binding thermodynamics. Such insights allow for the dissection of the dynamic structural components of functionally important intermolecular interactions.

Allosteric Breakage of the Hydrogen Bond within the Dual-Histidine Motif in the Active Site of Human Pin1 PPIase

Intimate cooperativity among active site residues in enzymes is a key factor for regulating elaborate reactions that would otherwise not occur readily. Peptidyl-prolyl cis-trans isomerase NIMA-interacting 1 (Pin1) is the phosphorylation-dependent cis-trans peptidyl-prolyl isomerase (PPIase) that specifically targets phosphorylated Ser/Thr-Pro motifs. Residues C113, H59, H157, and T152 form a hydrogen bond network in the active site, as in the noted connection. Theoretical studies have shown that protonation to thiolate C113 leads to rearrangement of this hydrogen bond network, with switching of the tautomeric states of adjacent histidines (H59 and H157) [Barman, A., and Hamelberg, D. (2014) Biochemistry 53, 3839-3850]. This is called the "dual-histidine motif". Here, C113A and C113S Pin1 mutants were found to alter the protonation states of H59 according to the respective residue type replaced at C113, and the mutations resulted in disruption of the hydrogen bond within the dual-histidine motif. In the C113A mutant, H59 was observed to be in exchange between ε- and δ-tautomers, which widened the entrance of the active site cavity, as seen by an increase in the distance between residues A113 and S154. The C113S mutant caused H59 to exchange between the ε-tautomer and imidazolium while not changing the active site structure. Moreover, the imidazole ring orientations of H59 and H157 were changed in the C113S mutant. These results demonstrated that a mutation at C113 modulates the hydrogen bond network dynamics. Thus, C113 acts as a pivot to drive the concerted function among the residues in the hydrogen bond network, as theoretically predicted.

NMR relaxation and structural elucidation of peptides in the presence and absence of trifluoroethanol illuminates the critical molecular nature of integrin αvβ6 ligand specificity

Integrin αvβ6 is an important emerging target for both imaging and therapy of cancer that requires specific ligands based on Arg-Gly-Asp (RGD) peptides. There remains little correlation between integrin-RGD ligand specificity despite studies suggesting an RGD-turn-helix ligand motif is required. Here, we describe the application of 15N NMR relaxation analyses and structure determination of αvβ6 peptide ligands in the presence and absence of trifluoroethanol (TFE) to identify their critical molecular nature that influences specificity, interaction and function. Two linear peptides; one known to demonstrate αvβ6 specificity (FMDV2) and the other based on a natural RGD ligand (LAP2), were compared to two additional peptides based on FMDV2 but cyclised in different positions using a disulphide bond (DBD1 and DBD2). The cyclic adaptation in DBD1 produces a significant alteration in backbone dynamic properties when compared to FMDV2; a potential driver for the loss in αvβ6 specificity by DBD1. The importance of ligand dynamics are highlighted through a comprehensive reduced spectral density and ModelFree analysis of peptide 15N NMR relaxation data and suggest αvβ6 specificity requires the formation of a structurally rigid helix preceded by a RGD motif exhibiting slow internal motion. Additional observations include the effect of TFE/water viscosity on global NMR dynamics and the advantages of using spectral density NMR relaxation data to estimate correlation times and motional time regimes for peptides in solution.

Influence of heme post-translational modification and distal ligation on the backbone dynamics of a monomeric hemoglobin

The cyanobacterium Synechococcus sp. PCC 7002 uses a hemoglobin of the truncated lineage (GlbN) in the detoxification of reactive species generated in the assimilation of nitrate. In view of a sensing or enzymatic role, several states of GlbN are of interest with respect to its structure-activity relationship. Nuclear magnetic resonance spectroscopy was applied to compare the structure and backbone dynamics of six GlbN forms differing in their oxidation state [Fe(II) or Fe(III)], distal ligand to the iron (histidine, carbon monoxide, or cyanide), or heme post-translational modification (b heme or covalently attached heme). Structural properties were assessed with pseudocontact shift calculations. (15)N relaxation data were analyzed by reduced spectral density mapping (picosecond to nanosecond motions) and by inspection of elevated R(2) values (microsecond to millisecond motions). On the picosecond to nanosecond time scale, GlbN exhibited little flexibility and was unresponsive to the differences among the various forms. Regions of slightly higher mobility were the CE turn, the EF loop, and the H-H' kink. In contrast, fluctuations on the microsecond to millisecond time scale depended on the form. Cyanide binding to the ferric state did not enhance motions, whereas reduction to the ferrous bis-histidine state resulted in elevated R(2) values for several amides. This response was attributed, at least in part, to a weakening of the distal histidine coordination. Carbon monoxide binding quenched some of these fluctuations. The results emphasized the role of the distal ligand in dictating backbone flexibility and illustrated the multiple ways in which motions are controlled by the hemoglobin fold.

Dynamics of backbone conformational heterogeneity in Bacillus subtilis ribonuclease P protein

Interconversion of protein conformations is imperative to function, as evidenced by conformational changes associated with enzyme catalytic cycles, ligand binding and post-translational modifications. In this study, we used 15N NMR relaxation experiments to probe the fast (i.e., ps-ns) and slow (i.e., micros-ms) conformational dynamics of Bacillus subtilis ribonuclease P protein (P protein) in its folded state, bound to two sulfate anions. Using the Lipari-Szabo mapping method [Andrec, M., Montelione, G. T., and Levy, R. M. (2000) J. Biomol. NMR 18, 83-100] to interpret the data, we find evidence for P protein dynamics on the mus-ms time scale in the ensemble. The residues that exhibit these slow internal motions are found in regions that have been previously identified as part of the P protein-P RNA interface. These results suggest that structural flexibility within the P protein ensemble may be important for proper RNase P holoenzyme assembly and/or catalysis.

Backbone dynamics of the monomeric lambda repressor denatured state ensemble under nondenaturing conditions

Oxidizing two native methionine residues predominantly populates the denatured state of monomeric lambda repressor (MetO-lambdaLS) under nondenaturing conditions. NMR was used to characterize the secondary structure and dynamics of MetO-lambdaLS in standard phosphate buffer. 13Calpha and 1Halpha chemical shift indices reveal a region of significant helicity between residues 9 and 29. This helical content is further supported by the observation of medium-range amide NOEs. The remaining residues do not exhibit significant helicity as determined by NMR. We determined 15N relaxation parameters for 64 of 85 residues at 600 and 800 MHz. There are two distinct regions of reduced flexibility, residues 8-32 in the N-terminal third and residues 50-83 in the C-terminal third. The middle third, residues 33-50, has greater flexibility. We have analyzed the amplitude of the backbone motions in terms of the physical properties of the amino acids and conclude that conformational restriction of the backbone MetO-lambdaLS is due to nascent helix formation in the region corresponding to native helix 1. The bulkiness of amino acid residues in the C-terminal third leads to the potential for hydrophobic interactions, which is suggested by chemical exchange detected by the difference in spectral density J(0) at the two static magnetic fields. The more flexible middle region is the result of a predominance of small side chains in this region.

Retracted article: In-cell protein dynamics

The native intracellular environment of proteins is crowded with metabolites and macromolecules. However, most biophysical information concerning proteins is acquired in dilute solution. To determine whether there are differences in dynamics, nuclear magnetic resonance spectroscopy can be used to measure 15N relaxation in uniformly 15N-enriched apocytochrome b5 inside living Escherichia coli and in dilute solution. Such data can then be used to compare the fast backbone dynamics of the partially folded protein in cells to its dynamics in dilute solution by using Lipari-Szabo analysis. It appears that the intracellular environment does not alter the protein's structure, or significantly change its fast dynamics. Specifically, the cytosol does not change the amplitude of fast backbone motions, but does increase the average timescale of these motions, most likely due to the increase in viscosity of the cytosol.

Solution Structure and Backbone Dynamics of the KH-QUA2 Region of the Xenopus STAR/GSG Quaking Protein

The Quaking protein belongs to the family of STAR/GSG domain RNA-binding proteins and is involved in multiple cell signalling and developmental processes in vertebrates, including the formation of myelin. Heteronuclear NMR methods were used to determine the solution structure of a 134 residue fragment spanning the KH and QUA2 homology regions of the Quaking protein from Xenopus laevis (pXqua) in the absence of RNA. The protein is shown to adopt an extended type I KH domain fold that is connected to a structured alpha-helix in the C-terminal QUA2 region by means of a highly flexible linker. A comparison with the solution structure of the related protein splicing factor 1 (SF1) indicates that most aspects of the RNA-binding interface are conserved in pXqua, although the "variable loop" region that follows the second beta-strand possesses two additional alpha-helices. The structure of pXqua provides an appropriate template for building models of important homologues, such as GLD-1 and Sam68. Measurements of the (15)N relaxation parameters of pXqua confirm that the polypeptide backbone of the QUA2 region is more dynamic than that of the KH portion, and that the C-terminal helix is partially structured in the absence of RNA. By comparison with a random coil reference state, the nascent structure in the QUA2 region is estimated to contribute 15.5kJmol(-1) to the change in conformational free energy that occurs on forming a complex with RNA. Since STAR/GSG proteins may regulate alternative splicing by competing with SF1 in the nucleus for specific branch-point sequences that signal intronic RNA, the formation of secondary structure in the QUA2 region in the unbound state of pXqua has important functional consequences.

Identification of a Residue Critical for Maintaining the Functional Conformation of BPTI

The effects of amino acid replacements on the backbone dynamics of bovine pancreatic trypsin inhibitor (BPTI) were examined using 15N NMR relaxation experiments. Previous studies have shown that backbone amide groups within the trypsin-binding region of the wild-type protein undergo conformational exchange processes on the micros time scale, and that replacement of Tyr35 with Gly greatly increases the number of backbone atoms involved in such motions. In order to determine whether these mutational effects are specific to the replacement of this residue with Gly, six additional replacements were examined in the present study. In two of these, Tyr35 was replaced with either Ala or Leu, and the other four were single replacements of Tyr23, Phe33, Asn43 or Asn44, all of which are highly buried in the native structure and conserved in homologous proteins. The Y35A and Y35L mutants displayed dynamic properties very similar to those of the Y35G mutant, with the backbone segments including residues 10-19 and 32-44 undergoing motions revealed by enhanced 15N transverse relaxation rates. On the other hand, the Y23L, N43G and N44A substitutions caused almost no detectable changes in backbone dynamics, on either the ns-ps or ms-micros time scales, even though each of these replacements significantly destabilizes the native conformation. Replacement of Phe33 with Leu caused intermediate effects, with several residues that have previously been implicated in motions in the wild-type protein displaying enhanced transverse relaxation rates. These results demonstrate that destabilizing amino acid replacements can be accommodated in a native protein with dramatically different effects on conformational dynamics and that Tyr35 plays a particularly important role in defining the conformation of the trypsin-binding site of BPTI.

Effects of temperature on the dynamic behaviour of the HIV-1 nucleocapsid NCp7 and its DNA complex

The nucleocapsid protein NCp7 of human immunodeficiency virus type 1 (HIV-1) contains two highly conserved CCHC zinc fingers and is involved in many crucial steps of the virus life-cycle. A large number of physiological rôles of NCp7 involve its binding to single-stranded nucleic acid chains. Several solution structures of NCp7 and its complex with single-stranded RNA or DNA have been reported. We have investigated the changes in the dynamic behaviour experienced by the (12-53)NCp7 peptide upon DNA binding using (15)N heteronuclear relaxation measurements at 293 K and 308 K, and fluorescence spectroscopy. The relaxation data were interpreted using the reduced spectral density approach, which allowed the high-frequency motion, overall tumbling rates and the conformational exchange contributions to be characterized for various states of the peptide without using a specific motional model. Analysis of the temperature-dependent correlation times derived from both NMR and fluorescence data indicated a co-operative change of the molecular shape of apo (12-53)NCp7 around 303 K, leading to an increased hydrodynamic radius at higher temperatures. The binding of (12-53)NCp7 to a single-stranded d(ACGCC) pentanucleotide DNA led to a reduction of the conformational flexibility that characterized the apo peptide. Translational diffusion experiments as well as rotational correlation times indicated that the (12-53)NCp7/d(ACGCC) complex tumbles as a rigid object. The amplitudes of high-frequency motions were restrained in the complex and the occurrence of conformational exchange was displaced from the second zinc finger to the linker residue Ala30.