Building reactive copper centers in human carbonic anhydrase II

Machine Learning-Enhanced Quantum Chemistry-Assisted Refinement of the Active Site Structure of Metalloproteins

Understanding the fine structural details of inhibitor binding at the active site of metalloenzymes can have a profound impact on the rational drug design targeted to this broad class of biomolecules. Structural techniques such as NMR, cryo-EM, and X-ray crystallography can provide bond lengths and angles, but the uncertainties in these measurements can be as large as the range of values that have been observed for these quantities in all the published structures. This uncertainty is far too large to allow for reliable calculations at the quantum chemical (QC) levels for developing precise structure-activity relationships or for improving the energetic considerations in protein-inhibitor studies. Therefore, the need arises to rely upon computational methods to refine the active site structures well beyond the resolution obtained with routine application of structural methods. In a recent paper, we have shown that it is possible to refine the active site of cobalt(II)-substituted MMP12, a metalloprotein that is a relevant drug target, by matching to the experimental pseudocontact shifts (PCS) those calculated using multireference ab initio QC methods. The computational cost of this methodology becomes a significant bottleneck when the starting structure is not sufficiently close to the final one, which is often the case with biomolecular structures. To tackle this problem, we have developed an approach based on a neural network (NN) and a support vector regression (SVR) and applied it to the refinement of the active site structure of oxalate-inhibited human carbonic anhydrase 2 (hCAII), another prototypical metalloprotein target. The refined structure gives a remarkably good agreement between the QC-calculated and the experimental PCS. This study not only contributes to the knowledge of CAII but also demonstrates the utility of combining machine learning (ML) algorithms with QC calculations, offering a promising avenue for investigating other drug targets and complex biological systems in general.

Elucidating the concentration-dependent effects of thiocyanate binding to carbonic anhydrase

Many proteins naturally carry metal centers, with a large share of them being in the active sites of several enzymes. Paramagnetic effects are a powerful source of structural information and, therefore, if the native metal is paramagnetic, or it can be functionally substituted with a paramagnetic one, paramagnetic effects can be used to study the metal sites, as well as the overall structure of the protein. One notable example is cobalt(II) substitution for zinc(II) in carbonic anhydrase. In this manuscript we investigate the effects of sodium thiocyanate on the chemical environment of the metal ion of the human carbonic anhydrase II. The electron paramagnetic resonance (EPR) titration of the cobalt(II) protein with thiocyanate shows that the EPR spectrum changes from A-type to C-type on passing from 1:1 to 1:1000-fold ligand excess. This indicates the occurrence of a change in the electronic structure, which may reflect a sizable change in the metal coordination environment in turn caused by a modification of the frozen solvent glass. However, paramagnetic nuclear magnetic resonance (NMR) data indicate that the metal coordination cage remains unperturbed even in 1:1000-fold ligand excess. This result proves that the C-type EPR spectrum observed at large ligand concentration should be ascribed to the low temperature at which EPR measurements are performed, which impacts on the structure of the protein when it is destabilized by a high concentration of a chaotropic agent.

Metal Ion Binding Induces Local Protein Unfolding and Destabilizes Human Carbonic Anhydrase II

Human carbonic anhydrase II (HCA) is a robust metalloprotein and an excellent biological model system to study the thermodynamics of metal ion coordination. Apo-HCA binds one zinc ion or two copper ions. We studied these binding processes at five temperatures (15-35 °C) using isothermal titration calorimetry, yielding thermodynamic parameters corrected for pH and buffer effects. We then sought to identify binding-induced structural changes. Our data suggest that binding at the active site organizes 6-8 residues; however, copper binding near the N-terminus results in a net unfolding of 6-7 residues. This surprising destabilization was confirmed using circular dichroism and protein stability measurements. Metal binding induced unfolding may represent an important regulatory mechanism, but it may be easily missed by NMR and X-ray crystallography. Thus, in addition to highlighting a highly novel binding-induced unfolding event, we demonstrate the value of calorimetry for studying the structural implications of metal binding.

Isothermal titration calorimetry and surface plasmon resonance analysis using the dynamic approach

