Insights into the Proton Transfer Mechanism of a Bilin Reductase PcyA Following Neutron Crystallography

A study on the catalytic activity of polypeptides toward the hydrolysis of glucoside compounds gastrodin, polydatin and esculin

The self-assembly of a series of catalytically active polypeptides toward hydrolysis of glucoside compounds, namely, gastrodin, polydatin and esculin was investigated. These active peptides are composed of two functional fragments: one is the hydrophobic sequence LHLHLRL, which forms assembling segments in the presence of Zn ions (Zn²⁺); another functional sequence of active peptides are catalytic sites such as Glu (E), Asp (D) and His (H), where carboxylic acids (-COOH) or imidazole groups act like scissors to cleave glucoside bonds of the compounds (according to the acid-base coupling mechanism). The effects of the amino acid sequence of the peptide, Zn²⁺ concentration, pH and the size or steric hindrance of glucoside compounds on the hydrolytic activity were studied. It was found that the crystalline structure of assembled peptides was crucial to provide the peptide with catalytic hydrolytic activity. Noncovalent interaction index was used to analyse the noncovalent interaction of PEs with glucoside compounds, including hydrogen bonds, van der Waals, and steric effect in the complexes. The binding energy of complexes, the direction and site of nucleophilic attack during deglycosylation processes were also investigated by molecular docking and the electron density Laplace function. This revealed that the differences in the hydrolytic activity of peptides toward glucoside compounds with different sizes originated from different hydrogen bond interactions between the peptides and substrates. These active peptides may find application in the preparation of drugs by de-glycosylation of natural compounds.

Neutron crystallography and quantum chemical analysis of bilin reductase PcyA mutants reveal substrate and catalytic residue protonation states

PcyA, a ferredoxin-dependent bilin pigment reductase, catalyzes the site-specific reduction of the two vinyl groups of biliverdin (BV), producing phycocyanobilin. Previous neutron crystallography detected both the neutral BV and its protonated form (BVH⁺) in the wildtype (WT) PcyA-BV complex, and a nearby catalytic residue Asp105 was found to have two conformations (protonated and deprotonated). Semiempirical calculations have suggested that the protonation states of BV are reflected in the absorption spectrum of the WT PcyA-BV complex. In the previously determined absorption spectra of the PcyA D105N and I86D mutants, complexed with BV, a peak at 730 nm, observed in the WT, disappeared and increased, respectively. Here, we performed neutron crystallography and quantum chemical analysis of the D105N-BV and I86D-BV complexes to determine the protonation states of BV and the surrounding residues and study the correlation between the absorption spectra and protonation states around BV. Neutron structures elucidated that BV in the D105N mutant is in a neutral state, whereas that in the I86D mutant is dominantly in a protonated state. Glu76 and His88 showed different hydrogen bonding with surrounding residues compared with WT PcyA, further explaining why D105N and I86D have much lower activities for phycocyanobilin synthesis than the WT PcyA. Our quantum mechanics/molecular mechanics calculations of the absorption spectra showed that the spectral change in D105N arises from Glu76 deprotonation, consistent with the neutron structure. Collectively, our findings reveal more mechanistic details of bilin pigment biosynthesis.

Re-evaluation of protein neutron crystallography with and without X-ray/neutron joint refinement

Protein neutron crystallography is a powerful technique to determine the positions of H atoms, providing crucial biochemical information such as the protonation states of catalytic groups and the geometry of hydrogen bonds. Recently, the crystal structure of a bacterial copper amine oxidase was determined by joint refinement using X-ray and neutron diffraction data sets at resolutions of 1.14 and 1.72 Å, respectively [Murakawa et al. (2020 ▸). Proc. Natl Acad. Sci. USA, 117, 10818-10824]. While joint refinement is effective for the determination of the accurate positions of heavy atoms on the basis of the electron density, the structural information on light atoms (hydrogen and deuterium) derived from the neutron diffraction data might be affected by the X-ray data. To unravel the information included in the neutron diffraction data, the structure determination was conducted again using only the neutron diffraction data at 1.72 Å resolution and the results were compared with those obtained in the previous study. Most H and D atoms were identified at essentially the same positions in both the neutron-only and the X-ray/neutron joint refinements. Nevertheless, neutron-only refinement was found to be less effective than joint refinement in providing very accurate heavy-atom coordinates that lead to significant improvement of the neutron scattering length density map, especially for the active-site cofactor. Consequently, it was confirmed that X-ray/neutron joint refinement is crucial for determination of the real chemical structure of the catalytic site of the enzyme.

