"To Be or Not to Be" Protonated: Atomic Details of Human Carbonic Anhydrase-Clinical Drug Complexes by Neutron Crystallography and Simulation

Multiple Binding Modes of Inhibitors to Human Carbonic Anhydrases: An Update on the Design of Isoform-Specific Modulators of Activity

Human carbonic anhydrases (hCAs) are widespread zinc enzymes that catalyze the hydration of CO2 to bicarbonate and a proton. Currently, 15 isoforms have been identified, of which only 12 are catalytically active. Given their involvement in numerous physiological and pathological processes, hCAs are recognized therapeutic targets for the development of inhibitors with biomedical applications. However, despite massive development efforts, very few of the presently available hCA inhibitors show selectivity for a specific isoform. X-ray crystallography is a very useful tool for the rational drug design of enzyme inhibitors. In 2012 we published in Chemical Reviews a highly cited review on hCA family (Alterio, V. et al. Chem Rev. 2012, 112, 4421-4468), analyzing about 300 crystallographic structures of hCA/inhibitor complexes and describing the different CA inhibition mechanisms existing up to that date. However, in the period 2012-2023, almost 700 new hCA/inhibitor complex structures have been deposited in the PDB and a large number of new inhibitor classes have been discovered. Based on these considerations, the aim of this Review is to give a comprehensive update of the structural aspects of hCA/inhibitor interactions covering the period 2012-2023 and to recapitulate how this information can be used for the rational design of more selective versions of such inhibitors.

Targeted anticancer pre-vinylsulfone covalent inhibitors of carbonic anhydrase IX

We designed novel pre-drug compounds that transform into an active form that covalently modifies particular His residue in the active site, a difficult task to achieve, and applied to carbonic anhydrase (CAIX), a transmembrane protein, highly overexpressed in hypoxic solid tumors, important for cancer cell survival and proliferation because it acidifies tumor microenvironment helping invasion and metastases processes. The designed compounds have several functionalities: (1) primary sulfonamide group recognizing carbonic anhydrases (CA), (2) high-affinity moieties specifically recognizing CAIX among all CA isozymes, and (3) forming a covalent bond with the His64 residue. Such targeted covalent compounds possess both high initial affinity and selectivity for the disease target protein followed by complete irreversible inactivation of the protein via covalent modification. Our designed prodrug candidates bearing moderately active pre-vinylsulfone esters or weakly active carbamates optimized for mild covalent modification activity to avoid toxic non-specific modifications and selectively target CAIX. The lead inhibitors reached 2 pM affinity, the highest among known CAIX inhibitors. The strategy could be used for any disease drug target protein bearing a His residue in the vicinity of the active site.

Acetazolamide and human carbonic anhydrases: retrospect, review and discussion of an intimate relationship

Acetazolamide (AZM) is a strong pharmacological sulphonamide-type (R-SO2-NH2, pKa 7.2) inhibitor of the activity of several carbonic anhydrase (CA) isoforms, notably of renal CA II (Ki, 12 nM) and CA IV (Ki, 74 nM). AZM is clinically used for about eighty years in various diseases including epilepsy and glaucoma. Pharmacological AZM increases temporarily the urinary excretion of bicarbonate (HCO3-) and sodium ions (Na+) and sustainably the urinary pH. AZM is excreted almost unchanged over several hours at high rates in the urine. Closely parallel concentrations of circulating and excretory AZM are observed upon administration of therapeutical doses of AZM. In a proof-of-principle study, we investigated the effects of the ingestion of a 250-mg AZM-containing tablet by a healthy volunteer on the urinary excretion of organic and inorganic substances over 5 h (range, 0, 0.5, 1, 1.5, 2, 3, 4, 5 h). Measured analytes included: AZM, amino acids and their metabolites such as guanidinoacetate, i.e. the precursor of creatine, of asymmetrically (ADMA) and symmetrically (SDMA) dimethylated arginine, nitrite (O = N-O-, pKa 3.4) and nitrate (O2N-O-, pKa -1.37), the major metabolites of nitric oxide (NO), the C-H acidic malondialdehyde (MDA; (CHO)2CH2, pKa 4.5), and creatinine for correction of analytes excretion. All analytes were measured by validated isotopologues using gas chromatography-mass spectrometry (GC-MS) methods. AZM excretion in the urine reached its maximum value after 2 h and was fairly stable for the next 3 h. Time series analysis by the ARIMA method was performed. AZM ingestion increased temporarily the urinary excretion of the amino acids Leu + Ile, nitrite and nitrate, decreased temporarily the urinary excretion of other amino acids. AZM decreased sustainably the urinary excretion of MDA, a biomarker of oxidative stress (i.e. lipid peroxidation). Whether this decrease is due to inhibition of the excretion of MDA or attenuation of oxidative stress by AZM is unknown. The acute and chronic effects of AZM on the urinary excretion of electrolytes and physiological substances reported in the literature are discussed in depth in the light of its extraordinary pharmacokinetics and pharmacodynamics. Tolerance development/drug resistance to AZM in chronic use and potential mechanisms are also addressed.

