Refinement of severely incomplete structures with maximum likelihood in BUSTER-TNT.

BUSTER-TNT is a maximum-likelihood macromolecular refinement package. BUSTER assembles the structural model, scales observed and calculated structure-factor amplitudes and computes the model likelihood, whilst TNT handles the stereochemistry and NCS restraints/constraints and shifts the atomic coordinates, B factors and occupancies. In real space, in addition to the traditional atomic and bulk-solvent models, BUSTER models the parts of the structure for which an atomic model is not yet available ('missing structure') as low-resolution probability distributions for the random positions of the missing atoms. In reciprocal space, the BUSTER structure-factor distribution in the complex plane is a two-dimensional Gaussian centred around the structure factor calculated from the atomic, bulk-solvent and missing-structure models. The errors associated with these three structural components are added to compute the overall spread of the Gaussian. When the atomic model is very incomplete, modelling of the missing structure and the consistency of the BUSTER statistical model help structure building and completion because (i) the accuracy of the overall scale factors is increased, (ii) the bias affecting atomic model refinement is reduced by accounting for some of the scattering from the missing structure, (iii) the addition of a spatial definition to the source of incompleteness improves on traditional Luzzati and sigmaA-based error models and (iv) the program can perform selective density modification in the regions of unbuilt structure alone.

Phasing in the presence of severe site-specific radiation damage through dose-dependent modelling of heavy atoms.

The case of a brominated RNA crystal structure determination in which standard three-wavelength MAD phasing was unsuccessful because of fast X-ray-induced debromination was reinvestigated [Ennifar et al. (2002), Acta Cryst. D58, 1262-1268]. It was found that if the data are kept unmerged and if a dose-stamp is associated with each reflection measurement, dose-dependent occupancies can be refined for the Br atoms. Such a parametrization has been implemented in the macromolecular phasing program SHARP. Refining such dose-dependent occupancies on an unmerged data set gave a dramatic improvement, even for SAD phases from only the first wavelength (peak), and resulted in a good electron-density map after solvent flattening. The adverse effect of radiation damage has been turned into a beneficial one. The crucial difference is made by the use of unmerged data: phasing power is generated through the intensity differences of symmetry-related reflections recorded at different doses, i.e. corresponding to different states of the X-ray-induced debromination. This approach should prove useful in all situations of experimental phasing where site-specific radiation damage occurs unavoidably and undesirably and not only in cases in which radiation damage is purposely being created in order to demonstrate its potential usefulness.

'MAD'ly phasing the extracellular domain of the LDL receptor: a medium-sized protein, large tungsten clusters and multiple non-isomorphous crystals.

The crystal structure of the extracellular domain of the LDL receptor (LDL-R) was determined in a MAD experiment using 12-tungstophosphate clusters as anomalous scatterers. While useful for phasing, the tungsten clusters rendered the crystals radiation-sensitive and non-isomorphous and profoundly altered the diffraction data, causing complications. The work is presented as a case study for phasing a medium-sized protein (700 residues) at low resolution (4 A) with multiple non-isomorphous crystals containing 31 W atoms in the asymmetric unit.

Generation, representation and flow of phase information in structure determination: recent developments in and around SHARP 2.0.

The methods for treating experimental data in the isomorphous replacement and anomalous scattering methods of macromolecular phase determination have undergone considerable evolution since their inception 50 years ago. The successive formulations used are reviewed, from the most simplistic viewpoint to the most advanced, including the exploration of some blind alleys. A new treatment is proposed and demonstrated for the improved encoding and subsequent exploitation of phase information in the complex plane. It is concluded that there is still considerable scope for further improvements in the statistical analysis of phase information, which touch upon numerous fundamental issues related to data processing and experimental design.

A novel combination of two classic catalytic schemes.

The crystal structure of an alkaline Bacillus cellulase catalytic core, from glucoside hydrolase family 5, reveals a novel combination of the catalytic machinery of two classic textbook enzymes. The enzyme has the expected two glutamate residues in close proximity to one another in the active-site that are typical of retaining cellulases. However, the proton donor, glutamate 139 is also unexpectedly a member of a serine-histidine-glutamate catalytic triad, forming a novel combination of catalytic machineries. Structure and sequence analysis of glucoside hydrolase family 5 reveal that the triad is highly conserved, but with variations at the equivalent of the serine position. We speculate that the purpose of this novel catalytic triad is to control the protonation of the acid/base glutamate, facilitating the first step of the catalytic reaction, protonation of the substrate, by the proton donor glutamate. If correct, this will be a novel use for a catalytic triad.

High-resolution experimental phases for tryptophanyl-tRNA synthetase (TrpRS) complexed with tryptophanyl-5'AMP.