Identification of significant residues for intermediate accumulation in phycocyanobilin synthesis

Phycocyanobilin, the primary pigment of both light perception and light-harvesting in cyanobacteria, is synthesized from biliverdin IXα (BV) through intermediate 18¹, 18²-dihydrobiliverdin (18¹, 18²-DHBV) by a phycocyanobilin:ferredoxin oxidoreductase (PcyA). In our previous study, we discovered two PcyA homologs (AmPcyAc and AmPcyAp) derived from Acaryochloris marina MBIC 11017 (A. marina) that exceptionally uses chlorophyll d as the primary photosynthetic pigment, absorbing longer wavelength far-red light than chlorophyll a, the photosynthetic pigment found in most cyanobacteria. Biochemical characterization of the two PcyA homologs identified functional diversification of these two enzymes: AmPcyAc provides 18¹, 18²-DHBV, and PCB to the cyanobacteriochrome (CBCR) photoreceptors, whereas, AmPcyAp specifically provides PCB to the light-harvesting phycobilisome subunit. In this study, we focused on the residues necessary for 18¹, 18²-DHBV supply to the CBCR photoreceptors by AmPcyAc. Based on the SyPcyA structure, we concentrated on the 30 residues that constitute the substrate-binding pocket. Among them, we discovered that Leu151 and Val225 in AmPcyAc were both substituted with isoleucine. During the enzymatic reaction, the SyPcyA variant molecule, possessing V225I and L151I replacements, accumulates the 18¹, 18²-DHBV and supplies it to a CBCR molecule derived from A. marina. It is worth noting that the substitution of Val225 with isoleucine was specifically conserved among the Acaryochloris genus. Collectively, we propose that the specific evolution of PcyA among the Acaryochloris genus may correlate with the acquisition of Chl. d synthetic ability and growth in long-wavelength far-red light environments.

Chlorophyllide a oxidoreductase Preferentially Catalyzes 8-Vinyl Reduction over B-Ring Reduction of 8-Vinyl Chlorophyllide a in the Late Steps of Bacteriochlorophyll Biosynthesis

Bacteriochlorophyll a (BChl) is an essential pigment for anoxygenic photosynthesis. In late steps of the BChl biosynthesis of Rhodobacter capsulatus, the C8 vinyl group and C7=C8 double bond of 8-vinyl chlorophyllide a (8 V-Chlide) are reduced by a C8 vinyl reductase (8VR), BciA, and a nitrogenase-like enzyme, chlorophyllide a oxidoreductase (COR), respectively, to produce 3-vinyl-bacteriochlorphyllide a. Recently, we discovered 8VR activity in COR. However, the kinetic parameters of the COR 8VR activity remain unknown, while those of the COR C7=C8 reductase activity and BciA have been reported. Here, we determined the kinetic parameters of COR 8VR activity by using 8 V-Chlide. The K_m value for 8 V-Chlide was 1.4 μM, which is much lower than the 6.2 μM determined for the C7=C8 reduction of Chlide. The kinetic parameters of the dual activities of COR suggest that COR catalyzes the reduction of the C8 vinyl group of 8 V-Chlide preferentially over C7=C8 reduction when both substrates are supplied during BChl biosynthesis.

Current status and near future plan of neutron protein crystallography at J-PARC

The IBARAKI Biological Crystal Diffractometer (iBIX) has been available for use at MLF (Material and Life Science Facility) in J-PARC (Japan Proton Accelerator Research Complex) since 2008. The development in state-of-the-art detector systems could enable iBIX to become one of the highest-performance neutron single-crystal diffractometers in the world. Here, together with other various developments, such as data reduction software, crystal growth, and new techniques in measurement coupled analysis, we provided new hydrogen and water structural data of several proteins and macromolecules. Although the proton power at MLF has not yet reached its planned maximum (1MW), a more powerful neutron source will be soon needed for neutron protein crystallography. A future idea is also proposed and discussed in this article.

