Quantum model of catalysis based on a mobile proton revealed by subatomic x-ray and neutron diffraction studies of h-aldose reductase

Applying Network Pharmacology and Molecular Docking in the Screening for Molecular Mechanisms of Ampalaya (Momordica charantia L.) and Banaba (Lagerstroemia speciosa L.) against Type 2 Diabetes Mellitus

Background and Objectives

Type 2 diabetes mellitus (T2DM) is a global health concern affecting more than 400 million people worldwide. Diabetic neuropathy, nephropathy, retinopathy, and cardiovascular complications lead to debilitating effects to patients. To prevent these, the treatment goal is to lower the blood sugar levels and maintain at a normal range which is achieved through conventional treatments like insulin and oral hypoglycemic agents. However, the high cost of these medications implicates patient treatment outcomes. Hence, alternatives are sought for including the use of herbal medicines. Momordica charantia (MC) and Lagerstroemia speciosa (LS) are common herbal medicines used to manage T2DM. In the Philippines, these herbal preparations are validated for their glucose lowering effects and are commonly found in combination in food supplements. The study aims to screen the possible mechanisms of compounds present in these herbal medicines which can offer possible explanations for their synergistic effects and rationalization of their combination in preparations.

Methods

Network pharmacology was employed to determine pivotal proteins that are targeted by MC and LS compounds. Molecular docking was then done to evaluate the favorability of the binding of these compounds toward their target proteins.

Results

Our results showed that TNF, HSP90AA1, MAPK3, ALDH2, GCK, AKR1B1, TTR and RBP4 are the possible pivotal targets of MC and LS compounds in T2DM.

Conclusion

Terpenoids from MC and decanoic acid from LS are the compounds which showed favorable binding towards pivotal protein targets in T2DM. By binding towards the different key proteins in T2DM, they may exhibit their synergistic effects. However, the results of this study are bound to the limitations of computational methods and experimental validation are needed to verify our findings.

Exploring potent aldose reductase inhibitors for anti-diabetic (anti-hyperglycemic) therapy: integrating structure-based drug design, and MMGBSA approaches

Aldose reductase (AR) is an important target in the development of therapeutics against hyper-glycemia-induced health complications such as retinopathy, etc. In this study, we employed a combination of structure-based drug design, molecular simulation, and free energy calculation approaches to identify potential hit molecules against anti-diabetic (anti-hyperglycemic)-induced health complications. The 3D structure of aldoreductase was screened for multiple compound libraries (1,00,000 compounds) and identified as ZINC35671852, ZINC78774792 from the ZINC database, Diamino-di nitro-methyl dioctyl phthalate, and Penta-o-galloyl-glucose from the South African natural compounds database, and Bisindolylmethane thiosemi-carbazides and Bisindolylme-thane-hydrazone from the Inhouse database for this study. The mode of binding interactions of the selected compounds later predicted their aldose reductase inhibitory potential. These com-pounds interact with the key active site residues through hydrogen bonds, salt bridges, and π-π interactions. The structural dynamics and binding free energy results further revealed that these compounds possess stable dynamics with excellent binding free energy scores. The structures of the lead inhibitors can serve as templates for developing novel inhibitors, and in vitro testing to confirm their anti-diabetic potential is warranted. The current study is the first to design small molecule inhibitors for the aldoreductase protein that can be used in the development of therapeutic agents to treat diabetes.

Aldose reductase inhibition alleviates diabetic cardiomyopathy and is associated with a decrease in myocardial fatty acid oxidation

BACKGROUND

Cardiovascular diseases, including diabetic cardiomyopathy, are major causes of death in people with type 2 diabetes. Aldose reductase activity is enhanced in hyperglycemic conditions, leading to altered cardiac energy metabolism and deterioration of cardiac function with adverse remodeling. Because disturbances in cardiac energy metabolism can promote cardiac inefficiency, we hypothesized that aldose reductase inhibition may mitigate diabetic cardiomyopathy via normalization of cardiac energy metabolism.

METHODS

Male C57BL/6J mice (8-week-old) were subjected to experimental type 2 diabetes/diabetic cardiomyopathy (high-fat diet [60% kcal from lard] for 10 weeks with a single intraperitoneal injection of streptozotocin (75 mg/kg) at 4 weeks), following which animals were randomized to treatment with either vehicle or AT-001, a next-generation aldose reductase inhibitor (40 mg/kg/day) for 3 weeks. At study completion, hearts were perfused in the isolated working mode to assess energy metabolism.

RESULTS

Aldose reductase inhibition by AT-001 treatment improved diastolic function and cardiac efficiency in mice subjected to experimental type 2 diabetes. This attenuation of diabetic cardiomyopathy was associated with decreased myocardial fatty acid oxidation rates (1.15 ± 0.19 vs 0.5 ± 0.1 µmol min-1 g dry wt-1 in the presence of insulin) but no change in glucose oxidation rates compared to the control group. In addition, cardiac fibrosis and hypertrophy were also mitigated via AT-001 treatment in mice with diabetic cardiomyopathy.

CONCLUSIONS

Inhibiting aldose reductase activity ameliorates diastolic dysfunction in mice with experimental type 2 diabetes, which may be due to the decline in myocardial fatty acid oxidation, indicating that treatment with AT-001 may be a novel approach to alleviate diabetic cardiomyopathy in patients with diabetes.

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.

Current status of neutron crystallography in structural biology