Benchmarking biomolecular force field-based Zn2+ for mono- and bimetallic ligand binding sites

Zn²⁺ is one of the most versatile biologically available metal ions, but accurate modeling of Zn²⁺ -containing metalloproteins at the biomolecular force field level can be challenging. Since most Zn²⁺ models are parameterized in bulk solvent, in-depth knowledge about their performance in a protein environment is limited. Thus, we systematically investigate here the behavior of non-polarizable Zn²⁺ models for their ability to reproduce experimentally determined metal coordination and ligand binding in metalloproteins. The benchmarking is performed in challenging environments, including mono- (carbonic anhydrase II) and bimetallic (metallo-β-lactamase VIM-2) ligand binding sites. We identify key differences in the performance between the Zn²⁺ models with regard to the preferred ligating atoms (charged/non-charged), attraction of water molecules, and the preferred coordination geometry. Based on these results, we suggest suitable simulation conditions for varying Zn²⁺ site geometries that could guide the further development of biomolecular Zn²⁺ models.

Rapid Evolution of a Fragment-like Molecule to Pan-Metallo-Beta-Lactamase Inhibitors: Initial Leads toward Clinical Candidates

With the emergence and rapid spreading of NDM-1 and existence of clinically relevant VIM-1 and IMP-1, discovery of pan inhibitors targeting metallo-beta-lactamases (MBLs) became critical in our battle against bacterial infection. Concurrent with our fragment and high-throughput screenings, we performed a knowledge-based search of known metallo-beta-lactamase inhibitors (MBLIs) to identify starting points for early engagement of medicinal chemistry. A class of compounds exemplified by 11, discovered earlier as B. fragilis metallo-beta-lactamase inhibitors, was selected for in silico virtual screening. From these efforts, compound 12 was identified with activity against NDM-1 only. Initial exploration on metal binding design followed by structure-guided optimization led to the discovery of a series of compounds represented by 23 with a pan MBL inhibition profile. In in vivo studies, compound 23 in combination with imipenem (IPM) robustly lowered the bacterial burden in a murine infection model and became the lead for the invention of MBLI clinical candidates.

Toward the Discovery of a Novel Class of Leads for High Altitude Disorders by Virtual Screening and Molecular Dynamics Approaches Targeting Carbonic Anhydrase

For decades, carbonic anhydrase (CA) inhibitors, most notably the acetazolamide-bearing 1,3,4-thiadiazole moiety, have been exploited at high altitudes to alleviate acute mountain sickness, a syndrome of symptomatic sensitivity to the altitude characterized by nausea, lethargy, headache, anorexia, and inadequate sleep. Therefore, inhibition of CA may be a promising therapeutic strategy for high-altitude disorders. In this study, co-crystallized inhibitors with 1,3,4-thiadiazole, 1,3-benzothiazole, and 1,2,5-oxadiazole scaffolds were employed for pharmacophore-based virtual screening of the ZINC database, followed by molecular docking and molecular dynamics simulation studies against CA to find possible ligands that may emerge as promising inhibitors. Compared to the co-crystal ligands of PDB-1YDB, 6BCC, and 6IC2, ZINC12336992, ZINC24751284, and ZINC58324738 had the highest docking scores of -9.0, -9.0, and -8.9 kcal/mol, respectively. A molecular dynamics (MD) simulation analysis of 100 ns was conducted to verify the interactions of the top-scoring molecules with CA. The system's backbone revealed minor fluctuations, indicating that the CA-ligand complex was stable during the simulation period. Simulated trajectories were used for the MM-GBSA analysis, showing free binding energies of -16.00 ± 0.19, -21.04 ± 0.17, and -19.70 ± 0.18 kcal/mol, respectively. In addition, study of the frontier molecular orbitals of these compounds by DFT-based optimization at the level of B3LYP and the 6-311G(d,p) basis set showed negative values of the HOMO and LUMO, indicating that the ligands are energetically stable, which is essential for forming a stable ligand-protein complex. These molecules may prove to be a promising therapy for high-altitude disorders, necessitating further investigations.

