Convergent views on disordered protein dynamics from NMR and computational approaches.

Intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs) is a class of biologically important proteins exhibiting specific biophysical characteristics. They lack a hydrophobic core, and their conformational behavior is strongly influenced by electrostatic interactions. IDPs and IDRs are highly dynamic, and a characterization of the motions of IDPs and IDRs is essential for their physically correct description. NMR together with molecular dynamics simulations are the methods best suited to such a task because they provide information about dynamics of proteins with atomistic resolution. Here, we present a study of motions of a disordered C-terminal domain of the delta subunit of RNA polymerase from Bacillus subtilis. Positively and negatively charged residues in the studied domain form transient electrostatic contacts critical for the biological function. Our study is focused on investigation of ps-ns dynamics of backbone of the delta subunit based on analysis of amide 15N NMR relaxation data and molecular dynamics simulations. In order to extend an informational content of NMR data to lower frequencies, which are more sensitive to slower motions, we combined standard (high-field) NMR relaxation experiments with high-resolution relaxometry. Altogether, we collected data reporting the relaxation at 12 different magnetic fields, resulting in an unprecedented data set. Our results document that the analysis of such data provides a consistent description of dynamics and confirms the validity of so far used protocols of the analysis of dynamics of IDPs also for a partially folded protein. In addition, the potential to access detailed description of motions at the timescale of tens of ns with the help of relaxometry data is discussed. Interestingly, in our case, it appears to be mostly relevant for a region involved in the formation of temporary contacts within the disordered region, which was previously proven to be biologically important.

Choice of Force Field for Proteins Containing Structured and Intrinsically Disordered Regions.

Biomolecular force fields optimized for globular proteins fail to properly reproduce properties of intrinsically disordered proteins. In particular, parameters of the water model need to be modified to improve applicability of the force fields to both ordered and disordered proteins. Here, we compared performance of force fields recommended for intrinsically disordered proteins in molecular dynamics simulations of three proteins differing in the content of ordered and disordered regions (two proteins consisting of a well-structured domain and of a disordered region with and without a transient helical motif and one disordered protein containing a region of increased helical propensity). The obtained molecular dynamics trajectories were used to predict measurable parameters, including radii of gyration of the proteins and chemical shifts, residual dipolar couplings, paramagnetic relaxation enhancement, and NMR relaxation data of their individual residues. The predicted quantities were compared with experimental data obtained within this study or published previously. The results showed that the NMR relaxation parameters, rarely used for benchmarking, are particularly sensitive to the choice of force-field parameters, especially those defining the water model. Interestingly, the TIP3P water model, leading to an artificial structural collapse, also resulted in unrealistic relaxation properties. The TIP4P-D water model, combined with three biomolecular force-field parameters for the protein part, significantly improved reliability of the simulations. Additional analysis revealed only one particular force field capable of retaining the transient helical motif observed in NMR experiments. The benchmarking protocol used in our study, being more sensitive to imperfections than the commonly used tests, is well suited to evaluate the performance of newly developed force fields.

Boosting the resolution of low-field [Formula: see text] relaxation experiments on intrinsically disordered proteins with triple-resonance NMR.

Improving our understanding of nanosecond motions in disordered proteins requires the enhanced sampling of the spectral density function obtained from relaxation at low magnetic fields. High-resolution relaxometry and two-field NMR measurements of relaxation have, so far, only been based on the recording of one- or two-dimensional spectra, which provide insufficient resolution for challenging disordered proteins. Here, we introduce a 3D-HNCO-based two-field NMR experiment for measurements of protein backbone [Formula: see text] amide longitudinal relaxation rates. The experiment provides accurate longitudinal relaxation rates at low field (0.33 T in our case) preserving the resolution and sensitivity typical for high-field NMR spectroscopy. Radiofrequency pulses applied on six different radiofrequency channels are used to manipulate the spin system at both fields. The experiment was demonstrated on the C-terminal domain of [Formula: see text] subunit of RNA polymerase from Bacillus subtilis, a protein with highly repetitive amino-acid sequence and very low dispersion of backbone chemical shifts.

Quantitative Conformational Analysis of Functionally Important Electrostatic Interactions in the Intrinsically Disordered Region of Delta Subunit of Bacterial RNA Polymerase.

Electrostatic interactions play important roles in the functional mechanisms exploited by intrinsically disordered proteins (IDPs). The atomic resolution description of long-range and local structural propensities that can both be crucial for the function of highly charged IDPs presents significant experimental challenges. Here, we investigate the conformational behavior of the δ subunit of RNA polymerase from Bacillus subtilis whose unfolded domain is highly charged, with 7 positively charged amino acids followed by 51 acidic amino acids. Using a specifically designed analytical strategy, we identify transient contacts between the two regions using a combination of NMR paramagnetic relaxation enhancements, residual dipolar couplings (RDCs), chemical shifts, and small-angle scattering. This strategy allows the resolution of long-range and local ensemble averaged structural contributions to the experimental RDCs, and reveals that the negatively charged segment folds back onto the positively charged strand, compacting the conformational sampling of the protein while remaining highly flexible in solution. Mutation of the positively charged region abrogates the long-range contact, leaving the disordered domain in an extended conformation, possibly due to local repulsion of like-charges along the chain. Remarkably, in vitro studies show that this mutation also has a significant effect on transcription activity, and results in diminished cell fitness of the mutated bacteria in vivo. This study highlights the importance of accurately describing electrostatic interactions for understanding the functional mechanisms of IDPs.

Conformational dynamics are a key factor in signaling mediated by the receiver domain of a sensor histidine kinase from Arabidopsis thaliana.

Multistep phosphorelay (MSP) cascades mediate responses to a wide spectrum of stimuli, including plant hormonal signaling, but several aspects of MSP await elucidation. Here, we provide first insight into the key step of MSP-mediated phosphotransfer in a eukaryotic system, the phosphorylation of the receiver domain of the histidine kinase CYTOKININ-INDEPENDENT 1 (CKI1RD) from Arabidopsis thaliana We observed that the crystal structures of free, Mg2+-bound, and beryllofluoridated CKI1RD (a stable analogue of the labile phosphorylated form) were identical and similar to the active state of receiver domains of bacterial response regulators. However, the three CKI1RD variants exhibited different conformational dynamics in solution. NMR studies revealed that Mg2+ binding and beryllofluoridation alter the conformational equilibrium of the ß3-α3 loop close to the phosphorylation site. Mutations that perturbed the conformational behavior of the ß3-α3 loop while keeping the active-site aspartate intact resulted in suppression of CKI1 function. Mechanistically, homology modeling indicated that the ß3-α3 loop directly interacts with the ATP-binding site of the CKI1 histidine kinase domain. The functional relevance of the conformational dynamics observed in the ß3-α3 loop of CKI1RD was supported by a comparison with another A. thaliana histidine kinase, ETR1. In contrast to the highly dynamic ß3-α3 loop of CKI1RD, the corresponding loop of the ETR1 receiver domain (ETR1RD) exhibited little conformational exchange and adopted a different orientation in crystals. Biochemical data indicated that ETR1RD is involved in phosphorylation-independent signaling, implying a direct link between conformational behavior and the ability of eukaryotic receiver domains to participate in MSP.