Neutron diffraction from a microgravity-grown crystal reveals the active site hydrogens of the internal aldimine form of tryptophan synthase

Pyridoxal 5'-phosphate (PLP), the biologically active form of vitamin B6, is an essential cofactor in many biosynthetic pathways. The emergence of PLP-dependent enzymes as drug targets and biocatalysts, such as tryptophan synthase (TS), has underlined the demand to understand PLP-dependent catalysis and reaction specificity. The ability of neutron diffraction to resolve the positions of hydrogen atoms makes it an ideal technique to understand how the electrostatic environment and selective protonation of PLP regulates PLP-dependent activities. Facilitated by microgravity crystallization of TS with the Toledo Crystallization Box, we report the 2.1 Å joint X-ray/neutron (XN) structure of TS with PLP in the internal aldimine form. Positions of hydrogens were directly determined in both the α- and β-active sites, including PLP cofactor. The joint XN structure thus provides insight into the selective protonation of the internal aldimine and the electrostatic environment of TS necessary to understand the overall catalytic mechanism.

Discovery of Novel Rhizoctonia solani DHFR Inhibitors as Fungicides Using Virtual Screening

Dihydrofolate reductase (DHFR) is an essential enzyme in the folate pathway and has been recognized as a well-known target for antibacterial and antifungal drugs. We discovered eight compounds from the ZINC database using virtual screening to inhibit Rhizoctonia solani (R. solani), a fungal pathogen in crops. These compounds were evaluated with in vitro assays for enzymatic and antifungal activity. Among these, compound Hit8 is the most active R. solani DHFR inhibitor, with the IC50 of 10.2 μM. The selectivity of inhibition is 22.3 against human DHFR with the IC50 of 227.7 μM. Moreover, Hit8 has higher antifungal activity against R. solani (EC50 of 38.2 mg L-1) compared with validamycin A (EC50 of 67.6 mg L-1), a well-documented fungicide. These results suggest that Hit8 may be a potential fungicide. Our study exemplifies a computer-aided method to discover novel inhibitors that could target plant pathogenic fungi.

Insights into the mechanism of SARS-CoV-2 main protease autocatalytic maturation from model precursors

A critical step for SARS-CoV-2 assembly and maturation involves the autoactivation of the main protease (MProWT) from precursor polyproteins. Upon expression, a model precursor of MProWT mediates its own release at its termini rapidly to yield a mature dimer. A construct with an E290A mutation within MPro exhibits time dependent autoprocessing of the accumulated precursor at the N-terminal nsp4/nsp5 site followed by the C-terminal nsp5/nsp6 cleavage. In contrast, a precursor containing E290A and R298A mutations (MProM) displays cleavage only at the nsp4/nsp5 site to yield an intermediate monomeric product, which is cleaved at the nsp5/nsp6 site only by MProWT. MProM and the catalytic domain (MPro1-199) fused to the truncated nsp4 region also show time-dependent conversion in vitro to produce MProM and MPro1-199, respectively. The reactions follow first-order kinetics indicating that the nsp4/nsp5 cleavage occurs via an intramolecular mechanism. These results support a mechanism involving an N-terminal intramolecular cleavage leading to an increase in the dimer population and followed by an intermolecular cleavage at the C-terminus. Thus, targeting the predominantly monomeric MPro precursor for inhibition may lead to the identification of potent drugs for treatment.

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.

Contribution of the catalytic dyad of SARS-CoV-2 main protease to binding covalent and noncovalent inhibitors

The effect of mutations of the catalytic dyad residues of SARS-CoV-2 main protease (MPro^WT) on the thermodynamics of binding of covalent inhibitors comprising nitrile [nirmatrelvir (NMV), NBH2], aldehyde (GC373) and ketone (BBH1) warheads to MPro is examined together with room temperature X-ray crystallography. When lacking the nucleophilic C145, NMV binding is ∼400-fold weaker corresponding to 3.5 kcal/mol and 13.3 °C decreases in free energy (ΔG) and thermal stability (T_m), respectively, relative to MPro^WT. The H41A mutation results in a 20-fold increase in the dissociation constant (K_d), and 1.7 kcal/mol and 1.4 °C decreases in ΔG and T_m, respectively. Increasing the pH from 7.2 to 8.2 enhances NMV binding to MPro^H41A, whereas no significant change is observed in binding to MPro^WT. Structures of the four inhibitor complexes with MPro^1-304/C145A show that the active site geometries of the complexes are nearly identical to that of MPro^WT with the nucleophilic sulfur of C145 positioned to react with the nitrile or the carbonyl carbon. These results support a two-step mechanism for the formation of the covalent complex involving an initial non-covalent binding followed by a nucleophilic attack by the thiolate anion of C145 on the warhead carbon. Noncovalent inhibitor ensitrelvir (ESV) exhibits a binding affinity to MPro^WT that is similar to NMV but differs in its thermodynamic signature from NMV. The binding of ESV to MPro^C145A also results in a significant, but smaller, increase in K_d and decrease in ΔG and T_m, relative to NMV.

AI-Accelerated Design of Targeted Covalent Inhibitors for SARS-CoV-2