Beta and Gamma Amino Acid-Substituted Benzenesulfonamides as Inhibitors of Human Carbonic Anhydrases

A series of novel benzenesulfonamide derivatives were synthesized bearing para-N β,γ-amino acid or para-N β-amino acid and thiazole moieties and their binding to the human carbonic anhydrase (CA) isozymes determined. These enzymes are involved in various illnesses, such as glaucoma, altitude sickness, epilepsy, obesity, and even cancer. There are numerous compounds that are inhibitors of CA and used as pharmaceuticals. However, most of them bind to most CA isozymes with little selectivity. The design of high affinity and selectivity towards one CA isozyme remains a significant challenge. The beta and gamma amino acid-substituted compound affinities were determined by the fluorescent thermal shift assay and isothermal titration calorimetry for all 12 catalytically active human carbonic anhydrase isozymes, showing the full affinity and selectivity profile. The structures of several compounds were determined by X-ray crystallography, and the binding mode in the active site of CA enzyme was shown.

In the footsteps of an inhibitor unbinding from the active site of human carbonic anhydrase II

The crystal structure of human carbonic anhydrase (HCA) II bound to an inhibitor molecule, 6-hydroxy-2-thioxocoumarin (FC5), shows FC5 to be located in a hydrophobic pocket at the active site. The present work employs classical molecular dynamics (MD) simulation to follow the FC5 molecule for 1 μs as it unbinds from its binding location, adopts the path of substrate/product diffusion (path 1) to leave the active site at around 75 ns. It is then found to undergo repeated binding and unbinding at different locations on the surface of the enzyme in water. Several transient excursions through different regions of the enzyme are also observed prior to its exit from the active site. These transient paths are combined with functionally relevant cavities/channels to enlist five additional pathways (path 2-6). Pathways 1-6 are subsequently explored using steered MD and umbrella sampling simulations. A free energy barrier of 0.969 kcal mol^-1 is encountered along path 1, while barriers in the range of 0.57-2.84 kcal mol^-1 are obtained along paths 2, 4 and 5. We also analyze in detail the interaction between FC5 and the enzyme along each path as the former leaves the active site of HCA II. Our results indicate path 1 to be the major exit pathway for FC5, although competing contributions may also come from the paths 2, 4 and 5.Communicated by Ramaswamy H. Sarma.

Structure and mechanism of secondary sulfonamide binding to carbonic anhydrases

Zinc-containing metalloenzyme carbonic anhydrase (CA) binds primary sulfonamides with extremely high, up to picomolar, affinity by forming a coordination bond between the negatively charged amino group and the zinc ion and making hydrogen bonds and hydrophobic contacts with other parts of the inhibitor molecule. However, N-methyl-substituted, secondary or tertiary sulfonamides bind CA with much lower affinity. In search for an explanation for this diminished affinity, a series of secondary sulfonamides were synthesized and, together with analogous primary sulfonamides, the affinities for 12 recombinant catalytically active human CA isoforms were determined by the fluorescent thermal shift assay, stopped-flow assay of the inhibition of enzymatic activity and isothermal titration calorimetry. The binding profile of secondary sulfonamides as a function of pH showed the same U-shape dependence seen for primary sulfonamides. This dependence demonstrated that there were protein binding-linked protonation reactions that should be dissected for the estimation of the intrinsic binding constants to perform structure-thermodynamics analysis. X-ray crystallographic structures of secondary sulfonamides and computational modeling dissected the atomic contributions to the binding energetics. Secondary sulfonamides bind to carbonic anhydrases via coordination bond between the negatively charged nitrogen of alkylated amino group and Zn(II) in the active site of CA. The binding reaction is linked to deprotonation of the amino group and protonation of the Zn(II)-bound hydroxide. To perform the structure-thermodynamics analysis, contributions of these linked reactions must be subtracted to determine the intrinsic energetics. In this aspect, the secondary sulfonamides are similar to primary sulfonamides as CA inhibitors.

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.

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.

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.

Isoform-Selective Enzyme Inhibitors by Exploring Pocket Size According to the Lock-and-Key Principle

