Carbonic Anhydrases: Different Active Sites, Same Metal Selectivity Rules.

Carbonic anhydrases are mononuclear metalloenzymes catalyzing the reversible hydration of carbon dioxide in organisms belonging to all three domains of life. Although the mechanism of the catalytic reaction is similar, different families of carbonic anhydrases do not have a common ancestor nor do they exhibit significant resemblance in the amino acid sequence or the structure and composition of the metal-binding sites. Little is known about the physical principles determining the metal affinity and selectivity of the catalytic centers, and how well the native metal is protected from being dislodged by other metal species from the local environment. Here, we endeavor to shed light on these issues by studying (via a combination of density functional theory calculations and polarizable continuum model computations) the thermodynamic outcome of the competition between the native metal cation and its noncognate competitor in various metal-binding sites. Typical representatives of the competing cations from the cellular environments of the respective classes of carbonic anhydrases are considered. The calculations reveal how the Gibbs energy of the metal competition changes when varying the metal type, structure, composition, and solvent exposure of the active center. Physical principles governing metal competition in different carbonic anhydrase metal-binding sites are delineated.

In Silico Analysis of the Ga3+/Fe3+ Competition for Binding the Iron-Scavenging Siderophores of P. aeruginosa-Implementation of Three Gallium-Based Complexes in the "Trojan Horse" Antibacterial Strategy.

The emergence of multidrug-resistant (MDR) microorganisms combined with the ever-draining antibiotic pipeline poses a disturbing and immensely growing public health challenge that requires a multidisciplinary approach and the application of novel therapies aimed at unconventional targets and/or applying innovative drug formulations. Hence, bacterial iron acquisition systems and bacterial Fe2+/3+-containing enzymes have been identified as a plausible target of great potential. The intriguing "Trojan horse" approach deprives microorganisms from the essential iron. Recently, gallium's potential in medicine as an iron mimicry species has attracted vast attention. Different Ga3+ formulations exhibit diverse effects upon entering the cell and thus supposedly have multiple targets. The aim of the current study is to specifically distinguish characteristics of great significance in regard to the initial gallium-based complex, allowing the alien cation to effectively compete with the native ferric ion for binding the siderophores pyochelin and pyoverdine secreted by the bacterium P. aeruginosa. Therefore, three gallium-based formulations were taken into consideration: the first-generation gallium nitrate, Ga(NO3)3, metabolized to Ga3+-hydrated forms, the second-generation gallium maltolate (tris(3-hydroxy-2-methyl-4-pyronato)gallium), and the experimentally proven Ga carrier in the bloodstream-the protein transferrin. We employed a reliable in silico approach based on DFT computations in order to understand the underlying biochemical processes that govern the Ga3+/Fe3+ rivalry for binding the two bacterial siderophores.

N-Methyl- and N-Phenylpiperazine Functionalized Styryl Dyes Inside Cucurbiturils: Theoretical Assessment of the Factors Governing the Host-Guest Recognition.

The family of cucurbiturils (CBs), the unique pumpkin-shaped macrocycles, has received great attention over the past four decades owing to their remarkable recognition properties. They have found diverse applications including biosensing and drug delivery technologies. The cucurbituril complexation of guest molecules can modulate their pKas, improve their solubility in aqueous solution, and reduce the adverse effects of the drugs, as well as enhance the stability and/or enable targeted delivery of the drug molecule. Employing twelve cationic styryl dyes with N-methyl- and N-phenylpiperazine functionality as probes, we attempted to understand the factors that govern the host-guest complexation of such molecules within CB[7] and CB[8] host systems. Various key factors determining the process were recognized, such as the pH and dielectric constant of the medium, the cavity size of the host, the chemical characteristics of the substituents in the guest entity, and the presence/absence of metal cations. The presented results add to our understanding (at the molecular level) of the mechanism of encapsulation of styryl dyes by cucurbiturils, thus shedding new light on various aspects of the intriguing complexation chemistry and the underlying recognition processes.

Theoretical Assessment of the Ligand/Metal/Quadruplex Recognition in the Non-Canonical Nucleic Acids Structures.