Functional diversification of two bilin reductases for light perception and harvesting in unique cyanobacterium Acaryochloris marina MBIC 11017

Bilin pigments play important roles for both light perception and harvesting in cyanobacteria by binding to cyanobacteriochromes (CBCRs) and phycobilisomes (PBS), respectively. Among various cyanobacteria, Acaryochloris marina MBIC 11017 (A. marina 11017) exceptionally uses chlorophyll d as the main photosynthetic pigment absorbing longer wavelength light than the canonical pigment, chlorophyll a, indicating existence of a system to sense longer wavelength light than others. On the other hand, A. marina 11017 has the PBS apparatus to harvest short-wavelength orange light, similar to most cyanobacteria. Thus, A. marina 11017 might sense longer wavelength light and harvest shorter wavelength light by using bilin pigments. Phycocyanobilin (PCB) is the main bilin pigment of both systems. Phycocyanobilin:ferredoxin oxidoreductase (PcyA) catalyzes PCB synthesis from biliverdin via the intermediate 18¹ ,18² -dihydrobiliverdin (18¹ ,18² -DHBV), resulting in the stepwise shortening of the absorbing wavelengths. In this study, we found that A. marina 11017 exceptionally encodes two PcyA homologs, AmPcyAc and AmPcyAp. AmPcyAc is encoded on the main chromosome with most photoreceptor genes, whereas AmPcyAp is encoded on a plasmid with PBS-related genes. High accumulation of 18¹ ,18² -DHBV for extended periods was observed during the reaction catalyzed by AmPcyAc, whereas 18¹ ,18² -DHBV was transiently accumulated for a short period during the reaction catalyzed by AmPcyAp. CBCRs could sense longer wavelength far-red light through 18¹ ,18² -DHBV incorporation, whereas PBS could only harvest orange light through PCB incorporation, suggesting functional diversification of PcyA as AmPcyAc and AmPcyAp to provide 18¹ ,18² -DHBV and PCB to the light perception and harvesting systems, respectively.

Crystal structure of phytochromobilin synthase in complex with biliverdin IXα, a key enzyme in the biosynthesis of phytochrome

Phytochromobilin (PΦB) is a red/far-red light sensory pigment in plant phytochrome. PΦB synthase is a ferredoxin-dependent bilin reductase (FDBR) that catalyzes the site-specific reduction of bilins, which are sensory and photosynthesis pigments, and produces PΦB from biliverdin, a heme-derived linear tetrapyrrole pigment. Here, we determined the crystal structure of tomato PΦB synthase in complex with biliverdin at 1.95 Å resolution. The overall structure of tomato PΦB synthase was similar to those of other FDBRs, except for the addition of a long C-terminal loop and short helices. The structure further revealed that the C-terminal loop is part of the biliverdin-binding pocket and that two basic residues in the C-terminal loop form salt bridges with the propionate groups of biliverdin. This suggested that the C-terminal loop is involved in the interaction with ferredoxin and biliverdin. The configuration of biliverdin bound to tomato PΦB synthase differed from that of biliverdin bound to other FDBRs, and its orientation in PΦB synthase was inverted relative to its orientation in the other FDBRs. Structural and enzymatic analyses disclosed that two aspartic acid residues, Asp-123 and Asp-263, form hydrogen bonds with water molecules and are essential for the site-specific A-ring reduction of biliverdin. On the basis of these observations and enzymatic assays with a V121A PΦB synthase variant, we propose the following mechanistic product release mechanism: PΦB synthase-catalyzed stereospecific reduction produces 2(R)-PΦB, which when bound to PΦB synthase collides with the side chain of Val-121, releasing 2(R)-PΦB from the synthase.

Direct Observation of the Protonation States in the Mutant Green Fluorescent Protein