In the design of high-affinity and enzyme isoform-selective inhibitors, we applied an approach of augmenting the substituents attached to the benzenesulfonamide scaffold in three ways, namely, substitutions at the 3,5- or 2,4,6-positions or expansion of the condensed ring system. The increased size of the substituents determined the spatial limitations of the active sites of the 12 catalytically active human carbonic anhydrase (CA) isoforms until no binding was observed because of the inability of the compounds to fit in the active site. This approach led to the discovery of high-affinity and high-selectivity compounds for the anticancer target CA IX and antiobesity target CA VB. The x-ray crystallographic structures of compounds bound to CA IX showed the positions of the bound compounds, whereas computational modeling confirmed that steric clashes prevent the binding of these compounds to other isoforms and thus avoid undesired side effects. Such an approach, based on the Lock-and-Key principle, could be used for the development of enzyme-specific drug candidate compounds.

Native mass spectrometry of human carbonic anhydrase I and its inhibitor complexes

Abstract

Native mass spectrometry is a potent technique to study and characterize biomacromolecules in their native state. Here, we have applied this method to explore the solution chemistry of human carbonic anhydrase I (hCA I) and its interactions with four different inhibitors, namely three sulfonamide inhibitors (AAZ, MZA, SLC-0111) and the dithiocarbamate derivative of morpholine (DTC). Through high-resolution ESI-Q-TOF measurements, the native state of hCA I and the binding of the above inhibitors were characterized in the molecular detail. Native mass spectrometry was also exploited to assess the direct competition in solution among the various inhibitors in relation to their affinity constants. Additional studies were conducted on the interaction of hCA I with the metallodrug auranofin, under various solution and instrumental conditions. Auranofin is a selective reagent for solvent-accessible free cysteine residues, and its reactivity was analyzed also in the presence of CA inhibitors. Overall, our investigation reveals that native mass spectrometry represents an excellent tool to characterize the solution behavior of carbonic anhydrase.

Graphic abstract

Electronic supplementary material

The online version of this article (10.1007/s00775-020-01818-8) contains supplementary material, which is available to authorized users.

Visualizing drug binding interactions using microcrystal electron diffraction

Visualizing ligand binding interactions is important for structure-based drug design and fragment-based screening methods. Rapid and uniform soaking with potentially reduced lattice defects make small macromolecular crystals attractive targets for studying drug binding using microcrystal electron diffraction (MicroED). However, so far no drug binding interactions could unambiguously be resolved by electron diffraction alone. Here, we use MicroED to study the binding of a sulfonamide inhibitor to human carbonic anhydrase isoform II (HCA II). We show that MicroED data can efficiently be collected on a conventional transmission electron microscope from thin hydrated microcrystals soaked with the clinical drug acetazolamide (AZM). The data are of high enough quality to unequivocally fit and resolve the bound inhibitor. We anticipate MicroED can play an important role in facilitating in-house fragment screening for drug discovery, complementing existing methods in structural biology such as X-ray and neutron diffraction.

Clabbers et al. utilize MicroED to present the structure of both apo and inhibitor-bound human carbonic anhydrase II at a high resolution to clearly identify the interaction of the inhibitor, acetazolamide. This method eases the difficulty of both crystallizing the protein and soaking the inhibitor in a smaller protein crystal.

Sulfamoyl Heteroarylcarboxylic Acids as Promising Metallo-β-Lactamase Inhibitors for Controlling Bacterial Carbapenem Resistance

Carbapenem antibiotics are the last resort for control of severe infectious diseases, bloodstream infections, and pneumonia caused by Gram-negative bacteria, including Enterobacteriaceae. However, carbapenem-resistant Enterobacteriaceae (CRE) strains have spread globally and are a critical concern in clinical settings because CRE infections are recognized as a leading cause of increased mortality among hospitalized patients. Most CRE produce certain kinds of serine carbapenemases (e.g., KPC- and GES-type β-lactamases) or metallo-β-lactamases (MBLs), which can hydrolyze carbapenems. Although effective MBL inhibitors are expected to restore carbapenem efficacy against MBL-producing CRE, no MBL inhibitor is currently clinically available. Here, we synthesized 2,5-diethyl-1-methyl-4-sulfamoylpyrrole-3-carboxylic acid (SPC), which is a potent inhibitor of MBLs. SPC is a remarkable lead compound for clinically useful MBL inhibitors and can potentially provide a considerable benefit to patients receiving treatment for lethal infectious diseases caused by MBL-producing CRE.

