Comparison of the backbone dynamics of a folded and an unfolded SH3 domain existing in equilibrium in aqueous buffer

SARS-CoV-2 Nsp8 N-terminal domain folds autonomously and binds dsRNA

The SARS-CoV-2 Nsp8 protein is a critical component of the RNA replicase, as its N-terminal domain (NTD) anchors Nsp12, the RNA, and Nsp13. Whereas its C-terminal domain (CTD) structure is well resolved, there is an open debate regarding the conformation adopted by the NTD as it is predicted as disordered but found in a variety of complex-dependent conformations or missing from many other structures. Using NMR spectroscopy, we show that the SARS CoV-2 Nsp8 NTD features both well folded secondary structure and disordered segments. Our results suggest that while part of this domain corresponding to two long α-helices forms autonomously, the folding of other segments would require interaction with other replicase components. When isolated, the α-helix population progressively declines towards the C-termini but surprisingly binds dsRNA while preserving structural disorder.

An autoinhibited state of 53BP1 revealed by small molecule antagonists and protein engineering

The recruitment of 53BP1 to chromatin, mediated by its recognition of histone H4 dimethylated at lysine 20 (H4K20me2), is important for DNA double-strand break repair. Using a series of small molecule antagonists, we demonstrate a conformational equilibrium between an open and a pre-existing lowly populated closed state of 53BP1 in which the H4K20me2 binding surface is buried at the interface between two interacting 53BP1 molecules. In cells, these antagonists inhibit the chromatin recruitment of wild type 53BP1, but do not affect 53BP1 variants unable to access the closed conformation despite preservation of the H4K20me2 binding site. Thus, this inhibition operates by shifting the conformational equilibrium toward the closed state. Our work therefore identifies an auto-associated form of 53BP1-autoinhibited for chromatin binding-that can be stabilized by small molecule ligands encapsulated between two 53BP1 protomers. Such ligands are valuable research tools to study the function of 53BP1 and have the potential to facilitate the development of new drugs for cancer therapy.

Studying micro to millisecond protein dynamics using simple amide 15N CEST experiments supplemented with major-state R2 and visible peak-position constraints

Over the last decade amide 15N CEST experiments have emerged as a popular tool to study protein dynamics that involves exchange between a 'visible' major state and sparsely populated 'invisible' minor states. Although initially introduced to study exchange between states that are in slow exchange with each other (typical exchange rates of, 10 to 400 s-1), they are now used to study interconversion between states on the intermediate to fast exchange timescale while still using low to moderate (5 to 350 Hz) 'saturating' B1 fields. The 15N CEST experiment is very sensitive to exchange as the exchange delay TEX can be quite long (~0.5 s) allowing for a large number of exchange events to occur making it a very powerful tool to detect minor sates populated ([Formula: see text]) to as low as 1%. When systems are in fast exchange and the 15N CEST data has to be described using a model that contains exchange, the exchange parameters are often poorly defined because the [Formula: see text] versus [Formula: see text] and [Formula: see text] versus exchange rate ([Formula: see text]) plots can be quite flat with shallow or no minima and the analysis of such 15N CEST data can lead to wrong estimates of the exchange parameters due to the presence of 'spurious' minima. Here we show that the inclusion of experimentally derived constraints on the intrinsic transverse relaxation rates and the inclusion of visible state peak-positions during the analysis of amide 15N CEST data acquired with moderate B1 values (~50 to ~350 Hz) results in convincing minima in the [Formula: see text] versus [Formula: see text] and the [Formula: see text] versus [Formula: see text] plots even when exchange occurs on the 100 μs timescale. The utility of this strategy is demonstrated on the fast-folding Bacillus stearothermophilus peripheral subunit binding domain that folds with a rate constant ~104 s-1. Here the analysis of 15N CEST data alone results in [Formula: see text] versus [Formula: see text] and [Formula: see text] versus [Formula: see text] plots that contain shallow minima, but the inclusion of visible-state peak positions and restraints on the intrinsic transverse relaxation rates of both states during the analysis of the 15N CEST data results in pronounced minima in the [Formula: see text] versus [Formula: see text] and [Formula: see text] versus [Formula: see text] plots and precise exchange parameters even in the fast exchange regime ([Formula: see text]~5). Using this strategy we find that the folding rate constant of PSBD is invariant (~10,500 s-1) from 33.2 to 42.9 °C while the unfolding rates (~70 to ~500 s-1) and unfolded state populations (~0.7 to ~4.3%) increase with temperature. The results presented here show that protein dynamics occurring on the 10 to 104 s-1 timescale can be studied using amide 15N CEST experiments.

An autoinhibited state of 53BP1 revealed by small molecule antagonists and protein engineering

The recruitment of 53BP1 to chromatin, mediated by its recognition of histone H4 dimethylated at lysine 20 (H4K20me2), is important for DNA double-strand break repair. Using a series of small molecule antagonists, we demonstrate a conformational equilibrium between an open and a pre-existing lowly populated closed state of 53BP1 in which the H4K20me2 binding surface is buried at the interface between two interacting 53BP1 molecules. In cells, these antagonists inhibit the chromatin recruitment of wild type 53BP1, but do not affect 53BP1 variants unable to access the closed conformation despite preservation of the H4K20me2 binding site. Thus, this inhibition operates by shifting the conformational equilibrium toward the closed state. Our work therefore identifies an auto-associated form of 53BP1 - autoinhibited for chromatin binding - that can be stabilized by small molecule ligands encapsulated between two 53BP1 protomers. Such ligands are valuable research tools to study the function of 53BP1 and have the potential to facilitate the development of new drugs for cancer therapy.

Mechanistic insights into the aggregation pathway of the patient-derived immunoglobulin light chain variable domain protein FOR005

Systemic antibody light chain (AL) amyloidosis is characterized by deposition of amyloid fibrils. Prior to fibril formation, soluble oligomeric AL protein has a direct cytotoxic effect on cardiomyocytes. We focus on the patient derived λ-III AL variable domain FOR005 which is mutated at five positions with respect to the closest germline protein. Using solution-state NMR spectroscopy, we follow the individual steps involved in protein misfolding from the native to the amyloid fibril state. Unfavorable mutations in the complementary determining regions introduce a strain in the native protein structure which yields partial unfolding. Driven by electrostatic interactions, the protein converts into a high molecular weight, oligomeric, molten globule. The high local concentration of aggregation prone regions in the oligomer finally catalyzes the conversion into fibrils. The topology is determined by balanced electrostatic interactions in the fibril core implying a 180° rotational switch of the beta-sheets around the conserved disulfide bond.

Retrospective study for the universal applicability of the residue-based linear free energy relationship in the two-state exchange of protein molecules

Multiprobe measurements, such as NMR and hydrogen exchange studies, can provide the equilibrium constant, K, and rate constants for forward and backward processes, k and k′, of the two-state structural changes of a polypeptide on a per-residue basis. We previously found a linear relationship between log K and log k and between log K and log k′ for the topological exchange of a 27-residue bioactive peptide. To test the general applicability of the residue-based linear free energy relationship (rbLEFR), we performed a literature search to collect residue-specific K, k, and k′ values in various exchange processes, including folding-unfolding equilibrium, coupled folding and binding of intrinsically disordered peptides, and structural fluctuations of folded proteins. The good linearity in a substantial number of the log–log plots proved that the rbLFER holds for the structural changes in a wide variety of protein-related phenomena. Among the successful cases, the hydrogen exchange study of apomyoglobin folding intermediates is particularly interesting. We found that the residues that deviated from the linear relationship corresponded to the α-helix, for which transient translocation had been identified by other experiments. Thus, the rbLFER is useful for studying the structures and energetics of the dynamic states of protein molecules.

