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

Revealing protonation states and tracking substrate in serine hydroxymethyltransferase with room-temperature X-ray and neutron crystallography

Pyridoxal 5'-phosphate (PLP)-dependent enzymes utilize a vitamin B6-derived cofactor to perform a myriad of chemical transformations on amino acids and other small molecules. Some PLP-dependent enzymes, such as serine hydroxymethyltransferase (SHMT), are promising drug targets for the design of small-molecule antimicrobials and anticancer therapeutics, while others have been used to synthesize pharmaceutical building blocks. Understanding PLP-dependent catalysis and the reaction specificity is crucial to advance structure-assisted drug design and enzyme engineering. Here we report the direct determination of the protonation states in the active site of Thermus thermophilus SHMT (TthSHMT) in the internal aldimine state using room-temperature joint X-ray/neutron crystallography. Conserved active site architecture of the model enzyme TthSHMT and of human mitochondrial SHMT (hSHMT2) were compared by obtaining a room-temperature X-ray structure of hSHMT2, suggesting identical protonation states in the human enzyme. The amino acid substrate serine pathway through the TthSHMT active site cavity was tracked, revealing the peripheral and cationic binding sites that correspond to the pre-Michaelis and pseudo-Michaelis complexes, respectively. At the peripheral binding site, the substrate is bound in the zwitterionic form. By analyzing the observed protonation states, Glu53, but not His residues, is proposed as the general base catalyst, orchestrating the retro-aldol transformation of L-serine into glycine.

Capture of activated dioxygen intermediates at the copper-active site of a lytic polysaccharide monooxygenase

Metalloproteins perform a diverse array of redox-related reactions facilitated by the increased chemical functionality afforded by their metallocofactors. Lytic polysaccharide monooxygenases (LPMOs) are a class of copper-dependent enzymes that are responsible for the breakdown of recalcitrant polysaccharides via oxidative cleavage at the glycosidic bond. The activated copper-oxygen intermediates and their mechanism of formation remains to be established. Neutron protein crystallography which permits direct visualization of protonation states was used to investigate the initial steps of oxygen activation directly following active site copper reduction in Neurospora crassa LPMO9D. Herein, we cryo-trap an activated dioxygen intermediate in a mixture of superoxo and hydroperoxo states, and we identify the conserved second coordination shell residue His157 as the proton donor. Density functional theory calculations indicate that both superoxo and hydroperoxo active site states are stable. The hydroperoxo formed is potentially an early LPMO catalytic reaction intermediate or the first step in the mechanism of hydrogen peroxide formation in the absence of substrate. We observe that the N-terminal amino group of the copper coordinating His1 remains doubly protonated directly following molecular oxygen reduction by copper. Aided by molecular dynamics and mining minima free energy calculations we establish that the conserved second-shell His161 in MtPMO3* displays conformational flexibility in solution and that this flexibility is also observed, though to a lesser extent, in His157 of NcLPMO9D. The imidazolate form of His157 observed in our structure following oxygen intermediate protonation can be attributed to abolished His157 flexibility due steric hindrance in the crystal as well as the solvent-occluded active site environment due to crystal packing. A neutron crystal structure of NcLPMO9D at low pH further supports occlusion of the active site since His157 remains singly protonated even at acidic conditions.

Hybrid QM/MM Free-Energy Evaluation of Drug-Resistant Mutational Effect on the Binding of an Inhibitor Indinavir to HIV-1 Protease

A human immunodeficiency virus-1 (HIV-1) protease is a homodimeric aspartic protease essential for the replication of HIV. The HIV-1 protease is a target protein in drug discovery for antiretroviral therapy, and various inhibitor molecules of transition state analogues have been developed. However, serious drug-resistant mutants have emerged. For understanding the molecular mechanism of the drug resistance, an accurate examination of the impacts of the mutations on ligand binding and enzymatic activity is necessary. Here, we present a molecular simulation study on the ligand binding of indinavir, a potent transition state analogue inhibitor, to the wild-type protein and a V82T/I84V drug-resistant mutant of the HIV-1 protease. We employed a hybrid ab initio quantum mechanical/molecular mechanical (QM/MM) free-energy optimization technique which combines a highly accurate QM description of the ligand molecule and its interaction with statistically ample conformational sampling of the MM protein environment by long-time molecular dynamics simulations. Through the free-energy calculations of protonation states of catalytic groups at the binding pocket and of the ligand-binding affinity changes upon the mutations, we successfully reproduced the experimentally observed significant reduction of the binding affinity upon the drug-resistant mutations and elucidated the underlying molecular mechanism. The present study opens the way for understanding the molecular mechanism of drug resistance through the direct quantitative comparison of ligand binding and enzymatic reaction with the same accuracy.