Production of metallo-β-lactamases (MBLs), which hydrolyze carbapenems, is a cause of carbapenem resistance in Enterobacteriaceae. Development of effective inhibitors for MBLs is one approach to restore carbapenem efficacy in carbapenem-resistant Enterobacteriaceae (CRE). We report here that sulfamoyl heteroarylcarboxylic acids (SHCs) can competitively inhibit the globally spreading and clinically relevant MBLs (i.e., IMP-, NDM-, and VIM-type MBLs) at nanomolar to micromolar orders of magnitude. Addition of SHCs restored meropenem efficacy against 17/19 IMP-type and 7/14 NDM-type MBL-producing Enterobacteriaceae to satisfactory clinical levels. SHCs were also effective against IMP-type MBL-producing Acinetobacter spp. and engineered Escherichia coli strains overproducing individual minor MBLs (i.e., TMB-2, SPM-1, DIM-1, SIM-1, and KHM-1). However, SHCs were less effective against MBL-producing Pseudomonas aeruginosa. Combination therapy with meropenem and SHCs successfully cured mice infected with IMP-1-producing E. coli and dually NDM-1/VIM-1-producing Klebsiella pneumoniae clinical isolates. X-ray crystallographic analyses revealed the inhibition mode of SHCs against MBLs; the sulfamoyl group of SHCs coordinated to two zinc ions, and the carboxylate group coordinated to one zinc ion and bound to positively charged amino acids Lys224/Arg228 conserved in MBLs. Preclinical testing revealed that the SHCs showed low toxicity in cell lines and mice and high stability in human liver microsomes. Our results indicate that SHCs are promising lead compounds for inhibitors of MBLs to combat MBL-producing CRE.

Protein X-ray Crystallography and Drug Discovery

With the advent of structural biology in the drug discovery process, medicinal chemists gained the opportunity to use detailed structural information in order to progress screening hits into leads or drug candidates. X-ray crystallography has proven to be an invaluable tool in this respect, as it is able to provide exquisitely comprehensive structural information about the interaction of a ligand with a pharmacological target. As fragment-based drug discovery emerged in the recent years, X-ray crystallography has also become a powerful screening technology, able to provide structural information on complexes involving low-molecular weight compounds, despite weak binding affinities. Given the low numbers of compounds needed in a fragment library, compared to the hundreds of thousand usually present in drug-like compound libraries, it now becomes feasible to screen a whole fragment library using X-ray crystallography, providing a wealth of structural details that will fuel the fragment to drug process. Here, we review theoretical and practical aspects as well as the pros and cons of using X-ray crystallography in the drug discovery process.

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

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

Large crystal growth for neutron protein crystallography

The use of neutron protein crystallography (NPX) is expanding rapidly, with most structures determined in the last decade. This growth is stimulated by a number of developments, spanning from the building of new NPX beamlines to the availability of improved software for structure refinement. The main bottleneck preventing structural biologists from adding NPX to the suite of methods commonly used is the large volume of the individual crystals required for a successful experiment. A survey of deposited NPX structures in the Protein Data Bank shows that about two-thirds came from crystals prepared using vapor diffusion, while batch and dialysis-based methods all-together contribute to most of the remaining one-third. This chapter explains the underlying principles of these protein crystallization methods and provides practical examples that may help others to successfully prepare large crystals for NPX.

Structural comparison of protiated, H/D-exchanged and deuterated human carbonic anhydrase IX

Human carbonic anhydrase IX (CA IX) expression is upregulated in hypoxic solid tumours, promoting cell survival and metastasis. This observation has made CA IX a target for the development of CA isoform-selective inhibitors. To enable structural studies of CA IX-inhibitor complexes using X-ray and neutron crystallography, a CA IX surface variant (CA IXSV; the catalytic domain with six surface amino-acid substitutions) has been developed that can be routinely crystallized. Here, the preparation of protiated (H/H), H/D-exchanged (H/D) and deuterated (D/D) CA IXSV for crystallographic studies and their structural comparison are described. Four CA IXSV X-ray crystal structures are compared: two H/H crystal forms, an H/D crystal form and a D/D crystal form. The overall active-site organization in each version is essentially the same, with only minor positional changes in active-site solvent, which may be owing to deuteration and/or resolution differences. Analysis of the crystal contacts and packing reveals different arrangements of CA IXSV compared with previous reports. To our knowledge, this is the first report comparing three different deuterium-labelled crystal structures of the same protein, marking an important step in validating the active-site structure of CA IXSV for neutron protein crystallography.

Engineered Carbonic Anhydrase VI-Mimic Enzyme Switched the Structure and Affinities of Inhibitors