Biophysical techniques such as isothermal titration calorimetry (ITC) and surface plasmon resonance (SPR) are routinely used to ascertain the global binding mechanisms of protein-protein or protein-ligand interaction. Recently, Dumas etal, have explicitly modelled the instrument response of the ligand dilution and analysed the ITC thermogram to obtain kinetic rate constants. Adopting a similar approach, we have integrated the dynamic instrument response with the binding mechanism to simulate the ITC profiles of equivalent and independent binding sites, equivalent and sequential binding sites and aggregating systems. The results were benchmarked against the standard commercial software Origin-ITC. Further, the experimental ITC chromatograms of 2′-CMP + RNASE and BH3I-1 + hBCL_XL interactions were analysed and shown to be comparable with that of the conventional analysis. Dynamic approach was applied to simulate the SPR profiles of a two-state model, and could reproduce the experimental profile accurately.

•

Incorporated instrument response within the kinetic framework using dynamic approach to analyse ITC and SPR data.

•

Different modelling approaches for instrument response such as lumped and kinetic modelling were compared and their equivalence were shown.

•

(1) equivalent single site, (2) equivalent sequential sites, (3) equivalent parallel sites and (4) aggregating molecular system were modelled using dynamic approach.

Carbonic anhydrase modification for carbon management

Carbonic anhydrase modification (chemical and biological) is an attractive strategy for its diverse application to accelerate the absorption of CO₂ from a flue gas with improved activity and stability. This article reports various possibilities of CA modification using metal-ligand homologous chemistry, cross-linking agents, and residue- and group-specific and genetic modifications, and assesses their role in carbon management. Chemically modified carbonic anhydrase is able to improve the absorption of carbon dioxide from a gas stream into mediation compounds with enhanced sequestration and mineral formation. Genetically modified CA polypeptide can also increase carbon dioxide conversion. Chemical modification of CA can be categorized in terms of (i) residue-specific modification (involves protein-ligand interaction in terms of substitution/addition) and group-specific modifications (based on the functional groups of the target CA). For every sustainable change, there should be no/limited toxic or immunological response. In this review, several CA modification pathways and biocompatibility rules are proposed as a theoretical support for emerging research in this area.

Biophysical approaches for the study of metal-protein interactions

Protein-protein interactions play important roles for a variety of cell functions, often involving metal ions; in fact, metal-ion binding mediates and regulates the activity of a wide range of biomolecules. Enlightening all of the specific features of metal-protein and metal-mediated protein-protein interactions can be a very challenging task; a detailed knowledge of the thermodynamic and spectroscopic parameters and the structural changes of the protein is normally required. For this purpose, many experimental techniques are employed, embracing all fields of Analytical and Bioinorganic Chemistry. In addition, the use of peptide models, reproducing the primary sequence of the metal-binding sites, is also proved to be useful. In this paper, a review of the most useful techniques for studying ligand-protein interactions with a special emphasis on metal-protein interactions is provided, with a critical summary of their strengths and limitations.

How Do Nuclei Couple to the Magnetic Moment of a Paramagnetic Center? A New Theory at the Gauntlet of the Experiments

The recent derivation, based on pure quantum chemistry (QC) first-principles, of the pseudocontact shifts (PCSs) caused by a paramagnetic metal center on far away nuclei has cast doubts on the validity of the semiempirical (SE) theory, predicting PCSs to arise from the metal magnetic susceptibility anisotropy. The SE theory has been used and applied countless times, especially in the last 2 decades, to obtain structural information on proteins containing paramagnetic metal ions. We show here that the QC and SE predictions can be directly tested against experiments, provided a suitable macromolecular system is used. The SE approach yields a good prediction of the experimental PCSs while the QC one does not. It appears that the classic theory is able to grasp satisfactorily the underlying physics.

The copper(II) binding centres of carbonic anhydrase are differently affected by reductants that ensure the redox intracellular environment

Copper is involved in several biological processes. The static and labile copper pools are controlled by means of a network of influx and efflux transporters, storage proteins, chaperones, transcription factors and small molecules as glutathione (GSH), which contributes to the cell reducing environment. To follow the fate of intracellular copper labile pool, a variant of human apocarbonic anhydrase has been proposed as fluorescent probe to monitor cytoplasmic Cu²⁺. Aware that in this cellular compartment copper ion is present as Cu⁺, electron spin resonance technique (ESR) was used to ascertain whether (bovine or human) carbonic anhydrase (CA) was able to accommodate Cu⁺ in the same sites occupied by Cu²⁺, in the presence of naturally occurring reducing agents such as ascorbate and GSH. Our ESR results on Cu²⁺ complexes with CA allow for a complete characterization of the two metal binding sites of the protein in solution. The use of the reported affinity constants of zinc in the catalytic site and of Cu²⁺ in the peripheral and catalytic site, allow us to obtain the speciation of copper species mimicking the spectroscopic study conditions. The different Cu²⁺ coordination features in the catalytic and the peripheral (the N-terminus cleft mouth) binding sites influence the chemical reduction effect of the two main naturally occurring reductants. Ascorbate reversibly reduces the Cu²⁺ complex with CA, while glutathione irreversibly induces the formation of Cu²⁺ complex with its oxidized form (GSSG). Our results questioned the use of CA as intracellular Cu²⁺ sensor. Furthermore, translating these findings to intracellular environment, the conversion of GSH in GSSG can significantly alter the metallostasis.