Hydrogen atoms and hydration water molecules in proteins are essential for many biochemical processes, especially enzyme catalysis. Neutron crystallography enables direct observation of hydrogen atoms, and reveals molecular recognition through hydrogen bonding and catalytic reactions involving proton-coupled electron transfer. The use of neutron crystallography is still limited for proteins, but its popularity is increasing owing to an increase in the number of diffractometers for structural biology at neutron facilities and advances in sample preparation. According to the characteristics of the neutrons, monochromatic or quasi-Laue methods and the time-of-flight method are used in nuclear reactors and pulsed spallation sources, respectively, to collect diffraction data. Growing large crystals is an inevitable problem in neutron crystallography for structural biology, but sample deuteration, especially protein perdeuteration, is effective in reducing background levels, which shortens data collection time and decreases the crystal size required. This review also introduces our recent neutron structure analyses of copper amine oxidase and copper-containing nitrite reductase. The neutron structure of copper amine oxidase gives detailed information on the protonation state of dissociable groups, such as the quinone cofactor, which are critical for catalytic reactions. Electron transfer via a hydrogen-bond jump and a hydroxide ion ligation in copper-containing nitrite reductase are clarified, and these observations are consistent with the results from the quantum chemical calculations. This review article is an extended version of the Japanese article, Elucidation of Enzymatic Reaction Mechanism by Neutron Crystallography, published in SEIBUTSU-BUTSURI Vol. 61, p.216–222 (2021).

Perspective on the Structural Basis for Human Aldo-Keto Reductase 1B10 Inhibition

Human aldo-keto reductase 1B10 (AKR1B10) is overexpressed in many cancer types and is involved in chemoresistance. This makes AKR1B10 to be an interesting drug target and thus many enzyme inhibitors have been investigated. High-resolution crystallographic structures of AKR1B10 with various reversible inhibitors were deeply analyzed and compared to those of analogous complexes with aldose reductase (AR). In both enzymes, the active site included an anion-binding pocket and, in some cases, inhibitor binding caused the opening of a transient specificity pocket. Different structural conformers were revealed upon inhibitor binding, emphasizing the importance of the highly variable loops, which participate in the transient opening of additional binding subpockets. Two key differences between AKR1B10 and AR were observed regarding the role of external loops in inhibitor binding. The first corresponded to the alternative conformation of Trp112 (Trp111 in AR). The second difference dealt with loop A mobility, which defined a larger and more loosely packed subpocket in AKR1B10. From this analysis, the general features that a selective AKR1B10 inhibitor should comply with are the following: an anchoring moiety to the anion-binding pocket, keeping Trp112 in its native conformation (AKR1B10-like), and not opening the specificity pocket in AR.

Physiological and Pathological Roles of Aldose Reductase

Aldose reductase (AR) is an aldo-keto reductase that catalyzes the first step in the polyol pathway which converts glucose to sorbitol. Under normal glucose homeostasis the pathway represents a minor route of glucose metabolism that operates in parallel with glycolysis. However, during hyperglycemia the flux of glucose via the polyol pathway increases significantly, leading to excessive formation of sorbitol. The polyol pathway-driven accumulation of osmotically active sorbitol has been implicated in the development of secondary diabetic complications such as retinopathy, nephropathy, and neuropathy. Based on the notion that inhibition of AR could prevent these complications a range of AR inhibitors have been developed and tested; however, their clinical efficacy has been found to be marginal at best. Moreover, recent work has shown that AR participates in the detoxification of aldehydes that are derived from lipid peroxidation and their glutathione conjugates. Although in some contexts this antioxidant function of AR helps protect against tissue injury and dysfunction, the metabolic transformation of the glutathione conjugates of lipid peroxidation-derived aldehydes could also lead to the generation of reactive metabolites that can stimulate mitogenic or inflammatory signaling events. Thus, inhibition of AR could have both salutary and injurious outcomes. Nevertheless, accumulating evidence suggests that inhibition of AR could modify the effects of cardiovascular disease, asthma, neuropathy, sepsis, and cancer; therefore, additional work is required to selectively target AR inhibitors to specific disease states. Despite past challenges, we opine that a more gainful consideration of therapeutic modulation of AR activity awaits clearer identification of the specific role(s) of the AR enzyme in health and disease.

The role of hydrogen atoms in redox catalysis by the flavoenzyme cholesterol oxidase

Flavoenzymes comprise a large class of proteins that carry out a diverse range of important redox chemistry. Although X-ray crystal structures of many flavoenzymes have been determined, there are still unresolved questions regarding the actual oxidation state of the flavin cofactors in these structures due to photoreduction by the ionizing radiation of the X-ray beam during the diffraction experiment. Additionally, the ability to visualize hydrogen atoms in X-ray structures is difficult due to the weak scattering capability of these atoms. Since hydrogen atoms affect the electrostatic nature of enzyme active sites and play important roles in the chemistry of key amino acid residues, visualizing the precise positions of these atoms provides a more detailed understanding of their role in enzyme catalysis. Single crystal neutron diffraction is an alternative method to structure determination, circumventing problems associated with photoreduction of the sample thus providing a clearer view of the structural features of a flavoenzyme in different redox states. Additionally, the larger neutron scattering factors for hydrogen and deuterium atoms enables one to visualize these atoms much more easily than from X-ray scattering measurements. In this chapter we give an overview of neutron and X-ray crystallography studies on the flavoenzyme, cholesterol oxidase and how the observations of unusual hydrogen atom positions provide insight into the redox chemistry of the flavin cofactor.

Dynamic nuclear polarization enhanced neutron crystallography: Amplifying hydrogen in biological crystals

Dynamic nuclear polarization (DNP) can provide a powerful means to amplify neutron diffraction from biological crystals by 10-100-fold, while simultaneously enhancing the visibility of hydrogen by an order of magnitude. Polarizing the neutron beam and aligning the proton spins in a polarized sample modulates the coherent and incoherent neutron scattering cross-sections of hydrogen, in ideal cases amplifying the coherent scattering by almost an order of magnitude and suppressing the incoherent background to zero. This chapter describes current efforts to develop and apply DNP techniques for spin polarized neutron protein crystallography, highlighting concepts, experimental design, labeling strategies and recent results, as well as considering new strategies for data collection and analysis that these techniques could enable.

What are the current limits on determination of protonation state using neutron macromolecular crystallography?

