Three-dimensional structure of the histidine-containing phosphocarrier protein (HPr) from Enterococcus faecalis in solution

The histidine-containing phosphocarrier protein (HPr) transfers a phosphate group between components of the prokaryotic phosphoenolpyruvate-dependent phosphotransferase system (PTS), which is finally used to phosphorylate the carbohydrate transported by the PTS through the cell membrane. Recently it has also been found to act as an intermediate in the signaling cascade that regulates transcription of genes related to the carbohydrate-response system. Both functions involve phosphorylation/dephosphorylation reactions, but at different sites. Using multidimensional (1)H-NMR spectroscopy and angular space simulated annealing calculations, we determined the structure of HPr from Enterococcus faecalis in aqueous solution using 1469 distance and 44 angle constraints derived from homonuclear NMR data. It has a similar overall fold to that found in HPrs from other organisms. Four beta strands, A, B, C, D, encompassing residues 2-7, 32-37, 40-42 and 60-66, form an antiparallel beta sheet lying opposite the two antiparallel alpha helices, a and c (residues 16-26 and 70-83). A short alpha helix, b, from residues 47-53 is also observed. The pairwise root mean square displacement for the backbone heavy atoms of the mean of the 16 NMR structures to the crystal structure is 0.164 nm. In contrast with the crystalline state, in which a torsion angle strain in the active-center loop has been described [Jia, Z., Vandonselaar, M., Quail, J.W. & Delbaere, L.T.J. (1993) Nature (London) 361, 94-97], in the solution structure, the active-site His15 rests on top of helix a, and the phosphorylation site N(delta 1) of the histidine ring is oriented towards the surface, making it easily accessible to the solvent. Back calculation of the 2D NOESY NMR spectra from both the NMR and X-ray structures shows that the active-center structure derived by X-ray crystallography is not compatible with experimental data recorded in solution. The observed torsional strain must either be a crystallization artefact or represents a conformational state that exists only to a small extent in solution.

RFAC, a program for automated NMR R-factor estimation

A computer program (RFAC) has been developed, which allows the automated estimation of residual indices (R-factors) for protein NMR structures and gives a reliable measure for the quality of the structures. The R-factor calculation is based on the comparison of experimental and simulated 1H NOESY NMR spectra. The approach comprises an automatic peak picking and a Bayesian analysis of the data, followed by an automated structure based assignment of the NOESY spectra and the calculation of the R-factor. The major difference to previously published R-factor definitions is that we take the non-assigned experimental peaks into account as well. The number and the intensities of the non-assigned signals are an important measure for the quality of an NMR structure. It turns out that for different problems optimally adapted R-factors should be used which are defined in the paper. The program allows to compute a global R-factor, different R-factors for the intra residual NOEs, the inter residual NOEs, sequential NOEs, medium range NOEs and long range NOEs. Furthermore, R-factors can be calculated for various user defined parts of the molecule or it is possible to obtain a residue-by-residue R-factor. Another possibility is to sort the R-factors according to their corresponding distances. The summary of all these different R-factors should allow the user to judge the structure in detail. The new program has been successfully tested on two medium sized proteins, the cold shock protein (TmCsp) from Termotoga maritima and the histidine containing protein (HPr) from Staphylococcus carnosus. A comparison with a previously published R-factor definition shows that our approach is more sensitive to errors in the calculated structure.

15N and 1H NMR study of histidine containing protein (HPr) from Staphylococcus carnosus at high pressure

The pressure-induced changes in 15N enriched HPr from Staphylococcus carnosus were investigated by two-dimensional (2D) heteronuclear NMR spectroscopy at pressures ranging from atmospheric pressure up to 200 MPa. The NMR experiments allowed the simultaneous observation of the backbone and side-chain amide protons and nitrogens. Most of the resonances shift downfield with increasing pressure indicating generalized pressure-induced conformational changes. The average pressure-induced shifts for amide protons and nitrogens are 0.285 ppm GPa(-1) at 278 K and 2.20 ppm GPa(-1), respectively. At 298 K the corresponding values are 0.275 and 2.41 ppm GPa(-1). Proton and nitrogen pressure coefficients show a significant but rather small correlation (0.31) if determined for all amide resonances. When restricting the analysis to amide groups in the beta-pleated sheet, the correlation between these coefficients is with 0.59 significantly higher. As already described for other proteins, the amide proton pressure coefficients are strongly correlated to the corresponding hydrogen bond distances, and thus are indicators for the pressure-induced changes of the hydrogen bond lengths. The nitrogen shift changes appear to sense other physical phenomena such as changes of the local backbone conformation as well. Interpretation of the pressure-induced shifts in terms of structural changes in the HPr protein suggests the following picture: the four-stranded beta-pleated sheet of HPr protein is the least compressible part of the structure showing only small pressure effects. The two long helices a and c show intermediary effects that could be explained by a higher compressibility and a concomitant bending of the helices. The largest pressure coefficients are found in the active center region around His15 and in the regulatory helix b which includes the phosphorylation site Ser46 for the HPr kinase. This suggests that this part of the structure occurs in a number of different structural states whose equilibrium populations are shifted by pressure. In contrast to the surrounding residues of the active center loop that show large pressure effects, Ile14 has a very small proton and nitrogen pressure coefficient. It could represent some kind of anchoring point of the active center loop that holds it in the right place in space, whereas other parts of the loop adapt themselves to changing external conditions.