Neutron crystallography has been used to elucidate the protonation states for the enhanced green fluorescent protein, which has revolutionized imaging technologies. The structure has a deprotonated hydroxyl group in the fluorescent chromophore. Also, the protonation states of His148 and Thr203, as well as the orientation of a critical water molecule in direct contact with the chromophore, could be determined. The results demonstrate that the deprotonated hydroxyl group in the chromophore and the nitrogen atom ND1 in His148 are charged negatively and positively, respectively, forming an ion pair. The position of the two deuterium atoms in the critical water molecule appears to be displaced slightly toward the acceptor oxygen atoms according to their omit maps. This displacement implies the formation of an intriguing electrostatic potential realized inside of the protein. Our findings provide new insights into future protein design strategies along with developments in quantum chemical calculations.

Crystal structure of the first eukaryotic bilin reductase GtPEBB reveals a flipped binding mode of dihydrobiliverdin

Phycobilins are light-harvesting pigments of cyanobacteria, red algae, and cryptophytes. The biosynthesis of phycoerythrobilin (PEB) is catalyzed by the subsequent action of two ferredoxin-dependent bilin reductases (FDBRs). Although 15,16-dihydrobiliverdin (DHBV):ferredoxin oxidoreductase (PebA) catalyzes the two-electron reduction of biliverdin IXα to 15,16-DHBV, PEB:ferredoxin oxidoreductase (PebB) reduces this intermediate further to PEB. Interestingly, marine viruses encode the FDBR PebS combining both activities within one enzyme. Although PebA and PebS share a canonical fold with similar substrate-binding pockets, the structural determinants for the stereo- and regiospecific modification of their tetrapyrrole substrates are incompletely understood, also because of the lack of a PebB structure. Here, we solved the X-ray crystal structures of both substrate-free and -bound PEBB from the cryptophyte Guillardia theta at 1.90 and 1.65 Å, respectively. The structures of PEBB exhibit the typical α/β/α-sandwich fold. Interestingly, the open-chain tetrapyrrole substrate DHBV is bound in an unexpected flipped orientation within the canonical FDBR active site. Biochemical analyses of the WT enzyme and active site variants identified two central aspartate residues Asp-99 and Asp-219 as essential for catalytic activity. In addition, the conserved Arg-215 plays a critical role in substrate specificity, binding orientation, and active site integrity. Because these critical residues are conserved within certain FDBRs displaying A-ring reduction activity, we propose that they present a conserved mechanism for this reaction. The flipped substrate-binding mode indicates that two-electron reducing FDBRs utilize the same primary site within the binding pocket and that substrate orientation is the determinant for A- or D-ring regiospecificity.

Bilin-metabolizing enzymes: site-specific reductions catalyzed by two different type of enzymes

In mammals, the green heme metabolite biliverdin is converted to a yellow anti-oxidant by NAD(P)H-dependent biliverdin reductase (BVR), whereas in O₂-dependent photosynthetic organisms it is converted to photosynthetic or light-sensing pigments by ferredoxin-dependent bilin reductases (FDBRs). In NADP⁺-bound and biliverdin-bound BVR-A, two biliverdins are stacked at the binding cleft; one is positioned to accept hydride from NADPH, and the other appears to donate a proton to the first biliverdin through a neighboring arginine residue. During the FDBR-catalyzed reaction, electrons and protons are supplied to bilins from ferredoxin and from FDBRs and waters bound within FDBRs, respectively. Thus, the protonation sites of bilin and catalytic residues are important for the analysis of site-specific reduction. The neutron structure of FDBR sheds light on this issue.

Neutron scattering in the biological sciences: progress and prospects

The scattering of neutrons can be used to provide information on the structure and dynamics of biological systems on multiple length and time scales. Pursuant to a National Science Foundation-funded workshop in February 2018, recent developments in this field are reviewed here, as well as future prospects that can be expected given recent advances in sources, instrumentation and computational power and methods. Crystallography, solution scattering, dynamics, membranes, labeling and imaging are examined. For the extraction of maximum information, the incorporation of judicious specific deuterium labeling, the integration of several types of experiment, and interpretation using high-performance computer simulation models are often found to be particularly powerful.

Status of the neutron time-of-flight single-crystal diffraction data-processing software STARGazer