Secretory human carbonic anhydrase VI (CA VI) has emerged as a potential drug target due to its role in pathological states, such as excess acidity-caused dental caries and injuries of gastric epithelium. Currently, there are no available CA VI-selective inhibitors or crystallographic structures of inhibitors bound to CA VI. The present study focuses on the site-directed CA II mutant mimicking the active site of CA VI for inhibitor screening. The interactions between CA VI-mimic and a series of benzenesulfonamides were evaluated by fluorescent thermal shift assay, stopped-flow CO₂ hydration assay, isothermal titration calorimetry, and X-ray crystallography. Kinetic parameters showed that A65T, N67Q, F130Y, V134Q, L203T mutations did not influence catalytic properties of CA II, but inhibitor affinities resembled CA VI, exhibiting up to 0.16 nM intrinsic affinity for CA VI-mimic. Structurally, binding site of CA VI-mimic was found to be similar to CA VI. The ligand interactions with mutated side chains observed in three crystallographic structures allowed to rationalize observed variation of binding modes and experimental binding affinities to CA VI. This integrative set of kinetic, thermodynamic, and structural data revealed CA VI-mimic as a useful model to design CA VI-specific inhibitors which could be beneficial for novel therapeutic applications.

The upcoming subatomic resolution revolution

Subatomic resolution macromolecular crystallography has been revealing the most fascinating details of macromolecular structures for many years. This most extreme form of macromolecular crystallography is going through rapid changes. A new generation of superbrilliant X-ray sources and detectors is facilitating the rapid acquisition of high-quality datasets. Equally important, a new breed of methods and highly integrated advanced computational tools for structure refinement and analysis is poised to change the way we use subatomic resolution data and reposition high-resolution macromolecular crystallography in medicinal chemistry studies. Subatomic resolution macromolecular crystallography may soon be a routine source of detailed molecular information besides precise geometries, including binding energies and other chemical descriptors, opening new possibilities of application.

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.

Neutron scattering in the biological sciences: progress and prospects

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

Temperature-Induced Replacement of Phosphate Proton with Metal Ion Captured in Neutron Structures of A-DNA

Nucleic acids can fold into well-defined 3D structures that help determine their function. Knowing precise nucleic acid structures can also be used for the design of nucleic acid-based therapeutics. However, locations of hydrogen atoms, which are key players of nucleic acid function, are normally not determined with X-ray crystallography. Accurate determination of hydrogen atom positions can provide indispensable information on protonation states, hydrogen bonding, and water architecture in nucleic acids. Here, we used neutron crystallography in combination with X-ray diffraction to obtain joint X-ray/neutron structures at both room and cryo temperatures of a self-complementary A-DNA oligonucleotide d[GTGG(CSe)CAC]2 containing 2'-SeCH3 modification on Cyt5 (CSe) at pH 5.6. We directly observed protonation of a backbone phosphate oxygen of Ade7 at room temperature. The proton is replaced with hydrated Mg2+ upon cooling the crystal to 100 K, indicating that metal binding is favored at low temperature, whereas proton binding is dominant at room temperature.

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

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

Biophysical, Biochemical, and Cell Based Approaches Used to Decipher the Role of Carbonic Anhydrases in Cancer and to Evaluate the Potency of Targeted Inhibitors

Carbonic anhydrases (CAs) are thought to be important for regulating pH in the tumor microenvironment. A few of the CA isoforms are upregulated in cancer cells, with only limited expression in normal cells. For these reasons, there is interest in developing inhibitors that target these tumor-associated CA isoforms, with increased efficacy but limited nonspecific cytotoxicity. Here we present some of the biophysical, biochemical, and cell based techniques and approaches that can be used to evaluate the potency of CA targeted inhibitors and decipher the role of CAs in tumorigenesis, cancer progression, and metastatic processes. These techniques include esterase activity assays, stop flow kinetics, and mass inlet mass spectroscopy (MIMS), all of which measure enzymatic activity of purified protein, in the presence or absence of inhibitors. Also discussed is the application of X-ray crystallography and Cryo-EM as well as other structure-based techniques and thermal shift assays to the studies of CA structure and function. Further, large-scale genomic and proteomic analytical methods, as well as cell based techniques like those that measure cell growth, apoptosis, clonogenicity, and cell migration and invasion, are discussed. We conclude by reviewing approaches that test the metastatic potential of CAs and how the aforementioned techniques have contributed to the field of CA cancer research.

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.