Direct-acting antivirals for the treatment of the COVID-19 pandemic caused by the SARS-CoV-2 virus are needed to complement vaccination efforts. Given the ongoing emergence of new variants, automated experimentation, and active learning based fast workflows for antiviral lead discovery remain critical to our ability to address the pandemic's evolution in a timely manner. While several such pipelines have been introduced to discover candidates with noncovalent interactions with the main protease (Mpro), here we developed a closed-loop artificial intelligence pipeline to design electrophilic warhead-based covalent candidates. This work introduces a deep learning-assisted automated computational workflow to introduce linkers and an electrophilic "warhead" to design covalent candidates and incorporates cutting-edge experimental techniques for validation. Using this process, promising candidates in the library were screened, and several potential hits were identified and tested experimentally using native mass spectrometry and fluorescence resonance energy transfer (FRET)-based screening assays. We identified four chloroacetamide-based covalent inhibitors of Mpro with micromolar affinities (KI of 5.27 μM) using our pipeline. Experimentally resolved binding modes for each compound were determined using room-temperature X-ray crystallography, which is consistent with the predicted poses. The induced conformational changes based on molecular dynamics simulations further suggest that the dynamics may be an important factor to further improve selectivity, thereby effectively lowering KI and reducing toxicity. These results demonstrate the utility of our modular and data-driven approach for potent and selective covalent inhibitor discovery and provide a platform to apply it to other emerging targets.

Structural and dynamic effects of paraoxon binding to human acetylcholinesterase by X-ray crystallography and inelastic neutron scattering

Organophosphorus (OP) compounds, including nerve agents and some pesticides, covalently bind to the catalytic serine of human acetylcholinesterase (hAChE), thereby inhibiting acetylcholine hydrolysis necessary for efficient neurotransmission. Oxime antidotes can reactivate the OP-conjugated hAChE, but reactivation efficiency can be low for pesticides, such as paraoxon (POX). Understanding structural and dynamic determinants of OP inhibition and reactivation can provide insights to design improved reactivators. Here, X-ray structures of hAChE with unaged POX, with POX and oximes MMB4 and RS170B, and with MMB4 are reported. A significant conformational distortion of the acyl loop was observed upon POX binding, being partially restored to the native conformation by oximes. Neutron vibrational spectroscopy combined with molecular dynamics simulations showed that picosecond vibrational dynamics of the acyl loop soften in the ∼20-50 cm-1 frequency range. The acyl loop structural perturbations may be correlated with its picosecond vibrational dynamics to yield more comprehensive template for structure-based reactivator design.

Autoprocessing and oxyanion loop reorganization upon GC373 and nirmatrelvir binding of monomeric SARS-CoV-2 main protease catalytic domain

The monomeric catalytic domain (residues 1–199) of SARS-CoV-2 main protease (MPro^1-199) fused to 25 amino acids of its flanking nsp4 region mediates its autoprocessing at the nsp4-MPro^1-199 junction. We report the catalytic activity and the dissociation constants of MPro^1-199 and its analogs with the covalent inhibitors GC373 and nirmatrelvir (NMV), and the estimated monomer-dimer equilibrium constants of these complexes. Mass spectrometry indicates the presence of the accumulated adduct of NMV bound to MPro^WT and MPro^1-199 and not of GC373. A room temperature crystal structure reveals a native-like fold of the catalytic domain with an unwound oxyanion loop (E state). In contrast, the structure of a covalent complex of the catalytic domain-GC373 or NMV shows an oxyanion loop conformation (E* state) resembling the full-length mature dimer. These results suggest that the E-E* equilibrium modulates autoprocessing of the main protease when converting from a monomeric polyprotein precursor to the mature dimer.

Structural characterization and catalytic activity of SARS-CoV-2 main protease reveal minimal interface regions enabling dimer formation driven by inhibitor-induced conformational changes of the oxyanion loop.

Structural and functional characterization of NEMO cleavage by SARS-CoV-2 3CLpro

In addition to its essential role in viral polyprotein processing, the SARS-CoV-2 3C-like protease (3CLpro) can cleave human immune signaling proteins, like NF-κB Essential Modulator (NEMO) and deregulate the host immune response. Here, in vitro assays show that SARS-CoV-2 3CLpro cleaves NEMO with fine-tuned efficiency. Analysis of the 2.50 Å resolution crystal structure of 3CLpro C145S bound to NEMO226-234 reveals subsites that tolerate a range of viral and host substrates through main chain hydrogen bonds while also enforcing specificity using side chain hydrogen bonds and hydrophobic contacts. Machine learning- and physics-based computational methods predict that variation in key binding residues of 3CLpro-NEMO helps explain the high fitness of SARS-CoV-2 in humans. We posit that cleavage of NEMO is an important piece of information to be accounted for, in the pathology of COVID-19.

EWALD: A macromolecular diffractometer for the second target station

Revealing the positions of all the atoms in large macromolecules is powerful but only possible with neutron macromolecular crystallography (NMC). Neutrons provide a sensitive and gentle probe for the direct detection of protonation states at near-physiological temperatures and clean of artifacts caused by x rays or electrons. Currently, NMC use is restricted by the requirement for large crystal volumes even at state-of-the-art instruments such as the macromolecular neutron diffractometer at the Spallation Neutron Source. EWALD's design will break the crystal volume barrier and, thus, open the door for new types of experiments, the study of grand challenge systems, and the more routine use of NMC in biology. EWALD is a single crystal diffractometer capable of collecting data from macromolecular crystals on orders of magnitude smaller than what is currently feasible. The construction of EWALD at the Second Target Station will cause a revolution in NMC by enabling key discoveries in the biological, biomedical, and bioenergy sciences.

