Protonation-state determination in proteins using high-resolution X-ray crystallography: effects of resolution and completeness.

A bond-distance analysis has been undertaken to determine the protonation states of ionizable amino acids in trypsin, subtilisin and lysozyme. The diffraction resolutions were 1.2 Šfor trypsin (97% complete, 12% H-atom visibility at 2.5σ), 1.26 Šfor subtilisin (100% complete, 11% H-atom visibility at 2.5σ) and 0.65 Šfor lysozyme (PDB entry 2vb1; 98% complete, 30% H-atom visibility at 3σ). These studies provide a wide diffraction resolution range for assessment. The bond-length e.s.d.s obtained are as small as 0.008 Šand thus provide an exceptional opportunity for bond-length analyses. The results indicate that useful information can be obtained from diffraction data at around 1.2-1.3 Šresolution and that minor increases in resolution can have significant effects on reducing the associated bond-length standard deviations. The protonation states in histidine residues were also considered; however, owing to the smaller differences between the protonated and deprotonated forms it is much more difficult to infer the protonation states of these residues. Not even the 0.65 Šresolution lysozyme structure provided the necessary accuracy to determine the protonation states of histidine.

IUCrData launches Raw Data Letters.

A new category of articles - Raw Data Letters - is introduced to IUCrData.

Determination of zinc incorporation in the Zn-substituted gallophosphate ZnULM-5 by multiple wavelength anomalous dispersion techniques.

The location of isomorphously substituted zinc over eight crystallographically different gallium sites has been determined in a single-crystal study of the gallophosphate ZnULM-5, Ga((16-x))Zn(x)(PO(4))(14)(HPO(4))(2)(OH)(2)F(7), [H(3)N{CH(2)}(6)NH(3)](4), 6H(2)O, in an 11 wavelength experiment, using data from Station 9.8, SRS Daresbury. The measurement of datasets around the K edges of both Ga and Zn, as well as two reference datasets away from each absorption edge, was utilized to selectively exploit dispersive differences of each metal atom type in turn, which allowed the major sites of Zn incorporation to be identified as the metal 1 and 3 sites, M1 and M3. The preferential substitution of Zn at these sites probably arises because they are located in double four-ring (D4R) building units which can relax to accommodate the incorporation of hetero atoms. As the crystal is non-centrosymmetric, with space group P2(1)2(1)2, it was also possible to use anomalous differences to corroborate the results obtained from the dispersive differences. These results were obtained firstly from difference Fourier maps, calculated using a phase set from the refined structure from data measured at the Zr K edge. Also, refined dispersive and anomalous occupancies, on an absolute scale, could be obtained using the program MLPHARE, allowing estimates for the Zn incorporation of approximately 22 and 18 at. % at the M1 and M3 sites to be obtained. In addition, f' and f'' values for Ga and Zn at each wavelength could be estimated both from MLPHARE results, and by refinement in JANA2006. The fully quantitative determinations of the dispersive and anomalous coefficients for Ga and Zn at each wavelength, as well as metal atom occupancies over the eight metal atom sites made use of the CCP4's MLPHARE program as well as SHELXL and JANA2006. The results by these methods agree closely, and JANA2006 allowed the ready determination of standard uncertainties on the occupancy parameters, which were for M1 and M3, 20.6 (3) and 17.2 (3) at %, respectively.

An investigation into structural changes due to deuteration.

Perdeuteration of proteins is becoming more commonplace and the assumption is in general that deuteration does not affect protein structure. In this work, the effect of deuteration on structure is examined by data mining, largely of the Cambridge Structural Database but also of the Inorganic Crystal Structure Database, for deuterated and hydrogenated pairs of small-molecule structures analysed by neutron and X-ray crystallography. Differences between these small-molecule structures have been calculated and the results thus far follow the initial assumption. However, functional changes are known, e.g. D(2)O is toxic to living systems but H(2)O is not, kinetics change, small pH to pD changes occur, proteins stiffen in D(2)O and ferroelectrics alter their properties.

The emergent rhizosphere: imaging the development of the porous architecture at the root-soil interface.

The rhizosphere is the zone of soil influenced by a plant root and is critical for plant health and nutrient acquisition. All below ground resources must pass through this dynamic zone prior to their capture by plant roots. However, researching the undisturbed rhizosphere has proved very challenging. Here we compare the temporal changes to the intact rhizosphere pore structure during the emergence of a developing root system in different soils. High resolution X-ray Computed Tomography (CT) was used to quantify the impact of root development on soil structural change, at scales relevant to individual micro-pores and aggregates (µm). A comparison of micro-scale structural evolution in homogenously packed soils highlighted the impacts of a penetrating root system in changing the surrounding porous architecture and morphology. Results indicate the structural zone of influence of a root can be more localised than previously reported (µm scale rather than mm scale). With time, growing roots significantly alter the soil physical environment in their immediate vicinity through reducing root-soil contact and crucially increasing porosity at the root-soil interface and not the converse as has often been postulated. This 'rhizosphere pore structure' and its impact on associated dynamics are discussed.

