MolProbity: More and better reference data for improved all-atom structure validation.

This paper describes the current update on macromolecular model validation services that are provided at the MolProbity website, emphasizing changes and additions since the previous review in 2010. There have been many infrastructure improvements, including rewrite of previous Java utilities to now use existing or newly written Python utilities in the open-source CCTBX portion of the Phenix software system. This improves long-term maintainability and enhances the thorough integration of MolProbity-style validation within Phenix. There is now a complete MolProbity mirror site at http://molprobity.manchester.ac.uk. GitHub serves our open-source code, reference datasets, and the resulting multi-dimensional distributions that define most validation criteria. Coordinate output after Asn/Gln/His "flip" correction is now more idealized, since the post-refinement step has apparently often been skipped in the past. Two distinct sets of heavy-atom-to-hydrogen distances and accompanying van der Waals radii have been researched and improved in accuracy, one for the electron-cloud-center positions suitable for X-ray crystallography and one for nuclear positions. New validations include messages at input about problem-causing format irregularities, updates of Ramachandran and rotamer criteria from the million quality-filtered residues in a new reference dataset, the CaBLAM Cα-CO virtual-angle analysis of backbone and secondary structure for cryoEM or low-resolution X-ray, and flagging of the very rare cis-nonProline and twisted peptides which have recently been greatly overused. Due to wide application of MolProbity validation and corrections by the research community, in Phenix, and at the worldwide Protein Data Bank, newly deposited structures have continued to improve greatly as measured by MolProbity's unique all-atom clashscore.

A new generation of crystallographic validation tools for the protein data bank.

This report presents the conclusions of the X-ray Validation Task Force of the worldwide Protein Data Bank (PDB). The PDB has expanded massively since current criteria for validation of deposited structures were adopted, allowing a much more sophisticated understanding of all the components of macromolecular crystals. The size of the PDB creates new opportunities to validate structures by comparison with the existing database, and the now-mandatory deposition of structure factors creates new opportunities to validate the underlying diffraction data. These developments highlighted the need for a new assessment of validation criteria. The Task Force recommends that a small set of validation data be presented in an easily understood format, relative to both the full PDB and the applicable resolution class, with greater detail available to interested users. Most importantly, we recommend that referees and editors judging the quality of structural experiments have access to a concise summary of well-established quality indicators.

MolProbity: all-atom structure validation for macromolecular crystallography.

MolProbity is a structure-validation web service that provides broad-spectrum solidly based evaluation of model quality at both the global and local levels for both proteins and nucleic acids. It relies heavily on the power and sensitivity provided by optimized hydrogen placement and all-atom contact analysis, complemented by updated versions of covalent-geometry and torsion-angle criteria. Some of the local corrections can be performed automatically in MolProbity and all of the diagnostics are presented in chart and graphical forms that help guide manual rebuilding. X-ray crystallography provides a wealth of biologically important molecular data in the form of atomic three-dimensional structures of proteins, nucleic acids and increasingly large complexes in multiple forms and states. Advances in automation, in everything from crystallization to data collection to phasing to model building to refinement, have made solving a structure using crystallography easier than ever. However, despite these improvements, local errors that can affect biological interpretation are widespread at low resolution and even high-resolution structures nearly all contain at least a few local errors such as Ramachandran outliers, flipped branched protein side chains and incorrect sugar puckers. It is critical both for the crystallographer and for the end user that there are easy and reliable methods to diagnose and correct these sorts of errors in structures. MolProbity is the authors' contribution to helping solve this problem and this article reviews its general capabilities, reports on recent enhancements and usage, and presents evidence that the resulting improvements are now beneficially affecting the global database.

The other 90% of the protein: assessment beyond the Calphas for CASP8 template-based and high-accuracy models.

For template-based modeling in the CASP8 Critical Assessment of Techniques for Protein Structure Prediction, this work develops and applies six new full-model metrics. They are designed to complement and add value to the traditional template-based assessment by the global distance test (GDT) and related scores (based on multiple superpositions of Calpha atoms between target structure and predictions labeled "Model 1"). The new metrics evaluate each predictor group on each target, using all atoms of their best model with above-average GDT. Two metrics evaluate how "protein-like" the predicted model is: the MolProbity score used for validating experimental structures, and a mainchain reality score using all-atom steric clashes, bond length and angle outliers, and backbone dihedrals. Four other new metrics evaluate match of model to target for mainchain and sidechain hydrogen bonds, sidechain end positioning, and sidechain rotamers. Group-average Z-score across the six full-model measures is averaged with group-average GDT Z-score to produce the overall ranking for full-model, high-accuracy performance. Separate assessments are reported for specific aspects of predictor-group performance, such as robustness of approximately correct template or fold identification, and self-scoring ability at identifying the best of their models. Fold identification is distinct from but correlated with group-average GDT Z-score if target difficulty is taken into account, whereas self-scoring is done best by servers and is uncorrelated with GDT performance. Outstanding individual models on specific targets are identified and discussed. Predictor groups excelled at different aspects, highlighting the diversity of current methodologies. However, good full-model scores correlate robustly with high Calpha accuracy.