The rate of deposition of models determined by neutron diffraction, or a hybrid approach that combines X-ray and neutron diffraction, has increased in recent years. The benefit of neutron diffraction is that hydrogen atom (H) positions are detectable, allowing for the determination of protonation state and water molecule orientation. This study analyses all neutron models deposited in the Protein Data Bank to date, focusing on protonation state and properties of H (or deuterium, D) atoms as well as the details of water molecules. In particular, clashes and hydrogen bonds involving H or D atoms are investigated. As water molecules are typically the least reproducible part of a structural model, their positions in neutron models were compared to those in homologous high-resolution X-ray structures. For models determined by joint refinement against X-ray and neutron data, the water structure comparison was also carried out for models re-refined against the X-ray data alone. The homologues have generally fewer conserved water molecules where X-ray only was used and the positions of equivalent waters vary more than in the case of the hybrid X-ray model. As neutron diffraction data are generally less complete than X-ray data, the influence of neutron data completeness on nuclear density maps was also analyzed. We observe and discuss systematic map quality deterioration as result of data incompleteness.

Implementation of the riding hydrogen model in CCTBX to support the next generation of X-ray and neutron joint refinement in Phenix

A fundamental prerequisite for implementing new procedures of atomic model refinement against neutron diffraction data is the efficient handling of hydrogen atoms. The riding hydrogen model, which constrains hydrogen atom parameters to those of the non-hydrogen atoms, is a plausible parameterization for refinements. This work describes the implementation of the riding hydrogen model in the Computational Crystallography Toolbox and in Phenix. Riding hydrogen atoms can be found in several different configurations that are characterized by specific geometries. For each configuration, the hydrogen atom parameterization and the expressions for the gradients of refinement target function with respect to non-hydrogen parameters are described.

X-ray crystallographic studies on the hydrogen isotope effects of green fluorescent protein at sub-ångström resolutions

Hydrogen atoms are critical to the nature and properties of proteins, and thus deuteration has the potential to influence protein function. In fact, it has been reported that some deuterated proteins show different physical and chemical properties to their protiated counterparts. Consequently, it is important to investigate protonation states around the active site when using deuterated proteins. Here, hydrogen isotope effects on the S65T/F99S/M153T/V163A variant of green fluorescent protein (GFP), in which the deprotonated B form is dominant at pH 8.5, were investigated. The pH/pD dependence of the absorption and fluorescence spectra indicates that the protonation state of the chromophore is the same in protiated GFP in H₂O and protiated GFP in D₂O at pH/pD 8.5, while the pK_a of the chromophore became higher in D₂O. Indeed, X-ray crystallographic analyses at sub-ångström resolution revealed no apparent changes in the protonation state of the chromophore between the two samples. However, detailed comparisons of the hydrogen OMIT maps revealed that the protonation state of His148 in the vicinity of the chromophore differed between the two samples. This indicates that protonation states around the active site should be carefully adjusted to be the same as those of the protiated protein when neutron crystallographic analyses of proteins are performed.

Improving the accuracy and resolution of neutron crystallographic data by three-dimensional profile fitting of Bragg peaks in reciprocal space

Neutron crystallography is a powerful technique for directly visualizing the locations of H atoms in biological macromolecules. This information has provided key new insights into enzyme mechanisms, ligand binding and hydration. However, despite the importance of this information, the application of neutron crystallography in biology has been limited by the relatively low flux of available neutron beams and the large incoherent neutron scattering from hydrogen, both of which contribute to weak diffraction data with relatively low signal-to-background ratios. A method has been developed to fit weak data based on three-dimensional profile fitting of Bragg peaks in reciprocal space by an Ikeda-Carpenter function with a bivariate Gaussian. When applied to data collected from three different proteins, three-dimensional profile fitting yields intensities with higher correlation coefficients (CC1/2) at high resolutions, decreased Rfree factors, extended resolutions and improved nuclear density maps. Importantly, additional features are revealed in nuclear density maps that may provide additional scientific information. These results suggest that three-dimensional profile fitting will help to extend the capabilities of neutron macromolecular crystallography.

Crystal structure of yeast xylose reductase in complex with a novel NADP-DTT adduct provides insights into substrate recognition and catalysis

Aldose reductases (ARs) belonging to the aldo-keto reductase (AKR) superfamily catalyze the conversion of carbonyl substrates into their respective alcohols. Here we report the crystal structures of the yeast Debaryomyces nepalensis xylose reductase (DnXR, AKR2B10) in the apo form and as a ternary complex with a novel NADP-DTT adduct. Xylose reductase, a key enzyme in the conversion of xylose to xylitol, has several industrial applications. The enzyme displayed the highest catalytic efficiency for l-threose (138 ± 7 mm^-1 ·s^-1 ) followed by d-erythrose (30 ± 3 mm^-1 ·s^-1 ). The crystal structure of the complex reveals a covalent linkage between the C4N atom of the nicotinamide ring of the cosubstrate and the S1 sulfur atom of DTT and provides the first structural evidence for a protein mediated NADP-low-molecular-mass thiol adduct. We hypothesize that the formation of the adduct is facilitated by an in-crystallo Michael addition of the DTT thiolate to the specific conformation of bound NADPH in the active site of DnXR. The interactions between DTT, a four-carbon sugar alcohol analog, and the enzyme are representative of a near-cognate product ternary complex and provide significant insights into the structural basis of aldose binding and specificity and the catalytic mechanism of ARs. DATABASE: Structural data are available in the PDB under the accession numbers 5ZCI and 5ZCM.

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.

New Insights into the Catalytic Mechanism of Aldose Reductase: A QM/MM Study