Self-association of a DNA loop creates a quadruplex: crystal structure of d(GCATGCT) at 1.8 A resolution.

BACKGROUND: The flexibility of DNA enables it to adopt three interconvertible types of duplex termed the A-, B- and Z-forms. It can also produce hairpin loops, triplex structures and guanine-rich quadruplex structures. Conformational flexibility assists in the tight packaging of DNA, for example in chromosomes. This is important given the large quantity of genetic information that must be packaged efficiently. Moreover, the ability of DNA to specifically self-associate or interact with complementary sequences is fundamental to many biological processes. Structural studies provide information about DNA conformation and DNA-DNA interactions and suggest features that might be relevant to how the molecule performs its biological role. RESULTS: We have characterized the structure of a synthetic heptanucleotide that folds into a novel loop structure. The loop is stabilized by association with a cation, by intra-strand hydrogen bonds between guanine and cytosine that are distinct from the normal Watson-Crick hydrogen bonds, and by van der Waals interactions. Two loops associate through the formation of four G.C pairs that exhibit pronounced base-stacking interactions. The formation of a symmetric A.A base pair further stabilizes loop dimerization. Stacking of the A.A pair on a symmetry-related A.A pairing assists the formation of a four-stranded assembly. A T.T pairing is also observed between symmetry-related loops. CONCLUSIONS: This analysis provides a rare example of an experimentally determined non-duplex DNA structure. It provides conformational detail relevant to the tight packaging or folding of a DNA strand and illustrates how a cation might modulate phosphate-phosphate repulsion in a tightly packed structure. The observation of base quartets involving G.C base pairs suggests a further structure to be considered in DNA-DNA interactions. The structure also provides detailed geometries for A.A and T.T base pairs.

The 1.2 A resolution structure of the Con A-dimannose complex.

The complex between concanavalin A (Con A) and alpha1-2 mannobiose (mannose alpha1-2 mannose) has been refined to 1.2 A resolution. This is the highest resolution structure reported for any sugar-lectin complex. As the native structure of Con A to 0.94 A resolution is already in the database, this gives us a unique opportunity to examine sugar-protein binding at high resolution. These data have allowed us to model a number of hydrogen atoms involved in the binding of the sugar to Con A, using the difference density map to place the hydrogen atoms. This map reveals the presence of the protonated form of Asp208 involved in binding. Asp208 is not protonated in the 0.94 A native structure. Our results clearly show that this residue is protonated and hydrogen bonds to the sugar. The structure accounts for the higher affinity of the alpha1-2 linked sugar when compared to other disaccharides. This structure identifies different interactions to those predicted by previous modelling studies. We believe that the additional data presented here will enable significant improvements to be made to the sugar-protein modelling algorithms.

Properties of a new crystal form of the complex of concanavalin A with methyl alpha-D-glucopyranoside.

The complex of concanavalin A with methyl alpha-D-glucopyranoside crystallizes as regular rhombic dodecahedra containing 35% protein by weight. The crystal is of space group I23 with a = 167.8 A (1 A = 0.1 nm) and contains one concanavalin A dimer per asymmetric unit. It diffracts to a resolution of 1.9 A and is suitable for crystallographic investigation of the structure of the saccharide-binding site.

Initiating a crystallographic study of trypanothione reductase.

We have obtained well-ordered single crystals of the flavoenzyme trypanothione reductase from Crithidia fasciculata. The crystals are tetragonal rods with unit cell dimensions a = 128.6 A, c = 92.5 A. The diffraction pattern corresponds to a primitive lattice. Laue class 4/m. Diffraction to better than 2.4 A has been recorded at the Daresbury Synchrotron. The accurate elucidation of the three-dimensional structure of this enzyme is required to support the rational design of compounds active against a variety of tropical diseases caused by trypanosomal parasites.

Active site of trypanothione reductase. A target for rational drug design.