Explicit models of motions to analyze NMR relaxation data in proteins

Nuclear Magnetic Resonance (NMR) is a tool of choice to characterize molecular motions. In biological macromolecules, pico- to nano-second motions, in particular, can be probed by nuclear spin relaxation rates which depend on the time fluctuations of the orientations of spin interaction frames. For the past 40 years, relaxation rates have been successfully analyzed using the Model Free (MF) approach which makes no assumption on the nature of motions and reports on the effective amplitude and time-scale of the motions. However, obtaining a mechanistic picture of motions from this type of analysis is difficult at best, unless complemented with molecular dynamics (MD) simulations. In spite of their limited accuracy, such simulations can be used to obtain the information necessary to build explicit models of motions designed to analyze NMR relaxation data. Here, we present how to build such models, suited in particular to describe motions of methyl-bearing protein side-chains and compare them with the MF approach. We show on synthetic data that explicit models of motions are more robust in the presence of rotamer jumps which dominate the relaxation in methyl groups of protein side-chains. We expect this work to motivate the use of explicit models of motion to analyze MD and NMR data.

The time-zero HSQC method improves the linear free energy relationship of a polypeptide chain through the accurate measurement of residue-specific equilibrium constants

EXSY (exchange spectroscopy) NMR provides the residue-specific equilibrium constants, K, and residue-specific kinetic rate constants, k, of a polypeptide chain in a two-state exchange in the slow exchange regime. A linear free energy relationship (LFER) discovered in a log k versus log K plot is considered to be a physicochemical basis for smooth folding and conformational changes of protein molecules. For accurate determination of the thermodynamic and kinetic parameters, the measurement bias arising from state-specific differences in the R₁ and R₂ relaxation rates of ¹H and other nuclei in HSQC and EXSY experiments must be minimized. Here, we showed that the time-zero HSQC acquisition scheme (HSQC0) is effective for this purpose, in combination with a special analytical method (Π analysis) for EXSY. As an example, we applied the HSQC0 + Π method to the two-state exchange of nukacin ISK-1 in an aqueous solution. Nukacin ISK-1 is a 27-residue lantibiotic peptide containing three mono-sulfide linkages. The resultant bias-free residue-based LFER provided valuable insights into the transition state of the topological interconversion of nukacin ISK-1. We found that two amino acid residues were exceptions in the residue-based LFER relationship. We inferred that the two residues could adopt special conformations in the transition state, to allow the threading of some side chains through a ring structure formed by one of the mono-sulfide linkages. In this context, the two residues are a useful target for the manipulation of the physicochemical properties and biological activities of nukacin ISK-1.

Structural basis of protein substrate processing by human mitochondrial high-temperature requirement A2 protease

Protein aggregates are often toxic, leading to impaired cellular activities and disease. The human HtrA2 trimeric enzyme cleaves such aggregates, and mutations in HtrA2 are causative for various neurodegenerative disorders, such as Parkinson’s disease and essential tremor. The mechanism by which cleavage occurs has been studied using small peptides, but little information is available as to how HtrA2 protects cells from the pathologic effects of aggregation involving protein molecules that can form well-folded structures. Using solution NMR spectroscopy, we investigated the structural dynamics of the interaction between HtrA2 and a model protein substrate, demonstrating that HtrA2 preferentially binds to an unfolded substrate ensemble and providing insights into how HtrA2 function is regulated.

The human high-temperature requirement A2 (HtrA2) protein is a trimeric protease that cleaves misfolded proteins to protect cells from stresses caused by toxic, proteinaceous aggregates, and the aberrant function of HtrA2 is closely related to the onset of neurodegenerative disorders. Our methyl-transverse relaxation optimized spectroscopy (TROSY)–based NMR studies using small-peptide ligands have previously revealed a stepwise activation mechanism involving multiple distinct conformational states. However, very little is known about how HtrA2 binds to protein substrates and if the distinct conformational states observed in previous peptide studies might be involved in the processing of protein clients. Herein, we use solution-based NMR spectroscopy to investigate the interaction between the N-terminal Src homology 3 domain from downstream of receptor kinase (drk) with an added C-terminal HtrA2-binding motif (drkN SH3-PDZbm) that exhibits marginal folding stability and serves as a mimic of a physiological protein substrate. We show that drkN SH3-PDZbm binds to HtrA2 via a two-pronged interaction, involving both its C-terminal PDZ-domain binding motif and a central hydrophobic region, with binding occurring preferentially via an unfolded ensemble of substrate molecules. Multivalent interactions between several clients and a single HtrA2 trimer significantly stimulate the catalytic activity of HtrA2, suggesting that binding avidity plays an important role in regulating substrate processing. Our results provide a thermodynamic, kinetic, and structural description of the interaction of HtrA2 with protein substrates and highlight the importance of a trimeric architecture for function as a stress-protective protease that mitigates aggregation.

NMR Provides Unique Insight into the Functional Dynamics and Interactions of Intrinsically Disordered Proteins

Intrinsically disordered proteins are ubiquitous throughout all known proteomes, playing essential roles in all aspects of cellular and extracellular biochemistry. To understand their function, it is necessary to determine their structural and dynamic behavior and to describe the physical chemistry of their interaction trajectories. Nuclear magnetic resonance is perfectly adapted to this task, providing ensemble averaged structural and dynamic parameters that report on each assigned resonance in the molecule, unveiling otherwise inaccessible insight into the reaction kinetics and thermodynamics that are essential for function. In this review, we describe recent applications of NMR-based approaches to understanding the conformational energy landscape, the nature and time scales of local and long-range dynamics and how they depend on the environment, even in the cell. Finally, we illustrate the ability of NMR to uncover the mechanistic basis of functional disordered molecular assemblies that are important for human health.

Protein Dynamics to Define and Refine Disordered Protein Ensembles

Intrinsically disordered proteins and unfolded proteins have fluctuating conformational ensembles that are fundamental to their biological function and impact protein folding, stability, and misfolding. Despite the importance of protein dynamics and conformational sampling, time-dependent data types are not fully exploited when defining and refining disordered protein ensembles. Here we introduce a computational framework using an elastic network model and normal-mode displacements to generate a dynamic disordered ensemble consistent with NMR-derived dynamics parameters, including transverse R₂ relaxation rates and Lipari-Szabo order parameters (S² values). We illustrate our approach using the unfolded state of the drkN SH3 domain to show that the dynamical ensembles give better agreement than a static ensemble for a wide range of experimental validation data including NMR chemical shifts, J-couplings, nuclear Overhauser effects, paramagnetic relaxation enhancements, residual dipolar couplings, hydrodynamic radii, single-molecule fluorescence Förster resonance energy transfer, and small-angle X-ray scattering.

Molecular Dynamics-Assisted Optimization of Protein NMR Relaxation Analysis