MolProbity: all-atom contacts and structure validation for proteins and nucleic acids.

MolProbity is a general-purpose web server offering quality validation for 3D structures of proteins, nucleic acids and complexes. It provides detailed all-atom contact analysis of any steric problems within the molecules as well as updated dihedral-angle diagnostics, and it can calculate and display the H-bond and van der Waals contacts in the interfaces between components. An integral step in the process is the addition and full optimization of all hydrogen atoms, both polar and nonpolar. New analysis functions have been added for RNA, for interfaces, and for NMR ensembles. Additionally, both the web site and major component programs have been rewritten to improve speed, convenience, clarity and integration with other resources. MolProbity results are reported in multiple forms: as overall numeric scores, as lists or charts of local problems, as downloadable PDB and graphics files, and most notably as informative, manipulable 3D kinemage graphics shown online in the KiNG viewer. This service is available free to all users at http://molprobity.biochem.duke.edu.

The backrub motion: how protein backbone shrugs when a sidechain dances.

Surprisingly, the frozen structures from ultra-high-resolution protein crystallography reveal a prevalent, but subtle, mode of local backbone motion coupled to much larger, two-state changes of sidechain conformation. This "backrub" motion provides an influential and common type of local plasticity in protein backbone. Concerted reorientation of two adjacent peptides swings the central sidechain perpendicular to the chain direction, changing accessible sidechain conformations while leaving flanking structure undisturbed. Alternate conformations in sub-1 angstroms crystal structures show backrub motions for two-thirds of the significant Cbeta shifts and 3% of the total residues in these proteins (126/3882), accompanied by two-state changes in sidechain rotamer. The Backrub modeling tool is effective in crystallographic rebuilding. For homology modeling or protein redesign, backrubs can provide realistic, small perturbations to rigid backbones. For large sidechain changes in protein dynamics or for single mutations, backrubs allow backbone accommodation while maintaining H bonds and ideal geometry.

Away from the edge II: in-house Se-SAS phasing with chromium radiation.

Recently, the demands of high-throughput macromolecular crystallography have driven continuous improvements in phasing methods, data-collection protocols and many other technologies. Single-wavelength anomalous scattering (SAS) phasing with chromium X-ray radiation opens a new possibility for phasing a protein with data collected in-house and has led to several successful examples of de novo structure solution using only weak anomalous scatterers such as sulfur. To further reduce data-collection time and make SAS phasing more robust, it is natural to combine selenomethionine-derivatized protein (SeMet protein) with Cr Kalpha radiation to take advantage of the larger anomalous scattering signal from selenium (f'' = 2.28 e(-)) compared with sulfur (f'' = 1.14 e(-)). As reported herein, the crystal structure of a putative chorismate mutase from Clostridium thermocellum was determined using Se-SAS with Cr Kalpha radiation. Each protein molecule contains eight selenomethionine residues in 148 amino-acid residues, providing a calculated Bijvoet ratio of about 3.5% at the Cr Kalpha wavelength. A single data set to 2.2 A resolution with approximately ninefold redundancy was collected using an imaging-plate detector coupled with a Cr source. Structure solution, refinement and deposition to the Protein Data Bank were performed within 9 h of the availability of the scaled diffraction data. The procedure used here is applicable to many other proteins and promises to become a routine pathway for in-house high-throughput crystallography.

A test of enhancing model accuracy in high-throughput crystallography.

The high throughput of structure determination pipelines relies on increased automation and, consequently, a reduction of time spent on interactive quality control. In order to meet and exceed current standards in model accuracy, new approaches are needed for the facile identification and correction of model errors during refinement. One such approach is provided by the validation and structure-improvement tools of the MOL: PROBITY: web service. To test their effectiveness in high-throughput mode, a large subset of the crystal structures from the SouthEast Collaboratory for Structural Genomics (SECSG) has used protocols based on the MOL: PROBITY: tools. Comparison of 29 working-set and 19 control-set SECSG structures shows that working-set outlier scores for updated Ramachandran-plot, sidechain rotamer, and all-atom steric criteria have been improved by factors of 5- to 10-fold (relative to the control set or to a Protein Data Bank sample), while quality of covalent geometry, R(work), R(free), electron density and difference density are maintained or improved. Some parts of this correction process are already fully automated; other parts involve manual rebuilding of conformations flagged by the tests as trapped in the wrong local minimum, often altering features of functional significance. The ease and effectiveness of this technique shows that macromolecular crystal structures from either traditional or high-throughput determinations can feasibly reach a new level of excellence in conformational accuracy and reliability.

