Characterization of a transitionally occupied state and thermal unfolding of domain 1.1 of σ A factor of RNA polymerase from Bacillus subtilis.

σ factors are essential parts of bacterial RNA polymerase (RNAP) as they allow to recognize promotor sequences and initiate transcription. Domain 1.1 of vegetative σ factors occupies the primary channel of RNAP and also prevents binding of the σ factor to promoter DNA alone. Here, we show that domain 1.1 of Bacillus subtilis σ A exists in more structurally distinct variants in dynamic equilibrium. The major conformation at room temperature is represented by a previously reported well-folded structure solved by nuclear magnetic resonance (NMR), but 4% of the protein molecules are present in a less thermodynamically favorable state. We show that this population increases with temperature and we predict its significant elevation at higher but still biologically relevant temperatures. We characterized the minor state of the domain 1.1 using specialized methods of NMR. We found that, in contrast to the major state, the detected minor state is partially unfolded. Its propensity to form secondary structure elements is especially decreased for the first and third α helices, while the second α helix and ß strand close to the C-terminus are more stable. We also analyzed thermal unfolding of the domain 1.1 and performed functional experiments with full length σ A and its shortened version lacking domain 1.1 ( σ A _ Δ 1.1 ). The results revealed that while full length σ A increases transcription activity of RNAP with increasing temperature, transcription with σ A _ Δ 1.1 remains constant. In summary, this study reveals conformational dynamics of domain 1.1 and provides a basis for studies of its interaction with RNAP and effects on transcription regulation.

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.

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.

Solution structure of domain 1.1 of the σA factor from Bacillus subtilis is preformed for binding to the RNA polymerase core.

Bacterial RNA polymerase (RNAP) requires σ factors to recognize promoter sequences. Domain 1.1 of primary σ factors (σ1.1) prevents their binding to promoter DNA in the absence of RNAP, and when in complex with RNAP, it occupies the DNA-binding channel of RNAP. Currently, two 3D structures of σ1.1 are available: from Escherichia coli in complex with RNAP and from T. maritima solved free in solution. However, these two structures significantly differ, and it is unclear whether this difference is due to an altered conformation upon RNAP binding or to differences in intrinsic properties between the proteins from these two distantly related species. Here, we report the solution structure of σ1.1 from the Gram-positive bacterium Bacillus subtilis We found that B. subtilis σ1.1 is highly compact because of additional stabilization not present in σ1.1 from the other two species and that it is more similar to E. coli σ1.1. Moreover, modeling studies suggested that B. subtilis σ1.1 requires minimal conformational changes for accommodating RNAP in the DNA channel, whereas T. maritima σ1.1 must be rearranged to fit therein. Thus, the mesophilic species B. subtilis and E. coli share the same σ1.1 fold, whereas the fold of σ1.1 from the thermophile T. maritima is distinctly different. Finally, we describe an intriguing similarity between σ1.1 and δ, an RNAP-associated protein in B. subtilis, bearing implications for the so-far unknown binding site of δ on RNAP. In conclusion, our results shed light on the conformational changes of σ1.1 required for its accommodation within bacterial RNAP.

Spectral density mapping at multiple magnetic fields suitable for (13)C NMR relaxation studies.

Standard spectral density mapping protocols, well suited for the analysis of (15)N relaxation rates, introduce significant systematic errors when applied to (13)C relaxation data, especially if the dynamics is dominated by motions with short correlation times (small molecules, dynamic residues of macromolecules). A possibility to improve the accuracy by employing cross-correlated relaxation rates and on measurements taken at several magnetic fields has been examined. A suite of protocols for analyzing such data has been developed and their performance tested. Applicability of the proposed protocols is documented in two case studies, spectral density mapping of a uniformly labeled RNA hairpin and of a selectively labeled disaccharide exhibiting highly anisotropic tumbling. Combination of auto- and cross-correlated relaxation data acquired at three magnetic fields was applied in the former case in order to separate effects of fast motions and conformational or chemical exchange. An approach using auto-correlated relaxation rates acquired at five magnetic fields, applicable to anisotropically moving molecules, was used in the latter case. The results were compared with a more advanced analysis of data obtained by interpolation of auto-correlated relaxation rates measured at seven magnetic fields, and with the spectral density mapping of cross-correlated relaxation rates. The results showed that sufficiently accurate values of auto- and cross-correlated spectral density functions at zero and (13)C frequencies can be obtained from data acquired at three magnetic fields for uniformly (13)C-labeled molecules with a moderate anisotropy of the rotational diffusion tensor. Analysis of auto-correlated relaxation rates at five magnetic fields represents an alternative for molecules undergoing highly anisotropic motions.

S3EPY: a Sparky extension for determination of small scalar couplings from spin-state-selective excitation NMR experiments.

S3EPY is a Python extension to the program Sparky written to facilitate the assessment of coupling constants from in-phase/antiphase and spin-state-selective excitation (S(3)E) experiments. It enables the routine use of small scalar couplings by automating the coupling evaluation procedure. S3EPY provides an integrated graphical user interface to programs which outputs graphs and the table of determined couplings.

NMR structure of the N-terminal domain of capsid protein from the mason-pfizer monkey virus.

