Jens P. Linge

Refinement of protein structures in explicit solvent.

We present a CPU efficient protocol for refinement of protein structures in a thin layer of explicit solvent and energy parameters with completely revised dihedral angle terms. Our approach is suitable for protein structures determined by theoretical (e.g., homology modeling or threading) or experimental methods (e.g., NMR). In contrast to other recently proposed refinement protocols, we put a strong emphasis on consistency with widely accepted covalent parameters and computational efficiency. We illustrate the method for NMR structure calculations of three proteins: interleukin‐4, ubiquitin, and crambin. We show a comparison of their structure ensembles before and after refinement in water with and without a force field energy term for the dihedral angles; crambin was also refined in DMSO. Our results demonstrate the significant improvement of structure quality by a short refinement in a thin layer of solvent. Further, they show that a dihedral angle energy term in the force field is beneficial for structure calculation and refinement. We discuss the optimal weight for the energy constant for the backbone angle omega and include an extensive discussion of meaning and relevance of the calculated validation criteria, in particular root mean square Z scores for covalent parameters such as bond lengths. Proteins 2003;50:496–506.

ARIA: automated NOE assignment and NMR structure calculation

MOTIVATION In the light of several ongoing structural genomics projects, faster and more reliable methods for structure calculation from NMR data are in great demand. The major bottleneck in the determination of solution NMR structures is the assignment of NOE peaks (nuclear Overhauser effect). Due to the high complexity of the assignment problem, most NOEs cannot be directly converted into unambiguous inter-proton distance restraints. RESULTS We present version 1.2 of our program ARIA (Ambiguous Restraints for Iterative Assignment) for automated assignment of NOE data and NMR structure calculation. We summarize recent progress in correcting for spin diffusion with a relaxation matrix approach, representing non-bonded interactions in the force field and refining final structures in explicit solvent. We also discuss book-keeping, data exchange with spectra assignment programs and deposition of the analysed experimental data to the databases. AVAILABILITY ARIA 1.2 is available from: http://www.pasteur.fr/recherche/unites/Binfs/aria/. SUPPLEMENTARY INFORMATION XML DTDs (for chemical shifts and NOE crosspeaks), Python scripts for the conversion of various NMR data formats and the results of example calculations using data from the S. cerevisiae HRDC domain are available from: http://www.pasteur.fr/recherche/unites/Binfs/aria/

The three-dimensional structure of the HRDC domain and implications for the Werner and Bloom syndrome proteins.

BACKGROUND The HRDC (helicase and RNaseD C-terminal) domain is found at the C terminus of many RecQ helicases, including the human Werner and Bloom syndrome proteins. RecQ helicases have been shown to unwind DNA in an ATP-dependent manner. However, the specific functional roles of these proteins in DNA recombination and replication are not known. An HRDC domain exists in both of the human RecQ homologues that are implicated in human disease and may have an important role in their function. RESULTS We have determined the three-dimensional structure of the HRDC domain in the Saccharomyces cerevisiae RecQ helicase Sgs1p by nuclear magnetic resonance (NMR) spectroscopy. The structure resembles auxiliary domains in bacterial DNA helicases and other proteins that interact with nucleic acids. We show that a positively charged region on the surface of the Sgs1p HRDC domain can interact with DNA. Structural similarities to bacterial DNA helicases suggest that the HRDC domain functions as an auxiliary domain in RecQ helicases. Homology models of the Werner and Bloom HRDC domains show different surface properties when compared with Sgs1p. CONCLUSIONS The HRDC domain represents a structural scaffold that resembles auxiliary domains in proteins that are involved in nucleic acid metabolism. In Sgs1p, the HRDC domain could modulate the helicase function via auxiliary contacts to DNA. However, in the Werner and Bloom syndrome helicases the HRDC domain may have a role in their functional differences by mediating diverse molecular interactions.

The CCPN project: an interim report on a data model for the NMR community.

A recent workshop discusses the progress toward integrating NMR data into a unifying data model.

NOE assignment with ARIA 2.0: the nuts and bolts.