The X-ray crystal structure of the enzyme trypanothione reductase, isolated from the trypanosomatid organism Crithidia fasciculata, has been solved by molecular replacement. The search model was the crystal structure of human glutathione reductase that shares approximately 40% sequence identity. The trypanosomal enzyme crystallizes in the tetragonal space group P4(1) with unit cell lengths of a = 128.9 A and c = 92.3 A. The asymmetric unit consists of a homodimer of approximate molecular mass 108 kDa. We present the structural detail of the active site as derived from the crystallographic model obtained at an intermediate stage of the analysis using diffraction data to 2.8 A resolution with an R-factor of 23.2%. This model has root-mean-square deviations from ideal geometry of 0.026 A for bond lengths and 4.7 degrees for bond angles. The trypanosomid enzyme assumes a similar biological function to glutathione reductase and, although similar in topology to human glutathione reductase, has an enlarged active site and a number of amino acid differences, steric and electrostatic, which allows it to process only the unique substrate trypanothione and not glutathione. This protein represents a prime target for chemotherapy of several debilitating tropical diseases caused by protozoan parasites belonging to the genera Trypanosoma and Leishmania. The structural differences between the parasite and host enzymes and their substrates thus provides a rational basis for the design of new drugs active against trypanosomes. In addition, our model explains the results of site-directed mutagenesis experiments, carried out on recombinant trypanothione reductase and glutathione reductases, designed by consideration of the crystal structure of human glutathione reductase.

An investigation into the protonation states of the C1 domain of cardiac myosin-binding protein C.

Myosin-binding protein C (MyBP-C) is a myofibril-associated protein found in cardiac and skeletal muscle. The cardiac isoform (cMyBP-C) is subject to reversible phosphorylation and the surface-charge state of the protein is of keen interest with regard to understanding the inter-protein interactions that are implicated in its function. Diffraction data from the C1 domain of cMyBP-C were extended to 1.30 A resolution, where the of the diffraction data crosses 2.0, using intense synchrotron radiation. The protonation-state determinations were not above 2sigma (the best was 1.81sigma) and therefore an extrapolation is given, based on 100% data completeness and the average DPI, that a 3sigma determination could be possible if X-ray data could be measured to 1.02 A resolution. This might be possible via improved crystallization or multiple sample evaluation, e.g. using robotics or a yet more intense/collimated X-ray beam or combinations thereof. An alternative would be neutron protein crystallography at 2 A resolution, where it is estimated that for the unit-cell volume of the cMyBP-C C1 domain crystal a crystal volume of 0.10 mm3 would be needed with fully deuterated protein on LADI III. These efforts would optimally be combined in a joint X-ray and neutron model refinement.

Science experiments via telepresence at a synchrotron radiation source facility.

Station 9.8 is one of the most oversubscribed and high-throughput stations at the Synchrotron Radiation Source, Daresbury, whereby awarded experimental time is limited, data collections last normally no longer than an hour, user changeover is normally every 24 h, and familiarity with the station systems can be low. Therefore time lost owing to technical failures on the station has a dramatic impact on productivity. To provide 24 h support, the application of a turnkey communication system has been implemented, and is described along with additional applications including its use for inter-continental classroom instruction, user training and remote participation.

The determination of protonation states in proteins.

The protonation states of aspartic acids and glutamic acids as well as histidine are investigated in four X-ray cases: Ni,Ca concanavalin A at 0.94 A, a thrombin-hirugen binary complex at 1.26 A resolution and two thrombin-hirugen-inhibitor ternary complexes at 1.32 and 1.39 A resolution. The truncation of the Ni,Ca concanavalin A data at various test resolutions between 0.94 and 1.50 A provided a test comparator for the ;unknown' thrombin-hirugen carboxylate bond lengths. The protonation states of aspartic acids and glutamic acids can be determined (on the basis of convincing evidence) even to the modest resolution of 1.20 A as exemplified by our X-ray crystal structure refinements of Ni and Mn concanavalin A and also as indicated in the 1.26 A structure of thrombin, both of which are reported here. The protonation-state indication of an Asp or a Glu is valid provided that the following criteria are met (in order of importance). (i) The acidic residue must have a single occupancy. (ii) Anisotropic refinement at a minimum diffraction resolution of 1.20 A (X-ray data-to-parameter ratio of approximately 3.5:1) is required. (iii) Both of the bond lengths must agree with the expectation (i.e. dictionary values), thus allowing some relaxation of the bond-distance standard uncertainties required to approximately 0.025 A for a '3sigma' determination or approximately 0.04 A for a '2sigma' determination, although some variation of the expected bond-distance values must be allowed according to the microenvironment of the hydrogen of interest. (iv) Although the F(o) - F(c) map peaks are most likely to be unreliable at the resolution range around 1.20 A, if admitted as evidence the peak at the hydrogen position must be greater than or equal to 2.5 sigma and in the correct geometry. (v) The atomic B factors need to be less than 10 A(2) for bond-length differentiation; furthermore, the C=O bond can also be expected to be observed with continuous 2F(o) - F(c) electron density and the C-OH bond with discontinuous electron density provided that the atomic B factors are less than approximately 20 A(2) and the contour level is increased. The final decisive option is to carry out more than one experiment, e.g. multiple X-ray crystallography experiments and ideally neutron crystallography. The complementary technique of neutron protein crystallography has provided evidence of the protonation states of histidine and acidic residues in concanavalin A and also the correct orientations of asparagine and glutamine side chains. Again, the truncation of the neutron data at various test resolutions between 2.5 and 3.0 A, even 3.25 and 3.75 A resolution, examines the limits of the neutron probe. These various studies indicate a widening of the scope of both X-ray and neutron probes in certain circumstances to elucidate the protonation states in proteins.