Aldose reductase is the first enzyme of the polyol pathway in which glucose is converted to fructose via sorbitol. The understanding of this key enzyme is important as it has been linked to some diabetes mellitus complications. The mechanism of the enzyme was investigated using a hybrid quantum mechanics/molecular mechanics (QM/MM) method. It was found that depending on the protonation state of His110 the mechanism can be concerted or stepwise and the proton donor can be either Tyr48 or His110. These findings are different from the previous theoretical studies based on QM/MM calculations using either AM1 or HF/4-31G, in which the reduction is, respectively, a stepwise or one-step process. The QM/MM energy barriers for the reduction of d-glyceraldehyde were evaluated at a B3LYP/6-31G* level for both HIP and HIE protonation states of His110. These were, respectively, 6.5 ± 2.2 and 16.7 ± 1.0 kcal/mol, which makes only the HIE protonation state consistent with the experimental value of 14.8 kcal/mol derived from kinetics experiments and makes Tyr48 the most probable proton donor.

Measurement of Enzyme Isotope Effects

Enzyme isotope effects, or the kinetic effects of "heavy" enzymes, refer to the effect of isotopically labeled protein residues on the enzyme's activity or physical properties. These effects are increasingly employed in the examination of the possible contributions of protein dynamics to enzyme catalysis. One hypothesis assumed that isotopic substitution of all ¹²C, ¹⁴N, and nonexchangeable ¹H by ¹³C, ¹⁵N, and ²H, would slow down protein picosecond to femtosecond dynamics without any effect on the system's electrostatics following the Born-Oppenheimer approximation. It was suggested that reduced reaction rates reported for several "heavy" enzymes accords with that hypothesis. However, numerous deviations from the predictions of that hypothesis were also reported. Current studies also attempt to test the role of individual residues by site-specific labeling or by labeling a pattern of residues on activity. It appears that in several systems the protein's fast dynamics are indeed reduced in "heavy" enzymes in a way that reduces the probability of barrier crossing of its chemical step. Other observations, however, indicated that slower protein dynamics are electrostatically altered in isotopically labeled enzymes. Interestingly, these effects appear to be system dependent, thus it might be premature to suggest a general role of "heavy" enzymes' effect on catalysis.

Characterization of erythrose reductase from Yarrowia lipolytica and its influence on erythritol synthesis

BACKGROUND

Erythritol is a natural sweetener that is used in the food industry. It is produced as an osmoprotectant by bacteria and yeast. Due to its chemical properties, it does not change the insulin level in the blood, and therefore it can be safely used by diabetics. Previously, it has been shown that erythrose reductase (ER), which catalyzes the final step, plays a crucial role in erythritol synthesis. ER reduces erythrose to erythritol with NAD(P)H as a cofactor. Despite many studies on erythritol synthesis by Yarrowia lipolytica, the enzymes involved in this metabolic pathway have ever been described.

RESULTS

The gene YALI0F18590g encoding the predicted erythrose reductase from Y. lipolytica was overexpressed, and its influence on erythritol synthesis was studied. The amino acid sequence of the Y. lipolytica ER showed a high degree of similarity to the previously described erythrose reductases from known erythritol producers, such as Candida magnoliae and Moniliella megachiliensis. Here, we found that the gene overexpression results in an enhanced titer of erythritol of 44.44 g/L (20% over the control), a yield of 0.44 g/g and productivity of 0.77 g/L/h. Moreover, on purification and characterization of the enzyme we found that it displays the highest activity at 37 °C and pH 3.0. The effects of various metal ions (Zn²⁺, Cu²⁺, Mn²⁺, Fe²⁺) on erythrose reductase were investigated. The addition of Zn²⁺ ions at 0.25 mM had a positive effect on the activity of erythrose reductase from Y. lipolytica, as well as on the erythritol production.

CONCLUSIONS

In this study we identified, overexpressed and characterized a native erythrose reductase in Y. lipolytica. Further optimizations of this strain via metabolic pathway engineering and media optimization strategies enabled 54 g/L to be produced in a shake-flask experiment. To date, this is the first reported study employing metabolic engineering of the native gene involved in the erythritol pathway to result in a high titer of the polyol. Moreover, it indicates the importance of environmental conditions for genetic targets in metabolic engineering.

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.

Structural Determinants of the Selectivity of 3-Benzyluracil-1-acetic Acids toward Human Enzymes Aldose Reductase and AKR1B10

The human enzymes aldose reductase (AR) and AKR1B10 have been thoroughly explored in terms of their roles in diabetes, inflammatory disorders, and cancer. In this study we identified two new lead compounds, 2-(3-(4-chloro-3-nitrobenzyl)-2,4-dioxo-3,4-dihydropyrimidin-1(2H)-yl)acetic acid (JF0048, 3) and 2-(2,4-dioxo-3-(2,3,4,5-tetrabromo-6-methoxybenzyl)-3,4-dihydropyrimidin-1(2H)-yl)acetic acid (JF0049, 4), which selectively target these enzymes. Although 3 and 4 share the 3-benzyluracil-1-acetic acid scaffold, they have different substituents in their aryl moieties. Inhibition studies along with thermodynamic and structural characterizations of both enzymes revealed that the chloronitrobenzyl moiety of compound 3 can open the AR specificity pocket but not that of the AKR1B10 cognate. In contrast, the larger atoms at the ortho and/or meta positions of compound 4 prevent the AR specificity pocket from opening due to steric hindrance and provide a tighter fit to the AKR1B10 inhibitor binding pocket, probably enhanced by the displacement of a disordered water molecule trapped in a hydrophobic subpocket, creating an enthalpic signature. Furthermore, this selectivity also occurs in the cell, which enables the development of a more efficient drug design strategy: compound 3 prevents sorbitol accumulation in human retinal ARPE-19 cells, whereas 4 stops proliferation in human lung cancer NCI-H460 cells.

The Macromolecular Neutron Diffractometer MaNDi at the Spallation Neutron Source