NMR relaxation analysis of the mobile residues in globular proteins is sensitive to the form of the experimentally fitted internal autocorrelation function, which is used to represent that motion. Different order parameter representations can precisely fit the same set of ¹⁵N R₁, R₂, and heteronuclear NOE measurements while yielding significantly divergent predictions of the underlying autocorrelation functions, indicating the insufficiency of these experimental relaxation data for assessing which order parameter representation provides the most physically realistic predictions. Molecular dynamics simulations offer an unparalleled capability for discriminating among different order parameter representations to assess which representation can most accurately model a wide range of physically realistic autocorrelation functions. Six currently utilized AMBER and CHARMM force fields were applied to calculate autocorrelation functions for the backbone H–N bond vectors of ubiquitin as an operational test set. An optimized time constant-constrained triexponential (TCCT) representation was shown to markedly outperform the widely used (S_f²,τ_s,S²) extended Lipari–Szabo representation and the more closely related (S_f²,S_H², S_N²) Larmor frequency-selective representation. Optimization of the TCCT representation at both 600 and 900 MHz ¹H converged to the same parameterization. The higher magnetic field yielded systematically larger deviations in the back-prediction of the autocorrelation functions for the mobile amides, indicating little added benefit from multiple field measurements in analyzing amides that lack slower (∼ms) exchange line-broadening effects. Experimental ¹⁵N relaxation data efficiently distinguished among the different force fields with regard to their prediction of ubiquitin backbone conformational dynamics in the ps–ns time frame. While the earlier AMBER 99SB and CHARMM27 force fields underestimate the scale of backbone dynamics, which occur in this time frame, AMBER 14SB provided the most consistent predictions for the well-averaged highly mobile C-terminal residues of ubiquitin.

Structure, dynamics, and function of SrnR, a transcription factor for nickel-dependent gene expression

Streptomyces griseus, a bacterium producing antibacterial drugs and featuring possible application in phytoremediation, expresses two metal-dependent superoxide dismutase (SOD) enzymes, containing either Fe(II) or Ni(II) in their active site. In particular, the alternative expression of the two proteins occurs in a metal-dependent mode, with the Fe(II)-enzyme gene (sodF) repressed at high intracellular Ni(II) concentrations by a two-component system (TCS). This complex involves two proteins, namely SgSrnR and SgSrnQ, which represent the transcriptional regulator and the Ni(II) sensor of the system, respectively. SgSrnR belongs to the ArsR/SmtB family of metal-dependent transcription factors; in the apo-form and in the absence of SgSrnQ, it can bind the DNA operator of sodF, upregulating gene transcription. According to a recently proposed hypothesis, Ni(II) binding to SgSrnQ would promote its interaction with SgSrnR, causing the release of the complex from DNA and the consequent downregulation of the sodF expression. SgSrnQ is predicted to be highly disordered, thus the understanding, at the molecular level, of how the SgSrnR/SgSrnQ TCS specifically responds to Ni(II) requires the knowledge of the structural, dynamic, and functional features of SgSrnR. These were investigated synergistically in this work using X-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy, atomistic molecular dynamics calculations, isothermal titration calorimetry, and in silico molecular docking. The results reveal that the homodimeric apo-SgSrnR binds to its operator in a two-step process that involves the more rigid globular portion of the protein and leaves its largely disordered regions available to possibly interact with the disordered SgSrnQ in a Ni-dependent process.

Fluorine NMR study of proline-rich sequences using fluoroprolines

Proline homopolymer motifs are found in many proteins; their peculiar conformational and dynamic properties are often directly involved in those proteins' functions. However, the dynamics of proline homopolymers is hard to study by NMR due to a lack of amide protons and small chemical shift dispersion. Exploiting the spectroscopic properties of fluorinated prolines opens interesting perspectives to address these issues. Fluorinated prolines are already widely used in protein structure engineering - they introduce conformational and dynamical biases - but their use as 19 F NMR reporters of proline conformation has not yet been explored. In this work, we look at model peptides where Cγ -fluorinated prolines with opposite configurations of the chiral Cγ centre have been introduced at two positions in distinct polyproline segments. By looking at the effects of swapping these (4R )-fluoroproline and (4S )-fluoroproline within the polyproline segments, we were able to separate the intrinsic conformational properties of the polyproline sequence from the conformational alterations instilled by fluorination. We assess the fluoroproline 19 F relaxation properties, and we exploit the latter in elucidating binding kinetics to the SH3 (Src homology 3) domain.

Residue-Specific Kinetic Insights into the Transition State in Slow Polypeptide Topological Isomerization by NMR Exchange Spectroscopy

The characterization of the transition state is a central issue in biophysical studies of protein folding. NMR is a multiprobe measurement technique that provides residue-specific information. Here, we used exchange spectroscopy to characterize the transition state of the two-state slow topological isomerization of a 27-residue lantibiotic peptide. The exchange kinetic rates varied on a per-residue basis, indicating the reduced kinetic cooperativity of the two-state exchange, as well as the previously observed reduced thermodynamic cooperativity. Furthermore, temperature-dependent measurements revealed large variations in the activation enthalpy and entropy terms among residues. Interestingly, we found a linear relationship between the logarithm of the equilibrium constants and that of the exchange rates. Because the data points are derived from amino acid residues in one polypeptide chain, we refer to the linear relationship as the residue-based linear free energy relationship (rbLFER). The rbLFER offers information about the transition state of the two-state exchange.

Rapid measurement of heteronuclear transverse relaxation rates using non-uniformly sampled R1ρ accordion experiments

Multidimensional, heteronuclear NMR relaxation methods are used extensively to characterize the dynamics of biological macromolecules. Acquisition of relaxation datasets on proteins typically requires significant measurement time, often several days. Accordion spectroscopy offers a powerful means to shorten relaxation rate measurements by encoding the "relaxation dimension" into the indirect evolution period in multidimensional experiments. Time savings can also be achieved by non-uniform sampling (NUS) of multidimensional NMR data, which is used increasingly to improve spectral resolution or increase sensitivity per unit time. However, NUS is not commonly implemented in relaxation experiments, because most reconstruction algorithms are inherently nonlinear, leading to problems when estimating signal intensities, relaxation rate constants and their error bounds. We have previously shown how to avoid these shortcomings by combining accordion spectroscopy with NUS, followed by data reconstruction using sparse exponential mode analysis, thereby achieving a dramatic decrease in the total length of longitudinal relaxation experiments. Here, we present the corresponding transverse relaxation experiment, taking into account the special considerations required for its successful implementation in the framework of the accordion-NUS approach. We attain the highest possible precision in the relaxation rate constants by optimizing the NUS scheme with respect to the Cramér-Rao lower bound of the variance of the estimated parameter, given the total number of sampling points and the spectrum-specific signal characteristics. The resulting accordion-NUS R 1 ρ relaxation experiment achieves comparable precision in the parameter estimates compared to conventional CPMG (Carr-Purcell-Meiboom-Gill) R 2 or spin-lock R 1 ρ experiments while saving an order of magnitude in experiment time.

Unveiling protein dynamics in solution with field-cycling NMR relaxometry

Field-cycling NMR relaxometry is a well-established technique that can give information on molecular structure and dynamics of biological systems. It provides the nuclear relaxation rates as a function of the applied magnetic field, starting from fields as low as ~ 10^-4 T up to about 1-3 T. The profiles so collected, called nuclear magnetic relaxation dispersion (NMRD) profiles, can be extended to include the relaxation rates at the largest fields achievable with high resolution NMR spectrometers. By exploiting this wide range of frequencies, the NMRD profiles can provide information on motions occurring on time scales from 10^-6 to 10^-9 s. ¹H NMRD measurements have proved very useful also for the characterization of paramagnetic proteins, because they can help characterise a number of parameters including the number, distance and residence time of water molecules coordinated to the paramagnetic center, the reorientation correlation times and the electron spin relaxation time, and the electronic structure at the metal site.

Danio rerio Oocytes for Eukaryotic In-Cell NMR