The mechanisms of catalysis and ligand binding for the SARS-CoV-2 NSP3 macrodomain from neutron and x-ray diffraction at room temperature

The nonstructural protein 3 (NSP3) macrodomain of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (Mac1) removes adenosine diphosphate (ADP) ribosylation posttranslational modifications, playing a key role in the immune evasion capabilities of the virus responsible for the coronavirus disease 2019 pandemic. Here, we determined neutron and x-ray crystal structures of the SARS-CoV-2 NSP3 macrodomain using multiple crystal forms, temperatures, and pHs, across the apo and ADP-ribose-bound states. We characterize extensive solvation in the Mac1 active site and visualize how water networks reorganize upon binding of ADP-ribose and non-native ligands, inspiring strategies for displacing waters to increase the potency of Mac1 inhibitors. Determining the precise orientations of active site water molecules and the protonation states of key catalytic site residues by neutron crystallography suggests a catalytic mechanism for coronavirus macrodomains distinct from the substrate-assisted mechanism proposed for human MacroD2. These data provoke a reevaluation of macrodomain catalytic mechanisms and will guide the optimization of Mac1 inhibitors.

Hit Expansion of a Noncovalent SARS-CoV-2 Main Protease Inhibitor

Inhibition of the SARS-CoV-2 main protease (M^pro) is a major focus of drug discovery efforts against COVID-19. Here we report a hit expansion of non-covalent inhibitors of M^pro. Starting from a recently discovered scaffold (The COVID Moonshot Consortium. Open Science Discovery of Oral Non-Covalent SARS-CoV-2 Main Protease Inhibitor Therapeutics. bioRxiv 2020.10.29.339317) represented by an isoquinoline series, we searched a database of over a billion compounds using a cheminformatics molecular fingerprinting approach. We identified and tested 48 compounds in enzyme inhibition assays, of which 21 exhibited inhibitory activity above 50% at 20 μM. Among these, four compounds with IC₅₀ values around 1 μM were found. Interestingly, despite the large search space, the isoquinolone motif was conserved in each of these four strongest binders. Room-temperature X-ray structures of co-crystallized protein–inhibitor complexes were determined up to 1.9 Å resolution for two of these compounds as well as one of the stronger inhibitors in the original isoquinoline series, revealing essential interactions with the binding site and water molecules. Molecular dynamics simulations and quantum chemical calculations further elucidate the binding interactions as well as electrostatic effects on ligand binding. The results help explain the strength of this new non-covalent scaffold for M^pro inhibition and inform lead optimization efforts for this series, while demonstrating the effectiveness of a high-throughput computational approach to expanding a pharmacophore library.

Covalent narlaprevir- and boceprevir-derived hybrid inhibitors of SARS-CoV-2 main protease: room-temperature X-ray and neutron crystallography, binding thermodynamics, and antiviral activity.. [PMID: 35169792 PMCID: PMC8845512 DOI: 10.21203/rs.3.rs-1318037/v1]

The COVID-19 pandemic continues to disrupt everyday life, with constantly emerging SARS-CoV-2 variants threatening to render current vaccines ineffective. Small-molecule antivirals can provide an important therapeutic treatment option that is subject to challenges caused by the virus variants. The viral main protease (M^pro) is critical for the virus replication and thus is considered an attractive drug target for specific protease inhibitors. We performed the design and characterization of three reversible covalent hybrid inhibitors BBH-1, BBH-2 and NBH-2, whose structures were derived from those of hepatitis C protease inhibitors boceprevir and narlaprevir. A joint X-ray/neutron structure of the M^pro/BBH-1 complex demonstrated that a Cys145 thiolate reaction with the inhibitor’s keto-warhead creates a negatively charged oxyanion, similar to that proposed for the M^pro-catalyzed peptide bond hydrolysis. Protonation states of the ionizable residues in the M^pro active site adapt to the inhibitor, which appears to be an intrinsic property of M^pro. Structural comparisons of the hybrid inhibitors with PF-07321332 revealed unconventional interactions of PF-07321332 with M^pro which may explain its more favorable enthalpy of binding and consequently higher potency. BBH-1, BBH-2 and NBH-2 demonstrated comparable antiviral properties in vitro relative to PF-07321332, making them good candidates for further design of improved antivirals.

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.

The mechanisms of catalysis and ligand binding for the SARS-CoV-2 NSP3 macrodomain from neutron and X-ray diffraction at room temperature

The NSP3 macrodomain of SARS CoV 2 (Mac1) removes ADP-ribosylation post-translational modifications, playing a key role in the immune evasion capabilities of the virus responsible for the COVID-19 pandemic. Here, we determined neutron and X-ray crystal structures of the SARS-CoV-2 NSP3 macrodomain using multiple crystal forms, temperatures, and pHs, across the apo and ADP-ribose-bound states. We characterize extensive solvation in the Mac1 active site, and visualize how water networks reorganize upon binding of ADP-ribose and non-native ligands, inspiring strategies for displacing waters to increase potency of Mac1 inhibitors. Determining the precise orientations of active site water molecules and the protonation states of key catalytic site residues by neutron crystallography suggests a catalytic mechanism for coronavirus macrodomains distinct from the substrate-assisted mechanism proposed for human MacroD2. These data provoke a re-evaluation of macrodomain catalytic mechanisms and will guide the optimization of Mac1 inhibitors.