Quadruplexes (GQs), peculiar DNA/RNA motifs concentrated in specific genomic regions, play a vital role in biological processes including telomere stability and, hence, represent promising targets for anticancer therapy. GQs are formed by folding guanine-rich sequences into square planar G-tetrads which stack onto one another. Metal cations, most often potassium, further stabilize the architecture by coordinating the lone electron pairs of the O atoms. The presence of additional nucleic acid bases, however, has been recently observed experimentally and contributes substantially to the structural heterogeneity of quadruplexes. Therefore, it is of paramount significance to understand the factors governing the underlying complex processes in these structures. The current study employs DFT calculations to model the interactions between metal cations (K+, Na+, Sr2+) and diverse tetrads composed of a guanine layer in combination with a guanine (G)-, adenine (A)-, cytosine (C)-, thymine (T)-, or uracil (U)-based tetrad layer. Moreover, the addition of 4-(3,4-dihydroisoquinolin-2-yl)-2-(quinolin-2-yl)quinazoline to the modeled quadruplexes as a possible mechanism of its well-exerted antitumor effect is assessed. The calculations imply that the metal cation competition and ligand complexation are influenced by the balance between electronic and implicit/explicit solvation effects, the composition of the tetrad layers, as well as by the solvent exposure to the surrounding environment expressed in terms of different dielectric constant values. The provided results significantly enhance our understanding of quadruplex diversity, ligand recognition, and the underlying mechanisms of stabilization at an atomic level.

Electrostatic interactions - key determinants of the metal selectivity in La3+ and Ca2+ binding proteins.

Nearly half of all known proteins contain metal co-factors. In the course of evolution two dozen metal cations (mostly monovalent and divalent species) have been selected to participate in processes of vital importance for living organisms. Trivalent metal cations have also been selected, although to a lesser extent as compared with their mono- and divalent counterparts. Notably, factors governing the metal selectivity in trivalent metal centers in proteins are less well understood than those in the respective divalent metal centers. Thus, the source of high La3+/Ca2+ selectivity in lanthanum-binding proteins, as compared with that of calcium-binding proteins (i.e., calmodulin), is still shrouded in mystery. The well-calibrated thermochemical calculations, performed here, reveal the dominating role of electrostatic interactions in shaping the metal selectivity in La3+-binding centers. The calculations also disclose other (second-order) determinants of metal selectivity in these systems, such as the rigidity and extent of solvent exposure of the binding site. All these factors are also implicated in shaping the metal selectivity in Ca2+-binding proteins.

Cu+/Ag+ Competition in Type I Copper Proteins (T1Cu).

Due to the similarity in the basic coordination behavior of their mono-charged cations, silver biochemistry is known to be linked to that of copper in biological systems. Still, Cu+/2+ is an essential micronutrient in many organisms, while no known biological process requires silver. In human cells, copper regulation and trafficking is strictly controlled by complex systems including many cytosolic copper chaperones, whereas some bacteria exploit the so-called "blue copper" proteins. Therefore, evaluating the controlling factors of the competition between these two metal cations is of enormous interest. By employing the tools of computational chemistry, we aim to delineate the extent to which Ag+ might be able to compete with the endogenous copper in its Type I (T1Cu) proteins, and where and if, alternatively, it is handled uniquely. The effect of the surrounding media (dielectric constant) and the type, number, and composition of amino acid residues are taken into account when modelling the reactions in the present study. The obtained results clearly indicate the susceptibility of the T1Cu proteins to a silver attack due to the favorable composition and geometry of the metal-binding centers, along with the similarity between the Ag+/Cu+-containing structures. Furthermore, by exploring intriguing questions of both metals' coordination chemistry, an important background for understanding the metabolism and biotransformation of silver in organisms is provided.

Lanthanides as Calcium Mimetic Species in Calcium-Signaling/Buffering Proteins: The Effect of Lanthanide Type on the Ca2+/Ln3+ Competition.

Lanthanides, the 14 4f-block elements plus Lanthanum, have been extensively used to study the structure and biochemical properties of metalloproteins. The characteristics of lanthanides within the lanthanide series are similar, but not identical. The present research offers a systematic investigation of the ability of the entire Ln3+ series to substitute for Ca2+ in biological systems. A well-calibrated DFT/PCM protocol is employed in studying the factors that control the metal selectivity in biological systems by modeling typical calcium signaling/buffering binding sites and elucidating the thermodynamic outcome of the competition between the "alien" La3+/Ln3+ and "native" Ca2+, and La3+ - Ln3+ within the lanthanide series. The calculations performed reveal that the major determinant of the Ca2+/Ln3+ selectivity in calcium proteins is the net charge of the calcium binding pocket; the more negative the charge, the higher the competitiveness of the trivalent Ln3+ with respect to its Ca2+ contender. Solvent exposure of the binding site also influences the process; buried active centers with net charge of -4 or -3 are characterized by higher Ln3+ over Ca2+ selectivity, whereas it is the opposite for sites with overall charge of -1. Within the series, the competition between La3+ and its fellow lanthanides is determined by the balance between two competing effects: electronic (favoring heavier lanthanides) and solvation (generally favoring the lighter lanthanides).