Understanding how the crowded and complex cellular milieu affects protein stability and dynamics has only recently become possible by using techniques such as in-cell nuclear magnetic resonance. However, the combination of stabilizing and destabilizing interactions makes simple predictions difficult. Here we show the potential of Danio rerio oocytes as an in-cell nuclear magnetic resonance model that can be widely used to measure protein stability and dynamics. We demonstrate that in eukaryotic oocytes, which are 3-6-fold less crowded than other cell types, attractive chemical interactions still dominate effects on protein stability and slow tumbling times, compared to the effects of dilute buffer.

Connecting the conformational behavior of cyclic octadepsipeptides with their ionophoric property and membrane permeability

Cyclic octadepsipeptides such as PF1022A and its synthetic derivative emodepside exhibit anthelmintic activity with the latter sold as a commercial drug treatment against gastrointestinal nematodes for animal health use. The structure-permeability relationship of these cyclic depsipeptides that could ultimately provide insights into the compound bioavailability is not yet well understood. The fully N-methylated amide backbone and apolar sidechain residues do not allow for the formation of intramolecular hydrogen bonds, normally observed in the membrane-permeable conformations of cyclic peptides. Hence, any understanding gained on these depsipeptides would serve as a prototype for future design strategies. In previous nuclear magnetic resonance (NMR) studies, two macrocyclic core conformers of emodepside were detected, one with all backbone amides in trans-configuration (hereon referred as the symmetric conformer) and the other with one amide in cis-configuration (hereon referred as the asymmetric conformer). In addition, these depsipeptides were also reported to be ionophores with a preference of potassium over sodium. In this study, we relate the conformational behavior of PF1022A, emodepside, and closely related analogs with their ionophoric characteristic probed using NMR and molecular dynamics (MD) simulations and finally evaluated their passive membrane permeability using PAMPA. We find that the equilibrium between the two core conformers shifts more towards the symmetric conformer upon addition of monovalent cations with selectivity for potassium over sodium. Both the NMR experiments and the theoretical Markov state models based on extensive MD simulations indicate a more rigid backbone for the asymmetric conformation, whereas the symmetric conformation shows greater flexibility. The experimental results further advocate for the symmetric conformation binding the cation. The PAMPA results suggest that the investigated depsipeptides are retained in the membrane, which may be advantageous for the likely target, a membrane-bound potassium channel.

Assessing Site-Specific Enhancements Imparted by Hyperpolarized Water in Folded and Unfolded Proteins by 2D HMQC NMR

Hyperpolarized water can be a valuable aid in protein NMR, leading to amide group ¹H polarizations that are orders of magnitude larger than their thermal counterparts. Suitable procedures can exploit this to deliver 2D ¹H–¹⁵N correlations with good resolution and enhanced sensitivity. These enhancements depend on the exchange rates between the amides and the water, thereby yielding diagnostic information about solvent accessibility. This study applied this “HyperW” method to four proteins exhibiting a gamut of exchange behaviors: PhoA^(350–471), an unfolded 122-residue fragment; barstar, a fully folded ribonuclease inhibitor; R17, a 13.3 kDa system possessing folded and unfolded forms under slow interconversion; and drkN SH3, a protein domain whose folded and unfolded forms interchange rapidly and with temperature-dependent population ratios. For PhoA4^(350–471) HyperW sensitivity enhancements were ≥300×, as expected for an unfolded protein sequence. Though fully folded, barstar also exhibited substantial enhancements; these, however, were not uniform and, according to CLEANEX experiments, reflected the solvent-exposed residues. R17 showed the expected superposition of ≥100-fold enhancements for its unfolded form, coexisting with more modest enhancements for their folded counterparts. Unexpected, however, was the behavior of drkN SH3, for which HyperW enhanced the unfolded but, surprisingly, enhanced even more certain folded protein sites. These preferential enhancements were repeatedly and reproducibly observed. A number of explanations—including three-site exchange magnetization transfers between water and the unfolded and folded states; cross-correlated relaxation processes from hyperpolarized “structural” waters and labile side-chain protons; and the possibility that faster solvent exchange rates characterize certain folded sites over their unfolded counterparts—are considered to account for them.

13C NMR Relaxation Analysis of Protein GB3 for the Assessment of Side Chain Dynamics Predictions by Current AMBER and CHARMM Force Fields

Molecular simulations with seven current AMBER- and CHARMM-based force fields yield markedly differing internal bond vector autocorrelation function predictions for many of the 223 methine and methylene H-C bonds of the 56-residue protein GB3. To enable quantification of accuracy, ¹³C R₁, R₂, and heteronuclear NOE relaxation rates have been determined for the methine and stereochemically assigned methylene C^α and C^β positions. With only three experimental relaxation values for each bond vector, central to this analysis is the accuracy with which MD-derived autocorrelation curves can be represented by a 3-parameter equation which, in turn, maps onto the NMR relaxation values. In contrast to the more widely used extended Lipari-Szabo order parameter representation, 95% of these MD-derived internal autocorrelation curves for GB3 can be fitted to within 1.0% rmsd over the time frame from 30 ps to 4 ns by a biexponential Larmor frequency-selective representation (LF-S²). Applying the LF-S² representation to the experimental relaxation rates and uncertainties serves to determine the boundary range for the autocorrelation function of each bond vector consistent with the experimental data. Not surprisingly, all seven force fields predict the autocorrelation functions for the more motionally restricted ¹H^α-¹³C^α and ¹H^β-¹³C^β bond vectors with reasonable accuracy. However, for the ¹H^β-¹³C^β bond vectors exhibiting aggregate order parameter S² values less than 0.85, only 1% of the MD-derived predictions lie with 1 σ of the experimentally determined autocorrelation functions and only 7% within 2 σ. On the other hand, substantial residue type-specific improvements in predictive performance were observed among the recent AMBER force fields. This analysis indicates considerable potential for the use of ¹³C relaxation measurements in guiding the optimization of the side chain dynamics characteristics of protein molecular simulations.

An autoinhibitory intramolecular interaction proof-reads RNA recognition by the essential splicing factor U2AF2

The recognition of cis-regulatory RNA motifs in human transcripts by RNA binding proteins (RBPs) is essential for gene regulation. The molecular features that determine RBP specificity are often poorly understood. Here, we combined NMR structural biology with high-throughput iCLIP approaches to identify a regulatory mechanism for U2AF2 RNA recognition. We found that the intrinsically disordered linker region connecting the two RNA recognition motif (RRM) domains of U2AF2 mediates autoinhibitory intramolecular interactions to reduce nonproductive binding to weak Py-tract RNAs. This proofreading favors binding of U2AF2 at stronger Py-tracts, as required to define 3' splice sites at early stages of spliceosome assembly. Mutations that impair the linker autoinhibition enhance the affinity for weak Py-tracts result in promiscuous binding of U2AF2 along mRNAs and impact on splicing fidelity. Our findings highlight an important role of intrinsically disordered linkers to modulate RNA interactions of multidomain RBPs.

Intrinsic disorder and amino acid specificity modulate binding of the WW2 domain in kidney and brain protein (KIBRA) to synaptopodin

