research papers
Lowresolution REFMAC5
tools in^{a}Structural Studies Division, MRC Laboratory of Molecular Biology, Cambridge CB2 0QH, England
^{*}Correspondence email: garib@mrclmb.cam.ac.uk
Two aspects of lowresolution macromolecular crystal structure analysis are considered: (i) the use of reference structures and structural units for provision of structural prior information and (ii) map sharpening in the presence of noise and the effects of Fourier series termination. The generation of interatomic distance restraints by ProSMART and their subsequent application in REFMAC5 is described. It is shown that the use of such external structural information can enhance the reliability of derived atomic models and stabilize The problem of map sharpening is considered as an inverse deblurring problem and is solved using Tikhonov regularizers. It is demonstrated that this type of map sharpening can automatically produce a map with more structural features whilst maintaining connectivity. Tests show that both of these directions are promising, although more work needs to be performed in order to further exploit structural information and to address the problem of reliable electrondensity calculation.
Keywords: lowresolution ; REFMAC5.
1. Introduction
Heterogeneous organization of molecules in the crystal lattice can lead to the diffraction data being of poor quality. Such heterogeneities may arise from effects such as crystal mosaicity, molecule/chain flexibility and localized disorder. This results in weak diffraction intensities, causing the data to be collected using a lower resolution threshold and thus with low information content. This behaviour is often observed for large complexes. However, the structures of the individual components of a complex may have been independently determined at higher resolution. Such information might then be used to aid the
of the lower resolution structure.There are other factors that can reduce the information content of macromolecular crystallographic (MX) data, thus reducing effective resolution. These include crystalgrowth peculiarities such as ^{1/3} = 1.26. Therefore, in the limit, the quality of the electrondensity map in the presence of perfect hemihedral at 2 Å would correspond to that at 2.52 Å in the singlecrystal case. The of models against data from twinned crystals is now routine (Murshudov et al., 2011; Adams et al., 2010; Sheldrick, 2008). However, statistics after against such data should be interpreted with care (Murshudov, 2011). It is important to remember that R factors and other overall statistics are dependent on the statistical properties of the data and therefore comparison of R factors from different crystals may give the wrong impression about the comparative quality of the models.
and OD order–disorder. In these cases, although the nominal resolution may be high, not all of the observations are independent. For example, in the case of perfect hemihedral the number of independent observations is decreased by a factor of two, corresponding to a resolution reduction by a factor of 2There are many problems that need to be tackled in order to make lowresolution structure analysis routine, two of which are considered in this paper.
Many tools have been developed to aid crystallographic et al., 2002). Structural information may be utilized in various forms, such as secondarystructure restraints, homologous reference structures and homology models, by various modern software packages including REFMAC5 (Murshudov et al., 1997, 2011) from CCP4 (Winn et al., 2011), BUSTERTNT (Blanc et al., 2004), phenix.refine (Adams et al., 2010; Afonine et al., 2005), SHELX (Sheldrick, 2008) and CNS (Schröder et al., 2010; Brünger et al., 1998).
at medium and higher resolutions over the past few decades. One of the current challenges is to develop complementary approaches for dealing with cases where only lowresolution data are available (lower than around 3 Å). One of the sources of available information is the threedimensional structures of macromolecules deposited in the Protein Data Bank (BermanThe concept of calculating an electrondensity map showing more features, e.g. side chains, has been proposed by many authors. Notably, Brunger and coworkers (Brunger et al., 2009; DeLaBarre & Brunger, 2006) have suggested a procedure known in the field of image processing (Gonzalez & Woods, 2002) as inverse filtering. However, it is known that such filters can amplify noise, thus masking out real signal. Unfortunately, the electron density always contains noise, which stems from several sources.

In this paper, we first describe the use of external structural information, specifically interatomic distance restraints obtained from reference homologous structures and general fragmentbased restraints (including quasisecondarystructure restraints, in particular to helical conformations). We then describe anisotropic regularized map sharpening. For each method, we provide examples of usage.
2. Structural restraints
2.1. Application of external structural restraints in crystallographic refinement
Information from external sources can be incorporated during
using a Bayesian framework where distribution of interatomic distances serves as prior knowledge. Thus, restraints generated using external structural information should help the macromolecule under to adopt a conformation that is more consistent with previous observations. This is similar to the use of geometry terms in which helps local structure adopt chemically reasonable conformations.The minus log posterior distribution target used in REFMAC5 (Murshudov et al., 1997) may be expressed
where f_{geom} and f_{xray} are the contributions of geometry terms (−prior distribution) and experimental data (−log likelihood) and w weights their relative contributions. The geometry component is a linear combination of various factors (effectively equivalent to the assumption that these contributors are independent), including any available external structural information,
where f_{ext}(, κ) is the unweighted contribution of an external interatomic distance restraint (d, r, σ) ∈ R to the target function, where R is the list of external restraints and the function f_{ext}(, κ) depends on the normalized residual and parameter κ. The parameter w_{ext} adjusts the weight of the external restraints relative to the other geometry components and f_{other} represents the contribution of all other prior information (Murshudov et al., 2011). An interatomic distance restraint comprises the current distance d between two atomic positions, the objective value r and the σ. The mechanism used for application of external restraints in REFMAC5 is described by Mooij et al. (2009).
Here, we stipulate that the objective value r of an external restraint should be lower than some threshold r_{max}, so that only reasonably shortrange restraints are utilized, thus allowing resistance to differences in global conformation. In REFMAC5, the Geman–McClure robust estimation function (Geman & McClure, 1987) is used for external restraints,
This function, which is equivalent to least squares for κ = 0, helps to reduce the influence of outliers and sensitivity to conformational changes.
Various criteria have been used for optimization of the Xray weight w, notably R_{free} (Brünger, 1997) and −LL_{free} (Bricogne, 1997; Tickle, 2007). Similarly, the appropriate selection of the external weight w_{ext} is not obvious. One potential optimization criterion might be to minimize R_{free}. However, it should be acknowledged that this may reduce the usefulness of R_{free} as an independent indicator of quality. Therefore, the weight for external structural information w_{ext} requires careful consideration.
2.1.1. Selection of external structural information
External structures should be selected on the basis of their reliability and similarity to the current model. For example, suitable reference structures may include sequenceidentical, homologous or structurally similar models solved at a higher resolution or generically derived structural information from nonhomologous sources (e.g. secondarystructure restraints obtained from an ideal αhelix).
The use of external restraints may in some cases be justified by any resultant increase in the reliability of atomic positions. However, it should be acknowledged that such an approach introduces bias; the influence of such bias may result in the model adopting a conformation that is less consistent with the observed data. The use of external restraints might make a particular model adopt a conformation very similar to a highresolution homologue, assuming it is appropriate to do so, and ideally result in improved
statistics.We suggest that external restraints should only be used if the benefits of any improvement in reliability are deemed to outweigh the negative effects. Indeed, this may well be the case for data of poor quality collected at low resolution. For example, r_{max}; see below).
of a model might cause some regions of very poor electron density to adopt an incorrect conformation. Increasing the weight of geometry terms may help the structure to adopt a more chemically reasonable configuration, but the region may still be incorrectly modelled owing to the effect of the misleading density; geometry restraints operate at a very high level of structural resolution. However, external restraints can operate at a lower level of structural resolution, as desired (by increasing the maximum restraint distance2.2. External restraint generation
Here, we shall refer to the chain that is to be refined as the target chain and to the chain that is to be used to generate the restraints as the external or reference chain.
External restraints for use in REFMAC5 may be generated using the ProSMART tool developed by Nicholls (2011). Amongst various other functionalities, ProSMART can generate restraints on interatomic distances utilizing structural information. Whilst reference structures would generally be nearidentical in sequence, the approach allows the alignment of, and subsequent restraint generation using, any reference chain(s). However, it is not implied that there would be any utility in using external restraints based on dissimilar structures; a high degree of local structural conservation would generally be required for the successful application of external restraints. In general, we assume that the target and external reference structures are sufficiently similar, although such decisions should ultimately be made by the user.
byThe alignment approach adopted by ProSMART is independent of global conformation, instead being concerned with the net conservation of local structure at a chosen level of structural resolution. Indeed, the restraints generated by ProSMART allow great global flexibility rather than rigidly pulling the target structure towards the same global conformation. Consequently, it is not necessary for the external reference chain to adopt the same global conformation as the target, e.g. structures in different bound states can be used. It is, however, necessary for local structure to be sufficiently well conserved along the chain so that the effect on is positive.
2.2.1. General approach
Suppose that we want to generate an interatomic distance restraint between two atoms in the target structure given knowledge of their positions and thus the distance d between them. Given an external reference structure and a residue alignment between the target and reference structures, it is possible to find the distance r between the corresponding atoms in the reference structure. The distance r is the objective value of the restraint.
If the target and external chains share a high degree of structural similarity, then we might expect d to be approximately equal to r, with some error. Consequently, the restraint distances r, with appropriate distributional assumptions, can be used as prior information during crystallographic Since we want to maintain a degree of global conformational independence between the target and reference structures, it is undesirable to generate restraints between atoms that are far apart. Therefore, restraints are only generated whose objective values are less than some threshold r_{max}. This parameter represents the structural resolution of the restraints; lower thresholds allow better conformational independence, whilst higher thresholds provide more information about the surrounding structural environment.
The adopted procedure of external restraint generation thus involves the identification of lists of corresponding intrachain atom pairs in both the target and reference structures, filtering these lists in order to identify only those atom pairs suitable for restraint generation (based on interatomic distance criteria), the identification of corresponding atom pairs between the target and reference structures and finally estimation of restraint distributions.
2.2.2. Identification of atom pairs to be restrained
Knowledge of an alignment between residues allows the direct inference of an atomic correspondence between target and reference structures. Such a correspondence may include both mainchain and sidechain atoms (providing the aligned amino acids are the same) or only mainchain N, C^{α}, C and O atoms (allowing mainchain structural restraints to be generated even for residues of different aminoacid type). The alignment may also be filtered according to conservation of local mainchain and/or sidechain structure in an attempt to avoid the generation of potentially unsuitable restraints.
Given an alignment of atoms, it is then necessary to identify a list of sufficiently close atom pairs independently for each of the two structures. Various methods of nearneighbour searching have been developed. Here, in order to efficiently reduce the search space, we use a cell technique (Bentley, 1975) previously used in biology (Levinthal, 1966), which involves the uniform partitioning (voxelization) of space into cubic cells with edge length r_{max} (the value of r_{max} is chosen to be 1.5 times greater for the target structure than for the reference). This approach allows the efficient identification of all atoms with positions x_{i} and x_{j} such that their interatomic distance satisfies the criteron r_{min} ≤ x_{i} − x_{j} ≤ r_{max}.
Using the achieved atomic correspondence, we may then calculate the list of all pairs of corresponding atom pairs, only considering those identified as being sufficiently close. The quantities of interest directly follow, namely the interatomic distance d_{ij} = x_{i}^{target} − x_{j}^{target} between atoms i and j in the target structure and the distance r_{ij} = x_{i}^{ref} − x_{j}^{ref} between corresponding atoms in the reference structure.
Finally, distances between atom pairs that are already tightly restrained by standard geometry terms are removed from the list of external restraints. In particular, we remove any short restraints separated by only one or two chemical bonds (see Figs. 1 and 2).
2.3. estimation of structural restraint distributions
Removal of distances restrained by standard geometry is vital for estimation of restraint distributions. It is reasonable to assume that the variability of longer range restraints would be very different to that of short restraints separated by only few bonds (see Figs. 1 and 2).
2.3.1. Form of the restraint distributions
Suppose the distributions of the positions of two atoms in the target structure are x_{1} ≃ N(c_{1}, σ_{1}^{2}) and x_{2} ≃ N(c_{2}, σ_{2}^{2}), where c_{i} is the coordinate corresponding to atom i. Since we are generally interested in lowresolution structures, we assume spherical normality; the variance terms are scalar to emphasize this point. Note that B factors are closely related to the variabilities of these distributions, which are usually chosen to be isotropic for lowresolution structures.
The distribution of vectors from the first atom to the second is given by
If the atoms are close then their positions are likely to be positively correlated, which will reduce the variability of the distance between them. Conversely, if the atoms are far apart then it is reasonable to surmise that their positions would be more independent and thus the variability of their interatomic distance would be larger. This is supported by Fig. 2, which demonstrates lower variability for atom pairs that are closer together.
Given that, under assumption of independence of atomic positions, [ (Δx_{i}/σ)^{2}]^{1/2} follows a noncentral χ distribution with three with noncentrality parameter {[E(Δx_{i})/σ]^{2}}^{1/2}, we deduce that the interatomic distance D = [(Δx_{i})^{2}]^{1/2} is related to the noncentral χ distribution; specifically, Dσ^{−1} ≃ χ′_{3}, where σ^{2} = var(Δx). However, for consistency between ProSMART and REFMAC5 we use the assumption of a normal distribution of external distances,
which constitutes the restraint to be used in μ = r, the distance between the corresponding atoms in the reference structure. Appropriate selection of the standard deviation σ is less obvious; currently used approaches are described below.
Given knowledge of external structural information, we estimate the mean as2.3.2. Estimating restraint standard deviations (SDs)
The observed distribution P(dr) of interatomic distances in the target structure given corresponding distances r in the reference structure may be used for selection of restraint SDs. For example, estimation of uniform SDs would allow restraints to be automatically weighted according to the overall agreement between interatomic distances in the two structures. In this trivial case, all SDs are estimated using
This would result in higher SDs (lower weights) being assigned to all external restraints when the reference structure is less similar to the target. Owing to the distancedependence of the variability of d − r, using a higher distance threshold r_{max} would also result in higher SDs. It follows that the restraints would naturally be weighted down if the target and reference structures exhibit conformational differences.
Another choice would be to allow the SDs to increase with the mean in order to account for the distancedependence of the observed distribution of restraints. This would allow restraints with small objective values (r) to have higher weights, whilst naturally weighting down the longerrange restraints. For example, the restraint variance could be allowed to increase linearly with restraint distance, that is
where the parameters k = [k_{1}, k_{2}] depend on the particular chain pair. This could be justified by the observation that any signalling causing correlation in atomic position would generally become weaker as the restraint distance increases. Furthermore, peculiar behaviour may be observed when there are multiple rigid substructures (e.g. domains) present, the effect of which would be exacerbated when the maximum restraint threshold r_{max} is large. The presence of multiple domains, or indeed any deterministic conformational changes, would tend to cause a systematic increase in observed restraint error for higher distances r.
The B values or reliability of atomic position if available could be used in the derivation of restraint SDs.
approach is general enough to allow estimation of parameters from other more complicated functional forms. For example, information such asAlternatively, attempts to sensibly estimate restraint SDs may be bypassed, instead allowing the weight of external restraints to be controlled using only the weighting terms in the
program.2.3.3. estimation of parameters
Given a functional form for the restraint variances σ^{2}(k), we then use estimation to optimize the parameters k. The optimization problem amounts to searching for parameter values k such that the constraints −∂log(L)/∂k_{j} = 0 are satisfied for all j within some acceptable error margin, so that the likelihood function L(k) is sufficiently maximized. In general, the probability density function of D is given by
Parameters of SDs are estimated by minimizing the −log likelihood, the gradient of which is given by
Note that other distributional forms could be considered and handled using this method. Minimization is performed using a quasiNewton method, in which an approximation of the Hessian matrix is updated after each procedural iteration. Specifically, we use the BFGS formula for updating the inverse Hessian approximation and a linesearch algorithm for selecting trial parameter values as described by Nocedal & Wright (1999).
2.4. Fragmentbased restraints
Further to using a reference structure, ProSMART is able to generate restraints based on individual structural units. This functionality may have broad application, including in the generation of restraints for secondarystructural elements. In particular, external restraints may be generated using an nresidue fragment representing an `ideal' αhelix, which may be used to keep helical structures intact. Such restraints might be used when a suitable reference structure is not available or when the reference chain is itself not sufficiently well refined. However, the suitability of other general insequence fragmentbased restraints, such as for βstrands, is less obvious owing to their comparatively high degrees of flexibility and the fact that hydrogen bonding occurs between, not within, βstrands. Another potential application would be when it is desired for a particular region to adopt a known conformation (e.g. if a specific small portion of conserved structure is found between the target and a reference chain); the suitability of such an approach would have to be carefully considered for the particular case.
Since aligned fragments may overlap (e.g. consecutive helical fragments), it is possible for a particular atom pair to be restrained to several atom pairs in the reference fragment. For example, in a helical fragment the distances between mainchain atoms in residues i and j may be very similar to those in residues i + 1 and j + 1. In such cases, restraints for a target atom pair in a helix might be generated using corresponding atoms from residues i and j or those from residues i + 1 and j + 1. Consequently, it is necessary to decide which residues to use for restraint generation. More generally, any restraint between atoms from residues i and j may result from several fragment alignments. Specifically, the reference fragment which has residue range [1, n] may be aligned to any of the residue ranges [j − n + 1, j], …, [i, i + n − 1] in the target structure, whilst still implying correspondences for residues i and j (where j − n < i < j). Therefore, ignoring heterogeneities and boundary conditions, there may be up to i − j + n potential alignments of residues i and j with some residues in the reference fragment.
The list of potential residue correspondences is reduced by fragment score criteria, since we only want to generate fragmentbased restraints for regions of structure sufficiently similar to the reference fragment; only configurations with associated Procrustes dissimilarity (local r.m.s.d.) scores below some threshold are included. Of the remaining potential residuepair alignments, if any, the one with the most favourable associated fragment Procrustes score is selected for restraint generation.
2.5. Examples of usage
Here, we present examples of the rerefinement of models previously deposited in the Protein Data Bank (Berman et al., 2002). Where appropriate, models were rerefined using 30 iterations of by REFMAC5 v.5.7 using external structural restraints generated by ProSMART; distancedependent SDs were used for all external restraints. Illustrations of protein structures were generated using CCP4mg (McNicholas et al., 2011), with comparative structural analyses performed by ProSMART. Model validation (geometry and contact analysis) was performed using the MolProbity server (Chen et al., 2010; Davis et al., 2007). Graphs were generated using R (R Development Core Team, 2011).
The consideration of altering some major parameters (w_{ext}, κ, r_{max}) demonstrates typical behaviour that might be expected in simple cases. Here, parameters were selected by optimizing R_{free} by trial and error, although it should be noted that other criteria may be chosen (e.g. −LL_{free}).
In our examples,
is automated, largely using default settings, and no attempt is made to achieve `good' final models. quality of local regions is not considered given the present purpose, in which we are interested only in the qualitative effect of external restraints on global statistics. Indeed, the examples shown here are neither further improved nor manually inspected following better models, and thus statistics, would be achieved by optimizing other parameters and by subsequent iterations of manual and automated model building and refinement.It should be noted that the examples of the rerefinement of deposited models presented here may not represent typical application, since external restraints may also be applied during earlier stages of the
process in order to help models to adopt more reliable conformations.2.5.1. Application of external structural restraints
Rerefinement of the 3.4 Å resolution structure with PDB code 2jha (Sutton et al., 2007) was attempted using both mainchain and sidechain external restraints from the sequenceidentical 2.5 Å resolution structure 2jhp . Both the target (2jha ) and external reference (2jhp ) structures comprise one chain crystallized in P6_{5}22. As can be seen in Fig. 3, they share very similar global conformations. However, the backbone trace is not identical. At a local level, differences in backbone conformation can be detected in a few regions and many residues have different sidechain conformations. For the purpose of this example, it is unknown/unassumed whether these differences are real, i.e. a consequence of suboptimal (of target or reference structures), or actual conformational differences. In practice, the reference structure would be manually inspected for poorly built/refined regions of the model.
As can be seen in Table 1, the unoptimized rerefinement of 2jha without external restraints results in a greatly reduced R factor. However, ΔR becomes large, suggesting that the rerefined model suffers from overfitting. The use of external restraints from 2jhp results in a considerable decrease in R_{free} and −LL_{free} and also a reduced ΔR, suggesting that the external structural information stabilizes and increases model reliability in this case. The parameters w_{ext} = 7.6, κ = 1.0 and r_{max} = 4.2 were chosen so as to minimize R_{free}. Since alone are not sufficient to unambiguously deduce model improvement, we also consider statistics provided by the MolProbity validation server. with external restraints results in improved geometry and a reduced number of steric clashes, suggesting an improved model, in agreement with our previous assertions based on statistics.

In order to understand the influence of external restraints, it is of interest to perform a structural comparison of the rerefined model with the target and reference models. It is important to confirm that the target structure is not restrained too tightly to the conformation of the reference; regions of structure that are actually different between the target and reference structures should be allowed to differ. In this particular case, the global r.m.s.d. of mainchain atoms between the externally restrained rerefined structure and the target 2jha (and the reference 2jhp ) is 0.442 (and 0.375), indicating that the rerefined structure has diverged from the conformation of 2jha whilst not being restrained too tightly to 2jhp . Furthermore, it is possible to see in Fig. 4 that some regions of backbone structure have remained close to the conformation of the original structure 2jha .
Robust estimation using the Geman–McClure function for the contribution of external restraints to the likelihood function, which reduces the effect of outliers, helps to ensure that regions of structure that correspond to actual differences between the target and reference structures are not restrained too tightly. As illustrated in Fig. 4, there are multiple residues whose side chains (and also some regions of backbone) have not been pulled into the of 2jhp (these residues are coloured red). These differences may represent actual differences between the crystals or errors in one of the models; manual inspection of such regions may reveal opportunities for further model improvement.
In contrast, the side chains of many residues have adopted very similar conformations to those in the reference structure (coloured yellow), whilst departing from the conformations of those in the original structure 2jha and the model rerefined without external restraints. This may be reasonable: the data corresponding to the target structure may be so poor that the bias introduced by the reference model appropriately increases model reliability (as in this example). In other cases where the restraints are deemed too tight it may be appropriate to alter relevant parameters, exclude certain residues from having external restraints or restrain only mainchain atoms.
It is important to acknowledge that the selection of certain parameters, most notably the external restraints weight w_{ext}, the Geman–McClure parameter κ and the maximum restraint length r_{max}, can have a very large effect on Furthermore, the appropriate choice of parameter values seems to be different for different cases, meaning that trial and error is currently required in order to produce reasonable results when using external structural restraints. Unfortunately, these parameters are highly correlated. The appropriate choice of parameters may depend on various factors such as data quality, resolution, other parameters, the suitability of the external reference structure and whether or not sidechain atoms are to be restrained. Consequently, careful consideration should be made in the application of external restraints.
Fig. 5 shows the range of achieved using different values of the parameters w_{ext}, κ and r_{max}, demonstrating that the choice of these parameters is important for the successful application of external restraints.
When the influence of external restraints is weak (small w_{ext} or large κ) the external restraints have little or no positive effect on If the external restraints do not contain sufficient positive information then the may worsen (R and R_{free} will rise). This negative effect can be countered by introducing restraints that have a high useful information content (requiring a suitable choice of reference structure), as is observed for intermediate parameter values in this case.
Conversely, selecting very high weights (high w_{ext} or low κ) has the effect of restraining the model too tightly to the reference structure, with the restraints behaving more like constraints. This often results in greatly increased values of R and R_{free} depending on the structural similarity of the target and reference models.
Refinement statistics arising from a variety of maximum restraint lengths r_{max} are shown in Fig. 5(c). Using a low r_{max} results in relatively few restraints being generated, thus having little effect on As r_{max} increases more restraints are generated and the external restraints have a greater impact on Note that longer restraints are less tolerant of conformational change, influencing tighter globally rigid structural agreement with the reference structure. Therefore, using larger values of r_{max} has strong negative effects in cases where conformational changes are present between target and reference structures. However, the effect of this is not dramatic in this example, since the structures are well conserved at the global level.
2.5.2. Fragmentbased αhelix restraints
We now consider the rerefinement of the 3.3 Å resolution model of human haemoglobin with PDB code 1ydz (Kavanaugh et al., 2005). The use of external restraints is demonstrated using fragment restraints from an ideal helix (using a fragment length of five residues and a Procrustes dissimilarity score threshold of 0.3 Å) and also external restraints from a nearsequenceidentical 1.07 Å resolution reference structure 2w72 (Savino et al., 2009). Both the target (1ydz ; P2_{1}2_{1}2) and the reference (2w72 ; P2_{1}) structures comprise four subunits.
As can be seen in Table 2, rerefinement of 1ydz without external restraints results in improved particularly when using local NCS restraints. Application of helix restraints results in a further decrease in R_{free} and ΔR (although an increased −LL_{free}; this may or may not have been the case were −LL_{free} used as the optimization criteria). More substantial improvements in are achieved by using external restraints from 2w72 both when using only mainchain restraints and when using both mainchain and sidechain restraints. Such situations, where the use of external structural information from both mainchain and sidechain atoms appears to improve the model, arise owing to local relative atomic positions being highly conserved between the reference and (actual) target structures. Thus, in such cases information contained in the reference structure can successfully be used to improve the model of the lowresolution structure. However, in other cases where the reference structure is less similar to the target the use of external sidechain restraints may not be appropriate.

It may be argued that both the homologous structures and the αhelical fragment appear to provide useful information that has a positive effect on crystallographic in this case. These assertions are also supported by general improvements in MolProbity validation scores, although one would anticipate subsequent manual to result in substantial further improvements.
Evidence suggests that external restraints on the homologous structure have a greater positive impact than restraints on the presumed helical conformation. This makes sense, since the helical restraints generically pull local backbone structure towards the αhelix attractor, whilst external restraints from 2w72 contain information specific to the particular protein class. The use of a library of helical fragments may result in further improvements. We conclude that helix restraints may be useful in some cases, particularly when appropriate highresolution reference structures are not available. More generally, given an appropriate fragment library, it may be possible to generate restraints to attractors in fragmentconformation space using this method.
Fig. 6 demonstrates the effect of altering the external restraint parameters w_{ext}, κ and r_{max} from their optimal values (according to R_{free} minimization criteria). These parameters are found to be optimal at different values depending on whether helical restraints, mainchain restraints or mainchain and sidechain restraints are applied. We conclude that the naive application of external restraints without parameter optimization can result in reduced model quality, even in cases where the reference structure is an appropriate choice. The appropriate choice of parameters may be very different in different cases. The successful application of external restraints requires the suitable selection of reference structure(s), atom pairs, estimation of SDs and parameters.
3. Map sharpening
The mapsharpening problem can be written in the general form
where ρ_{0} is the underlying signal that we would like to observe (actual electron density), ρ is the observed signal (model of electron density from observation), g is a process through which blurring operates on the signal, k is a blurring function that changes the signal (ρ_{0}) before observation is carried out and n is noise. However, this formulation is too general to be practical. In order to make the problem manageable, we must make assumptions regarding the functional forms of g and k and assume a model for the noise n. Therefore, for simplicity, we assume that noise is additive and the blurring function is linear,
If there were no noise then the problem would be a linear equation. This problem is illposed, especially when k is near singular, i.e. small perturbations in input parameters may cause large variations in output. For example, the effects of small noise addition, an incorrectly defined blurring function or Fourier series termination may result in an uninterpretable `deblurred' electrondensity map. It should be noted that in crystallography we always deal with limited noisy data and that Fourier series termination is always present. Even if there were no noise and we had knowledge of the exact blurring function k(x, y), solving (11) would still not be straightforward. The numbers of equations and parameters to be estimated are equal to the number of grid points in the electron density, which can be very large.
The problem becomes manageable, whilst not completely reflecting reality, when we make the further assumption that the blurring function is independent of position. This simplification essentially means that the whole content of the k(x, y) = k(x − y, 0). Using the notation k(x) = k(x, 0), (11) becomes
oscillates as a unit with no rotational component, resulting in the blurring function having the propertySince the problem is illposed, we can approach its solution utilizing ideas from the field of regularization (Tikhonov & Arsenin, 1977). Under the assumption of white noise, our illposed problem may be replaced by the minimization problem
where . denotes the L_{2} norm, f is a regularization function and γ is a regularization parameter to be selected. Usually, regularizers are chosen so that the resultant function obeys certain conditions. For example purposes, we shall consider two popular conditions: (i) the function should be small and (ii) the first derivatives of the function should be small (i.e. the function should vary slowly). For the first case we have
and for the second case
which is known as a firstorder Sobolev norm. Since ρ is a periodic function, we can write
where Δ is the Laplace operator [Δ = (∂^{2}/∂x_{i}^{2})] and (.,.) denotes the scalar product in Hilbert space.
Now the problem is reduced to finding the minimum of the functional
where L = I (identity operator) for L_{2}type regularizers (first case) and L = −Δ for Sobolevtype regularizers (second case).
Using Plancherel's theorem, the convolution theorem and the fact that the Fourier transformation of the Laplacian is proportional to the negative squared length of the reciprocalspace vector, we can rewrite the problem as
where F_{hkl} is the before sharpening (e.g. 2mF_{o} − DF_{c}type maps), F_{0hkl} is that after sharpening and s = 2sinθ/λ is the length of the reciprocalspace vector, with t(s) = 1, α = γ for regularization function f_{1} and t(s) = s^{2}, α = (2π)^{2}γ for f_{2}. This minimization problem has a very simple solution,
When k(x) is Gaussian then the equation has an especially simple form, since K(s) = [k(x)] = exp(−s^{T}B_{deblur}s/4), where B_{deblur} is an anisotropic deblurring B value.
Unfortunately, in reality neither B values nor α are known. Whilst there are several techniques to find an `optimal' value for α (Vogel, 2002) when the blurring function is known, in our implementation such an approach did not give consistent results. Therefore, we used the following ad hoc procedure for selection of the regularization parameter. Denoting K_{α}(s) = K(s)/[K^{2}(s) + αt(s)] and A_{α}(s) = K_{α}(s)K(s), we see that A_{α} is similar to the hat function used in regression analysis (Stuart et al., 2009). We can define the of errors (the number of observations minus the effective number of parameters) as^{1}
Note that when α = 0 then n_{df} = 0 and when α → ∞ then n_{df} is equal to the number of observations. We select α so that n_{df} is equal to 10–20% of the number of observations. Since we do not know the exact values of B and α, we also perform ad hoc integration using an empirically derived distribution of these parameters. The necessary integral may then be written
where B_{aniso} reflects the anisotropy of the data and is calculated during scaling of the calculated structure factors relative to the observed structure factors, under the conditions that it obeys crystal symmetry, and tr(B_{aniso}) = 0.
The joint probability distribution of B and α can be written
The mean value of the isotropic part is taken to be equal to the median value (B) of the coordinates (although it may be better to use Wilson's B value estimated using intensity curves derived by Popov & Bourenkov (2003). We approximate the distribution of the isotropic part of the B values using a Gaussian distribution centred at B_{sharp} with standard deviation equal to B_{sharp}/10. For each B value, we select α so that n_{df} is 10–20% of the number of observations and the standard deviation of the distribution of α is taken to be α_{B}/10.
Note that (17) and (18) suggest a class of regularizers. They can be selected to use particular knowledge about the electron density in real and For example, if it is desired to suppress the effect of ice rings then one can select t(s) so that the corresponding reflections are weighted down.
3.1. Implementations and an example
We have implemented anisotropic sharpening with L_{2} and Tikhonov–Sobolev regularizers with and without integration over the ad hoc joint probability distribution of B and α using probability distribution (24). We have also implemented the regularization function t(s) = 1 + s^{2}. These are available from REFMAC5 v.5.7. In our tests, all regularization functions gave similar results. This is not surprising, as the major problem is that the blurring function is not positionindependent. Before finding accurate regularizers, the problem of modelling positiondependent blurring functions should be dealt with. All results presented here were achieved using the L_{2}type regularizer.
Map sharpening was tested for many cases using data sets from the PDB (Berman et al., 2002) with resolution below 3 Å. The best results were obtained for PDB entry 2r6c (Bailey et al., 2007). For any lowresolution data taken from the PDB, before map calculation we generally try jellybody, local NCS (if present) and external reference structure (if applicable) and take the best refined results for further analysis. For 2r6c , the original R/R_{free} statistics reported in the PDB were 0.321/0.344. After these values became 0.240/0.300. Fig. 7 shows an illustration of the maps after with and without unregularized and regularized map sharpening. It is apparent that in this case using regularized mapsharpening coefficients shows more features (possibly side chains) and connectivity. Whilst this example shows regularization using the L_{2}type regularizer, it should be noted that the Sobolevtype regularizer gave similar results.
4. Conclusions and future directions
In this paper, we have presented two tools to aid in lowresolution
namely external structural restraints and regularized map sharpening.The use of external restraints based on homologous reference structures and/or structural fragments gives promising results. In particular, we have demonstrated how improved models can be achieved by the externally restrained rerefinement of deposited models.
Since the use of external restraints will alter global geometry validation statistics, such results should be interpreted accordingly and the integrity of local structure should always be considered. Indeed, local regions should still be manually inspected in order to ensure local suitability of the use of external restraints, despite any apparent improvements in overall statistics. If there are any serious artefacts that arise owing to bias towards the reference structure, it may be appropriate to exclude particular residues from external restraint generation.
In some cases, better results can be achieved by utilizing information from multiple reference structures, the difficulty often being that this requires the existence and availability of multiple structures suitable as references. Our implementation allows the generation of external restraints based on multiple reference structures; currently, the restraints most consistent with the target model are selected for use during refinement.
For practical application, we anticipate external restraints to also be of particular use during intermediate stages of model building/refinement, for stabilizing local structure and in helping to achieve sensible model geometry. Of course, the degree of any improvement owing to external restraints will be limited by the quality of the reference structure. For example, the MolProbity statistics for the reference structure 2jhp (used in one of our examples) are not particularly low given its resolution (see Table 1). Nevertheless, structural information contained in this reference model was able to improve the lower resolution target. However, the use of a more reliable reference model may have resulted in further improvements to the target structure 2jha . This scenario highlights the immediate need for ways to automatically validate the suitability of reference structures, most importantly at the local level, so that destructive restraints are not generated or are appropriately weighted down (whilst downweighting is already effectively performed by using robust estimators in our implementation, other complementary approaches would be desirable). In application, it may be sensible to attempt rerefinement of any reference structures before restraint generation in an attempt to improve the quality of the prior information. For example, this might be achieved automatically by using the PDB_REDO protocol (Joosten et al., 2009). In some cases, manual model rebuilding and of reference structures may be necessary/appropriate and thus should ideally always be considered. Such approaches may reduce any error propagation from reference to target models.
There is much room for improvement and future exploration in the generation and application of external structural information.
We have also implemented DNA/RNA basepair restraints based on interatomic distances, torsion angles and chirality; testing is currently in progress. Parameters for these restraints have been taken from Neidle (2008). For accurate of DNA/RNA it is necessary to use sugarpuckering as well as basestacking restraints. Whilst it is relatively simple to implement sugarpuckering restraints, e.g. using the elegant method presented by Cremer & Pople (1975), determining appropriate distributional parameters will take some time. Designing restraints for base stacking is a much more challenging problem, for which we do not currently have any satisfactory approaches.
The implemented method of regularized map sharpening uses the assumption that the blurring function is positionindependent. However, this assumption may not always be valid: it is expected that the oscillation of molecules within a crystal will be more complex and crystal disorder will be more anisotropic. One natural extension to map sharpening would be to use TLS parameters (Winn et al., 2001) as a blurring function. However, we are not aware of a simple solution to this problem. Future work will include deblurring using TLS parameters. Another problem with the current approach is that we assume that noise and signal are uncorrelated and that the noise is white noise. This may not reflect reality, especially when atomic model errors are dominating contributors to the noise. For density modification, the problem may become even more complicated.
Footnotes
^{1}This is the simplest way of defining effective of errors. Another formula is n_{df} = tr(I − A_{α})^{2}. Equation (20) has a simple form in when blurring is positionindependent then the sharpening matrices are diagonal.
Acknowledgements
Part of this work was carried out whilst the authors were at the Structural Biology Laboratory, Department of Chemistry, University of York, during which time RAN was funded by a BBSRC PhD Studentship and FL and GNM were funded by the Wellcome Trust. Part of this work was supported by the Medical Research Council (grant No. MC_US_A025_0104).
References
Adams, P. D. et al. (2010). Acta Cryst. D66, 213–221. Web of Science CrossRef CAS IUCr Journals
Afonine, P. V., GrosseKunstleve, R. W. & Adams, P. D. (2005). Acta Cryst. D61, 850–855. Web of Science CrossRef CAS IUCr Journals
Bailey, S., Eliason, W. & Steitz, T. (2007). Science, 318, 459–463. CrossRef PubMed CAS
Bentley, J. L. (1975). A Survey of Techniques for Fixed Radius Near Neighbor Searching. Technical Report SLAC186 and STANCS75513. Stanford Linear Accelerator Center.
Berman, H. M. et al. (2002). Acta Cryst. D58, 899–907. Web of Science CrossRef CAS IUCr Journals
Blanc, E., Roversi, P., Vonrhein, C., Flensburg, C., Lea, S. M. & Bricogne, G. (2004). Acta Cryst. D60, 2210–2221. Web of Science CrossRef CAS IUCr Journals
Bricogne, G. (1997). Methods Enzymol. 276, 361–423. CrossRef CAS Web of Science
Brünger, A. T. (1997). Methods Enzymol. 277, 366–396. CrossRef PubMed CAS Web of Science
Brünger, A. T., Adams, P. D., Clore, G. M., DeLano, W. L., Gros, P., GrosseKunstleve, R. W., Jiang, J.S., Kuszewski, J., Nilges, M., Pannu, N. S., Read, R. J., Rice, L. M., Simonson, T. & Warren, G. L. (1998). Acta Cryst. D54, 905–921. Web of Science CrossRef IUCr Journals
Brunger, A. T., DeLaBarre, B., Davies, J. M. & Weis, W. I. (2009). Acta Cryst. D65, 128–133. Web of Science CrossRef CAS IUCr Journals
Chen, V. B., Arendall, W. B., Headd, J. J., Keedy, D. A., Immormino, R. M., Kapral, G. J., Murray, L. W., Richardson, J. S. & Richardson, D. C. (2010). Acta Cryst. D66, 12–21. Web of Science CrossRef CAS IUCr Journals
Cremer, D. & Pople, J. (1975). J. Am. Chem. Soc. 97, 1354–1358. CrossRef CAS Web of Science
Davis, I., LeaverFay, A., Chen, V., Block, J., Kapral, G., Wang, X., Murray, L., Arendall, W., Snoeyink, J., Richardson, J. & Richardson, D. (2007). Nucleic Acids Res. 35, W375–W383. CrossRef PubMed
DeLaBarre, B. & Brunger, A. T. (2006). Acta Cryst. D62, 923–932. Web of Science CrossRef CAS IUCr Journals
Geman, S. & McClure, D. (1987). Bull. Int. Stat. Inst. 52, 5–21.
Gonzalez, R. C. & Woods, R. E. (2002). Digital Image Processing. New Jersey: Prentice Hall.
Joosten, R. P. et al. (2009). J. Appl. Cryst. 42, 376–384. Web of Science CrossRef CAS IUCr Journals
Kavanaugh, J., Rogers, P. & Arnone, A. (2005). Biochemistry, 44, 6101–6121. CrossRef PubMed CAS
Levinthal, C. (1966). Sci. Am. 214, 42–52. CrossRef CAS PubMed
Long, F., Vagin, A. A., Young, P. & Murshudov, G. N. (2008). Acta Cryst. D64, 125–132. Web of Science CrossRef CAS IUCr Journals
McNicholas, S., Potterton, E., Wilson, K. S. & Noble, M. E. M. (2011). Acta Cryst. D67, 386–394. Web of Science CrossRef CAS IUCr Journals
Mooij, W., Cohen, S., Joosten, K., Murshudov, G. & Perrakis, A. (2009). Structure, 17, 183–189. Web of Science CrossRef PubMed CAS
Murshudov, G. N. (2011). Appl. Comput. Math. 10, 250–261.
Murshudov, G. N., Skubák, P., Lebedev, A. A., Pannu, N. S., Steiner, R. A., Nicholls, R. A., Winn, M. D., Long, F. & Vagin, A. A. (2011). Acta Cryst. D67, 355–367. Web of Science CrossRef CAS IUCr Journals
Murshudov, G. N., Vagin, A. A. & Dodson, E. J. (1997). Acta Cryst. D53, 240–255. CrossRef CAS Web of Science IUCr Journals
Neidle, S. (2008). Principles of Nucleic Acid Structures. London: Academic Press.
Nicholls, R. (2011). PhD thesis. Department of Chemistry, University of York, UK.
Nocedal, J. & Wright, S. (1999). Numerical Optimization. New York: Springer.
Popov, A. N. & Bourenkov, G. P. (2003). Acta Cryst. D59, 1145–1153. Web of Science CrossRef CAS IUCr Journals
R Development Core Team (2011). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing. http://www.Rproject.org .
Savino, C., Miele, A., Draghi, F., Johnson, K., Sciara, G., Brunori, M. & Vallone, B. (2009). Biopolymers, 91, 1097–1107. CrossRef PubMed CAS
Schröder, G., Brünger, A. & Levitt, M. (2007). Structure, 15, 1630–1641. PubMed
Schröder, G., Levitt, M. & Brünger, A. (2010). Nature (London), 464, 1218–1222. PubMed
Sheldrick, G. M. (2008). Acta Cryst. A64, 112–122. Web of Science CrossRef CAS IUCr Journals
Smart, O. S., Womack, T. O., Flensburg, C., Keller, P., Sharff, A., Paciorek, W., Vonrhein, C. & Bricogne, G. (2012). Acta Cryst. D68, 368–380. CrossRef IUCr Journals
Stuart, A., Ord, K. & Arnold, S. (2009). Kendall's Advanced Theory of Statistics, Vol. 2A, Classical Inference. Hoboken: John Wiley & Sons.
Sutton, G., Grimes, J., Stuart, D. & Roy, P. (2007). Nature Struct. Mol. Biol. 14, 449–451. CrossRef CAS
Tickle, I. J. (2007). Acta Cryst. D63, 1274–1281. Web of Science CrossRef CAS IUCr Journals
Tikhonov, A. N. & Arsenin, V. Y. (1977). Solution of Illposed Problems. Washington: Winston & Sons.
Vagin, A. A., Steiner, R. A., Lebedev, A. A., Potterton, L., McNicholas, S., Long, F. & Murshudov, G. N. (2004). Acta Cryst. D60, 2184–2195. Web of Science CrossRef CAS IUCr Journals
Vogel, C. R. (2002). Computational Methods for Inverse Problems. Philadephia: SIAM.
Winn, M. D. et al. (2011). Acta Cryst. D67, 235–242. Web of Science CrossRef CAS IUCr Journals
Winn, M. D., Isupov, M. N. & Murshudov, G. N. (2001). Acta Cryst. D57, 122–133. Web of Science CrossRef CAS IUCr Journals
This is an openaccess article distributed under the terms of the Creative Commons Attribution (CCBY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.