Joint neutron/molecular dynamics vibrational spectroscopy reveals softening of HIV-1 protease upon binding of a tight inhibitor

Biomacromolecules are inherently dynamic, and their dynamics are interwoven into function. The fast collective vibrational dynamics in proteins occurs in the low picosecond timescale corresponding to frequencies of ∼5-50 cm^-1. This sub-to-low THz frequency regime covers the low-amplitude collective breathing motions of a whole protein and vibrations of the constituent secondary structure elements, such as α-helices, β-sheets and loops. We have used inelastic neutron scattering experiments in combination with molecular dynamics simulations to demonstrate the vibrational dynamics softening of HIV-1 protease, a target of HIV/AIDS antivirals, upon binding of a tight clinical inhibitor darunavir. Changes in the vibrational density of states of matching structural elements in the two monomers of the homodimeric protein are not identical, indicating asymmetric effects of darunavir on the vibrational dynamics. Three of the 11 major secondary structure elements contribute over 40% to the overall changes in the vibrational density of states upon darunavir binding. Molecular dynamics simulations informed by experiments allowed us to estimate that the altered vibrational dynamics of the protease would contribute -3.6 kcal mol^-1 at 300 K, or 25%, to the free energy of darunavir binding. As HIV-1 protease drug resistance remains a concern, our results open a new avenue to help establish a direct quantitative link between protein vibrational dynamics and drug resistance.

An N⋯H⋯N low-barrier hydrogen bond preorganizes the catalytic site of aspartate aminotransferase to facilitate the second half-reaction

Pyridoxal 5′-phosphate (PLP)-dependent enzymes have been extensively studied for their ability to fine-tune PLP cofactor electronics to promote a wide array of chemistries. Neutron crystallography offers a straightforward approach to studying the electronic states of PLP and the electrostatics of enzyme active sites, responsible for the reaction specificities, by enabling direct visualization of hydrogen atom positions. Here we report a room-temperature joint X-ray/neutron structure of aspartate aminotransferase (AAT) with pyridoxamine 5′-phosphate (PMP), the cofactor product of the first half reaction catalyzed by the enzyme. Between PMP N_SB and catalytic Lys258 Nζ amino groups an equally shared deuterium is observed in an apparent low-barrier hydrogen bond (LBHB). Density functional theory calculations were performed to provide further evidence of this LBHB interaction. The structural arrangement and the juxtaposition of PMP and Lys258, facilitated by the LBHB, suggests active site preorganization for the incoming ketoacid substrate that initiates the second half-reaction.

The neutron structure of pyridoxal 5′-phosphate-dependent enzyme aspartate aminotransferase with pyridoxamine 5′-phosphate (PMP) reveals a low-barrier hydrogen bond between the amino groups of PMP and catalytic Lys258, preorganizing the active site for catalysis

Capturing the Catalytic Proton of Dihydrofolate Reductase: Implications for General Acid-Base Catalysis

Acid–base catalysis, which involves one or more proton transfer reactions, is a chemical mechanism commonly employed by many enzymes. The molecular basis for catalysis is often derived from structures determined at the optimal pH for enzyme activity. However, direct observation of protons from experimental structures is quite difficult; thus, a complete mechanistic description for most enzymes remains lacking. Dihydrofolate reductase (DHFR) exemplifies general acid–base catalysis, requiring hydride transfer and protonation of its substrate, DHF, to form the product, tetrahydrofolate (THF). Previous X-ray and neutron crystal structures coupled with theoretical calculations have proposed that solvent mediates the protonation step. However, visualization of a proton transfer has been elusive. Based on a 2.1 Å resolution neutron structure of a pseudo-Michaelis complex of E. coli DHFR determined at acidic pH, we report the direct observation of the catalytic proton and its parent solvent molecule. Comparison of X-ray and neutron structures elucidated at acidic and neutral pH reveals dampened dynamics at acidic pH, even for the regulatory Met20 loop. Guided by the structures and calculations, we propose a mechanism where dynamics are crucial for solvent entry and protonation of substrate. This mechanism invokes the release of a sole proton from a hydronium (H₃O⁺) ion, its pathway through a narrow channel that sterically hinders the passage of water, and the ultimate protonation of DHF at the N5 atom.