The second WW domain (WW2) of the kidney and brain scaffolding protein, KIBRA, has an isoleucine (Ile-81) rather than a second conserved tryptophan and is primarily unstructured. However, it adopts the canonical triple-stranded antiparallel β-sheet structure of WW domains when bound to a two-PPXY motif peptide of the synaptic protein Dendrin. Here, using a series of biophysical experiments, we demonstrate that the WW2 domain remains largely disordered when bound to a 69-residue two-PPXY motif polypeptide of the synaptic and podocyte protein synaptopodin (SYNPO). Isothermal titration calorimetry and CD experiments revealed that the interactions of the disordered WW2 domain with SYNPO are significantly weaker than SYNPO's interactions with the well-folded WW1 domain and that an I81W substitution in the WW2 domain neither enhances binding affinity nor induces substantial WW2 domain folding. In the tandem polypeptide, the two WW domains synergized, enhancing the overall binding affinity with the I81W variant tandem polypeptide 2-fold compared with the WT polypeptide. Solution NMR results showed that SYNPO binding induces small but definite chemical shift perturbations in the WW2 domain, confirming the disordered state of the WW2 domain in this complex. These analyses also disclosed that SYNPO binds the tandem WW domain polypeptide in an antiparallel manner, that is, the WW1 domain binds the second PPXY motif of SYNPO. We propose a binding model consisting of a bipartite interaction mode in which the largely disordered WW2 forms a "fuzzy" complex with SYNPO. This binding mode may be important for specific cellular functions.

Observation of Fast Two-Dimensional NMR Spectra during Protein Folding Using Polarization Transfer from Hyperpolarized Water

Nuclear spin hyperpolarized water is utilized to obtain protein spectra not only in the folded state but also during the refolding process. Polarization transfer to Ribonuclease Sa through proton exchange and the nuclear Overhauser effect (NOE) results in NMR signal enhancements of amide protons by up to 24-fold. These enhancements enable the measurement of fast two-dimensional NMR spectra on the same time scale as the folding. Resolved amide proton signals corresponding to the folded protein are observed both under folded and refolding conditions, whereby the refolding protein shows smaller transferred signals. Residue-specific evaluation of contributions to the polarization transfer indicates that signals attributed to a relayed intramolecular NOE are not observable in the refolding experiment. These differences are explained by the absence of long-range contacts and faster molecular motions in the unfolded protein. Applications of this method include accessing residue-specific information on structure and dynamics during multistate protein folding.

The Cytoplasm-Entry Domain of Antibacterial CdiA Is a Dynamic α-Helical Bundle with Disulfide-Dependent Structural Features

Many Gram-negative bacterial species use contact-dependent growth inhibition (CDI) systems to compete with neighboring cells. CDI+ strains express cell-surface CdiA effector proteins, which carry a toxic C-terminal region (CdiA-CT) that is cleaved from the effector upon transfer into the periplasm of target bacteria. The released CdiA-CT consists of two domains. The C-terminal domain is typically a nuclease that inhibits cell growth, and the N-terminal "cytoplasm-entry" domain mediates toxin translocation into the target-cell cytosol. Here, we use NMR and circular dichroism spectroscopic approaches to probe the structure, stability, and dynamics of the cytoplasm-entry domain from Escherichia coli STEC_MHI813. Chemical shift analysis reveals that the CdiA-CTMHI813 entry domain is composed of a C-terminal helical bundle and a dynamic N-terminal region containing two disulfide linkages. Disruption of the disulfides by mutagenesis or chemical reduction destabilizes secondary structure over the N-terminus, but has no effect on the C-terminal helices. Although critical for N-terminal structure, the disulfides have only modest effects on global thermodynamic stability, and the entry domain exhibits characteristics of a molten globule. We find that the disulfides form in vivo as the entry domain dwells in the periplasm of inhibitor cells prior to target-cell recognition. CdiA-CTMHI813 variants lacking either disulfide still kill target bacteria, but disruption of both bonds abrogates growth inhibition activity. We propose that the entry domain's dynamic structural features are critical for function. In its molten globule-like state, the domain resists degradation after delivery, yet remains pliable enough to unfold for membrane translocation.

Solution Structure of the Carboxy-Terminal Tandem Repeat Domain of Eukaryotic Elongation Factor 2 Kinase and Its Role in Substrate Recognition

Eukaryotic elongation factor 2 kinase (eEF-2K), an atypical calmodulin-activated protein kinase, regulates translational elongation by phosphorylating its substrate, eukaryotic elongation factor 2 (eEF-2), thereby reducing its affinity for the ribosome. The activation and activity of eEF-2K are critical for survival under energy-deprived conditions and is implicated in a variety of essential physiological processes. Previous biochemical experiments have indicated that the binding site for the substrate eEF-2 is located in the C-terminal domain of eEF-2K, a region predicted to harbor several α-helical repeats. Here, using NMR methodology, we have determined the solution structure of a C-terminal fragment of eEF-2K, eEF-2K562-725 that encodes two α-helical repeats. The structure of eEF-2K562-725 shows signatures characteristic of TPR domains and of their SEL1-like sub-family. Furthermore, using the analyses of NMR spectral perturbations and ITC measurements, we have localized the eEF-2 binding site on eEF-2K562-725. We find that eEF-2K562-725 engages eEF-2 with an affinity comparable to that of the full-length enzyme. Furthermore, eEF-2K562-725 is able to inhibit the phosphorylation of eEF-2 by full-length eEF-2K in trans. Our present studies establish that eEF-2K562-725 encodes the major elements necessary to enable the eEF-2K/eEF-2 interactions.

Reorientational Dynamics of Amyloid-β from NMR Spin Relaxation and Molecular Simulation

Amyloid-β (Aβ) aggregation is a hallmark of Alzheimer's disease. As an intrinsically disordered protein, Aβ undergoes extensive dynamics on multiple length and time scales. Access to a comprehensive picture of the reorientational dynamics in Aβ requires therefore the combination of complementary techniques. Here, we integrate ¹⁵N spin relaxation rates at three magnetic fields with microseconds-long molecular dynamics simulation, ensemble-based hydrodynamic calculations, and previously published nanosecond fluorescence correlation spectroscopy to investigate the reorientational dynamics of Aβ1-40 (Aβ40) at single-residue resolution. The integrative analysis shows that librational and dihedral angle fluctuations occurring at fast and intermediate time scales are not sufficient to decorrelate orientational memory in Aβ40. Instead, slow segmental motions occurring at ∼5 ns are detected throughout the Aβ40 sequence and reach up to ∼10 ns for selected residues. We propose that the modulation of time scales of reorientational dynamics with respect to intra- and intermolecular diffusion plays an important role in disease-related Aβ aggregation.

Insights into Allosteric Control of Human Blood Group A and B Glycosyltransferases from Dynamic NMR

Human blood group A and B glycosyltransferases (GTA, GTB) are retaining glycosyltransferases, requiring a catalytic mechanism that conserves the anomeric configuration of the hexopyranose moiety of the donor substrate (UDP-GalNAc, UDP-Gal). Previous studies have shown that GTA and GTB cycle through structurally distinct states during catalysis. Here, we link binding and release of substrates, substrate-analogs, and products to transitions between open, semi-closed, and closed states of the enzymes. Methyl TROSY based titration experiments in combination with zz-exchange experiments uncover dramatic changes of binding kinetics associated with allosteric interactions between donor-type and acceptor-type ligands. Taken together, this highlights how allosteric control of on- and off-rates correlates with conformational changes, driving catalysis to completion.

Dynamic Studies on Intrinsically Disordered Regions of Two Paralogous Transcription Factors Reveal Rigid Segments with Important Biological Functions