Metal-Assisted Complexation of Fluorogenic Dyes by Cucurbit[7]uril and Cucurbit[8]uril: A DFT Evaluation of the Key Factors Governing the Host-Guest Recognition.

With the emergence of host-guest systems, a novel branch of complexation chemistry has found wide application in industries such as food, pharmacy, medicine, environmental protection and cosmetics. Along with the extensively studied cyclodextrins and calixarenes, the innovative cucurbiturils (CB) have enjoyed increased popularity among the scientific community as they possess even better qualities as cavitands as compared to the former molecules. Moreover, their complexation abilities could further be enhanced with the assistance of metal cations, which can interestingly exert a dual effect on the complexation process: either by competitively binding to the host entity or cooperatively associating with the CB@guest structures. In our previous work, two metal species (Mg2+ and Ga3+) have been found to bind to CB molecules in the strongest fashion upon the formation of host-guest complexes. The current study focuses on their role in the complex formation with three dye molecules: thiazole orange, neutral red, and thioflavin T. Various key factors influencing the process have been recognized, such as pH and the dielectric constant of the medium, the cavity size of the host, Mn+ charge, and the presence/absence of hydration shell around the metal cation. A well-calibrated DFT methodology, solidly based and validated and presented in the literature experimental data, is applied. The obtained results shed new light on several aspects of the cucurbituril complexation chemistry.

Host-guest complexation of cucurbit[7]uril and cucurbit[8]uril with the antimuscarinic drugs tropicamide and atropine.

Cucurbiturils are useful excipients in eye drop formulations: they can increase the water solubility of the drug, enhance drug absorption into the eye, improve aqueous stability and reduce local irritation. Effective and safe drug delivery is, however, a challenge and the information on the host (CBs)/guest (tropicamide and atropine) interactions can help improving the existing treatments and develop novel therapies not limited only to eye diseases/conditions. Since this carrier system can easily modify the properties of the drug and ensure its delivery at the targeted ocular tissue, further insight into the intimate mechanism of the host-guest recognition is crucial. The present DFT/SMD study focuses on the role of numerous factors governing this process, namely the specific position of the guest molecule in the cavity of the cucurbituril, the ionization form (non/protonated) of the antimuscarinic drug, the dielectric constant of the medium, and the size of the cavitant pore. The obtained results are in line with experimental observations and shed light on the mechanism, at atomic resolution, of recognition between the CBs and the two parasympatholytic drugs.

Competition between Ag+ and Ni2+ in nickel enzymes: Implications for the Ag+ antibacterial activity.

Silver's antimicrobial properties have been known for centuries, but exactly how it kills bacteria is still a mystery. Information on the competition between the native Ni2+ and abiogenic Ag+ cations in bacterial systems is also critically lacking. For example, urease, a famous nickel-containing enzyme that hydrolyzes urea into carbon dioxide and ammonia (a key step in the biogeochemical nitrogen cycle on Earth), is inhibited by Ag+ cations, but the molecular mechanism of silver's action is poorly understood. By employing density functional theory (DFT) calculations combined with the polarizable continuum model (PCM) computations we assess the susceptibility of the mono/binuclear Ni2+ binding sites in the nickel enzymatic centers to Ni2+→Ag+ substitution. The active centers in the mononuclear glyoxalase I and acireductone dioxygenase enzymes appear to be well protected against Ag+ attack and, presumably, stay functional even in its presence. On the other hand, the binuclear nickel binding site in urease appears vulnerable to silver attack - the results obtained are in line with available experimental data and give reason to assume a possible substitution of the essential Ni2+ cation from the urease metal center by Ag+.

Holo-chromodulin: competition between the native Cr3+ and other biogenic cations (Fe3+, Fe2+, Mg2+, and Zn2+) for the binding sites.

