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

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.

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.

Direct Effects of Toxic Divalent Cations on Contractile Proteins with Implications for the Heart: Unraveling Mechanisms of Dysfunction

The binding of calcium and magnesium ions to proteins is crucial for regulating heart contraction. However, other divalent cations, including xenobiotics, can accumulate in the myocardium and enter cardiomyocytes, where they can bind to proteins. In this article, we summarized the impact of these cations on myosin ATPase activity and EF-hand proteins, with special attention given to toxic cations. Optimal binding to EF-hand proteins occurs at an ionic radius close to that of Mg²⁺ and Ca²⁺. In skeletal Troponin C, Cd²⁺, Sr²⁺, Pb²⁺, Mn²⁺, Co²⁺, Ni²⁺, Ba²⁺, Mg²⁺, Zn²⁺, and trivalent lanthanides can substitute for Ca²⁺. As myosin ATPase is not a specific MgATPase, Ca²⁺, Fe²⁺, Mn²⁺, Ni²⁺, and Sr²⁺ could support myosin ATPase activity. On the other hand, Zn²⁺ and Cu² significantly inhibit ATPase activity. The affinity to various divalent cations depends on certain proteins or their isoforms and can alter with amino acid substitution and post-translational modification. Cardiac EF-hand proteins and the myosin ATP-binding pocket are potential molecular targets for toxic cations, which could significantly alter the mechanical characteristics of the heart muscle at the molecular level.

Protein Ca2+-Sites Prone to Sr2+ Substitution: Implications for Strontium Therapy

Strontium (Sr), an alkali metal with properties similar to calcium, in the form of soluble salts is used to treat osteoporosis. Despite the information accumulated on the role of Sr²⁺ as a Ca²⁺ mimetic in biology and medicine, there is no systematic study of how the outcome of the competition between the two dications depends on the physicochemical properties of (i) the metal ions, (ii) the first- and second-shell ligands, and (iii) the protein matrix. Specifically, the key features of a Ca²⁺-binding protein that enable Sr²⁺ to displace Ca²⁺ remain unclear. To address this, we studied the competition between Ca²⁺ and Sr²⁺ in protein Ca²⁺-binding sites using density functional theory combined with the polarizable continuum model. Our findings indicate that Ca²⁺-sites with multiple strong charge-donating protein ligands, including one or more bidentately bound Asp^-/Glu^- that are relatively buried and rigid are protected against Sr²⁺ attack. On the other hand, Ca²⁺-sites crowded with multiple protein ligands may be prone to Sr²⁺ displacement if they are solvent-exposed and flexible enough so that an extra backbone ligand from the outer shell can bind to Sr²⁺. In addition, solvent-exposed Ca²⁺ sites with only a few weak charge-donating ligands that can rearrange to fit the strontium's coordination requirements are susceptible to Sr²⁺ displacement. We provide the physical basis of these results and discuss potential novel protein targets of therapeutic Sr²⁺.

A Network Pharmacology and Multi-Omics Combination Approach to Reveal the Effect of Strontium on Ca2+ Metabolism in Bovine Rumen Epithelial Cells

Strontium (Sr) belongs to the same group in the periodic table as calcium (Ca). Sr level can serve as an index of rumen Ca absorption capacity; however, the effects of Sr on Ca²⁺ metabolism are unclear. This study aims to investigate the effect of Sr on Ca²⁺ metabolism in bovine rumen epithelial cells. The bovine rumen epithelial cells were isolated from the rumen of newborn Holstein male calves (n = 3, 1 day old, 38.0 ± 2.8 kg, fasting). The half maximal inhibitory concentration (IC50) of Sr-treated bovine rumen epithelial cells and cell cycle were used to establish the Sr treatment model. Transcriptomics, proteomics, and network pharmacology were conducted to investigate the core targets of Sr-mediated regulation of Ca²⁺ metabolism in bovine rumen epithelial cells. The data of transcriptomics and proteomics were analyzed using bioinformatic analysis (Gene Ontology and Kyoto Encyclopedia of genes/protein). Quantitative data were analyzed using one-way ANOVA in GraphPad Prism 8.4.3 and the Shapiro-Wilk test was used for the normality test. Results presented that the IC50 of Sr treatment bovine rumen epithelial cells for 24 h was 43.21 mmol/L, and Sr increased intracellular Ca²⁺ levels. Multi-omics results demonstrated the differential expression of 770 mRNAs and 2436 proteins after Sr treatment; network pharmacology and reverse transcriptase polymerase chain reaction (RT-PCR) revealed Adenosylhomocysteine hydrolase-like protein 2 (AHCYL2), Semaphoring 3A (SEMA3A), Parathyroid hormone-related protein (PTHLH), Transforming growth factor β2 (TGF-β2), and Cholesterol side-chain cleavage enzyme (CYP11A1) as potential targets for Sr-mediated Ca²⁺ metabolism regulation. Together these results will improve the current comprehension of the regulatory effect of Sr on Ca²⁺ metabolism and pave a theoretical basis for Sr application in bovine hypocalcemia.

How Theoretical Evaluations Can Generate Guidelines for Designing/Engineering Metalloproteins with Desired Metal Affinity and Selectivity