Long stretches of intrinsically disordered regions (IDRs) are abundantly present in eukaryotic transcription factors. Although their biological significance is well appreciated, the underlying structural and dynamic mechanisms of their function are still not clear. Using solution NMR spectroscopy, we have studied the structural and dynamic features of two paralogous HOX transcription factors, SCR and DFD, from Drosophila. Both proteins have a conserved DNA-binding homeodomain and a long stretch of functionally important IDR. Using NMR dynamics, we determined flexibility of each residue in these proteins. The flexibility of the residues in the disordered region is not uniform. In both proteins, the IDRs have short stretches of consecutive residues with relatively less flexibility, that is, higher rigidity. We show that one such rigid segment is specifically recognized by another co-transcription factor, thus highlighting the importance of these rigid segments in IDR-mediated protein-protein interactions. Using molecular dynamics simulation, we further show that the rigid segments sample less conformations compared to the rest of the residues in the disordered region. The restrained conformational sampling of these rigid residues should lower the loss in conformational entropy during their interactions with binding partners resulting in sequence specific binding. This work provides experimental evidence of a "rigid-segment" model of IDRs, where functionally important rigid segments are connected by highly flexible linkers. Furthermore, a comparative study of IDRs in paralogous proteins reveals that in spite of low-sequence conservation, the rigid and flexible segments are sequentially maintained to preserve related functions and regulations of these proteins.

Chemical Exchange

The phenomenon of chemical or conformational exchange in NMR spectroscopy has enabled detailed characterization of time-dependent aspects of biomolecular function, including folding, molecular recognition, allostery, and catalysis, on timescales from microsecond to second. Importantly, NMR methods based on a variety of spin relaxation parameters have been developed that provide quantitative information on interconversion kinetics, thermodynamic properties, and structural features of molecular states populated to a fraction of a percent at equilibrium and otherwise unobservable by other NMR approaches. The ongoing development of more sophisticated experimental techniques and the necessity to apply these methods to larger and more complex molecular systems engenders a corresponding need for theoretical advances describing such techniques and facilitating data analysis in applications. This review surveys current aspects of the theory of chemical exchange, as utilized in ZZ-exchange; Hahn and Carr-Purcell-Meiboom-Gill (CPMG) spin-echo; and R_1ρ, chemical exchange saturation transfer (CEST), and dark state saturation transfer (DEST) spin-locking experiments. The review emphasizes theoretical results for kinetic topologies with more than two interconverting states, both to obtain compact analytical forms suitable for data analysis and to establish conditions for distinguishability between alternative kinetic schemes.

Local and Global Dynamics in Intrinsically Disordered Synuclein

Intrinsically disordered proteins (IDPs) experience a diverse spectrum of motions that are difficult to characterize with a single experimental technique. Herein we combine high- and low-field nuclear spin relaxation, nanosecond fluorescence correlation spectroscopy (nsFCS), and long molecular dynamics simulations of alpha-synuclein, an IDP involved in Parkinson disease, to obtain a comprehensive picture of its conformational dynamics. The combined analysis shows that fast motions below 2 ns caused by local dihedral angle fluctuations and conformational sampling within and between Ramachandran substates decorrelate most of the backbone N-H orientational memory. However, slow motions with correlation times of up to ca. 13 ns from segmental dynamics are present throughout the alpha-synuclein chain, in particular in its C-terminal domain, and global chain reconfiguration occurs on a timescale of ca. 60 ns. Our study demonstrates a powerful strategy to determine residue-specific protein dynamics in IDPs at different time and length scales.

Magnesium Activates Microsecond Dynamics to Regulate Integrin-Collagen Recognition

Integrin receptors bind collagen via metal-mediated interactions that are modulated by magnesium (Mg²⁺) levels in the extracellular matrix. Nuclear magnetic resonance-based relaxation experiments, isothermal titration calorimetry, and adhesion assays reveal that Mg²⁺ functions as both a structural anchor and dynamic switch of the α₁β₁ integrin I domain (α₁I). Specifically, Mg²⁺ binding activates micro- to millisecond timescale motions of residues distal to the binding site, particularly those surrounding the salt bridge at helix 7 and near the metal ion-dependent adhesion site. Mutagenesis of these residues impacts α₁I functional activity, thereby suggesting that Mg-bound α₁I dynamics are important for collagen binding and consequent allosteric rearrangement of the low-affinity closed to high-affinity open conformation. We propose a multistep recognition mechanism for α₁I-Mg-collagen interactions involving both conformational selection and induced-fit processes. Our findings unravel the multifaceted role of Mg²⁺ in integrin-collagen recognition and assist in elucidating the molecular mechanisms by which metals regulate protein-protein interactions.

NMR Structure, Dynamics and Interactions of the Integrin β2 Cytoplasmic Tail with Filamin Domain IgFLNa21

Integrins are transmembrane proteins that mediate cell adhesion and migration. Each integrin is a heterodimer formed by an α and a β subunit. A large number of cytoplasmic proteins interact with the cytoplasmic tails (CTs) of integrins. The actin-binding cytoskeletal protein filamin A is a negative regulator of integrin activation. The IgFLNa21 domain of filamin A binds to the C-terminus of β2 CT that contains a TTT-motif. Based on x-ray crystallography, it has been reported that the integrin β2 CT forms a β strand that docks into the β strands C and D of IgFLNa21. In this study, we performed solution NMR analyses of IgFLNa21 in the presence of integrin β2 CT peptides, and hybrid IgFLNa21, a construct of covalently linked IgFLNa21 and β2 CT. The atomic resolution structure of the hybrid IgFLNa21 demonstrated conserved binding mode with β2 CT. Although, ¹⁵N relaxation, model free analyses and H-D exchange studies have uncovered important insights into the conformational dynamics and stability of β2 CT in complex with IgFLNa21. Such dynamical characteristics are likely to be necessary for the TTT-motif to serve as a phosphorylation switch that regulates filamin A binding to integrin β2 CT.

Atomic resolution conformational dynamics of intrinsically disordered proteins from NMR spin relaxation

Nuclear magnetic resonance (NMR) spectroscopy is one of the most powerful experimental approaches for investigating the conformational behaviour of intrinsically disordered proteins (IDPs). IDPs represent a significant fraction of all proteomes, and, despite their importance for understanding fundamental biological processes, the molecular basis of their activity still remains largely unknown. The functional mechanisms exploited by IDPs in their interactions with other biomolecules are defined by their intrinsic dynamic modes and associated timescales, justifying the considerable interest over recent years in the development of technologies adapted to measure and describe this behaviour. NMR spin relaxation delivers information-rich, site-specific data reporting on conformational fluctuations occurring throughout the molecule. Here we review recent progress in the use of ¹⁵N relaxation to identify local backbone dynamics and long-range chain-like motions in unfolded proteins.

The PipX Protein, When Not Bound to Its Targets, Has Its Signaling C-Terminal Helix in a Flexed Conformation

PipX, an 89-residue protein, acts as a coactivator of the global nitrogen regulator NtcA in cyanobacteria. NtcA-PipX interactions are regulated by 2-oxoglutarate (2-OG), an inverse indicator of the ammonia abundance, and by P_II, a protein that binds to PipX at low 2-OG concentrations. The structure of PipX, when bound to NtcA or P_II, consists of an N-terminal, five-stranded β-sheet (conforming a Tudor-like domain), and two long α-helices. These helices adopt either a flexed conformation, where they are in close contact and in an antiparallel mutual orientation, also packing against the β-sheet, or an open conformation (observed only in the P_II-PipX complex) where the last α-helix moves apart from the rest of the protein. The aim of this work was to study the structure and dynamics of isolated PipX in solution by NMR. The backbone chemical shifts, the hydrogen-exchange, and the NOE patterns indicated that the isolated, monomeric PipX structure was formed by an N-terminal five-stranded β-sheet and two C-terminal α-helices. Furthermore, the observed NOEs between the two helices, and of α-helix2 with β-strand2 suggested that PipX adopted a flexed conformation. The β-strands 1 and 5 were highly flexible, as shown by the lack of interstrand backbone-backbone NOEs; in addition, the ¹⁵N-dynamics indicated that the C terminus of β-strand4 and the following β-turn (Phe42-Thr47), and the C-cap of α-helix1 (Arg70-Asn71) were particularly mobile. These two regions could act as hinges, allowing PipX to interact with its partners, including PlmA in the newly recognized P_II-PipX-PlmA ternary complex.