The STARGazer data-processing software is used for neutron time-of-flight (TOF) single-crystal diffraction data collected using the IBARAKI Biological Crystal Diffractometer (iBIX) at the Japan Proton Accelerator Research Complex (J-PARC). This software creates hkl intensity data from three-dimensional (x, y, TOF) diffraction data. STARGazer is composed of a data-processing component and a data-visualization component. The former is used to calculate the hkl intensity data. The latter displays the three-dimensional diffraction data with searched or predicted peak positions and is used to determine and confirm integration regions. STARGazer has been developed to make it easier to use and to obtain more accurate intensity data. For example, a profile-fitting method for peak integration was developed and the data statistics were improved. STARGazer and its manual, containing installation and data-processing components, have been prepared and provided to iBIX users. This article describes the status of the STARGazer data-processing software and its data-processing algorithms.

Hydronium Ions Accompanying Buried Acidic Residues Lead to High Apparent Dielectric Constants in the Interior of Proteins

Internal ionizable groups are known to play important roles in protein functions. A mystery that has attracted decades of extensive experimental and theoretical studies is the apparent dielectric constants experienced by buried ionizable groups, which are much higher than values expected for protein interiors. Many interpretations have been proposed, such as water penetration, conformational relaxation, local unfolding, protein intrinsic backbone fluctuations, etc. However, these interpretations conflict with many experimental observations. The virtual mixture of multiple states (VMMS) simulation method developed in our lab provides a direct approach for studying the equilibrium of multiple chemical states and can monitor p K_a values along simulation trajectories. Through VMMS simulations of staphylococcal nuclease (SNase) variants with internal Asp or Glu residues, we discovered that cations were attracted to buried deprotonated acidic groups and the presence of the nearby cations were essential to reproduce experimentally measured p K_a values. This finding, combined with structural analysis and validation simulations, suggests that the proton released from a deprotonation process stays near the deprotonated group inside proteins, possibly in the form of a hydronium ion. The existence of a proton near a buried charge has many implications in our understanding of protein functions.

QM/MM Investigation for Protonation States in a Bilin Reductase PcyA-Biliverdin IXα Complex

Herein we report quantum mechanical/molecular mechanical (QM/MM) studies to investigate the most probable protonation states of active site amino acids and bound substrate based on a recently reported neutron diffraction structure of phycocyanobilin:ferredoxin oxidoreductase (PcyA) by Unno et al. This structure was considered to be bound in its initial state of biliverdin IXα (BV), which has the C-pyrrole ring in the deprotonated state. The protonation state of BV suggested by neutron and spectroscopic studies is a stable, two-electron reduced complex with a bound hydronium ion. Several ambiguities in the neutron structure were observed which prompted a further theoretical analysis of the structure. This structural investigation provides new understanding of the PcyA and BV protonation states not previously reported in the literature. Our calculations suggest that the hydronium ion (H3 O+ ) is energetically unfavorable, preferentially protonating the neighboring His88 residue and that the C-ring of BV is not protonated.

Neutron macromolecular crystallography

Neutron diffraction techniques permit direct determination of the hydrogen (H) and deuterium (D) positions in crystal structures of biological macromolecules at resolutions of ∼1.5 and 2.5 Å, respectively. In addition, neutron diffraction data can be collected from a single crystal at room temperature without radiation damage issues. By locating the positions of H/D-atoms, protonation states and water molecule orientations can be determined, leading to a more complete understanding of many biological processes and drug-binding. In the last ca. 5 years, new beamlines have come online at reactor neutron sources, such as BIODIFF at Heinz Maier-Leibnitz Zentrum and IMAGINE at Oak Ridge National Laboratory (ORNL), and at spallation neutron sources, such as MaNDi at ORNL and iBIX at the Japan Proton Accelerator Research Complex. In addition, significant improvements have been made to existing beamlines, such as LADI-III at the Institut Laue-Langevin. The new and improved instrumentations are allowing sub-mm3 crystals to be regularly used for data collection and permitting the study of larger systems (unit-cell edges >100 Å). Owing to this increase in capacity and capability, many more studies have been performed and for a wider range of macromolecules, including enzymes, signalling proteins, transport proteins, sugar-binding proteins, fluorescent proteins, hormones and oligonucleotides; of the 126 structures deposited in the Protein Data Bank, more than half have been released since 2013 (65/126, 52%). Although the overall number is still relatively small, there are a growing number of examples for which neutron macromolecular crystallography has provided the answers to questions that otherwise remained elusive.