High-Throughput Virtual Screening and Validation of a SARS-CoV-2 Main Protease Noncovalent Inhibitor

Despite the recent availability of vaccines against the acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the search for inhibitory therapeutic agents has assumed importance especially in the context of emerging new viral variants. In this paper, we describe the discovery of a novel noncovalent small-molecule inhibitor, MCULE-5948770040, that binds to and inhibits the SARS-Cov-2 main protease (Mpro) by employing a scalable high-throughput virtual screening (HTVS) framework and a targeted compound library of over 6.5 million molecules that could be readily ordered and purchased. Our HTVS framework leverages the U.S. supercomputing infrastructure achieving nearly 91% resource utilization and nearly 126 million docking calculations per hour. Downstream biochemical assays validate this Mpro inhibitor with an inhibition constant (Ki) of 2.9 μM (95% CI 2.2, 4.0). Furthermore, using room-temperature X-ray crystallography, we show that MCULE-5948770040 binds to a cleft in the primary binding site of Mpro forming stable hydrogen bond and hydrophobic interactions. We then used multiple μs-time scale molecular dynamics (MD) simulations and machine learning (ML) techniques to elucidate how the bound ligand alters the conformational states accessed by Mpro, involving motions both proximal and distal to the binding site. Together, our results demonstrate how MCULE-5948770040 inhibits Mpro and offers a springboard for further therapeutic design.

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

Structural, Electronic, and Electrostatic Determinants for Inhibitor Binding to Subsites S1 and S2 in SARS-CoV-2 Main Protease

Creating small-molecule antivirals specific for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) proteins is crucial to battle coronavirus disease 2019 (COVID-19). SARS-CoV-2 main protease (Mpro) is an established drug target for the design of protease inhibitors. We performed a structure-activity relationship (SAR) study of noncovalent compounds that bind in the enzyme's substrate-binding subsites S1 and S2, revealing structural, electronic, and electrostatic determinants of these sites. The study was guided by the X-ray/neutron structure of Mpro complexed with Mcule-5948770040 (compound 1), in which protonation states were directly visualized. Virtual reality-assisted structure analysis and small-molecule building were employed to generate analogues of 1. In vitro enzyme inhibition assays and room-temperature X-ray structures demonstrated the effect of chemical modifications on Mpro inhibition, showing that (1) maintaining correct geometry of an inhibitor's P1 group is essential to preserve the hydrogen bond with the protonated His163; (2) a positively charged linker is preferred; and (3) subsite S2 prefers nonbulky modestly electronegative groups.

Structural and functional characterization of NEMO cleavage by SARS-CoV-2 3CLpro.. [PMID: 34816264 PMCID: PMC8609902 DOI: 10.1101/2021.11.11.468228] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

In addition to its essential role in viral polyprotein processing, the SARS-CoV-2 3C-like (3CLpro) protease can cleave human immune signaling proteins, like NF-κB Essential Modulator (NEMO) and deregulate the host immune response. Here, in vitro assays show that SARS-CoV-2 3CLpro cleaves NEMO with fine-tuned efficiency. Analysis of the 2.14 Å resolution crystal structure of 3CLpro C145S bound to NEMO_226–235 reveals subsites that tolerate a range of viral and host substrates through main chain hydrogen bonds while also enforcing specificity using side chain hydrogen bonds and hydrophobic contacts. Machine learning- and physics-based computational methods predict that variation in key binding residues of 3CLpro-NEMO helps explain the high fitness of SARS-CoV-2 in humans. We posit that cleavage of NEMO is an important piece of information to be accounted for in the pathology of COVID-19.

Michaelis-like complex of SARS-CoV-2 main protease visualized by room-temperature X-ray crystallography

SARS-CoV-2 emerged at the end of 2019 to cause an unprecedented pandemic of the deadly respiratory disease COVID-19 that continues to date. The viral main protease (Mpro) is essential for SARS-CoV-2 replication and is therefore an important drug target. Understanding the catalytic mechanism of Mpro, a cysteine protease with a catalytic site comprising the noncanonical Cys145-His41 dyad, can help in guiding drug design. Here, a 2.0 Å resolution room-temperature X-ray crystal structure is reported of a Michaelis-like complex of Mpro harboring a single inactivating mutation C145A bound to the octapeptide Ac-SAVLQSGF-CONH2 corresponding to the nsp4/nsp5 autocleavage site. The peptide substrate is unambiguously defined in subsites S5 to S3' by strong electron density. Superposition of the Michaelis-like complex with the neutron structure of substrate-free Mpro demonstrates that the catalytic site is inherently pre-organized for catalysis prior to substrate binding. Induced fit to the substrate is driven by P1 Gln binding in the predetermined subsite S1 and rearrangement of subsite S2 to accommodate P2 Leu. The Michaelis-like complex structure is ideal for in silico modeling of the SARS-CoV-2 Mpro catalytic mechanism.