The Macromolecular Neutron Diffractometer (MaNDi) is located on beamline 11B of the Spallation Neutron Source at Oak Ridge National Laboratory. The instrument is a neutron time-of-flight wavelength-resolved Laue diffractometer optimized to collect diffraction data from single crystals. The instrument has been designed to provide flexibility in several instrumental parameters, such as beam divergence and wavelength bandwidth, to allow data collection from a range of macromolecular systems.

Deuterium Labeling Strategies for Creating Contrast in Structure-Function Studies of Model Bacterial Outer Membranes Using Neutron Reflectometry

Studying the outer membrane of Gram-negative bacteria is challenging due to the complex nature of its structure. Therefore, simplified models are required to undertake structure-function studies of processes that occur at the outer membrane/fluid interface. Model membranes can be created by immobilizing bilayers to solid supports such as gold or silicon surfaces, or as monolayers on a liquid support where the surface pressure and fluidity of the lipids can be controlled. Both model systems are amenable to having their structure probed by neutron reflectometry, a technique that provides a one-dimensional depth profile through a membrane detailing its thickness and composition. One of the strengths of neutron scattering is the ability to use contrast matching, allowing molecules containing hydrogen and those enriched with deuterium to be highlighted or matched out against the bulk isotopic composition of the solvent. Lipopolysaccharides, a major component of the outer membrane, can be isolated for incorporation into model membranes. Here, we describe the deuteration of lipopolysaccharides from rough strains of Escherichia coli for incorporation into model outer membranes, and how the use of deuterated materials enhances structural analysis of model membranes by neutron reflectometry.

Production and characterization of recombinant perdeuterated cholesterol oxidase

Cholesterol oxidase (CO) is a FAD (flavin adenine dinucleotide) containing enzyme that catalyzes the oxidization and isomerization of cholesterol. Studies directed toward elucidating the catalytic mechanism of CO will provide an important general understanding of Flavin-assisted redox catalysis. Hydrogen atoms play an important role in enzyme catalysis; however, they are not readily visualized in protein X-ray diffraction structures. Neutron crystallography is an ideal method for directly visualizing hydrogen positions at moderate resolutions because hydrogen and deuterium have comparable neutron scattering lengths to other heavy atoms present in proteins. The negative coherent and large incoherent scattering lengths of hydrogen atoms in neutron diffraction experiments can be circumvented by replacing hydrogen atoms with its isotope, deuterium. The perdeuterated form of CO was successfully expressed from minimal medium, purified, and crystallized. X-ray crystallographic structures of the enzyme in the perdeuterated and hydrogenated states confirm that there are no apparent structural differences between the two enzyme forms. Kinetic assays demonstrate that perdeuterated and hydrogenated enzymes are functionally identical. Together, structural and functional studies indicate that the perdeuterated protein is suitable for structural studies by neutron crystallography directed at understanding the role of hydrogen atoms in enzyme catalysis.

Perdeuteration: improved visualization of solvent structure in neutron macromolecular crystallography

The 1.8 Å resolution neutron structure of deuterated type III antifreeze protein in which the methyl groups of leucine and valine residues are selectively protonated is presented. Comparison between this and the 1.85 Å resolution neutron structure of perdeuterated type III antifreeze protein indicates that perdeuteration improves the visibility of solvent molecules located in close vicinity to hydrophobic residues, as cancellation effects between H atoms of the methyl groups and nearby heavy-water molecules (D2O) are avoided.

B-factor Analysis and Conformational Rearrangement of Aldose Reductase

The NADPH-dependent reduction of glucose reaction that is catalyzed by Aldose Reductase (AR) follows a sequential ordered kinetic mechanism in which the co-factor NADPH binds to the enzyme prior to the aldehyde substrate. The kinetic/structural experiments have found a conformational change involving a hinge-like movement of a surface loop (residues 213-224) which is anticipated to take place upon the binding of the diphosphate moiety of NADPH. The reorientation of this loop, expected to permit the release of NADP⁺, represents the rate-limiting step of the catalytic mechanism. This study reveals: 1) The Translation/Libration/Screw (TLS) analysis of absolute B-factors of apo AR crystal structures indicates that the 212-224 loop might move as a rigid group. 2) Residues that make the flexible loop slide in the AR binary and ternary complexes. 3) The normalized B-factors separate this segment into three different clusters with fewer residues.

Induced fit docking, pharmacophore modeling, and molecular dynamic simulations on thiazolidinedione derivatives to explore key interactions with Tyr48 in polyol pathway

To obtain a scientific thought and expedition to explore key interactions with Tyr48 in aldose reductase (ALR), combined study of pharmacophore modeling, induced fit docking, and dynamics studies were performed on ALR. A stereo chemically and energetically valid model of ALR-NADP+ complex was developed using homology modeling technique. Statistically a significant five point pharmacophore model was designed on a set of 54 thiazolidinedione derivatives with good external and internal predictive ability. Rigid and induced fit docking protocols were applied on ALR protein for both with and without NADP+ cofactor to identify a suitable binding mode that facilitates the key hydrogen bond interactions with Tyr48. Docking of thiazolidinedione derivatives into ALR-NADP+ complex gave more promising results by reducing false positive binding of inhibitors into the co-factor binding site. Structural changes within Try48 and Asp43 during the binding process in enzyme inhibitor complex were analyzed using molecular dynamics (MD) simulations. The results obtained from dynamic simulations emphasized the role of Tyr48 in maintaining inter or intra molecular hydrogen bond interaction with the protein or inhibitor respectively. New molecules were designed and checked for their binding interactions and showed improved results compared to existing thiazolidinediones derivatives. Hence, these combined protocols will be helpful and cooperative to design and optimize molecules with better inhibitory activity against the biologically active target.

Looking for Hydrogen Atoms: Neutron Crystallography Provides Novel Insights Into Protein Structure and Function

