A general method for the unbiased improvement of solution NMR structures by the use of related X-ray data, the AUREMOL-ISIC algorithm

Accurate protein structure modeling using sparse NMR data and homologous structure information

While information from homologous structures plays a central role in X-ray structure determination by molecular replacement, such information is rarely used in NMR structure determination because it can be incorrect, both locally and globally, when evolutionary relationships are inferred incorrectly or there has been considerable evolutionary structural divergence. Here we describe a method that allows robust modeling of protein structures of up to 225 residues by combining (1)H(N), (13)C, and (15)N backbone and (13)Cβ chemical shift data, distance restraints derived from homologous structures, and a physically realistic all-atom energy function. Accurate models are distinguished from inaccurate models generated using incorrect sequence alignments by requiring that (i) the all-atom energies of models generated using the restraints are lower than models generated in unrestrained calculations and (ii) the low-energy structures converge to within 2.0 Å backbone rmsd over 75% of the protein. Benchmark calculations on known structures and blind targets show that the method can accurately model protein structures, even with very remote homology information, to a backbone rmsd of 1.2-1.9 Å relative to the conventional determined NMR ensembles and of 0.9-1.6 Å relative to X-ray structures for well-defined regions of the protein structures. This approach facilitates the accurate modeling of protein structures using backbone chemical shift data without need for side-chain resonance assignments and extensive analysis of NOESY cross-peak assignments.

Automated protein NMR structure determination in solution

The main drawback of protein NMR spectroscopy today is still the extensive amount of time required for solving a single structure. The main bottleneck in this respect is the manual evaluation of the experimental spectra. A clear solution to this challenge is the development of automated methods for this purpose. At the current stage of development, this goal has been almost or in a few cases fully reached for favorable cases such as well-behaved, stably folding smaller proteins below the 25 kDa range. For larger and/or more difficult molecules, the input of a human expert is still required. However, even here, automated routines will substantially speed up the structure determination process. In this report, we will summarize recent developments in this field and especially emphasize practical aspects important for a successful automated protein structure determination in solution. An important aspect closely related to structure determination is structure validation. Therefore, we devote a section to automated approaches for this topic.

Protein structure calculation with data imputation: the use of substitute restraints

The amount of experimental restraints e.g., NOEs is often too small for calculating high quality three-dimensional structures by restrained molecular dynamics. Considering this as a typical missing value problem we propose here a model based data imputation technique that should lead to an improved estimation of the correct structure. The novel automated method implemented in AUREMOL makes a more efficient use of the experimental information to obtain NMR structures with higher accuracy. It creates a large set of substitute restraints that are used either alone or together with the experimental restraints. The new approach was successfully tested on three examples: firstly, the Ras-binding domain of Byr2 from Schizosaccharomyces pombe, the mutant HPr (H15A) from Staphylococcus aureus, and a X-ray structure of human ubiquitin. In all three examples, the quality of the resulting final bundles was improved considerably by the use of additional substitute restraints, as assessed quantitatively by the calculation of RMSD values to the "true" structure and NMR R-factors directly calculated from the original NOESY spectra or the published diffraction data.

Ca2+-dependent conformational changes in a C-terminal cytosolic domain of polycystin-2

The PKD1 and PKD2 genes are the genes that are mutated in patients suffering from autosomal dominant polycystic kidney disease. The human PKD2 gene codes for a 968-amino acid long membrane protein called polycystin-2 that represents a cation channel whose activity can be regulated by Ca(2+) ions. By CD, fluorescence, and NMR spectroscopy, we have studied a 117-amino acid-long fragment of the cytoplasmic domain of polycystin-2, polycystin-2-(680-796) that was proposed to contain a Ca(2+)-binding site. NMR structure determination reveals the existence of two Ca(2+)-binding sites in polycystin-2-(680-796) arranged in a typical and an atypical EF-hand motif. In the absence of Ca(2+) the protein forms a dimer that is dissociated by Ca(2+) binding. This dissociation may be related to the Ca(2+) inactivation observed earlier. The calcium affinity of the protein was determined by fluorescence and NMR spectroscopy. At 293 K, the K(D) values for the high and low affinity sites are 55 mum and 179 mum, respectively.