Identification and structural analysis of a thermophilic β-1,3-glucanase from compost

β-1,3-glucanase can specifically hydrolyze glucans to oligosaccharides and has potential applications in biotechnology. We used the metatranscriptomic technology to discover a thermophilic β-1,3-glucanase from compost. The phylogenetic study shows that it belongs to the family 16 glycoside hydrolase (GH16) and is most homologous with an enzyme from Streptomyces sioyaensis, an actinobacterium. It has the activity of 146.9 U/mg in the optimal reaction condition (75 °C and pH 5.5). Its catalytic domain was crystallized and diffracted to 1.14 Å resolution. The crystal structure shows a sandwich-like β-jelly-roll fold with two disulfide bonds. After analyzing the occurring frequencies of these cysteine residues, we designed two mutants (C160G and C180I) to study the role of these disulfide bonds. Both mutants have decreased their optimal temperature from 75 to 70 °C, which indicate that the disulfide bonds are important to maintain thermostability. Interestingly, the activity of C160G has increased ~ 17% to reach 171.4 U/mg. We speculate that the increased activity of C160G mutant is due to increased dynamics near the active site. Our studies give a good example of balancing the rigidity and flexibility for enzyme activity, which is helpful for protein engineering.

Characterization and structural analysis of a thermophilic GH11 xylanase from compost metatranscriptome

Xylanase is efficient for xylan degradation and widely applied in industries. We found a GH11 family xylanase (Xyn11A) with high thermostability and catalytic activity from compost metatranscriptome. This xylanase has the optimal reaction temperature at 80 °C with the activity of 2907.3 U/mg. The X-ray crystallographic structure shows a typical "right hand" architecture, which is the characteristics of the GH11 family enzymes. Comparing it with the mesophilic XYN II, a well-studied GH11 xylanase from Trichoderma reesei, Xyn11A is more compact with more H-bonds. Our mutagenic results show that the electrostatic interactions in the thumb and palm region of Xyn11A could result in its high thermostability and activity. Introducing a disulfide bond at the N-terminus further increased its optimal reaction temperature to 90 °C with augmented activity. KEY POINTS: • A hyperthermophilic xylanase with high activity was discovered using the metatranscriptomic method. • The mechanisms of thermophilicity and high activity were revealed using X-ray crystallography, mutagenesis, and molecular dynamics simulations. • The thermostability and activity were further improved by introducing a disulfide bond.

Room temperature crystallography of human acetylcholinesterase bound to a substrate analogue 4K-TMA: Towards a neutron structure

Acetylcholinesterase (AChE) catalyzes hydrolysis of acetylcholine thereby terminating cholinergic nerve impulses for efficient neurotransmission. Human AChE (hAChE) is a target of nerve agent and pesticide organophosphorus compounds that covalently attach to the catalytic Ser203 residue. Reactivation of inhibited hAChE can be achieved with nucleophilic antidotes, such as oximes. Understanding structural and electrostatic (i.e. protonation states) determinants of the catalytic and reactivation processes is crucial to improve design of oxime reactivators. Here we report X-ray structures of hAChE conjugated with a reversible covalent inhibitor 4K-TMA (4K-TMA:hAChE) at 2.8 ​Å resolution and of 4K-TMA:hAChE conjugate with oxime reactivator methoxime, MMB4 (4K-TMA:hAChE:MMB4) at 2.6 ​Å resolution, both at physiologically relevant room temperature, as well as cryo-crystallographic structure of 4K-TMA:hAChE at 2.4 ​Å resolution. 4K-TMA acts as a substrate analogue reacting with the hydroxyl of Ser203 and generating a reversible tetrahedral hemiketal intermediate that closely resembles the first tetrahedral intermediate state during hAChE-catalyzed acetylcholine hydrolysis. Structural comparisons of room temperature with cryo-crystallographic structures of 4K-TMA:hAChE and published mAChE complexes with 4K-TMA, as well as the effect of MMB4 binding to the peripheral anionic site (PAS) of the 4K-TMA:hAChE complex, revealed only discrete, minor differences. The active center geometry of AChE, already highly evolved for the efficient catalysis, was thus indicative of only minor conformational adjustments to accommodate the tetrahedral intermediate in the hydrolysis of the neurotransmitter acetylcholine (ACh). To map protonation states in the hAChE active site gorge we collected 3.5 ​Å neutron diffraction data paving the way for obtaining higher resolution datasets that will be needed to determine locations of individual hydrogen atoms.

Revertant mutation V48G alters conformational dynamics of highly drug resistant HIV protease PRS17

