Blind prediction of interfacial water positions in CAPRI.

We report the first assessment of blind predictions of water positions at protein-protein interfaces, performed as part of the critical assessment of predicted interactions (CAPRI) community-wide experiment. Groups submitting docking predictions for the complex of the DNase domain of colicin E2 and Im2 immunity protein (CAPRI Target 47), were invited to predict the positions of interfacial water molecules using the method of their choice. The predictions-20 groups submitted a total of 195 models-were assessed by measuring the recall fraction of water-mediated protein contacts. Of the 176 high- or medium-quality docking models-a very good docking performance per se-only 44% had a recall fraction above 0.3, and a mere 6% above 0.5. The actual water positions were in general predicted to an accuracy level no better than 1.5 Å, and even in good models about half of the contacts represented false positives. This notwithstanding, three hotspot interface water positions were quite well predicted, and so was one of the water positions that is believed to stabilize the loop that confers specificity in these complexes. Overall the best interface water predictions was achieved by groups that also produced high-quality docking models, indicating that accurate modelling of the protein portion is a determinant factor. The use of established molecular mechanics force fields, coupled to sampling and optimization procedures also seemed to confer an advantage. Insights gained from this analysis should help improve the prediction of protein-water interactions and their role in stabilizing protein complexes.

Clustering biomolecular complexes by residue contacts similarity.

Inaccuracies in computational molecular modeling methods are often counterweighed by brute-force generation of a plethora of putative solutions. These are then typically sieved via structural clustering based on similarity measures such as the root mean square deviation (RMSD) of atomic positions. Albeit widely used, these measures suffer from several theoretical and technical limitations (e.g., choice of regions for fitting) that impair their application in multicomponent systems (N > 2), large-scale studies (e.g., interactomes), and other time-critical scenarios. We present here a simple similarity measure for structural clustering based on atomic contacts--the fraction of common contacts--and compare it with the most used similarity measure of the protein docking community--interface backbone RMSD. We show that this method produces very compact clusters in remarkably short time when applied to a collection of binary and multicomponent protein-protein and protein-DNA complexes. Furthermore, it allows easy clustering of similar conformations of multicomponent symmetrical assemblies in which chain permutations can occur. Simple contact-based metrics should be applicable to other structural biology clustering problems, in particular for time-critical or large-scale endeavors.

Protein structure determination from pseudocontact shifts using ROSETTA.

Paramagnetic metal ions generate pseudocontact shifts (PCSs) in nuclear magnetic resonance spectra that are manifested as easily measurable changes in chemical shifts. Metals can be incorporated into proteins through metal binding tags, and PCS data constitute powerful long-range restraints on the positions of nuclear spins relative to the coordinate system of the magnetic susceptibility anisotropy tensor (Δχ-tensor) of the metal ion. We show that three-dimensional structures of proteins can reliably be determined using PCS data from a single metal binding site combined with backbone chemical shifts. The program PCS-ROSETTA automatically determines the Δχ-tensor and metal position from the PCS data during the structure calculations, without any prior knowledge of the protein structure. The program can determine structures accurately for proteins of up to 150 residues, offering a powerful new approach to protein structure determination that relies exclusively on readily measurable backbone chemical shifts and easily discriminates between correctly and incorrectly folded conformations.

Protein-protein HADDocking using exclusively pseudocontact shifts.

In order to enhance the structure determination process of macromolecular assemblies by NMR, we have implemented long-range pseudocontact shift (PCS) restraints into the data-driven protein docking package HADDOCK. We demonstrate the efficiency of the method on a synthetic, yet realistic case based on the lanthanide-labeled N-terminal ε domain of the E. coli DNA polymerase III (ε186) in complex with the HOT domain. Docking from the bound form of the two partners is swiftly executed (interface RMSDs < 1 Å) even with addition of very large amount of noise, while the conformational changes of the free form still present some challenges (interface RMSDs in a 3.1-3.9 Å range for the ten lowest energy complexes). Finally, using exclusively PCS as experimental information, we determine the structure of ε186 in complex with the HOT-homologue θ subunit of the E. coli DNA polymerase III.