Evolution and molecular mechanism of four-electron reducing ferredoxin-dependent bilin reductases from oceanic phages

Ferredoxin-dependent bilin reductases (FDBRs) are a class of enzymes reducing the heme metabolite biliverdin IXα (BV) to form open-chain tetrapyrroles used for light-perception and light-harvesting in photosynthetic organisms. Thus far, seven FDBR families have been identified, each catalysing a distinct reaction and either transferring two or four electrons from ferredoxin onto the substrate. The newest addition to the family is PcyX, originally identified from metagenomics data derived from phage. Phylogenetically, PcyA is the closest relative catalysing the reduction of BV to phycocyanobilin. PcyX, however, converts the same substrate to phycoerythrobilin, resembling the reaction catalysed by cyanophage PebS. Within this study, we aimed at understanding the evolution of catalytic activities within FDBRs using PcyX as an example. Additional members of the PcyX clade and a remote member of the PcyA family were investigated to gain insights into catalysis. Biochemical data in combination with the PcyX crystal structure revealed that a conserved aspartate-histidine pair is critical for activity. Interestingly, the same residues are part of a catalytic Asp-His-Glu triad in PcyA, including an additional Glu. While this Glu residue is replaced by Asp in PcyX, it is not involved in catalysis. Substitution back to a Glu failed to convert PcyX to a PcyA. Therefore, the change in regiospecificity is not only caused by individual catalytic amino acid residues. Rather the combination of the architecture of the active site with the positioning of the substrate triggers specific proton transfer yielding the individual phycobilin products.

ENZYMES

Suggested EC number for PcyX: 1.3.7.6 DATABASES: The PcyX X-ray structure was deposited in the PDB with the accession code 5OWG.

The Existence of an Isolated Hydronium Ion in the Interior of Proteins

Neutron diffraction analysis studies reported an isolated hydronium ion (H₃O⁺) in the interior of d‐xylose isomerase (XI) and phycocyanobilin‐ferredoxin oxidoreductase (PcyA). H₃O⁺ forms hydrogen bonds (H‐bonds) with two histidine side‐chains and a backbone carbonyl group in PcyA, whereas H₃O⁺ forms H‐bonds with three acidic residues in XI. Using a quantum mechanical/molecular mechanical (QM/MM) approach, we analyzed stabilization of H₃O⁺ by the protein environment. QM/MM calculations indicated that H₃O⁺ was unstable in the PcyA crystal structure, releasing a proton to an H‐bond partner His88, producing H₂O and protonated His88. On the other hand, H₃O⁺ was stable in the XI crystal structure. H‐bond partners of isolated H₃O⁺ would be practically limited to acidic residues such as aspartic and glutamic acids in the protein environment.

A substrate-bound structure of cyanobacterial biliverdin reductase identifies stacked substrates as critical for activity

Biliverdin reductase catalyses the last step in haem degradation and produces the major lipophilic antioxidant bilirubin via reduction of biliverdin, using NAD(P)H as a cofactor. Despite the importance of biliverdin reductase in maintaining the redox balance, the molecular details of the reaction it catalyses remain unknown. Here we present the crystal structure of biliverdin reductase in complex with biliverdin and NADP+. Unexpectedly, two biliverdin molecules, which we designated the proximal and distal biliverdins, bind with stacked geometry in the active site. The nicotinamide ring of the NADP+ is located close to the reaction site on the proximal biliverdin, supporting that the hydride directly attacks this position of the proximal biliverdin. The results of mutagenesis studies suggest that a conserved Arg185 is essential for the catalysis. The distal biliverdin probably acts as a conduit to deliver the proton from Arg185 to the proximal biliverdin, thus yielding bilirubin.

Application of profile fitting method to neutron time-of-flight protein single crystal diffraction data collected at the iBIX