Chromodulin is an oligopeptide that has an essential role for the flawless functioning of insulin. Although the precise sequence of the constituent amino acid residues and the 3D structure of the molecule has not yet been deciphered, it is known that chromodulin contains only four amino acids in the ratio of Glu-: Gly: Cys: Asp- = 4: 2: 2: 2. An indispensable part for the integrity of the molecule in its active (holo-) form are four chromium cations (hence the name) in the oxidation state of 3+, positioned in two metal binding sites containing one and three Cr3+ ions. Experimental works provide some hints/clues concerning the structure of the metal centers, although their exact composition, type, and arrangement of metal ligating entities remain enigmatic. In the current study, we endeavor to unveil possible structure(s) of the Cr3+ loaded binding sites by strictly following the evidence provided by the experimental data. Well-calibrated in silico methodology for optimization and evaluation of Gibbs free energies is applied and gives strong premises for reliably deciphering the composition/structure of chromodulin metal binding sites. Additional computations reveal the advantage of choosing Cr3+ over other tri- (Fe3+) and divalent (Fe2+, Mg2+, and Zn2+) biogenic ions for securing maximum stability of the metal-occupied binding sites.

Theoretical Insight into the Phosphate-Targeted Silver's Antibacterial Action: Differentiation between Gram (+) and Gram (-) Bacteria.

Although silver is one of the first metals finding broad applications in everyday life, specific key points of the intimate mechanism of its bacteriostatic/bactericidal activity lack explanation. It is widely accepted that the antimicrobial potential of the silver cation depends on the composition and thickness of the bacterial external envelope: the outer membrane in Gram-negative bacteria is more prone to Ag+ attack than the cell wall in Gram-positive bacteria. The major cellular components able to interact strongly with Ag+ (teichoic acids, phospholipids, and lipopolysaccharides) contain mono/diesterified phosphate moieties. By applying a reliable DFT/SMD methodology, we modeled the reactions between the aforementioned constituents in typical Gram-positive and Gram-negative bacteria and hydrated Ag+ species, thus disclosing the factors that govern the process of metal-model ligand complexation. The conducted research indicates thermodynamically possible reactions in all cases but still a greater preference of the Ag+ toward the constituents in Gram-negative bacteria in comparison with their counterparts in Gram-positive bacteria. The observed tendencies shed light on the specific interactions of the silver cation with the modeled phosphate-containing units at the atomic level.

Complexation of trivalent metal cations (Al3+, Ga3+, In3+, La3+, Lu3+) to cucurbiturils: a DFT/SMD evaluation of the key factors governing the host-guest recognition.

Cucurbiturils (CBs), the pumpkin-shaped macrocycles, are suitable hosts for an array of neutral and cationic species. A plethora of host-guest complexes between CBs and a variety of guest molecules has been studied. However, much remains unknown, even in the complexation of very simple guests such as metal cations. In the computational study herein, DFT molecular modeling has been employed to investigate the interactions of a series of trivalent metal cations (Al3+, Ga3+, In3+, La3+, Lu3+) to cucurbit[n]urils and to evaluate the main factors controlling the host-guest complexation. The thermodynamic descriptors (Gibbs energies in the gas phase and in a water environment) of the corresponding complexation reactions have been estimated. This research is a logical continuation of an earlier study on the interaction between CB[n]s and a series of biologically essential mono- and divalent metal cations (Na+/K+ and Mg2+/Ca2+, respectively).

Strontium Binding to α-Parvalbumin, a Canonical Calcium-Binding Protein of the "EF-Hand" Family.

