research papers
Coping with strong translational Phaser: human Rab27a
and extreme anisotropy in withaFaculty of Natural Sciences, Department of Life Sciences, Imperial College London, Exhibition Road, South Kensington, London SW7 2AZ, England, bDepartment of Chemistry, Molecular Sciences Research Hub, Imperial College London, Wood Lane, London W12 0BZ, England, and cCambridge Institute for Medical Research and Department of Haematology, University of Cambridge, Wellcome Trust/MRC Building, Hills Road, Cambridge CB2 0XY, England
*Correspondence e-mail: e.cota@imperial.ac.uk, rjr27@cam.ac.uk
Data pathologies caused by effects such as diffraction anisotropy and translational Mut, crystallized within 24 h and diffracted to 2.82 Å resolution, with a possessing room for a large number of protein copies. Initial efforts to solve the structure using by Phaser were not successful. Analysis of the data set revealed that the crystals suffered from both extreme anisotropy and strong tNCS. As a result, large numbers of reflections had estimated standard deviations that were much larger than their measured intensities and their expected intensities, revealing problems with the use of such data at the time in Phaser. By eliminating extremely weak reflections with the largest combined effects of anisotropy and tNCS, these problems could be avoided, allowing a molecular-replacement solution to be found. The lessons that were learned in solving this structure have guided improvements in the numerical analysis used in Phaser, particularly in identifying diffraction measurements that convey very little information content. The calculation of information content could also be applied as an alternative to ellipsoidal truncation. The post-mortem analysis also revealed an oversight in accounting for measurement errors in the fast rotation function. While the crystal of mutant Rab27a is not amenable to drug screening, the structure can guide new modifications to obtain more suitable crystal forms.
(tNCS) can dramatically complicate the solution of the crystal structures of macromolecules. Such problems were encountered in determining the structure of a mutant form of Rab27a, a member of the Rab GTPases. Mutant Rab27a constructs that crystallize in the free form were designed for use in the discovery of drugs to reduce primary tumour invasiveness and metastasis. One construct, hRab27aKeywords: molecular replacement; Phaser; anisotropy; translational noncrystallography symmetry; Rab27a; information content.
PDB reference: human Rab27a, 6huf
1. Introduction
Accounting rigorously for the effects of errors in a statistical model can dramatically enhance the sensitivity of likelihood-based methods. For instance, in molecular-replacement (MR) calculations, Phaser (McCoy et al., 2007) is able to account for the effects of errors in both the search model and in the measured diffraction data; this is difficult to achieve with methods based on the properties of the or on the computation of correlation coefficients. In addition, information obtained from already placed search components significantly improves the signal in rotation and translation searches for subsequent components, as measured by the log-likelihood gain (LLG) and Z-scores (McCoy, 2007; Storoni et al., 2004; McCoy et al., 2005).
This sensitivity is a double-edged sword, as likelihood-based methods are also highly sensitive to defects in their statistical models. For this reason, in crystallographic applications it is essential to account for the statistical effects of anisotropy (McCoy et al., 2007) and translational (tNCS; Sliwiak et al., 2014). The likelihood targets in versions of Phaser since v.2.5.4 account for the statistical effects of tNCS arising from translations combined with small changes in conformation and orientation differences up to 10°. These yield tNCS correction parameters describing changes in the expected intensities (and their probability distribution). Automated algorithms in Phaser can deal with simple cases of tNCS, for instance a single tNCS vector between two groups of molecules, but manual intervention by the user can be required for more complex situations, which includes a complete understanding of the cell content and identifying the tNCS vectors between the molecules (Sliwiak et al., 2014).
One consequence of the intensity modulations introduced by significant anisotropy and/or tNCS is that there are bound to be systematically weak intensities with relatively large measurement errors, regardless of any overall resolution limit applied to the data. In these circumstances, it is particularly important to account rigorously for the effects of intensity-measurement error, for instance with the log-likelihood gain on intensities (LLGI) target (Read & McCoy, 2016). Problems encountered in solving the structure of Rab27a have highlighted the importance of these issues.
Rab27a is a small GTPase belonging to the large family of Ras-related in brain (Rab) proteins. Rab27a is part of the secretory pathway involved in the transport of melanosomes (Strom et al., 2002) and the secretion of vesicles containing insulin (Yamaoka et al., 2015), histamine (Goishi et al., 2004), chemokines, matrix metalloproteases (MMPs) and exosomes (Fukuda, 2013; Brozzi et al., 2012; Ostrowski et al., 2010). In humans, Rab27a is overexpressed in multiple types of cancer, including breast (Wang et al., 2008), lung (Li et al., 2014), pancreatic (Wang et al., 2015) and liver (Dong et al., 2012) cancers.
Evidence supporting the role of human Rab27a (hRab27a) in multiple cancer types suggests that the inhibition of this GTPase could be a target for cancer therapy. Therefore, structural characterization of Rab27a is required for the development of specific inhibitors. Crystallographic structures of mouse Rab27a and Rab27b (mRab27a and mRab27b) in complex with the human Slp2a and Slac2a (hSlp2a and hSlac2a) effectors have been reported (Kukimoto-Niino et al., 2008; Chavas et al., 2008). Potential ligandable sites are located at or near the mRab27–hSlp2a and mRab27–hSlac2a interfaces, and therefore these complexes cannot be used for the characterization of Rab27a–ligand complexes. While the crystallization of Rab27a on its own would be the ideal solution to this problem, this has been unsuccessful for the human and mouse homologues (Chavas et al., 2010). We therefore generated hRab27a mutants that were capable of crystallizing in the absence of effectors and were suitable for ligand-binding studies. Point mutations in hRab27a were made based on the crystal packing of mouse Rab3, the highest identity hRab27a homologue with known structure (Dumas et al., 1999). This led to a construct, referred to as hRab27aMut, that is able to form crystals that diffract to a maximum resolution of 2.82 Å and with the potential ligand-binding sites exposed. A complete description of the design of these mutants will be reported elsewhere.
Initial attempts to solve the structure by MR using Phaser (McCoy et al., 2007) were unsuccessful. Inspection of the X-ray data showed that these crystals were highly anisotropic, the native indicated strong translational noncrystallographic symmetry (tNCS) and a high copy number was predicted per asymmetric unit.
Here, we describe the solution of this difficult MR problem, as well as the improvements that the experience has inspired in Phaser. Moreover, the has given us directions for further improvements in the design of Rab27a constructs that crystallize in the free form suitable for ligand discovery, which will be reported in detail elsewhere.
2. Materials and methods
2.1. Protein production
The cDNA template for hRab27a (UniProt code P51159) was kindly provided by Dr Miguel Seabra (Imperial College London). A gene corresponding to residues 1–192 was amplified from this cDNA and cloned into the pET-15b plasmid, generating the pET-15b-rab27a construct. The construct contains an N-terminal His tag followed by a Tobacco etch virus (TEV) protease cleavage site. PCR amplification was performed using Q5 High-Fidelity DNA Polymerase (New England Biolabs; NEB); the 5′-CGGCTCATATGTCTGATGGAGATTATGATTAC-3′ and 5′-CGGCTGGATCCTCAGGACTTGTCCACACTCC-3′ were used as the forward and reverse primers, respectively. A Q5 Site-Directed Mutagenesis Kit (NEB) was used to introduce several mutations (Q105E, Q118K, M119T, Q140E, K144A, E145A, E146A, I149R, A150Q and K154H; the Arg50–His69 loop was replaced with the sequence TIYRNDKRIK) in the pET-15b-rab27a construct to generate the pET-15b-hrab27amut construct. A Q78L mutation was introduced to decrease the GTPase activity of the protein, and C123S and C188S mutations were used to avoid aggregation during protein preparation. A glycine would remain as the initial residue after tag removal using TEV protease.
For the production of hRab27aMut, the pET-15b-hrab27amut construct was transformed into Escherichia coli BL21 (DE3) cells (NEB). The bacteria were grown in lysogenic broth (LB) at 37°C to an OD at 600 nm of 0.6–0.8, and protein expression was then induced with 0.5 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) at 37°C for 3 h. The cells were harvested by centrifugation at 4000 rev min−1 for 10 min at room temperature. The cell pellets were resuspended in 50 mM Tris–HCl pH 8.0, 500 mM NaCl, 5 mM MgCl2 (buffer A) supplemented with 10 mM imidazole. The cells were lysed with a cell disruptor (Constant Systems) at 172 MPa and centrifuged at 15 000 rev min−1 for 45 min at 4°C. The supernatant was loaded onto an Ni–NTA affinity column (Qiagen) equilibrated in buffer A supplemented with 10 mM imidazole. The resin was washed with 20 volumes of buffer A with 10 mM imidazole, and the protein was then eluted in buffer A with 300 mM imidazole. The protein was dialyzed against buffer B (50 mM Tris–HCl pH 8.0, 100 mM NaCl, 5 mM MgCl2) and the His tag was removed by overnight incubation with TEV protease (His-tagged) at a molar ratio of 1:20 in buffer B supplemented with 1 mM DTT at 4°C. DTT was removed by dialysis against buffer B and the protein was reloaded onto an Ni–NTA column to remove TEV protease and uncleaved protein. The purity was assessed by SDS–PAGE. The protein concentration was determined by UV–Vis absorption at 280 nm using a Nanodrop spectrophotometer (ThermoFisher).
The locked-active (GTP-bound) form of hRab27aMut was obtained by loading the protein with the nonhydrolysable GTP analogue GppNHp (Jena Bioscience). GppNHp was loaded by overnight incubation of 10 mg hRab27aMut with 25 units of Antarctic Phosphatase (NEB) in buffer B with 1 mM zinc chloride, 0.2 M ammonium sulfate and a fourfold molar excess of GppNHp in a final reaction volume of 2 ml at 4°C. The GTPase was further purified by with a Superdex 75 HiLoad (10/30) column (GE Healthcare) equilibrated in 20 mM Tris–HCl pH 8.0, 150 mM NaCl, 5 mM MgCl2. The eluted protein was concentrated to 25 mg ml−1 and flash-frozen in liquid nitrogen for storage.
2.2. Crystallization and X-ray data collection
Sitting-drop vapour-diffusion crystallization experiments with hRab27aMut(GppNHp) were set up using a Mosquito robot (TTP Labtech) at 20°C. A search for crystallization conditions was performed using ∼1000 commercial conditions. Drops consisting of 400 nl were formed by mixing equal volumes of protein solution and precipitant solution. The best crystals were obtained in 20%(v/v) ethylene glycol, 10%(w/v) PEG 8000, 30 mM MgCl2, 30 mM CaCl2, 100 mM HEPES pH 7.5 after 3–4 days at 20°C. Crystals were cryoprotected in the crystallization-condition solution supplemented with 30%(v/v) ethylene glycol and were flash-cooled in a nylon loop in liquid nitrogen. A complete X-ray data set to 2.82 Å resolution was collected at 100 K on beamline I02 at Diamond Light Source (DLS), Oxford, England. The data were processed and scaled with DIALS (Waterman et al., 2016; Winter et al., 2018), POINTLESS (Evans, 2011) and AIMLESS (Evans & Murshudov, 2013) using the xia2 pipeline (Winter, 2010). Statistics for the data collection are presented in Table 1. An initial model generated by with Phaser was refined through an iterative cycle using Coot (Emsley et al., 2010) and REFMAC5 (Winn et al., 2003). The final model structures were validated using the MolProbity server (Chen et al., 2010) at https://molprobity.biochem.duke.edu. All structure images were prepared using PyMOL (Schrödinger).
‡. §. |
A self-rotation function was calculated with MOLREP (Vagin & Teplyakov, 2010). Native Patterson maps were calculated with the FFT program (Ten Eyck, 1973) from the CCP4 package (Winn et al., 2011). Anisotropic atomic displacement parameters, including the anisotropic delta-B, were calculated using the ANO (anisotropy) mode and tNCS expected intensity factors using the TNCS mode in Phaser. SFTOOLS from the CCP4 package (B. Hazes, unpublished results) was used to combine the anisotropy and tNCS factors, to select a subset of data for the initial structure solution and to compute the equivalent resolution corresponding to a full data set with a specified number of reflections. The Matthews coefficient (Matthews, 1968) and solvent-content calculations for different possible compositions of the were carried out with MATTHEWS_COEF from the CCP4 package (Winn et al., 2011).
3. Results and discussion
3.1. composition and translational noncrystallographic symmetry
The Mut(GppNHp) crystal was estimated to contain a large number of GTPase molecules (between 16 and 24; see Table 2; Kantardjieff & Rupp, 2003; Matthews, 1968; McCoy, 2007). With high NCS, the contribution of each component is small, making structure solution by MR much more challenging.
of the hRab27a
|
The self-rotation function reveals the angular relationship between two or more identical molecules in the ω, φ and κ spherical polar angles. Self-rotation function peaks often correspond to rotational NCS in the crystal (Drenth, 2007). There is a κ = 90° (ω = [90°], φ = [54°]) peak in the self-rotation function (Fig. 1a), corresponding to a fourfold rotation axis. There are also 13 κ = 180° peaks corresponding to twofold rotation axes. One interpretation of this is that there are two assemblies with dihedral D4 point-group symmetry in the crystal, with the two fourfold axes parallel.
This function measures the correlation of the native with a rotated copy, often calculated usingTranslational et al., 2013). This changes the usual Wilson distribution of structure-factor intensities (Read et al., 2013; Wilson, 1949). The calculation of a native for the hRab27aMut(GppNHp) data reveals a peak at fractional coordinates (0.000, 0.022, 0.500) of 45% of the height of the origin peak (Fig. 1b), showing strong tNCS that broadens the intensity distribution; because of the half-unit-cell component of the translation along the c axis, reflections with l odd will tend to be very weak, although this will be modulated by the size of the k index (because of the small but not insignificant translation along the b axis).
(tNCS) occurs when two or more independent copies of a molecule have similar orientations in the tNCS-related molecules would contribute with the same or similar amplitudes to a However, their relative phases are determined by the projection of the translation vector on the diffraction vector, resulting in systematic interference that generates stronger and weaker reflections (Read3.2. Extreme diffraction anisotropy
The hRab27aMut diffraction pattern was extremely anisotropic (Fig. 2), with the diffraction intensity falling off at different rates in different reciprocal-lattice directions. This is potentially owing to the pattern of lattice contacts in the crystal, which can give variations in the relative ordering of molecules along different directions. If not accounted for, the presence of significant anisotropy in the data will affect the likelihood functions used by Phaser, so it is important to refine and apply anisotropic correction factors. The degree of anisotropy of an X-ray data set can be described using the anisotropic delta-B, which is the difference between the two most extreme principal components of the anisotropic atomic displacement parameter along different directions in Delta-B values of above 50 Å2 are considered to indicate extreme anisotropy. The diffraction anisotropy of the hRab27aMut crystals was estimated with the ANO mode of Phaser to be 207.3 Å2. This indicates a case of severe diffraction anisotropy (Fig. 2a), with an effective resolution of 2.82 Å in the strongest direction and 5.0 Å in the weakest direction (Fig. 2b).
3.3. Solving the molecular-replacement problem
After failed attempts to solve the structure with Phaser using the structure of mRab27a as a model, we used Sculptor (Bunkóczi & Read, 2011a) and Ensembler (Bunkóczi & Read, 2011b) to generate an optimized ensemble model. This ensemble was generated on the basis of the closest homologue structures reported for hRab27a: mRab27a(GppNHp) (PDB entry 3bc1; 87% identical in amino-acid sequence; Chavas et al., 2008), mRab27b(GDP) (PDB entry 2iey; 68% identical; Chavas et al., 2007), mRab27b(GppNHp) (PDB entry 2zet; Kukimoto-Niino et al., 2008) and human Rab8a(GppNHp) (PDB entry 4lhw; 49% identical; Guo et al., 2013). Regions with different conformations among the input models were removed using the `trim' option of Ensembler (Fig. 3a). A MR calculation with Phaser using this ensemble failed in the first attempt, where a solution was found for only one pair of tNCS-related copies.
It appears that the combination of strong tNCS and extremely high anisotropy led to a very wide distribution of expected intensities, with many reflections expected to have extremely weak intensities based on these systematic effects. In addition, the high number of molecules in the Phaser (Read & McCoy, 2016) should compensate for the effects of anisotropy and tNCS by downweighting the systematically weak reflections with standard deviations that are large compared with their expected intensities. However, there could potentially be significant errors in the estimates of the standard deviations, as well as in the anisotropy and/or tNCS correction factors applied to the expected intensities. In addition, the presence of reflections with standard deviations much larger than their expected intensities could lead to numerical instabilities in the evaluation of the intensity-based likelihood target. To avoid these potential problems, the systematically weakest reflections with the largest relative errors were omitted from the molecular-replacement calculations. The anisotropic scale factors and tNCS scale factors were calculated using the ANO (anisotropy) and TNCS modes, respectively, in Phaser. Using SFTOOLS, these correction factors were multiplied together and then used to discard the systematically weakest intensities. In the initial calculation with the pruned data, any reflection for which the combined correction factor was greater than 10 was discarded; as a result, around 40% of the data were discarded (Fig. 4). Although both tNCS and anisotropy are present, for this data set by far the largest corrections arise primarily from the effects of anisotropy. The correction factors for anisotropy vary by a factor of nearly 330 000, while those for tNCS vary by a factor of less than 700, combining to give an overall variation by a factor of about 900 000 (Fig. 4). Note that the largest effects of tNCS are seen at low resolution, where small rotations and conformational differences have less effect on the correlations between the structure-factor contributions of tNCS-related molecules, while the largest effects of anisotropy are seen at high resolution; as a result, the range of the combined effects of tNCS and anisotropy is smaller than one would expect if the two effects varied independently.
is likely to complicate the rotation and translation search functions. In principle, the new intensity-based likelihood target inUsing the trimmed data, a clear and correct molecular-replacement solution could be found with a TFZ score of 12.8 for the final copy, placing 16 copies of the trimmed ensemble model in a physically plausible crystal-packing arrangement (Fig. 3b); solutions with a TFZ of greater than 8 are almost always correct (Oeffner et al., 2013). Testing different thresholds for the scaling-factor cutoff suggested that a 50× scaling-factor cutoff still gave an equivalent MR solution, enabling us to cut only 20% of the reflections. Density for the nucleotide, which was not included in the model, was observed in the NCS-averaged 2Fo − Fc and Fo − Fc electron-density maps (Fig. 3c), strongly suggesting that the molecular-replacement solution was correct.
The solution is also consistent with the self-rotation function. The D4 assemblies that superpose with very low r.m.s.d. values (<0.1 Å) using molecules A and I of T1 and T3, indicating that they have the same conformation/structure (Figs. 5a and 5b). The fourfold axis of the octamer correlates with the peak in the self-rotation function at κ = 90° (κ = 90°, φ = ±180°, ϕ = 54°) and κ = 180° (κ = 180°, φ = ±180°, ϕ = 54° for the twofold axis within the same tetramer) (Fig. 1). The twofold axes relating molecules in one tetramer to molecules within other tetramers explain the peaks observed in the self-rotation function at κ = 180°. The peaks labelled 1–13 correlate to twofold axes between molecules in T1–T3, T1–T4, T2–T3 and T3–T4 (Fig. 1). A full description of the relationships is given in Table 3. In agreement with the prominent off-origin peak in the native translational symmetry between the two octamers is observed in the structure (Fig. 5e).
consists of two octamers, giving two
|
The structure was completed and refined using Coot for manual rebuilding and REFMAC5 for during which restraints were applied. Most residues in all 16 molecules were modelled, apart from flexible residues at the N-terminus of the construct. Residues with poor side-chain density (930 out of a total of 2736 in the model) were truncated at the Cβ atom. The final used a pruned data set from which reflections conveying less than 0.05 bits of information (24% of the data set) were removed, as discussed below. The agreement with the measured data (Rfree = 0.342 and Rwork = 0.312) is consistent with what one might expect from a data set containing 69 568 reflections; this corresponds to the number of reflections that would be contained in a complete isotropic data set at a resolution of 3.09 Å. The coordinates and structure factors have been deposited in the wwPDB (Berman et al., 2007) as PDB entry 6huf.
In the Rab27a structure, the SF4 pocket, formed by the α3–β5 loop (a highly variable region among Ras superfamily members) and the C-terminal region of the α5 helix, is of particular interest, as it is fundamental to the interaction of Rab27a with the WF motif of Slp2a. A model was built for the SF4 pocket in all 16 molecules of the solution structure. Interestingly, the pocket is free from contacts with neighbouring symmetry-related molecules (Fig. 6), making it suitable for protein–ligand interaction studies if the problems with anisotropy in the data could be resolved.
3.4. Excluding systematically weak data based on information content
Subsequent to, and inspired by, this structure solution, an automated method to exclude the systematically weakest reflections from the MR likelihood calculations has been implemented in Phaser. The method applied in the initial structure solution was chosen to eliminate the reflections that would suffer most from the combined effects of anisotropy and tNCS, but it did not account for the precision of the individual measurements.
The new method considers the precision of the measurement relative to the intensity expected for the particular reflection when the effects of anisotropy and tNCS are taken into account. One way to evaluate the precision of a measurement is to consider how much information that measurement conveys; in other words, how much more is known after making the measurement than before. This information gain can be evaluated by considering the loss of relative ) distribution of true intensities] to the posterior probability distribution. In information theory, this quantity is known as the Kullback–Leibler divergence or KL-divergence (Kullback & Leibler, 1951), which is defined in (1) and is represented subsequently as simply DKL:
in going from the prior probability distribution [the null hypothesis, in this case the Wilson (1949If the natural logarithm is used in this expression, the information content is expressed in units of nats, whereas the equivalent expression using the base 2 logarithm gives information in terms of bits, which can therefore be obtained from that in nats by dividing by ln(2). The KL-divergence is always non-negative, but because the integral is weighted by only one of the two probability distributions it is not symmetric and is therefore not strictly a distance metric.
This information-based measure is a natural choice in the context of likelihood-based optimization methods. If in the KL-divergence in (1) the prior probability is replaced by a prior probability conditional on a model, then it can be shown that maximizing a likelihood function (i.e. the probability of the data given the model) is equivalent to minimizing this KL-divergence (Bishop, 2006). In other words, maximizing the likelihood minimizes the divergence between the probability of the true value of the data given the model and the probability of the true value of the data given the measurements of the data.
For diffraction data measured in terms of intensities and their estimated standard deviations, the expressions are simpler if cast in terms of normalized intensity values, for which the expected true intensity is 1, i.e. E2. For clarity, we will represent the normalized intensity as Z (= E2). The prior probability is simply the Wilson (1949) distribution of normalized intensities, given in (2a) for the acentric case and (2b) for the centric case:
In computing the KL-divergence for diffraction intensities, the posterior probability of the true intensity given the measured intensity, which plays a key role in the procedures of French & Wilson (1978), can be defined in terms of other probabilities using Bayes' theorem (3), yielding (4):
In this equation, the probability distribution for the observed intensity given the true intensity is taken as the Gaussian distribution in (5),
The probability distribution for the observed normalized intensity is given by (6a) for acentric reflections and by (6b) for centric reflections, which are reproduced from equations (9a) and (9b) of Read & McCoy (2016):
In (6b), erfc is the complement of the error function and D is a parabolic cylinder function (Whittaker & Watson, 1990).
The integral in (4) could be used to evaluate the information content of individual reflections, and a minimum information content could be defined for reflections that are accepted for further calculations. We chose instead to evaluate and use the expected value of the information content, based only on the estimated standard deviation and ignoring the particular value found for the measured intensity. The primary argument for this choice is that outlier observations are probably more likely to be encountered for the systematically weak intensities, partly because of inaccuracies in the determination of the correction factors; outliers that are substantially larger than expected will be evaluated, according to (4), as conveying more information and would thus be more likely to be kept in the data set. An additional advantage to using the expected information content is that this is a function of only the standard deviation of the normalized intensity, so a simple threshold can be set. In contrast, evaluating the integral in (4) is considerably more difficult, but in the future we will test whether there is a practical difference in outcome.
The expected information content is a probability-weighted average over all possible values of the measured intensity, given in (7):
The derivation of (7) implicitly assumes that the standard deviation of the intensity is independent of the measured intensity, which would not be valid for well measured intensities. However, the information thresholds are only applied in practice to observations in which the uncertainty of the measurement is at least several times larger than the expected intensity itself (see below); in these circumstances the uncertainty comes primarily from the counting statistics of the background rather than the peak.
To construct lookup tables for normalized intensity standard deviation thresholds, (7) was evaluated by numerical integration in Mathematica v.10 (Wolfram Research, Champaign, Illinois, USA) for a variety of expected information-content thresholds. Information-content filtering based on these thresholds was implemented in Phaser (McCoy et al., 2007), with the feature being available in v.2.7.17 (November, 2016) or newer. Note that the systematically weak reflections contribute to the of parameters describing the anisotropy and tNCS, and are only excluded for subsequent MR likelihood calculations; for this reason, it is better to provide the full, unpruned set of data to Phaser.
An examination of (7) gives further insight into the connection between the KL-divergence and likelihood. The form of this equation is highly reminiscent of the expected log-likelihood gain (eLLG) used to predict the outcome of molecular-replacement calculations, as defined in equation (3) of McCoy et al. (2017). This equation can be recast in terms of observed intensities rather than effective amplitudes, yielding (8),
For the case of a perfect model, where the calculated
is identical to the true this equation for the eLLG is equivalent to the expected KL-divergence. In other words, the expected KL-divergence corresponds to the estimated maximum contribution of an observation to the total likelihood that could be achieved with a perfect model.3.5. Accounting for measurement error in the likelihood-based fast rotation function
Inspection of the log files obtained in the initial structure solution before and after pruning the data with the largest anisotropy and tNCS correction factors suggested that the greatest improvements from omitting systematically weak data were in the results of the fast rotation function. This revealed an oversight in the implementation of the intensity-based LLGI target function in Phaser (Read & McCoy, 2016). In almost all cases, implementing this target simply involves replacing the structure-factor amplitude with an effective amplitude, Feff, and applying an additional factor Dobs to any σA values in the likelihood targets; both Feff and Dobs are derived from the intensity and its standard deviation (Read & McCoy, 2016).
Applying this to the likelihood-based fast rotation function, LERF1 (Storoni et al., 2004), requires a slightly different approach. LERF1 is based on a first-order series expansion of the log of the rotation likelihood function, given in (9) (adapted from equation 17 of Storoni et al., 2004),
where χΩ is the Fourier transform of the sphere inside of which Patterson-like functions of the observed intensities and contributions of the fixed and rotating components of the model are compared as a function of rotation. (Note that the post-multiplication of k by R−1 corresponds in to rotating the calculated Patterson in by pre-multiplying the coordinates by R.) The Patterson-like functions I1t and I1s are defined in (10a)–(10c), which are adapted from equations (18) and (19) of Storoni et al. (2004):
In (10b) and (10c) D is the Luzzati factor (Luzzati, 1952), which is proportional to σA. In the initial adaptation of LERF1 to the LLGI intensity-based likelihood target, any instances of D in the variance term ΣN′ in (10a) were multiplied by Dobs. However, the Luzzati factor D in (10b) was not modified, because rotation of the model associates different indices k with the observed reflections indexed by h. To compensate in (9) for this omission, the expression for I1t has to be multiplied by D2obs. This correction was introduced into Phaser at the same time as the filtering on information content.
3.6. Tests of modified Phaser
As described above, eliminating the systematically weakest reflections from the data set was sufficient to give a clear solution to the hRab27a structure, even before the fast rotation function was modified to properly account for intensity-measurement errors.
With the new algorithms, the hRab27a structure and others suffering from severe anisotropy and/or tNCS can now be solved more easily and without manual intervention. Table 4 illustrates the effect of applying different information-content thresholds on the course of the molecular-replacement calculation. With the corrected fast rotation function, it is no longer necessary to prune the systematically weak reflections in order to obtain a solution. Pruning up to about 19% of the weakest reflections in this data set (those conveying less than 0.01 bits of information each) has very little effect on the signal; if anything, the final LLG value increases very slightly. For this case at least, there is very little disadvantage to including even exceptionally weak data as long as the effects of measurement errors are accounted for properly. The main effect is a tendency for the total computing time to increase with the number of reflections included. (Note that there is a stochastic element to the total computing time, which is influenced by the number of potential partial solutions identified at any point in the calculation.) For other cases, where the estimates of measurement errors might be poorer or where the effects of anisotropy and/or tNCS might be modelled less accurately, omitting the weakest reflections might still improve the outcome of the calculation.
‡Reflections rejected as Wilson distribution outliers. |
However, our experience with the oversight in the implementation of the fast rotation function shows that when an algorithm fails to account properly for the effects of measurement error, there is a real advantage to pruning the weakest data. In the uncorrected fast rotation function, terms corresponding to weak observations with little information content were being included at a higher weight than they should have been given. The same general effect will apply in any other calculation in which weak data are not appropriately downweighted. For instance, the use of amplitudes and their standard deviations obtained through the French & Wilson (1978) algorithm in amplitude-based likelihood targets will overweight extremely weak data because the French and Wilson amplitude standard deviation has a finite value even in the limit of intensities with infinite measurement error (Read & McCoy, 2016).
The relationship between the expected LLG and the expected KL-divergence (equations 7 and 8) shows that even for a model approaching perfection, the omission of data with low information content will have very little effect on a properly calculated likelihood function, indicating that such observations should have very little leverage. For instance, measurements contributing 0.01 bits of information will contribute at most 0.01ln(2) to the likelihood score, so it would take over 140 such observations to change the likelihood score by a single unit. If such observations are omitted from algorithms in which the effects of errors are not properly accounted for, this will remove a potential source of systematic bias or noise.
The expected information content could therefore potentially be used as an alternative to ellipsoidal truncation to prune weak data (Strong et al., 2006). The initial approach, that of pruning the reflections with the highest combined anisotropy and tNCS correction factors, led to a successful structure solution but does not work nearly as well. For instance, if the 23 629 reflections with a combined intensity-correction factor of greater than 60 are omitted, the final LLG decreases from 3667.3 to 3560.2, whereas if the 23 868 reflections conveying less than 0.1 bits of information are omitted the final LLG only decreases to 3646.8. As a less extreme example, 17 457 reflections have a combined correction factor of greater than 160; if these are omitted the final LLG decreases to 3659.3, whereas setting the information-content threshold to 0.01 bits actually gives a slight increase in LLG while omitting a very similar number of reflections (Table 4).
Based on these data and similar tests on other systems (results not shown), the default threshold chosen for likelihood calculations in Phaser is 0.01 bits of information per reflection; note that all data should be used in the data-preparation calculations in Phaser that characterize anisotropy and tNCS effects. Optimal thresholds for computations in other software are likely to differ from this. In addition, the information calculations depend on the accuracy of the parameters describing anisotropy and tNCS, and do not yet account for other effects on intensities such as those from or order–disorder structures. The full data set should therefore always be maintained without permanently excluding data at any information threshold.
4. Conclusions
The hRab27aMut(GppNHp) data show how difficult cases of can be solved using Phaser if anisotropy and tNCS are properly accounted for using strategies that are applied automatically in Phaser v.2.7.17 or newer. Moreover, the structure of the hRab27aMut(GppNHp) crystals shows that the SF4 pocket, which is the primary target for ligand-binding studies, is unoccupied and could be used to study the structure of ligands binding to Rab27a. The only major drawback is the data quality, specifically the overall resolution and severe anisotropy, which would be problematic for weak binding ligands with low occupancy. Optimization of crystallization conditions, additive screens and the structure of hRab27aMut(GppNHp) reported here will guide further construct design to obtain a more tractable crystal form for ligand-binding studies.
Acknowledgements
We thank Diamond Light Source for access to beamline I02 (proposal No. mx9424), which contributed to the results presented here, and Dr Marc Morgan of the Imperial College London X-ray facility for his assistance throughout the project.
Funding information
The research was supported by Cancer Research UK (Drug Discovery Committee grant C29637/A20781 to EC and EWT) and the award of a Principal Research Fellowship to RJR by the Wellcome Trust (grant 082961/Z/07/Z).
References
Berman, H., Henrick, K., Nakamura, H. & Markley, J. L. (2007). Nucleic Acids Res. 35, D301–D303. Web of Science CrossRef PubMed CAS Google Scholar
Bishop, C. M. (2006). Pattern Recognition and Machine Learning. New York: Springer. Google Scholar
Brozzi, F., Diraison, F., Lajus, S., Rajatileka, S., Philips, T., Regazzi, R., Fukuda, M., Verkade, P., Molnár, E. & Váradi, A. (2012). Traffic, 13, 54–69. CrossRef CAS Google Scholar
Bunkóczi, G. & Read, R. J. (2011a). Acta Cryst. D67, 303–312. CrossRef IUCr Journals Google Scholar
Bunkóczi, G. & Read, R. J. (2011b). Comput. Crystallogr. Newsl. 2, 8–9. Google Scholar
Chavas, L. M., Ihara, K., Kawasaki, M., Torii, S., Uejima, T., Kato, R., Izumi, T. & Wakatsuki, S. (2008). Structure, 16, 1468–1477. CrossRef CAS Google Scholar
Chavas, L. M. G., Ihara, K., Kawasaki, M. & Wakatsuki, S. (2010). Nihon Kessho Gakkaishi, 51, 320–326. Google Scholar
Chavas, L. M. G., Torii, S., Kamikubo, H., Kawasaki, M., Ihara, K., Kato, R., Kataoka, M., Izumi, T. & Wakatsuki, S. (2007). Acta Cryst. D63, 769–779. Web of Science CrossRef IUCr Journals Google Scholar
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 Google Scholar
Dong, W.-W., Mou, Q., Chen, J., Cui, J.-T., Li, W.-M. & Xiao, W.-H. (2012). World J. Gastroenterol. 18, 1806–1813. CrossRef CAS Google Scholar
Drenth, J. (2007). Principles of Protein X-ray Crystallography, 3rd ed. New York: Springer. Google Scholar
Dumas, J. J., Zhu, Z., Connolly, J. L. & Lambright, D. G. (1999). Structure, 7, 413–423. Web of Science CrossRef PubMed CAS Google Scholar
Emsley, P., Lohkamp, B., Scott, W. G. & Cowtan, K. (2010). Acta Cryst. D66, 486–501. Web of Science CrossRef CAS IUCr Journals Google Scholar
Evans, P. R. (2011). Acta Cryst. D67, 282–292. Web of Science CrossRef CAS IUCr Journals Google Scholar
Evans, P. R. & Murshudov, G. N. (2013). Acta Cryst. D69, 1204–1214. Web of Science CrossRef CAS IUCr Journals Google Scholar
French, S. & Wilson, K. (1978). Acta Cryst. A34, 517–525. CrossRef CAS IUCr Journals Web of Science Google Scholar
Fukuda, M. (2013). Traffic, 14, 949–963. CrossRef CAS Google Scholar
Goishi, K., Mizuno, K., Nakanishi, H. & Sasaki, T. (2004). Biochem. Biophys. Res. Commun. 324, 294–301. CrossRef CAS Google Scholar
Guo, Z., Hou, X., Goody, R. S. & Itzen, A. (2013). J. Biol. Chem. 288, 32466–32474. CrossRef CAS Google Scholar
Kantardjieff, K. A. & Rupp, B. (2003). Protein Sci. 12, 1865–1871. Web of Science CrossRef PubMed CAS Google Scholar
Kukimoto-Niino, M., Sakamoto, A., Kanno, E., Hanawa-Suetsugu, K., Terada, T., Shirouzu, M., Fukuda, M. & Yokoyama, S. (2008). Structure, 16, 1478–1490. CAS Google Scholar
Kullback, S. & Leibler, R. A. (1951). Ann. Math. Stat. 22, 79–86. CrossRef Web of Science Google Scholar
Li, W., Hu, Y., Jiang, T., Han, Y., Han, G., Chen, J. & Li, X. (2014). APMIS, 122, 1080–1087. CAS Google Scholar
Luzzati, V. (1952). Acta Cryst. 5, 802–810. CrossRef IUCr Journals Web of Science Google Scholar
Matthews, B. W. (1968). J. Mol. Biol. 33, 491–497. CrossRef CAS PubMed Web of Science Google Scholar
McCoy, A. J. (2007). Acta Cryst. D63, 32–41. Web of Science CrossRef CAS IUCr Journals Google Scholar
McCoy, A. J., Grosse-Kunstleve, R. W., Adams, P. D., Winn, M. D., Storoni, L. C. & Read, R. J. (2007). J. Appl. Cryst. 40, 658–674. Web of Science CrossRef CAS IUCr Journals Google Scholar
McCoy, A. J., Grosse-Kunstleve, R. W., Storoni, L. C. & Read, R. J. (2005). Acta Cryst. D61, 458–464. Web of Science CrossRef CAS IUCr Journals Google Scholar
McCoy, A. J., Oeffner, R. D., Wrobel, A. G., Ojala, J. R. M., Tryggvason, K., Lohkamp, B. & Read, R. J. (2017). Proc. Natl Acad. Sci. USA, 114, 3637–3641. Web of Science CrossRef CAS PubMed Google Scholar
Oeffner, R. D., Bunkóczi, G., McCoy, A. J. & Read, R. J. (2013). Acta Cryst. D69, 2209–2215. Web of Science CrossRef CAS IUCr Journals Google Scholar
Ostrowski, M., Carmo, N. B., Krumeich, S., Fanget, I., Raposo, G., Savina, A., Moita, C. F., Schauer, K., Hume, A. N., Freitas, R. P., Goud, B., Benaroch, P., Hacohen, N., Fukuda, M., Desnos, C., Seabra, M. C., Darchen, F., Amigorena, S., Moita, L. F. & Thery, C. (2010). Nature Cell Biol. 12, 19–30. CrossRef CAS Google Scholar
Read, R. J., Adams, P. D. & McCoy, A. J. (2013). Acta Cryst. D69, 176–183. Web of Science CrossRef CAS IUCr Journals Google Scholar
Read, R. J. & McCoy, A. J. (2016). Acta Cryst. D72, 375–387. Web of Science CrossRef IUCr Journals Google Scholar
Sliwiak, J., Jaskolski, M., Dauter, Z., McCoy, A. J. & Read, R. J. (2014). Acta Cryst. D70, 471–480. Web of Science CrossRef CAS IUCr Journals Google Scholar
Storoni, L. C., McCoy, A. J. & Read, R. J. (2004). Acta Cryst. D60, 432–438. Web of Science CrossRef CAS IUCr Journals Google Scholar
Strom, M., Hume, A. N., Tarafder, A. K., Barkagianni, E. & Seabra, M. C. (2002). J. Biol. Chem. 277, 25423–25430. CrossRef CAS Google Scholar
Strong, M., Sawaya, M. R., Wang, S., Phillips, M., Cascio, D. & Eisenberg, D. (2006). Proc. Natl Acad. Sci. USA, 103, 8060–8065. Web of Science CrossRef PubMed CAS Google Scholar
Ten Eyck, L. F. (1973). Acta Cryst. A29, 183–191. CrossRef CAS IUCr Journals Web of Science Google Scholar
Vagin, A. & Teplyakov, A. (2010). Acta Cryst. D66, 22–25. Web of Science CrossRef CAS IUCr Journals Google Scholar
Wang, J.-S., Wang, F.-B., Zhang, Q.-G., Shen, Z.-Z. & Shao, Z.-M. (2008). Mol. Cancer Res. 6, 372–382. CrossRef CAS Google Scholar
Wang, Q., Ni, Q., Wang, X., Zhu, H., Wang, Z. & Huang, J. (2015). Med. Oncol. 32, 372. CrossRef Google Scholar
Waterman, D. G., Winter, G., Gildea, R. J., Parkhurst, J. M., Brewster, A. S., Sauter, N. K. & Evans, G. (2016). Acta Cryst. D72, 558–575. Web of Science CrossRef IUCr Journals Google Scholar
Whittaker, E. T. & Watson, G. N. (1990). A Course in Modern Analysis, 4th ed., pp. 347–348. Cambridge University Press. Google Scholar
Wilson, A. J. C. (1949). Acta Cryst. 2, 318–321. CrossRef IUCr Journals Web of Science Google Scholar
Winn, M. D., Ballard, C. C., Cowtan, K. D., Dodson, E. J., Emsley, P., Evans, P. R., Keegan, R. M., Krissinel, E. B., Leslie, A. G. W., McCoy, A., McNicholas, S. J., Murshudov, G. N., Pannu, N. S., Potterton, E. A., Powell, H. R., Read, R. J., Vagin, A. & Wilson, K. S. (2011). Acta Cryst. D67, 235–242. Web of Science CrossRef CAS IUCr Journals Google Scholar
Winn, M. D., Murshudov, G. N. & Papiz, M. Z. (2003). Methods Enzymol. 374, 300–321. Web of Science CrossRef PubMed CAS Google Scholar
Winter, G. (2010). J. Appl. Cryst. 43, 186–190. Web of Science CrossRef CAS IUCr Journals Google Scholar
Winter, G., Waterman, D. G., Parkhurst, J. M., Brewster, A. S., Gildea, R. J., Gerstel, M., Fuentes-Montero, L., Vollmar, M., Michels-Clark, T., Young, I. D., Sauter, N. K. & Evans, G. (2018). Acta Cryst. D74, 85–97. Web of Science CrossRef IUCr Journals Google Scholar
Yamaoka, M., Ishizaki, T. & Kimura, T. (2015). World J. Diabetes, 6, 508–516. CrossRef Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.