Drug resistance is a serious problem for controlling the HIV/AIDS pandemic. Current antiviral drugs show several orders of magnitude worse inhibition of highly resistant clinical variant PRS17 of HIV-1 protease compared with wild-type protease. We have analyzed the effects of a common resistance mutation G48V in the flexible flaps of the protease by assessing the revertant PRS17_V48G for changes in enzyme kinetics, inhibition, structure, and dynamics. Both PRS17 and the revertant showed about 10-fold poorer catalytic efficiency than wild-type enzyme (0.55 and 0.39 μM^-1min^-1 compared to 6.3 μM^-1min^-1). Clinical inhibitors, amprenavir and darunavir, showed 2-fold and 8-fold better inhibition, respectively, of the revertant than of PRS17, although the inhibition constants for PRS17_V48G were still 25 to 1,200-fold worse than for wild-type protease. Crystal structures of inhibitor-free revertant and amprenavir complexes with revertant and PRS17 were solved at 1.3-1.5 Å resolution. The amprenavir complexes of PRS17_V48G and PRS17 showed no significant differences in the interactions with inhibitor, although changes were observed in the conformation of Phe53 and the interactions of the flaps. The inhibitor-free structure of the revertant showed flaps in an open conformation, however, the flap tips do not have the unusual curled conformation seen in inhibitor-free PRS17. Molecular dynamics simulations were run for 1 μs on the two inhibitor-free mutants and wild-type protease. PRS17 exhibited higher conformational fluctuations than the revertant, while the wild-type protease adopted the closed conformation and showed the least variation. The second half of the simulations captured the transition of the flaps of PRS17 from a closed to a semi-open state, whereas the flaps of PRS17_V48G tucked into the active site and the wild-type protease retained the closed conformation. These results suggest that mutation G48V contributes to drug resistance by altering the conformational dynamics of the flaps.

Covalent inhibition of hAChE by organophosphates causes homodimer dissociation through long-range allosteric effects

Acetylcholinesterase (EC 3.1.1.7), a key acetylcholine-hydrolyzing enzyme in cholinergic neurotransmission, is present in a variety of states in situ, including monomers, C-terminally disulfide-linked homodimers, homotetramers, and up to three tetramers covalently attached to structural subunits. Could oligomerization that ensures high local concentrations of catalytic sites necessary for efficient neurotransmission be affected by environmental factors? Using small-angle X-ray scattering (SAXS) and cryo-EM, we demonstrate that homodimerization of recombinant monomeric human acetylcholinesterase (hAChE) in solution occurs through a C-terminal four-helix bundle at micromolar concentrations. We show that diethylphosphorylation of the active serine in the catalytic gorge or isopropylmethylphosphonylation by the R_P enantiomer of sarin promotes a 10-fold increase in homodimer dissociation. We also demonstrate the dissociation of organophosphate (OP)-conjugated dimers is reversed by structurally diverse oximes 2PAM, HI6, or RS194B, as demonstrated by SAXS of diethylphosphoryl-hAChE. However, binding of oximes to the native ligand-free hAChE, binding of high-affinity reversible ligands, or formation of an S_P-sarin-hAChE conjugate had no effect on homodimerization. Dissociation monitored by time-resolved SAXS occurs in milliseconds, consistent with rates of hAChE covalent inhibition. OP-induced dissociation was not observed in the SAXS profiles of the double-mutant Y337A/F338A, where the active center gorge volume is larger than in wildtype hAChE. These observations suggest a key role of the tightly packed acyl pocket in allosterically triggered OP-induced dimer dissociation, with the potential for local reduction of acetylcholine-hydrolytic power in situ. Computational models predict allosteric correlated motions extending from the acyl pocket toward the four-helix bundle dimerization interface 25 Å away.

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.

Direct Observation of Protonation State Modulation in SARS-CoV-2 Main Protease upon Inhibitor Binding with Neutron Crystallography

The main protease (3CL Mpro) from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, is an essential enzyme for viral replication with no human counterpart, making it an attractive drug target. To date, no small-molecule clinical drugs are available that specifically inhibit SARS-CoV-2 Mpro. To aid rational drug design, we determined a neutron structure of Mpro in complex with the α-ketoamide inhibitor telaprevir at near-physiological (22 °C) temperature. We directly observed protonation states in the inhibitor complex and compared them with those in the ligand-free Mpro, revealing modulation of the active-site protonation states upon telaprevir binding. We suggest that binding of other α-ketoamide covalent inhibitors can lead to the same protonation state changes in the Mpro active site. Thus, by studying the protonation state changes induced by inhibitors, we provide crucial insights to help guide rational drug design, allowing precise tailoring of inhibitors to manipulate the electrostatic environment of SARS-CoV-2 Mpro.

Inhibitor binding influences the protonation states of histidines in SARS-CoV-2 main protease