Native data, anomalous data at three wavelengths and an independent peak-wavelength data set for SeMet-substituted protein have been collected from cryoprotected crystals of the TrpRS-adenylate product (TAM) complex to a resolution limit of 1.7 A. Independent phase sets were developed using SHARP and improved by solvent flipping with SOLOMON using molecular envelopes derived from experimental densities for, respectively, peak-wavelength SAD data from four different crystals, MAD data and their M(S)IRAS combinations with native data. Hendrickson-Lattman phase-probability coefficients from each phase set were used in BUSTER to drive maximum-likelihood refinements of well defined parts of the previously refined room-temperature 2.9 A structure. Maximum-entropy completion followed by manual rebuilding was then used to generate a model for the missing segments, bound ligand and solvent molecules. Surprisingly, peak-wavelength SAD experiments produced the smallest phase errors relative to the refined structures. Selenomethionylated models deviate from one another by 0.25 A and from the native model by 0.38 A, but all have r.m.s. deviations of approximately 1.0 A from the 2.9 A model. Difference Fourier calculations between amplitudes from the 300 K experiment and the new amplitudes at 100 K using 1.7 A model phases show no significant structural changes arising from temperature variation or addition of cryoprotectant. The main differences between low- and high-resolution structures arise from correcting side-chain rotamers in the core of the protein as well as on the surface. These changes improve various structure-validation criteria.

Crystal structure of MalK, the ATPase subunit of the trehalose/maltose ABC transporter of the archaeon Thermococcus litoralis.

The members of the ABC transporter family transport a wide variety of molecules into or out of cells and cellular compartments. Apart from a translocation pore, each member possesses two similar nucleoside triphosphate-binding subunits or domains in order to couple the energy-providing reaction with transport. In the maltose transporter of several Gram-negative bacteria and the archaeon Thermo coccus litoralis, the nucleoside triphosphate-binding subunit contains a C-terminal regulatory domain. A dimer of the subunit is attached cytoplasmically to the translocation pore. Here we report the crystal structure of this dimer showing two bound pyrophosphate molecules at 1.9 A resolution. The dimer forms by association of the ATPase domains, with the two regulatory domains attached at opposite poles. Significant deviation from 2-fold symmetry is seen at the interface of the dimer and in the regions corresponding to those residues known to be in contact with the translocation pore. The structure and its relationship to function are discussed in the light of known mutations from the homologous Escherichia coli and Salmonella typhimurium proteins.

Structure of the 30S ribosomal subunit.

Genetic information encoded in messenger RNA is translated into protein by the ribosome, which is a large nucleoprotein complex comprising two subunits, denoted 30S and 50S in bacteria. Here we report the crystal structure of the 30S subunit from Thermus thermophilus, refined to 3 A resolution. The final atomic model rationalizes over four decades of biochemical data on the ribosome, and provides a wealth of information about RNA and protein structure, protein-RNA interactions and ribosome assembly. It is also a structural basis for analysis of the functions of the 30S subunit, such as decoding, and for understanding the action of antibiotics. The structure will facilitate the interpretation in molecular terms of lower resolution structural data on several functional states of the ribosome from electron microscopy and crystallography.

Modelling prior distributions of atoms for macromolecular refinement and completion.

Until modelling is complete, macromolecular structures are refined in the absence of a model for some of the atoms in the crystal. Techniques for defining positional probability distributions of atoms, and using them to model the missing part of a macromolecular crystal structure and the bulk solvent, are described. The starting information may consist of either a tentative structural model for the missing atoms or an electron-density map. During structure completion and refinement, the use of probability distributions enables the retention of low-resolution phase information while avoiding premature commitment to uncertain higher resolution features. Homographic exponential modelling is proposed as a flexible, compact and robust parametrization that proves to be superior to a traditional Fourier expansion in approximating a model protein envelope. The homographic exponential model also has potential applications to ab initio phasing of Fourier amplitudes associated with macromolecular envelopes.

Applications of single-wavelength anomalous dispersion at high and atomic resolution.

Two examples of the application of single-wavelength anomalous dispersion (SAD) in macromolecular structure determination are described, both using the statistical phasing program SHARP. For the holmium-substituted calcium-binding protein psoriasin (22.7 kDa), a set of accurate phases has been obtained to a resolution of 1.05 A without recourse to an atomic model of the molecule. The accuracy of the phases resulted in an electron-density map of a quality comparable to sigma(A)-weighted 2mF(o) - DF(c) maps derived from the final model refined with SHELX97. Comparison of the refined and SAD electron-density maps showed significant discrepancies resulting from the iterative refinement in reciprocal space. Additionally, it is shown that the structure of psoriasin can be determined from native data extending to 2.0 A alone by exploiting the minute anomalous signal from a bound zinc ion.