Almost half of all known proteins contain metal co-factors. Crucial for the flawless performance of a metalloprotein is the selection with high fidelity of the cognate metal cation from the surrounding biological fluids. Therefore, elucidating the factors controlling the metal binding and selectivity in metalloproteins is of particular significance. The knowledge thus acquired not only contributes to better understanding of the intimate mechanism of these events but, also, significantly enriches the researcher's toolbox that could be used in designing/engineering novel metalloprotein structures with pre-programmed properties. A powerful tool in aid of deciphering the physical principles behind the processes of metal recognition and selectivity is theoretical modeling of metal-containing biological structures. This review summarizes recent findings in the field with an emphasis on elucidating the major factors governing these processes. The results from theoretical evaluations are discussed. It is the hope that the physical principles evaluated can serve as guidelines in designing/engineering of novel metalloproteins of interest to both science and industry.

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.

What Is Parvalbumin for?

Parvalbumin (PA) is a small, acidic, mostly cytosolic Ca2+-binding protein of the EF-hand superfamily. Structural and physical properties of PA are well studied but recently two highly conserved structural motifs consisting of three amino acids each (clusters I and II), which contribute to the hydrophobic core of the EF-hand domains, have been revealed. Despite several decades of studies, physiological functions of PA are still poorly known. Since no target proteins have been revealed for PA so far, it is believed that PA acts as a slow calcium buffer. Numerous experiments on various muscle systems have shown that PA accelerates the relaxation of fast skeletal muscles. It has been found that oxidation of PA by reactive oxygen species (ROS) is conformation-dependent and one more physiological function of PA in fast muscles could be a protection of these cells from ROS. PA is thought to regulate calcium-dependent metabolic and electric processes within the population of gamma-aminobutyric acid (GABA) neurons. Genetic elimination of PA results in changes in GABAergic synaptic transmission. Mammalian oncomodulin (OM), the β isoform of PA, is expressed mostly in cochlear outer hair cells and in vestibular hair cells. OM knockout mice lose their hearing after 3–4 months. It was suggested that, in sensory cells, OM maintains auditory function, most likely affecting outer hair cells’ motility mechanisms.

Mechanism of Zn2+ and Ca2+ Binding to Human S100A1

S100A1 is a member of the S100 family of small ubiquitous Ca²⁺-binding proteins, which participates in the regulation of cell differentiation, motility, and survival. It exists as homo- or heterodimers. S100A1 has also been shown to bind Zn²⁺, but the molecular mechanisms of this binding are not yet known. In this work, using ESI-MS and ITC, we demonstrate that S100A1 can coordinate 4 zinc ions per monomer, with two high affinity (K_D~4 and 770 nm) and two low affinity sites. Using competitive binding experiments between Ca²⁺ and Zn²⁺ and QM/MM molecular modeling we conclude that Zn²⁺ high affinity sites are located in the EF-hand motifs of S100A1. In addition, two lower affinity sites can bind Zn²⁺ even when the EF-hands are saturated by Ca²⁺, resulting in a 2Ca²⁺:S100A1:2Zn²⁺ conformer. Finally, we show that, in contrast to calcium, an excess of Zn²⁺ produces a destabilizing effect on S100A1 structure and leads to its aggregation. We also determined a higher affinity to Ca²⁺ (K_D~0.16 and 24 μm) than was previously reported for S100A1, which would allow this protein to function as a Ca²⁺/Zn²⁺-sensor both inside and outside cells, participating in diverse signaling pathways under normal and pathological conditions.

Ca2+/Sr2+ Selectivity in Calcium-Sensing Receptor (CaSR): Implications for Strontium's Anti-Osteoporosis Effect

The extracellular calcium-sensing receptor (CaSR) controls vital bone cell functions such as cell growth, differentiation and apoptosis. The binding of the native agonist (Ca²⁺) to CaSR activates the receptor, which undergoes structural changes that trigger a cascade of events along the cellular signaling pathways. Strontium (in the form of soluble salts) has been found to also be a CaSR agonist. The activation of the receptor by Sr²⁺ is considered to be the major mechanism through which strontium exerts its anti-osteoporosis effect, mostly in postmenopausal women. Strontium-activated CaSR initiates a series of signal transduction events resulting in both osteoclast apoptosis and osteoblast differentiation, thus strengthening the bone tissue. The intimate mechanism of Sr²⁺ activation of CaSR is still enigmatic. Herewith, by employing a combination of density functional theory (DFT) calculations and polarizable continuum model (PCM) computations, we have found that the Ca²⁺ binding sites 1, 3, and 4 in the activated CaSR, although possessing a different number and type of protein ligands, overall structure and charge state, are all selective for Ca²⁺ over Sr²⁺. The three binding sites, regardless of their structural differences, exhibit almost equal metal selectivity if they are flexible and have no geometrical constraints on the incoming Sr²⁺. In contrast to Ca²⁺ and Sr²⁺, Mg²⁺ constructs, when allowed to fully relax during the optimization process, adopt their stringent six-coordinated octahedral structure at the expense of detaching a one-backbone carbonyl ligand and shifting it to the second coordination layer of the metal. The binding of Mg²⁺ and Sr²⁺ to a rigid/inflexible calcium-designed binding pocket requires an additional energy penalty for the binding ion; however, the price for doing so (to be paid by Sr²⁺) is much less than that of Mg²⁺. The results obtained delineate the key factors controlling the competition between metal cations for the receptor and shed light on some aspects of strontium’s therapeutic effects.