Neutron crystallography allows direct localization of hydrogen positions in biological macromolecules. Within enzymes, hydrogen atoms play a pivotal role in catalysis. Recent advances in instrumentation and sample preparation have helped to overcome the difficulties of performing neutron diffraction experiments on protein crystals. The application of neutron macromolecular crystallography to a growing number of proteins has yielded novel structural insights. The ability to accurately position water molecules, hydronium ions, and hydrogen atoms within protein structures has helped in the study of low-barrier hydrogen bonds and hydrogen-bonding networks. The determination of protonation states of protein side chains, substrates, and inhibitors in the context of the macromolecule has provided important insights into enzyme chemistry and ligand binding affinities, which can assist in the design of potent therapeutic agents. In this review, we give an overview of the method and highlight advances in knowledge attained through the application of neutron protein crystallography.

Modulation of aldose reductase inhibition by halogen bond tuning

In this paper, we studied a designed series of aldose reductase (AR) inhibitors. The series was derived from a known AR binder, which had previously been shown to form a halogen bond between its bromine atom and the oxygen atom of the Thr-113 side chain of AR. In the series, the strength of the halogen bond was modulated by two factors, namely bromine-iodine substitution and the fluorination of the aromatic ring in several positions. The role of the single halogen bond in AR-ligand binding was elucidated by advanced binding free energy calculations involving the semiempirical quantum chemical Hamiltonian. The results were complemented with ultrahigh-resolution X-ray crystallography and IC50 measurements. All of the AR inhibitors studied were shown by X-ray crystallography to bind in an identical manner. Further, it was demonstrated that it was possible to decrease the IC50 value by about 1 order of magnitude by tuning the strength of the halogen bond by a monoatomic substitution. The calculations revealed that the protein-ligand interaction energy increased upon the substitution of iodine for bromine or upon the addition of electron-withdrawing fluorine atoms to the ring. However, the effect on the binding affinity was found to be more complex due to the change of the solvation/desolvation properties within the ligand series. The study shows that it is possible to modulate the strength of a halogen bond in a protein-ligand complex as was designed based on the previous studies of low-molecular-weight complexes.

The IMAGINE instrument: first neutron protein structure and new capabilities for neutron macromolecular crystallography

The first high-resolution neutron protein structure of perdeuterated rubredoxin fromPyrococcus furiosus(PfRd) determined using the new IMAGINE macromolecular neutron crystallography instrument at the Oak Ridge National Laboratory is reported. Neutron diffraction data extending to 1.65 Å resolution were collected from a relatively small 0.7 mm3PfRd crystal using 2.5 d (60 h) of beam time. The refined structure contains 371 out of 391, or 95%, of the D atoms of the protein and 58 solvent molecules. The IMAGINE instrument is designed to provide neutron data at or near atomic resolution (1.5 Å) from crystals with volume <1.0 mm3and with unit-cell edges <100 Å. Beamline features include novel elliptical focusing mirrors that deliver neutrons into a 2.0 × 3.2 mm focal spot at the sample position with full-width vertical and horizontal divergences of 0.5 and 0.6°, respectively. Variable short- and long-wavelength cutoff optics provide automated exchange between multiple-wavelength configurations (λmin= 2.0, 2.8, 3.3 Å to λmax= 3.0, 4.0, 4.5, ∼20 Å). These optics produce a more than 20-fold increase in the flux density at the sample and should help to enable more routine collection of high-resolution data from submillimetre-cubed crystals. Notably, the crystal used to collect thesePfRd data was 5–10 times smaller than those previously reported.

Neutron diffraction studies towards deciphering the protonation state of catalytic residues in the bacterial KDN9P phosphatase

The enzyme 2-keto-3-deoxy-9-O-phosphonononic acid phosphatase (KDN9P phosphatase) functions in the pathway for the production of 2-keto-3-deoxy-D-glycero-D-galacto-nononic acid, a sialic acid that is important for the survival of commensal bacteria in the human intestine. The enzyme is a member of the haloalkanoate dehalogenase superfamily and represents a good model for the active-site protonation state of family members. Crystals of approximate dimensions 1.5 × 1.0 × 1.0 mm were obtained in space group P2(1)2(1)2, with unit-cell parameters a = 83.1, b = 108.9, c = 75.7 Å. A complete neutron data set was collected from a medium-sized H/D-exchanged crystal at BIODIFF at the Heinz Maier-Leibnitz Zentrum (MLZ), Garching, Germany in 18 d. Initial refinement to 2.3 Å resolution using only neutron data showed significant density for catalytically important residues.

Neutron scattering techniques and applications in structural biology

Neutron scattering is exquisitely sensitive to the position, concentration, and dynamics of hydrogen atoms in materials and is a powerful tool for the characterization of structure-function and interfacial relationships in biological systems. Modern neutron scattering facilities offer access to a sophisticated, nondestructive suite of instruments for biophysical characterization that provides spatial and dynamic information spanning from Ångstroms to microns and from picoseconds to microseconds, respectively. Applications in structural biology range from the atomic-resolution analysis of individual hydrogen atoms in enzymes through to meso- and macro-scale analysis of complex biological structures, membranes, and assemblies. The large difference in neutron scattering length between hydrogen and deuterium allows contrast variation experiments to be performed and enables H/D isotopic labeling to be used for selective and systematic analysis of the local structure, dynamics, and interactions of multi-component systems. This overview describes the available techniques and summarizes their practical application to the study of biomolecular systems.

Joint X-ray/neutron crystallographic study of HIV-1 protease with clinical inhibitor amprenavir: insights for drug design

HIV-1 protease is an important target for the development of antiviral inhibitors to treat AIDS. A room-temperature joint X-ray/neutron structure of the protease in complex with clinical drug amprenavir has been determined at 2.0 Å resolution. The structure provides direct determination of hydrogen atom positions in the enzyme active site. Analysis of the enzyme-drug interactions suggests that some hydrogen bonds may be weaker than deduced from the non-hydrogen interatomic distances. This information may be valuable for the design of improved protease inhibitors.