The cis-Pro touch-turn: a rare motif preferred at functional sites.

A new motif of three-dimensional (3D) protein structure is described, called the cis-Pro touch-turn. In this four-residue, three-peptide motif, the central peptide is cis. Residue 2, which precedes the proline, has phi, psi values either in the "prePro region" of the Ramachandran plot near -130 degrees, 75 degrees or in the Lalpha region near +60 degrees, +60 degrees. The Calpha(1)-Calpha(4) distance is 4-5 A and the two flanking peptides lie parallel to one another, making van der Waals contact rather than a hydrogen bond. Apparently, this arrangement is locally unfavorable and therefore rare, usually occurring only if needed for biological function. Of the 12 examples in a 500-protein database, cis-Pro touch-turns are found at the catalytic sites of pectate lyase, Ni-Fe hydrogenase, glucoamylase, xylanase, and opine dehydrogenase and at the primary binding sites of ribonuclease H, type I DNA polymerase, ribotoxin, and phage gene 3 protein. In each of these protein families, the touch-turns serve different roles; their functional importance is supported by conservation and mutagenesis data. In analyzing the conservation patterns of these 3D motifs, new methods for in-depth quality evaluation of the structural bioinformatic data are employed to distinguish between significant exceptions and errors

RNA backbone is rotameric.

Despite the importance of local structural detail to a mechanistic understanding of RNA catalysis and binding functions, RNA backbone conformation has been quite recalcitrant to analysis. There are too many variable torsion angles per residue, and their raw empirical distributions are poorly clustered. This study applies quality-filtering techniques (using resolution, crystallographic B factor, and all-atom steric clashes) to the backbone torsion angle distributions from an 8,636-residue RNA database. With noise levels greatly reduced, clear signal appears for the underlying angle preferences. Half-residue torsion angle distributions for alpha-beta-gamma and for delta-epsilon-zeta are plotted and contoured in 3D; each shows about a dozen distinct peaks, which can then be combined in pairs to define complete RNA backbone conformers. Traditional nucleic acid residues are defined from phosphate to phosphate, but here we use a base-to-base (or sugar-to-sugar) division into "suites" to parse the RNA backbone repeats, both because most backbone steric clashes are within suites and because the relationship of successive bases is both reliably determined and conformationally important. A suite conformer has seven variables, with sugar pucker specified at both ends. Potential suite conformers were omitted if not represented by at least a small cluster of convincing data points after application of quality filters. The final result is a small library of 42 RNA backbone conformers, which should provide valid conformations for nearly all RNA backbone encountered in experimental structures.

Structure validation by Calpha geometry: phi,psi and Cbeta deviation.

Geometrical validation around the Calpha is described, with a new Cbeta measure and updated Ramachandran plot. Deviation of the observed Cbeta atom from ideal position provides a single measure encapsulating the major structure-validation information contained in bond angle distortions. Cbeta deviation is sensitive to incompatibilities between sidechain and backbone caused by misfit conformations or inappropriate refinement restraints. A new phi,psi plot using density-dependent smoothing for 81,234 non-Gly, non-Pro, and non-prePro residues with B < 30 from 500 high-resolution proteins shows sharp boundaries at critical edges and clear delineation between large empty areas and regions that are allowed but disfavored. One such region is the gamma-turn conformation near +75 degrees,-60 degrees, counted as forbidden by common structure-validation programs; however, it occurs in well-ordered parts of good structures, it is overrepresented near functional sites, and strain is partly compensated by the gamma-turn H-bond. Favored and allowed phi,psi regions are also defined for Pro, pre-Pro, and Gly (important because Gly phi,psi angles are more permissive but less accurately determined). Details of these accurate empirical distributions are poorly predicted by previous theoretical calculations, including a region left of alpha-helix, which rates as favorable in energy yet rarely occurs. A proposed factor explaining this discrepancy is that crowding of the two-peptide NHs permits donating only a single H-bond. New calculations by Hu et al. [Proteins 2002 (this issue)] for Ala and Gly dipeptides, using mixed quantum mechanics and molecular mechanics, fit our nonrepetitive data in excellent detail. To run our geometrical evaluations on a user-uploaded file, see MOLPROBITY (http://kinemage.biochem.duke.edu) or RAMPAGE (http://www-cryst.bioc.cam.ac.uk/rampage).