The sensitivity of the commonly used progressive multiple sequence alignment method has been greatly improved for the alignment of divergent protein sequences. Firstly, individual weights are assigned to each sequence in a partial alignment in order to down-weight near-duplicate sequences and up-weight the most divergent ones. Secondly, amino acid substitution matrices are varied at different alignment stages according to the divergence of the sequences to be aligned. Thirdly, residue-specific gap penalties and locally reduced gap penalties in hydrophilic regions encourage new gaps in potential loop regions rather than regular secondary structure. Fourthly, positions in early alignments where gaps have been opened receive locally reduced gap penalties to encourage the opening up of new gaps at these positions. These modifications are incorporated into a new program, CLUSTAL W which is freely available.

/pdf/clustal-w-improving-the-sensitivity-of-progressive-multiple-2na2s9sq2h.pdf

Clustal w: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice

The design, implementation, and capabilities of an extensible visualization system, UCSF Chimera, are discussed. Chimera is segmented into a core that provides basic services and visualization, and extensions that provide most higher level functionality. This architecture ensures that the extension mechanism satisfies the demands of outside developers who wish to incorporate new features. Two unusual extensions are presented: Multiscale, which adds the ability to visualize large-scale molecular assemblies such as viral coats, and Collaboratory, which allows researchers to share a Chimera session interactively despite being at separate locales. Other extensions include Multalign Viewer, for showing multiple sequence alignments and associated structures; ViewDock, for screening docked ligand orientations; Movie, for replaying molecular dynamics trajectories; and Volume Viewer, for display and analysis of volumetric data. A discussion of the usage of Chimera in real-world situations is given, along with anticipated future directions. Chimera includes full user documentation, is free to academic and nonprofit users, and is available for Microsoft Windows, Linux, Apple Mac OS X, SGI IRIX, and HP Tru64 Unix from http://www.cgl.ucsf.edu/chimera/.

/pdf/ucsf-chimera-a-visualization-system-for-exploratory-research-59xu54fhxg.pdf

UCSF Chimera--a visualization system for exploratory research and analysis.

The Protein Data Bank (PDB; http://www.rcsb.org/pdb/ ) is the single worldwide archive of structural data of biological macromolecules. This paper describes the goals of the PDB, the systems in place for data deposition and access, how to obtain further information, and near-term plans for the future development of the resource.

/pdf/the-protein-data-bank-4lbxxcio6v.pdf

The Protein Data Bank

Publisher Summary X-ray data can be collected with zero-, one-, and two-dimensional detectors, zero-dimensional (single counter) being the simplest and two-dimensional the most efficient in terms of measuring diffracted X-rays in all directions. To analyze the single-crystal diffraction data collected with these detectors, several computer programs have been developed. Two-dimensional detectors and related software are now predominantly used to measure and integrate diffraction from single crystals of biological macromolecules. Macromolecular crystallography is an iterative process. To monitor the progress, the HKL package provides two tools: (1) statistics, both weighted (χ2) and unweighted (R-merge), where the Bayesian reasoning and multicomponent error model helps obtain proper error estimates and (2) visualization of the process, which helps an operator to confirm that the process of data reduction, including the resulting statistics, is correct and allows the evaluation of the problems for which there are no good statistical criteria. Visualization also provides confidence that the point of diminishing returns in data collection and reduction has been reached. At that point, the effort should be directed to solving the structure. The methods presented in the chapter have been applied to solve a large variety of problems, from inorganic molecules with 5 A unit cell to rotavirus of 700 A diameters crystallized in 700 × 1000 × 1400 A cell.

Processing of X-ray diffraction data collected in oscillation mode

An N⋅log(N) method for evaluating electrostatic energies and forces of large periodic systems is presented. The method is based on interpolation of the reciprocal space Ewald sums and evaluation of the resulting convolutions using fast Fourier transforms. Timings and accuracies are presented for three large crystalline ionic systems.

Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems

For a successful analysis of the relation between amino acid sequence and protein structure, an unambiguous and physically meaningful definition of secondary structure is essential. We have developed a set of simple and physically motivated criteria for secondary structure, programmed as a pattern-recognition process of hydrogen-bonded and geometrical features extracted from x-ray coordinates. Cooperative secondary structure is recognized as repeats of the elementary hydrogen-bonding patterns “turn” and “bridge.” Repeating turns are “helices,” repeating bridges are “ladders,” connected ladders are “sheets.” Geometric structure is defined in terms of the concepts torsion and curvature of differential geometry. Local chain “chirality” is the torsional handedness of four consecutive Cα positions and is positive for right-handed helices and negative for ideal twisted β-sheets. Curved pieces are defined as “bends.” Solvent “exposure” is given as the number of water molecules in possible contact with a residue. The end result is a compilation of the primary structure, including SS bonds, secondary structure, and solvent exposure of 62 different globular proteins. The presentation is in linear form: strip graphs for an overall view and strip tables for the details of each of 10.925 residues. The dictionary is also available in computer-readable form for protein structure prediction work.

Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features

An algorithm has been developed for the automatic interpretation of a given set of observed reciprocal-lattice points. It extracts a reduced cell and assigns indices to each reflection by a graph-theoretical implementation of the local indexing method. All possible symmetries of the observed lattice compatible with the metric of the reduced cell are recognized and reported, together with the unit-cell constants and the linear index transformation relating the conventional to the reduced cell. This algorithm has been incorporated into the program XDS [Kabsch (1988). J. Appl. Cryst. 21, 916–924], which is now able to process single-crystal area-detector data without prior knowledge of the symmetry and the unit-cell constants.

Automatic processing of rotation diffraction data from crystals of initially unknown symmetry and cell constants

A simple procedure is derived which determines a best rotation of a given vector set into a second vector set by minimizing the weighted sum of squared deviations. The method is generalized for any given metric constraint on the transformation.

A solution for the best rotation to relate two sets of vectors

Important steps in the processing of rotation data are described that are common to most software packages. These programs differ in the details and in the methods implemented to carry out the tasks. Here, the working principles underlying the data-reduction package XDS are explained, including the new features of automatic determination of spot size and reflecting range, recognition and assignment of crystal sym­metry and a highly efficient algorithm for the determination of correction/scaling factors.

/pdf/integration-scaling-space-group-assignment-and-post-1xors09pmu.pdf

Integration, scaling, space‐group assignment and post‐refinement

The atomic models of the complex between rabbit skeletal muscle actin and bovine pancreatic deoxyribonuclease I both in the ATP and ADP forms have been determined by X-ray analysis at an effective resolution of 2.8 A and 3A, respectively. The two structures are very similar. The actin molecule consists of two domains which can be further subdivided into two subdomains. ADP or ATP is located in the cleft between the domains with a calcium ion bound to the beta- or beta- and gamma-phosphates, respectively. The motif of a five-stranded beta sheet consisting of a beta meander and a right handed beta alpha beta unit appears in each domain suggesting that gene duplication might have occurred. These sheets have the same topology as that found in hexokinase.

Wolfgang Kabsch

Papers

Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features

Automatic processing of rotation diffraction data from crystals of initially unknown symmetry and cell constants

A solution for the best rotation to relate two sets of vectors

Integration, scaling, space‐group assignment and post‐refinement

Atomic structure of the actin:DNase I complex.