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

Improved joint X-ray and neutron refinement procedure in Phenix

Neutron diffraction is one of the three crystallographic techniques (X-ray, neutron and electron diffraction) used to determine the atomic structures of molecules. Its particular strengths derive from the fact that H (and D) atoms are strong neutron scatterers, meaning that their positions, and thus protonation states, can be derived from crystallographic maps. However, because of technical limitations and experimental obstacles, the quality of neutron diffraction data is typically much poorer (completeness, resolution and signal to noise) than that of X-ray diffraction data for the same sample. Further, refinement is more complex as it usually requires additional parameters to describe the H (and D) atoms. The increase in the number of parameters may be mitigated by using the `riding hydrogen' refinement strategy, in which the positions of H atoms without a rotational degree of freedom are inferred from their neighboring heavy atoms. However, this does not address the issues related to poor data quality. Therefore, neutron structure determination often relies on the presence of an X-ray data set for joint X-ray and neutron (XN) refinement. In this approach, the X-ray data serve to compensate for the deficiencies of the neutron diffraction data by refining one model simultaneously against the X-ray and neutron data sets. To be applicable, it is assumed that both data sets are highly isomorphous, and preferably collected from the same crystals and at the same temperature. However, the approach has a number of limitations that are discussed in this work by comparing four separately re-refined neutron models. To address the limitations, a new method for joint XN refinement is introduced that optimizes two different models against the different data sets. This approach is tested using neutron models and data deposited in the Protein Data Bank. The efficacy of refining models with H atoms as riding or as individual atoms is also investigated.

Structural and Synthetic Aspects of Small Ring Oxa- and Aza-Heterocyclic Ring Systems as Antiviral Activities

Antiviral properties of different oxa- and aza-heterocycles are identified and properly correlated with their structural features and discussed in this review article. The primary objective is to explore the activity of such ring systems as antiviral agents, as well as their synthetic routes and biological significance. Eventually, the structure-activity relationship (SAR) of the heterocyclic compounds, along with their salient characteristics are exhibited to build a suitable platform for medicinal chemists and biotechnologists. The synergistic conclusions are extremely important for the introduction of a newer tool for the future drug discovery program.

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.

Understanding the impact of halogen functional group (Br, Cl, F, OH) in amprenavir ligand of the HIV protease

We focused our attention towards the most dreadful disease that threatens the mankind of 20th century - Acquired immunodeficiency syndrome (AIDS), caused through the human immunodeficiency virus (HIV) and a sexually transmitted infection (STI). In this study, our foremost interest was to identify the potency and stability of HIV ligand- Amprenavir (APV) and its modelled functional group (Br, Cl, F, CF3, CH3, NH2) ligands through halogen and hydrogen bond contact, which will have a clear portrait on the structure activity of protein ligand interactions. This will assist chemist in synthesizing novel APV ligands, which are expected to inhibit the activity of HIV-1 protease enzyme. The binding strength of Amprenavir ligand with interacting hinge region amino acid side chains: Isoleucine (ILE 147, 150, 184), Valine (VAL 82), Alanine (ALA 28), Aspartic acid (25, 30, 125, 130) and Glycine (GLY 127, 149) were understood through interaction energy calculations at HF, B3LYP, M052X, MP2 level of theories for different basis set (6-311 G**, LANL2DZ). The present work will reveal an understandable picture about the halogen and hydrogen bond interaction that grip the contact of ligand and amino acids in the hinge region. Overall the Halogen atom (Br, Cl, F) functional groups improved the binding strength of APV in HIV protease; which provide a new novel path for the functional group preference on the ligand that enclose perfectly with the amino acid in the hinge region.Communicated by Ramaswamy H. Sarma.

Complementarity of neutron, XFEL and synchrotron crystallography for defining the structures of metalloenzymes at room temperature

Room-temperature macromolecular crystallography allows protein structures to be determined under close-to-physiological conditions, permits dynamic freedom in protein motions and enables time-resolved studies. In the case of metalloenzymes that are highly sensitive to radiation damage, such room-temperature experiments can present challenges, including increased rates of X-ray reduction of metal centres and site-specific radiation-damage artefacts, as well as in devising appropriate sample-delivery and data-collection methods. It can also be problematic to compare structures measured using different crystal sizes and light sources. In this study, structures of a multifunctional globin, dehaloperoxidase B (DHP-B), obtained using several methods of room-temperature crystallographic structure determination are described and compared. Here, data were measured from large single crystals and multiple microcrystals using neutrons, X-ray free-electron laser pulses, monochromatic synchrotron radiation and polychromatic (Laue) radiation light sources. These approaches span a range of 18 orders of magnitude in measurement time per diffraction pattern and four orders of magnitude in crystal volume. The first room-temperature neutron structures of DHP-B are also presented, allowing the explicit identification of the hydrogen positions. The neutron data proved to be complementary to the serial femtosecond crystallography data, with both methods providing structures free of the effects of X-ray radiation damage when compared with standard cryo-crystallography. Comparison of these room-temperature methods demonstrated the large differences in sample requirements, data-collection time and the potential for radiation damage between them. With regard to the structure and function of DHP-B, despite the results being partly limited by differences in the underlying structures, new information was gained on the protonation states of active-site residues which may guide future studies of DHP-B.

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

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

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

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

A protocol for production of perdeuterated OmpF porin for neutron crystallography

Hydrogen atoms are at the limit of visibility in X-ray structures even at high resolution. Neutron macromolecular crystallography (NMX) is an unambiguous method to locate hydrogens and study the significance of hydrogen bonding interactions in biological systems. Since NMX requires very large crystals, very few neutron structures of proteins have been determined yet. In addition, the most common hydrogen isotope ¹H gives rise to significant background due to its large incoherent scattering cross-section. Therefore, it is advantageous to substitute as many hydrogens as possible with the heavier isotope ²H (deuterium) to reduce the sample volume requirement. While the solvent exchangeable hydrogens can be substituted by dissolving the protein in heavy water, complete deuterium labelling - perdeuteration - requires the protein to be expressed in heavy water with a deuterated carbon source. In this work, we developed an optimized method for large scale production of deuterium-labelled bacterial outer membrane protein F (OmpF) for NMX. OmpF was produced using deuterated media with different carbon sources. Mass spectrometry verified the integrity and level of deuteration of purified OmpF. Perdeuterated OmpF crystals diffracted X-rays to a resolution of 1.9 Å. This work lays the foundation for structural studies of membrane protein by neutron diffraction in future.

Insights into effect of the Asp25/Asp25' protonation states on binding of inhibitors Amprenavir and MKP97 to HIV-1 protease using molecular dynamics simulations and MM-GBSA calculations

The protonation states of two aspartic acids in the catalytic strands of HIV-1 protease (PR) remarkably affect bindings of inhibitors to PR. It is requisite for the design of potent inhibitors towards PR to investigate the influences of Asp25/Asp25' protonated states on dynamics behaviour of PR and binding mechanism of inhibitors to PR. In this work, molecular dynamics (MD) simulations, MM-GBSA method and principal component (PC) analysis were coupled to explore the effect of Asp25/Asp25' protonation states on conformational changes of PR and bindings of Amprenavir and MKP97 to PR. The results show that the Asp25/Asp25' protonation states exert different impacts on structural fluctuations, flexibility and motion modes of PR. Dynamics analysis verifies that Asp25/Asp25' protonated states highly affect conformational dynamics of two flaps in PR. The binding free energy calculations results suggest that the Asp25/Asp25' protonated states obviously strengthen bindings of inhibitors to PR compared to the non-protonation state. Calculations of residue-based free energy decomposition indicate that the Asp25/Asp25' protonation not only disturbs the interaction network of inhibitors with PR but also stabilizes bindings of inhibitors to PR by cancelling the electrostatic repulsive interaction. Therefore, special attentions should be paid to the Asp25/Asp25' protonation in the design of potent inhibitors towards PR.

On the Case of the Misplaced Hydrogens

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

Opportunities and challenges in neutron crystallography

Neutron and X-ray crystallography are complementary to each other. While X-ray scattering is directly proportional to the number of electrons of an atom, neutrons interact with the atomic nuclei themselves. Neutron crystallography therefore provides an excellent alternative in determining the positions of hydrogens in a biological molecule. In particular, since highly polarized hydrogen atoms (H+) do not have electrons, they cannot be observed by X-rays. Neutron crystallography has its own limitations, mainly due to inherent low flux of neutrons sources, and as a consequence, the need for much larger crystals and for different data collection and analysis strategies. These technical challenges can however be overcome to yield crucial structural insights about protonation states in enzyme catalysis, ligand recognition, as well as the presence of unusual hydrogen bonds in proteins.

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.

Substrate-Enzyme Interactions in Intramembrane Proteolysis: γ-Secretase as the Prototype

Intramembrane-cleaving proteases (I-CLiPs) catalyze the hydrolysis of peptide bonds within the transmembrane regions of membrane protein substrates, releasing bioactive fragments that play roles in many physiological and pathological processes. Based on their catalytic mechanism and nucleophile, I-CLiPs are classified into metallo, serine, aspartyl, and glutamyl proteases. Presenilin is the most prominent among I-CLiPs, as the catalytic subunit of γ-secretase (GS) complex responsible for cleaving the amyloid precursor protein (APP) and Notch, as well as many other membrane substrates. Recent cryo-electron microscopy (cryo-EM) structures of GS provide new details on how presenilin recognizes and cleaves APP and Notch. First, presenilin transmembrane helix (TM) 2 and 6 are dynamic. Second, upon binding to GS, the substrate TM helix is unwound from the C-terminus, resulting in an intermolecular β-sheet between the substrate and presenilin. The transition of the substrate C-terminus from α-helix to β-sheet is proposed to expose the scissile peptide bond in an extended conformation, leaving it susceptible to protease cleavage. Despite the astounding new insights in recent years, many crucial questions remain unanswered regarding the inner workings of γ-secretase, however. Key unanswered questions include how the enzyme recognizes and recruits substrates, how substrates are translocated from an initial docking site to the active site, how active site aspartates recruit and coordinate catalytic water, and the nature of the mechanisms of processive trimming of the substrate and product release. Answering these questions will have important implications for drug discovery aimed at selectively reducing the amyloid load in Alzheimer's disease (AD) with minimal side effects.

Inhibition of the activity of HIV-1 protease through antibody binding and mutations probed by molecular dynamics simulations

HIV-1 protease is an essential enzyme in the life cycle of the HIV-1 virus. The conformational dynamics of the flap region of the protease is critical for the ligand binding mechanism, as well as for the catalytic activity. The monoclonal antibody F11.2.32 raised against HIV-1 protease inhibits its activity on binding. We have studied the conformational dynamics of protease in its free, inhibitor ritonavir and antibody bound forms using molecular dynamics simulations. We find that upon Ab binding to the epitope region (residues 36-46) of protease, the overall flexibility of the protease is decreased including the flap region and the active site, which is similar to the decrease in flexibility observed by inhibitor binding to the protease. This suggests an allosteric mechanism to inhibit protease activity. Further, the protease mutants G40E and G40R are known to have decreased activity and were also subjected to MD simulations. We find that the loss of flexibility in the mutants is similar to that observed in the protease bound to the Ab/inhibitor. These insights highlight the role played by dynamics in the function of the protease and how control of flexibility through Ab binding and site specific mutations can inhibit protease activity.

Performance of radial distribution function-based descriptors in the chemoinformatic studies of HIV-1 protease

Aim: This letter investigates the role of radial distribution function-based descriptors for in silico design of new drugs. Methodology: The multiple linear regression models for HIV-1 protease and its complexes with a series of inhibitors were constructed. A detailed analysis of major atomic contributions to the radial distribution function descriptor weighted by the number of valence shell electrons identified residues Arg8, Asp29 and residues of the catalytic triad as crucial for the correlation with the inhibition constant, together with residues Asp30 and Ile50, whose mutations are known to cause an emergence of drug resistant variants. Conclusion: This study demonstrates an easy and fast assessment of the activity of potential drugs and the derivation of structural information of their complexes with the receptor or enzyme.

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.

The Influence of Hydrogen Atoms on the Performance of Radial Distribution Function-Based Descriptors in the Chemoinformatic Studies of HIV-1 Protease Complexes with Inhibitors

AIMS

The aim of this letter is to explore the influence of adding hydrogen atoms to the crystallographic structures of HIV-1 protease complexes with a series of inhibitors on the performance of radial distribution function based descriptors recently introduced in chemoinformatic studies.

BACKGROUND

Quite recently the successful application of molecular descriptors based on a radial distribution function to correlate it with biologically interesting properties of a ligand - enzyme complex was demonstrated. Except its predictive power, the analysis of atoms with dominant contributions to the RDFs can be used to identify relevant atoms and interactions. Since original paper was published on dataset consisting of the X-ray structures of complexes without hydrogen atoms, we wonder weather addition of light atoms can provide us new piece of information.

OBJECTIVE

The primarily objective is to create the model correlating the RDF based descriptors and physicochemical properties of the HIV-1 protease complexes with inhibitors with hydrogen atoms. Then, we will compare the performance of new model with previous one, where the hydrogen atoms were discarded. Information about interactions between the enzyme and the inhibitors will be extracted from the analysis of the RDF.

METHODS

The radial distribution function descriptor weighted by the number of valence shell electrons has proven to be sensitive to the changes in the structure of the enzyme and enzyme-ligand complexes. For each structure in our data set, RDF will be calculated and using multiple linear regression method the mathematical model will be designed correlating RDF based descriptors and the physicochemical properties. Statistical analysis of the atom's contribution to the total RDF will reveal relevant interactions.

RESULTS

The applicability of RDF based descriptor for the correlation with pKi and EC50 values is demonstrated, while simple models containing only two or three parameters are able to explain 78 and 86 % of the variance, respectively. The models with explicitly included hydrogens are of comparable quality with the previous models without hydrogens. The analysis of the atom's dominant contributions highlighted the importance of the hydroxyl groups of the inhibitor near the Asp25 and Asp25' residues when it is bounded to the protease.

CONCLUSION

Models based on the RDF weighted by the number of valence shell electrons for correlating small number of molecular descriptors and physicocehmical properties for structures with and without hydrogens are of comparable quality and both can be used for identification of relevant functional groups and interactions. Other: Our approach can be integrated to the next generation virtual screening methods, because is fast, reliable with high predictability potential.

Productive reorientation of a bound oxime reactivator revealed in room temperature X-ray structures of native and VX-inhibited human acetylcholinesterase

Exposure to organophosphorus compounds (OPs) may be fatal if untreated, and a clear and present danger posed by nerve agent OPs has become palpable in recent years. OPs inactivate acetylcholinesterase (AChE) by covalently modifying its catalytic serine. Inhibited AChE cannot hydrolyze the neurotransmitter acetylcholine leading to its build-up at the cholinergic synapses and creating an acute cholinergic crisis. Current antidotes, including oxime reactivators that attack the OP-AChE conjugate to free the active enzyme, are inefficient. Better reactivators are sought, but their design is hampered by a conformationally rigid portrait of AChE extracted exclusively from 100K X-ray crystallography and scarcity of structural knowledge on human AChE (hAChE). Here, we present room temperature X-ray structures of native and VX-phosphonylated hAChE with an imidazole-based oxime reactivator, RS-170B. We discovered that inhibition with VX triggers substantial conformational changes in bound RS-170B from a "nonproductive" pose (the reactive aldoxime group points away from the VX-bound serine) in the reactivator-only complex to a "semi-productive" orientation in the VX-modified complex. This observation, supported by concurrent molecular simulations, suggested that the narrow active-site gorge of hAChE may be significantly more dynamic than previously thought, allowing RS-170B to reorient inside the gorge. Furthermore, we found that small molecules can bind in the choline-binding site hindering approach to the phosphorous of VX-bound serine. Our results provide structural and mechanistic perspectives on the reactivation of OP-inhibited hAChE and demonstrate that structural studies at physiologically relevant temperatures can deliver previously overlooked insights applicable for designing next-generation antidotes.

Membrane-protein crystals for neutron diffraction

Neutron macromolecular crystallography (NMX) has the potential to provide the experimental input to address unresolved aspects of transport mechanisms and protonation in membrane proteins. However, despite this clear scientific motivation, the practical challenges of obtaining crystals that are large enough to make NMX feasible have so far been prohibitive. Here, the potential impact on feasibility of a more powerful neutron source is reviewed and a strategy for obtaining larger crystals is formulated, exemplified by the calcium-transporting ATPase SERCA1. The challenges encountered at the various steps in the process from crystal nucleation and growth to crystal mounting are explored, and it is demonstrated that NMX-compatible membrane-protein crystals can indeed be obtained.

Conformational and dynamical basis for cross-reactivity observed between anti HIV-1 protease antibody with protease and an epitope peptide from it

F11.2.32 is a monoclonal antibody raised against HIV-1 protease and it inhibits protease activity. While the structure of the epitope peptide in complex with the antibody is known, how protease interacts with the antibody is not known. In this study, we model the conformational features of the free and bound epitope peptide and protease-antibody interactions. We find through our simulations, that the free epitope peptide P36-46 samples conformations akin to the bound conformation of the peptide in complex with the Ab, with a β-turn conformation sampled by the ³⁸LPGR⁴¹ sequence highlighting the role of inherent conformational preferences of the peptide. Further, to determine the interactions present between the protease and antibody, we docked the protease in its conformation observed in the crystal structure, onto the antibody and simulated the dynamics of the complex in explicit water. We have identified the key residues involved in hydrogen-bond interactions and salt-bridges in Ag-Ab complex and examined the role of CDR flexibility in binding different conformations of the same epitope sequence in peptide and protein antigens. Thus, our results provide the basis for understanding the cross-reactivity observed between the antibody with protease and the epitope peptide from it.

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

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

Evaluation of models determined by neutron diffraction and proposed improvements to their validation and deposition

The Protein Data Bank (PDB) contains a growing number of models that have been determined using neutron diffraction or a hybrid method that combines X-ray and neutron diffraction. The advantage of neutron diffraction experiments is that the positions of all atoms can be determined, including H atoms, which are hardly detectable by X-ray diffraction. This allows the determination of protonation states and the assignment of H atoms to water molecules. Because neutrons are scattered differently by hydrogen and its isotope deuterium, neutron diffraction in combination with H/D exchange can provide information on accessibility, dynamics and chemical lability. In this study, the deposited data, models and model-to-data fit for all PDB entries that used neutron diffraction as the source of experimental data have been analysed. In many cases, the reported Rwork and Rfree values were not reproducible. In such cases, the model and data files were analysed to identify the reasons for this mismatch. The issues responsible for the discrepancies are summarized and explained. The analysis unveiled limitations to the annotation, deposition and validation of models and data, and a lack of community-wide accepted standards for the description of neutron models and data, as well as deficiencies in current model refinement tools. Most of the issues identified concern the handling of H atoms. Since the primary use of neutron macromolecular crystallography is to locate and directly visualize H atoms, it is important to address these issues, so that the deposited neutron models allow the retrieval of the maximum amount of information with the smallest effort of manual intervention. A path forward to improving the annotation, validation and deposition of neutron models and hybrid X-ray and neutron models is suggested.

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

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

Neutron Crystallography Detects Differences in Protein Dynamics: Structure of the PKG II Cyclic Nucleotide Binding Domain in Complex with an Activator

As one of the main receptors of a second messenger, cGMP, cGMP-dependent protein kinase (PKG) isoforms I and II regulate distinct physiological processes. The design of isoform-specific activators is thus of great biomedical importance and requires detailed structural information about PKG isoforms bound with activators, including accurate positions of hydrogen atoms and a description of the hydrogen bonding and water architecture. Here, we determined a 2.2 Å room-temperature joint X-ray/neutron (XN) structure of the human PKG II carboxyl cyclic nucleotide binding (CNB-B) domain bound with a potent PKG II activator, 8-pCPT-cGMP. The XN structure directly visualizes intermolecular interactions and reveals changes in hydrogen bonding patterns upon comparison to the X-ray structure determined at cryo-temperatures. Comparative analysis of the backbone hydrogen/deuterium exchange patterns in PKG II:8-pCPT-cGMP and previously reported PKG Iβ:cGMP XN structures suggests that the ability of these agonists to activate PKG is related to how effectively they quench dynamics of the cyclic nucleotide binding pocket and the surrounding regions.

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.

Design and synthesis of potent HIV-1 protease inhibitors with (S)-tetrahydrofuran-tertiary amine-acetamide as P2-ligand: Structure-activity studies and biological evaluation

The design, synthesis, and SAR study of a new series of HIV-1 protease inhibitors incorporating stereochemically defined tetrahydrofuran-tertiary amine-acetamide P2-ligand are described. Various substituent effects on the tertiary amine P2-ligand and phenylsulfonamide P2'-ligand were investigated to maximize the ligand-binding site interactions in the protease active site. Most of inhibitors displayed low nanomolar to subnanomolar inhibitory potency. Inhibitor 20e containing N-(S-tetrahydrofuran)-N-(2-methoxyethyl)acetamide as P2-ligand along with 4-methoxylphenylsulfonamide as P2'-ligand displayed the most potent enzyme inhibitory activity (IC₅₀ = 0.35 nM) and remarkably low cytotoxicity (CC₅₀ = 305 μM).

Neutron Crystallography for the Study of Hydrogen Bonds in Macromolecules

The hydrogen bond (H bond) is one of the most important interactions that form the foundation of secondary and tertiary protein structure. Beyond holding protein structures together, H bonds are also intimately involved in solvent coordination, ligand binding, and enzyme catalysis. The H bond by definition involves the light atom, H, and it is very difficult to study directly, especially with X-ray crystallographic techniques, due to the poor scattering power of H atoms. Neutron protein crystallography provides a powerful, complementary tool that can give unambiguous information to structural biologists on solvent organization and coordination, the electrostatics of ligand binding, the protonation states of amino acid side chains and catalytic water species. The method is complementary to X-ray crystallography and the dynamic data obtainable with NMR spectroscopy. Also, as it gives explicit H atom positions, it can be very valuable to computational chemistry where exact knowledge of protonation and solvent orientation can make a large difference in modeling. This article gives general information about neutron crystallography and shows specific examples of how the method has contributed to structural biology, structure-based drug design; and the understanding of fundamental questions of reaction mechanisms.

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

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

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

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

Combining X-ray and neutron crystallography with spectroscopy

X-ray protein crystallography has, through the determination of the three-dimensional structures of enzymes and their complexes, been essential to the understanding of biological chemistry. However, as X-rays are scattered by electrons, the technique has difficulty locating the presence and position of H atoms (and cannot locate H+ ions), knowledge of which is often crucially important for the understanding of enzyme mechanism. Furthermore, X-ray irradiation, through photoelectronic effects, will perturb the redox state in the crystal. By using single-crystal spectrophotometry, reactions taking place in the crystal can be monitored, either to trap intermediates or follow photoreduction during X-ray data collection. By using neutron crystallography, the positions of H atoms can be located, as it is the nuclei rather than the electrons that scatter neutrons, and the scattering length is not determined by the atomic number. Combining the two techniques allows much greater insight into both reaction mechanism and X-ray-induced photoreduction.

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

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

Neutron crystallography aids in drug design

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

Neutron structure of human carbonic anhydrase II in complex with methazolamide: mapping the solvent and hydrogen-bonding patterns of an effective clinical drug

Carbonic anhydrases (CAs; EC 4.2.1.1) catalyze the interconversion of CO2 and HCO3-, and their inhibitors have long been used as diuretics and as a therapeutic treatment for many disorders such as glaucoma and epilepsy. Acetazolamide (AZM) and methazolamide (MZM, a methyl derivative of AZM) are two of the classical CA inhibitory drugs that have been used clinically for decades. The jointly refined X-ray/neutron structure of MZM in complex with human CA isoform II (hCA II) has been determined to a resolution of 2.2 Å with an Rcryst of ∼16.0%. Presented in this article, along with only the second neutron structure of a clinical drug-bound hCA, is an in-depth structural comparison and analyses of differences in hydrogen-bonding network, water-molecule orientation and solvent displacement that take place upon the binding of AZM and MZM in the active site of hCA II. Even though MZM is slightly more hydrophobic and displaces more waters than AZM, the overall binding affinity (Ki) for both of the drugs against hCA II is similar (∼10 nM). The plausible reasons behind this finding have also been discussed using molecular dynamics and X-ray crystal structures of hCA II-MZM determined at cryotemperature and room temperature. This study not only allows a direct comparison of the hydrogen bonding, protonation states and solvent orientation/displacement of AZM and MZM, but also shows the significant effect that the methyl derivative has on the solvent organization in the hCA II active site.

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

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

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

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

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.

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

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

Large-volume protein crystal growth for neutron macromolecular crystallography

Neutron macromolecular crystallography (NMC) is the prevailing method for the accurate determination of the positions of H atoms in macromolecules. As neutron sources are becoming more available to general users, finding means to optimize the growth of protein crystals to sizes suitable for NMC is extremely important. Historically, much has been learned about growing crystals for X-ray diffraction. However, owing to new-generation synchrotron X-ray facilities and sensitive detectors, protein crystal sizes as small as in the nano-range have become adequate for structure determination, lessening the necessity to grow large crystals. Here, some of the approaches, techniques and considerations for the growth of crystals to significant dimensions that are now relevant to NMC are revisited. These include experimental strategies utilizing solubility diagrams, ripening effects, classical crystallization techniques, microgravity and theoretical considerations.

Synthesis and biological evaluation of novel HIV-1 protease inhibitors using tertiary amine as P2-ligands

A series of tertiary amine derivatives exhibiting potent HIV-1 protease inhibiting properties were identified. These novel inhibitors were designed based on the structure of Darunavir with modification on the P2 and P2' position. This effort led to discovery of 35e and 38e, which exhibited excellent HIV-1 protease inhibition with IC50 values of 15 nM and 64 nM, respectively.

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.

Binding site asymmetry in human transthyretin: insights from a joint neutron and X-ray crystallographic analysis using perdeuterated protein

Human transthyretin has an intrinsic tendency to form amyloid fibrils and is heavily implicated in senile systemic amyloidosis. Here, detailed neutron structural studies of perdeuterated transthyretin are described. The analyses, which fully exploit the enhanced visibility of isotopically replaced hydrogen atoms, yield new information on the stability of the protein and the possible mechanisms of amyloid formation. Residue Ser117 may play a pivotal role in that a single water molecule is closely associated with the γ-hydrogen atoms in one of the binding pockets, and could be important in determining which of the two sites is available to the substrate. The hydrogen-bond network at the monomer-monomer interface is more extensive than that at the dimer-dimer interface. Additionally, the edge strands of the primary dimer are seen to be favourable for continuation of the β-sheet and the formation of an extended cross-β structure through sequential dimer couplings. It is argued that the precursor to fibril formation is the dimeric form of the protein.