Hyperinsulinism-Causing Mutations Cause Multiple Molecular Defects in SUR1 NBD1

The sulfonylurea receptor 1 (SUR1) protein forms the regulatory subunit in ATP sensitive K⁺ (K_ATP) channels in the pancreas. SUR proteins are members of the ATP binding cassette (ABC) superfamily of proteins. Binding and hydrolysis of MgATP at the SUR nucleotide binding domains (NBDs) lead to channel opening. Pancreatic K_ATP channels play an important role in insulin secretion. SUR1 mutations that result in increased levels of channel opening ultimately inhibit insulin secretion and lead to neonatal diabetes. In contrast, SUR1 mutations that disrupt trafficking and/or decrease gating of K_ATP channels cause congenital hyperinsulinism, where oversecretion of insulin occurs even in the presence of low glucose levels. Here, we present data on the effects of specific congenital hyperinsulinism-causing mutations (G716V, R842G, and K890T) located in different regions of the first nucleotide binding domain (NBD1). Nuclear magnetic resonance (NMR) and fluorescence data indicate that the K890T mutation affects residues throughout NBD1, including residues that bind MgATP, NBD2, and coupling helices. The mutations also decrease the MgATP binding affinity of NBD1. Size exclusion and NMR data indicate that the G716V and R842G mutations cause aggregation of NBD1 in vitro, possibly because of destabilization of the domain. These data describe structural characterization of SUR1 NBD1 and shed light on the underlying molecular basis of mutations that cause congenital hyperinsulinism.

Understanding Protein Structure Deformation on the Surface of Gold Nanoparticles of Varying Size

Gold nanoparticles (AuNPs) have been of recent interest due to their unique optical properties and their biocompatibility. Biomolecules spontaneously adsorb to their surface, a trait that could potentially be exploited for drug targeting. Currently, it is unclear whether protein-AuNP interactions at the nanoparticle surface are dependent on nanoparticle size. In this work, we investigate whether varying surface curvature can induce protein unfolding and multilayer binding in citrate-coated AuNPs of various sizes. A recently developed NMR-based approach was utilized to determine the adsorption capacity, and protein NMR spectra were compared to determine whether nanoparticle size influences protein interactions at the surface. In addition, transmission electron microscopy (TEM) and dynamic light scattering (DLS) were employed to corroborate the NMR studies. Over a broad range of AuNP sizes (14-86 nm), we show that adsorption capacity can be predicted by assuming that proteins are compact and globular on the nanoparticle surface. Additionally, roughly one layer of proteins is adsorbed regardless of AuNP size. Our results hold for two proteins of significantly different sizes, GB3 (6 kDa) and bovine carbonic anhydrase (BCA, 29 kDa). However, the unstable drkN SH3 domain (ΔḠ0 ≈ 0, 7 kDa) does not appear to follow the same trend seen for stable, globular proteins. This observation suggests that unstable proteins can deform significantly when bound to AuNP surfaces. Taken together, the results of this work can be used to improve our knowledge of the mechanism of protein-AuNP interactions to optimize their use in the biomedical field.

Quantifying protein dynamics in the ps-ns time regime by NMR relaxation

Both ¹⁵N chemical shift anisotropy (CSA) and sufficiently rapid exchange linebroadening transitions exhibit relaxation contributions that are proportional to the square of the magnetic field. Deconvoluting these contributions is further complicated by residue-dependent variations in protein amide ¹⁵N CSA values which have proven difficult to accurately measure. Exploiting recently reported improvements for the implementation of T₁ and T_1ρ experiments, field strength-dependent studies have been carried out on the B3 domain of protein G (GB3) as well as on the immunophilin FKBP12 and a H87V variant of that protein in which the major conformational exchange linebroadening transition is suppressed. By applying a zero frequency spectral density rescaling analysis to the relaxation data collected at magnetic fields from 500 to 900 MHz ¹H, differential residue-specific ¹⁵N CSA values have been obtained for GB3 which correlate with those derived from solid state and liquid crystalline NMR measurements to a level similar to the correlation among those previously reported studies. Application of this analysis protocol to FKBP12 demonstrated an efficient quantitation of both weak exchange linebroadening contributions and differential residue-specific ¹⁵N CSA values. Experimental access to such differential residue-specific ¹⁵N CSA values should significantly facilitate more accurate comparisons with molecular dynamics simulations of protein motion that occurs within the timeframe of global molecular tumbling.

(13)C-NMR studies on disulfide bond isomerization in bovine pancreatic trypsin inhibitor (BPTI)

Conformational isomerization of disulfide bonds is associated with the dynamics and thus the functional aspects of proteins. However, our understanding of the isomerization is limited by experimental difficulties in probing it. We explored the disulfide conformational isomerization of the Cys14-Cys38 disulfide bond in bovine pancreatic trypsin inhibitor (BPTI), by performing an NMR line-shape analysis of its Cys carbon peaks. In this approach, 1D (13)C spectra were recorded at small temperature intervals for BPTI samples selectively labeled with site-specifically (13)C-enriched Cys, and the recorded peaks were displayed in the order of the temperature after the spectral scales were normalized to a carbon peak. Over the profile of the line-shape, exchange broadening that altered with temperature was manifested for the carbon peaks of Cys14 and Cys38. The Cys14-Cys38 disulfide bond reportedly exists in equilibrium between a high-populated (M) and two low-populated states (m c14 and m c38). Consistent with the three-site exchange model, biphasic exchange broadening arising from the two processes was observed for the peak of the Cys14 α-carbon. As the exchange broadening is maximized when the exchange rate equals the chemical shift difference in Hz between equilibrating sites, semi-quantitative information that was useful for establishing conditions for (13)C relaxation dispersion experiments was obtained through the carbon line-shape profile. With respect to the m c38 isomerization, the (1)H-(13)C signals at the β-position of the minor state were resolved from the major peaks and detected by exchange experiments at a low temperature.

Identification of Dynamic Modes in an Intrinsically Disordered Protein Using Temperature-Dependent NMR Relaxation

The dynamic modes and time scales sampled by intrinsically disordered proteins (IDPs) define their function. Nuclear magnetic resonance (NMR) spin relaxation is probably the most powerful tool for investigating these motions delivering site-specific descriptions of conformational fluctuations from throughout the molecule. Despite the abundance of experimental measurement of relaxation in IDPs, the physical origin of the measured relaxation rates remains poorly understood. Here we measure an extensive range of auto- and cross-correlated spin relaxation rates at multiple magnetic field strengths on the C-terminal domain of the nucleoprotein of Sendai virus, over a large range of temperatures (268-298 K), and combine these data to describe the dynamic behavior of this archetypal IDP. An Arrhenius-type relationship is used to simultaneously analyze up to 61 relaxation rates per amino acid over the entire temperature range, allowing the measurement of local activation energies along the chain, and the assignment of physically distinct dynamic modes. Fast (τ ≤ 50 ps) components report on librational motions, a dominant mode occurs on time scales around 1 ns, apparently reporting on backbone sampling within Ramachandran substates, while a slower component (5-25 ns) reports on segmental dynamics dominated by the chain-like nature of the protein. Extending the study to three protein constructs of different lengths (59, 81, and 124 amino acids) substantiates the assignment of these contributions. The analysis is shown to be remarkably robust, accurately predicting a broad range of relaxation data measured at different magnetic field strengths and temperatures. The ability to delineate intrinsic modes and time scales from NMR spin relaxation will improve our understanding of the behavior and function of IDPs, adding a new and essential dimension to the description of this biologically important and ubiquitous class of proteins.