Strontium salts are used for treatment of osteoporosis and bone cancer, but their impact on calcium-mediated physiological processes remains obscure. To explore Sr2+ interference with Ca2+ binding to proteins of the EF-hand family, we studied Sr2+/Ca2+ interaction with a canonical EF-hand protein, α-parvalbumin (α-PA). Evaluation of the equilibrium metal association constants for the active Ca2+ binding sites of recombinant human α-PA ('CD' and 'EF' sites) from fluorimetric titration experiments and isothermal titration calorimetry data gave 4 × 109 M-1 and 4 × 109 M-1 for Ca2+, and 2 × 107 M-1 and 2 × 106 M-1 for Sr2+. Inactivation of the EF site by homologous substitution of the Ca2+-coordinating Glu in position 12 of the EF-loop by Gln decreased Ca2+/Sr2+ affinity of the protein by an order of magnitude, whereas the analogous inactivation of the CD site induced much deeper suppression of the Ca2+/Sr2+ affinity. These results suggest that Sr2+ and Ca2+ bind to CD/EF sites of α-PA and the Ca2+/Sr2+ binding are sequential processes with the CD site being occupied first. Spectrofluorimetric Sr2+ titration of the Ca2+-loaded α-PA revealed presence of secondary Sr2+ binding site(s) with an apparent equilibrium association constant of 4 × 105 M-1. Fourier-transform infrared spectroscopy data evidence that Ca2+/Sr2+-loaded forms of α-PA exhibit similar states of their COO- groups. Near-UV circular dichroism (CD) data show that Ca2+/Sr2+ binding to α-PA induce similar changes in symmetry of microenvironment of its Phe residues. Far-UV CD experiments reveal that Ca2+/Sr2+ binding are accompanied by nearly identical changes in secondary structure of α-PA. Meanwhile, scanning calorimetry measurements show markedly lower Sr2+-induced increase in stability of tertiary structure of α-PA, compared to the Ca2+-induced effect. Theoretical modeling using Density Functional Theory computations with Polarizable Continuum Model calculations confirms that Ca2+-binding sites of α-PA are well protected against exchange of Ca2+ for Sr2+ regardless of coordination number of Sr2+, solvent exposure or rigidity of sites. The latter appears to be a key determinant of the Ca2+/Sr2+ selectivity. Overall, despite lowered affinity of α-PA to Sr2+, the latter competes with Ca2+ for the same EF-hands and induces similar structural rearrangements. The presence of a secondary Sr2+ binding site(s) could be a factor contributing to Sr2+ impact on the functional activity of proteins.

Molecular insights into the interaction of angiotensin I-converting enzyme (ACE) inhibitors and HEXXH motif.

Nutraceuticals and functional foods garner a lot of attention as potential alternative therapies for treatment of (pre)hypertension. Food-derived proteins release large variety of bioactive peptides which are similar in structure to peptide sequences acting in the organism and therefore can modulate their physiological functions. Val-Pro-Pro (VPP) is a milk-derived tripeptide with assumed mild inhibitory activity against angiotensin-converting enzyme (ACE). Computational (DFT) methods are applied on simplified models of Zn2+-HEXXH binding motif without/with bound inhibitors in order to assess the ability of two pharmaceutical drugs (Captopril and Lisinopril) and Val-Pro-Pro to coordinate with Zn2+-HEXXH binding motif of ACE. Both drugs have significant affinity towards the active site, while the Val-Pro-Pro tripeptide has weaker affinity. The obtained results shed light on the thermodynamic aspects of the inhibitors coordination to the Zn2+-HEXXH binding motif of ACE.

Host-Guest Complexation of Cucurbit[7]Uril and Cucurbit[8]Uril with the Antineoplastic and Multiple Sclerosis Agent Mitoxantrone (Novantrone).

The nature of interactions between the neutral/protonated mitoxantrone and the cucurbit[n]uril (n = 7, 8) host system was analyzed by employing density functional theory calculations. A comparison between the inclusion complexes of CB[7] and CB[8] shows various subtle differences in the complexation thermodynamics, given as changes in the Gibbs energy. Doubly and quadruply charged mitoxantrone (MX) molecules spontaneously form complexes in a water solvent, which are modeled using the polarizable continuum model approach. Both CB[7] and CB[8] complexes are stable as the geometry of the cavity allows for electrostatic interactions between the charged MX arms and the rim of the CB cavity. CB[8] also forms a stable complex with two mitoxantrone molecules with their aromatic rings stacked inside the cavity. Both CB[7] and CB[8] show properties that can be utilized in drug delivery.

Competition between abiogenic and biogenic metal cations in biological systems: Mechanisms of gallium's anticancer and antibacterial effect.

Metal cations are key players in a plethora of essential biological processes. Over the course of evolution specific biological functions have been bestowed upon two dozen of (biogenic) metal species, some of the most frequently found being sodium, potassium, magnesium, calcium, zinc, manganese, iron, and copper. On the other hand, there is a group of less studied abiogenic metals like lithium, strontium and gallium that possess not known functions in living organisms, but, by mimicking the native ions and/or competing with them for binding to key metalloenzymes, may exert beneficial effect on humans in particular medical conditions. This review summarizes and critically examines the mechanisms of gallium's therapeutic action in anticancer and antibacterial therapies by exploiting the tools of molecular modeling and experimental biochemistry. These approaches allow for identifying key factors for Ga3+ beneficial effect such as the electrostatic interactions with the protein ligands, substrates or bacterial siderophores, intramolecular hydrogen bond formation, and pH and dielectric properties of the medium. Several intriguing questions concerning the gallium competition with the native ferric ion have found their answers.

