Ultrahigh-resolution backbone structure of perdeuterated protein GB1 using residual dipolar couplings from two alignment media.
Explore the source record for details and available documents.
Bacterial Proteins↗
Biomedical subjects
Martin Blackledge
Publications and source records attributed to Martin Blackledge.
Simultaneous determination of protein backbone structure and dynamics from residual dipolar couplings.
Determination of protein structure classically results in a single average configuration that takes no account of conformational fluctuation. Dynamics are, however, inherently linked to structure and crucial to our understanding of biological function. In this study we have used analytical descriptions of dynamic averaging of residual dipolar couplings (RDCs) to simultaneously determine the backbone structure and dynamics of protein GB3. RDCs alone are used to determine an ultrahigh-resolution structure that compares very closely with a refined X-ray structure (rmsd of 0.34 A overall backbone residues). Dynamic amplitudes reporting on motions up to the millisecond time scale reproduce the main characteristics of dynamics previously determined in conjunction with the crystal structure. The use of RDCs alone allows a bias-free comparison with a purely static approach to structure determination. Extensive cross validation clearly demonstrates that the dynamic description is superior to the static approximation. The demonstration that this level of structural resolution and dynamic detail can be extracted from RDCs supports previous indications that these parameters contain extremely precise information about biomolecular conformational sampling.
Models, Molecular↗
Charge-induced molecular alignment of intrinsically disordered proteins.
Explore the source record for details and available documents.
Amino Acids↗
Residual dipolar couplings in short peptides reveal systematic conformational preferences of individual amino acids.
Residual dipolar couplings (RDCs) observed by NMR in solution under weak alignment conditions can monitor average net orientations and order parameters of individual bonds. By their simple geometrical dependence, RDCs bear particular promise for the quantitative characterization of conformations in partially folded or unfolded proteins. We have systematically investigated the influence of amino acid substitutions X on the conformation of unfolded model peptides EGAAXAASS as monitored by their (1)H(Nu)-(15)N and (1)H(alpha)-(13)C(alpha) RDCs detected at natural abundance of (15)N and (13)C in strained polyacrylamide gels. In total, 14 single amino acid substitutions were investigated. The RDCs show a specific dependence on the substitution X that correlates to steric or hydrophobic interactions with adjacent amino acids. In particular, the RDCs for the glycine and proline substitutions indicate less or more order, respectively, than the other amino acids. The RDCs for aromatic substitutions tryptophane and tyrosine give evidence of a kink in the peptide backbone. This effect is also observable for orientation by Pf1 phages and corroborated by variations in (13)C(alpha) secondary shifts and (3)J(HNH)(alpha) scalar couplings in isotropic samples. RDCs for a substitution with the beta-turn sequence KNGE differ from single amino acid substitutions. Terminal effects and next neighbor effects could be demonstrated by further specific substitutions. The results were compared to statistical models of unfolded peptide conformations derived from PDB coil subsets, which reproduce overall trends for (1)H(Nu)-(15)N RDCs for most substitutions, but deviate more strongly for (1)H(alpha)-(13)C(alpha) RDCs. The outlined approach opens the possibility to obtain a systematic experimental characterization of the influence of individual amino acid/amino acid interactions on orientational preferences in polypeptides.
Amino Acid Sequence↗
The influence of nitrogen-15 proton-driven spin diffusion on the measurement of nitrogen-15 longitudinal relaxation times.
The effect of nitrogen-15 proton-driven spin diffusion on quantitative (15)N T(1) measurements in solid proteins is investigated, and the impact on the measurement of dynamic parameters is assessed. A simple model of exchange between neighboring nitrogens is used to reproduce the evolution of (15)N spin systems whose longitudinal relaxation rates and exchange rates are compatible with experimental measurements. We show that the induced error in the measured T(1) and its effect on the determination of dynamics parameters is likely to be less than the current experimental error. The use of deuterated protein samples is shown to have a small but sometimes visible effect, and may also considerably slow down or even suppress the exchange of magnetization due to spin diffusion.
Artifacts↗
Sequence-specific solvent accessibilities of protein residues in unfolded protein ensembles.
Protein stability cannot be understood without the correct description of the unfolded state. We present here an efficient method for accurate calculation of atomic solvent exposures for denatured protein ensembles. The method used to generate the ensembles has been shown to reproduce diverse biophysical experimental data corresponding to natively and chemically unfolded proteins. Using a data set of 19 nonhomologous proteins containing from 98 to 579 residues, we report average accessibilities for all residue types. These averaged accessibilities are considerably lower than those previously reported for tripeptides and close to the lower limit reported by Creamer and co-workers. Of importance, we observe remarkable sequence dependence for the exposure to solvent of all residue types, which indicates that average residue solvent exposures can be inappropriate to interpret mutational studies. In addition, we observe smaller influences of both protein size and protein amino acid composition in the averaged residue solvent exposures for individual proteins. Calculating residue-specific solvent accessibilities within the context of real sequences is thus necessary and feasible. The approach presented here may allow a more precise parameterization of protein energetics as a function of polar- and apolar-area burial and opens new ways to investigate the energetics of the unfolded state of proteins.
Computer Simulation↗
Observation of heteronuclear overhauser effects confirms the 15N-1H dipolar relaxation mechanism in a crystalline protein.
The observation of proton to nitrogen-15 heteronuclear Overhauser effects in the microcrystalline protein Crh is used to confirm that the principal mechanism of relaxation of amide nitrogens is due to the fluctuation of the N-H dipolar couplings caused by N-H bond dynamics. Our observations reveal the central role of water as the main source of proton magnetization, and we provide an analysis of the different pathways that could lead to the observed results.
Bacterial Proteins↗
Effect of the solvent on the conformation of a depsipeptide: NMR-derived solution structure of hormaomycin in DMSO from residual dipolar couplings in a novel DMSO-compatible alignment medium.
The macrocyclic compound hormaomycin has been investigated by NMR spectroscopy and by restrained molecular-dynamics simulations. Measurement of residual dipolar couplings induced by dissolving the depsipeptide in a polyacrylamide gel compatible with DMSO and their incorporation into the structure calculation of the title compound improved the precision of the family of structures. In DMSO the macrocyclic ring shows two beta-turns, whose positions in the sequence differ from those found in the CDCl3 solution structure and in the crystal structure obtained from hexylene glycol/H2O (50:50). The bulky side chain consisting of a 3-(2-nitrocyclopropyl)alanine and a chlorinated N-hydroxypyrrole moiety is flexible in DMSO.
Bacterial Proteins↗
Quantitative analysis of backbone dynamics in a crystalline protein from nitrogen-15 spin-lattice relaxation.
A detailed analysis of nitrogen-15 longitudinal relaxation times in microcrystalline proteins is presented. A theoretical model to quantitatively interpret relaxation times is developed in terms of motional amplitude and characteristic time scale. Different averaging schemes are examined in order to propose an analysis of relaxation curves that takes into account the specificity of MAS experiments. In particular, it is shown that magic angle spinning averages the relaxation rate experienced by a single spin over one rotor period, resulting in individual relaxation curves that are dependent on the orientation of their corresponding carousel with respect to the rotor axis. Powder averaging thus leads to a nonexponential behavior in the observed decay curves. We extract dynamic information from experimental decay curves, using a diffusion in a cone model. We apply this study to the analysis of spin-lattice relaxation rates of the microcrystalline protein Crh at two different fields and determine differential dynamic parameters for several residues in the protein.
Bacterial Proteins↗
Defining long-range order and local disorder in native alpha-synuclein using residual dipolar couplings.
Intrinsically unstructured proteins play key biochemical roles in a vast range of normal and pathological processes. To study these systems, it is necessary to invoke an ensemble of rapidly interconverting conformations. Residual dipolar couplings (RDCs) are particularly powerful probes of the behavior of unfolded proteins, reporting on time and ensemble-averaged conformations up to and beyond the millisecond time scale. In this study, we present a novel interpretation of RDCs in unfolded systems that simultaneously defines long-range structural order and local conformational sampling. This approach is used to describe the structure and dynamics of alpha-Synuclein (alphaS), a protein that is strongly implicated in the development of Parkinson's disease (PD), allowing unambiguous detection of strongly populated conformers containing long-range contacts between the N- and C-terminal domains. The structural model combines two features required for the description of alphaS in solution: local conformational fluctuation based on random sampling of residue-specific phi/psi distributions, and long-range contacts induced by the presence of nonbonding interactions between domains that are distant in primary sequence. Both aspects are found to be necessary for the reproduction of the nonaveraged RDCs from alphaS. Although RDCs have previously been shown to report on local conformational preferences in unstructured proteins, this study shows the additional sensitivity of these measurements to the presence of long-range order in highly flexible systems.
Protein Conformation↗
A structural model for unfolded proteins from residual dipolar couplings and small-angle x-ray scattering.
Natively unfolded proteins play key roles in normal and pathological biochemical processes. Despite their importance for function, this category of proteins remains beyond the reach of classical structural biology because of their inherent conformational heterogeneity. We present a description of the intrinsic conformational sampling of unfolded proteins based on residue-specific /Psi propensities from loop regions of a folded protein database and simple volume exclusion. This approach is used to propose a structural model of the 57-aa, natively disordered region of the nucleocapsid-binding domain of Sendai virus phosphoprotein. Structural ensembles obeying these simple rules of conformational sampling are used to simulate averaged residual dipolar couplings (RDCs) and small-angle x-ray scattering data. This protein is particularly informative because RDC data from the equally sized folded and unfolded domains both report on the unstructured region, allowing a quantitative analysis of the degree of order present in this part of the protein. Close agreement between experimental and simulated RDC and small-angle x-ray scattering data validates this simple model of conformational sampling, providing a precise description of local structure and dynamics and average dimensions of the ensemble of sampled structures. RDC data from two urea-unfolded systems are also closely reproduced. The demonstration that conformational behavior of unfolded proteins can be accurately predicted from the primary sequence by using a simple set of rules has important consequences for our understanding of the structure and dynamics of the unstructured state.
Models, Molecular↗
Identification of slow correlated motions in proteins using residual dipolar and hydrogen-bond scalar couplings.
Despite their importance for biological activity, slower molecular motions beyond the nanosecond range remain poorly understood. We have assembled an unprecedented set of experimental NMR data, comprising up to 27 residual dipolar couplings per amino acid, to define the nature and amplitude of backbone motion in protein G using the Gaussian axial fluctuation model in three dimensions. Slower motions occur in the loops, and in the beta-sheet, and are absent in other regions of the molecule, including the alpha-helix. In the beta-sheet an alternating pattern of dynamics along the peptide sequence is found to form a long-range network of slow motion in the form of a standing wave extending across the beta-sheet, resulting in maximal conformational sampling at the interaction site. The alternating nodes along the sequence match the alternation of strongly hydrophobic side chains buried in the protein core. Confirmation of the motion is provided through extensive cross-validation and by independent hydrogen-bond scalar coupling analysis that shows this motion to be correlated. These observations strongly suggest that dynamical information can be transmitted across hydrogen bonds and have important implications for understanding collective motions and long-range information transfer in proteins.
Amino Acids↗
Solution structure of the sulfite reductase flavodoxin-like domain from Escherichia coli.
The flavoprotein moiety of Escherichia coli sulfite reductase (SiR-FP) is homologous to electron transfer proteins such as cytochrome-P450 reductase (CPR) or nitric oxide synthase (NOS). We report on the three-dimensional structure of SiR-FP18, the flavodoxin-like domain of SiR-FP, which has been determined by NMR. In the holoenzyme, this domain plays an important role by shuttling electrons from the FAD to the hemoprotein (the beta-subunit). The structure presented here was determined using distance and torsion angle information in combination with residual dipolar couplings determined in two different alignment media. Several protein-FMN NOEs allowed us to place the prosthetic group in its binding pocket. The structure is well-resolved, and (15)N relaxation data indicate that SiR-FP18 is a compact domain. The binding interface with cytochrome c, a nonphysiological electron acceptor, has been determined using chemical shift mapping. Comparison of the SiR-FP18 structure with the corresponding domains from CPR and NOS shows that the fold of the protein core is highly conserved, but the analysis of the electrostatic surfaces reveals significant differences between the three domains. These observations are placed in the physiological context so they can contribute to the understanding of the electron transfer mechanism in the SiR holoenzyme.
Binding Sites↗
Protein backbone dynamics from N-HN dipolar couplings in partially aligned systems: a comparison of motional models in the presence of structural noise.
Residual dipolar couplings (RDCs) provide excellent probes for the exploration of dynamics in biomolecules on biologically relevant time-scales. Applying geometric motional models in combination with high-resolution structures to fit experimental RDCs allows the extraction of local dynamic amplitudes of peptide planes in proteins using only a limited number of data points. Here we compare the behaviour of three simple and intuitive dynamic modes: the Gaussian axial fluctuation model (1D-GAF), the two-site jump model, and a model supposing axially symmetric motion about a mean orientation. The requirement of a structural model makes this kind of methodology potentially very sensitive to structural imprecision. Numerical simulations of RDC dynamic averaging under different regimes show that the anisotropic motional models are more geometrically stringent than the axially symmetric model making it more difficult to alias structural noise as artificial dynamic amplitudes. Indeed, it appears that the model extracts accurate motional amplitudes even in the presence of significant structural error. We also show that a two-site jump model, also assuming the (alpha)C(i-1)-(alpha)C(i) as rotation axis, can only be distinguished from the previously developed GAF model beyond amplitude/jumps of around 40 degrees. The importance of appropriate estimation of the molecular alignment tensor for determination of local motional amplitudes is also illustrated here. We demonstrate a systematic scaling of extracted dynamic amplitudes if a static structure is assumed when determining the alignment tensor from dynamically averaged RDCs. As an example an artificial increase of 0.14 (0.85 compared to the expected 0.71) is observed in the extracted S2 if a pervasive 20 degrees GAF motion is present that is ignored in the tensor determination. Finally we apply a combined approach using the most appropriate motional model, to complete the analysis of dynamic motions from protein G.
Amino Acids↗
Site-specific backbone dynamics from a crystalline protein by solid-state NMR spectroscopy.
Site-specific nitrogen-15 longitudinal relaxation rates are measured for the microcrystalline dimeric form of the protein Crh using multidimensional high-resolution solid-state NMR methods. The measured rates are used to provide a qualitative description of the site-specific internal mobility of the protein present in the solid state.
Bacterial Proteins↗
Local dynamic amplitudes on the protein backbone from dipolar couplings: toward the elucidation of slower motions in biomolecules.
Local protein backbone dynamics on time scales reaching up to milliseconds have been investigated using residual dipolar couplings (RDC) using an analytical description of conformational averaging under the influence of anisotropic peptide plane dynamics. We have applied this technique to RDC from protein G and find that sites in the alpha-helix exhibit overall higher-order parameters than loops, suggesting a high degree of conformational integrity even over this extended time period. The approach is shown to be stable when using data from a smaller number of alignment media. Order parameters derived from combinations of independent subsets of two and three of the five alignment media from protein G reveal results essentially identical to those from the complete data set. Structures of lysozyme determined at different crystal diffraction resolutions ranging from 0.9 to 2.1 A give similar dynamic parameters using this method, demonstrating robustness with respect to structural error.
Models, Molecular↗
Anisotropic small amplitude Peptide plane dynamics in proteins from residual dipolar couplings.
The effect of small amplitude anisotropic peptide plane motion on residual dipolar couplings (RDC) measured in proteins has been investigated. RDC averaging effects in the presence of GAF (Gaussian axial fluctuation) motions are found to vary strongly depending on the peptide plane orientation. Even low amplitude dynamics can significantly affect derived alignment tensor parameters if this motion is not taken into account. An analytical description of averaged N-(N)H RDCs is introduced that includes basic GAF-like motion. The averaging depends on the orientation of the peptide plane (alpha, beta, gamma) in the alignment frame and on the motional amplitude (sigma). This expression is used to investigate the presence of anisotropic reorientational dynamics in proteins by incorporating sigma as an additional parameter into the alignment tensor analysis. Average GAF amplitudes (sigma(av)) are determined for secondary structural elements from single experimental N-(N)H RDC data sets from five different proteins, in combination with high-resolution structural models. This yields statistically significant improvement over the static description, and detects sigma(av) values ranging from 14.4 to 17.0 degrees for the different proteins. A higher value of sigma(av) = 20 degrees from loop regions was found using two independent sets of N-(N)H RDC in the protein lysozyme, for which a very high-resolution structure is available. Comparison of fitting behavior over 13 structures from lysozyme of crystal diffraction resolution ranging from 0.9 to 2.1A indicates a small spread of motional amplitudes, demonstrating that the method is robust up to this level of resolution. A combined definition of (alpha)C-C' and N-(N)H RDC under the influence of GAF motions allows simultaneous fitting of both RDC. Application to three proteins leads to similar sigma(av) values and a more significant improvement with respect to the static model. Using the GAF model to describe conformationally averaged RDC is important for two reasons: a more accurate definition of the alignment tensor magnitude can be derived, and the method can be used to detect average small amplitude motions in protein backbones from readily accessible data, on time scales not easily sampled by other NMR techniques.
Anisotropy↗
Structure and dynamics of the nucleocapsid-binding domain of the Sendai virus phosphoprotein in solution.
The RNA-dependent RNA polymerase of the Sendai virus (SeV) consists of the large protein (L) and the phosphoprotein (P). P plays a crucial role in the enzyme by positioning L (which carries the polymerase activity) onto the matrix for transcription and replication formed by the RNA and the nucleoprotein, the N-RNA. P has a modular structure with distinct functional domains: an N-terminal domain involved in binding to N degrees (N that is not yet bound to RNA) and a C-terminal domain that carries the oligomerisation domain, the N-RNA binding domain and the L binding domain and that, combined with L, is active in transcription. Structural data have previously been obtained on the N-terminal domain and on the oligomerisation domain of P, but not yet on its N-RNA binding domain (also-called the X protein). Here we present an NMR and a small angle neutron scattering study of the SeV X protein. We show that this molecule presents two subdomains linked by an 11-residue linker, with the N-subdomain lacking a well-defined conformation. The 3D structure of the C-subdomain consists of three alpha-helices revealing an asymmetric charge distribution that may be important for binding to RNA-bound nucleoprotein. The structure of the entire C-terminal domain of P is modelled from its constituent parts in combination with small angle scattering data on this domain.
Amino Acid Sequence↗