We developed and employed a profile fitting method for the peak integration of neutron time-of-flight diffraction data collected by the IBARAKI Biological Crystal Diffractometer (iBIX) at the Japan Proton Accelerator Research Complex (J-PARC) for protein ribonuclease A and α-thrombin single crystals. In order to determine proper fitting functions, four asymmetric functions were evaluated using strong intensity peaks. A Gaussian convolved with two back-to-back exponentials was selected as the most suitable fitting function, and a profile fitting algorithm for the integration method was developed. The intensity and structure refinement data statistics of the profile fitting method were compared to those of the summation integration method. It was clearly demonstrated that the profile fitting method provides more accurate integrated intensities and model structures than the summation integration method at higher resolution shells. The integration component with the profile fitting method has already been implemented in the iBIX data processing software STARGazer and its user manual has been prepared.

The use of neutron scattering to determine the functional structure of glycoside hydrolase

Neutron diffraction provides different information from X-ray diffraction, because neutrons are scattered by atomic nuclei, whereas X-rays are scattered by electrons. One of the key advantages of neutron crystallography is the ability to visualize hydrogen and deuterium atoms, making it possible to observe the protonation state of amino acid residues, hydrogen bonds, networks of water molecules and proton relay pathways in enzymes. But, because of technical difficulties, less than 100 enzyme structures have been evaluated by neutron crystallography to date. In this review, we discuss the advantages and disadvantages of neutron crystallography as a tool to investigate the functional structure of glycoside hydrolases, with some examples.

Long-Range Electrostatics-Induced Two-Proton Transfer Captured by Neutron Crystallography in an Enzyme Catalytic Site

Neutron crystallography was used to directly locate two protons before and after a pH-induced two-proton transfer between catalytic aspartic acid residues and the hydroxy group of the bound clinical drug darunavir, located in the catalytic site of enzyme HIV-1 protease. The two-proton transfer is triggered by electrostatic effects arising from protonation state changes of surface residues far from the active site. The mechanism and pH effect are supported by quantum mechanics/molecular mechanics (QM/MM) calculations. The low-pH proton configuration in the catalytic site is deemed critical for the catalytic action of this enzyme and may apply more generally to other aspartic proteases. Neutrons therefore represent a superb probe to obtain structural details for proton transfer reactions in biological systems at a truly atomic level.

Sub-atomic resolution X-ray crystallography and neutron crystallography: promise, challenges and potential

The International Year of Crystallography saw the number of macromolecular structures deposited in the Protein Data Bank cross the 100000 mark, with more than 90000 of these provided by X-ray crystallography. The number of X-ray structures determined to sub-atomic resolution (i.e. ≤1 Å) has passed 600 and this is likely to continue to grow rapidly with diffraction-limited synchrotron radiation sources such as MAX-IV (Sweden) and Sirius (Brazil) under construction. A dozen X-ray structures have been deposited to ultra-high resolution (i.e. ≤0.7 Å), for which precise electron density can be exploited to obtain charge density and provide information on the bonding character of catalytic or electron transfer sites. Although the development of neutron macromolecular crystallography over the years has been far less pronounced, and its application much less widespread, the availability of new and improved instrumentation, combined with dedicated deuteration facilities, are beginning to transform the field. Of the 83 macromolecular structures deposited with neutron diffraction data, more than half (49/83, 59%) were released since 2010. Sub-mm(3) crystals are now regularly being used for data collection, structures have been determined to atomic resolution for a few small proteins, and much larger unit-cell systems (cell edges >100 Å) are being successfully studied. While some details relating to H-atom positions are tractable with X-ray crystallography at sub-atomic resolution, the mobility of certain H atoms precludes them from being located. In addition, highly polarized H atoms and protons (H(+)) remain invisible with X-rays. Moreover, the majority of X-ray structures are determined from cryo-cooled crystals at 100 K, and, although radiation damage can be strongly controlled, especially since the advent of shutterless fast detectors, and by using limited doses and crystal translation at micro-focus beams, radiation damage can still take place. Neutron crystallography therefore remains the only approach where diffraction data can be collected at room temperature without radiation damage issues and the only approach to locate mobile or highly polarized H atoms and protons. Here a review of the current status of sub-atomic X-ray and neutron macromolecular crystallography is given and future prospects for combined approaches are outlined. New results from two metalloproteins, copper nitrite reductase and cytochrome c', are also included, which illustrate the type of information that can be obtained from sub-atomic-resolution (∼0.8 Å) X-ray structures, while also highlighting the need for complementary neutron studies that can provide details of H atoms not provided by X-ray crystallography.