NDM-1, the ultimate promiscuous enzyme: substrate recognition and catalytic mechanism

The specter of a return to an era in which infectious disease looms as a significant threat to human health is not just hyperbole; there are serious concerns about the widespread overuse and misuse of antibiotics contributing to increased antibiotic resistance in pathogens. The recent discovery of a new enzyme, first identified in Klebsiella pneumoniae from a patient from New Delhi and denoted as NDM-1, represents an example of extreme promiscuity: It hydrolyzes and inactivates nearly all known β-lactam-based antibiotics with startling efficiency. NDM-1 can utilize different metal cofactors and seems to exploit an alternative mechanism based on the reaction conditions. Here we report the results of a combined experimental and theoretical study that examines the substrate, metal binding, and catalytic mechanism of the enzyme. We utilize structures obtained through X-ray crystallography, biochemical assays, and numerical simulation to construct a model of the enzyme catalytic pathway. The NDM-1 enzyme interacts with the substrate solely through zinc, or other metals, bound in the active site, explaining the observed lack of specificity against a broad range of β-lactam antibiotic agents. The zinc ions also serve to activate a water molecule that hydrolyzes the β-lactam ring through a proton shuttle.

Solution conformations of early intermediates in Mos1 transposition

DNA transposases facilitate genome rearrangements by moving DNA transposons around and between genomes by a cut-and-paste mechanism. DNA transposition proceeds in an ordered series of nucleoprotein complexes that coordinate pairing and cleavage of the transposon ends and integration of the cleaved ends at a new genomic site. Transposition is initiated by transposase recognition of the inverted repeat sequences marking each transposon end. Using a combination of solution scattering and biochemical techniques, we have determined the solution conformations and stoichiometries of DNA-free Mos1 transposase and of the transposase bound to a single transposon end. We show that Mos1 transposase is an elongated homodimer in the absence of DNA and that the N-terminal 55 residues, containing the first helix-turn-helix motif, are required for dimerization. This arrangement is remarkably different from the compact, crossed architecture of the dimer in the Mos1 paired-end complex (PEC). The transposase remains elongated when bound to a single-transposon end in a pre-cleavage complex, and the DNA is bound predominantly to one transposase monomer. We propose that a conformational change in the single-end complex, involving rotation of one half of the transposase along with binding of a second transposon end, could facilitate PEC assembly.

Neutron and X-ray crystal structures of a perdeuterated enzyme inhibitor complex reveal the catalytic proton network of the Toho-1 β-lactamase for the acylation reaction

The mechanism by which class A β-lactamases hydrolyze β-lactam antibiotics has been the subject of intensive investigation using many different experimental techniques. Here, we report on the novel use of both neutron and high resolution x-ray diffraction to help elucidate the identity of the catalytic base in the acylation part of the catalytic cycle, wherein the β-lactam ring is opened and an acyl-enzyme intermediate forms. To generate protein crystals optimized for neutron diffraction, we produced a perdeuterated form of the Toho-1 β-lactamase R274N/R276N mutant. Protein perdeuteration, which involves replacing all of the hydrogen atoms in a protein with deuterium, gives a much stronger signal in neutron diffraction and enables the positions of individual deuterium atoms to be located. We also synthesized a perdeuterated acylation transition state analog, benzothiophene-2-boronic acid, which was also isotopically enriched with (11)B, as (10)B is a known neutron absorber. Using the neutron diffraction data from the perdeuterated enzyme-inhibitor complex, we were able to determine the positions of deuterium atoms in the active site directly rather than by inference. The neutron diffraction results, along with supporting bond-length analysis from high resolution x-ray diffraction, strongly suggest that Glu-166 acts as the general base during the acylation reaction.

Protonation-state determination in proteins using high-resolution X-ray crystallography: effects of resolution and completeness

A bond-distance analysis has been undertaken to determine the protonation states of ionizable amino acids in trypsin, subtilisin and lysozyme. The diffraction resolutions were 1.2 Å for trypsin (97% complete, 12% H-atom visibility at 2.5σ), 1.26 Å for subtilisin (100% complete, 11% H-atom visibility at 2.5σ) and 0.65 Å for lysozyme (PDB entry 2vb1; 98% complete, 30% H-atom visibility at 3σ). These studies provide a wide diffraction resolution range for assessment. The bond-length e.s.d.s obtained are as small as 0.008 Å and thus provide an exceptional opportunity for bond-length analyses. The results indicate that useful information can be obtained from diffraction data at around 1.2-1.3 Å resolution and that minor increases in resolution can have significant effects on reducing the associated bond-length standard deviations. The protonation states in histidine residues were also considered; however, owing to the smaller differences between the protonated and deprotonated forms it is much more difficult to infer the protonation states of these residues. Not even the 0.65 Å resolution lysozyme structure provided the necessary accuracy to determine the protonation states of histidine.

Room-temperature ultrahigh-resolution time-of-flight neutron and X-ray diffraction studies of H/D-exchanged crambin

The room-temperature (RT) X-ray structure of H/D-exchanged crambin is reported at 0.85 Å resolution. As one of the very few proteins refined with anisotropic atomic displacement parameters at two temperatures, the dynamics of atoms in the RT and 100 K structures are compared. Neutron diffraction data from an H/D-exchanged crambin crystal collected at the Protein Crystallography Station (PCS) showed diffraction beyond 1.1 Å resolution. This is the highest resolution neutron diffraction reported to date for a protein crystal and will reveal important details of the anisotropic motions of H and D atoms in protein structures.

Rapid visualization of hydrogen positions in protein neutron crystallographic structures