Preliminary results of neutron and X-ray diffraction data collection on a lytic polysaccharide monooxygenase under reduced and acidic conditions

Lytic polysaccharide monooxygenases (LPMOs) are copper-center enzymes that are involved in the oxidative cleavage of the glycosidic bond in crystalline cellulose and other polysaccharides. The LPMO reaction is initiated by the addition of a reductant and oxygen to ultimately form an unknown activated copper-oxygen species that is responsible for polysaccharide-substrate H-atom abstraction. Given the sensitivity of metalloproteins to radiation damage, neutron protein crystallography provides a nondestructive technique for structural characterization while also informing on the positions of H atoms. Neutron cryo-crystallography permits the trapping of catalytic intermediates, thereby providing insight into the protonation states and chemical nature of otherwise short-lived species in the reaction mechanism. To characterize the reaction-mechanism intermediates of LPMO9D from Neurospora crassa, a cryo-neutron diffraction data set was collected from an ascorbate-reduced crystal. A second neutron diffraction data set was collected at room temperature from an LPMO9D crystal exposed to low-pH conditions to probe the protonation states of ionizable groups involved in catalysis under acidic conditions.

Mechanisms of Proteolytic Enzymes and Their Inhibition in QM/MM Studies

Experimental evidence for enzymatic mechanisms is often scarce, and in many cases inadvertently biased by the employed methods. Thus, apparently contradictory model mechanisms can result in decade long discussions about the correct interpretation of data and the true theory behind it. However, often such opposing views turn out to be special cases of a more comprehensive and superior concept. Molecular dynamics (MD) and the more advanced molecular mechanical and quantum mechanical approach (QM/MM) provide a relatively consistent framework to treat enzymatic mechanisms, in particular, the activity of proteolytic enzymes. In line with this, computational chemistry based on experimental structures came up with studies on all major protease classes in recent years; examples of aspartic, metallo-, cysteine, serine, and threonine protease mechanisms are well founded on corresponding standards. In addition, experimental evidence from enzyme kinetics, structural research, and various other methods supports the described calculated mechanisms. One step beyond is the application of this information to the design of new and powerful inhibitors of disease-related enzymes, such as the HIV protease. In this overview, a few examples demonstrate the high potential of the QM/MM approach for sophisticated pharmaceutical compound design and supporting functions in the analysis of biomolecular structures.

On the Case of the Misplaced Hydrogens

Few other elements play a more central role in biology than hydrogen. The interactions, bonding and movement of hydrogen atoms are central to biological catalysis, structure and function. Yet owing to the elusive nature of a single hydrogen atom few experimental and computational techniques can precisely determine its location. This is exemplified in short hydrogen bonds (SHBs) where the location of the hydrogen atom is indicative of the underlying strength of the bonds, which can vary from 1-5 kcal/mol in canonical hydrogen bonds, to an almost covalent nature in single-well hydrogen bonds. Owing to the often-times inferred position of hydrogen, the role of SHBs in biology has remained highly contested and debated. This has also led to discrepancies in computational, biochemical and structural studies of proteins thought to use SHBs in performing chemistry and stabilizing interactions. Herein, we discuss in detail two distinct examples, namely the conserved catalytic triad and the photoreceptor, photoactive yellow protein, where studies of these SHB-containing systems have permitted contextualization of the role these unique hydrogen bonds play in biology.

Room-temperature neutron and X-ray data collection of 3CL Mpro from SARS-CoV-2

The replication of SARS-CoV-2 produces two large polyproteins, pp1a and pp1ab, that are inactive until cleavage by the viral chymotrypsin-like cysteine protease enzyme (3CL Mpro) into a series of smaller functional proteins. At the heart of 3CL Mpro is an unusual catalytic dyad formed by the side chains of His41 and Cys145 and a coordinated water molecule. The catalytic mechanism by which the enzyme operates is still unknown, as crucial information on the protonation states within the active site is unclear. To experimentally determine the protonation states of the catalytic site and of the other residues in the substrate-binding cavity, and to visualize the hydrogen-bonding networks throughout the enzyme, room-temperature neutron and X-ray data were collected from a large H/D-exchanged crystal of ligand-free (apo) 3CL Mpro.

Structural Analysis of Potent Hybrid HIV-1 Protease Inhibitors Containing Bis-tetrahydrofuran in a Pseudosymmetric Dipeptide Isostere