Determination of the populations and structures of multiple conformers in an ensemble from NMR data: multiple-copy refinement of nucleic acid structures using floating weights

A new algorithm is presented for determination of structural conformers and their populations based on NMR data. Restrained Metropolis Monte Carlo simulations or restrained energy minimizations are performed for several copies of a molecule simultaneously. The calculations are restrained with dipolar relaxation rates derived from measured NOE intensities via complete relaxation matrix analysis. The novel feature of the algorithm is that the weights of individual conformers are determined in every refinement step, by the quadratic programming algorithm, in such a way that the restraint energy is minimized. Its design ensures that the calculated populations of the individual conformers are based only on experimental restraints. Presence of internally inconsistent restraints is the driving force for determination of distinct multiple conformers. The method is applied to various simulated test systems. Conformational calculations on nucleic acids are carried out using generalized helical parameters with the program DNAminiCarlo. From different mixtures of A- and B-DNA, minor fractions as low as 10% could be determined with restrained energy minimization. For B-DNA with three local conformers (C2'-endo, O4'-exo, C3'-endo), the minor O4'-exo conformer could not be reliably determined using NOE data typically measured for DNA. The other two conformers, C2'-endo and C3'-endo, could be reproduced by Metropolis Monte Carlo simulated annealing. The behavior of the algorithm in various situations is analyzed, and a number of refinement protocols are discussed. Prior to application of this algorithm to each experimental system, it is suggested that the presence of internal inconsistencies in experimental data be ascertained. In addition, because the performance of the algorithm depends on the type of conformers involved and experimental data available, it is advisable to carry out test calculations with simulated data modeling each experimental system studied.

Computer assisted assignment of 13C or 15N edited 3D-NOESY-HSQC spectra using back calculated and experimental spectra

A new tool, for the simulation of 15N or 13C edited 3D-NOESY-HSQC spectra using the complete relaxation matrix approach, has been developed and integrated in the program AURELIA. This tool should be particularly useful for the fast and reliable computer assisted assignment of 3D-NOESY-HSQC spectra by comparing back-calculated and experimental spectra in an iterative process. Folded spectra are sometimes used to enhance the digital resolution in the indirect dimensions of multidimensional spectra. However, these spectra are usually difficult to analyze. To simplify this assignment process we have incorporated the simulation and automated annotation of folded peaks into the program. It is hereby possible to simulate multiple folding in all three dimensions of 3D 15N- or 13C-NOESY-HSQC spectra. By comparing experimental 3D-NOESY-HSQC spectra with spectra back calculated from a single trial structure or a set of trial structures, a user can easily check if the final structures explain all experimental NOEs. The new feature has been successfully tested with the histidine-containing phosphocarrier protein HPr from Staphylococcus carnosus.

Relax, a flexible program for the back calculation of NOESY spectra based on complete-relaxation-matrix formalism

RELAX is a flexible program for the quantitative analysis of NOESY spectra. It allows the simultaneous application of different models describing the internal and overall motion of the molecule under investigation for individual spin pairs or groups of spins. A correction for anisotropy effects due to the deviation of the molecule from a spherical shape is calculated automatically from the trial structure. The program can deal with completely relaxed spectra as well as spectra recorded with a short relaxation delay. An execution-time-controlled splitting of the relaxation matrix reduces the computation time significantly without any loss of accuracy. This is especially important for large molecules or medium distance cutoffs.

AURELIA, a program for computer-aided analysis of multidimensional NMR spectra

AURELIA is an advanced program for the computer-aided evaluation of two-, three- and four-dimensional NMR spectra of any type of molecule. It can be used for the analysis of spectra of small molecules as well as for evaluation of complicated spectra of biological macromolecules such as proteins. AURELIA is highly interactive and offers a large number of tools, such as artefact reduction, cluster and multiplet analysis, spin system searches, resonance assignments, automated calculation of volumes in multidimensional spectra, calculation of distances with different approaches, including the full relaxation matrix approach, Bayesian analysis of peak features, correlation of molecular structures with NMR data, comparison of spectra via spectral algebra and pattern match techniques, automated sequential assignments on the basis of triple resonance spectra, and automatic strip calculation. In contrast to most other programs, many tasks are performed automatically.