Zinc and Its Critical Role in Retinitis pigmentosa: Insights from DFT/SMD Calculations.

Metal cations are required for the proper function of a great amount of biological processes, as they are indispensable cofactors participating in up to 40% of the active sites of the proteins. In the case of some diseases, however, metal cations could exhibit a dual function. As an example, the role of the zinc cation in the development of Retinitis pigmentosa could be given. Experimental works indicate the loss of thermostability of the rhodopsin protein, subjected to the combination of-typical for the disease-mutations and increased quantity of Zn2+. Two structural networks in the intradiscal domain surrounding His100 and His195 are supposed to be susceptible to pathophysiological changes in trace metal concentrations. From a thermodynamic point of view, it is of particular interest to decipher the foundations of the observed outcome, as well as to closely characterize the intimate interactions between the "native" cation and the building amino acid residues of the studied centers. Therefore, the powerful, but fundamentally limited, tools of computational chemistry were applied on simplified models of rhodopsin metal centers in order to shed light on the following aspects: (1) what is the preferred geometry of the Zn2+-containing complexes with the amino acid ligands from the binding pockets; (2) what is the role of the mutations for the interactions between Zn2+ and the examined centers; (3) could other divalent cations such as Ca2+ and Cu2+ substitute for the native zinc; (4) how does the dielectric constant of the environment affect the processes? The obtained results illuminate some aspects of the zinc coordination to amino acid residues and zinc biochemistry related to the presumed pathogenesis of Retinitis pigmentosa.

Gallium as an Antibacterial Agent: A DFT/SMD Study of the Ga3+/Fe3+ Competition for Binding Bacterial Siderophores.

The urgency of finding novel antibacterial drugs (not only antibiotics), exhibiting different mechanisms of therapeutic action, is significant and has served as a premise for recognizing bacteria's siderophores as a plausible drug target. Bacteria secrete siderophores in order to sequester iron(III) from the surrounding medium by binding the essential metal with high affinity. Gallium, on the other hand, is an "abiogenic" ion, known for its anticancer, antibacterial, and anti-inflammatory action. The rationale behind its therapeutic effect lies in its close mimicry of the ferric ion. Since both ions share various physicochemical characteristics, it is of particular interest to understand if gallium could compete with the native ferric ion for binding siderophores and to decipher which molecular characteristics favor Ga3+ binding over Fe3+ binding. It is also well-known that some bacteria are susceptible to gallium-based therapy, while others are not. Therefore, many questions arise such as the following: (1) Which main group/groups building the siderophores promote gallium's attack? (2) Does the combination of the building blocks affect the preference toward a metal? (3) Does the environment play a crucial role? (4) Could the pH of the medium influence the balance between the ions? We try to address these questions by evaluating the free energy of the competition between Ga3+ and Fe3+ ions for siderophore ligands of various structures, denticities, and charge states by employing the tools of the computational chemistry at the DFT/SMD level. Our results not only fall in line with recent experimental data but also complement our knowledge about "Trojan horse" gallium-based therapy.

Complexation of biologically essential (mono- and divalent) metal cations to cucurbiturils: a DFT/SMD evaluation of the key factors governing the host-guest recognition.

Supramolecular complexes based on classical synthetic macrocyclic host molecules such as cyclodextrins and calixarenes have received much attention recently due to their broad applications as biological and chemical sensors, bioimaging agents, drug delivery carriers, light-emitting materials, etc. Cucurbit[n]urils comprise another group of cavitands known for their high affinity for various guest molecules. Nonetheless, some aspects of their coordination chemistry remain enigmatic. Although they are recognized as potential biomimetic scaffolds, they are still not tested as metalloenzyme models and not much is known about their metal-binding properties. Furthermore, there is no systematic study on the key factors controlling the processes of metal coordination to these systems. In the computational study herein, DFT molecular modeling has been employed in order to investigate the interactions of biologically essential (mono- and divalent) metal cations to cucurbit[n]urils and evaluate the major determinants shaping the process. The thermodynamic descriptors (Gibbs energies in the gas phase and in a water medium) of the corresponding complexation reactions have been estimated. The results obtained shed light on the mechanism of host-guest recognition and disclose which factors more specifically affect the metal binding process.