The design, synthesis, and X-ray structural analysis of hybrid HIV-1 protease inhibitors (PIs) containing bis-tetrahydrofuran (bis-THF) in a pseudo-C₂-symmetric dipeptide isostere are described. A series of PIs were synthesized by incorporating bis-THF of darunavir on either side of the Phe-Phe isostere of lopinavir in combination with hydrophobic amino acids on the opposite P2/P2' position. Structure-activity relationship studies indicated that the bis-THF moiety can be attached at either the P2 or P2' position without significantly affecting potency. However, the group on the opposite P2/P2' position had a dramatic effect on potency depending on the size and shape of the side chain. Cocrystal structures of inhibitors with wild-type HIV-1 protease revealed that the bis-THF moiety retained similar interactions as observed in the darunavir-protease complex regardless of the position on the Phe-Phe isostere. Analyses of cocrystal structures and molecular dynamics simulations provide insights into optimizing HIV-1 PIs containing bis-THF in non-sulfonamide dipeptide isosteres.

Highly drug-resistant HIV-1 protease reveals decreased intra-subunit interactions due to clusters of mutations

Drug-resistance is a serious problem for treatment of the HIV/AIDS pandemic. Potent clinical inhibitors of HIV-1 protease show several orders of magnitude worse inhibition of highly drug-resistant variants. Hence, the structure and enzyme activities were analyzed for HIV protease mutant HIV-1 protease (EC 3.4.23.16) (PR) with 22 mutations (PRS5B) from a clinical isolate that was selected by machine learning to represent high-level drug-resistance. PRS5B has 22 mutations including only one (I84V) in the inhibitor binding site; however, clinical inhibitors had poor inhibition of PRS5B activity with kinetic inhibition value (K_i ) values of 4-1000 nm or 18- to 8000-fold worse than for wild-type PR. High-resolution crystal structures of PRS5B complexes with the best inhibitors, amprenavir (APV) and darunavir (DRV) (K_i ~ 4 nm), revealed only minor changes in protease-inhibitor interactions. Instead, two distinct clusters of mutations in distal regions induce coordinated conformational changes that decrease favorable internal interactions across the entire protein subunit. The largest structural rearrangements are described and compared to other characterized resistant mutants. In the protease hinge region, the N83D mutation eliminates a hydrogen bond connecting the hinge and core of the protease and increases disorder compared to highly resistant mutants PR with 17 mutations and PR with 20 mutations with similar hinge mutations. In a distal β-sheet, mutations G73T and A71V coordinate with accessory mutations to bring about shifts that propagate throughout the subunit. Molecular dynamics simulations of ligand-free dimers show differences consistent with loss of interactions in mutant compared to wild-type PR. Clusters of mutations exhibit both coordinated and antagonistic effects, suggesting PRS5B may represent an intermediate stage in the evolution of more highly resistant variants. DATABASES: Structural data are available in Protein Data Bank under the accession codes 6P9A and 6P9B for PRS5B/DRV and PRS5B/APV, respectively.

Visualizing Tetrahedral Oxyanion Bound in HIV-1 Protease Using Neutrons: Implications for the Catalytic Mechanism and Drug Design

HIV-1 protease is indispensable for virus propagation and an important therapeutic target for antiviral inhibitors to treat AIDS. As such inhibitors are transition-state mimics, a detailed understanding of the enzyme mechanism is crucial for the development of better anti-HIV drugs. Here, we used room-temperature joint X-ray/neutron crystallography to directly visualize hydrogen atoms and map hydrogen bonding interactions in a protease complex with peptidomimetic inhibitor KVS-1 containing a reactive nonhydrolyzable ketomethylene isostere, which, upon reacting with the catalytic water molecule, is converted into a tetrahedral intermediate state, KVS-1_TI. We unambiguously determined that the resulting tetrahedral intermediate is an oxyanion, rather than the gem-diol, and both catalytic aspartic acid residues are protonated. The oxyanion tetrahedral intermediate appears to be unstable, even though the negative charge on the oxyanion is delocalized through a strong n → π* hyperconjugative interaction into the nearby peptidic carbonyl group of the inhibitor. To better understand the influence of the ketomethylene isostere as a protease inhibitor, we have also examined the protease structure and binding affinity with keto-darunavir (keto-DRV), which similar to KVS-1 includes the ketomethylene isostere. We show that keto-DRV is a significantly less potent protease inhibitor than DRV. These findings shed light on the reaction mechanism of peptide hydrolysis catalyzed by HIV-1 protease and provide valuable insights into further improvements in the design of protease inhibitors.

Neutron crystallography of copper amine oxidase reveals keto/enolate interconversion of the quinone cofactor and unusual proton sharing