Neutron crystallography is a powerful technique for experimental visualization of the positions of light atoms, including hydrogen and its isotope deuterium. In recent years, structural biologists have shown increasing interest in the technique as it uniquely complements X-ray crystallographic data by revealing the positions of D atoms in macromolecules. With this regained interest, access to macromolecular neutron crystallography beamlines is becoming a limiting step. In this report, it is shown that a rapid data-collection strategy can be a valuable alternative to longer data-collection times in appropriate cases. Comparison of perdeuterated rubredoxin structures refined against neutron data sets collected over hours and up to 5 d shows that rapid neutron data collection in just 14 h is sufficient to provide the positions of 269 D atoms without ambiguity.

Preliminary neutron crystallographic study of human transthyretin

Preliminary studies of perdeuterated crystals of human transthyretin (TTR) have been carried out using the LADI-III and D19 diffractometers at the Institut Laue-Langevin in Grenoble. The results demonstrate the feasibility of a full crystallographic analysis to a resolution of 2.0 Å using Laue diffraction and also illustrate the potential of using monochromatic instruments such as D19 for higher resolution studies where larger crystals having smaller unit cells are available. This study will yield important information on hydrogen bonding, amino-acid protonation states and hydration in the protein. Such information will be of general interest for an understanding of the factors that stabilize/destabilize TTR and for the design of ligands that may be used to counter TTR amyloid fibrillogenesis.

Time-of-flight neutron diffraction study of bovine γ-chymotrypsin at the Protein Crystallography Station

The overarching goal of this research project is to determine, for a subset of proteins, exact hydrogen positions using neutron diffraction, thereby improving H-atom placement in proteins so that they may be better used in various computational methods that are critically dependent upon said placement. In order to be considered applicable for neutron diffraction studies, the protein of choice must be amenable to ultrahigh-resolution X-ray crystallography, be able to form large crystals (1 mm(3) or greater) and have a modestly sized unit cell (no dimension longer than 100 Å). As such, γ-chymotrypsin is a perfect candidate for neutron diffraction. To understand and probe the role of specific active-site residues and hydrogen-bonding patterns in γ-chymotrypsin, neutron diffraction studies were initiated at the Protein Crystallography Station (PCS) at Los Alamos Neutron Science Center (LANSCE). A large single crystal was subjected to H/D exchange prior to data collection. Time-of-flight neutron diffraction data were collected to 2.0 Å resolution at the PCS with ~85% completeness. Here, the first time-of-flight neutron data collection from γ-chymotrypsin is reported.

Neutron fibre diffraction studies of amyloid using H2O/D2O isotopic replacement

The first neutron fibre diffraction studies of an amyloid system are presented. The techniques used to prepare the large samples needed are described, as well as the procedures used to isotopically replace H2O in the sample by D2O. The results demonstrate the feasibility of this type of approach for the pursuit of novel structural analyses that will strongly complement X-ray fibre diffraction studies and probe aspects of amyloid structure that to date have remained obscure. The approach is demonstrated using an amyloid form of the peptide NSGAITIG, but is equally applicable for the study of other systems such as Alzheimer's Aβ peptide.

Preliminary joint X-ray and neutron protein crystallographic studies of endoxylanase II from the fungus Trichoderma longibrachiatum

Room-temperature X-ray and neutron diffraction data were measured from a family 11 endoxylanase holoenzyme (XynII) originating from the filamentous fungus Trichoderma longibrachiatum to 1.55 Å resolution using a home source and to 1.80 Å resolution using the Protein Crystallography Station at LANSCE. Crystals of XynII, which is an important enzyme for biofuel production, were grown at pH 8.5 in order to examine the effect of basic conditions on the protonation-state distribution in the active site and throughout the protein molecule and to provide insights for rational engineering of catalytically improved XynII for industrial applications.

Preliminary joint neutron time-of-flight and X-ray crystallographic study of human ABO(H) blood group A glycosyltransferase

The biosyntheses of oligosaccharides and glycoconjugates are conducted by glycosyltransferases. These extraordinarily diverse and widespread enzymes catalyze the formation of glycosidic bonds through the transfer of a monosaccharide from a donor molecule to an acceptor molecule, with the stereochemistry about the anomeric carbon being either inverted or retained. Human ABO(H) blood group A α-1,3-N-acetylgalactosaminyltransferase (GTA) generates the corresponding antigen by the transfer of N-acetylgalactosamine from UDP-GalNAc to the blood group H antigen. To understand better how specific active-site-residue protons and hydrogen-bonding patterns affect substrate recognition and catalysis, neutron diffraction studies were initiated at the Protein Crystallography Station (PCS) at Los Alamos Neutron Science Center (LANSCE). A large single crystal was subjected to H/D exchange prior to data collection and time-of-flight neutron diffraction data were collected to 2.5 Å resolution at the PCS to ∼85% overall completeness, with complementary X-ray diffraction data collected from a crystal from the same drop and extending to 1.85 Å resolution. Here, the first successful neutron data collection from a glycosyltransferase is reported.

Neutron structure and mechanistic studies of diisopropyl fluorophosphatase (DFPase)

Diisopropyl fluorophosphatase (DFPase) is a calcium-dependent phosphotriesterase that acts on a variety of highly toxic organophosphorus compounds that act as inhibitors of acetylcholinesterase. The mechanism of DFPase has been probed using a variety of methods, including isotopic labelling, which demonstrated the presence of a phosphoenzyme intermediate in the reaction mechanism. In order to further elucidate the mechanism of DFPase and to ascertain the protonation states of the residues and solvent molecules in the active site, the neutron structure of DFPase was solved at 2.2 Å resolution. The proposed nucleophile Asp229 is deprotonated, while the active-site solvent molecule W33 was identified as water and not hydroxide. These data support a mechanism involving direct nucleophilic attack by Asp229 on the substrate and rule out a mechanism involving metal-assisted water activation. These data also allowed for the re-engineering of DFPase through rational design to bind and productively orient the more toxic S(P) stereoisomers of the nerve agents sarin and cyclosarin, creating a modified enzyme with enhanced overall activity and significantly increased detoxification properties.