The main protease (M^pro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of M^pro, a cysteine protease, have been determined, facilitating structure-based drug design. M^pro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41–Cys145, M^pro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nucleophile Cys145 have been debated in previous studies of SARS-CoV M^pro, but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 M^pro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of M^pro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α-ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored N_δ (HD) and N_ϵ (HE) protonation of His41 and His164, respectively, the α-ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 M^pro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.

The main protease (M^pro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics.

Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19

We present a supercomputer-driven pipeline for in silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. Ensemble docking makes use of MD results by docking compound databases into representative protein binding-site conformations, thus taking into account the dynamic properties of the binding sites. We also describe preliminary results obtained for 24 systems involving eight proteins of the proteome of SARS-CoV-2. The MD involves temperature replica exchange enhanced sampling, making use of massively parallel supercomputing to quickly sample the configurational space of protein drug targets. Using the Summit supercomputer at the Oak Ridge National Laboratory, more than 1 ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to 10 configurations of each of the 24 SARS-CoV-2 systems using AutoDock Vina. Comparison to experiment demonstrates remarkably high hit rates for the top scoring tranches of compounds identified by our ensemble approach. We also demonstrate that, using Autodock-GPU on Summit, it is possible to perform exhaustive docking of one billion compounds in under 24 h. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and artificial intelligence (AI) methods to cluster MD trajectories and rescore docking poses.

Unusual zwitterionic catalytic site of SARS-CoV-2 main protease revealed by neutron crystallography

The main protease (3CL Mpro) from SARS-CoV-2, the etiological agent of COVID-19, is an essential enzyme for viral replication. 3CL Mpro possesses an unusual catalytic dyad composed of Cys145 and His41 residues. A critical question in the field has been what the protonation states of the ionizable residues in the substrate-binding active-site cavity are; resolving this point would help understand the catalytic details of the enzyme and inform rational drug development against this pernicious virus. Here, we present the room-temperature neutron structure of 3CL Mpro, which allowed direct determination of hydrogen atom positions and, hence, protonation states in the protease. We observe that the catalytic site natively adopts a zwitterionic reactive form in which Cys145 is in the negatively charged thiolate state and His41 is doubly protonated and positively charged, instead of the neutral unreactive state usually envisaged. The neutron structure also identified the protonation states, and thus electrical charges, of all other amino acid residues and revealed intricate hydrogen-bonding networks in the active-site cavity and at the dimer interface. The fine atomic details present in this structure were made possible by the unique scattering properties of the neutron, which is an ideal probe for locating hydrogen positions and experimentally determining protonation states at near-physiological temperature. Our observations provide critical information for structure-assisted and computational drug design, allowing precise tailoring of inhibitors to the enzyme's electrostatic environment.

Malleability of the SARS-CoV-2 3CL Mpro Active-Site Cavity Facilitates Binding of Clinical Antivirals

The COVID-19 pandemic caused by SARS-CoV-2 requires rapid development of specific therapeutics and vaccines. The main protease of SARS-CoV-2, 3CL Mpro, is an established drug target for the design of inhibitors to stop the virus replication. Repurposing existing clinical drugs can offer a faster route to treatments. Here, we report on the binding mode and inhibition properties of several inhibitors using room temperature X-ray crystallography and in vitro enzyme kinetics. The enzyme active-site cavity reveals a high degree of malleability, allowing aldehyde leupeptin and hepatitis C clinical protease inhibitors (telaprevir, narlaprevir, and boceprevir) to bind and inhibit SARS-CoV-2 3CL Mpro. Narlaprevir, boceprevir, and telaprevir are low-micromolar inhibitors, whereas the binding affinity of leupeptin is substantially weaker. Repurposing hepatitis C clinical drugs as COVID-19 treatments may be a useful option to pursue. The observed malleability of the enzyme active-site cavity should be considered for the successful design of specific protease inhibitors.

Studying the Role of a Single Mutation of a Family 11 Glycoside Hydrolase Using High-Resolution X-ray Crystallography

XynII is a family 11 glycoside hydrolase that uses the retaining mechanism for catalysis. In the active site, E177 works as the acid/base and E86 works as the nucleophile. Mutating an uncharged residue (N44) to an acidic residue (D) near E177 decreases the enzyme's optimal pH by ~ 1.0 unit. D44 was previously suggested to be a second proton carrier for catalysis. To test this hypothesis, we abolished the activity of E177 by mutating it to be Q, and mutated N44 to be D or E. These double mutants have dramatically decreased activities. Our high-resolution crystallographic structures and the microscopic pK_a calculations show that D44 has similar position and pK_a value during catalysis, indicating that D44 changes electrostatics around E177, which makes it prone to rotate as the acid/base in acidic conditions, thus decreases the pH optimum. Our results could be helpful to design enzymes with different pH optimum.

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.

Room-temperature X-ray crystallography reveals the oxidation and reactivity of cysteine residues in SARS-CoV-2 3CL Mpro: insights into enzyme mechanism and drug design

The emergence of the novel coronavirus SARS-CoV-2 has resulted in a worldwide pandemic not seen in generations. Creating treatments and vaccines to battle COVID-19, the disease caused by the virus, is of paramount importance in order to stop its spread and save lives. The viral main protease, 3CL Mpro, is indispensable for the replication of SARS-CoV-2 and is therefore an important target for the design of specific protease inhibitors. Detailed knowledge of the structure and function of 3CL Mpro is crucial to guide structure-aided and computational drug-design efforts. Here, the oxidation and reactivity of the cysteine residues of the protease are reported using room-temperature X-ray crystallography, revealing that the catalytic Cys145 can be trapped in the peroxysulfenic acid oxidation state at physiological pH, while the other surface cysteines remain reduced. Only Cys145 and Cys156 react with the alkylating agent N-ethylmaleimide. It is suggested that the zwitterionic Cys145-His45 catalytic dyad is the reactive species that initiates catalysis, rather than Cys145-to-His41 proton transfer via the general acid-base mechanism upon substrate binding. The structures also provide insight into the design of improved 3CL Mpro inhibitors.

Inhibitor binding influences the protonation states of histidines in SARS-CoV-2 main protease

The main protease (M ^pro ) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of M ^pro , a cysteine protease, have been determined, facilitating structure-based drug design. M ^pro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41-Cys145, M ^pro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nu-cleophile Cys145 have been debated in previous studies of SARS-CoV M ^pro , but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 M ^pro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of M ^pro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α -ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored N _δ (HD) and N _ϵ (HE) protonation of His41 and His164, respectively, the α -ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 M ^pro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.

Substrate Binding Stiffens Aspartate Aminotransferase by Altering the Enzyme Picosecond Vibrational Dynamics

Protein dynamics on various time scales from femtoseconds to milliseconds impacts biological function by driving proteins to conformations conducive to ligand binding and creating functional states in enzyme catalysis. Neutron vibrational spectroscopy carried out by measuring inelastic neutron scattering from protein molecules in combination with molecular simulations has the unique ability of detecting and visualizing changes in the picosecond protein vibrational dynamics due to ligand binding. Here we present neutron vibrational spectra of a homodimeric pyridoxal 5'-phosphate-dependent enzyme, aspartate aminotransferase, obtained from the open internal aldimine and closed external aldimine conformational states. We observe that in the external aldimine state the protein structure stiffens relative to the internal aldimine state, indicating rigidified vibrational dynamics on the picosecond time scale in the low-frequency regime of 5-50 cm^-1. Our molecular dynamics simulations indicate substantial changes in the picosecond dynamics of the enzyme secondary structure elements upon substrate binding, with the largest contributions from just two helices and the β-sheet.

Structural plasticity of SARS-CoV-2 3CL Mpro active site cavity revealed by room temperature X-ray crystallography

The COVID-19 disease caused by the SARS-CoV-2 coronavirus has become a pandemic health crisis. An attractive target for antiviral inhibitors is the main protease 3CL M^pro due to its essential role in processing the polyproteins translated from viral RNA. Here we report the room temperature X-ray structure of unliganded SARS-CoV-2 3CL M^pro, revealing the ligand-free structure of the active site and the conformation of the catalytic site cavity at near-physiological temperature. Comparison with previously reported low-temperature ligand-free and inhibitor-bound structures suggest that the room temperature structure may provide more relevant information at physiological temperatures for aiding in molecular docking studies.

The SARS-CoV-2 3CL main protease (3CL M^pro) is a chymotrypsin-like protease that facilitates the production of non-structural proteins, which are essential for viral replication and is therefore of great interest as a drug target. Here, the authors present the 2.30 Å room temperature crystal structure of ligand-free 3CL M^pro and compare it with the earlier determined low-temperature ligand-free and inhibitor-bound crystal structures.

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.

Rational design, synthesis, and evaluation of uncharged, "smart" bis-oxime antidotes of organophosphate-inhibited human acetylcholinesterase

Organophosphate (OP) intoxications from nerve agent and OP pesticide exposures are managed with pyridinium aldoxime-based therapies whose success rates are currently limited. The pyridinium cation hampers uptake of OPs into the central nervous system (CNS). Furthermore, it frequently binds to aromatic residues of OP-inhibited acetylcholinesterase (AChE) in orientations that are nonproductive for AChE reactivation, and the structural diversity of OPs impedes efficient reactivation. Improvements of OP antidotes need to include much better access of AChE reactivators to the CNS and optimized orientation of the antidotes' nucleophile within the AChE active-center gorge. On the basis of X-ray structures of a CNS-penetrating reactivator, monoxime RS194B, reversibly bound to native and venomous agent X (VX)-inhibited human AChE, here we created seven uncharged acetamido bis-oximes as candidate antidotes. Both oxime groups in these bis-oximes were attached to the same central, saturated heterocyclic core. Diverse protonation of the heterocyclic amines and oxime groups of the bis-oximes resulted in equilibration among up to 16 distinct ionization forms, including uncharged forms capable of diffusing into the CNS and multiple zwitterionic forms optimal for reactivation reactions. Conformationally diverse zwitterions that could act as structural antidote variants significantly improved in vitro reactivation of diverse OP-human AChE conjugates. Oxime group reorientation of one of the bis-oximes, forcing it to point into the active center for reactivation, was confirmed by X-ray structural analysis. Our findings provide detailed structure-activity properties of several CNS-directed, uncharged aliphatic bis-oximes holding promise for use as protonation-dependent, conformationally adaptive, "smart" accelerated antidotes against OP toxicity.

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.

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.

IMAGINE: The neutron protein crystallography beamline at the high flux isotope reactor

IMAGINE is a high intensity, quasi-Laue neutron crystallography beamline developed at the 85MW High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (ORNL). This state-of-the-art facility for neutron-diffraction enables neutron protein structures to be determined at or near atomic resolutions from crystals with volumes of <1mm³ and unit cell edges of <150Å. The beamline features include elliptical focusing mirrors that deliver neutrons into a 2.0×3.2mm² focal spot at the sample position, and variable short and long wavelength cutoff optics that provide automated exchange between multiple wavelength configurations. The beamline is equipped with a single-axis goniometer, neutron-sensitive cylindrical image plate detector and room temperature and cryogenic sample environments. This article describes the beamline components, the diffractometer and the data collection and data analysis protocols that are used, and outlines the protein deuteration, crystallization and conventional crystallography capabilities that are available to users at ORNL's neutron facilities. We also present examples of the scientific questions being addressed at this beamline and highlight important findings in enzyme chemistry that have been made possible by IMAGINE.