Recent advances in neutron crystallographic studies have provided structural bases for quantum behaviors of protons observed in enzymatic reactions. Thus, we resolved the neutron crystal structure of a bacterial copper (Cu) amine oxidase (CAO), which contains a prosthetic Cu ion and a protein-derived redox cofactor, topa quinone (TPQ). We solved hitherto unknown structures of the active site, including a keto/enolate equilibrium of the cofactor with a nonplanar quinone ring, unusual proton sharing between the cofactor and the catalytic base, and metal-induced deprotonation of a histidine residue that coordinates to the Cu. Our findings show a refined active-site structure that gives detailed information on the protonation state of dissociable groups, such as the quinone cofactor, which are critical for catalytic reactions.

Pyridoxal 5'-phosphate dependent reactions: Analyzing the mechanism of aspartate aminotransferase

Enzyme catalysis is the primary activity in energy and information metabolism and enzyme cofactors are key to the catalytic ability of most enzymes. Pyridoxal 5'-phosphate (PLP) cofactor, derived from Vitamin B6, is widely distributed in nature and has significant latitude in catalytic diversity. X-ray crystallography has revealed the structures of diverse PLP dependent enzymes from multiple families. But these structures are incomplete, lacking the positions of protons essential for understanding enzymatic mechanisms. Here, we review the diversity of PLP and discuss the use of neutron crystallography and joint X-ray/neutron refinement of Fold Type I aspartate aminotransferase to visualize the positions of protons in both the internal and external aldimine forms. Strategies used to prepare extremely large crystals required for neutron diffraction and the approach to data refinement including the PLP cofactor are discussed. The observed positions of protons, including one located in a previously unknown low-barrier hydrogen bond, have been used to create more accurate models for computational analysis. The results revealed a new mechanism for the transaminase reaction where hyperconjugation is key to reducing the energy barrier which finally provides a clear explanation of the Dunathan alignment.

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.

Proton transfer and drug binding details revealed in neutron diffraction studies of wild-type and drug resistant HIV-1 protease

HIV-1 protease is an essential therapeutic target for the design and development of antiviral inhibitors to treat AIDS. We used room temperature neutron crystallography to accurately determine hydrogen atom positions in several protease complexes with clinical drugs, amprenavir and darunavir. Hydrogen bonding interactions were carefully mapped to provide an unprecedented picture of drug binding to the protease target. We demonstrate that hydrogen atom positions within the enzyme catalytic site can be altered by introducing drug resistant mutations and by protonating surface residues that trigger proton transfer reactions between the catalytic Asp residues and the hydroxyl group of darunavir. When protein perdeuteration is not feasible, we validate the use of initial H/D exchange with unfolded protein and partial deuteration in pure D₂O with hydrogenous glycerol to maximize deuterium incorporation into the protein, with no detrimental effects on the growth of quality crystals suitable for neutron diffraction experiments.

Titration of ionizable groups in proteins using multiple neutron data sets from a single crystal: application to the small GTPase Ras

Neutron protein crystallography (NPC) reveals the three-dimensional structures of proteins, including the positions of H atoms. The technique is particularly suited to elucidate ambiguous catalytic steps in complex biochemical reactions. While NPC uniquely complements biochemical assays and X-ray structural analyses by revealing the protonation states of ionizable groups at and around the active site of enzymes, the technique suffers from a major drawback: large single crystals must be grown to compensate for the relatively low flux of neutron beams. However, in addition to revealing the positions of hydrogens involved in enzyme catalysis, NPC has the advantage over X-ray crystallography that the crystals do not suffer radiation damage. The lack of radiation damage can be exploited to conduct in crystallo parametric studies. Here, the use of a single crystal of the small GTPase Ras to collect three neutron data sets at pD 8.4, 9.0 and 9.4 is reported, enabling an in crystallo titration study using NPC. In addition to revealing the behavior of titratable groups in the active site, the data sets will allow the analysis of allosteric water-mediated communication networks across the molecule, particularly regarding Cys118 and three tyrosine residues central to these networks, Tyr32, Tyr96 and Tyr137, with pK_a values expected to be in the range sampled in our experiments.

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.

The Neutron Macromolecular Crystallography Instruments at Oak Ridge National Laboratory: Advances, Challenges, and Opportunities