Effect of protein structure and/or conformation on the dityrosine cross-linking induced by haem-hydrogen peroxide

BACKGROUND

Haem, an essential cofactor in aerobic organisms, can cause oxidative stress and impose toxic effects on tissues and organs. It can induce aggregation of proteins via dityrosine cross-linking and cause neurodegenerative diseases. Although dityrosine cross-linking in many proteins induced by haem has been reported, not all the proteins have the same effect or the efficiency of cross-linking varies, while the reason has not been clarified.

METHODS

The correlation of protein structure/conformation with its aggregation tendency via dityrosine induced by hematin (oxidized form of haem) in the presence of hydrogen peroxide (H2O2) was studied through reducing sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), fluorescence and circular dichroism (CD) measurements, and the mechanism was investigated by performing UV-Vis absorbance, Raman spectroscopy and low-temperature electron spin resonance (ESR) experiments.

RESULTS

It was found that proteins in unstructured state are more readily to be cross-linked via dityrosine formation by hematin-H2O2. The unstructured protein without steric effect can coordinate with hematin to form six-coordinated protein-hematin complex, in which the produced tyrosyl radicals by H2O2 are with high tendency to dimerize to form dityrosine.

CONCLUSIONS

Our results demonstrate that protein structure/conformation can affect its coordination state with haem, and the tendency of reaction of two tyrosyl radicals, further influencing the yield and efficiency of dityrosine cross-linking in the presence of H2O2.

GENERAL SIGNIFICANCE

This research can help to deepen our understanding of the protein aggregation and inactivation mechanisms in varied sophisticated conditions, and especially give us the new insight into the toxic effects under haem stress.

Dynamics of GCN4 facilitate DNA interaction: a model-free analysis of an intrinsically disordered region

Intrinsically disordered proteins (IDPs) and proteins with intrinsically disordered regions (IDRs) are known to play important roles in regulatory and signaling pathways. A critical aspect of these functions is the ability of IDP/IDRs to form highly specific complexes with target molecules. However, elucidation of the contributions of conformational dynamics to function has been limited by challenges associated with structural heterogeneity of IDP/IDRs. Using NMR spin relaxation parameters ((15)N R1, (15)N R2, and {(1)H}-(15)N heteronuclear NOE) collected at four static magnetic fields ranging from 14.1 to 21.1 T, we have analyzed the backbone dynamics of the basic leucine-zipper (bZip) domain of the Saccharomyces cerevisiae transcription factor GCN4, whose DNA binding domain is intrinsically disordered in the absence of DNA substrate. We demonstrate that the extended model-free analysis can be applied to proteins with IDRs such as apo GCN4 and that these results significantly extend previous NMR studies of GCN4 dynamics performed using a single static magnetic field of 11.74 T [Bracken, et al., J. Mol. Biol., 1999, 285, 2133-2146] and correlate well with molecular dynamics simulations [Robustelli, et al., J. Chem. Theory Comput., 2013, 9, 5190-5200]. In contrast to the earlier work, data at multiple static fields allows the time scales of internal dynamics of GCN4 to be reliably quantified. Large amplitude dynamic fluctuations in the DNA-binding region have correlation times (τs ≈ 1.4-2.5 ns) consistent with a two-step mechanism in which partially ordered bZip conformations of GCN4 form initial encounter complexes with DNA and then rapidly rearrange to the high affinity state with fully formed basic region recognition helices.

In-cell thermodynamics and a new role for protein surfaces

There is abundant, physiologically relevant knowledge about protein cores; they are hydrophobic, exquisitely well packed, and nearly all hydrogen bonds are satisfied. An equivalent understanding of protein surfaces has remained elusive because proteins are almost exclusively studied in vitro in simple aqueous solutions. Here, we establish the essential physiological roles played by protein surfaces by measuring the equilibrium thermodynamics and kinetics of protein folding in the complex environment of living Escherichia coli cells, and under physiologically relevant in vitro conditions. Fluorine NMR data on the 7-kDa globular N-terminal SH3 domain of Drosophila signal transduction protein drk (SH3) show that charge-charge interactions are fundamental to protein stability and folding kinetics in cells. Our results contradict predictions from accepted theories of macromolecular crowding and show that cosolutes commonly used to mimic the cellular interior do not yield physiologically relevant information. As such, we provide the foundation for a complete picture of protein chemistry in cells.

Investigating protein conformational energy landscapes and atomic resolution dynamics from NMR dipolar couplings: a review

Nuclear magnetic resonance spectroscopy is exquisitely sensitive to protein dynamics. In particular inter-nuclear dipolar couplings, that become measurable in solution when the protein is dissolved in a dilute liquid crystalline solution, report on all conformations sampled up to millisecond timescales. As such they provide the opportunity to describe the Boltzmann distribution present in solution at atomic resolution, and thereby to map the conformational energy landscape in unprecedented detail. The development of analytical methods and approaches based on numerical simulation and their application to numerous biologically important systems is presented.

Long range Trp-Trp interaction initiates the folding pathway of a pro-angiogenic β-hairpin peptide

HPLW, a designed VEGF (Vascular Endothelium Growth Factor) receptor-binding peptide, assumes a well folded β-hairpin conformation in water and is able to induce angiogenesis in vivo. In this study, we investigated at atomic resolution the thermal folding/unfolding pathway of HPLW by means of an original multi-technique approach combining DSC, NMR, MD and mutagenesis analyses. In particular, careful NMR investigation of the single proton melting temperatures together with DSC analysis accurately delineate the peptide folding mechanism, which is corroborated by computational folding/unfolding simulations. The HPLW folding process consists of two main events, which are successive but do not superimpose. The first folding step initiates at 320 K upon the hydrophobic collapse of the Trp5 and Trp13 side-chains which stabilizes the concurrent β-turn formation, whose COi-HNi + 3 hydrogen bond (Asp10 → Arg7) appears particularly stable. At 316 K, once the β-turn is completely formed, the two β-strands pair, very likely starting by Trp5 and Trp13, which thus play a key role also in the final step of the β-hairpin folding. Overall, here we describe a multi-state hierarchical folding pathway of a highly structured β-hairpin, which can be classified as a broken-zipper mechanism.

Long-range correlated dynamics in intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are involved in a wide variety of physiological and pathological processes and are best described by ensembles of rapidly interconverting conformers. Using fast field cycling relaxation measurements we here show that the IDP α-synuclein as well as a variety of other IDPs undergoes slow reorientations at time scales comparable to folded proteins. The slow motions are not perturbed by mutations in α-synuclein, which are related to genetic forms of Parkinson's disease, and do not depend on secondary and tertiary structural propensities. Ensemble-based hydrodynamic calculations suggest that the time scale of the underlying correlated motion is largely determined by hydrodynamic coupling between locally rigid segments. Our study indicates that long-range correlated dynamics are an intrinsic property of IDPs and offers a general physical mechanism of correlated motions in highly flexible biomolecular systems.