Dioxygen, an unexpected carbonic anhydrase ligand

Carbonic anhydrases (CAs, EC 4.2.1.1) are ubiquitous metalloenzymes, grouped into seven different classes, which catalyze the reaction of CO₂ hydration to bicarbonate and protons. All of the fifteen human isoforms reported to date belong to the α-class and contain zinc as a cofactor. The structure of human Zn,Cu-CA II has been solved which contains a copper ion bound at its N-terminal, coordinated to His4 and His64. In the active site a dioxygen molecule is coordinated to the zinc ion. Since dioxygen is a rather unexpected CA ligand, molecular dynamics (MD) simulations were performed which suggested a superoxide character of the zinc bound O₂.

The Complex Relationship between Metals and Carbonic Anhydrase: New Insights and Perspectives

Carbonic anhydrase is a ubiquitous metalloenzyme, which catalyzes the reversible hydration of CO₂ to HCO₃⁻ and H⁺. Metals play a key role in the bioactivity of this metalloenzyme, although their relationships with CA have not been completely clarified to date. The aim of this review is to explore the complexity and multi-aspect nature of these relationships, since metals can be cofactors of CA, but also inhibitors of CA activity and modulators of CA expression. Moreover, this work analyzes new insights and perspectives that allow translating new advances in basic science on the interaction between CA and metals to applications in several fields of research, ranging from biotechnology to environmental sciences.

Isothermal Titration Calorimetry Measurements of Metal Ions Binding to Proteins

ITC measurements involving metal ions are susceptible to a number of competing reactions (oxidation, precipitation, and hydrolysis) and coupled reactions involving the buffer and protons. Stabilization and delivery of the metal ion as a well-defined and well-characterized complex with the buffer, or a specific ligand, can suppress undesired solution chemistry and, depending on the stability of the metal complex, allow accurate measurements of higher affinity protein-binding sites. This requires, however, knowledge of the thermodynamics of formation of the metal complex and accounting for its contribution to the experimentally measured values (KITC and ΔHITC) through a post hoc analysis that provides the condition-independent binding thermodynamics (K, ΔG(o), ΔH, ΔS, and ΔCP). This analysis also quantifies the number of protons that are displaced when the metal ion binds to the protein.

Characterization of the Copper(II) Binding Sites in Human Carbonic Anhydrase II

Human carbonic anhydrase (CA) is a well-studied, robust, mononuclear Zn-containing metalloprotein that serves as an excellent biological ligand system to study the thermodynamics associated with metal ion coordination chemistry in aqueous solution. The apo form of human carbonic anhydrase II (CA) binds 2 equiv of copper(II) with high affinity. The Cu(2+) ions bind independently forming two noncoupled type II copper centers in CA (CuA and CuB). However, the location and coordination mode of the CuA site in solution is unclear, compared to the CuB site that has been well-characterized. Using paramagnetic NMR techniques and X-ray absorption spectroscopy we identified an N-terminal Cu(2+) binding location and collected information on the coordination mode of the CuA site in CA, which is consistent with a four- to five-coordinate N-terminal Cu(2+) binding site reminiscent to a number of N-terminal copper(II) binding sites including the copper(II)-amino terminal Cu(2+) and Ni(2+) and copper(II)-β-amyloid complexes. Additionally, we report a more detailed analysis of the thermodynamics associated with copper(II) binding to CA. Although we are still unable to fully deconvolute Cu(2+) binding data to the high-affinity CuA site, we derived pH- and buffer-independent values for the thermodynamics parameters K and ΔH associated with Cu(2+) binding to the CuB site of CA to be 2 × 10(9) and -17.4 kcal/mol, respectively.