The IMAGINE and MaNDi instruments, located at Oak Ridge National Laboratory High Flux Isotope Reactor and Spallation Neutron Source, respectively, are powerful tools for determining the positions of hydrogen atoms in biological macromolecules and their ligands, orienting water molecules, and for differentiating chemical states in macromolecular structures. The possibility to model hydrogen and deuterium atoms in neutron structures arises from the strong interaction of neutrons with the nuclei of these isotopes. Positions can be unambiguously assigned from diffraction studies at the 1.5–2.5 Å resolutions, which are typical for protein crystals. Neutrons have the additional benefit for structural biology of not inducing radiation damage to protein crystals, which can be critical in the study of metalloproteins. Here we review the specifications of the IMAGINE and MaNDi beamlines and illustrate their complementarity. IMAGINE is suitable for crystals with unit cell edges up to 150 Å using a quasi-Laue technique, whereas MaNDi provides neutron crystallography resources for large unit cell samples with unit cell edges up to 300 Å using the time of flight (TOF) Laue technique. The microbial culture and crystal growth facilities which support the IMAGINE and MaNDi user programs are also described.

Drug Resistance Mutation L76V Alters Nonpolar Interactions at the Flap-Core Interface of HIV-1 Protease

Four HIV-1 protease (PR) inhibitors, clinical inhibitors lopinavir and tipranavir, and two investigational compounds 4 and 5, were studied for their effect on the structure and activity of PR with drug-resistant mutation L76V (PR_L76V). Compound 5 exhibited the best K _i value of 1.9 nM for PR_L76V, whereas the other three inhibitors had K _i values of 4.5-7.6 nM, 2-3 orders of magnitude worse than for wild-type enzymes. Crystal structures showed only minor differences in interactions of inhibitors with PR_L76V compared to wild-type complexes. The shorter side chain of Val76 in the mutant lost hydrophobic interactions with Lys45 and Ile47 in the flap, and with Asp30 and Thr74 in the protein core, consistent with decreased stability. Inhibitors forming additional polar interactions with the flaps or dimer interface of PR_L76V were unable to compensate for the decrease in internal hydrophobic contacts. These structures provide insights for inhibitor design.

IMAGINE: neutrons reveal enzyme chemistry

Neutron diffraction is exquisitely sensitive to the positions of H atoms in protein crystal structures. IMAGINE is a high-intensity, quasi-Laue neutron crystallography beamline developed at the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory. This state-of-the-art facility for neutron diffraction has enabled detailed structural analysis of macromolecules. IMAGINE is especially suited to resolve individual H atoms in protein structures, enabling neutron protein structures to be determined at or near atomic resolutions from crystals with volumes of less than 1 mm3 and unit-cell edges of less than 150 Å. Beamline features include elliptical focusing mirrors that deliver neutrons into a 2.0 × 3.2 mm focal spot at the sample position, and variable short- and long-wavelength cutoff optics that provide automated exchange between multiple wavelength configurations. This review gives an overview of the IMAGINE beamline at the HFIR, presents examples of the scientific questions being addressed at this beamline, and highlights important findings in enzyme chemistry that have been made using the neutron diffraction capabilities offered by IMAGINE.

Consistent Prediction of Mutation Effect on Drug Binding in HIV-1 Protease Using Alchemical Calculations

Despite a large number of antiretroviral drugs targeting HIV-1 protease for inhibition, mutations in this protein during the course of patient treatment can render them inefficient. This emerging resistance inspired numerous computational studies of the HIV-1 protease aimed at predicting the effect of mutations on drug binding in terms of free binding energy Δ G, as well as in mechanistic terms. In this study, we analyze ten different protease-inhibitor complexes carrying major resistance-associated mutations (RAMs) G48V, I50V, and L90M using molecular dynamics simulations. We demonstrate that alchemical free energy calculations can consistently predict the effect of mutations on drug binding. By explicitly probing different protonation states of the catalytic aspartic dyad, we reveal the importance of the correct choice of protonation state for the accuracy of the result. We also provide insight into how different mutations affect drug binding in their specific ways, with the unifying theme of how all of them affect the crucial drug binding regions of the protease.

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.

Proton in a Confined Space: Structural Studies of H+ ⊂Crypt-111 Iodide and Some Halogen-Bonded Derivatives