A dipicolinic acid tag for rigid lanthanide tagging of proteins and paramagnetic NMR spectroscopy.

A new lanthanide tag was designed for site-specific labeling of proteins with paramagnetic lanthanide ions. The tag, 4-mercaptomethyl-dipicolinic acid, binds lanthanide ions with nanomolar affinity, is readily attached to proteins via a disulfide bond, and avoids the problems of diastereomer formation associated with most of the conventional lanthanide tags. The high lanthanide affinity of the tag opens the possibility to measure residual dipolar couplings in a single sample containing a mixture of paramagnetic and diamagnetic lanthanides. Using the DNA-binding domain of the E. coli arginine repressor as an example, it is demonstrated that the tag allows immobilization of the lanthanide ion in close proximity of the protein by additional coordination of the lanthanide by a carboxyl group of the protein. The close proximity of the lanthanide ion promotes accurate determinations of magnetic susceptibility anisotropy tensors. In addition, the small size of the tag makes it highly suitable for studies of intermolecular interactions.

Numbat: an interactive software tool for fitting Deltachi-tensors to molecular coordinates using pseudocontact shifts.

Pseudocontact shift (PCS) effects induced by a paramagnetic lanthanide bound to a protein have become increasingly popular in NMR spectroscopy as they yield a complementary set of orientational and long-range structural restraints. PCS are a manifestation of the chi-tensor anisotropy, the Deltachi-tensor, which in turn can be determined from the PCS. Once the Deltachi-tensor has been determined, PCS become powerful long-range restraints for the study of protein structure and protein-ligand complexes. Here we present the newly developed package Numbat (New User-friendly Method Built for Automatic Deltachi-Tensor determination). With a Graphical User Interface (GUI) that allows a high degree of interactivity, Numbat is specifically designed for the computation of the complete set of Deltachi-tensor parameters (including shape, location and orientation with respect to the protein) from a set of experimentally measured PCS and the protein structure coordinates. Use of the program for Linux and Windows operating systems is illustrated by building a model of the complex between the E. coli DNA polymerase III subunits epsilon186 and theta using PCS.

Sequence-specific and stereospecific assignment of methyl groups using paramagnetic lanthanides.

Pseudocontact shifts (PCSs) induced by a site-specifically bound paramagnetic lanthanide ion are shown to provide fast access to sequence-specific resonance assignments of methyl groups in proteins of known three-dimensional structure. Stereospecific assignments of Val and Leu methyls are obtained as well as resonance assignments of all other methyls, including Met epsilonCH3 groups. No prior assignments of the diamagnetic protein are required nor are experiments that transfer magnetization between the methyl groups and the protein backbone. Methyl Cz-exchange experiments were designed to provide convenient access to PCS measurements in situations where a paramagnetic lanthanide is in exchange with a diamagnetic lanthanide. In the absence of exchange, simultaneous 13C-HSQC assignments and PCS measurements are delivered by the newly developed program Possum. The approaches are demonstrated with the complex between the N-terminal domain of the subunit epsilon and the subunit theta of the Escherichia coli DNA polymerase III.

Efficient chi-tensor determination and NH assignment of paramagnetic proteins.

Anisotropic magnetic susceptibility tensors chi of paramagnetic metal ions are manifested in pseudocontact shifts, residual dipolar couplings, and other paramagnetic observables that present valuable long-range information for structure determinations of protein-ligand complexes. A program was developed for automatic determination of the chi-tensor anisotropy parameters and amide resonance assignments in proteins labeled with paramagnetic metal ions. The program requires knowledge of the three-dimensional structure of the protein, the backbone resonance assignments of the diamagnetic protein, and a pair of 2D 15N-HSQC or 3D HNCO spectra recorded with and without paramagnetic metal ion. It allows the determination of reliable chi-tensor anisotropy parameters from 2D spectra of uniformly 15N-labeled proteins of fairly high molecular weight. Examples are shown for the 185-residue N-terminal domain of the subunit epsilon from E. coli DNA polymerase III in complex with the subunit theta and La3+ in its diamagnetic and Dy3+, Tb3+, and Er3+ in its paramagnetic form.