The assignment of nuclear Overhauser effect (NOE) resonances is the crucial step in determining the three-dimensional structure of biomolecules from nuclear magnetic resonance (NMR) data. Our program, Ambiguous Restraints for Iterative Assignment (ARIA), treats Noe assignment as an integral part of the structure determination process. This chapter briefly outlines the method and discusses how to carry out a complete structure determination project with the new version 2.0 of ARIA. Two new features greatly streamline the procedure: a new graphical user interface (GUI) and the incorporation of the data model of the Collaborative Computing Project for the NMR community (CCPN). The GUI supports the user in setting up and managing a project. The CCPN data model facilitates data exchange with a great variety of other programs. We give practical guidelines for how to use ARIA and how to analyze results.

Improving the quality of protein structures derived by NMR spectroscopy

Biomolecular structures provide the basis for many studies in several research areas such as homology modelling, structure-based drug design and functional genomics. It is an important prerequisite that the structure is reliable in terms of accurate description of the experimental data, and in terms of good quality of local- and overall geometry. Recent surveys indicate that structures solved by NMR-spectroscopy normally are of lower precision than high-resolution X-ray structures. Here, we present a refinement protocol that improves the quality of protein structures determined by NMR-spectroscopy to the level of those determined by high resolution X-ray crystallography in terms of local geometry. The protocol was tested on experimental data of the proteins IL4 and Ubiquitin and on simulated data of the protein Crambin. In almost all aspects, the protocol yielded better results in terms of accuracy and precision. Independent validation of the results for Ubiquitin, using residual dipolar couplings, indicates that the ensemble of NMR structure is substantially improved by the protocol.

Comparative analysis of structural and dynamic properties of the loaded and unloaded hemophore HasA: functional implications.

A heme-acquisition system present in several Gram-negative bacteria requires the secretion of hemophores. These extracellular carrier proteins capture heme and deliver it to specific outer membrane receptors. The Serratia marcescens HasA hemophore is a monodomain protein that binds heme with a very high affinity. Its alpha/beta structure, as that of its binding pocket, has no common features with other iron- or heme-binding proteins. Heme is held by two loops L1 and L2 and coordinated to iron by an unusual ligand pair, H32/Y75. Two independent regions of the hemophore beta-sheet are involved in HasA-HasR receptor interaction. Here, we report the 3-D NMR structure of apoHasA and the backbone dynamics of both loaded and unloaded hemophore. While the overall structure of HasA is very similar in the apo and holo forms, the hemophore presents a transition from an open to a closed form upon ligand binding, through a large movement, of up to 30 A, of loop L1 bearing H32. Comparison of loaded and unloaded HasA dynamics on different time scales reveals striking flexibility changes in the binding pocket. We propose a mechanism by which these structural and dynamic features provide the dual function of heme binding and release to the HasR receptor.

Refinement of the protein backbone angle psi in NMR structure calculations.

Cross-correlated relaxation rates involving the Cα-Hα dipolar interaction and the carbonyl (C′) chemical shift anisotropy (CSA) have been measured using two complementary 3D experiments. We show that the protein backbone angle ψ can be directly refined against such cross-correlated relaxation rates (ΓHα Cα,C′) and the three-bond H/D isotope effect on the Cα chemical shifts (3ΔCα(ND)). By simultaneously using both experimental parameters as restraints during NMR structure calculations, a unique value for the backbone angle ψ is defined. We have applied the new refinement method to the α-Spectrin SH3 domain (a β-sheet protein) and to the Sgs1p HRDC domain (an α-helical protein) and show that the quality of the NMR structures is substantially improved, judging from the atomic coordinate precision and the Ramachandran map. In addition, the ψ-refined NMR structures of the SH3 domain deviate less from the 1.8 Å crystal structure, suggesting an improved accuracy. The proposed refinement method can be used to significantly improve the quality of NMR structures and will be applicable to larger proteins.

StarDOM: From STAR format to XML

StarDOM is a software package for the representation of STAR files as document object models and the conversion of STAR files into XML. This allows interactive navigation by using the Document Object Model representation of the data as well as easy access by XML query languages. As an example application, the entire BioMagResBank has been transformed into XML format. Using an XML query language, statistical queries on the collected NMR data sets can be constructed with very little effort. The BioMagResBank/XML data and the software can be obtained at http://www.nmr.embl-heidelberg.de/nmr/StarDOM/