Experimental observations and modeling data are reported on the solid-state structural features of crypt- 111⋅HI (1) and the three-component co-crystals that 1 forms with α,ω-diiodoperfluoroalkanes 2 a-d. X-ray analyses indicate that, in all five systems and at low temperature, the caged proton is covalently bonded to a single nitrogen atom and is involved in a network of intramolecular hydrogen bonds. In contrast, room-temperature, solid-state ¹⁵ N NMR spectroscopy suggests magnetic equivalency of the two N atoms of crypt-111 in both 1 and co-crystals of 1 with diiodoperfluoroalkanes. Computational modelling confirms that the acidic hydrogen inside the cavity preferentially sits along the internitrogen axis and is covalently bonded to one nitrogen. The computed energy barriers suggest that the hopping of the encapsulated proton between the two N atoms of the cage can occur in the halogen-bonded co-crystals of 1⋅2, but it is hardly possible in the pure H⁺ ⊂crypt-111 iodide 1. These different pictures of the proton position and dynamics obtained by using different techniques and conditions confirm the unique characteristics of the confined space within the cavity of crypr-111 and the distinctive features of processes occurring therein.

Direct visualization of critical hydrogen atoms in a pyridoxal 5'-phosphate enzyme

Enzymes dependent on pyridoxal 5′-phosphate (PLP, the active form of vitamin B₆) perform a myriad of diverse chemical transformations. They promote various reactions by modulating the electronic states of PLP through weak interactions in the active site. Neutron crystallography has the unique ability of visualizing the nuclear positions of hydrogen atoms in macromolecules. Here we present a room-temperature neutron structure of a homodimeric PLP-dependent enzyme, aspartate aminotransferase, which was reacted in situ with α-methylaspartate. In one monomer, the PLP remained as an internal aldimine with a deprotonated Schiff base. In the second monomer, the external aldimine formed with the substrate analog. We observe a deuterium equidistant between the Schiff base and the C-terminal carboxylate of the substrate, a position indicative of a low-barrier hydrogen bond. Quantum chemical calculations and a low-pH room-temperature X-ray structure provide insight into the physical phenomena that control the electronic modulation in aspartate aminotransferase.

Pyridoxal 5’-phosphate (PLP) is a ubiquitous co factor for diverse enzymes, among them aspartate aminotransferase. Here the authors use neutron crystallography, which allows the visualization of the positions of hydrogen atoms, and computation to characterize the catalytic mechanism of the enzyme.

High Pressure Spectroscopic Investigation on Proton Transfer in Squaric Acid and 4,4'-Bipyridine Co-crystal

In attempt to the obtain detailed geometric information of proton transfer compound (subsequently denote as SQBP) formed between squaric acid (SQ)and 4,4′-bipyridine(BP), and to investigate the mechanisms of pressure-induced double proton transfer and related structural phase transition, we carried out in-situ high pressure Raman spectroscopy of SQBP up to 20 GPa. A solid-solid phase transition together with double proton transfer phenomenon was confirmed by Raman spectroscopy at about 1.5 GPa, and the activation of C = O stretching mode in Raman spectra indicates a square-ring structure of SQ with four symmetric C = O bond formation. These results are further supported by first-principals calculations and in-situ high pressure infrared absorption spectroscopy. Additionally, Raman intensity analysis suggests that a higher-order phase transition with planar BP molecular structure occurred in the pressure range of 3~6 GPa. As a result, the π electron delocalization effect in BP dominated the intensity enhancement of C = O stretching mode in SQ. To the best of our knowledge, this is the first time observation of the intensity enhancement of proton donor’s normal modes induced by proton acceptor’s π electron delocalization.

Back-exchange of deuterium in neutron crystallography: characterization by IR spectroscopy

The application of IR spectroscopy to the characterization and quality control of samples used in neutron crystallography is described. While neutron crystallography is a growing field, the limited availability of neutron beamtime means that there may be a delay between crystallogenesis and data collection. Since essentially all neutron crystallographic work is carried out using D2O-based solvent buffers, a particular concern for these experiments is the possibility of H2O back-exchange across reservoir or capillary sealants. This may limit the quality of neutron scattering length density maps and of the associated analysis. Given the expense of central facility beamtime and the effort that goes into the production of suitably sized (usually perdeuterated) crystals, a systematic method of exploiting IR spectroscopy for the analysis of back-exchange phenomena in the reservoirs used for crystal growth is valuable. Examples are given in which the characterization of D2O/H2O back-exchange in transthyretin crystals is described.

Room Temperature Neutron Crystallography of Drug Resistant HIV-1 Protease Uncovers Limitations of X-ray Structural Analysis at 100 K