The high-resolution structure of the N-terminal domain (NTD) of the retroviral capsid protein (CA) of Mason-Pfizer monkey virus (M-PMV), a member of the betaretrovirus family, has been determined by NMR. The M-PMV NTD CA structure is similar to the other retroviral capsid structures and is characterized by a six alpha-helix bundle and an N-terminal beta-hairpin, stabilized by an interaction of highly conserved residues, Pro1 and Asp57. Since the role of the beta-hairpin has been shown to be critical for formation of infectious viral core, we also investigated the functional role of M-PMV beta-hairpin in two mutants (i.e., DeltaP1NTDCA and D57ANTDCA) where the salt bridge stabilizing the wild-type structure was disrupted. NMR data obtained for these mutants were compared with those obtained for the wild type. The main structural changes were observed within the beta-hairpin structure; within helices 2, 3, and 5; and in the loop connecting helices 2 and 3. This observation is supported by biochemical data showing different cleavage patterns of the wild-type and the mutated capsid-nucleocapsid fusion protein (CANC) by M-PMV protease. Despite these structural changes, the mutants with disrupted salt bridge are still able to assemble into immature, spherical particles. This confirms that the mutual interaction and topology within the beta-hairpin and helix 3 might correlate with the changes in interaction between immature and mature lattices.

31P chemical shift tensors for canonical and non-canonical conformations of nucleic acids: a DFT study and NMR implications.

31P chemical shift anisotropy (CSA) tensors have been calculated for a set of selected DNA and RNA backbone conformations using density functional theory. The set includes canonical A-RNA, A-DNA, BI-DNA, BII-DNA, ZI-DNA, and ZII-DNA as well as four A-RNA-type, seven non-A-RNA-type, and three non-canonical DNA conformations. Hexahydrated dimethyl phosphate has been employed as a model. The 31P chemical shift tensors obtained are discussed in terms of similarities in the behavior observed for gauche-gauche (gg) and gauche-trans (gt) conformations around the P-O bonds. We show that torsion angles alpha and zeta are major determinants of the isotropic chemical shift deltaiso and of the deltaCSA11 component of the traceless chemical shift tensor, which is revealed in separate ranges of both deltaiso and deltaCSA11 for gg- and gt-conformers, respectively. A clear distinction between the two conformation types has not been found for the deltaCSA22 and deltaCSA33 components, which is attributed to their different directional properties. The 31P CSA tensors exhibit considerable variations resulting in large spans of approximately 16 ppm for deltaCSA11 and approximately 22 ppm for deltaCSA22 and deltaCSA33. We examine the consequences of the CSA variations for predicting the chemical shift changes upon partial alignment deltacsa and for the values of CSA order parameters extracted from the analysis of 31P NMR relaxation data. The theoretical 31P CSA tensors as well as the experimental 31P CSA tensor of barium diethyl phosphate (BDEP) are used to calculate deltacsa for two eclipsed orientations of the CSA and molecular alignment tensors. Percentage differences between the CSA order parameters obtained using the theoretical 31P CSA tensors and the experimental 31P CSA tensor of BDEP, respectively, are also determined.

Structure of Bombyx mori chemosensory protein 1 in solution.

Chemosensory Proteins (CSPs) represent a family of conserved proteins found in insects that may be involved in chemosensory functions. BmorCSP1 is expressed mainly in antennae and legs of the silkworm moth Bombyx mori and was cloned from antennal cDNA. Here we report the determination of the structure of Bombyx mori CSP1 (BmorCSP1) by NMR. The overall fold of BmorCSP1 is globular and comprises six alpha-helices. These helices span residues 10-14, 17-27, 35-49, 57-72, 75-85, and 92-100. The internal hydrophobic sides of the helices are formed mostly by leucine and isoleucine residues and, therefore, well suited to constitute a binding site for hydrophobic ligands.

Internal consistency of NMR data obtained in partially aligned biomacromolecules.

A method of testing structure-related NMR data prior to structure calculations is presented. The test is applicable to second rank tensor interactions (dipolar coupling, anisotropic chemical shielding, and quadrupolar interaction) observed in partially aligned samples of biomacromolecules. The method utilizes the fact that only limited number of frequencies corresponding to the mentioned interactions can be measured independently in a rigid fragment of the macromolecule. Additional values can be predicted as linear combinations of the set of independent frequencies. Internal consistency of sufficiently large sets of frequencies measured in individual molecular fragments is tested by comparing the experimental data with their predicted values. The method requires only knowledge of local geometry (i.e., definition of the interaction tensors in the local coordinate frames of the fragments). No information about the alignment or shape of the molecule is required. The test is best suited for planar fragments. Application to peptide bonds and nucleic acid bases is demonstrated.

Refinement of d(GCGAAGC) hairpin structure using one- and two-bond residual dipolar couplings.

The structure of the 13C,15N-labeled d(GCGAAGC) hairpin, as determined by NMR spectroscopy and refined using molecular dynamics with NOE-derived distances, torsion angles, and residual dipolar couplings (RDCs), is presented. Although the studied molecule is of small size, it is demonstrated that the incorporation of diminutive RDCs can significantly improve local structure determination of regions undefined by the conventional restraints. Very good correlation between the experimental and back-calculated small one- and two-bond 1H-13C, 1H-15N, 13C-13C and 13C-15N coupling constants has been attained. The final structures clearly show typical features of the miniloop architecture. The structure is discussed in context of the extraordinary stability of the d(GCGAAGC) hairpin, which originates from a complex interplay between the aromatic base stacking and hydrogen bonding interactions.

Program MULDER -- a tool for extracting torsion angles from NMR data.

MULDER (Mostly UniversaL Dihedral angle ExtractoR) is a program for extraction of torsion angle information from NMR data. Currently, it can analyze two types of input data: The torsion angle data, where several (3)J-coupling constants and/or interatomic distances are combined in order to reduce the torsion angle ambiguity arising from solving the isolated Karplus (or distance) equation, and the sugar pucker data, where the dynamics of five-membered sugar rings is evaluated by postprocessing the results calculated from (3)J((HH)) coupling constants by program PSEUROT. Program MULDER can be used either as an alternative to r-MD programs in situations where only specific structural features are studied, or as a preparatory tool in connection with full r-MD structure calculation for extraction of unambiguous torsion angle restraints.