research papers
Joint neutron/X-ray
of a mechanistically relevant complex of perdeuterated urate oxidase and simulations provide insight into the hydration step of catalysisaRandall Centre for Cell and Molecular Biophysics, King's College London, London SE1 1UL, United Kingdom, bLarge Scale Structures Group, Institut Laue-Langevin, 71 avenue des Martyrs, 38042 Cedex 9, Grenoble, France, cDepartment of Chemistry, King's College London, London SE1 1DB, United Kingdom, dDepartment of Physics and Astronomy, University College, London WC1E 6BT, United Kingdom, and eLife Sciences Group, Institut Laue-Langevin, 71 avenue des Martyrs, 38042 Cedex 9, Grenoble, France
*Correspondence e-mail: roberto.steiner@kcl.ac.uk
Cofactor-independent urate oxidase (UOX) is an ∼137 kDa tetrameric enzyme essential for uric acid (UA) catabolism in many organisms. UA is first oxidized by O2 to dehydroisourate (DHU) via a peroxo intermediate. DHU then undergoes hydration to 5-hydroxyisourate (5HIU). At different stages of the reaction both catalytic O2 and water occupy the `peroxo hole' above the organic substrate. Here, high-resolution neutron/X-ray crystallographic analysis at room temperature has been integrated with simulations to investigate the hydration step of the reaction. The joint neutron/X-ray structure of perdeuterated Aspergillus flavus UOX in complex with its 8-azaxanthine (8AZA) inhibitor shows that the catalytic water molecule (W1) is present in the peroxo hole as neutral H2O, oriented at 45° with respect to the ligand. It is stabilized by Thr57 and Asn254 on different UOX protomers as well as by an O—H⋯π interaction with 8AZA. The active site Lys10–Thr57 dyad features a charged Lys10–NH3+ side chain engaged in a strong hydrogen bond with Thr57OG1, while the Thr57OG1–HG1 bond is rotationally dynamic and oriented toward the π system of the ligand, on average. Our analysis offers support for a mechanism in which W1 performs a nucleophilic attack on DHUC5 with Thr57HG1 central to a Lys10-assisted proton-relay system. Room-temperature crystallography and simulations also reveal conformational heterogeneity for Asn254 that modulates W1 stability in the peroxo hole. This is proposed to be an active mechanism to facilitate W1/O2 exchange during catalysis.
Keywords: neutron/X-ray diffraction; urate oxidase; protein perdeuteration; cofactor-independent oxidase; biomolecular simulations.
1. Introduction
In many organisms, urate oxidase (UOX) is an essential enzyme that catalyses the O2-dependent degradation of uric acid (UA) to 5-hydroxyisourate (5HIU) (Kahn et al., 1997). Humans, however, lack UOX, as it was lost by pseudogenization during primate evolution resulting in inactivation of the uricolytic pathway (Oda et al., 2002). Therefore, recombinant UOX finds a therapeutic application for the enzymatic treatment of severe hyperuricemia and prevents the UA burst accompanying tumour lysis after certain chemotherapy treatments (Terkeltaub, 2010).
Mechanistically, UOX-mediated degradation of UA follows two sequential steps: (1) an initial oxidation step whereby UA reacts with O2 to yield dehydroisourate (DHU) via a 5-peroxoisourate (5PIU) intermediate and (2) a hydration step in which DHU is hydroxylated to 5HIU [Fig. 1(a)] (Kahn, 1999; Wei et al., 2017). The latter is ultimately transformed to the more soluble allantoin either enzymatically or non-enzymatically (Ramazzina et al., 2006). High-resolution X-ray crystallographic studies of the anaerobic complex of Aspergillus flavus UOX (hereafter simply UOX) in complex with its natural substrate UA (Bui et al., 2014) as well as with several other similar scaffolds including the 8-azaxanthine (8AZA) inhibitor employed in this study [Fig. 1(b)] (Colloc'h et al., 1997; Gabison et al., 2010; Retailleau et al., 2004), reveal that ligands bind at the interface between two protomers of the ∼137 kDa UOX tetramer stabilized by several hydrogen bonds with the additional contribution of a π–π stacking interaction [Figs. 1(c) and 1(d)]. Directly above the ligand plane, a water molecule (W1) suggested to be responsible for the hydration step of catalysis (Wei et al., 2017) is coordinated by the side chains of Thr57* and Asn254 (an asterisk indicates residues belonging to the second UOX molecule contributing to the active site). A comparison with the X-ray structure of UOX in complex with the 5PIU intermediate trapped under cryo-conditions (Bui et al., 2014) shows that W1 occupies the same position of the peroxo group [Fig. S1(a) of the supporting information]. This location (the `peroxo hole') also provides stabilization to O2 during the formation of the 5PIU intermediate (Bui et al., 2014) [Fig. S1(b)]. Small negative ions like cyanide and chloride have been visualized in the peroxo hole in the presence of intact UA under aerobic conditions (Gabison et al., 2008, 2010; Oksanen et al., 2014) offering structural support for their behaviour as competitive inhibitors.
UOX belongs to a group of intriguing cofactor-independent oxidases and oxygenases that rely on their organic substrates for O2 activation (Fetzner & Steiner, 2010). At optimal catalytic pH (∼8.0) (Pfrimer et al., 2010), UA is a monoanion in solution and an early mechanistic proposal suggested the existence of a UOX general base dyad formed by Lys10* and Thr57* which further deprotonate UA upon binding, generating the dianion (Imhoff et al., 2003). A signature for the latter has been identified spectroscopically during catalysis (Kahn & Tipton, 1998). As deduced by its tendency to undergo spontaneous oxidation in the gas phase (Kahn & Tipton, 1997), the UA dianion has been proposed to transfer one electron to O2 to generate the superoxide anion (). This would result in a [] pair that, following radical recombination, yields 5PIU. However, a recent theoretical study questioned the involvement of the UA dianion (Wei et al., 2017). Instead, it proposed a in which the Lys10*–Thr57* dyad abstracts the H8 proton from the UA monoanion bound in its C8-lactim form [tautomer shown in Fig. 1(a)] with a concomitant one-electron transfer to O2 resulting directly in the formation of the [] pair. The latter mechanism is supported by the neutron crystallographic structure of the UOX–UA complex solved under aerobic conditions in the presence of inhibitory chloride in the peroxo hole (UOX–UA–Cl−, PDB entry 4n9m; Oksanen et al., 2014) which revealed UA bound in its monoanion form deprotonated at N3 with a deuteron bound to O8. Deprotonation at N3 was also observed for 8AZA in the neutron structure of the UOX–8AZA–Cl− complex (PDB entry 4n3m; Oksanen et al., 2014).
Displacement of W1 by chloride in the neutron structures currently available precludes obtaining direct information on the catalytic water and on the protonation state of its environment relevant to the hydration step of catalysis. Here, we have integrated high-resolution room-temperature neutron/X-ray crystallographic structural analysis with DUOX) instead of relying on protium (1H)/deuterium (2H or D) replacement of exchangeable hydrogens as done in previous studies on this system (Oksanen et al., 2014). Protein perdeuteration, which replaces all hydrogen atoms with deuterium, is an effective method of improving S/N as it dramatically lowers the incoherent background while enhancing the coherent scattering signal (Fisher et al., 2014; Shu et al., 2000). In addition, perdeuteration avoids map cancellation issues, allowing all D atoms to be visualized at ∼2.5 Å resolution or better (Blakeley & Podjarny, 2018; Logan, 2020).
(MD) simulations to study the UOX–8AZA complex in the presence of W1 in the peroxo hole. To maximize the signal-to-noise ratio (S/N) of the neutron diffraction data we employed perdeuterated UOX (The high-resolution joint neutron/X-ray structure of DUOX–8AZA–W1 presented here provides experimental support for a mechanism in which W1 attacks DHUC5 (atom C5 of DHU) with the assistance of the Lys10*–NH3+⋯HO–Thr57* dyad in which Thr57*HG1 is central to a proton-relay system. Experiment and simulations highlight a degree of rotational freedom for the Thr57* hydroxyl group. We suggest that this helps Thr57* to fulfil a dual role of stabilizing W1 as well as acting as a proton donor to the DHUN7 π system to kickstart hydroxylation. Moreover, experiment and simulations also reveal an unforeseen conformational heterogeneity for Asn254 that modulates W1 stability in the peroxo hole. This might represent an active mechanism that facilitates W1/O2 exchange during catalysis.
2. Materials and methods
2.1. Perdeuterated A. flavus UOX (DUOX)
A. flavus urate oxidase (UOX) cDNA codon-optimized for bacterial expression was purchased from Genscript (Piscataway, NJ, USA) and cloned using the NdeI/XhoI restriction sites in a pET24b(+) expression vector (Novagen). Recombinant untagged perdeuterated UOX (DUOX) was expressed in E. coli BL21 (DE3) cells. Cells were adapted for growth under deuterated conditions via a stepwise adaptation process using a minimal medium with glycerol-d8 as the carbon source. DUOX was produced in a large quantity via batch-fed (Haertlein et al., 2016). DUOX was purified using a multi-step protocol previously developed for protiated UOX (Bui et al., 2014). Briefly, pelleted cells were re-suspended in 50 mM Tris–HCl (pH 8.0), 250 mM NaCl, supplemented with lysozyme, DNAse and a protease inhibitor cocktail, and lysed by sonication. Protein purification was performed using a combination of ammonium sulfate precipitation, DEAE and Resource Q ion exchange, Phenyl Sepharose hydrophobic interaction, and Superdex 75 size-exclusion chromatographic steps. After each step samples were analysed by SDS-PAGE and the purest fractions were pooled for further work. Both hydrophobic interaction and size-exclusion purification steps were carried out in the absence of NaCl in the buffers.
2.2. Mass spectrometry
M with solution A (0.03% TFA in water). Protein samples were firstly desalted on a reverse-phase C8 cartridge (Zorbax 300SB-C8, 5 µm, 300 µm ID × 5 mm, Agilent Technologies) for 3 min at a flow rate of 50 µl min−1 with 100% solvent A and then eluted with 70% solvent B at flow rate of 50 µl min−1 for MS detection. MS acquisition was carried out in positive ion mode in the 300–3200 m/z range. MS spectra were acquired and the data processed with the MassHunter workstation software (v.B.02.00, Agilent Technologies).
electrospray ionization (LC/ESI-MS) was performed on a 6210 LC-TOF spectrometer coupled to an HPLC system (Agilent Technologies). All solvents used were HPLC grade (Chromasolv, Sigma–Aldrich), trifluoroacetic acid (TFA) was from Acros Organics (puriss p.a.). Solvent A was 0.03% TFA in water, solvent B was 95% acetonitrile, 5% water and 0.03% TFA. Just before analysis protein samples were diluted under acidic denaturing conditions to a final concentration of 5 µ2.3. Crystallization, neutron/X-ray data collection and structural analysis
DUOX eluted from in 50 mM Tris–acetate, 30 mM sodium acetate at pH 8.0 was thoroughly exchanged against 50 mM Tris–acetate, 30 mM sodium acetate in D2O, pD 7.59, before concentrating to >20 mg ml−1. DUOX–8AZA crystals were grown at 291 K using the batch crystallization method. Deuterated crystallization buffer consisted of 50 mM Tris–acetate, 8% PEG 8000, 2 mg ml−1 8AZA, pD 7.59. Crystallization drops of 30 µl total volume were prepared with protein solution and crystallization buffer in a 2:1 ratio. After ∼10 d, a single crystal of approximate volume 1.5 mm3 (1.4 × 1.2 × 0.9 mm) was mounted in a quartz capillary and sealed with wax for data collection.
Neutron quasi-Laue diffraction data to 2.10 Å resolution were collected at 298 K on the LADI-III instrument at the Institut Laue-Langevin (ILL), Grenoble, France (Blakeley et al., 2010). Data were collected using the wavelength range 3.12–4.20 Å, centred at 3.68 Å, and an exposure time of 1.5 h per image. Neutron diffraction data were indexed and integrated with the program LAUEGEN (Campbell, 1995). LSCALE (Arzt et al., 1999) was used in the wavelength normalization step, followed by data scaling and merging using SCALA (Evans, 2006), and conversion to structure factors using the program TRUNCATE (French & Wilson, 1978; Winn et al., 2011). Neutron data collection statistics are given in Table 1. Subsequent to the neutron diffraction experiment, the same crystal was used for X-ray diffraction data collection at 298 K on the BM30A (FIP) beamline at the European Synchrotron Radiation Facility (ESRF), Grenoble, France. Data to 1.33 Å resolution were indexed and integrated using the xia2 pipeline (Winter, 2010). Statistics for the X-ray dataset are also reported in Table 1.
|
Joint neutron/X-ray phenix.refine (Afonine et al., 2010) package of the Phenix suite (Liebschner et al., 2019) using the anisotropic treatment of ADPs for all non-hydrogen atoms. Deuterium atoms, added to the model at hydrogen positions using the READYSET program (Liebschner et al., 2019) were refined isotropically. Model fitting to X-ray and neutron maps was carried out with the package Coot (Emsley et al., 2010). can be found in Table 2. The final model was deposited in the PDB (PDB entry 7a0l).
was performed with the
|
2.4. simulations
NAMD 2.12 software (Phillips et al., 2005, 2020). Parametrization of the system, including the 8AZA ligand, was performed through CHARMM-GUI (Jo et al., 2008). The electrostatic description of the ligand was enhanced by ESP partial atomic charges, calculated separately with Gaussian09 (Frisch et al., 2013) at the B3LYP-BJD3/6-311G(d,p) level of theory (Becke, 1993; Grimme et al., 2011; McLean & Chandler, 1980; Krishnan et al., 1980) with the CHELPG algorithm (Breneman & Wiberg, 1990). The simulated systems were solvated in a TIP3P water box with the box size based on the dimensions of the system and giving 10 Å extension in all directions. Four K+ counterions were added to neutralize the system. For each system, our MD protocol consisted of the following steps: (1) energy minimization over 10 000 steps using a timestep of 2 fs; (2) equilibration over 1 ns in the NVT ensemble (T = 303.15 K) with the RMSD of heavy atoms of the host and guest molecules constrained to their initial position using a force constant of 1 kcal mol−1 Å−2; (3) production run in NPT ensemble with a time step of 2 fs, where the temperature and pressure were held at constant at 303.15 K and 1.0 atm, respectively. Constant temperature was set by a Langevin thermostat with a of 1 ps−1, whereas the target pressure was reached using the Berendsen barostat. The Shake algorithm (Andersen, 1983) was used for bonds and angles involving hydrogen atoms. We used the particle mesh Ewald method (Darden et al., 1993) for the long-range electrostatics in combination with a cutoff of 12 Å for the evaluation of the non-bonded interactions.
(MD) simulations were performed with the2.5. Quantum mechanics/molecular mechanics calculations
Quantum mechanics/molecular mechanics (QM/MM) calculations, including geometry minimizations and CHARMM molecular mechanics simulation software (Brooks et al., 2009). Parametrization of the system was performed through CHARMM-GUI (Jo et al., 2008). A cutoff distance of 12 Å was used for the Lennard–Jones and electrostatic interactions. In order to keep to zero, shifting was applied to the pairwise interactions over 10 Å distances. Residues further than 25 Å from the ligand were deleted from the structures. We applied constraints on the atoms lying further than 15 Å from a ligand to keep them at their original positions. These geometry restrictions were defined for the starting structure and kept unaltered through the calculations. The forces acting on the QM region of the model were calculated externally with the Q-Chem 4.3 software suite (Shao et al., 2015), which also involved the surrounding atoms as point charges. DFT calculations were performed at the B3LYP/6–31+G* level of theory. For geometry minimizations a limit of 1000 steps and a threshold for the RMSD of the gradient of 0.005 kcal mol−1 was applied. The connection between the QM and the MM regions was maintained with hydrogen link atoms.
simulations, were performed with3. Results
3.1. Perdeuterated A. flavus UOX (DUOX) for neutron crystallography
A. flavus UOX expressed recombinantly in Saccharomyces cerevisiae is a therapeutic commercial product of Sanofi–Aventis under the brand name Fasturtec. Several high-resolution structural studies on this enzyme relied on protein material kindly gifted by the company (Gabison et al., 2010; Oksanen et al., 2014; Retailleau et al., 2004). To circumvent this limitation we previously set up a heterologous expression system in E. coli (Bui et al., 2014). Using codon-optimized cDNA, untagged UOX can be expressed in good quantities and purified to allowing us to obtain crystals that diffract X-rays to atomic resolution (Bui et al., 2014). Bacterial expression of UOX from other organisms has also been reported enabling X-ray crystallographic studies to be performed, albeit at lower resolution (Hibi et al., 2016; Marchetti et al., 2016; Imhoff et al., 2003).
As neutron crystallographic studies benefit substantially from data collection on fully deuterated crystals we took advantage of the D-LAB facility (ILL, Grenoble, France) to produce perdeuterated UOX (DUOX) using our expression construct. Stepwise adaptation of E. coli to 100% D2O and the use of deuterated glycerol as a carbon source allowed us to obtain high levels of DUOX after optimization of conditions (Haertlein et al., 2016). Purification of DUOX mirrored the strategy we established previously for (protiated) UOX affording more than 50 mg of pure DUOX per litre of culture [Fig. 2(a)]. The degree of deuteration was estimated by MS at >99% [Figs. 2(b) and 2(c)].
3.2. Joint neutron/X-ray structure of the DUOX–8AZA–W1 complex
A single crystal (volume ≃ 1.5 mm3) of DUOX in complex with 8AZA obtained under chloride-free conditions was utilized to collect sequentially room-temperature neutron and X-ray diffraction data at resolutions of 2.10 and 1.33 Å, respectively. Data collection statistics are given in Table 1. As the crystal was both perdeuterated and had a large volume, relatively short exposure times could be used for neutron data collection, allowing the complete dataset to be collected in only 16.5 h at the LADI-III beamline at the ILL (Blakeley et al., 2010). Typically, neutron data collections take several days (Blakeley & Podjarny, 2018) and so the data collection for our study described here represents (as far as we are aware) the fastest neutron data collection for a protein with a relatively large For X-ray data collection we did not attempt to measure data at the highest possible resolution. Instead, we opted to collect a complete dataset with negligible radiation damage. Electron density maps display none of the specific signatures typically associated with radiation damage, such as decarboxylation of acidic residues (Ravelli & McSweeney, 2000; Weik et al., 2000). In addition, neutron and X-ray datasets are isomorphous ruling out any obvious unit-cell expansion during X-ray data collection, a global indicator of radiation damage (Garman, 2010).
As typically observed for protiated UOX crystals, the DUOX–8AZA–W1 complex crystallizes in the orthorhombic I222. The contains a single DUOX chain that generates the biological tetramer of D2 symmetry following the application of operations. Joint neutron/X-ray delivered a final model characterized by R/Rfree values of 21.47%/23.03% and 10.16%/10.93% for neutron and X-ray data, respectively. DUOX residues are visible for residues 1–295, whereas the six C-terminal residues are flexible. In addition to the 8AZA ligand, 259 solvent molecules are modelled as D2O and 92 additional ones are modelled as oxygen atoms only owing to rotational disorder. are summarized in Table 2. The use of perdeuterated UOX results in nuclear scattering length density free from protium-induced cancellation effects which allows the visualization of both exchangeable and non-exchangeable hydrogens. Representative examples of the excellent quality of X-ray and neutron maps are given in Figs. 2(d), 2(e) and S2.
A superposition of the final room-temperature DUOX–8AZA–W1 model with that of the 1.75 Å resolution UOX–8AZA–W1 X-ray structure solved at 4°C (277 K) (PDB entry 1r51; Retailleau et al., 2004) or with that of H/D-exchanged 1.1 Å resolution H/DUOX–8AZA–Cl− X-ray structure solved at 100 K (PDB entry 4n9v; Oksanen et al., 2014) reveals essentially no differences, with RMSD values of 0.18 and 0.22 Å, respectively, for all common main-chain atoms. Thus, in line with what has been reported for other proteins (Koruza et al., 2019), partial or complete deuterium-labelling does not induce major structural changes in UOX.
3.3. AZA and W1
As observed in previous UOX complexes (Colloc'h et al., 1997, 2007; Retailleau et al., 2004; Oksanen et al., 2014), 8AZA binds in the active site at the interface between two UOX protomers with direct hydrogen-bond stabilizing contributions from Thr57*, Asp58*, Arg176, Val227 and Gln228, while residues Tyr8*, Ile54*, Ala56*, Phe159, Leu170 and Ser226 further line the binding pocket with Phe159 engaged in π–π stacking with the pyrimidine ring of 8AZA [Fig. 3(a)]. We observe dual conformations for the side chains of Tyr8* and Ile54*, as well as for Ile288 located directly above Gln228.
8AZA binds as the N3 monoanion with clear nuclear scattering length density for deuterons linked to N1 and N9. Although nuclear scattering length density in omit maps appeared somewhat stronger for the hydrogen bound to N1, occupancy D9 (8AZAD9⋯W2O = 1.84 Å), though there is no evidence for a solvent molecule interacting with 8AZAN8 as seen in the anaerobic complexes with UA and its 9-methyl derivative (9-MUA), both of which can react with O2 to form C5-peroxides (Bui et al., 2014).
gives unitary values for both deuterons. A water molecule (W2) forms a strong hydrogen bond with 8AZAW1 is visualized in the electrostatically positive peroxo hole above the 8AZA plane and with its oxygen atom facing Thr57* [Fig. 3(b)]. Electron density maps clearly define the position of the W1 oxygen (W1O), whereas omit neutron maps calculated by modelling W1 as W1O only show elongated positive difference density next to the oxygen, which is consistent with the presence of two deuterium atoms [Fig. 3(c)]. In the final model W1 refines to an overall occupancy of 0.93. It is sandwiched between Thr57* and Asn254 with its plane oriented at ∼45° with respect to that of 8AZA with its O—D1 bond almost parallel to it and directed into a small cavity lined by 8AZA, Lys10*, Thr57*, Asn254 and His256 [Fig. 3(d)]. W1O is hydrogen bonded to Thr57*DG1 at a distance of 2.16 Å. Similar distances (2.23–2.25 Å) are also observed between W1O and the amide deuterons of Asn254, which, similar to other residues lining the active site (Tyr8*, Ile54*, Ile288), also exhibit alternative conformations. The visualization of W1 hydrogens allows us to also define an unanticipated O—H⋯π interaction with the ligand. Following the classical parametrization for this type of bond (Steiner & Koellner, 2001), which uses the distance (dO⋯m) between the donor atom (O) and the midpoint of the ring (m), the ω(O) angle between the O⋯m line and the ring normal as well as the O—H⋯m angle [Fig. 3(d)], we find values of dO⋯m = 3.64 Å, ω(O) = 23.8° and ∠(O—H⋯m) = 141° that are consistent with those most commonly observed for the peptide n → (n − 2) N—H⋯π interaction for which there are abundant data in proteins (Steiner & Koellner, 2001).
3.4. Lys10*–Thr57* dyad
The active-site dyad formed by Lys10*–Thr57* has been implicated in various stages of the catalytic cycle (Gabison et al., 2010; Imhoff et al., 2003; Wei et al., 2017), thus its protonation state was carefully investigated. While difference neutron maps for Lys10* showed a clear `tri-lobe' nuclear density distribution around NZ that is consistent with a positively charged (Lys10*–ND3+) group [Fig. 4(a)], the identification of the protonation state for Thr57* was not so straightforward. Difference neutron maps did not provide evidence for a positive peak adjacent to Thr57*OG1 at the +3.0σ level but positive density that can be attributed to DG1 was seen at the +2.0σ level [Fig. 4(a)]. As neutron omit maps offered only a mild suggestion for the presence of a deuteron bound to Thr57*OG1 and difference map interpretation at low threshold levels need to be considered very carefully, we decided to investigate the unlikely possibility of a Lys10*–NH3+⋯O−–Thr57* ionic pair using theoretical methods. However, this scenario was quickly ruled out as already during the first few frames of QM/MM minimization (order of a few femtoseconds) a spontaneous proton transfer occurs from the Lys10* amine to the Thr57* hydroxylate resulting in the establishment of a Lys10*–NH2⋯HO–Thr57* hydrogen bond (Fig. S3). This strongly argues against the deprotonation of Thr57*.
We then considered the possibility that Thr57* is in the neutral state and that the lack of clear density for DG1 in neutron difference maps can be attributed to positional heterogeneity which cannot be interpreted fully at the resolution of our neutron diffraction data (2.10 Å). To further explore this, we calculated the distribution of the CA–CB–OG1–HG1 torsion angle during a 100 ns classical MD simulation. As shown in Fig. 4(b), this dihedral exhibits a rather broad asymmetric profile ranging from −10 to −190° with a maximum at around −72°, resulting in the OG1–HG1 bond pointing `down' towards 8AZA, almost perpendicular to it. A smaller `shoulder' is also seen around −135°, with OG1–HG1 parallel to 8AZA. We therefore modelled a deuteron atom bound to Thr57*OG1 as suggested by omit neutron maps. Occupancy for this atom gives a value of 0.78. In the final model the CA–CB–OG1–DG1 torsion angle displays a value of −110.4° resulting in a slight downward orientation of the OG1–DG1 bond towards the ligand plane with the deuteron at 2.77 Å from 8AZAN7, thus indicating another O—H⋯π interaction [Fig. 4(c)]. All deuterons bound to Lys10*NZ also refine to high occupancy (>0.90) with a hydrogen bond between Lys10*DZ1 and Thr57*OG1 at a distance of 1.97 Å. Lys10* is in turn hydrogen bonded to His256. The latter is neutral, with its hydrogen residing on ND1. As a result, His256NE2 accepts a hydrogen bond from the Lys10*DZ3, positioned at a distance of 2.24 Å. Lys10*DZ2 is hydrogen bonded to water W4 that, together with W5 and W2, completes a chain of solvent molecules connecting the Lys10*–Thr57* dyad to AZAD9. Overall, experimental data and theoretical calculations are consistent with the existence of a Lys10*–ND3+⋯DO–Thr57* dyad in which the Thr57*OG1–DG1 bond is directed towards 8AZA, on average.
3.5. Heterogeneity of Asn254 and its flanking residues
During crystallographic β7 and β7′ strands. This indicated the existence of an alternative conformation in addition to the main one typically observed, which we could model [Fig. 5(a)]. Occupancy for this short stretch that hosts Asn254 suggests an approximate 75:25 ratio for the major [yellow in Fig. 5(a)] and minor (dark grey) conformations. In the latter, Pro253 is shifted towards β8 resulting in a movement of Asn254 in the opposite direction. This leads to an Asn254CA displacement of ∼1.4 Å and the consequent movement of its side chain that nevertheless remains at hydrogen-bond distance from W1. Modelling of this region in an alternative conformation made a residual peak more obvious in mFo−DFc electron density maps at the +4.0σ level (∼0.25 e− Å−3) near the amide group of Asn254 in its major conformation (Fig. S4). We interpreted this peak as a low-occupancy water molecule (oxygen atom only, refined occupancy = 0.21). Inclusion of this solvent molecule in the model fully accounted for the positive difference density. Moreover, we calculated 2mFo−DFc feature enhanced maps (FEM) that can strengthen weak signals, if present, and can reduce model bias and noise (Afonine et al., 2015). The FEM map for this region shows clearly the alternative conformation for Pro253–Lys255 peptide as well as a `bulge' on the side of the amide plane that supports our modelling of this low-occupancy water (W′) [Fig. 5(b)].
inspection of near-atomic resolution X-ray difference maps showed a pattern of positive and negative peaks for the tripeptide linker Pro253–Lys255 connecting theTo further validate our experimental observation of flexibility for this region we looked at the distribution of Asn254 torsion angles using MD simulations [Fig. 5(c)]. Two independent 100 ns simulations were carried out using a UOX dimer model (thus allowing the analysis of two active sites in each run) (Fig. S5) with the Pro253–Lys255 stretch either in the major or minor conformation. The initial 25 ns of each simulation were removed from the analysis to avoid bias caused by system conditioning. As shown in Fig. 5(c), the MD simulations allow us to define four main cluster areas (green, yellow, aqua, dark grey) in the (CA–CB–CG–ND2)–(N–CA–CB–CG) 2D contour plot. The cluster highlighted in yellow is representative of the main conformation seen in the crystal [red square at 179°, 192° in Fig. 5(c)] and confirms that this rotamer is sampled often during the simulation. A snapshot from the simulation illustrative of this cluster is given in Fig. 5(d) and shows an arrangement similar to that seen in the crystal. The dark grey cluster representing the minor crystallographic conformation (red square at 272°, 81°) is also occasionally visited, thus lending support to our structural observation. In many structures of this cluster W1 moves away from the peroxo hole with the hydrogen of the Asn254 amide replacing W1 in an N—H⋯π interaction with 8AZA [Fig. 5(e)]. This novel W1 position appears consistent with that of that of the low-occupancy W′ in our model.
The simulations also reveal additional torsion angles frequently adopted by Asn254 that define almost a continuum (aqua and green clusters). Representative structures are shown in Figs. 5(f), and 5(g), respectively. Structures belonging to the aqua cluster display the Asn254 side-chain carbonyl group pointing towards the peroxo hole typically resulting in the displacement of W1 and the establishment of a solvent chain organized above the ligand [Fig. 5(f)]. The transition from the aqua to the green cluster involves a 90° rotation of the amide plane resulting in the carbonyl `pointing away' from the peroxo hole and thus establishing a hydrogen bond between one of the HD hydrogens and W1 [Fig. 4(g)]. Taken together, experiment and theory indicate that Asn254 is rather dynamic and able to explore a range of conformations.
3.6. Asn254–W1 coupling
Together with Thr57*, Asn254 is involved in W1 binding and also in the stabilization of O2 and 5PIU during the oxidative stage of catalysis (Bui et al., 2014). As both H2O and O2 must access the peroxo hole during the reaction cycle we considered the possiblity that Asn254 dynamics might offer a mechanism to modulate W1 stability at this location. To investigate this further, we performed QM/MM MD simulations. Similar to the approach employed in our classical MD simulations, we performed independent 8 ps-long runs using initial structures in which the Pro253–Lys255 region was either in the major or minor conformation. Asn254 torsion angles remain stable during the simulations indicating energetic minima (Fig. S6). We then inspected W1 movements using the AZAC5–W1O distance and the AZAC4–AZAC5–W1O–W1average (H1,H2) dihedral angle as coordinate parameters. For the major conformation [Fig. 6(a)], W1 mostly explores distances between 3.1 and 4.0 Å with a maximum at 3.6 Å and a narrow dihedral distribution centred around 10°. This is in excellent agreement with the experimental values in our structure (3.47 Å, 6°). The snapshots at different time points in Figs. 6(c), 6(d) and 6(e) provide further visual confirmation that W1 is relatively stable within the peroxo hole when Asn254 is in the major conformation (see Supplementary Movie S1 in the supporting information). The situation is different for the minor conformation [Fig. 6(b)]. Under these conditions W1 moves at longer distances from AZAC5 (maximum at 4.2–4.4 Å) exhibiting a broader distribution of different dihedral angles (225–255°). This is the result of W1 leaving the peroxo hole (we observe this after ∼500 fs) stabilized by the amide group of Asn254 as seen in the snapshots in Figs. 6(f) 6(g) and 6(h) (see Supplementary Movie S2). In our final crystallographic model W′ refines at a distance of 4.4 Å from AZAC5 and at hydrogen-bonding distance of 2.75 Å from Asn254OD1 (minor). Taken together, experiment and theory (both standard MD and QM/MM-MD) are in excellent agreement and support the view that the minor crystallographic conformation favours the displacement of W1 from the peroxo hole.
4. Discussion
The synergistic use of experimental and theoretical approaches provides a powerful way to gain insight into the mechanisms of enzyme catalysis (Bui & Steiner, 2016; Huggins et al., 2019). Here, we have integrated high-resolution room-temperature neutron/X-ray crystallographic analysis with biomolecular simulations to study co-factor-independent UOX in complex with its competitive inhibitor 8AZA in the presence of the catalytic water molecule W1 bound in the peroxo hole.
Our room-temperature analysis has uncovered structural heterogeneity for the Pro253–Lys255 tripeptide. This involves residue Asn254 which interacts directly with W1. High-resolution crystallographic studies at room temperature are extremely valuable in the detection of functional multiconformers or rare states in enzymes restricted under cryogenic conditions (Fraser et al., 2011). Early studies involving crystals of myoglobin (Frauenfelder et al., 1987, 1979) and ribonuclease (Rasmussen et al., 1992; Tilton et al., 1992) have shown that conformational distributions are restricted at low temperatures and, for human proline isomerase cyclophilin A, only collection of X-ray diffraction data at ambient temperature allowed an agreement between the crystallographic and NMR functional conformational substates to be obtained (Fraser et al., 2009). The two Asn254 conformers observed here map onto high-probability clusters of torsion-angle distributions obtained from the simulations [yellow and dark grey clusters in Fig. 5(b)]. Asn254 side-chain variability has also been observed under ligand-free conditions (Retailleau et al., 2004). In the absence of a ligand, UOX has been crystallized in a different (P21) with two complete tetramers in the (PDB entry 1r56; Retailleau et al., 2004). In the 2.6 Å resolution structure at 288 K, Asn254 samples a range of states that map onto the aqua–green clusters in Fig. 5(b). Conformational heterogeneity for Asn254 therefore appears to be a constitutive property of UOX. In cholesterol oxidase flavoenzyme, an arginine and a glutamate residue at the active site, found in two conformations, have been proposed to control oxygen access to the active-site cavity from a neighbouring channel (Coulombe et al., 2001). It is tempting to suggest that Asn254 dynamics might play a similar role considering that this residue marks the boundary between the active site and an adjacent hydrophobic cavity proposed to act as a transient O2 reservoir (Colloc'h & Prangé, 2014).
A mechanism for O2 gating facilitated by Asn254 dynamics probably involves W1 destabilization. MD simulations show that, though the main crystallographic conformation favours long W1 residence times in the peroxo hole, the minor conformation is associated with a higher probability of escaping from it. In fact, this conformation stabilizes W1 outside the peroxo hole at a position consistent with that of the low-occupancy W′ in our structure. It is therefore possible that the marginally lower than unitary occupancy for W1 (occupancy 0.93) arises from its fractional mobilization to the W′ site promoted by the minor conformation of Asn254. Coupling of Asn254 conformational heterogeneity and W1 stability is further supported by the of ligand-free UOX in which Asn254 is visualized adopting multiple alternative rotamers (Retailleau et al., 2004). Here, W1 is observed in only three of the eight independent active sites of the two UOX tetramers present in the A mechanism that takes advantage of protein dynamics to `actively' promote W1/O2 exchange for the oxidation step of catalysis appears particularly useful for cofactor-independent UOX as O2 affinity for the electrostatically positive active site is presumably rather low. Therefore, similarly to other cofactor-independent O2-metabolizing enzymes such as DpgC, poor ability to concentrate dioxygen in the active site, for example by metal coordination, must be compensated by mechanisms that increase the frequency of O2 entries (Di Russo et al., 2015). Asn254 dynamics is therefore suggested to decrease W1 residence times favouring O2 exchange with the neighbouring hydrophobic pocket.
The availability of the joint neutron/X-ray structures of DUOX–8AZA–W1 allows for an interesting comparison with that of the UOX–8AZA–Cl− complex (PDB entry 4n3m) solved in previous work (Oksanen et al., 2014). Functionally, the key difference between these structures resides in the nature of the molecule present in the peroxo hole (water versus chloride). Thus, comparative analysis is particularly meaningful for active-site residues. While both structures show that 8AZA is bound as a monoanion deprotonated at N3, differences exist for the important Lys10*–Thr57* dyad (Fig. S7). In UOX–8AZA–Cl−, Lys10* is modelled in the neutral NH2 state accepting a hydrogen bond from Thr57* whose DG1 atom points toward Lys10*NZ. The authors, however, state that a second conformation (not modelled) might exist in which Lys10* is positively charged (–NH3+) and the protonated Thr57* hydroxyl group points towards the Cl− ion. In UOX–8AZA–W1, the catalytic water W1 is present in the peroxo hole as a neutral H2O (here we will use H for both 1H and D) oriented at 45° with respect to the ligand plane and stabilized by a combination of hydrogen bonds. These involve not only the side chains of Thr57* and Asn254 (both conformations) but also an O—H⋯π interaction with the 8AZA monoanion. We find clear evidence for a charged Lys10*–NH3+ side chain whilst the Thr57*HG1 atom does not appear to be strongly localized. However, our experiment and simulation agree that the Thr57*OG1–HG1 bond points, on average, in the direction of the ligand, i.e. rotated by almost 180° compared with what was reported in the presence of chloride. We suggest that differences in active-site protonation states between these structures most likely reflect the nature and charge of the different molecules bound in the peroxo hole. In future work it will therefore be important to determine the neutron structure of the UOX–UA–W1 complex under anaerobic conditions as this has direct implications on the mechanism of substrate-induced O2 activation.
The protonation state observed in our UOX–8AZA–W1 complex is directly relevant for the hydration stage of the et al., 2017) proposed the existence of a strong hydrogen-bonding network between charged Lys10*, Thr57*, W1 and DHUN7 (note that in their paper the authors use protein numbering for B. subtilis UOX in which Lys9* and Thr69* correspond to A. flavus Lys10* and Thr57*, respectively). As W1O approaches DHUC5 in the nucleophilic attack, W1H1 comes close to Thr57*OG1 while HG1 gradually forms a bond with DHUN7 changing the C5=N7 double bond to a single bond (Fig. 7). In this mechanism the positively charged Lys10* stabilizes the negative charge on Thr57* during the transition state. Our structural data and simulation results are fully consistent with this proposal. Specifically, the charged Lys10* side chain favours the orientation of the Thr57*OG1–HG1 bond toward the π system of the ligand, thus facilitating proton transfer to N7. As W1O is oriented toward C5, nucleophilic attack is expected to occur on this atom (instead of C4) while Thr57* regains neutrality by abstracting the H1 proton from W1. Thr57HG1 is therefore central to a Lys10–NH3+-assisted proton-relay system.
A recent theoretical study (WeiSupporting information
Supporting movie 1. DOI: https://doi.org/10.1107/S2052252520013615/jt5052sup1.mpg
Supplorting movie 2. DOI: https://doi.org/10.1107/S2052252520013615/jt5052sup2.mpg
Supporting figures. DOI: https://doi.org/10.1107/S2052252520013615/jt5052sup3.pdf
Footnotes
‡These authors contributed equally to this work.
Acknowledgements
We thank the ESRF for access to the BM30A beamline and Dr Jean-Luc Ferrer for assisting LMcG during data collection. Dr Luca Signor (IBS, Grenoble) is acknowledged for his help with the
analysis.Funding information
LMcG is supported by an ILL-sponsored PhD studentship awarded to RAS and MPB. SB is funded by BBSRC grant BB/P000169/1 awarded to RAS and ER. The authors thank the ILL for access to the Deuteration Laboratory (proposal No. DL-03-189) and for the provision of beam time on the LADI-III diffractometer (proposal Nos. 8-01-490, 8-01-522, 8-01-549). This work used the platforms of the Grenoble Instruct-ERIC center (ISBG; UMS 3518 CNRS-CEA-UGA-EMBL) within the Grenoble Partnership for Structural Biology (PSB), supported by FRISBI (grant No. ANR-10-INBS-05-02) and GRAL, financed within the University Grenoble Alpes graduate school (Ecoles Universitaires de Recherche) CBH-EUR-GS (grant No. ANR-17-EURE-0003). ER also acknowledges EPSRC (EP/R013012/1) and ERC (project 757850 BioNet) grants. The authors acknowledge use of the research computing facility at King's College London, Rosalind.
References
Afonine, P. V., Moriarty, N. W., Mustyakimov, M., Sobolev, O. V., Terwilliger, T. C., Turk, D., Urzhumtsev, A. & Adams, P. D. (2015). Acta Cryst. D71, 646–666. Web of Science CrossRef IUCr Journals Google Scholar
Afonine, P. V., Mustyakimov, M., Grosse-Kunstleve, R. W., Moriarty, N. W., Langan, P. & Adams, P. D. (2010). Acta Cryst. D66, 1153–1163. Web of Science CrossRef CAS IUCr Journals Google Scholar
Andersen, H. C. (1983). J. Comput. Phys. 52, 24–34. CrossRef CAS Web of Science Google Scholar
Arzt, S., Campbell, J. W., Harding, M. M., Hao, Q. & Helliwell, J. R. (1999). J. Appl. Cryst. 32, 554–562. Web of Science CrossRef CAS IUCr Journals Google Scholar
Becke, A. D. (1993). J. Chem. Phys. 98, 5648–5652. CrossRef CAS Web of Science Google Scholar
Blakeley, M. P. & Podjarny, A. (2018). Emerg. Top. Life Sci. 2, 39–55. CAS Google Scholar
Blakeley, M. P., Teixeira, S. C. M., Petit-Haertlein, I., Hazemann, I., Mitschler, A., Haertlein, M., Howard, E. & Podjarny, A. D. (2010). Acta Cryst. D66, 1198–1205. Web of Science CrossRef CAS IUCr Journals Google Scholar
Breneman, C. M. & Wiberg, K. B. (1990). J. Comput. Chem. 11, 361–373. CrossRef CAS Web of Science Google Scholar
Brooks, B. R., Brooks, C. L., Mackerell, A. D. Jr, Nilsson, L., Petrella, R. J., Roux, B., Won, Y., Archontis, G., Bartels, C., Boresch, S., Caflisch, A., Caves, L., Cui, Q., Dinner, A. R., Feig, M., Fischer, S., Gao, J., Hodoscek, M., Im, W., Kuczera, K., Lazaridis, T., Ma, J., Ovchinnikov, V., Paci, E., Pastor, R. W., Post, C. B., Pu, J. Z., Schaefer, M., Tidor, B., Venable, R. M., Woodcock, H. L., Wu, X., Yang, W., York, D. M. & Karplus, M. (2009). J. Comput. Chem. 30, 1545–1614. Web of Science CrossRef PubMed CAS Google Scholar
Bui, S. & Steiner, R. A. (2016). Curr. Opin. Struct. Biol. 41, 109–118. Web of Science CrossRef CAS PubMed Google Scholar
Bui, S., von Stetten, D., Jambrina, P. G., Prangé, T., Colloc'h, N., de Sanctis, D., Royant, A., Rosta, E. & Steiner, R. A. (2014). Angew. Chem. Int. Ed. 53, 13710–13714. Web of Science CrossRef CAS Google Scholar
Campbell, J. W. (1995). J. Appl. Cryst. 28, 228–236. CrossRef CAS Web of Science IUCr Journals Google Scholar
Colloc'h, N., el Hajji, M., Bachet, B., L'Hermite, G., Schiltz, M., Prangé, T., Castro, B. & Mornon, J. P. (1997). Nat. Struct. Biol. 4, 947–952. CAS PubMed Web of Science Google Scholar
Colloc'h, N. & Prangé, T. (2014). FEBS Lett. 588, 1715–1719. Web of Science CAS PubMed Google Scholar
Colloc'h, N., Sopkova-de Oliveira Santos, J., Retailleau, P., Vivarès, D., Bonneté, F., Langlois d'Estainto, B., Gallois, B., Brisson, A., Risso, J. J., Lemaire, M., Prangé, T. & Abraini, J. H. (2007). Biophys. J. 92, 217–224. Web of Science PubMed CAS Google Scholar
Coulombe, R., Yue, K. Q., Ghisla, S. & Vrielink, A. (2001). J. Biol. Chem. 276, 30435–30441. Web of Science CrossRef PubMed CAS Google Scholar
Darden, T., York, D. & Pedersen, L. (1993). J. Chem. Phys. 98, 10089–10092. CrossRef CAS Web of Science Google Scholar
DeLano, W. L. (2002). CCP4 Newsl. Protein Crystallogr. 40, 82–92. Google Scholar
Di Russo, N. V., Condurso, H. L., Li, K., Bruner, S. D. & Roitberg, A. E. (2015). Chem. Sci. 6, 6341–6348. Web of Science CrossRef CAS PubMed 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. (2006). Acta Cryst. D62, 72–82. Web of Science CrossRef CAS IUCr Journals Google Scholar
Fetzner, S. & Steiner, R. A. (2010). Appl. Microbiol. Biotechnol. 86, 791–804. Web of Science CrossRef CAS PubMed Google Scholar
Fisher, S. J., Blakeley, M. P., Howard, E. I., Petit-Haertlein, I., Haertlein, M., Mitschler, A., Cousido-Siah, A., Salvay, A. G., Popov, A., Muller-Dieckmann, C., Petrova, T. & Podjarny, A. (2014). Acta Cryst. D70, 3266–3272. Web of Science CrossRef IUCr Journals Google Scholar
Fraser, J. S., Clarkson, M. W., Degnan, S. C., Erion, R., Kern, D. & Alber, T. (2009). Nature, 462, 669–673. Web of Science CrossRef PubMed CAS Google Scholar
Fraser, J. S., van den Bedem, H., Samelson, A. J., Lang, P. T., Holton, J. M., Echols, N. & Alber, T. (2011). Proc. Natl Acad. Sci. USA, 108, 16247–16252. Web of Science CrossRef CAS PubMed Google Scholar
Frauenfelder, H., Hartmann, H., Karplus, M., Kuntz, I. D. Jr, Kuriyan, J., Parak, F., Petsko, G. A., Ringe, D., Tilton, R. F. Jr & Connolly, M. L. (1987). Biochemistry, 26, 254–261. CrossRef CAS PubMed Web of Science Google Scholar
Frauenfelder, H., Petsko, G. A. & Tsernoglou, D. (1979). Nature, 280, 558–563. CrossRef CAS PubMed Web of Science Google Scholar
French, S. & Wilson, K. (1978). Acta Cryst. A34, 517–525. CrossRef CAS IUCr Journals Web of Science Google Scholar
Frisch, M. J., Trucks, G. W., Schlegel, H. B., Scuseria, G. M., Robb, M. A., Cheeseman, J. R., Scalmani, G., Barone, V., Mennucci, B., Petersson, G. A., Nakatsuji, H., Caricato, M., Li, X., Hratchian, H. P., Izmaylov, A. F., Bloino, J., Zheng, G., Sonnenberg, J. L., Hada, M., Ehara, M., Toyota, K., Fukuda, R., Hasegawa, J., Ishida, M., Nakajima, T., Honda, Y., Kitao, O., Nakai, H., Vreven, T., Montgomery, J. A. Jr, Peralta, J. E., Ogliaro, F., Bearpark, M., Heyd, J. J., Brothers, E., Kudin, K. N., Staroverov, V. N., Keith, T., Kobayashi, R., Normand, J., Raghavachari, K., Rendell, A., Burant, J. C., Iyengar, S. S., Tomasi, J., Cossi, M., Rega, N., Millam, J. M., Klene, M., Knox, J. E., Cross, J. B., Bakken, V., Adamo, C., Jaramillo, J., Gomperts, R., Stratmann, R. E., Yazyev, O., Austin, A. J., Cammi, R., Pomelli, C., Ochterski, J. W., Martin, R. L., Morokuma, K., Zakrzewski, V. G., Voth, G. A., Salvador, P., Dannenberg, J. J., Dapprich, S., Daniels, A. D., Farkas, O., Foresman, J. B., Ortiz, J. V., Cioslowski, J. & Fox, D. J. (2013). Gaussian09. Gaussian, Inc., Wallingford CT, USA. Google Scholar
Gabison, L., Chiadmi, M., El Hajji, M., Castro, B., Colloc'h, N. & Prangé, T. (2010). Acta Cryst. D66, 714–724. Web of Science CrossRef CAS IUCr Journals Google Scholar
Gabison, L., Prangé, T., Colloc'h, N., El Hajji, M., Castro, B. & Chiadmi, M. (2008). BMC Struct. Biol. 8, 32. Google Scholar
Garman, E. F. (2010). Acta Cryst. D66, 339–351. Web of Science CrossRef CAS IUCr Journals Google Scholar
Grimme, S., Ehrlich, S. & Goerigk, L. (2011). J. Comput. Chem. 32, 1456–1465. Web of Science CrossRef CAS PubMed Google Scholar
Haertlein, M., Moulin, M., Devos, J. M., Laux, V., Dunne, O. & Forsyth, V. T. (2016). Methods Enzymol. 566, 113–157. Web of Science CrossRef CAS PubMed Google Scholar
Hibi, T., Kume, A., Kawamura, A., Itoh, T., Fukada, H. & Nishiya, Y. (2016). Biochemistry, 55, 724–732. Web of Science CrossRef CAS PubMed Google Scholar
Huggins, D. J., Biggin, P. C., Dämgen, M. A., Essex, J. W., Harris, S. A., Henchman, R. H., Khalid, S., Kuzmanic, A., Laughton, C. A., Michel, J., Mulholland, A. J., Rosta, E., Sansom, M. S. P. & van der Kamp, M. W. (2019). WIREs Comput. Mol. Sci. 9, e1393. Google Scholar
Imhoff, R. D., Power, N. P., Borrok, M. J. & Tipton, P. A. (2003). Biochemistry, 42, 4094–4100. Web of Science CrossRef PubMed CAS Google Scholar
Jo, S., Kim, T., Iyer, V. G. & Im, W. (2008). J. Comput. Chem. 29, 1859–1865. Web of Science CrossRef PubMed CAS Google Scholar
Kahn, K. (1999). Bioorg. Chem. 27, 351–362. Web of Science CrossRef CAS Google Scholar
Kahn, K., Serfozo, P. & Tipton, P. A. (1997). J. Am. Chem. Soc. 119, 5435–5442. CrossRef CAS Web of Science Google Scholar
Kahn, K. & Tipton, P. A. (1997). Biochemistry, 36, 4731–4738. CrossRef CAS PubMed Web of Science Google Scholar
Kahn, K. & Tipton, P. A. (1998). Biochemistry, 37, 11651–11659. Web of Science CrossRef CAS PubMed Google Scholar
Koruza, K., Lafumat, B., Nyblom, M., Mahon, B. P., Knecht, W., McKenna, R. & Fisher, S. Z. (2019). Acta Cryst. D75, 895–903. Web of Science CrossRef IUCr Journals Google Scholar
Krishnan, R., Binkley, J. S., Seeger, R. & Pople, J. A. (1980). J. Chem. Phys. 72, 650–654. CrossRef CAS Web of Science Google Scholar
Liebschner, D., Afonine, P. V., Baker, M. L., Bunkóczi, G., Chen, V. B., Croll, T. I., Hintze, B., Hung, L.-W., Jain, S., McCoy, A. J., Moriarty, N. W., Oeffner, R. D., Poon, B. K., Prisant, M. G., Read, R. J., Richardson, J. S., Richardson, D. C., Sammito, M. D., Sobolev, O. V., Stockwell, D. H., Terwilliger, T. C., Urzhumtsev, A. G., Videau, L. L., Williams, C. J. & Adams, P. D. (2019). Acta Cryst. D75, 861–877. Web of Science CrossRef IUCr Journals Google Scholar
Logan, D. T. (2020). Methods Enzymol. 634, 201–224. Web of Science CrossRef PubMed Google Scholar
Marchetti, M., Liuzzi, A., Fermi, B., Corsini, R., Folli, C., Speranzini, V., Gandolfi, F., Bettati, S., Ronda, L., Cendron, L., Berni, R., Zanotti, G. & Percudani, R. (2016). Sci. Rep. 6, 38302. Web of Science CrossRef PubMed Google Scholar
McLean, A. D. & Chandler, G. S. (1980). J. Chem. Phys. 72, 5639–5648. CrossRef CAS Web of Science Google Scholar
Oda, M., Satta, Y., Takenaka, O. & Takahata, N. (2002). Mol. Biol. Evol. 19, 640–653. Web of Science CrossRef PubMed CAS Google Scholar
Oksanen, E., Blakeley, M. P., El-Hajji, M., Ryde, U. & Budayova-Spano, M. (2014). PLoS One, 9, e86651. Web of Science CrossRef PubMed Google Scholar
Pfrimer, P., de Moraes, L. M., Galdino, A. S., Salles, L. P., Reis, V. C., De Marco, J. L., Prates, M. V., Bloch, C. Jr & Torres, F. A. (2010). J. Biomed. Biotechnol. 2010, 674908. Web of Science CrossRef PubMed Google Scholar
Phillips, J. C., Braun, R., Wang, W., Gumbart, J., Tajkhorshid, E., Villa, E., Chipot, C., Skeel, R. D., Kalé, L. & Schulten, K. (2005). J. Comput. Chem. 26, 1781–1802. Web of Science CrossRef PubMed CAS Google Scholar
Phillips, J. C., Hardy, D. J., Maia, J. D. C., Stone, J. E., Ribeiro, J. V., Bernardi, R. C., Buch, R., Fiorin, G., Hénin, J., Jiang, W., McGreevy, R., Melo, M. C. R., Radak, B. K., Skeel, R. D., Singharoy, A., Wang, Y., Roux, B., Aksimentiev, A., Luthey-Schulten, Z., Kalé, L. V., Schulten, K., Chipot, C. & Tajkhorshid, E. (2020). J. Chem. Phys. 153, 044130. Web of Science CrossRef PubMed Google Scholar
Ramazzina, I., Folli, C., Secchi, A., Berni, R. & Percudani, R. (2006). Nat. Chem. Biol. 2, 144–148. Web of Science CrossRef PubMed CAS Google Scholar
Rasmussen, B. F., Stock, A. M., Ringe, D. & Petsko, G. A. (1992). Nature, 357, 423–424. CrossRef PubMed CAS Web of Science Google Scholar
Ravelli, R. B. & McSweeney, S. M. (2000). Structure, 8, 315–328. Web of Science CrossRef PubMed CAS Google Scholar
Retailleau, P., Colloc'h, N., Vivarès, D., Bonneté, F., Castro, B., El Hajji, M., Mornon, J.-P., Monard, G. & Prangé, T. (2004). Acta Cryst. D60, 453–462. Web of Science CrossRef CAS IUCr Journals Google Scholar
Shao, Y., Gan, Z., Epifanovsky, E., Gilbert, A. T. B., Wormit, M., Kussmann, J., Lange, A. W., Behn, A., Deng, J., Feng, X., Ghosh, D., Goldey, M., Horn, P. R., Jacobson, L. D., Kaliman, I., Khaliullin, R. Z., Kuś, T., Landau, A., Liu, J., Proynov, E. I., Rhee, Y. M., Richard, R. M., Rohrdanz, M. A., Steele, R. P., Sundstrom, E. J., Woodcock, H. L. III, Zimmerman, P. M., Zuev, D., Albrecht, B., Alguire, E., Austin, B., Beran, G. J. O., Bernard, Y. A., Berquist, E., Brandhorst, K., Bravaya, K. B., Brown, S. T., Casanova, D., Chang, C.., Chen, Y., Chien, S. H., Closser, K. D., Crittenden, D. L., Diedenhofen, M., DiStasio, R. A. Jr., Do, H., Dutoi, A. D., Edgar, R. G., Fatehi, S., Fusti-Molnar, L., Ghysels, A., Golubeva-Zadorozhnaya, A., Gomes, J., Hanson-Heine, M. W. D., Harbach, P. H. P., Hauser, A. W., Hohenstein, E. G., Holden, Z. C., Jagau, T.., Ji, H., Kaduk, B., Khistyaev, K., Kim, J., Kim, J., King, R. A., Klunzinger, P., Kosenkov, D., Kowalczyk, T., Krauter, C. M., Lao, K. U., Laurent, A. D., Lawler, K. V., Levchenko, S. V., Lin, C. Y., Liu, F., Livshits, E., Lochan, R. C., Luenser, A., Manohar, P., Manzer, S. F., Mao, S.., Mardirossian, N., Marenich, A. V., Maurer, S. A., Mayhall, N. J., Neuscamman, E., Oana, C. M., Olivares-Amaya, R., O'Neill, D. P., Parkhill, J. A., Perrine, T. M., Peverati, R., Prociuk, A., Rehn, D. R., Rosta, E., Russ, N. J., Sharada, S. M., Sharma, S., Small, D. W., Sodt, A., Stein, T., Stück, D., Su, Y.., Thom, A. J. W., Tsuchimochi, T., Vanovschi, V., Vogt, L., Vydrov, O., Wang, T., Watson, M. A., Wenzel, J., White, A., Williams, C. F., Yang, J., Yeganeh, S., Yost, S. R., You, Z.., Zhang, I. Y., Zhang, X., Zhao, Y., Brooks, B. R., Chan, G. K. L., Chipman, D. M., Cramer, C. J., Goddard, W. A. III, Gordon, M. S., Hehre, W. J., Klamt, A., Schaefer, H. F. III, Schmidt, M. W., Sherrill, C. D., Truhlar, D. G., Warshel, A., Xu, X., Aspuru-Guzik, A., Baer, R., Bell, A. T., Besley, N. A., Chai, J.., Dreuw, A., Dunietz, B. D., Furlani, T. R., Gwaltney, S. R., Hsu, C.., Jung, Y., Kong, J., Lambrecht, D. S., Liang, W., Ochsenfeld, C., Rassolov, V. A., Slipchenko, L. V., Subotnik, J. E., Van Voorhis, T., Herbert, J. M., Krylov, A. I., Gill, P. M. W. & Head-Gordon, M. (2015). Mol. Phys. 113, 184–215. Web of Science CrossRef CAS Google Scholar
Shu, F., Ramakrishnan, V. & Schoenborn, B. P. (2000). Proc. Natl Acad. Sci. USA, 97, 3872–3877. Web of Science CrossRef PubMed CAS Google Scholar
Steiner, T. & Koellner, G. (2001). J. Mol. Biol. 305, 535–557. Web of Science CrossRef PubMed CAS Google Scholar
Terkeltaub, R. (2010). Nat. Rev. Rheumatol. 6, 30–38. Web of Science CrossRef CAS PubMed Google Scholar
Tilton, R. F. Jr, Dewan, J. C. & Petsko, G. A. (1992). Biochemistry, 31, 2469–2481. CrossRef PubMed CAS Web of Science Google Scholar
Wei, D., Huang, X., Qiao, Y., Rao, J., Wang, L., Liao, F. & Zhan, C. G. (2017). ACS Catal. 7, 4623–4636. Web of Science CrossRef CAS PubMed Google Scholar
Weik, M., Ravelli, R. B., Kryger, G., McSweeney, S., Raves, M. L., Harel, M., Gros, P., Silman, I., Kroon, J. & Sussman, J. L. (2000). Proc. Natl Acad. Sci. USA, 97, 623–628. Web of Science CrossRef PubMed CAS 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
Winter, G. (2010). J. Appl. Cryst. 43, 186–190. Web of Science CrossRef CAS IUCr Journals 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.