HIV-1 protease inhibitors are crucial for treatment of HIV-1/AIDS, but their effectiveness is thwarted by rapid emergence of drug resistance. To better understand binding of clinical inhibitors to resistant HIV-1 protease, we used room-temperature joint X-ray/neutron (XN) crystallography to obtain an atomic-resolution structure of the protease triple mutant (V32I/I47V/V82I) in complex with amprenavir. The XN structure reveals a D⁺ ion located midway between the inner Oδ1 oxygen atoms of the catalytic aspartic acid residues. Comparison of the current XN structure with our previous XN structure of the wild-type HIV-1 protease-amprenavir complex suggests that the three mutations do not significantly alter the drug-enzyme interactions. This is in contrast to the observations in previous 100 K X-ray structures of these complexes that indicated loss of interactions by the drug with the triple mutant protease. These findings, thus, uncover limitations of structural analysis of drug binding using X-ray structures obtained at 100 K.

Fifteen years of the Protein Crystallography Station: the coming of age of macromolecular neutron crystallography

The Protein Crystallography Station (PCS), located at the Los Alamos Neutron Scattering Center (LANSCE), was the first macromolecular crystallography beamline to be built at a spallation neutron source. Following testing and commissioning, the PCS user program was funded by the Biology and Environmental Research program of the Department of Energy Office of Science (DOE-OBER) for 13 years (2002-2014). The PCS remained the only dedicated macromolecular neutron crystallography station in North America until the construction and commissioning of the MaNDi and IMAGINE instruments at Oak Ridge National Laboratory, which started in 2012. The instrument produced a number of research and technical outcomes that have contributed to the field, clearly demonstrating the power of neutron crystallo-graphy in helping scientists to understand enzyme reaction mechanisms, hydrogen bonding and visualization of H-atom positions, which are critical to nearly all chemical reactions. During this period, neutron crystallography became a technique that increasingly gained traction, and became more integrated into macromolecular crystallography through software developments led by investigators at the PCS. This review highlights the contributions of the PCS to macromolecular neutron crystallography, and gives an overview of the history of neutron crystallography and the development of macromolecular neutron crystallography from the 1960s to the 1990s and onwards through the 2000s.

Dangerous liaisons: anion-induced protonation in phosphate–polyamine interactions and their implications for the charge states of biologically relevant surfaces

The specific binding of phosphates on surface amino groups has complex consequences on the protonation equilibriums and effective surface charges with serious implications for further biochemically relevant interactions.

Neutron structures of the Helicobacter pylori 5'-methylthioadenosine nucleosidase highlight proton sharing and protonation states

MTAN (5'-methylthioadenosine nucleosidase) catalyzes the hydrolysis of the N-ribosidic bond of a variety of adenosine-containing metabolites. The Helicobacter pylori MTAN (HpMTAN) hydrolyzes 6-amino-6-deoxyfutalosine in the second step of the alternative menaquinone biosynthetic pathway. Substrate binding of the adenine moiety is mediated almost exclusively by hydrogen bonds, and the proposed catalytic mechanism requires multiple proton-transfer events. Of particular interest is the protonation state of residue D198, which possesses a pK_a above 8 and functions as a general acid to initiate the enzymatic reaction. In this study we present three corefined neutron/X-ray crystal structures of wild-type HpMTAN cocrystallized with S-adenosylhomocysteine (SAH), Formycin A (FMA), and (3R,4S)-4-(4-Chlorophenylthiomethyl)-1-[(9-deaza-adenin-9-yl)methyl]-3-hydroxypyrrolidine (p-ClPh-Thio-DADMe-ImmA) as well as one neutron/X-ray crystal structure of an inactive variant (HpMTAN-D198N) cocrystallized with SAH. These results support a mechanism of D198 pKa elevation through the unexpected sharing of a proton with atom N7 of the adenine moiety possessing unconventional hydrogen-bond geometry. Additionally, the neutron structures also highlight active site features that promote the stabilization of the transition state and slight variations in these interactions that result in 100-fold difference in binding affinities between the DADMe-ImmA and ImmA analogs.

Neutron crystallography aids in drug design

Since drugs bind to their targets through directional H bonding and non-directional hydrophobic and electrostatic interactions, neutron crystallography can help guide structure-based drug design. This is illustrated by McKenna and co-workers [Aggarwal et al. (2016), IUCrJ, 3, 319-325] who describe the room-temperature neutron structure of human carbonic anyhydrase II in complex with the clinical inhibitor methazolamide to 2.2 Å resolution, and compare this with the previously determined room-temperature neutron structure of human carbonic anyhydrase II in complex with the clinical inhibitor acetazolamide to 2.0 Å resolution [Fisher et al. (2012). J. Am. Chem. Soc.134, 14726-14729].