Time differences between friedel reflections: accuracy of crystal setting and requirements on beam stability.

In synchrotron radiation data collections where the wavelength is carefully set to optimize f''-derived crystallography intensity differences (Friedel pairs), careful alignment of the crystal is useful to minimize the time differences of stimulation of the reflections in the pairs. This paper quantifies these time differences as a function of crystal misorientation, with typical parameters, using the angular velocity of the crystal. Likewise, the time spent in the diffraction condition is also calculated via the angular reflecting range for a common synchrotron beam geometry. These times offer the user direct insight into the time-dependent aspects of the diffraction measurements. This therefore allows the optimum conditions to be set up so as to extract as accurate anomalous differences as possible in the context of synchrotron radiation beam stability and lifetime.

The properties of (2Fo - Fc) and (Fo - Fc) electron-density maps at medium-to-high resolutions.

This paper reports on the efficacy of (F(o) - F(c)) versus (2F(o) - F(c)) electron-density maps at 3.2 A resolution. Firstly, a study is reported of a simple truncation at 2.3 and 3.2 A of the 1.6 A resolution crystal structure of concanavalin A at room temperature [Emmerich et al. (1994), Acta Cryst. D50, 749-756] with 149 known bound water molecules. Secondly, the concanavalin A 1.6 A resolution model was re-refined but with the data truncated to 3.2 A. In a similar evaluation, these procedures were repeated for the apocrustacyanin A1 cryotemperature 1.4 A resolution model [Cianci et al. (2001), Acta Cryst. D57, 1219-1229]. Maps at 1.4, 2.3 and 3.2 A resolutions were first generated and the structure was then re-refined at 3.2 A and additionally at 2.3 A resolution. The results on concanavalin A show that the number of bound water molecules that are resolved decreases by two thirds from 1.6 to 3.2 A, but that key structural waters, for example at the transition metal and the calcium ion, are still resolved in the (F(o) - F(c)) map but not in the (2F(o) - F(c)) map. For apocrustacyanin A1, the results with these two difference maps were less clear-cut. Two key structural bound waters (w93 and w105) were selected that had been previously identified in beta-crustacyanin [Cianci et al. (2002), Proc. Natl Acad. Sci. USA, 99, 9795-9800] in protein-carotenoid interactions. The behaviour of w93 is similar to that of concanavalin A key waters, but that of w105 is not. These behaviours were therefore explored in finer resolution increments, namely 2.9, 2.7 and 2.5 A. Finally, further tests on "real" data sets for peanut lectin and concanavalin A at medium resolution confirm these map properties, namely that an (F(o) - F(c)) difference electron-density map is more effective than a (2F(o) - F(c)) map in showing bound water structure at lower resolutions ( approximately 3.2 A). This result is important since a growing number of protein crystal structure studies are concerned with multi-macromolecular complexes and are at such resolutions. Details of the bound solvent can still be revealed at 3.2 A via the (F(o) - F(c)) map calculation. The physical basis of the limitation of the (2F(o) - F(c)) map presumably lies in the series-termination error effect on such a map involving the first negative ripple from the protein atom to which a bound water oxygen is hydrogen bonded, sufficiently cancelling its peak. In addition, re-refinements at 3.2 A show distances that can agree with known values but B values that do not agree with known values.

Integration of macromolecular diffraction data using radial basis function networks.

This paper presents a novel approach for intensity calculation of X-ray diffraction spots based on a two-stage radial basis function (RBF) network. The first stage uses pre-determined reference profiles from a database as basis functions in order to locate the diffraction spots and identify any overlapping regions. The second-stage RBF network employs narrow basis functions capable of local modifications of the reference profiles leading to a more accurate observed diffraction spot approximation and therefore accurate determination of spot positions and integrated intensities.