The structure of adrenodoxin reductase of mitochondrial P450 systems: electron transfer for steroid biosynthesis.

Adrenodoxin reductase is a monomeric 51 kDa flavoenzyme that is involved in the biosynthesis of all steroid hormones. The structure of the native bovine enzyme was determined at 2.8 A resolution, and the structure of the respective recombinant enzyme at 1.7 A resolution. Adrenodoxin reductase receives a two-electron package from NADPH and converts it to two single electrons that are transferred via adrenodoxin to all mitochondrial cytochromes P 450. The structure suggests how the observed flavin semiquinone is stabilized. A striking feature is the asymmetric charge distribution, which most likely controls the approach of the electron carrier adrenodoxin. A model for the interaction is proposed. Adrenodoxin reductase shows clear sequence homology to half a dozen proteins identified in genome analysis projects, but neither sequence nor structural homology to established, functionally related electron transferases. Yet, the structure revealed a relationship to the disulfide oxidoreductases, permitting the assignment of the NADP-binding site.

Locating proper non-crystallographic symmetry in low-resolution electron-density maps with the program GETAX.

Non-crystallographic symmetry averaging for improving and extending an initial set of phases can be crucial at an early stage of a protein structure analysis. A method is described which detects the position of a proper rotation axis in a surprisingly poor electron-density map and is fast enough to run through a large number of axis orientations. It uses a simple multimer mask to define the searching unit, which is then shifted through the whole unit cell looking for the position with the highest correlation coefficient between the interrelated parts. Appropriate weighting and averaging enhances the signal-to-noise ratio. Examples of the application of this algorithm are given. The use of the local rotation axis for phasing is commented on. A search of the Protein Data Bank showed that 27% of the unique crystal forms contain proper local n-fold axes, which could have been located with the presented method.

Chaperone-assisted expression of authentic bovine adrenodoxin reductase in Escherichia coli.

Adrenodoxin reductase is an essential component of the mitochondrial monooxygenase systems that are involved in the synthesis of steroid hormones and related compounds. After removing by mutagenesis a secondary ribosome binding site and an mRNA loop formed between the gene and the vector, large amounts of the enzyme could be produced in Escherichia coli by coexpression with the HSP60-chaperone system. The purified protein was homogeneous enough for reproducible crystallization. The crystals diffracted X-rays isotropically beyond 1.7 A resolution permitting a structure analysis.

The structure of a trimeric archaeal adenylate kinase.

The adenylate kinase from the hyperthermophilic archaean species Sulfolobus acidocaldarius has been cloned, expressed in Escherichia coli, purified and crystallized. The crystal structure was elucidated by multiple isomorphous replacement and non-crystallographic density averaging. The structure was refined at 2.6 A (1 A=0.1 nm) resolution. The enzyme is trimeric, in contrast to previous solution measurements that suggested a dimeric structure, and in contrast to the vast majority of adenylate kinases, which are monomeric. In large parts of each subunit the chain fold resembles the known enzyme structure from eubacteria and eukaryotes although the sequence homology is negligible. Since the asymmetric unit contains two trimers with and without bound AMP at the AMP sites and with an ADP at one of the six ATP sites, the analysis shows the enzyme in several states. The conformational differences between these states resemble those of other adenylate kinases. Because of sequence homology, the structure presented provides a good model for the methanococcal adenylate kinases.

Movie of the structural changes during a catalytic cycle of nucleoside monophosphate kinases.

BACKGROUND: There are 17 crystal structures of nucleoside monophosphate kinases known. As expected for kinases, they show large conformational changes upon binding of substrates. These are concentrated in two chain segments, or domains, of 30 and 38 residues that are involved in binding of the substrates N1TP and N2MP (nucleoside tri- and monophosphates with bases N1 and N2), respectively. RESULTS: After aligning the 17 structures on the main parts of their polypeptide chains, two domains in various conformational states were revealed. These states were caused by bound substrate (or analogues) and by crystal-packing forces, and ranged between a 'closed' conformation and a less well defined 'open' conformation. The structures were visually sorted yielding an approximately evenly spaced series of domain states that outlines the closing motions when the substrates bind. The packing forces in the crystals are weak, leaving the natural domain trajectories essentially intact. Packing is necessary, however, to produce stable intermediates. The ordered experimental structures were then recorded as still pictures of a movie and animated to represent the motions of the molecule during a catalytic cycle. The motions were smoothed out by adding interpolated structures to the observed ones. The resulting movies are available through the World Wide Web (http:@bio5.chemie.uni-freiburg.de/ak movie.html). CONCLUSIONS: Given the proliferating number of homologous proteins known to exist in different conformational states, it is becoming possible to outline the motions of chain segments and combine them into a movie, which can then represent protein action much more effectively than static pictures alone are able to do.