Next Article in Journal
In Silico Screening of Natural Compounds for Candidates 5HT6 Receptor Antagonists against Alzheimer’s Disease
Previous Article in Journal
The Influence of High-Intensity Ultrasonication on Properties of Cellulose Produced from the Hop Stems, the Byproduct of the Hop Cones Production
Previous Article in Special Issue
Synthesis of New Thiourea-Metal Complexes with Promising Anticancer Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Coordination Chemistry of Nucleotides and Antivirally Active Acyclic Nucleoside Phosphonates, including Mechanistic Considerations †

1
Department of Chemistry, University of Basel, St. Johannsring 19, CH-4056 Basel, Switzerland
2
Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
*
Authors to whom correspondence should be addressed.
This review is dedicated to Professor Luigi Messory from the University of Florence, Italy, with the very best wishes of the authors for all his future endeavors.
Molecules 2022, 27(9), 2625; https://doi.org/10.3390/molecules27092625
Submission received: 17 March 2022 / Revised: 8 April 2022 / Accepted: 10 April 2022 / Published: 19 April 2022

Abstract

:
Considering that practically all reactions that involve nucleotides also involve metal ions, it is evident that the coordination chemistry of nucleotides and their derivatives is an essential corner stone of biological inorganic chemistry. Nucleotides are either directly or indirectly involved in all processes occurring in Nature. It is therefore no surprise that the constituents of nucleotides have been chemically altered—that is, at the nucleobase residue, the sugar moiety, and also at the phosphate group, often with the aim of discovering medically useful compounds. Among such derivatives are acyclic nucleoside phosphonates (ANPs), where the sugar moiety has been replaced by an aliphatic chain (often also containing an ether oxygen atom) and the phosphate group has been replaced by a phosphonate carrying a carbon–phosphorus bond to make the compounds less hydrolysis-sensitive. Several of these ANPs show antiviral activity, and some of them are nowadays used as drugs. The antiviral activity results from the incorporation of the ANPs into the growing nucleic acid chain—i.e., polymerases accept the ANPs as substrates, leading to chain termination because of the missing 3′-hydroxyl group. We have tried in this review to describe the coordination chemistry (mainly) of the adenine nucleotides AMP and ATP and whenever possible to compare it with that of the dianion of 9-[2-(phosphonomethoxy)ethyl]adenine (PMEA2− = adenine(N9)-CH2-CH2-O-CH2- PO 3 2 ) [or its diphosphate (PMEApp4−)] as a representative of the ANPs. Why is PMEApp4− a better substrate for polymerases than ATP4−? There are three reasons: (i) PMEA2− with its anti-like conformation (like AMP2−) fits well into the active site of the enzyme. (ii) The phosphonate group has an enhanced metal ion affinity because of its increased basicity. (iii) The ether oxygen forms a 5-membered chelate with the neighboring phosphonate and favors thus coordination at the Pα group. Research on ANPs containing a purine residue revealed that the kind and position of the substituent at C2 or C6 has a significant influence on the biological activity. For example, the shift of the (C6)NH2 group in PMEA to the C2 position leads to 9-[2-(phosphonomethoxy)ethyl]-2-aminopurine (PME2AP), an isomer with only a moderate antiviral activity. Removal of (C6)NH2 favors N7 coordination, e.g., of Cu2+, whereas the ether O atom binding of Cu2+ in PMEA facilitates N3 coordination via adjacent 5- and 7-membered chelates, giving rise to a Cu(PMEA)cl/O/N3 isomer. If the metal ions (M2+) are M(α,β)-M(γ)-coordinated at a triphosphate chain, transphosphorylation occurs (kinases, etc.), whereas metal ion binding in a M(α)-M(β,γ)-type fashion is relevant for polymerases. It may be noted that with diphosphorylated PMEA, (PMEApp4−), the M(α)-M(β,γ) binding is favored because of the formation of the 5-membered chelate involving the ether O atom (see above). The self-association tendency of purines leads to the formation of dimeric [M2(ATP)]2(OH) stacks, which occur in low concentration and where one half of the molecule undergoes the dephosphorylation reaction and the other half stabilizes the structure—i.e., acts as the “enzyme” by bridging the two ATPs. In accord herewith, one may enhance the reaction rate by adding AMP2− to the [Cu2(ATP)]2(OH) solution, as this leads to the formation of mixed stacked Cu3(ATP)(AMP)(OH) species, in which AMP2− takes over the structuring role, while the other “half” of the molecule undergoes dephosphorylation. It may be added that Cu3(ATP)(PMEA) or better Cu3(ATP)(PMEA)(OH) is even a more reactive species than Cu3(ATP)(AMP)(OH). – The matrix-assisted self-association and its significance for cell organelles with high ATP concentrations is summarized and discussed, as is, e.g., the effect of tryptophanate (Trp), which leads to the formation of intramolecular stacks in M(ATP)(Trp)3− complexes (formation degree about 75%). Furthermore, it is well-known that in the active-site cavities of enzymes the dielectric constant, compared with bulk water, is reduced; therefore, we have summarized and discussed the effect of a change in solvent polarity on the stability and structure of binary and ternary complexes: Opposite effects on charged O sites and neutral N sites are observed, and this leads to interesting insights.

Graphical Abstract

Table of Contents
1. Some General Considerations
2. Why Is the Antiviral PMEApp4− a Better Substrate for Nucleic Acid Polymerases than (2’-Deoxy)Adenosine 5′-Triphosphate (dATP4−/ATP4−)?
3. Modelling the Interactions of Metal Ions with the Ether Oxygen of PMEA2−
4. The N3 versus N7 Metal Ion-Binding Mode in Acyclic Nucleoside Phosphonates (ANPs) Containing a Purine Moiety
    4.1. The Metal Ion-Coordinating Properties of 9-[2-(Phosphonomethoxy)ethyl]-2-aminopurine (PME2AP)
    4.2. The Metal Ion-Coordinating Properties of 9-[2-(Phosphonomethoxy)ethyl]adenine (PMEA)
    4.3. The Metal Ion-Coordinating Properties of 9-[2-(Phosphonomethoxy)ethyl]-2,6-diaminopurine (PMEDAP)
5. The Various Inhibiting Substituents in 9-[2-(Phosphonomethoxy)ethyl]-2-amino-6-dimethylaminopurine (PME2A6DMAP) Lead to PME-Like Metal Ion-Binding Properties
6. Mechanistic Considerations on the Metal Ion-Promoted Dephosphorylation of ATP
    6.1. Effects of Increasing Amounts of Metal Ions on the Dephosphorylation Rate
    6.2. Oligo Formation and the Effect of Self-Association on the Dephosphorylation Rate
    6.3. Promotion of ATP Hydrolysis by AMP and PMEA, and Inhibition of the Reactivity by Adenine-Altered AMP Relatives
    6.4. Matrix-Assisted Self-Association and Its Significance for Cell Organelles with High ATP Concentrations
7. Solution Structures of Mixed-Ligand Complexes Containing ATP4− and Related Ligands
    7.1. Ligands Containing a Ribose Residue with a Phosphate Group
    7.2. Acyclic Nucleoside Phosphonates. Cu(arm)(PMEA) as an Example with Intramolecular Stack Formation
8. The Effect of a Change in Solvent Polarity on the Stability and Structure of Binary and Ternary Complexes
    8.1. The Dianion of (Phosphonomethoxy)ethane (PME2−) as an Example of a Polar O-Ligand. The Properties of Cu(PME) Are Largely as Expected
    8.2. Nitrogen Donor Sites Are Especially Sensitive to Polarity Changes of the Solvent
    8.3. Competing Solvent Effects on N- versus O-Sites in the Same Ligand and Its Bearing on Complex Stability
9. Discussion
10. Conclusions and Outlook
Abbreviations and Definitions
References

1. Some General Considerations

Nucleobases (Nb), nucleosides (Ns), and nucleotides, in the form of nucleoside monophosphates (NMP2−), nucleoside diphosphates (NDP3−), and nucleoside triphosphates (NTP4−) (Figure 1; upper part), participate in many biologically relevant reactions, most of them being metal ion-dependent [1,2,3]. Due to the increasing charge and the number of available phosphate groups (denticity), the stability of complexes formed with divalent metal ions increases in the order M(NMP) < M(NDP) < M(NTP)2−. This order is independent of the kind of metal ion considered, though of course, e.g., complexes with Zn2+ are considerably more stable than those with Na+. In addition, and depending on the circumstances, the nucleobase may also leave its mark (e.g., Section 4 and Section 5), especially if the focus is on selectivity (e.g., Section 2 and Section 7).
The structures of the important nucleobases are depictured in the lower part of Figure 1 [4,5,6,7]. Note, the nucleosides present in ribonucleic acid (RNA) and in deoxyribonucleic acid (DNA) differ: The ribose in RNA carries on OH group in position 2′ of the ribosyl moiety, whereas in the case of DNA at this position only a hydrogen occurs (see Figure 1; bottom part). Moreover, the important nucleobases in RNA are adenine, cytosine, guanine, and uracil, whereas in DNA uracil is replaced by thymine. Next to these five nucleobases there is one more (Figure 1; middle part)—namely, hypoxanthine with its nucleoside, inosine; these compounds (including their phosphates) play important roles in the metabolism of purines [8,9]—i.e., in the derivatives of adenine and guanine.
The importance of divalent metal ions for the biological transfer of phosphoryl or nucleotidyl groups has long been recognized [10,11,12,13]. It is evident that there are various ways in which a triphosphate chain can coordinate to a metal ion: All three phosphate groups, depending on the size of the metal ion, could possibly directly bind, but also outersphere interactions are possible, etc. In the upper part of Figure 2 [11,14,15,16,17,18,19,20,21,22,23] two structures are indicated that are of relevance with regard to the reactions that are known to occur. The structure at the top shows two metal ions coordinated in such a manner that the terminal γ-phosphate group is activated for transfer [11], thus giving rise to a kinase-type reaction; indeed, such a transphosphorylation mechanism was confirmed for Escherichia coli phosphoenolpyruvate carboxykinase by an X-ray structure [24]. Similarly, in the M(α)-M(β,γ) coordinated complex, the break between Pα and Pβ is activated leading to a nucleotidyl transfer—that is, a polymerase-type reaction, as was concluded from solution [11,22,25] as well as solid state studies [15,16,17,18,19,26]. Evidently, the means that favor the one or other activated species need to be discussed (see Section 2 and Section 6).
Of special interest is the structure at the bottom of Figure 2 ([11,14,15,16,17,18,19,20,21,22,23]), where two metal ions are coordinated to diphosphorylated PMEA, (9-[2-(phosphonomethoxy)ethyl]adenine [20,21,22,23]), giving rise to an acyclic ATP analogue. A nucleophilic attack at Pα is clearly favored by formation of a 5-membered chelate with the ether oxygen of the aliphatic chain [25]. The effect of a phosphonate Pα compared with a phosphate Pα needs to be considered in detail (see Section 2).
Figure 1. In the upper part are shown the chemical structures of the adenosine and uridine 5′-phosphates—namely, adenosine 5′-monophosphate (AMP2−; n = 1), adenosine 5′-diphosphate (ADP3−; n = 2), and adenosine 5′-triphosphate (ATP4−; n = 3), as well as those of uridine 5′-monophosphate (UMP2−; n = 1), uridine 5′-diphosphate (UDP3−; n = 2), and uridine 5′-triphosphate (UTP4−; n = 3). The phosphate groups are named α, β, γ; the γ group being the terminal one. The nucleotides are depictured in their dominating anti conformation [4,5,6,7], which means that the adenine residue is pointing away from the ribose plane as is the (C2)O group in the uridine 5′-phosphates. It is obvious that the substitution with other nucleobases will also lead to the anti conformation. In the lower part of the figure, the structures are shown of the nucleobases (Nb), of the nucleosides (Ns), and of the 2′-deoxynucleosides (dNs). The Abbreviations employed (following the order given in the figure) are as follows: Ade = adenine, Ado = adenosine, and dAdo = 2′-deoxyadenosine; Hyp = hypoxanthine, Ino = inosine, and dIno = 2′-deoxyinosine; Gua = guanine, Guo = guanosine, and dGuo = 2′-deoxyguanosine; Cyt = cytosine, Cyd = cytidine, and dCyd = 2′-deoxycytidine; Ura = uracil, Urd = uridine, and dUrd = 2′-deoxyuridine; Thy = thymine; Thd = thymidine, and dThd = 2′-deoxythymidine = 2′-deoxy-5-methyluridine.
Figure 1. In the upper part are shown the chemical structures of the adenosine and uridine 5′-phosphates—namely, adenosine 5′-monophosphate (AMP2−; n = 1), adenosine 5′-diphosphate (ADP3−; n = 2), and adenosine 5′-triphosphate (ATP4−; n = 3), as well as those of uridine 5′-monophosphate (UMP2−; n = 1), uridine 5′-diphosphate (UDP3−; n = 2), and uridine 5′-triphosphate (UTP4−; n = 3). The phosphate groups are named α, β, γ; the γ group being the terminal one. The nucleotides are depictured in their dominating anti conformation [4,5,6,7], which means that the adenine residue is pointing away from the ribose plane as is the (C2)O group in the uridine 5′-phosphates. It is obvious that the substitution with other nucleobases will also lead to the anti conformation. In the lower part of the figure, the structures are shown of the nucleobases (Nb), of the nucleosides (Ns), and of the 2′-deoxynucleosides (dNs). The Abbreviations employed (following the order given in the figure) are as follows: Ade = adenine, Ado = adenosine, and dAdo = 2′-deoxyadenosine; Hyp = hypoxanthine, Ino = inosine, and dIno = 2′-deoxyinosine; Gua = guanine, Guo = guanosine, and dGuo = 2′-deoxyguanosine; Cyt = cytosine, Cyd = cytidine, and dCyd = 2′-deoxycytidine; Ura = uracil, Urd = uridine, and dUrd = 2′-deoxyuridine; Thy = thymine; Thd = thymidine, and dThd = 2′-deoxythymidine = 2′-deoxy-5-methyluridine.
Molecules 27 02625 g001
Figure 2. The upper and middle parts of the figure show structures of two M2(NTP) complexes, where NTP4− = nucleoside 5′-triphosphate: In one case, a M(α,β)-M(γ) coordination is depictured (upper part) indicating the structure relevant for transphosphorylations (kinase, etc.). In the other case, the metal ions are bound in a M(α)-M(β,γ)-type fashion (middle part) that is relevant for nucleic acid polymerases, which catalyze the transfer of a nucleotidyl unit. For the latter binding mode, the structure needs to be enforced by the enzyme; this means, the two metal ions need to be anchored [11,12,13,14] to amino acid side chains, often carboxylate groups of aspartate or glutamate residues of the enzyme [15,16,17,18,19]. The (lower part) shows a complex formed between two metal ions and diphosphorylated PMEA (9-[2-(phosphonomethoxy)ethyl]adenine)—that is, M2(PMEApp). Note, the M(α)-M(β,γ) binding mode, crucial for the polymerase reaction, is favored by the formation of a 5-membered chelate with the ether oxygen of the aliphatic chain (see Section 3). Of course, the adenine residue may be replaced by any other nucleobase moiety and the nucleophile (N) may in addition interact with the M2+ at the α-phosphonate group. Altered versions of similar situations are depictured in Refs [20,21,22,23].
Figure 2. The upper and middle parts of the figure show structures of two M2(NTP) complexes, where NTP4− = nucleoside 5′-triphosphate: In one case, a M(α,β)-M(γ) coordination is depictured (upper part) indicating the structure relevant for transphosphorylations (kinase, etc.). In the other case, the metal ions are bound in a M(α)-M(β,γ)-type fashion (middle part) that is relevant for nucleic acid polymerases, which catalyze the transfer of a nucleotidyl unit. For the latter binding mode, the structure needs to be enforced by the enzyme; this means, the two metal ions need to be anchored [11,12,13,14] to amino acid side chains, often carboxylate groups of aspartate or glutamate residues of the enzyme [15,16,17,18,19]. The (lower part) shows a complex formed between two metal ions and diphosphorylated PMEA (9-[2-(phosphonomethoxy)ethyl]adenine)—that is, M2(PMEApp). Note, the M(α)-M(β,γ) binding mode, crucial for the polymerase reaction, is favored by the formation of a 5-membered chelate with the ether oxygen of the aliphatic chain (see Section 3). Of course, the adenine residue may be replaced by any other nucleobase moiety and the nucleophile (N) may in addition interact with the M2+ at the α-phosphonate group. Altered versions of similar situations are depictured in Refs [20,21,22,23].
Molecules 27 02625 g002

2. Why Is the Antiviral PMEApp4− a Better Substrate for Nucleic Acid Polymerases than (2′-Deoxy)Adenosine 5′-Triphosphate (dATP4−/ATP4−)?

The advent of the human immunodeficiency viruses (HIV-1 and HIV-2) [27,28,29] spurred research on antivirals. As the virus uses the cellular machinery of the host for its own reproduction (for details see [22,30]), the transcription machinery including nucleotide derivatives of such compounds were also developed and studied in detail. Nucleotides were systematically altered at the nucleobase, at the sugar moiety, and at the phosphate group (e.g., thiophosphate [31,32,33]) with the aim of discovering compounds that inhibit the synthesis of viral DNA or RNA. Among early successful compounds were acyclic nucleoside phosphonates [30], which have the advantage that the phosphorus–carbon bond is not split by dephosphorylation enzymes and that they have no OH group available that allows the continuation of a growing nucleic acid chain—instead, the incorporation of such a derivative into the growing chain leads to its termination [34,35].
9-[2-(Phosphonomethoxy)ethyl]adenine (PMEA) is one of these acyclic nucleoside phosphonates [34], and we will concentrate here on this example. As pointed out already in the preceding Section, this adenine derivative needs to be diphosphorylated to PMEApp4−, which is an ATP/dATP analogue, as is indicated in the lower part of Figure 2. In fact, initially PMEApp4− is a better substrate for polymerases than (d)ATP [34] for three main reasons:
-
PMEA has an anti-like conformation just as AMP.
-
The phosphonate group owns an increased basicity and therefore also an enhanced metal ion affinity, and
-
formation of a 5-membered chelate involving the ether oxygen favors coordination at the Pα group.
  • These will be discussed below.
In Figure 3, the anti-like conformation of PMEA2− [36,37] is shown, and it is expected that this also holds for PMEApp4− ([4,5,6,36,37,38]). Hence, PMEApp4− is further expected to fit into the active site of a polymerase just as well as ATP4− or dATP4−.
An oxygen atom is more electronegative than a carbon atom. Therefore, a P–O bond is more polar than a P–C bond, and a phosphonyl group is more basic than a phosphoryl one. This conclusion is confirmed by the data of Table 1 [39,40,41], where the proton affinities (and thus the basicities) of methyl phosphate and methylphosphonate (=P2−) are compared [39,40,41,42,43] based on the following two deprotonation equilibria:
H 2 ( P )   H ( P ) + H + K H 2 ( P ) H = [ H ( P ) ] [ H + ] / [ H 2 ( P ) ]
H ( P )   P 2 + H + K H ( P ) H = [ H ( P ) ] [ H + ] / [ H 2 ( P ) ]
The average difference Δ pKa = p K H ( P ) H p K H 2 ( P ) H = 5.3 ± 0.2 encompasses the two individual values for CH3OPO(OH)2 (5.26 ± 0.2) and CH3PO(OH)2 (5.41 ± 0.03). Indeed, these values are in the expected order [41,42,43], and the differences for phosphates and phosphonates are alike, even though the actual acidity constants differ (Table 1).
An enhanced basicity should give rise to an enhanced metal ion affinity. This postulate can be probed by constructing plots of log K M ( R - PO 3 ) M (y-axis) versus p K H ( R - PO 3 ) H (x-axis), where R - PO 3 2 is a phosphate or phosphonate ligand. Indeed, if the residue R does not interfere with metal ion binding (neither in a positive nor negative way), all corresponding data points fit on straight lines, as has been shown for more than 10 metal ions [44,45,46,47]. These lines can be quantified with the straight-line Equation (3)
y = m · x + b
where x represents the pKa value of the phosphate monoester or phosphonate ligand and y the log stability constant of the complex. The slope m and the intercept b with the y-axis are characteristic for a given metal ion.
We consider now as an example the situation with zinc(II), an important bio-metal ion (Equation (4)):
log   K Zn ( R - PO 3 ) Zn = ( 0.345   ±   0.026 )   p K H ( R - PO 3 ) H ( 0.017   ±   0.171 )
This straight-line equation holds in the pKa range of about 5 to 8 with an error limit of ±0.060 log units (3σ) [45,46,47].
Application of p K CH 3 OP ( O ) 2 OH H = 6.36 (Table 1) to Equation (4) gives log K Zn ( CH 3 OPO 3 ) Zn = 2.18 ± 0.06, and application of p K CH 3 P ( O ) 2 OH H = 7.51 (Table 1) leads to log K Zn ( CH 3 PO 3 ) Zn = 2.57 ± 0.06; hence, the phosphonate complex Zn(CH3PO3) is by 0.39 ± 0.08 log units more stable than the phosphate complex Zn(CH3OPO3). Or, to express it differently, the Zn2+-phosphonate complex is by a factor of 2.5 more stable than the Zn2+-phosphate one.
In the structure at the bottom of Figure 2, it is indicated that metal ion coordination is at the Pα group—that is, next to the increased basicity of this group—apparently further favored by the formation of a 5-membered chelate involving the ether oxygen. This suspicion will be dealt with in the next Section 3.

3. Modelling the Interactions of Metal Ions with the Ether Oxygen of PMEA2−

To be unequivocal in the evaluations, we need a ligand that offers for metal ion binding only the phosph(on)ate group [48] and the ether oxygen and that should structurally be as close to PMEA2− as possible. These requests are fulfilled by the dianion of ethoxymethanephosphonate, also addressed as (phosphonomethoxy)ethane (=PME2−), which is PMEA2− that has lost its adenine residue. As outlined in the preceding Section 2, plots of log K M ( R - PO 3 ) M (y-axis) versus p K H ( R - PO 3 ) H (x-axis), where R of the ligand R - PO 3 2 is a residue that does in no way interfere with metal ions, give rise to straight-lines (Equation (3)). Figure 4 shows the straight-line examples for Mg2+, Mn2+, and Cu2+ [44,45].
The data points for the corresponding complexes of PMEA2− and PME2− show that at least in the case of Cu2+ the adenine residue participates in metal ion binding because Cu(PMEA) is more stable than Cu(PME) (see Figure 4). More important for the present is the observation that the M(PME) complexes are more stable than is predicted based on the basicity of the phosphonate group; hence, the following intramolecular Equilibrium (5) must be of relevance:
Molecules 27 02625 i001
Its position is determined by the intramolecular, and dimension-less, equilibrium constant KI and the ratio of the “closed” (cl) and “open” isomers (Equation (6) where for the present PME2− = PE2−):
K I = [ M ( PE ) cl ] [ M ( PE ) op ]
Consequently, the overall complex formation is described by Equilibrium (7)
M 2 + + PE 2     M ( PE ) op     M ( PE ) cl
and the experimentally measured stability of the M(PE) complex by Equation (8):
K M ( PE ) M = ( [ M ( PE ) op ] + [ M ( PE ) cl ] ) [ M 2 + ] [ PE 2 ] = [ M ( PE ) op ] [ M 2 + ] [ PE 2 ] + [ M ( PE ) cl ] [ M 2 + ] [ PE 2 ]
From Equations (6) and (8) follow [45,47,49,50] Equations (9) and (10):
K M ( PE ) M   = K M ( PE ) op M + K I   ·   K M ( PE ) op M = K M ( PE ) op M ( 1 + K I )
K I = K M ( PE ) M K M ( PE ) op M 1 = 10 log   Δ 1
The stability constant of the open isomer (Equation (11))
K M ( PE ) op M = [ M ( PE ) op ] / ( [ M 2 + ] [ PE 2 ] )
is not directly accessible by experiments, but it may be calculated by applying the acidity constant K H ( PE ) H (analogue to Equation (2)) to the straight-line equations, as defined by Equation (3). Now with the stability constant of the open isomer available, the stability-constant difference (= stability enhancement; Figure 4) according to Equation (12)
log   Δ = log   Δ M ( PE ) = log   K M ( PE ) M K M ( PE ) op M = log   K M ( PE ) exper M log   K M ( PE ) calc M
can be calculated; this also defines the second term in the above Equation (10). Furthermore, the log Δ values correspond to the vertical dotted lines in Figure 4. From Figure 4, it is evident that the stability enhancements (Equation (12)) for the M(PME) and M(PMEA) complexes can differ. To place this conclusion on solid grounds, the mentioned values are compared according to Equation (13), where for the present PMEA = PE:
Δ log ΔM/PE/PME = log ΔM/PE − log ΔM/PME
These results are listed in Table 2 (Column 4).
Evidently, with values for KI known, the percentages for the closed or chelated species in Equilibrium (5) can be calculated according to Equation (14):
% M(PE)cl = 100 · KI/(1 + KI)
Of course, for any calculation of this type, well-defined error limits are compulsory. The corresponding results are summarized in Table 2 [45].
In Columns 2 and 3 of Table 2, the stability enhancements (Equation (12)) for the M(PMEA) and M(PME) complexes are listed [44,45]. These data indicate for all complexes of the various metal ions that they are more stable than expected based on the basicity of the phosphonate group (see also Figure 4). As these values have the same dimension, they can be directly compared, as expressed by Equation (13); the results of these comparisons are listed in Column 4 of Table 2. Interestingly, for the metal ions Mg2+, Ca2+, Sr2+, Ba2+, Mn2+, Co2+, Zn2+, and Cd2+, the values for Δ log Δ are zero within the error limits, meaning that the M(PMEA) and M(PME) complexes for a given metal ion show the same stability enhancement and thus the same extent of metal ion–ether oxygen interaction in Equilibrium (5). The corresponding formation degrees (Equation (14)) of the M(PME)cl isomers are listed in the terminal Column of Table 2; the same formation degrees hold within the error limits for the M(PMEA)cl isomers.
The Cu(PMEA) complex is the sole significant exception in Column 4 of Table 2 that shows a further stability enhancement that is beyond that of Cu(PME)—namely, by about 0.3 log unit (see Table 2, Column 4). A corresponding, though much smaller, effect is observed for Ni(PMEA) with Δ log Δ = 0.11 ± 0.09, which is just at the edge of significance. However, one has to conclude for both instances that the adenine residue must be responsible for this effect [44,45,51]. A M2+–N1 interaction cannot be the reason because a phosphonate-coordinated metal ion cannot reach N1 [44,45]. However, interactions appear possible with the other two N atoms; indeed, macrochelate formation with N7 is well-known for many related systems [11,50,51,52,53,54,55,56]. The situation with N3 [44] is more complicated because N3 can only be reached if the metal ion coordinates simultaneously with the ether O atom (see the PMEA2− structure in the upper part of Figure 3), the result being that a 5-membered and a 7-membered chelate are linked to each other [44,45,51]. It appears that the latter isomer involving N3 dominates (for further details see Section 4). Of course, the total amount of chelated species for these two instances is rather large: % Cu(PMEA)cl/tot = 83 ± 3 and % Ni(PMEA)cl/tot = 50 ± 8 [25,44,45].

4. The N3 versus N7 Metal Ion-Binding Mode in Acyclic Nucleoside Phosphonates (ANPs) Containing a Purine Moiety

4.1. The Metal Ion-Coordinating Properties of 9-[2-(Phosphonomethoxy)ethyl]-2-aminopurine (PME2AP)

Research on acyclic nucleoside phosphonates (ANPs) containing a purine residue revealed that the kind and the position of the substituent (at C2 or C6) has a significant influence on the biological activity. For example, the shift of the (C6)NH2 group in PMEA to the C2 position (Figure 5) leads to an ANP isomer, i.e., PME2AP, which shows only a moderate antiviral activity compared with that of PMEA [57].
Of course, the removal of the (C6)NH2 substituent removes the steric effect of this group on N7, and consequently, macrochelate formation involving N7 is strongly enhanced in M(PME2AP) species [43,58], whereas no metal ion interaction with N3 could be discovered. This is because the NH2 substituent at C2 inhibits metal ion binding to N3. This observation supports the hypothesis that the C6 substituent plays an important role for the biological activity.
For the M(PME2AP) complexes, the stability increases between 0.13 and 1.06 log units (=log ΔM/PME2AP) [43]. For the complexes of Mg2+, Ca2+, (Sr2+, Ba2+ were not studied), and Mn2+, only Equilibrium (5) is responsible for the increased complex stabilities. However, for the complexes of Co2+, Ni2+, Cu2+, Zn2+, and Cd2+, macrochelate formation is of relevance as well [43], which means that the Equilibrium Scheme (15) in Figure 6 describes the situation.
The formation degrees of the two relevant and chelated species M(PME2AP)cl/O and M(PME2AP)cl/N7, plus the “open” species M(PME2AP)op are summarized in Table 3. The larger formation degree of Ni(PME2AP)cl/N7 compared with Cu(PME2AP)cl/N7 corresponds to the situation with the M(AMP) complexes [50] and can be explained by statistical considerations (see Section 9) based on the different geometries of the coordination spheres of Ni2+ and Cu2+ [59].

4.2. The Metal Ion-Coordinating Properties of 9-[2-(Phosphonomethoxy)ethyl]adenine (PMEA)

The best known ANP of the compounds shown in Figure 5 is PMEA, which exhibits various biological activities [59,60] and which in the form of its bis(pivaloyloxymethyl)ester [22] is used as an oral prodrug (adefovir dipivoxil) already for many years [59]. This neutral prodrug passes membranes more easily than PMEA2−, which is released inside the cell and then diphosphorylated giving PMEApp4− as the biologically active dATP4−/ATP4− analogue (Section 2).
For most of the studied M(PMEA) complexes, it is Equilibrium (5) that explains all observed properties (Table 2; Section 3). The significant exception is Cu(PMEA), which is about 0.3 log unit more stable beyond the ether O–metal ion interaction (Section 3). This is also seen in Figure 7 [44,61,62,63] if the stabilities of the Cu(PMEA) and Cu(PME-R) are compared; this extra stability must be attributed to a metal ion interaction with the adenine residue.
Based on steric and chemical considerations [45,64], it was concluded that the Equilibrium Scheme (17) (Figure 6) involving Cu(PMEA)cl/O/N3 is the relevant one and not Equilibrium Scheme (15) with Cu(PMEA)cl/N7. This conclusion was confirmed by 1H NMR line-broadening experiments [65]; hence, the (C6)NH2 substituent inhibits Cu2+ coordination at N7 and the ether O atom-binding of Cu2+ facilitates N3 coordination via adjacent 5- and 7-membered chelate rings.
Later studies with 9-(4-phosphonobutyl)adenine (=3′-deoxa-PMEA2− = dPMEA2−) [66] provided evidence for the formation of in total four isomers—that is, in addition to those appearing in the Equilibrium Scheme (17) (Figure 6), also a macrochelated species involving N7, Cu(PMEA)cl/N7 may form in traces, and the Equilibrium Scheme (16) is valid. The calculated percentages (in parentheses) [66] for the four isomers are Cu(PMEA)op (17 ± 3%), Cu(PMEA)cl/O (34 ± 10%), Cu(PMEA)cl/N7 (7.7 ± 5.3%), and Cu(PMEA)cl/O/N3 (41 ± 12%). These results are based on the assumption that the butyl chain of dPMEA2− and the one of PMEA containing the ether oxygen in their formation of the macrochelate with N7 behave alike; there is, however, a small caveat here: The macrochelate with dPMEA2− could be somewhat stabilized by a hydrophobic interaction between the butyl chain and the imidazole ring of the adenine residue; then the given estimation of 7.7 ± 5.3% for Cu(PMEA)cl/N7 would be somewhat too large.

4.3. The Metal Ion-Coordinating Properties of 9-[2-(Phosphonomethoxy)ethyl]-2,6-diaminopurine (PMEDAP)

We have already noted that PMEA (Figure 5) has antiviral properties and that those of PME2AP are smaller. The related ANP, PMEDAP (Figure 5), is also potentially a useful antiviral compound [34,67,68]. However, PMEDAP has a more pronounced cytostatic effect than PMEA; it affects human leukemia cell lines [69] and inhibits the cellular DNA polymerase δ [70].
PMEDAP is an interesting molecule with an amino group each at C2 and C6. In Figure 7, it is seen that its Cu2+ complex has about the same stability as Cu(PMEA) and that both complexes are more stable than Cu(PME-R) (R = non-interfering residue). This proves that the nucleobase residues must be involved in Cu2+ binding. Considering that the (C2)NH2 group inhibits metal ion binding at N3, this observation seems surprising because N3, like with PMEA (Section 4.2), cannot be the cause. However, due to the further NH2 group in PMEDAP, compared with PMEA, at the adenine residue the basicity of the 2,6-diaminopurine moiety is by Δ pKa of about 0.7 more basic than the adenine residue [62]. This basicity increase enhances the Cu2+ affinity for N7 so much that the steric inhibiting effect of (C6)NH2 on N7 is compensated. To summarize, Cu(PMEA) has an increased complex stability due to an N3 coordination and Cu(PMEDAP) due to an involvement of N7 in macrochelate formation. These different binding sites in the seemingly similar complexes are confirmed by 1H NMR line-broadening experiments [62].
In Table 4 in Columns 2 and 3, the stability enhancements (Equation (12)) for the M(PMEDAP) and M(PME-R) complexes are listed, and in Column 4, these stability enhancements are compared according to Equation (13). This comparison shows that for the M(PMEDAP) complexes of Mg2+, Ca2+, Sr2+, Ba2+, Mn2+, (Zn2+), and Cd2+, only Equilibrium (5) is of relevance; hence, the nucleobase moiety does not participate in these instances in metal ion binding. This is different for the complexes of Co2+, Ni2+, and Cu2+, where the 2,6-diaminopurine moiety is involved in metal ion binding (see Section 4.2). That apparently the Δ log Δ value for Ni(PMEDAP) is a bit larger than that for Cu(PMEDAP) (Table 4, Column 4) corresponds to the observation made (vide supra) with the corresponding M(PME2AP) complexes (Section 4.1), as well as with those for M(AMP), and can be explained, as said, by statistical considerations based on the different geometries of the coordination spheres of Ni2+ and Cu2+ [50] (see Section 9).
From Column 6 in Table 4, it follows that the formation degrees of the 5-membered chelates, only Equilibrium (5) being of relevance, varies for the six metal ions in question between about 15% (Sr2+) and 50% (Cd2+). In the case of Co2+, Ni2+, and Cu2+ complexes, as discussed, the macrochelate formed by the phosphonate-coordinated metal ions with N7 is of relevance. The formation degrees of these macrochelates, M(PMEDAP)cl/N7, are for Co2+, Ni2+, and Cu2+ 26 ± 14%, 54 ± 10%, and 43 ± 11%, respectively [62].
Interestingly, 9-[2-(phosphonomethoxy)ethyl]-2-amino-6-dimethylaminopurine (PME2A6DMAP) shows similar therapeutic effects as the parent PMEDAP (Figure 5) [60,71], although it exerts even a stronger cytostatic activity in mouse leukemia L1210 cells or in human cervix carcinoma [60]. In the indicated context of the biological activity, it is amazing to note that the metal ion-coordinating properties between these two ANPs differ considerably. Those for PMEDAP were discussed in this Section 4.3; in the next Section 5, those for PME2A6DMAP (Figure 5) will follow.

5. The Various Inhibiting Substituents in 9-[2-(Phosphonomethoxy)ethyl]-2-amino-6-dimethylaminopurine (PME2A6DMAP) Lead to PME-Like Metal Ion-Binding Properties

As far as the proton affinity of the N1 site is concerned, we have seen that the introduction of a second amino group into the adenine moiety makes the N1 site by Δ pKa ca. 0.7 more basic (Section 4.3). A similar situation is observed now for PMEA and PME2A6DMAP, with an enhanced proton affinity of N1 by Δ pKa ca. 0.9 ( p K H 2 ( PMEA ) H = 4.16 ± 0.02 and p K H 2 ( PME 2 ADMAP ) H = 5.09 ± 0.02 [61]). In contrast, the basicity of the phosphonate group is very similar for the two ANPs [61] (see also the p K H ( PE ) H values in Figure 7). However, this increase in basicity of the nucleic base moiety has no influence on the stability of the Cu2+ complexes of PME2A6DMAP2−, in contrast with that of PMEDAP2− (Figure 7): The stabilities of Cu(PMEA) and Cu(PMEDAP) are enhanced and very similar, whereas the stability of Cu(PME2A6DMAP) is close to the nucleic base-free Cu(PME-R) complex.
The reason lies in the structure of the M(PME2A6DMAP) complex. The ligand PME2A6DMAP2− offers in principle four metal ion binding sites, N1 not being accessible [44,45] (Section 3), leaving N3, N7, the ether oxygen of the –CH2–CH2–O–CH2 PO 3 2 chain, and the phosphonate group. From Figure 7, it follows that the phosphonate group is the primary metal ion-binding site, and therefore, the data need to be evaluated first for the occurrence of Equilibrium (5). For this reason, we have collected in Column 2 of Table 5 [61] the stability enhancements for log ΔM(PME2A6DMAP) = log ΔM(PE) (Equation (12)). These data are compared according to Equation (13) with the log ΔM(PME-R) values (Table 5, Column 3), which contain solely Equilibrium (5) as a source for an increased complex stability. It is revealing to see that all the values for Δ log Δ (Equation (13)) in Column 4 of Table 5 are zero within the error limits. From this result for Δ log Δ, it follows that the nucleobase of PME2A6DMAP2− does not participate in metal ion binding but that the coordination of the ether O atom is solely responsible for the increased complex stability.
Application of log Δ = log ΔM(PE) = log ΔM(PME2A6DMAP) to Equation (10) provides the intramolecular equilibrium constant KI = KI/PE for Equilibrium (5) (Table 5, Column 5), and from this the formation degree of the M(PME2A6DMAP)cl/O species follows with Equation (14) (Table 5, Column 6). This formation degree varies between about 15% to 70%, depending on the metal ion involved in Equilibrium (5). As one would expect, these formation degrees are within their error limits identical with those listed in the terminating Column of Table 5 for the % M(PME-R)cl/O species. In other words, PME2A6DMAP2− coordinates to metal ions in a PME-like manner, as seen in Equilibrium (5). However, this does not exclude, e.g., stacking interactions [72] between the purine residue of PME2A6DMAP2− in a protein or nucleic acid.

6. Mechanistic Considerations on the Metal Ion-Promoted Dephosphorylation of ATP

6.1. Effects of Increasing Amounts of Metal Ions on the Dephosphorylation Rate

The rate of hydrolysis of ATP [10] and other NTPs [14,73] increases with decreasing pH—i.e., with an increasing proton concentration [14]. This is confirmed in a quantitative manner by the following considerations: The pKa values for the deprotonation of H2(ATP)2− are p K H ( ATP ) H = 6.47 ± 0.01 and p K H 2 ( ATP ) H = 4.00 ± 0.01 (25 °C; aqueous solution; I = 0.1 M, NaNO3 or NaClO4) [74]—that is, for the deprotonation of the terminal γ-phosphate group and the (N1)H+ site, respectively. The initial rate, v0, of the dephosphorylation reaction of ATP in 10−3 M aqueous solution at 50 °C and pH0 ca. 9 amounts to v0 = 1.5 × 10−10 M s−1; at pH0 5.5 one obtains v0 = 10 × 10−10 M s−1. This means, roughly speaking, at 50 °C and I = 0.1 M (Na+) H(ATP)3− hydrolyzes by a factor of about 7 faster than ATP4−. The addition of one equivalent of Cu2+ to the 10−3 M ATP solution at 50 °C and pH0 5.5 has significant consequences because now v0 = 25 × 10−9 M s−1; in other words, the dephosphorylation rate is increased by a factor of about 25.
The mentioned observations stipulated the experiments shown in Figure 8, where the initial rate v0 at pH0 5.5 is plotted for the dephosphorylation reaction in dependence on the concentration of a second metal ion. In the case of Cu2+ as the second metal ion, the effect is most pronounced (Figure 8) because with one more equivalent of Cu2+ the system reaches already its maximum rate. The reason is that the Cu2(ATP) complex is rather stable [10] and that the second Cu2+ is present at pH0 5.5 as Cu(OH)+, and this allows an intramolecular OH attack at the terminal γ-phosphate group. This type of reactivity involving OH applies to some extent also for Zn2+, which is an extremely versatile metal ion [75,76].
However, Mg2+, Ni2+, and Cd2+ reach as the second metal ion and as promoters of the dephosphorylation reaction all a limiting v0 value of about 10·10−8 M s−1—that is, compared with the Cu2+:ATP = 1:1 system (v0 = 2.5 × 10−8 M s−1), the dephosphorylation rate is enhanced by a factor of about four. Interestingly, such a second metal ion is of relevance also for kinases [77,78,79]. Note, the only difference between the mentioned three metal ions is that the maximum rate of 10 × 10−8 M s−1 is reached at different concentrations; in other words, the stabilities of these ATP4− complexes decrease in the order Cd2+ > Ni2+ > Mg2+ [74], and the reaction proceeds with H2O as the nucleophile.
The situation with Cu(dien)2+ as promoter is most interesting because in Cu(dien) only a single equatorial binding site at Cu2+ is left. Despite this handicap, with a 15-fold excess of Cu(dien) over Cu(ATP)2−, the same maximum rate is obtained as with Mg2+, Ni2+, and Cd2+. Hence, generally speaking, the reactive complex has the composition M[Cu(ATP)]. This is in accord with the observations that for a given metal ion the reactive species has the composition M2(ATP), as is confirmed [10,14,80] for several metal ions at various pH values by Job’s plots [81]. In this context, one may mention that polymerases are known (for DNA [25,82,83,84] and RNA [2,25,83,84,85]) that contain an intrinsic divalent metal ion, often Zn2+ [82], but in order to be active with a NTP4− substrate, also an extrinsic second divalent metal ion, often Mg2+ (or Mn2+), is needed [26,84,86,87].

6.2. Oligo Formation and the Effect of Self-Association on the Dephosphorylation Rate

Purine nucleotide residues, as they occur in DNA or RNA, form intramolecular stacks and determine in this way the structure of these macromolecules. Such π stacks can also occur between simple purine units [88] leading to associates [89,90,91,92], and their formation is reflected in the 1H NMR spectrum that changes with increasing nucleotide concentration. Plots of the changes in dependence on concentration can be successfully interpreted with the isodesmic model of non-cooperative self-association [89]; i.e., it is assumed, e.g., for an adenine derivative (A) that for Equilibrium (18) the association constants for the various steps are all equal (Equation (19)):
( A ) n + A     ( A ) n + 1
K = [(A)n+1]/[(A)n][A]
In Table 6, the association constants for adenosine (Ado) and several of its derivatives are listed [89,93]. The comparison of the constants for adenosine and AMP2− shows that the negative charge leads to repulsion and thus inhibits the association; indeed, the constants decrease further for ADP3− and ATP4−. However, for Mg(ATP)2−, the association increases (K = 4.0 ± 0.5 M−1) beyond that of AMP2− (K = 2.1 ± 0.4 M−1), despite the equality of the charges. This indicates that Mg2+ does not only compensate part of the negative charge but that this metal ion links in addition the phosphate residues of different ATPs together and that this stabilizes the oligomers further [89].
Yet, for Zn(ATP)2− and Cd(ATP)2−, the constants are much larger than those for Mg(ATP)2− (Table 6), indicating that interactions with the adenine residue must now also occur. Indeed, the chemical shifts for (C8)H indicate that at low concentrations a downfield shift occurs, whereas all other shifts are upfield [90]; this suggests the preferred formation of Zn2+/N7 bridged [ Zn ( ATP ) ] 2 4 dimers, which then associate further in an unbridged manner [89]. The analogous observation is made with Cd2+/ATP4−. For M2+ = Zn2+ or Cd2+, one may write the following two equilibria:
2   M ( ATP ) 2   [ M ( ATP ) ] 2 4 K D * = [ [ M ( ATP ) ] 2 4 ] / [ M ( ATP ) 2 ] 2
[ M ( ATP ) ] 2 4 + M n ( ATP ) n 2 n   M 2 + n ( ATP ) 2 + n ( 4 + 2 n ) K st = [ M 2 + n ( ATP ) 2 + n ( 4 + 2 n ) ] / [ [ M ( ATP ) ] 2 4 ] [ M n ( ATP ) ] n 2 n ]
In Figure 9, we plotted the formation degrees of monomeric to octameric species in the ATP concentration range 0 to 0.4 M. The preferred formation of dimers in the Zn2+/ATP = 1:1 system is evident. These data allow definition of the concentrations one should use if one wants to describe the properties of monomers. For example, in 1 mM solutions of Mg(ATP)2−, about 99% are present in the monomeric form, whereas under the same conditions for Zn(ATP)2−, only about 96% are monomers [53]. It is thus advisable to work with [ATP] < 1 mM in studies aimed for the properties of monomers.
That for the dephosphorylation reaction of ATP in the presence of the mentioned metal ions, Ni2+, Cd2+, Cu2+, and Zn2+, the formation of dimeric species for reaching the reactive state is important (in the presence of one or two equivalents of M2+), follows from plots of log v0 versus log [ATP]tot; ATP being present with one or two equivalents of metal ions. These plots show slopes of two for Cd2+, e.g., at pH0 7.2, 8.2, and 9.0 [80], for Cu2+ at pH0 5.5 and 6.7, for Ni2+ at pH0 8.0, and for Zn2+ at pH0 7.2 [10], meaning that the reactive M(ATP)2− or M2(ATP) complex is a dimer. The latter species is the more reactive one, as is also confirmed by Job’s plots (Section 6.1); thus [M2(ATP)]2 or [M2(ATP)]2OH, depending on the metal ion involved, is the most reactive complex.
The tentative structure of such a dimer is shown in the upper part of Figure 10 [14,80,94]. It is evident that one of the two ATP molecules (left part) is used to organize with a metal ion (M2+) a bridge to N7 and to stabilize the stack, whereas the other ATP is in the reactive state for OH or H2O attack and thus for the release of the terminal γ-phosphate group. One could conclude that ATP is its own “enzyme” [55]. Indeed, this view is confirmed in the next section.

6.3. Promotion of ATP Hydrolysis by AMP and PMEA, and Inhibition of the Reactivity by Adenine-Altered AMP Relatives

The dimer shown in the upper part of Figure 10 consists of two M2(ATP) units. The left part of this dimer is the structuring part, and “removing” the β,γ-phosphate group leaves an AMP unit with one metal ion. Hence, the question arises whether the remainder—that is, M3(ATP)(AMP)—is a reactive complex. Indeed, with M2+ = Cu2+, this is the case, as is shown in Figure 11.
In a Cu2+/ATP 2:1 system with [Cu]tot = 2 × 10−3 M and [ATP]tot = 10−3 M, the addition of AMP enhances the dephosphorylation rate. At an AMP concentration of about 7 × 10−3 M, the rate enhancement amounts nearly to a factor of four, and this demonstrates that AMP can take over the role of a structuring ATP. What are the reasons that the mixed stacks are formed? They are twofold: (i) The concentration of the mixed stacks is favored by statistics simply because the AMP concentration is larger than that of ATP, and (ii) the charge repulsion in a [(AMP)(ATP)]6− stack is smaller than in an [(ATP)2]8− one.
That the M2+/N7 interaction in the stacks is crucial follows from experiments in which tubercidin 5′-monophosphate (TuMP2− = 7-deaza-AMP2−; see Figure 12), which lacks N7 as this is replaced by CH, is added to the Cu2+/ATP 2:1 system, and this reduces the dephosphorylation rate considerably [14,80]. This is despite the structural identity and the equality of self-stacking of TuMP2− and AMP2− [6].
The structural sensitivity of the reactive intermediate (Figure 10) is further demonstrated by the addition of 1,N6-ethenoadenosine 5′-monophosphate (ε-AMP; see Figure 12) to the Cu2+/ATP 2:1 system, which reduces the dephopsphorylation rate drastically: At [ε-AMP]tot = 4·10−3 M, the rate is practically zero (Figure 11). The stacking properties of ε-AMP2− are very similar to those of AMP [95,96], but the coordination of Cu2+ to the 1,10-phenenthroline-like (C6)N,N7 site [95,96] leads to a different orientation of the metal ion in space and thus inhibits the dephosphorylation.
Similarly, adenosine 5′-monophosphate N(1)-oxide (AMP·NO; Figure 12) loses under the experimental conditions (which are given in the legend of Figure 11) a proton, and Cu2+ binds thus to the ionized o-amino N-oxide group [97,98,99]. This means that the reactive stack involving N7 (Figure 10) cannot be formed; the Cu2+ ion is wrongly orientated in space, and AMP·NO is thus a strong inhibitor of the dephosphorylation reaction [14,80].
In fact, all ligands that have a larger affinity for M2+ than the metal ion has for N7 will destroy the reactive intermediates seen in Figure 10. Hence, if one adds to the Cu2+/ATP 2:1 system at pH 6.7 (see legend of Figure 11) ligands such as phosphate or ATP, the reactivity will decrease [14]. Ligands such as L-tryptophan (trp) or 2,2′-bipyridine (bpy) will form ternary complexes such as M(trp)(ATP)3− [100,101,102,103] or M(bpy)(ATP)2− [104,105,106,107,108] and thus inhibit the dephosphorylation as well; indeed, bpy can even protect ATP against hydrolysis [109]. Ligands such as dien (diethylenetriamine = 1,4,7-triazaheptane) or dpa [di(2-picolyl)amine] will even remove in part the metal ions such as Cu2+ from the reactive complexes and will thus inhibit dephosphorylation strongly (see Figure 9 in [14] or Figure 22 in [11]).
In the preceding paragraphs and in the context of discussing the impact on the dephosphorylation rate of the AMP derivatives seen in Figure 12, it became evident that the slightest alteration in the structure of the adenine residue reduces the dephosphorylation rate. Evidently, these altered AMP derivatives are unable to take over the structuring role in the dimers (left site of the dimers in Figure 10). It was therefore interesting to see whether the acyclic AMP analogue, i.e., PMEA (Figure 5), had any effect on the reactivity [94].
PMEA owns an unaltered adenine residue; hence, one condition for activity is fulfilled and the question is: Can the phosphonate group at the aliphatic ether chain form with a metal ion a bridge and take over the structuring role? The answer is YES, as Figure 11 shows. Addition of about eight equivalents of PMEA to the Cu2+/ATP 2:1 system with [Cu2+]tot = 2·10−3 M and [ATP]tot = 10−3 M increases the dephosphorylation rate by a factor of nearly seven. Hence, Cu3(ATP)(PMEA) or better Cu3(ATP)(PMEA)(OH) (see Section 6.1) is even a more reactive species than Cu3(ATP)(AMP)(OH), which has an enhancement factor of “only” about four. The reason for this observation is clearly the larger flexibility of the aliphatic ether chain, compared with the rigid ribose residue of AMP.

6.4. Matrix-Assisted Self-Association and Its Significance for Cell Organelles with High ATP Concentrations

This section (with small alterations) is reproduced with permission from our publication in Pure and Applied Chemistry [21]; copyright 2004, International Union of Pure and Applied Chemistry (IUPAC).
The results summarized in Section 6.2 show that the extent of aggregation (as described by the isodesmic model of non-cooperative self-association (Equations (18) and (19)) “is much affected by external conditions such as concentration, pH or the presence of metal ions ([11,23,80,91,92,93]), the neutralization of the negative charges at the phosphate groups being an important factor for stack formation. This has prompted a study of the effect of poly-α,l-lysine, p(Lys)n, on the self-stacking properties of ATP in D2O [110,111]. At pD 8.4 all the ε-amino groups of the side chains of p(Lys)n are protonated and carry a positive charge; i.e., ATP4− with its negative phosphate groups is expected to “line up” along the p ( H · Lys ) n n + matrix. Indeed, under the mentioned conditions with [p-Lys]side chains about 0.4 M and [ATP4−] ≤ 0.25 M, K = 11.5 ± 2.1 M−1 based on the isodesmic model [Equations (18) and (19)] [111]. This value should be compared with K = 1.3 ± 0.2 M−1 (Table 6) measured for ATP4− in the absence of any promoter. That indeed the positively charged side chains of p ( H · Lys ) n n + are responsible for the increased stacking tendency is proven by the fact that at pD 12, where the side chains of p(Lys)n are largely deprotonated, the stack-promoting effect of p(Lys)n has largely disappeared [111].”
“The indicated self-association of ATP is important for certain cell organelles; for example, for the 5-hydroxytryptamine organelles of mammalian blood platelets it was estimated that among other solutes [ATP] is 0.5 M, [ADP] 0.1 M, [Mg2+] 0.44 M, and [Ca2+] 0.11 M [112,113]. Similar ATP concentrations appear to hold, e.g., for cholinergic, neurosecretory or synaptic vesicles and granules (see Refs and Table 4 in [111]). Clearly, such vesicles should be osmotically unstable [114,115,116], yet based on the summarized results one may conclude that the high nucleotide concentrations in such vesicles can be handled by nature via self-association and aggregate formation [114,115,116]. Indeed, e.g., the chromaffin granules of bovine adrenal medulla contain next to nucleotides (0.2 M; mainly ATP) and metal ions (0.03 M; mainly Mg2+ and Ca2+) [117,118] also so-called chromogranines, i.e., low-molecular-weight proteins which might be used as a matrix for aggregate formation”. (For further “Conclusions Regarding Other Matrices” see [111]).
“In fact, one is tempted to speculate further: Considering that extracellular ATP is universally employed in cell-cell communication, particularly in synaptic transmission [119] (references therein) and if one recalls that electrons may migrate over long distances in DNA [120,121], one may propose that one way to achieve information transfer (and there are others) occurs in the following way. Assuming that, say six ATPs are lined up to form a stack covering a distance of approximately 20 Å one could imagine that at one end of the stack a metal ion (Fe2+, Mn2+, Cu+) is oxidized to a higher charged state (Fe3+, Mn3+, Cu2+) by a redox reaction and that this oxidized ion triggers hydrolysis of the triphosphate residue and that the electron travels through the stack to an acceptor at its other end [111]. This is depicted in a simplified fashion in Figure 13. That oxidation of a divalent metal ion (Mn2+) to a trivalent one (Mn3+) dramatically enhances the dephosphorylation rate of nucleoside 5′-triphosphates is known [11,73].”

7. Solution Structures of Mixed-Ligand Complexes Containing ATP4− and Related Ligands

7.1. Ligands Containing a Ribose Residue with a Phosphate Group

Nearly 50 years ago, it was discovered by UV-spectrophotometry that a charge-transfer band occurs in mixed ligand Cu2+ complexes formed by AMP2− or ATP4− and a heteroaromatic amine (arm)—i.e., 2,2′-bipyridine (bpy) or 1,10-phenanthroline (phen) [104,122]. This charge-transfer band was attributed to the formation of an intramolecular stack, as indicated in Equilibrium (22):
Molecules 27 02625 i002
From the fact that (i) the stability constant, e.g., of Cu(bpy)(ATP)2−, as determined from potentiomtric pH titrations on the one hand (log K Cu ( bpy ) ( ATP ) Cu ( bpy ) = 6.91 ± 0.15 [123]) and based on the charge-transfer band on the other (log K Cu ( bpy ) ( ATP ) Cu ( bpy ) = 6.96 ± 0.25 [104]) were practically identical and that (ii) also the extinction coefficient of the binary ATP4−/bpy adduct was similar to the one of the corresponding ternary complex indicated that the formation degree of the intramolecular stack in Cu(bpy)(ATP)2− is high; a tentative structure is shown in Figure 14 [104,124,125,126].
Of course, the stability of the binary ATP4−/bpy stacking adduct is low (log KCT = 0.91 ± 0.22) [104]; addition of a metal ion bridges the two ligands, and the stability increases, as is known from many examples [101,106,127,128,129,130]. Indeed, for example, the formation degree of the intramolecular stack in the ternary Cu(bpy)(AMP) and Cu(phen)(AMP) complexes is high—it amounts to 81 (±4) and 90 (±2)%, respectively [125]. Finally, solids of such complexes were isolated, and by X-ray crystal structure studies, also intramolecular stacking interactions were proven [107,108,131].
“The first mixed ligand complex containing ATP and an amino acid was one with tryptophanate, i.e., Zn(ATP)(Trp)3−. By 1H NMR shift experiments it was shown that an indole-adenine interaction takes place [132] which may be promoted by Zn2+ (see also [101]). Later, the position of the intramolecular Equilibrium (23)
Molecules 27 02625 i003
was determined with Zn2+ as metal ion and it was concluded [103] that the stacked species occurs with a formation degree of approximately 75%. Other metal ions were studied as well [100,103,133,134] and the occurrence of intramolecular stacks in M(ATP)(Trp)3− complexes” (for a tentative structure [101] see Figure 15) “was confirmed by several groups [135,136,137,138]. Related complexes containing phenyl and imidazole residues have also been investigated [124,139,140].”
“Moreover, it was shown that the isopropyl residue of leucinate is also able to undergo a hydrophobic interaction in M(ATP)(Leu)3− complexes; e.g., the formation degree of the “closed” species of the corresponding Mn2+ complex amounts to about 40% [103]. The occurrence of such species with a hydrophobic interaction was again proven by 1H NMR shift studies [103] and by comparisons of stability constants, which were mostly determined via potentiometric pH titrations [103]” (see also [126]).
“Quantification of the formation of intramolecular stacks in mixed ligand complexes containing a heteroaromatic amine (Arm) allowed the conclusion [126] that recognition between the adenine residue and amino acid side chains in mixed ligand complexes of the type M(ATP)(aa)3− (where aa = amino acetate derivative) decreases in the order (partly tentative): indole residue (tryptophan) > phenyl residue (phenylalanine) ≥ isopropyl residue (leucine) ≥ imidazole residue (histidine) > methyl residue (alanine).”
The preceding three paragraphs (with small alterations are reproduced with permission from our publication in Pure and Applied Chemistry [21]; copyright 2004, International Union of Pure and Applied Chemistry (IUPAC).

7.2. Acyclic Nucleoside Phosphonates: Cu(arm)(PMEA) as an Example with Intramolecular Stack Formation

In Section 6.2, we have seen that aromatic moieties can form stacks and that this holds quite generally [125]—e.g., also for 1,10-phenanthroline [141]. Therefore, it is important to point out that all results described below were obtained under conditions where self-association was negligible, and the results refer thus to monomeric species [45,141]. That in Cu(arm)(PMEA) complexes the adenine residue affects the stability of the complexes in a positive manner is immediately evident from Figure 16 [141].
With arm coordinated to Cu2+, there are only two equatorial coordination sites left for binding of a second ligand. Because in PMEA2− the phosphonate group is the primary binding site (Section 3), Equilibrium (5), which involves the ether O atom, is of relevance also for the ternary complexes, but a further direct equatorial coordination of the adenine residue is not possible. Hence, the only remaining way how the adenine residue can affect complex stability is via the formation of an intramolecular stack (Figure 17) [45,141].
Careful considerations of the steric situation [45,141] reveal that the formation of an intramolecular stack is only possible if the ether oxygen (Equation (5)) is released from the coordination sphere; only this provides enough flexibility. Aside from the open isomer, Cu(arm)(PMEA)op, there are two more isomers, one that involves the ether O atom, Cu(arm)(PMEA)cl/O and one with the intramolecular stack, Cu(arm)(PMEA)st; this then gives rise [45,141] to Equilibrium Scheme (24):
Molecules 27 02625 i004
To determine the formation degrees of the various species, one has to recall that the total stability increase is experimentally accessible and that the stability of the open isomers can be calculated via the reference lines seen in Figure 16; this leads to the following equations [45]:
K Cu ( arm ) ( PMEA ) op Cu ( arm ) = [ Cu ( arm ) ( PMEA ) op ] / ( [ Cu ( arm ) 2 + ] [ PMEA 2 ] )
KI/O = [Cu(arm)(PMEA)cl/O]/[Cu(arm)(PMEA)op]
KI/st = [Cu(arm)(PMEA)st]/[Cu(arm)(PMEA)op]
For the experimentally accessible equilibrium constant it holds,
K Cu ( arm ) ( PMEA ) Cu ( arm ) = [ Cu ( arm ) ( PMEA ) ] [ Cu ( arm ) 2 + ] [ PMEA 2 ] = ( [ Cu ( arm ) ( PMEA ) op ] + [ Cu ( arm ) ( PMEA ) cl / O ] + [ Cu ( arm ) ( PMEA ) st ] ) [ Cu ( arm ) 2 + ] [ PMEA 2 ] = K Cu ( arm ) ( PMEA ) op Cu ( arm ) + K I / O · K Cu ( arm ) ( PMEA ) op Cu ( arm ) + K I / st · K Cu ( arm ) ( PMEA ) op Cu ( arm ) = K Cu ( arm ) ( PMEA ) op Cu ( arm ) ( 1 + K I / O + K I / st )
and for the intramolecular interactions one obtains
K I = K I / tot = K Cu ( arm ) ( PMEA ) Cu ( arm ) K Cu ( arm ) ( PMEA ) op Cu ( arm ) = 10 log   Δ 1 = [ Cu ( arm ) ( PMEA ) int / tot ] [ Cu ( arm ) ( PMEA ) op ] = ( [ Cu ( arm ) ( PMEA ) cl / O ] + [ Cu ( arm ) ( PMEA ) st ] ) [ Cu ( arm ) ( PMEA ) op ] = K I / O + K I / st
Clearly, Cu(arm)(PMEA)int/tot in Equation (29) represents the sum of both species with an intramolecular interaction.
Since the overall stability constant as defined by Equation (28) is experimentally accessible (Table 7, Column 2), and because the stability constant of the open isomer (Equation (25)), in which only the phosphonate group is coordinated (Table 7, Column 3), can be calculated with the information given in Figure 16, the stability enhancement log Δ, as defined by Equation (12)), is obtained (Table 7, Column 4). These log Δ values can be used with Equation (29) to calculate KI/tot (Table 7, Column 5), and from these data follow (Equation (14)) the formation degrees in percentage for Cu(arm)(PMEA)int/tot (Table 7, Column 6).
Table 7. Quantification of the stability increase via log Δ (analogous to Equation (12)) for the Cu(arm)(PME) and Cu(arm)(PMEA) complexes and extent of the intramolecular chelate formation (Equation (14)) in the Cu(arm)(PME) complexes at 25 °C and I = 0.1 M (NaNO3); the results for the corresponding binary Cu(PME) complex are given for comparison.
Table 7. Quantification of the stability increase via log Δ (analogous to Equation (12)) for the Cu(arm)(PME) and Cu(arm)(PMEA) complexes and extent of the intramolecular chelate formation (Equation (14)) in the Cu(arm)(PME) complexes at 25 °C and I = 0.1 M (NaNO3); the results for the corresponding binary Cu(PME) complex are given for comparison.
M2+log   K M ( PME ) M
(Equations (8), and (28)) a
log   K M ( PME ) op M
(Analogous to Equation (25)) b
log ΔPME
(Equation (12)) c
KI
(Equations (6), (9) and (29)) c
% M(PME)cl
(Equation (14)) c,d
Cu2+3.73 ± 0.033.25 ± 0.060.48 ± 0.072.02 ± 0.4767 ± 5
Cu(bpy)2+3.86 ± 0.033.27 ± 0.070.59 ± 0.082.89 ± 0.6874 ± 5
Cu(phen)2+3.90 ± 0.043.28 ± 0.060.62 ± 0.073.17± 0.6976 ± 4
M2+log   K M ( PMEA ) M
[Equation (28)] a
log   K M ( PMEA ) op M
[Equation (25)] b
log ΔPMEA
[Equation (12)] c
KI = KI/tot
[Equation (29)] c
% Cu(arm)(PMEA)int/tot
[Equation (29); analog. to Equation (14)]
Cu(bpy)2+4.70 ± 0.023.22 ± 0.071.48 ± 0.0729.20 ± 4.8796.69 ± 0.53
Cu(phen)2+4.97 ± 0.033.23± 0.061.74 ± 0.0753.95 ± 8.8698.18 ± 0.29
a The constants were measured in aqueous solution via potentiometric pH titrations [44]. The constants are from [141]. For the error limits see Footnote ‘a’ of Table 2. b These constants were calculated with the pKa values 7.02 and 6.90 for H(PME) and H(PMEA), respectively (see footnotes in Table 2) and the straight-line equations given in Table 5 of Ref [141]; the corresponding error limits are 3 times the SD values (see Table 5 in Ref [141]). c The error limits for log Δ, KI, and % M(PME)cl were calculated according to the error propagation after Gauss; the same applies for the lower part of the table. d See Equilibrium (5). This table is reproduced (slight alterations were made) with permission from Table 8 in our publication in Coordination Chemistry Reviews [45]; copyright 1995, Elsevier Science S.A., Lausanne, Switzerland.
Based on the justified assumption that the stability of Cu(arm)(PME)cl equals that of Cu(arm)(PMEA)cl/O, one obtains KI/O for the latter species (Table 8, Column 5), and thus from Equation (29) KI/st follows (see Table 8, Column 6) because KI/tot is known. Therefore, the formation degrees of all three isomers that occur in Equilibrium Scheme (24) are known (Table 8, Columns 4, 7, and 8).
First, it should be noted that all three isomers according to Equilibrium Scheme (24) occur, but that the open isomer is a minority species for both Cu(arm)(PMEA)op complexes: The Cu(bpy)(arm)cl/O and the Cu(phen)(arm)cl/O species form to about 10% and 6%, respectively. The dominating species are clearly the one with the intramolecular stack, Cu(arm)(PMEA)st, which reach a formation degree of about 90% (Table 8, Column 8). This high formation degree of the stacked species suppresses the formation of the 5-membered chelate of Equilibrium (5); this contrasts with Cu(arm)(PME)cl, which reaches a formation degree of about 75% (Table 7, upper part, Column 6).
With regard to biological systems, one may note that it is expected that intramolecular stacks in ternary complexes do not play a significant role. However, what follows from the high formation degree of the Cu(arm)(PMEA)st species is that the adenine residue provides a good anchoring point. Such a point could be created by stack formation with another nucleobase, e.g., of a ribozyme, or by stack formation with an aromatic amino acid side chain of a protein. Such stacks can then give rise to the formation of reactive complexes, as seen in Figure 2; that is, the coordination of two metal ions to a triphosphate can occur either via M(α,β)-M(γ) or via M(α)-M(β,γ), leading to different products.

8. The Effect of a Change in Solvent Polarity on the Stability and Structure of Binary and Ternary Complexes

8.1. The Dianion of (Phosphonomethoxy)ethane (PME2−) as an Example of a Polar O-Ligand. The Properties of Cu(PME) Are Largely as Expected

A change in solvent polarity is expected to affect especially polar or ionic interactions. For PME2−, this means, as follows from Equilibrium (5), metal ion binding to the phosphonate group and the ether oxygen site is likely to be altered. Indeed, these interactions are expected to be affected by the so-called “effective” or “equivalent solution” dielectric constant. This constant is reduced in proteins [142,143,144] and in active site cavities of enzymes [145], if compared with bulk water; hence, the activity of water is reduced [146]. For example, the effective dielectric constant in the active sites of bovine carbonic anhydrase and carboxypeptidase A is 35 and <70, respectively [145], compared with about 78.5 for bulk water at 25 °C [147,148]. By applying 1,4-dioxane-water mixtures, one may expect to simulate to some extent the situation in an active site cavity; the dielectric constants are about 53 and 35 in water containing 30% or 50% (v/v) 1,4-dioxane, respectively [147,148].
To address the above situation, log K Cu ( R - PO 3 ) Cu versus p K H ( R - PO 3 ) H plots for water and water containing 30% or 50% (v/v) 1,4-dioxane were established [149,150]; the corresponding straight-line references are shown in Figure 18 [44,45,149,150].
It is evident that with increasing 1,4-dioxane concentration, the stability of the Cu2+-phosphonate complexes increases; at p K H ( R - PO 3 ) H = 7.0, this increase amounts to about 0.57 and 0.93 log units by the solvent change from water to water containing 30% or 50% (v/v) 1,4-dioxane; these are stability increases by factors of nearly 4 and 10, respectively. The properties of these solvents are listed in Columns 2 to 4 of Table 9. To conclude, the Cu2+-phosphonate binding is clearly facilitated by a decreasing solvent polarity, as one would have expected.
The data points for the Cu(PME) complex (●) are in all three solvents significantly above the reference lines. It should be recalled that the vertical differences to the straight-reference lines correspond to the stability enhancements log ΔCu(PME) (Equation (12)). These values are listed in Column 5 of Table 9; application of these data to Equation (10) provides KI (Table 9, Column 6) and thus with Equation (14) also the formation degrees of the closed species are obtained (Table 9, Column 7). It is interesting to note that the arithmetic means of log ΔCu(PME) (0.43 ± 0.07) and % Cu(PME)cl (62 ± 3) overlap within the error limits with the individual results and that this seems to indicate that the size of the stability enhancements log ΔCu(PME) and the formation degrees % Cu(PME)cl appear largely to be independent of the solvent. This may seem somewhat surprising, especially in the light of the overall stability increase in the Cu(PME) complexes (Figure 18) with the decreasing solvent polarity. On the other hand, this result implies that the solvent effect is about the same for both binding sites, the phosphonate group, and the ether oxygen of the aliphatic chain.

8.2. Nitrogen Donor Sites Are Especially Sensitive to Polarity Changes of the Solvent

With the observations of the oxygen-binding sites described in the preceding Section 8.1 in mind, it seems appropriate to consider now ligands with N-sites, as such sites are of relevance for biological systems as well. Some of such results are compiled in Table 10 [127].
For the planar 1,10-phenanthroline (phen) with its three fused aromatic 6-membered rings, one expects self-association. This is in aqueous solution very pronounced indeed (Table 10; Entry (1.a)), but it drops with 2,2′-bipyridine (bpy) to about one fourth due to the removal of the bridging aromatic ring (Table 10; Entry (1.b)). For both N-ligands, the self-stacking tendency decreases very significantly by the addition of 1,4-dioxane to the aqueous solution. This indicates that the ethylene units of dioxane solvate the aromatic rings of phen and bpy and suppress in this way self-stacking.
The results just indicated for the self-stacking process are reflected in the stability constants K ( arm ) ( ATP ) ( ATP ) observed for the binary stacking adducts (phen)(ATP)4− and (bpy)(ATP)4− (Table 10; Entries (2.a) and (2.b)). The adenine residue of ATP4− is well suited to form stacks (see Section 7). That the adduct with phen in D2O is about twice as stable as the one with bpy fits into the picture, as does the dramatic decrease in the stability constants K ( arm ) ( ATP ) ( ATP ) of the binary adducts by the addition of dioxane.
What happens if the two ligands ATP4− and arm (= phen or bpy) of Entries (2.a) and (2.b) are bridged by a metal ion? This question is best answered by considering the Δ log KCu/arm/ATP values that are defined in Equation (30) and that quantify Equilibrium (31) and therefore also the stability of the ternary complexes:
Δ   log   K Cu / arm / ATP = log   K Cu ( arm ) ( ATP ) Cu ( arm ) log   K Cu ( ATP ) Cu = log   K Cu ( ATP ) ( arm ) Cu ( ATP ) log   K Cu ( arm ) Cu
Cu ( arm ) 2 + + Cu ( ATP ) 2   Cu ( arm ) ( ATP ) 2 + Cu 2 +
For Δ log KCu/arm/ATP, one expects for statistical reasons a negative value: For the Jahn–Teller distorted Cu2+ and two bidentate ligands, Δ log KCu/statist ≅ –0.9 was estimated [127]. The results listed in Table 10 in Entries (3.a) and (3.b) for water are thus a surprise because they mean that Equilibrium (31) is shifted towards its right side. That the value for Cu(bpy)(ATP)2− is smaller than the one for Cu(phen)(ATP)2− is expected by having the results for Entries (2.a) and (2.b) in mind. Addition of dioxane affects the stability of the intramolecular stacks in Cu(arm)(ATP) in a negative manner, shifting Equilibrium (31) to the left, a result in accord with this given in Entries (2.a) and (2.b).
The two ligands of which Cu(arm)(ATP)2− is composed need to be considered also in an indirect manner. For phen and bpy, this has already been carried out (Table 10; Entries (1.a) and (1.b)), and for ATP4− (see also Table 6) it is made in the form of the Cu(ATP)2− complex (Table 10; Entry (4.b)). This result is compared with that for Cu(UTP)2−, UTP4− being known as having a very small tendency for self-association [127]. Addition of dioxane to Cu(UTP)2− (Table 10; Entry (4.a)) increases the stability of the complex, and the three available data sets fit in a log K Cu ( UTP ) Cu versus p K H ( UTP ) H plot on a straight line [127]. This means that the stability of Cu(UTP)2− increases, like the one of simple Cu2+-phosph(on)ate complexes [150,151], with decreasing solvent polarity and behaves thus in a normal manner.
This observation with Cu(UTP)2− contrasts strongly with the one for Cu(ATP)2−, the stability of which is practically independent of the solvent and also larger in all three solvents (Table 10; Entry (4.b)) than is expected on the basis of the basicity of the phosphate residue. This basicity is the same in water for UTP4− ( p K H ( UTP ) H = 6.46 ± 0.01) and ATP4− ( p K H ( ATP ) H = 6.49 ± 0.01) [127]; this identity of the basicity also holds for the mixed solvents [127]. Hence, the increased stability of Cu(ATP)2− in water must be due to the well-known [21,46,53,95] macrochelate formation with N7. This N7 interaction is inhibited by solvation of dioxane of the adenine residue hindering the coordination of Cu2+ to N7. Hence, the extent of the N7 interaction decreases with increasing amounts of dioxane, but in contrast, the Cu2+-phosphate binding strength increases. In other words, we have two opposing solvent effects that cancel each other and give rise to the observed (Table 10; Entry (4.b)) solvent-independent stability of the Cu(ATP)2− complex.

8.3. Competing Solvent Effects on N- versus O-Sites in the Same Ligand and Its Bearing on Complex Stability

After the conclusion in the preceding Section 8.2 that N- and O-sites face opposing solvent effects, it is of interest to see how a ligand with such ambivalent properties affects complex stability upon a change in solvent polarity. An ideal ligand in this respect is the well-studied AMP2− (Figure 3), which offers for coordination next to the primary binding site—hat is, the phosphate groups—the N7 site of the adenine residue that is involved in macrochelate formation (see Footnote ‘c’ in Table 11). The available data for the anion of the AMP analogue 9-[2-(phosphonomethoxy)ethyl]adenine (PMEA2−; Figure 3) are given for comparison. The values for the complexes of both ligands are listed in Table 11 [21,45,46,50,53,95,149,152,153].
The log stability constants of the Cu(AMP) complex (Table 11; Column 3) increase with increasing dioxane concentration, as expected for a ligand with a phosphate group (Section 3) as primary coordination site. The behavior of Cu(PMEA) is analogous: It is interesting to note that the stability enhancement by going from water to 50% (v/v) dioxane is identical for the two complexes—that is, it amounts for Cu(AMP) to 1.59 ± 0.04 and for Cu(PMEA) to 1.58 ± 0.08 log units. This confirms the close relationship between phosphate and phosphonate ligands. The corresponding differences in stability for the open isomers Cu(AD)op are within the error limits identical with the mentioned ones (Table 11; Column 4). Most interesting, however, are the percentages of the macrochelates, Cu(AMP)cl/N7 (Table 11; Column 7): They pass, beginning with a formation degree of 46 ± 10% in water, through a minimum of 9 ± 8% in 30 % (v/v) dioxane, and then go up again to 48 ± 5% in 50% dioxane. This reflects exactly the discussed opposing solvent effects and thus the ambivalent nature of the AMP2− ligand with its phosphate and N7 sites. The properties of Cu(PMEA) are analogous, though the stability differences are less pronounced (see the three bottom lines in Table 11).
As noted, weak interactions between metal ions and oxygen donor sites on the one hand and nitrogen donors on the other are very differently affected by a decreasing polarity of the solvent. This is nicely confirmed by the contents of Table 11 and illustrated in Figure 19, where the curves for % Cu(AMP)cl and % Cu(PMEA)cl pass through a minimum each, demonstrating again the ambivalent properties of these two ligands [41,47,154].
The ligand structures of the other Cu2+ complexes that appear in Figure 19 are pure O-donor ligands, as is seen in Figure 20 [155,156]. In all instances is the phosph(on)ate group the primary and thus the stability-determining binding site, but they all offer for weak interactions in their Cu2+ complexes either a carbonyl (DHAP2−, AnP2−, AcP2−) or a hydroxyl group (G1P2−).
In aqueous solution, the extent of this weak interaction varies strongly [47]: In Cu(AcP) the closed isomer amounts to 76 ± 4% and in Cu(AnP) to 56 ± 7% (also other M2+ reach high formation degrees), whereas in Cu(DHAP) and Cu(G1P) at best only traces (9 ± 13% for Cu(DHAP)) of the closed form exist. The reason for this different behavior is that in the case of Cu(AcP/AnP) 6-membered and in the case of Cu(DHAP/G1P) 7-membered chelates form. However, in the presence of dioxane, the formation degree of the chelates increases, and in 50% (v/v) 1,4-dioxane with Cu(DHAP) and Cu(G1P) 45 ± 4 and 42 ± 5%, respectively, of the complexes are present as 7-membered chelates [47]. For Cu(AnP) 80 ± 2% form 6-membered chelates [47]. In Figure 19, the corresponding plots are seen: In all three cases for which data exist, there is no minimum observed, but the curves increase with increasing dioxane concentration. This is as it should be for complexes of ligands that offer only O-donor sites.
To conclude, the discrimination that is achieved between N and O sites with a decreasing solvent polarity is fascinating and impressive.

9. Discussion

In Section 4.1 and Section 4.3, we have seen that macrochelate formation involving N7 and the phosphate group of an ANP is more pronounced in the Ni2+ complexes than in those containing Cu2+. The corresponding observation has been made earlier with Cu(AMP), which has a smaller stability enhancement compared with that of Ni(AMP); of course, the stability enhancement is a reflection of the N7 interaction—i.e., of the extent of macrochelate formation. These observations have been explained by statistical considerations [50]. Assuming Cu has a Jahn–Teller distorted octahedral coordination sphere [157], then there are only the two equatorial cis positions of the coordinated phosphate group left, which are able to form macrochelates by coordinating to N7; the trans position cannot be reached due to steric restrictions. In the case of Ni2+ with its regular octahedral coordination sphere and a bound phosphate group, four cis positions for N7 coordination remain; only the trans position is again not accessible. Hence, Ni2+ compared with Cu2+ is statistically favored by a factor of 2 (= 4/2) [50]; i.e., the stability enhancement is expected to be 0.3 log unit larger for Ni(AMP) than Cu(AMP). Indeed, the experimental results (see Equations (12) and (13)) are Δ log ΔNi/Cu/AMP = log ΔNi(AMP) − log ΔCu(AMP) = (0.54 ± 0.06) − (0.27 ± 0.08) = 0.27 ± 0.10 (from Column 4 in Table VI of [50]) and thus close to the statistical expectation.
In Section 4.1, we dealt with PME2AP, an ANP isomer in which the amino group is shifted to C2 making for steric reasons N3 less accessible but favoring N7. From the results listed in Table 3, we learned that macrochelate formation with N7 is more pronounced in Ni(PME2AP) (formation degree = % Ni(PME2AP)cl/N7 = 85 ± 4%) than in Cu(PME2AP) (73 ± 6%). Application of the above procedure gives Δ log ΔNi/Cu/PME2AP = Δ log ΔNi(PME2AP)Δ log ΔCu(PME2AP) = (0.82 ± 0.11) − (0.58 ± 0.11) = 0.24 ± 0.16 (data from Column 6 in Table 7 of [43]) indicating the correct trend. The complication is that in the present case not only two isomers (open closed N7) form, as with the M(AMP) complexes: Next to these, there is now the five-membered chelate involving the ether oxygen (see structure in Figure 5). To give the complete picture [43], for Ni(PME2AP) = 85 ± 4% (cl/N7), 4 ± 3% (cl/O), and 11 ± 3% (op), and also for Cu(PME2AP) = 73 ± 6% (cl/N7), 18 ± 6% (cl/O), and 9 ± 2% (op).
For PMEDAP2− (structure in Figure 5), the situation is similar, despite the fact that there are two amino groups: one at C2 and one at C6. Due to the increased basicity of the DAP residue [62], the steric hindrance is overcome, and N7 becomes a further binding site leading to macrochelate formation. This then gives rise again to three isomeric complexes: one with the macrochelate, one that involves the ether O atom, plus one open form [62]. The species distribution is as follows: Ni(PMEDAP) = 54 ± 10% (cl/N7), 13 ± 8% (cl/O), and 33 ± 6% (op), and for Cu(PMEDAP) = 43 ± 11% (cl/N7), 38 ± 11% (cl/O), and 19 ± 3% (op) (data from entries 2 and 3 of Table 6 in [62]). The larger formation degree of Ni(PMEDAP)cl/N7 confirms also here, despite the large errors, that Ni2+ is favored over Cu2+.
Above, we have discussed metal ion interactions with N7. In this context one should point out how sensitive these interactions are. For example, solvation of the adenine residue by 1,4-dioxane (in other words, in mixed solvents) inhibits the Cu2+–N7 interaction, and this also holds certainly for other metal ions. We assume that the ethylene groups of dioxane undergo hydrophobic interactions with the aromatic residue of the adenine parts and that this screens N7 (Section 8.2). On the other hand, addition of 1,4-dioxane to an aqueous solution reduces the dielectric constant of the solvent, and this favors polar or ionic interactions—e.g., between M2+ and a phosph(on)ate group. Such effects (at N and O sites) can give rise to opposing “reactions”, leading thus at a certain solvent mixture to a minimum of the interactions (see Table 11 and Figure 19). Such “solvent effects” can also be thought to be created by amino acid (aa) side chains of proteins. For example, for mixed ligand complexes of the type M(ATP)(aa)3−, it is known that the extent of the interaction decreases (in part tentatively) in the series aa = indole residue (tryptophan) > phenyl residue (phenylalanine) ≥ isopropyl residue (leucine) ≥ imidazole residue (histidine) > methyl residue (alanine) (see Section 7.1). Clearly, this can lead to selectivity via the binding strength but also via the orientation in space and therefore, one is in for surprises.
To the best of our knowledge, there is only one study of a larger RNA, where the Mg2+ ion binding affinity has been characterized in the presence of 20% 1,4-dioxane [158]. In domain 5 of a group II intron ribozyme, under such conditions, small structural changes were observed in a bulge known to bind Mg2+, a new coordination site was detected at a GU wobble pair, and in general the affinity of Mg2+ to the RNA increased by a factor of 1.5. More often than the application of 1,4-dioxane to reduce the solvent polarity, other small molecules, sugars as well as polymers such as dextran or PEG (poly-ethylene glycol) have been used to mimic in cellulo conditions in in vitro experiments [159]. These substances are called crowding agents, meaning they simulate the tight environment within a cell packed with proteins, biomolecules, and other “spectator” molecules, creating an effect of volume exclusion. However, these molecules also alter the solvent properties, such as the dielectric constant. Generally, such excluded volume effects drastically reduce the requirement for unnaturally high concentration of divalent metal ions needed under traditional in vitro conditions to stabilize the native fold of large RNAs. For example, in a group II intron, the addition of 10% PEG reduces the Mg2+ requirement for folding and catalytic activity by a factor of 2 and leads to an increase in catalytic activity by a factor of 4 when keeping the Mg2+ concentration constant [160]. Similar effects have been shown for many other classes of ribozymes [159]. In the case of the above-mentioned group II Intron, it was found that the increasing concentration of crowding agent reaches a maximum, after which catalytic activity is again reduced, probably because of too dense crowding becoming detrimental for a correct RNA fold [160].
Another such surprise occurs if one decides to compare the metal ion affinity of the phosphate residues of nucleotides; after all, these phosphate residues are the stability-determining binding sites of the nucleotides. To facilitate this overview, we plotted the logarithmic stability constants log K M ( N ) M , of the mono- (R-MP2−), di- (R-DP3−), and triphosphates (R-TP4−) (R always being a non-interacting residue) in an Irving–Williams-type fashion [161,162] (Figure 21).
Even though known for a long time, it is always a new surprise to observe that plotted data with phosphate ligands do not strictly follow the Irving–Williams sequence [161,162], the Mn2+ and the Zn2+ complexes are too stable [163,164] (see also pp. 364, 365 in [165]), or one could say that the complexes of the other metal ions are not stable enough; this may be a matter of outersphere complex formation. There is a large stability span of about 1.1 to 2.4 log units, the effect being especially pronounced for Cu2+, by going from the M(R-MP) complexes to those of M(R-DP), whereas the span from M(R-DP) to M(R-TP)2− amounts only to about 0.8 to 1.0 log unit (if the special case of Cu2+ is ignored); this indicates to us that in the M(R-MP) complexes, outersphere binding is of relevance [48,50].
That the phosphate residue plays the dominating part in the determination of complex stability of nucleotides is confirmed by the data that are inserted into Figure 21 for the adenine-nucleotide complexes; these are at least for the third transition metal ions always a touch more stable, which is the result of N7 binding. This small contribution to complex stability should not be misleading; as far as selectivity is concerned, the nucleobases play crucial roles, as we have repeatedly seen. Striking examples for such a selectivity are small ribozymes such as the TS [166] or hammerhead ribozyme: Both classes are most active in the presence of Mn2+, with the RzB and Schistosoma hammerhead ribozymes experiencing a 400-fold increase in activity over Mg2+ [167,168]. Analysis of the cleavage rates in the presence of metal ions along the Irving–Williams series showed a striking parallelism with the metal ion affinity towards the phosphate group [169]. This clearly indicates that the coordination to a phosphate moiety is the decisive factor in positioning the catalytic metal ion with the RNA structure. In contrast, for example in the glms ribozyme, named after its cofactor glucosamine-6-phosphate, highest activity is reached in the presence of Mg2+, whereas other earth alkaline metal ions and Mn2+, Co2+, Zn2+, and Cd2+ inhibit cleavage [170]. Inhibiting effects of divalent metal ions other than Mg2+ are also common in the large ribozymes like the group II Intron ribozyme Sc.ai5γ [171]. In such cases, the binding property of the metal ion to other functionals groups, such as the guanine or adenine N7, must play a crucial role too in correctly positioning the metal ion within the RNA structure or even in the catalytic core.

10. Conclusions and Outlook

In the preceding part of this review, we have repeatedly seen that metal ions interact with N7 of a purine residue and affect such the structure and function of DNA and RNA molecules [172]. A simple example is the coordination of Pt(II) to N7, which acidifies the (N1)H+ site in nucleoside phosphate derivates [173], though to varying extents bio-metals [174] are expected to have the same effect. However, metal ions bind in RNA and DNA not only to N7 but in rare cases also to N1 or N3—e.g., Mg(II) binds in an innersphere to the N3 position of a guanine residue within an adenine riboswitch aptamer (Mg105 to G6-N3) [175,176]. It is therefore of interest to quantify the basicity of the various N sites because commonly complex stability increases with basicity (Section 2). For adenosine (the structure is given in Figure 1), the micro acidity constants have been determined [177]; they are p k N 7 - N 1 · H N 7 - N 1 = 3.63 ± 0.02, p k H · N 7 - N 1 N 7 - N 1 = 2.15 ± 0.15, and p k H · N 3 - N 1 , N 7 N 3 - N 1 , N 7 = 1.5 ± 0.3. For example, the value for p k H · N 7 - N 1 N 7 - N 1 quantifies the basicity/acidity of the N7 site under conditions where the N1 site does not carry a proton. Analogously, the value of 1.5 for pkH·N3 holds for conditions when N1 and N7 are free. Hence, not only for N1 but also for N7 and N3 some metal ion affinity is expected, and this is of relevance for nucleic acids since adenine N7 is exposed in the major and N3 in the minor groove [178,179].
In this review, we have also dealt with acyclic nucleoside phosphonates, which have antiviral properties. An important representative is (S)-9-[3-hydroxy-2-(phosphonomethoxy)propyl]adenine (HPMPA), which is active against a range of DNA viruses [180]. The acid–base properties of this compound have recently been studied in detail [181]. The basicity of the adenine nitrogens decreases in the order N1 > N7 > N3, like it was observed for adenosine [177]. The micro acidity constant for +H·N7(HPMPA)N1 was estimated to be p k H · N 7 - N 1 N 7 - N 1 ca 3.5; the acid–base properties of PMEA are expected to be similar [181,182]. The metal ion-coordinating properties of HPMPA have also been studied [183], and it is interesting to compare these with those of (S)-1-[3-hydroxy-2-(phosphonomethoxy)propyl]cytosine (HPMPC); note, the cytosine residue does not affect complex stability [184]. It may be added that HPMPC, known as cidofovir [185], is used for the treatment of a variety of DNA virus infections [186]. One may also mention here Acyclovir, which after the addition of a phosphonate group gives the ANP2− = guanine(N9)–CH2–O–CH2–CH2 PO 3 2 ; one may consider this ANP2− as a relative of PMEA2− in which the adenine residue is replaced by a guanine one and where the "aliphatic" chain is also altered but the number of four chain links is kept. Acyclovir is used to treat herpes infections including outbreaks of genital herpes [187].
Coordination to N7 of purine moieties and associated macrochelate formation is also highly relevant and regularly observed in large nucleic acid structures—i.e., complex RNA molecules such as ribozymes, riboswitches, and other functional RNAs. In such cases, mostly Mg2+ is involved, with the Mg2+ being party dehydrated, forming innersphere as well as outersphere contacts to the nucleotide moieties. Depending on the local structure, the macrochelate is also formed, not necessarily within a given nucleotide—i.e., from the bridging phosphate to its own nucleobase—but more often between phosphate and N7 of two nucleotides far apart in the primary structure but close in space in the tertiary fold. Prominent examples of such phosphate-purine macrochelate formation by Mg2+ are structural metal ion binding motifs such as tandem GC and sheared GA base pairs. Large RNA structures such as the H. marismortui large ribosomal subunit [188] contain numerous of such Mg2+ binding sites, but they are also known from domains of the self-splicing yeast group II intron Sc.ai5γ [189,190,191], and the stem loop C of a Mg2+ riboswitch [192], as is summarized in Ref [163].
An interesting case is the self-activation of Cu2(ATP) via dimerization, [Cu2(ATP)]2 (Figure 10). In this dimer, one half of the complex provides the structure and the other half undergoes the reaction. It is especially fascinating that the structural part of one ATP4− can be taken over by AMP2− or even PMEA2− (Figure 10 and Figure 11). Nonetheless, the smallest structural alteration at the adenine residue (Figure 12) kills the activity. Though apparently simple, the situation is very sensitive, but it corresponds to the M(α,β)-M(γ) binding mode.
Furthermore, nucleoside 5′-triphosphates are in the form of their metal ion complexes important substrates for many reactions of biological relevance. There are two principal structural types—one with the metal ion coordinated as is indicated in the upper part of Figure 2, where the γ-phosphate group is transferred and we observe a kinase-type reaction that in the simplest case is a phosphoryl transfer to water; i.e., the NTP is hydrolyzed. If an NTP is exposed to an excess of metal ions, an M(α,β)-M(γ)-type coordination results. This is different for the structure shown in the middle part of Figure 2, which is of the M(α)-M(β,γ) type and leads to the transfer of a nucleotidyl group—that is, a polymerase-type reaction. In order to achieve this binding mode, the two metal ions need to be anchored by the enzyme; this fixation process can occur via amino acid side chains (such as carboxylate or imidazole groups) of the protein, but fixation via stacking is also possible.
This alternate positioning of two Mg2+ along the triphosphate chain is used by two prominent classes of enzymes: In a kinase-type reaction (i.e., the transfer of the γ-phosphate, a typical “phosphoryl-transfer” type of Mg2+ coordination is observed), for example, in a phosphoenolpyruvate carboxykinase, the two divalent metal ions (Mn2+ and Mg2+ in this case) coordinate to the ATP triphosphate chain in a M(α,β,γ)-M(γ)-fashion [24]. In contrast, a polymerase-type coordination is observed in T7 RNA polymerase [193] as well as in a high-fidelity DNA polymerase RB69 [26]. In addition, the employment of two metal ions is also known from ribozymes, where a classical two-metal ion mechanism is used to cleave the RNA chain at a specific phosphodiester—e.g., in group I introns [194,195] and in group II introns [196].
For a metal ion to be shifted along a triphosphate chain into the reactive position, Mg2+ appears to be ideal because in the initial 1:1 complex it seems to be present as a mixture of β,γ-bidentate and α,β,γ-tridentate species [197]. Mg2+ is by itself not very reactive, but in combination with another metal ion it is quite useful (Figure 8).

Funding

This research was funded by the Swiss National Science Foundation grant number 200020_192153 (to RKOS).

Acknowledgments

The support of our work (A.S. and H.S.) by the Department of Chemistry of the University of Basel, Switzerland, is gratefully acknowledged, as is the support for R.K.O.S. by the Swiss National Science Foundation and the University of Zürich.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations and Definitions

aaamino acid
AcP2−acetyl phosphate (Figure 20)
AIDSacquired immune deficiency syndrome
AMP·NO2−adenosine 5′-monophosphate N(1)-oxide (Figure 12)
AMP2−adenosine 5′-monophosphate (Figure 3)
ε-AMP2−1,N6-ethenoadenosine 5′-monophosphate (Figure 12)
AnP2−acetonylphosphonate (Figure 20)
armheteroaromatic nitrogen base, e.g., bpy or phen
ATP4−adenosine 5′-triphosphate
bpy2,2-bipyridine
DHAP2−dihydroxyacetone phosphate (Figure 20)
diendiethylenetriamine = 1,4,7-triazaheptane
DNAdeoxyribonucleic acid
dPMEA2−dianion of 9-(4-phosphonobutyl)adenine = 3′-deoxa-PMEA2−
G1P2−glycerol 1-phosphate (= α-glycerophosphate; in many biochemistry texts also designated as glycerol 3-phosphate) (Figure 20)
HIVhuman immunodeficiency virus
L2−any phosph(on)ate ligand (R- PO 3 2 )
M2+any divalent metal ion (in a few instances also Cu(bpy)2+ and Cu(phen)2+ are represented by this abbreviation)
Nnucleophile
NDP3−nucleoside 5′-diphosphate
NMP2−nucleoside 5′-monophosphate
NTP4−nucleoside 5′-triphosphate
PEEA2−dianion of 9-[2-(2-phosphonoethoxy)ethyl]adenine
phen1,10-phenanthroline
PME2−dianion of (phosphonomethoxy)ethane (=ethoxymethanephosphonate) [Equation (5)]
PMEA2−dianion of 9-[2-(phosphonomethoxy)ethyl]adenine (Figure 3)
RNAribonucleic acid
TuMP2−tubercidin 5′-monophosphate = 7-deaza-AMP2− (Figure 12)
Several additional definitions are given in the legends to Figure 1, Figure 2, Figure 4 and Figure 5.—Species that are given in the text without a charge either do not carry one or represent the species in general (i.e., independent of their degree of protonation); which of the two versions applies is always clear from the context.

References

  1. Sigel, A.; Sigel, H. (Eds.) Interactions of Metal Ions with Nucleotides, Nucleid Acids, and Their Constituents. In Metal Ions in Biological Systems; Marcel Dekker, Inc.: New York, NY, USA; Basel, Switzerland; Hong Kong, China, 1996; Volume 32, pp. 1–814. [Google Scholar]
  2. Sigel, A.; Sigel, H.; Sigel, R.K.O. (Eds.) Structural and Catalytic Roles of Metal Ions in RNA. In Metal Ions in Life Sciences; Royal Society of Chemistry: Cambridge, UK, 2011; Volume 9, pp. 1–391. [Google Scholar]
  3. Sigel, A.; Sigel, H.; Sigel, R.K.O. (Eds.) Interplay Between Metal Ions and Nucleic Acids. In Metal Ions in Life Sciences; Springer: Dordrecht, The Netherlands, 2012; Volume 10, pp. 1–351. [Google Scholar]
  4. Aoki, K. General Conclusions from Solid State Studies of Nucleotide-Metal Ion Complexes. Met. Ions Biol. Syst. 1996, 32, 91–134. [Google Scholar]
  5. Martin, R.B.; Mariam, Y.H. Interactions between Metal Ions and Nucleic Bases, Nucleosides, and Nucleotides in Solution. Met. Ions Biol. Syst. 1979, 8, 57–124. [Google Scholar]
  6. Tribolet, R.; Sigel, H. Self-Association and Protonation of Adenosine 5′-Monophosphate in Comparison with Its 2′- and 3′-Analogues and Tubercidin 5′-Monophosphate (7-deaza-AMP). Eur. J. Biochem. 1987, 163, 353–363. [Google Scholar] [CrossRef] [PubMed]
  7. Davies, D.B.; Rajani, P.; Sadikot, H. Determination of Glycosidic Bond Conformations of Pyrimidine Nucleosides and Nucleotides Using Vicinal Carbon–Proton Coupling Constants. J. Chem. Soc. Perkin Trans. 1985, 2, 279–285. [Google Scholar] [CrossRef]
  8. Pedley, A.M.; Benkovic, S.J. A New View into the Regulation of Purine Metabolism—The Purinosome. Trends Biochem. Sci. 2017, 42, 141–154. [Google Scholar] [CrossRef] [Green Version]
  9. Kanehisa, M.; Goto, S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  10. Sigel, H.; Amsler, P.E. Hydrolysis of Nucleoside Phosphates. VI. On the Mechanism of the Metal Ion Promoted Dephosphorylation of Purine Nucleoside 5′-Triphosphates. J. Am. Chem. Soc. 1976, 98, 7390–7400. [Google Scholar] [CrossRef]
  11. Sigel, H. Mechanistic Aspects of the Metal Ion Promoted Hydrolysis of Nucleoside 5′-Triphosphates (NTPs). Coord. Chem. Rev. 1990, 100, 453–539. [Google Scholar] [CrossRef]
  12. Mildvan, A.S. Role of Magnesium and Other Divalent Cations in ATP-Utilizing Enzymes. Magnesium 1987, 6, 28–33. [Google Scholar]
  13. Sträter, N.; Lipscomb, W.N.; Klabunde, T.; Krebs, B. Two-Metal Ion Catalysis in Enzymatic Acyl- and Phosphoryl-Transfer Reactions. Angew. Chemie Int. Ed. Engl. 1996, 35, 2024–2055. [Google Scholar] [CrossRef]
  14. Sigel, H.; Hofstetter, F.; Martin, R.B.; Milburn, R.M.; Scheller-Krattiger, V.; Scheller, K.H. General Considerations on Transphosphorylations: Mechanism of the Metal Ion Facilitated Dephosphorylation of Nucleoside 5′-Triphosphates, Including Promotion of ATP Dephosphorylation by Addition of Adenosine 5′-Monophosphate. (Hydrolysis of Nucleoside Phosphates. Part 8). J. Am. Chem. Soc. 1984, 106, 7935–7946. [Google Scholar]
  15. Steitz, T.A. DNA Polymerases: Structural Diversity and Common Mechanisms. J. Biol. Chem. 1999, 274, 17395–17398. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Brautigam, C.A.; Steitz, T.A. Structural and Functional Insights Provided by Crystal Structures of DNA Polymerases and Their Substrate Complexes. Curr. Opin. Struct. Biol. 1998, 8, 54–63. [Google Scholar] [CrossRef]
  17. Pelletier, H.; Sawaya, M.R.; Kumar, A.; Wilson, S.H.; Kraut, J. Structures of Ternary Complexes of Rat DNA Polymerase Beta, a DNA Template-Primer, and ddCTP. Science 1994, 264, 1891–1903. [Google Scholar] [CrossRef]
  18. Pelletier, H. Polymerase Structures and Mechanism. Science 1994, 266, 2025–2026. [Google Scholar] [CrossRef] [Green Version]
  19. Pelletier, H.; Sawaya, M.R.; Wolfle, W.; Wilson, S.H.; Kraut, J. Crystal Structures of Human DNA Polymerase Beta Complexed with DNA: Implications for Catalytic Mechanism, Processivity, and Fidelity. Biochemistry 1996, 35, 12762–12777. [Google Scholar] [CrossRef]
  20. Sigel, H.; Song, B.; Blindauer, C.A.; Kapinos, L.E.; Gregáň, F.; Prónayová, N. Why is the Antiviral Nucleotide Analogue 9-[2-(Phosphonomethoxy)ethyl]adenine in Its Diphosphorylated Form (PMEApp4−) Initially a Better Substrate for Polymerases than (2′-Deoxy)adenosine 5′-Triphosphate (dATP4−/ATP4−)? Considerations on the Mechanism of Polymerases. Chem. Commun. 1999, 743–744. [Google Scholar] [CrossRef]
  21. Sigel, H. Adenosine 5′-Triphosphate (ATP4−): Aspects of the Coordination Chemistry of a Multi-Talented Biological Substrate. Pure Appl. Chem. 2004, 76, 375–388. [Google Scholar] [CrossRef] [Green Version]
  22. Sigel, H. Metal Ion Complexes of Antivirally Active Nucleotide Analogues. Conclusions Regarding Their Biological Action. Chem. Soc. Rev. 2004, 33, 191–200. [Google Scholar] [CrossRef]
  23. Sigel, H.; Griesser, R. Nucleoside 5′-Triphosphates: Self-Association, Acid-Base, and Metal Ion-Binding Properties in Solution. Chem. Soc. Rev. 2005, 34, 875–900. [Google Scholar] [CrossRef]
  24. Tari, L.W.; Matte, A.; Goldie, H.; Delbaere, L.T.J. Mg2+–Mn2+ Clusters in Enzyme-Catalyzed Phosphoryl-Transfer Reactions. Nature Struct. Biol. 1997, 4, 990–994. [Google Scholar] [CrossRef] [PubMed]
  25. Sigel, H. Metal Ion-Binding Properties of the Antiviral Nucleotide Analogue 9-[2-(Phosphonomethoxy)ethyl]adenine (PMEA). Why Is Its Diphosphorylated Form, PMEApp4−, Initially a Better Substrate for Nucleic Acid Polymerases than (2′-Deoxy)-Adenosine 5′-Triphosphate (dATP4−/ATP4−)? Pure Appl. Chem. 1999, 71, 1727–1740. [Google Scholar]
  26. Freisinger, E.; Grollman, A.P.; Miller, H.; Kisker, C. Lesion (In)tolerance Reveals Insights into DNA Replication Fidelity. EMBO J. 2004, 23, 1494–1505. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Barré-Sinoussi, F.; Chermann, J.C.; Rey, F.; Nugeyre, M.T.; Chamaret, S.; Gruest, J.; Dauguet, C.; Axler-Blin, C.; Vézinet-Brun, F.; Rouzioux, C.; et al. Isolation of a T-Lymphotropic Retrovirus from a Patient at Risk for Acquired Immune Deficiency Syndrome (AIDS). Science 1983, 220, 868–871. [Google Scholar] [CrossRef] [Green Version]
  28. Popovic, M.; Sarngadharan, M.G.; Read, E.; Gallo, R.C. Detection, Isolation, and Continuous Production of Cytopathic Retroviruses (HTLV-III) from Patients with AIDS and Pre-AIDS. Science 1984, 224, 497–500. [Google Scholar] [CrossRef]
  29. Gallo, R.C.; Montagnier, L. AIDS in 1988. Sci. Am. 1988, 259, 41–48. [Google Scholar] [CrossRef]
  30. De Clercq, E. New Perspectives for the Treatment of HIV Infections. Coll. Czech. Chem. Commun. 1998, 63, 449–479. [Google Scholar] [CrossRef]
  31. Song, B.; Sigel, R.K.O.; Sigel, H. Acid-Base Properties of Adenosine 5′-O-Thiomonophosphate in Aqueous Solution. Chem. Eur. J. 1997, 3, 29–33. [Google Scholar] [CrossRef]
  32. Sigel, R.K.O.; Song, B.; Sigel, H. Stabilities and Structures of Metal Ion Complexes of Adenosine-5′-O-Thiomonophosphate (AMPS2−) in Comparison with Those of Its Parent Nucleotide (AMP2−) in Aqueous Solution. J. Am. Chem. Soc. 1997, 119, 744–755. [Google Scholar] [CrossRef]
  33. Sigel, A.; Operschall, B.P.; Sigel, R.K.O.; Sigel, H. Metal Ion Complexes of Nucleoside Phosphorothioates Reflecting the Ambivalent Properties of Lead(II). New J. Chem. 2018, 42, 7551–7559. [Google Scholar] [CrossRef] [Green Version]
  34. Holý, A.; De Clercq, E.; Votruba, I. Phosphonylmethyl Ethers of Nucleosides and Their Acyclic Analogues. ACS Symp. Ser. 1989, 401, 51–71. [Google Scholar]
  35. Holý, A. Acyclic Nucleoside and Nucleotide Analogues. A New Class of Antivirals. Il Farmaco 1991, 46, 141–146. [Google Scholar] [PubMed]
  36. Blindauer, C.A.; Holý, A.; Dvořáková, A.; Sigel, H. Solution Properties of Antiviral Adenine-Nucleotide Analogues. The Acid-Base Properties of 9-[2-(Phosphonomethoxy)ethyl]adenine (PMEA) and of Its N1, N3, and N7 Deaza Derivatives in Aqueous Solution. J. Chem. Soc. Perkin Trans. 2 1997, 11, 2353–2363. [Google Scholar] [CrossRef]
  37. Schwalbe, C.H.; Thomson, W.; Freeman, S. Structural Studies on Bio-Active Molecules. Part 17. Crystal Structure of 9-(2′-phosphonomethoxyethyl)adenine (PMEA). J. Chem. Soc. Perkin Trans. 1 1991, 1348–1349. [Google Scholar] [CrossRef]
  38. Aoki, K.; Murayama, K. Nucleic Acid-Metal Ion Interactions in the Solid State. Met. Ions Life Sci. 2021, 10, 43–102. [Google Scholar]
  39. Sigel, H.; Zuberbühler, A.D.; Yamauchi, O. Comments on Potentiometric pH Titrations and the Relationship between pH-Meter Reading and Hydrogen Ion Concentration. Anal. Chim. Acta 1991, 255, 63–72. [Google Scholar] [CrossRef]
  40. Saha, A.; Saha, N.; Ji, L.-N.; Zhao, J.; Gregáň, F.; Sajadi, S.A.A.; Song, B.; Sigel, H. Stability of Metal Ion Complexes Formed with Methyl Phosphate and Hydrogen Phosphate. J. Biol. Inorg. Chem. 1996, 1, 231–238. [Google Scholar] [CrossRef]
  41. Sigel, H.; Da Costa, C.P.; Song, B.; Carloni, P.; Gregáň, F. Stability and Structure of Metal Ion Complexes Formed in Solution with Acetyl Phosphate and Acetonylphosphonate: Quantification of Isomeric Equilibria. J. Am. Chem. Soc. 1999, 121, 6248–6257. [Google Scholar] [CrossRef]
  42. Fernández-Botello, A.; Griesser, R.; Holý, A.; Moreno, V.; Sigel, H. Acid-Base and Metal Ion-Binding Properties of 9-[2-(2-Phosphonoethoxy)ethyl]adenine (PEEA), a Relative of the Antiviral Nucleotide Analogue 9-[2-(Phosphonomethoxy)ethyl]adenine (PMEA). An Exercise on the Quantification of Isomeric Complex Equilibria in Solution. Inorg. Chem. 2005, 44, 5104–5117. [Google Scholar]
  43. Fernández-Botello, A.; Operschall, B.P.; Holý, A.; Moreno, V.; Sigel, H. Metal Ion-Binding Properties of 9-[2-(Phosphonomethoxy)ethyl]-2-aminopurine (PME2AP), an Isomer of the Antiviral Nucleotide Analogue 9-[2-(Phosphonomethoxy)ethyl]adenine (PMEA). Steric Guiding of Metal Ion-Coordination by the Purine-Amino Group. Dalton Trans. 2010, 39, 6344–6354. [Google Scholar] [CrossRef]
  44. Sigel, H.; Chen, D.; Corfù, N.A.; Gregáň, F.; Holý, A.; Strašák, M. Metal-Ion-Coordinating Properties of Various Phosphonate Derivatives, Including 9-[2-(Phosphonylmethoxy)ethyl]adenine (PMEA)—An Adenosine Monophosphate (AMP) Analogue with Antiviral Properties. Helv. Chim. Acta 1992, 75, 2634–2656. [Google Scholar] [CrossRef]
  45. Sigel, H. Intramolecular Equilibria in Metal Ion Complexes of Artificial Nucleotide Analogues with Antiviral Properties. A Case Study. Coord. Chem. Rev. 1995, 144, 287–319. [Google Scholar] [CrossRef]
  46. Sigel, H.; Song, B. Solution Structures of Nucleotide Metal Ion Complexes. Isomeric Equilibria. Met. Ions Biol. Syst. 1996, 32, 135–205. [Google Scholar]
  47. Sigel, H.; Kapinos, L.E. Quantification of Isomeric Equilibria for Metal Ion Complexes Formed in Solution by Phosphate or Phosphonate Ligands with a Weakly Coordinating Second Site. Coord. Chem. Rev. 2000, 200–202, 563–594. [Google Scholar] [CrossRef]
  48. Massoud, S.S.; Sigel, H. Metal Ion Coordinating Properties of Pyrimidine-Nucleoside 5′-Monophosphates (CMP, UMP, TMP) and of Simple Phosphate Monoesters, Including D-Ribose 5′-Monophosphate. Establishment of Relations Between Complex Stability and Phosphate Basicity. Inorg. Chem. 1988, 27, 1447–1453. [Google Scholar] [CrossRef]
  49. Martin, R.B.; Sigel, H. Quantification of Intramolecular Ligand Equilibria in Metal-Ion Complexes. Comments Inorg. Chem. 1988, 6, 285–314. [Google Scholar] [CrossRef]
  50. Sigel, H.; Massoud, S.S.; Tribolet, R. Comparison of the Metal Ion Coordinating Properties of Tubercidin 5′-Monophosphate (7-Deaza-AMP) with Those of Adenosine 5′-Monophosphate (AMP) and 1,N6-Ethenoadenosine 5′-Monophosphate (ε-AMP). Definite Evidence for Metal Ion-Base Backbinding to N-7 and Extent of Macrochelate Formation in M(AMP) and M(ε-AMP). J. Am. Chem. Soc. 1988, 110, 6857–6865. [Google Scholar]
  51. Blindauer, C.A. Untersuchungen zu Stabilität, Struktur und Reaktivität von Metallionenkomplexen Einiger Antiviraler Nukleotid-Analoga; Inauguraldissertation. Ph.D. Thesis, Universität Basel, Basel, Switzerland, 1998; pp. 1–235. [Google Scholar]
  52. Diebler, H. Stability and Structure of Complexes of Transition Metal Ions with Nucleotides and Related Compounds. J. Mol. Catal. 1984, 23, 209–217. [Google Scholar] [CrossRef]
  53. Sigel, H. Isomeric Equilibria in Complexes of Adenosine 5′-Triphosphate with Divalent Metal Ions. Solution Structures of M(ATP)2− Complexes. Eur. J. Biochem. 1987, 165, 65–72. [Google Scholar] [CrossRef]
  54. Sigel, H. Metal-Nucleotide Interactions. ACS Symp. Ser. 1989, 402, 159–204. [Google Scholar]
  55. Sigel, H. Have Adenosine 5′-Triphosphate (ATP4−) and Related Purine-Nucleotides Played a Role in Early Evolution? ATP, Its Own ‘Enzyme’ in Metal Ion Facilitated Hydrolysis! Inorg. Chim. Acta 1992, 198, 1–11. [Google Scholar] [CrossRef]
  56. Sigel, H. Interactions of Metal Ions with Nucleotides and Nucleic Acids and Their Constituents. Chem. Soc. Rev. 1993, 22, 255–267. [Google Scholar] [CrossRef]
  57. Holý, A.; Günter, J.; Dvořákova, H.; Masojídkova, M.; Andrei, G.; Snoeck, R.; Balzarini, J.; De Clercq, E. Structure−Antiviral Activity Relationship in the Series of Pyrimidine and Purine N-[2-(2-Phosphonomethoxy)ethyl] Nucleotide Analogues. 1. Derivatives Substituted at the Carbon Atoms of the Base. J. Med. Chem. 1999, 42, 2064–2086. [Google Scholar] [CrossRef] [PubMed]
  58. Fernández-Botello, A.; Holý, A.; Moreno, V.; Operschall, B.P.; Sigel, H. Intramolecular π-π Stacking Interactions in Aqueous Solution in Mixed-Ligand Copper(II) Complexes Formed by Heteroaromatic Amines and the Nucleotide Analogue 9-[2-(Phos phonomethoxy)ethyl]-2-aminopurine (PME2AP), an Isomer of the Antivirally Active 9-[2-(Phosphonomethoxy)ethyl]adenine (PMEA). (Ternary Complexes in Solution. Part 69). Inorg. Chim. Acta 2009, 362, 799–810. [Google Scholar]
  59. De Clercq, E. The Holý Trinity: The Acyclic Nucleoside Phosphonates. Adv. Pharmacol. 2013, 67, 293–316. [Google Scholar]
  60. Holý, A. Phosphonomethoxyalkyl Analogs of Nucleotides. Curr. Pharm. Des. 2003, 9, 2567–2592. [Google Scholar] [CrossRef]
  61. Gómez-Coca, R.B.; Sigel, A.; Operschall, B.P.; Holý, A.; Sigel, H. Solution Properties of Metal Ion Complexes Formed with the Antiviral and Cytostatic Nucleotide Analogue 9-[2-(Phosphonomethoxy)ethyl]-2-amino-6-dimethylaminopurine (PME2A6DMAP). Can. J. Chem. 2014, 92, 771–780. [Google Scholar] [CrossRef]
  62. Blindauer, C.A.; Sjåstad, T.I.; Holý, A.; Sletten, E.; Sigel, H. Aspects of the Co-ordination Chemistry of the Antiviral Nucleotide Analogue, 9-[2-(Phosphonomethoxy)ethyl]-2,6-diaminopurine (PMEDAP). J. Chem. Soc., Dalton Trans. 1999, 3661–3671. [Google Scholar] [CrossRef]
  63. Blindauer, C.A.; Holý, A.; Sigel, H. Metal Ion-Binding Properties of the Nucleotide Analogue 1-[2-(Phosphonomethoxy)ethyl]cytosine (PMEC) in Aqueous Solution. Coll. Czech. Chem. Commun. 1999, 64, 613–632. [Google Scholar] [CrossRef]
  64. Sigel, H. On the Metal Ion-Coordinating Properties of the Antiviral Nucleotide Analogue 9-[2-(Phosphonomethoxy)ethyl]adenine and Related Compounds. J. Indian Chem. Soc. 1997, 74, 261–271. [Google Scholar] [CrossRef]
  65. Blindauer, C.A.; Emwas, A.H.; Holý, A.; Dvořákova, H.; Sletten, E.; Sigel, H. Complex Formation of the Antiviral 9-[2-(Phosphonomethoxy)ethyl]adenine (PMEA) and of Its N1, N3, and N7 Deaza Derivatives with Copper(II) in Aqueous Solution. Chem. Eur. J. 1997, 3, 1526–1536. [Google Scholar] [CrossRef]
  66. Gómez-Coca, R.B.; Kapinos, L.E.; Holý, A.; Vilaplana, R.A.; González-Vílchez, F.; Sigel, H. Metal Ion-Binding Properties of 9-(4-Phosphonobutyl)adenine (dPMEA), a Sister Compound of the Antiviral Nucleotide Analogue 9-[2-(Phosphonomethoxy)ethyl]adenine (PMEA), and Quantification of the Equilibria Involving Four Cu(PMEA) Isomers. J. Chem. Soc. Dalton Trans. 2000, 13, 2077–2084. [Google Scholar] [CrossRef]
  67. Holý, A.; Votruba, I.; Merta, A.; Černý, J.; Veselý, J.; Vlach, J.; Šedivá, K.; Rosenberg, I.; Otmar, M.; Hřebabecký, H.; et al. Acyclic Nucleotide Analogues: Synthesis, Antiviral Activity and Inhibitory Effects on some Cellular and Virus-Encoded Enzymes in vitro. Antivir. Res. 1990, 13, 295–311. [Google Scholar] [CrossRef]
  68. Neyts, J.; Stals, F.; Bruggeman, C.; De Clercq, E. Activity of the Anti-HIV Agent 9-(2-phosphonyl-methoxyethyl)-2,6-diaminopurine against Cytomegalovirus in vitro and in vivo. Eur. J. Clin. Microbiol. Infect. Dis. 1993, 12, 437–446. [Google Scholar] [CrossRef]
  69. Franek, F.; Holý, A.; Votruba, I.; Eckschlager, T. Acyclic Nucleotide Analogues Suppress Growth and Induce Apoptosis in Human Leukemia Cell Lines. Int. J. Oncol. 1999, 14, 745–797. [Google Scholar] [CrossRef]
  70. Kramata, P.; Votruba, I.; Otová, B.; Holý, A. Different Inhibitory Potencies of Acyclic Phosphonomethoxyalkyl Nucleotide Analogs toward DNA Polymerases alpha, delta and epsilon. Mol. Pharmacol. 1996, 49, 1005–1011. [Google Scholar]
  71. Valeriánová, M.; Otová, B.; Bílá, V.; Hanzalová, J.; Votruba, I.; Holý, A.; Eckschlager, T.; Krejčí, O.; Trka, J. PMEDAP and Its N6-Substituted Derivatives: Genotoxic Effect and Apoptosis in in vitro Conditions. J. Anticancer Res. 2003, 23, 4933–4939. [Google Scholar]
  72. Gómez-Coca, R.B.; Blindauer, C.A.; Sigel, A.; Operschall, B.P.; Holý, A.; Sigel, H. Extent of Intramolecular π-Stacks in Aqueous Solution in Mixed-Ligand Copper(II) Complexes Formed by Heteroaromatic Amines and Several 2-Aminopurine Derivatives of the Antivirally Active Nucleotide Analogue 9-[2-(Phosphonomethoxy)ethyl]adenine (PMEA). (Ternary Complexes. Part 71). Chem. Biodivers. 2012, 9, 2008–2034. [Google Scholar]
  73. Amsler, P.E.; Sigel, H. Hydrolysis of Nucleoside Phosphates. V. Comparison of the Metal-Ion-Promoted Dephosphorylation of the 5′-Triphosphates of Adenosine, Inosine, Guanosine and Cytidine by Mn2+, Ni2+ and Zn2+ in Binary and Ternary Complexes. Eur. J. Biochem. 1976, 63, 569–581. [Google Scholar] [CrossRef]
  74. Sigel, H.; Tribolet, R.; Malini-Balakrishnan, R.; Martin, R.B. Comparison of the Stabilities of Monomeric Metal Ion Complexes Formed with Adenosine 5′-Triphosphate (ATP) and Pyrimidine-Nucleoside 5′-Triphosphates (CTP, UTP, TTP) and Evaluation of the Isomeric Equilibria in the Complexes of ATP and CTP. Inorg. Chem. 1987, 26, 2149–2157. [Google Scholar] [CrossRef]
  75. Sigel, H.; Sigel, A. (Eds.) Zinc and Its Role in Biology and Nutrition. In Metal Ions in Biological Systems; Marcel Dekker, Inc.: New York, NY, USA, 1983; Volume 15, pp. 1–493. [Google Scholar]
  76. Sigel, H.; Martin, R.B. The Colourless “Chameleon” or the Peculiar Properties of Zn2+ in Complexes in Solution. Quantification of Equilibria Involving a Change of the Coordination Number of the Metal Ion. Chem. Soc. Rev. 1994, 23, 83–91. [Google Scholar] [CrossRef]
  77. Sun, G.; Budde, R.J. Requirement for an Additional Divalent Metal Cation to Activate Protein Tyrosine Kinases. Biochemistry 1997, 36, 2139–2146. [Google Scholar] [CrossRef]
  78. Knape, M.J.; Ballez, M.; Berghardt, N.C.; Zimmermann, B.; Bertinetti, D.; Kornev, A.P.; Herberg, F.W. Divalent Metal Ions Control Activity and Inhibition of Protein Kinases. Metallomics 2017, 9, 1576–1584. [Google Scholar] [CrossRef]
  79. Jacobsen, D.M.; Bao, Z.-Q.; O’Brien, P.; Brooks, C.L., 3rd; Young, M.A. Price to be Paid for Two-Metal Catalysis. Magnesium Ions that Accelerate Chemistry Unavoidably Limit Product Release from a Protein Kinase. J. Am. Chem. Soc. 2012, 134, 15357–15370. [Google Scholar] [CrossRef]
  80. Sigel, H. Metal Ion-Assisted Stacking Interactions and the Facilitated Hydrolysis of Nucleoside 5′-Triphosphates. Pure Appl. Chem. 1998, 70, 969–976. [Google Scholar] [CrossRef]
  81. Job, M.P. Sur la Constitution des Solutions Chlorhydriques de Sels de Cobalt. C. R. Hebd. Séances Acad. Sci. 1933, 196, 181–183. [Google Scholar]
  82. Slater, J.P.; Mildvan, A.S.; Loeb, L.A. Zinc in DNA Polymerases. Biochem. Biophys. Res. Commun. 1971, 44, 37–43. [Google Scholar] [CrossRef]
  83. Palermo, G.; Cavalli, A.; Klein, M.L.; Alfonso-Prieto, M.; Dal Peraro, M.; De Vivo, M. Catalytic Metal Ions and Enzymatic Progressing of DNA and RNA. Accounts Chem. Res. 2015, 48, 220–228. [Google Scholar] [CrossRef]
  84. Wu, F.Y.-H.; Wu, C.-W. The Role of Zinc in DNA and RNA Polymerases. Met. Ions Biol. Syst. 1983, 15, 157–192. [Google Scholar]
  85. Scrutton, M.C.; Wu, C.W.; Goldthwait, D.A. The Presence and Possible Role of Zinc in RNA Polymerase Obtained from Escherichia coli. Proc. Natl. Acad. Sci. USA 1971, 68, 2497–2501. [Google Scholar] [CrossRef] [Green Version]
  86. Eichhorn, G.L. The Function of Metal Ions in Genetic Regulation. Met. Ions Biol. Syst. 1980, 10, 1–21. [Google Scholar]
  87. Genna, V.; Donati, E.; De Vivo, M. The Catalytic Mechanism of DNA and RNA Polymerases. ACS Catalysis 2018, 8, 11103–11118. [Google Scholar] [CrossRef]
  88. Broom, A.D.; Schweizer, M.P.; Ts’o, P.O.P. Interaction and Association of Bases and Nucleosides in Aqueous Solutions. V. Studies of the Association of Purine Nucleosides by Vapor Pressure Osmometry and by Proton Magnetic Resonance. J. Am. Chem. Soc. 1967, 89, 3612–3622. [Google Scholar] [CrossRef]
  89. Scheller, K.H.; Hofstetter, F.; Mitchell, P.R.; Prijs, B.; Sigel, H. Macrochelate Formation in Monomeric Metal Ion Complexes of Nucleoside 5′-Triphosphates and the Promotion of Stacking by Metal Ions. Comparison of the Self-Association of Purine and Pyrimidine 5′-Triphosphates Using Proton Nuclear Magnetic Resonance. J. Am. Chem. Soc. 1981, 103, 247–260. [Google Scholar] [CrossRef]
  90. Mitchell, P.R.; Sigel, H. A Proton Nuclear Magnetic Resonance Study of Self-Stacking in Purine and Pyrimidine Nucleosides and Nucleotides. Eur. J. Biochem. 1978, 88, 149–154. [Google Scholar] [CrossRef] [PubMed]
  91. Tribolet, R.; Sigel, H. Influence of the Protonation Degree on the Self-Association Properties of Adenosine 5′-Triphosphate (ATP). Eur. J. Biochem. 1988, 170, 617–626. [Google Scholar] [CrossRef] [PubMed]
  92. Corfù, N.A.; Tribolet, R.; Sigel, H. Comparison of the Self-Association Properties of the 5′-Triphosphates of Inosine (ITP), Guanosine (GTP), and Adenosine (ATP). Further Evidence for Ionic Interactions in the Highly Stable Dimeric [H2(ATP) ] 2 4 Stack. Eur. J. Biochem. 1990, 191, 721–735. [Google Scholar] [CrossRef]
  93. Scheller, K.H.; Sigel, H. A Proton Nuclear Magnetic Resonance Study of Purine and Pyrimidine Nucleoside 5′-Diphosphates. Extent of Macrochelate Formation in Monomeric Metal Ion Complexes and Promotion of Self-Stacking by Metal Ions. J. Am. Chem. Soc. 1983, 105, 5891–5900. [Google Scholar] [CrossRef]
  94. Sigel, H.; Blindauer, C.A.; Holý, A.; Dvořáková, H. Facilitation of the Copper(II)-Promoted Dephosphorylation of Adenosine 5′-Triphosphate (ATP4−) by the Antiviral Nucleotide Analogue, 9-[2-(Phosphonomethoxy)ethyl]adenine (PMEA). Chem. Commun. 1998, 1219–1220. [Google Scholar] [CrossRef]
  95. Sigel, H.; Scheller, K.H. On the Metal-Ion Coordinating Properties of the 5′-Monophosphates of 1,N6-Ethenoadenosine (ε-AMP), Adenosine, and Uridine. Comparison of the Macrochelate Formation in the Complexes of ε-AMP, AMP, ADP, and ATP. Eur. J. Biochem. 1984, 138, 291–299. [Google Scholar] [CrossRef]
  96. Sigel, H. Self-Association, Protonation, and Metal-Coordination of 1,N6-Ethenoadenine Derivatives in Comparison with Their Parent Compounds Adenosine, AMP and ATP. Chimia 1987, 41, 11–26. [Google Scholar] [CrossRef]
  97. Sigel, H.; Brintzinger, H. Adenosin-5′-monophosphat-N(1)-oxid. Ein ‘Ambivalentes’ Ligandsystem. Helv. Chim. Acta 1964, 47, 1701–1717. [Google Scholar] [CrossRef]
  98. Sigel, H.; Prijs, B. Ternäre Komplexe in Lösung. IV. Einfluss von 2,2’-Bipyridyl auf Stabilität und Acidität des Cu2+-Adenosin-5′-monophosphat-N(1)-oxid-1:1-Komplexes. Helv. Chim. Acta 1967, 50, 2357–2362. [Google Scholar] [CrossRef] [PubMed]
  99. Sigel, H. The Ambivalent Properties of Some Base-Modified Nucleotides. Met. Ions Biol. Syst. 1979, 8, 125–158. [Google Scholar]
  100. Sigel, H.; Naumann, C.F. Ternary Complexes in Solution. XXIV. Metal Ion Bridging of Stacked Purine-Indole Adducts. The Mixed-Ligand Complexes of Adenosine 5′-Triphosphate, Tryptophan, and Manganese(II), Copper(II), or Zinc(II). J. Am. Chem. Soc. 1976, 98, 730–739. [Google Scholar] [CrossRef]
  101. Mitchell, P.R.; Prijs, B.; Sigel, H. On the Position of the Intramolecular Equilibrium between Opened and Aromatic-Ring Stacked Forms in Ternary Complexes Composed of Adenosine 5′-Triphosphate, Mg2+ or Zn2+, and L-Tryptophanate, and in Related Ternary Systems. (Part 32). Helv. Chim. Acta 1979, 62, 1723–1735. [Google Scholar] [CrossRef]
  102. Sigel, H. Stability, Structure and Reactivity of Mixed Ligand Complexes in Solution. In Coordination Chemistry—20; Banerjea, D., Ed.; IUPAC through Pergamon Press: Oxford, UK; New York, NY, USA, 1980; pp. 27–45. [Google Scholar]
  103. Sigel, H.; Fischer, B.E.; Farkas, E. Metal Ion-Promoted Hydrophobic Interactions between Nucleotides and Amino Acids. Mixed-Ligand Adenosine 5′-Triphosphate/Metal Ion(II)/L-Leucinate Systems and Related Ternary Complexes. (Part 42). Inorg. Chem. 1983, 22, 925–934. [Google Scholar] [CrossRef]
  104. Naumann, C.F.; Sigel, H. Ternary Complexes in Solution. XVIII. The Stability-Enhancement of Nucleotide-Containing Charge-Transfer Adducts through the Formation of a Metal Ion Bridge. J. Am. Chem. Soc. 1974, 96, 2750–2756. [Google Scholar] [CrossRef]
  105. Chaudhuri, P.; Sigel, H. Ternary Complexes in Solution. 26. Stacking Interactions in the Mixed-Ligand Complexes Formed by Adenosine or Inosine 5′-Triphosphate, 2,2’-Bipyridyl, and Cobalt(II), Nickel(II), Copper(II), or Zinc(II). Evidence for Phosphate-Protonated Complexes. J. Am. Chem. Soc. 1977, 99, 3142–3150. [Google Scholar] [CrossRef]
  106. Mitchell, P.R.; Sigel, H. Enhanced Stability of Ternary Metal Ion/Adenosine 5′-Triphosphate Complexes. Cooperative Effects Caused by Stacking Interactions in Complexes Containing Adenosine Triphosphate, Phenanthroline, and Magnesium, Calcium, or Zinc Ions. (Ternary Complexes in Solution. Part 28). J. Am. Chem. Soc. 1978, 100, 1564–1570. [Google Scholar]
  107. Orioli, P.; Cini, R.; Donati, D.; Mangani, S. Crystal and Molecular Structure of the Ternary Complex Bis[(adenosine 5′-triphosphato)(2,2’-bipyridine)zinc(II)] Tetrahydrate. J. Am. Chem. Soc. 1981, 103, 4446–4452. [Google Scholar] [CrossRef]
  108. Sheldrick, W.S. Preparation and Structures of Ternary Copper(II) Complexes of ADP and ATP Models for Enzyme-Metal Ion-Nucleoside Polyphosphate Complexes. Z. Naturforsch. B Anorg. Chem. Org. Chem. 1982, 37, 863–871. [Google Scholar] [CrossRef]
  109. Buisson, D.H.; Sigel, H. Protection of Adenosine 5′-Triphosphate toward Hydrolysis by the Formation of a Mixed-Ligand Metal Ion Complex. Angew. Chem. Internat. Ed. Engl. 1972, 11, 1025. [Google Scholar] [CrossRef] [PubMed]
  110. Sigel, H.; Corfù, N.A. Facilitated Formation of High-Molecular-Weight Associates of Adenosine 5′-Triphosphate (ATP). Naturwissenschaften 1995, 82, 237–238. [Google Scholar] [CrossRef] [PubMed]
  111. Sigel, H.; Corfù, N.A. The Assisted Self-Association of ATP4− by a Poly(amino acid), [poly(Lys)], and Its Significance for Cell Organelles that Contain High Concentrations of Nucleotides. Eur. J. Biochem. 1996, 240, 508–517. [Google Scholar] [CrossRef] [PubMed]
  112. Pletscher, A.; Da Prada, M.; Berneis, K.H.; Steffen, H.; Lütold, B.; Weder, H.G. Molecular Organization of Amine Storage Organelles of Blood Platelets and Adrenal Medulla. Adv. Cytopharmacol. 1974, 2, 257–264. [Google Scholar]
  113. Da Prada, M.; Richards, J.G.; Kettler, R. Amine Storage Organelles in Platelets. In Platelets in Biology and Pathology-2; Gordon, J.L., Ed.; Elsevier: Amsterdam, The Netherlands, 1981; pp. 107–145. [Google Scholar]
  114. Berneis, K.H.; Da Prada, M.; Pletscher, A. Metal-Dependent Aggregation of Nucleotides with Formation of Biphasic Liquid Systems. Biochim. Biophys. Acta 1970, 215, 547–549. [Google Scholar] [CrossRef]
  115. Da Prada, M.; Berneis, K.H.; Pletscher, A. Storage of Catecholamines in Adrenal Medulary Granules: Formation of Aggregates with Nucleotides. Life Sci. 1971, 10, 639–646. [Google Scholar] [CrossRef]
  116. Berneis, K.H.; Goetz, U.; Da Prada, M.; Pletscher, A. Interaction of Aggregated Catecholamines and Nucleotides with Intragranular Proteins. Naunyn-Schmiedeberg’s Arch. Pharmacol. 1973, 277, 291–296. [Google Scholar] [CrossRef]
  117. Winkler, H.; Carmichael, S.W. The Chromaffin Granule. In The Secretory Granule; Poisner, A.M., Trifaró, J.M., Eds.; Elsevier Biomedical Press: Amsterdam, The Netherlands; New York, NY, USA; Oxford, UK, 1982; pp. 3–79. [Google Scholar]
  118. Winkler, H.; Apps, D.K.; Fischer-Colbrie, R. The Molecular Function of Adrenal Chromaffin Granules: Established Facts and Unresolved Topics. Neuroscience 1986, 18, 261–290. [Google Scholar] [CrossRef]
  119. Dzhandzhugazyan, K.; Bock, E. Demonstration of an Extracellular ATP-Binding Site in NCAM:  Functional Implications of Nucleotide Binding. Biochemistry 1997, 36, 15381–15395. [Google Scholar] [CrossRef]
  120. Kelley, S.O.; Holmlin, R.E.; Stemp, E.D.A.; Barton, J.K. Photoinduced Electron Transfer in Ethidium-Modified DNA Duplexes:  Dependence on Distance and Base Stacking. J. Am. Chem. Soc. 1997, 119, 9861–9870. [Google Scholar] [CrossRef]
  121. Giese, B. Electron Transfer in DNA. Curr. Opin. Chem. Biol. 2002, 6, 612–618. [Google Scholar] [CrossRef]
  122. Naumann, C.F.; Prijs, B.; Sigel, H. Adenosine and Inosine 5′-Triphosphates: Protonation, Metal Ion Coordination, and Charge-Transfer Interaction between Two Ligands within Ternary Complexes. Eur. J. Biochem. 1974, 41, 209–216. [Google Scholar] [CrossRef] [PubMed]
  123. Sigel, H.; Becker, K.; McCormick, D.B. Ternary Complexes in Solution. III. Influence of 2,2’-Bipyridyl on the Stability of 1:1 Complexes of Co2+, Ni2+, Cu2+, and Zn2+ with Hydrogen Phosphate, Adenosine 5′-Monophosphate, and Adenosine 5′-Triphosphate. Biochim. Biophys. Acta 1967, 148, 655–664. [Google Scholar] [CrossRef]
  124. Sigel, H. Hydrophobic Interactions in Biological Systems: Some Background Information Based on Ligand-Ligand Interactions in Metal Ion Complexes. Pure Appl. Chem. 1989, 61, 923–932. [Google Scholar] [CrossRef]
  125. Massoud, S.S.; Tribolet, R.; Sigel, H. Metal-Ion-Governed Molecular Recognition: Extent of Intramolecular Stack Formation in Mixed-Ligand Copper(II) Complexes Containing a Heteroaromatic N Base and an Adenosine Monophosphate (2’-AMP, 3′-AMP, or 5′-AMP). A Structuring Effect of the Metal-Ion Bridge. Eur. J. Biochem. 1990, 187, 387–393. [Google Scholar]
  126. Yamauchi, O.; Odani, A.; Masuda, H.; Sigel, H. Stacking Interactions Involving Nucleotides and Metal Ion Complexes. Met. Ions Biol. Syst. 1996, 32, 207–270. [Google Scholar] [CrossRef]
  127. Tribolet, R.; Malini-Balakrishnan, R.; Sigel, H. Influence of Decreasing Solvent-Polarity (Dioxane-Water Mixtures) on the Stability and Structure of Binary and Ternary Complexes of Adenosine 5′-Triphosphate and Uridine 5′-Triphosphate. J. Chem. Soc. Dalton Trans. 1985, 2291–2303. [Google Scholar] [CrossRef]
  128. Fischer, B.E.; Sigel, H. Ternary Complexes in Solution. 35. Intramolecular Hydrophobic Ligand-Ligand Interactions in Mixed Ligand Complexes Containing an Aliphatic Amino Acid. J. Am. Chem. Soc. 1980, 102, 2998–3008. [Google Scholar] [CrossRef]
  129. Sigel, H.; Malini-Balakrishnan, R.; Häring, U.K. Ternary Complexes in Solution (Part 48). Influence of Organic Solvents on Intramolecular Aromatic-Ring Stacks in Aqueous Mixed-Ligand Metal Ion Complexes. Opposing Solvent Effects. J. Am. Chem. Soc. 1985, 107, 5137–5148. [Google Scholar] [CrossRef]
  130. Bianchi, E.M.; Sajadi, S.A.A.; Song, B.; Sigel, H. Intramolecular Stacking Interactions in Mixed Ligand Complexes Formed by Copper(II), 2,2’-Bipyridine or 1,10-Phenanthroline, and Monoprotonated or Deprotonated Adenosine 5′-Diphosphate (ADP3−). Evaluation of Isomeric Equilibria. Inorg. Chim. Acta 2000, 300–302, 487–498. [Google Scholar] [CrossRef]
  131. Aoki, K. X-Ray Evidence for the Metal Ion Bridged Intra- and Intermolecular Stacking Interactions between Nucleotide Bases and Aromatic Heterocyclic Rings within the Ternary Complex [Cu(5′-AMP)(bpy)(H2O)]2.(NO3)2.6H2O. J. Am. Chem. Soc. 1978, 100, 7106–7108. [Google Scholar] [CrossRef]
  132. Naumann, C.F.; Sigel, H. Ternary Complexes in Solution. (Part 21). Bridging of the Stacked Adduct Between Tryptophan and Adenosine 5′-Triphosphate by Zinc(II). FEBS Lett. 1974, 47, 122–124. [Google Scholar] [CrossRef] [Green Version]
  133. Orenberg, J.B.; Fischer, B.E.; Sigel, H. Binary and Ternary Complexes of Metal Ions, Nucleoside 5′-Monophosphates, and Amino Acids. (Ternary Complexes in Solution. Part 33). J. Inorg. Nucl. Chem. 1980, 42, 785–792. [Google Scholar] [CrossRef]
  134. Sigel, H.; Scheller, K.H.; Scheller-Krattiger, V.; Prijs, B. Comparison of the Properties of Binary and Ternary Metal Ion Complexes of 1,N6-Ethenoadenosine 5′-Triphosphate (ε-ATP) and Adenosine 5′-Triphosphate (ATP), Including Macrochelate and Purine-Indole Stack Formation. J. Am. Chem. Soc. 1986, 108, 4171–4178. [Google Scholar] [CrossRef]
  135. Toulmé, J.-J. Stacking Interactions between Aromatic Amino Acids and the Adenine Ring of ATP in Zinc-Mediated Ternary Complexes. Bioinorg. Chem. 1978, 8, 319–329. [Google Scholar] [CrossRef]
  136. Basosi, R.; Gaggelli, E.; Tiezzi, E. Electron Spin Resonance Evidence of Ternary Complex Formation. J. Chem. Res. (S) 1977, 278–279. [Google Scholar]
  137. Basosi, R.; Gaggelli, E.; Tiezzi, E. Mixed Ligand Complexes in Solution of Manganous Ion. J. Chem. Soc. Faraday Trans. 1 1979, 75, 2433–2438. [Google Scholar] [CrossRef]
  138. Gaggelli, E.; Rossi, C.; Tiezzi, E. Noncovalent Interactions and Paramagnetic Relaxation Probes. Investigation of the Correlation Times. J. Phys. Chem. 1980, 84, 116–118. [Google Scholar] [CrossRef]
  139. Arena, G.; Cali, R.; Cucinotta, V.; Musumeci, S.; Rizzarelli, E.; Sammartano, S. Thermodynamics of Metal Complexes with Ligand–Ligand Interaction. Mixed Complexes of Copper(II) and Zinc(II) with Adenosine 5′-Triphosphate and L-Histidine or Histamine. J. Chem. Soc. Dalton Trans. 1984, 1651–1658. [Google Scholar] [CrossRef]
  140. Sigel, H.; Tribolet, R.; Yamauchi, O. The Imidazole Group and Its Stacking Properties in Mixed Ligand Metal Ion Complexes. Comments Inorg. Chem. 1990, 9, 305–330. [Google Scholar] [CrossRef]
  141. Chen, D.; Bastian, M.; Gregáň, F.; Holý, A.; Sigel, H. Ternary Complexes in Solution (Part 55) with Phosphonates as Ligands. Various Intramolecular Equilibria in Mixed-Ligand Complexes Containing the Antiviral 9-(2-Phosphonomethoxyethyl)adenine, an Adenosine Monophosphate Analogue. J. Chem. Soc. Dalton Trans. 1993, 1537–1546. [Google Scholar] [CrossRef]
  142. Rees, D.C. Experimental Evaluation of the Effective Dielectric Constant of Proteins. J. Mol. Biol. 1980, 141, 323–326. [Google Scholar] [CrossRef]
  143. Moore, G.R. Control of Redox Properties of Cytochrome c by Special Electrostatic Interactions. FEBS Lett. 1983, 161, 171–175. [Google Scholar] [CrossRef] [Green Version]
  144. Rogers, N.K.; Moore, G.R.; Sternberg, M.J.E. Electrostatic Interactions in Globular Proteins: Calculations of the pH Dependence of the Redox Potential of Cytochrome c551. J. Mol. Biol. 1985, 182, 613–616. [Google Scholar] [CrossRef]
  145. Sigel, H.; Martin, R.B.; Tribolet, R.; Häring, U.K.; Malini-Balakrishnan, R. An Estimation of the Equivalent Solution Dielectric Constant in the Active-Site Cavity of Metalloenzymes. Dependence of Carboxylate–Metal-Ion Complex Stabilities on the Polarity of Mixed Aqueous/Organic Solvents. Eur. J. Biochem. 1985, 152, 187–193, See also the comment in the footnote on page 258 in: Bastian, M.; Sigel, H. On the Metal Ion Binding Properties of Orotidine. Inorg. Chim. Acta 1990, 178, 249–259. [Google Scholar] [CrossRef]
  146. De Meis, L. Role of Water in the Energy of Hydrolysis of Phosphate Compounds—Energy Transduction in Biological Membranes. Biochim. Biophys. Acta 1989, 973, 333–349. [Google Scholar] [CrossRef]
  147. Åkerlöf, G.; Short, O. The Dielectric Constant of Dioxane-Water Mixtures between 0 and 80. J. Am. Chem. Soc. 1936, 58, 1241–1243, Correction in. J. Am. Chem. Soc. 1953, 75, 6357. [Google Scholar] [CrossRef]
  148. Critchfield, F.E.; Gibson, J.A., Jr.; Hall, J.L. Dielectric Constant for the Dioxane-Water System from 20 to 35°. J. Am. Chem. Soc. 1953, 75, 1991–1992. [Google Scholar] [CrossRef]
  149. Chen, D.; Gregáň, F.; Holy, A.; Sigel, H. Solvent-Dependent Metal Ion-Adenine Recognition. Quantification of the Intramolecular Equilibria between Various Isomers of the Cu2+ Complexes Formed in Water-Dioxane Mixtures with the Anions of the Antiviral 9-(2-(Phosphonomethoxy)ethyl)adenine (PMEA), an Adenosine Monophosphate (AMP) Analogue. Inorg. Chem. 1993, 32, 5377–5384. [Google Scholar]
  150. Magalhães, M.C.F.; Sigel, H. Establishment of Relations between Complex Stability and Phosphate Group Basicity for the Copper(II) Complexes of Simple Phosphate Monoesters in Dioxane-Water Mixtures. J. Indian Chem. Soc. 1992, 69, 437–441. [Google Scholar]
  151. Bastian, M.; Sigel, H. Stability and Structure of the Mg2+, Ca2+ and Cu2+ Complexes of Orotidinate 5′-Monophosphate (OMP)3− in Various Aqueous 1,4-Dioxane Mixtures. Inorg. Chim. Acta 1991, 187, 227–237. [Google Scholar] [CrossRef]
  152. Liang, G.; Sigel, H. Solvent-Dependent Metal Ion-Nucleic Base Recognition. Extent of Macrochelate Formation in the Binary Copper(II) Complexes of Adenosine 5′-Monophosphate (AMP) and Adenosine 5′-Triphosphate (ATP) in Water-Dioxane Mixtures. Inorg. Chem. 1990, 29, 3631–3632. [Google Scholar] [CrossRef]
  153. Bianchi, E.M.; Sajadi, S.A.A.; Song, B.; Sigel, H. Stabilities and Isomeric Equilibria in Aqueous Solution of Monomeric Metal Ion Complexes of Adenosine 5′-Diphosphate (ADP3−) in Comparison with Those of Adenosine 5′-Monophosphate (AMP2−). Chem. Eur. J. 2003, 9, 881–892. [Google Scholar] [CrossRef]
  154. Da Costa, C.P.; Song, B.; Gregáň, F.; Sigel, H. Intramolecular Chelate Formation Involving the Carbonyl Oxygen of Acetyl Phosphate or Acetonylphosphonate in Mixed Ligand Copper(II) Complexes Containing also 2,2’-Bipyridine or 1,10-Phenanthroline. A Decreased Solvent Polarity Favours Metal Ion-Carbonyl Oxygen Recognition. (Part 64). J. Chem. Soc. Dalton Trans. 2000, 6, 899–904. [Google Scholar]
  155. Rawn, J.D. Biochemistry; Patterson: Burlington, NC, USA, 1989. [Google Scholar]
  156. Stryer, L. Biochemistry, 4th ed.; Freeman: New York, NY, USA, 1997. [Google Scholar]
  157. Sigel, H.; Martin, R.B. Coordinating Properties of the Amide Bond. Stability and Structure of Metal Ion Complexes of Peptides and Related Ligands. Chem. Rev. 1982, 82, 385–426. [Google Scholar] [CrossRef]
  158. Furler, M.; Knobloch, B.; Sigel, R.K.O. Influence of decreased solvent permittivity on the structure and magnesium(II)-binding properties of the catalytic domain 5 of a group II intron ribozyme. Inorg. Chim. Acta 2009, 362, 771–776. [Google Scholar] [CrossRef] [Green Version]
  159. DasGupta, S. Molecular Crowding and RNA Catalysis. Org. Biomol. Chem. 2020, 18, 7724–7739. [Google Scholar] [CrossRef]
  160. Fiorini, E.; Paudel, B.; Börner, R.; Sigel, R.K.O.; Rueda, D. Optimal Molecular Crowding Accelerates Group II Intron Folding and Maximizes Catalysis. Proc. Natl. Acad. Sci. USA 2018, 115, 11917–11922. [Google Scholar]
  161. Irving, H.; Williams, R.J.P. Order of Stability of Metal Complexes. Nature 1948, 162, 746–747. [Google Scholar] [CrossRef]
  162. Irving, H.; Williams, R.J.P. The Stability of Transition-Metal Complexes. J. Chem. Soc. 1953, 3192–3210. [Google Scholar] [CrossRef]
  163. Sigel, R.K.O.; Sigel, H. Metal Ion Interactions with Nucleic Acids and Their Constituents. In Bioinorganic Fundamentals and Applications: Metals in Natural Living Systems and Metals in Toxicology and Medicine; Pecoraro, V.L., Hambley, T., Eds.; In Comprehensive Inorganic Chemistry II; Reedijk, J., Poeppelmeier, K., Eds.; Elsevier: Oxford, UK, 2013; Volume 3, pp. 623–660. [Google Scholar]
  164. Sajadi, S.A.A.; Song, B.; Gregáň, F.; Sigel, H. Acid-Base and Metal Ions-Coordinating Properties of Pyrimidine-Nucleoside 5′-Diphosphates (CDP, UDP, dTDP) and of Several Simple Diphosphate Monoesters. Establishment of Relations between Complex Stability and Diphosphate Basicity. Inorg. Chem. 1999, 38, 439–448. [Google Scholar] [CrossRef] [PubMed]
  165. Sigel, A.; Operschall, B.P.; Sigel, H. Complex Formation of Lead(II) with Nucleotides and Their Constituents. Met. Ions Life Sci. 2017, 17, 319–402. [Google Scholar]
  166. Liu, Y.; Wilson, T.J.; Lilley, D.M.J. The Structure of a Nucleolytic Ribozyme that Employs a Catalytic Metal Ion. Nature Chem. Biol. 2017, 13, 508–513. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  167. Roychowdhury-Saha, M.; Burke, D.H. Extraordinary Rates of Transition Metal Ion-Mediated Ribozyme Catalysis. RNA 2006, 12, 1846–1852. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  168. Boots, J.L.; Canny, M.D.; Azimi, E.; Pardi, A. Metal Ion Specificities for Folding and Cleavage Activity in the Schistosoma Hammerhead Ribozyme. RNA 2008, 14, 2212–2222. [Google Scholar] [CrossRef] [Green Version]
  169. Schnabl, J.; Sigel, R.K.O. Controlling Ribozyme Activity by Metal Ions. Curr. Opin. Chem. Biol. 2010, 14, 269–275. [Google Scholar] [CrossRef] [Green Version]
  170. Klawuhn, K.; Jansen, J.A.; Souchek, J.; Soukup, G.A.; Soukup, J.K. Analysis of Metal Ion-Dependence in glmS Ribozyme Self-Cleavage and Coenzyme Binding. ChemBioChem 2010, 11, 2567–2571. [Google Scholar] [CrossRef]
  171. Erat, M.C.; Sigel, R.K.O. Divalent Metal Ions Tune the Self-splicing Reaction of the Yeast Mitochondrial Group II Intron Sc.ai5γ. J. Biol. Inorg. Chem. 2008, 13, 1025–1036. [Google Scholar] [CrossRef] [Green Version]
  172. Sigel, R.K.O.; Pyle, A.M. Alternative Roles for Metal Ions in Enzyme Catalysis and the Implications for Ribozyme Chemistry. Chem. Rev. 2007, 107, 97–113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  173. Sigel, A.; Operschall, B.P.; Griesser, R.; Song, B.; Okruszek, A.; Odani, A.; Katsuta, T.; Lippert, B.; Sigel, H. (N-7)-Platination and Its Effect on (N1)H-Acidification in Nucleoside Phosphate Derivatives. Inorg. Chim. Acta 2016, 452, 137–141. [Google Scholar] [CrossRef]
  174. Sigel, H.; Sigel, A. The Bio-Relevant Metals of the Periodic Table of the Elements. Zeitschr. für Naturforschung. B 2019, 74b, 461–471. [Google Scholar] [CrossRef]
  175. Delfosse, V.; Bouchard, P.; Bonneau, E.; Dagenais, P.; Lemay, J.-F.; Lafontaine, D.A.; Legault, P. Riboswitch structure: An internal residue mimicking the purine ligand. Nucleic Acids Res. 2010, 38, 2057–2068. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  176. Schnabl, J.; Suter, P.; Sigel, R.K.O. MINAS—A Database of Metal Ions in Nucleic AcidS. Nucleic Acids Res. 2012, 40, D434–D438. [Google Scholar] [CrossRef] [Green Version]
  177. Kapinos, L.E.; Operschall, B.P.; Larsen, E.; Sigel, H. Understanding the Acid-Base Properties of Adenosine: The Intrinsic Basicities of N1, N3 and N7. Chem. Eur. J. 2011, 17, 8156–8164. [Google Scholar] [CrossRef]
  178. Saenger, W. Principles of Nucleic Acid Structure; Springer: New York, NY, USA, 1984. [Google Scholar]
  179. Erat, M.C.; Kovacs, H.; Sigel, R.K.O. Metal Ion-N7 Coordination in a Ribozyme Branch Domain by NMR. J. Inorg. Biochem. 2010, 104, 611–613. [Google Scholar] [CrossRef] [Green Version]
  180. De Clercq, E.; Holý, A.; Rosenberg, I.; Sakuma, T.; Balzarini, J.; Maudgal, P.C. A Novel Selective Broad-Spectrum Anti-DNA Virus Agent. Nature 1986, 323, 464–467. [Google Scholar] [CrossRef]
  181. Blindauer, C.A.; Holý, A.; Sigel, A.; Operschall, B.P.; Griesser, R.; Sigel, H. Acid–Base Properties of an Antivirally Active Acyclic Nucleoside Phosphonate: (S)-9-[3-Hydroxy-2-(phosphonomethoxy)propyl]adenine (HPMPA). New J. Chem. 2022, 46, 6484–6493. [Google Scholar] [CrossRef]
  182. Blindauer, C.A.; Griesser, R.; Holý, A.; Operschall, B.P.; Sigel, A.; Song, B.H.; Sigel, H. Intramolecular π-Stacks in Mixed-Ligand Copper(II) Complexes Formed by Heteroaromatic Amines and Antivirally Active Acyclic Nucleotide Analogues Carrying a Hydroxy-2-(phosphonomethoxy)propyl Residue. J. Coord. Chem. 2018, 71, 1910–1934. [Google Scholar] [CrossRef]
  183. Blindauer, C.A.; Holý, A.; Operschall, B.P.; Sigel, A.; Song, B.; Sigel, H. Metal Ion-Coordinating Properties in Aqueous Solution of the Antivirally Active Nucleotide Analogue (S)-9-[3-Hydroxy-2-(phosphonomethoxy)propyl]adenine (HPMPA). Quantification of Complex Isomeric Equilibria. Eur. J. Inorg. Chem. 2019, 2019, 3892–3903. [Google Scholar] [CrossRef]
  184. Blindauer, C.A.; Sigel, A.; Operschall, B.P.; Holý, A.; Sigel, H. Metal-Ion Binding Properties of (S)-1-[3-Hydroxy-2-(phosphonomethoxy)propyl]cytosine (HPMPC, Cidofovir). A Nucleotide Analogue with Activity Against DNA Viruses. Inorg. Chim. Acta 2018, 472, 283–294. [Google Scholar] [CrossRef]
  185. De Clercq, E. Towards an Effective Chemotherapy of Virus Infections: Therapeutic Potential of Cidofovir [(S)-1-[3-Hydroxy-2-(phosphonomethoxy)propyl]cytosine, HPMPC] for the Treatment of DNA Virus Infections. Collect. Czech. Chem. Commun. 1998, 63, 480–506. [Google Scholar] [CrossRef]
  186. De Clercq, E.; Li, G. Approved Antiviral Drugs over the Past 50 Years. Clin. Microbiol. Rev. 2016, 29, 695–747. [Google Scholar] [CrossRef] [Green Version]
  187. Edwards, A.M.; Baric, R.S.; Ollmann Saphire, E.; Ullmer, J.B. Stopping Pandemics before They Start: Lessons Learned from SARS-CoV-2. Science 2022, 375, 1133–1139. [Google Scholar] [CrossRef]
  188. Klein, D.J.; Moore, P.B.; Steitz, T.A. The Contribution of Metal Ions to the Structural Stability of the Large Ribosomal Subunit. RNA 2004, 10, 1366–1379. [Google Scholar] [CrossRef] [Green Version]
  189. Erat, M.C.; Sigel, R.K.O. Determination of the Intrinsic Affinities of Multiple Site-Specific Mg2+ Ions Coordinated to Domain 6 of a Group II Intron Ribozyme. Inorg. Chem. 2007, 46, 11224–11234. [Google Scholar] [CrossRef] [Green Version]
  190. Erat, M.C.; Zerbe, O.; Fox, T.; Sigel, R.K.O. Solution Structure of Domain 6 from a Self-Splicing Group II Intron Ribozyme: A Mg2+ Binding Site is Located Close to the Stacked Branch Adenosine. ChemBioChem 2007, 8, 306–314. [Google Scholar] [CrossRef] [Green Version]
  191. Sigel, R.K.O.; Vaidya, A.; Pyle, A.M. Metal Ion Binding Sites in a Group II Intron Core. Nat. Struct. Biol. 2000, 7, 1111–1116. [Google Scholar]
  192. Korth, M.M.T.; Sigel, R.K.O. Unusually High-Affinity Mg2+ Binding at the AU-Rich Sequence within the Antiterminator Hairpin of a Mg2+ Riboswitch. Chem. Biodivers. 2012, 9, 2035–2049. [Google Scholar] [CrossRef] [Green Version]
  193. Yin, Y.W.; Steitz, T.A. The Structural Mechanism of Translocation and Helicase Activity in T7 RNA Polymerase. Cell 2004, 116, 393–404. [Google Scholar] [CrossRef] [Green Version]
  194. Shan, S.-O.; Yoshida, A.; Sun, S.; Piccirilli, J.A.; Herschlag, D. Three Metal Ions at the Active Site of the Tetrahymena Group I Ribozyme. Proc. Natl. Acad. Sci. USA 1999, 96, 12299–12304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  195. Shan, S.-O.; Kravchuk, A.V.; Piccirilli, J.A.; Herschlag, D. Defining the Catalytic Metal Ion Interactions in the Tetrahymena Ribozyme Reaction. Biochemistry 2001, 40, 5161–5171. [Google Scholar] [CrossRef] [PubMed]
  196. Manigrasso, J.; Chillón, I.; Genna, V.; Vidossich, P.; Somarowthu, S.; Pyle, A.M.; De Vivo, M.; Marcia, M. Visualizing Group II Intron Dynamics between the First and Second Step of Splicing. Nature Commun. 2020, 11, 2837. [Google Scholar] [CrossRef]
  197. Takeuchi, H.; Murata, H.; Harada, I. Interaction of Adenosine 5′-Triphosphate with Mg2+: Vibrational Study of Coordination Sites by Use of 18O-Labeled Triphosphates. J. Am. Chem. Soc. 1988, 110, 392–397. [Google Scholar] [CrossRef]
Figure 3. Comparison of the chemical structure of the dianion of 9-[2-(phosphonomethoxy)ethyl]adenine (PMEA2−) with the one of adenosine 5′-monophosphate (AMP2−). The orientation of PMEA2− in solution [36] and in the solid state [37] resembles the anti conformation of AMP2− [4,5,6,38] and thus PMEA2− may be considered as an analogue of AMP2− or 2′-deoxy-AMP2−.
Figure 3. Comparison of the chemical structure of the dianion of 9-[2-(phosphonomethoxy)ethyl]adenine (PMEA2−) with the one of adenosine 5′-monophosphate (AMP2−). The orientation of PMEA2− in solution [36] and in the solid state [37] resembles the anti conformation of AMP2− [4,5,6,38] and thus PMEA2− may be considered as an analogue of AMP2− or 2′-deoxy-AMP2−.
Molecules 27 02625 g003
Figure 4. Evidence for enhanced stabilities of some M(PMEA) (●) complexes in comparison with those of the corresponding M(PME) (⊗) species, based on the relationship between log K M ( R - PO 3 ) M and p K H ( R - PO 3 ) H   for M(R-PO3) complexes of some simple phosphate monoester and phosphonate ligands ( R - PO 3 2 ) (◯): 4-nitrophenyl phosphate (NPhP2−), phenyl phosphate (PhP2−), uridine 5′-monophosphate (UMP2−), D-ribose 5-monophosphate (RibMP2−), thymidine [= 1-(2′-deoxy-β-D-ribofuranosyl)thymine] 5′-monophosphate (TMP2−), n-butyl phosphate (BuP2−), methanephosphonate (MeP2−), and ethanephosphonate (EtP2−) (from left to right). The least squares lines (Equation (3)) are drawn through the corresponding 8 data sets (◯) taken from Ref [48] for the phosphate monoesters and from Ref [44] for the phosphonates. The data due to the M2+/PMEA (●) and the M2+/PME (⊗) systems are from Ref [44]. The vertical broken lines emphasize the stability differences to the reference lines; they equal log Δ as defined in Equation (12) for the M(PME) complexes (for the M(PMEA) complexes the analogous formulation holds). All the plotted equilibrium constants refer to aqueous solutions at 25 °C and I = 0.1 M (NaNO3). Reproduced with permission from our publication in Coordination Chemistry Reviews [45]; copyright 1995, Elsevier Science S.A., Lausanne, Switzerland.
Figure 4. Evidence for enhanced stabilities of some M(PMEA) (●) complexes in comparison with those of the corresponding M(PME) (⊗) species, based on the relationship between log K M ( R - PO 3 ) M and p K H ( R - PO 3 ) H   for M(R-PO3) complexes of some simple phosphate monoester and phosphonate ligands ( R - PO 3 2 ) (◯): 4-nitrophenyl phosphate (NPhP2−), phenyl phosphate (PhP2−), uridine 5′-monophosphate (UMP2−), D-ribose 5-monophosphate (RibMP2−), thymidine [= 1-(2′-deoxy-β-D-ribofuranosyl)thymine] 5′-monophosphate (TMP2−), n-butyl phosphate (BuP2−), methanephosphonate (MeP2−), and ethanephosphonate (EtP2−) (from left to right). The least squares lines (Equation (3)) are drawn through the corresponding 8 data sets (◯) taken from Ref [48] for the phosphate monoesters and from Ref [44] for the phosphonates. The data due to the M2+/PMEA (●) and the M2+/PME (⊗) systems are from Ref [44]. The vertical broken lines emphasize the stability differences to the reference lines; they equal log Δ as defined in Equation (12) for the M(PME) complexes (for the M(PMEA) complexes the analogous formulation holds). All the plotted equilibrium constants refer to aqueous solutions at 25 °C and I = 0.1 M (NaNO3). Reproduced with permission from our publication in Coordination Chemistry Reviews [45]; copyright 1995, Elsevier Science S.A., Lausanne, Switzerland.
Molecules 27 02625 g004
Figure 5. Chemical structures of the dianions of 9-[2-(phosphonomethoxy)ethyl]adenine (PMEA2−), 9-[2-(phosphonomethoxy)ethyl]-2-aminopurine (PME2AP2−), 9-[2-(phosphonomethoxy)ethyl]-2,6-diaminopurine (PMEDAP2−), and 9-[2-(phosphonomethoxy)ethyl]-2-amino-6-dimethylaminopurine (PME2A6DMAP2−). These compounds are abbreviated as PE2−. For the solution structure of PMEA2−, see legend of Figure 3. Replacement of the adenine residue in PMEA2− by a cytosine residue (Figure 1) gives PMEC2− = 1-[2-(phosphonomethoxy)ethyl]cytosine.
Figure 5. Chemical structures of the dianions of 9-[2-(phosphonomethoxy)ethyl]adenine (PMEA2−), 9-[2-(phosphonomethoxy)ethyl]-2-aminopurine (PME2AP2−), 9-[2-(phosphonomethoxy)ethyl]-2,6-diaminopurine (PMEDAP2−), and 9-[2-(phosphonomethoxy)ethyl]-2-amino-6-dimethylaminopurine (PME2A6DMAP2−). These compounds are abbreviated as PE2−. For the solution structure of PMEA2−, see legend of Figure 3. Replacement of the adenine residue in PMEA2− by a cytosine residue (Figure 1) gives PMEC2− = 1-[2-(phosphonomethoxy)ethyl]cytosine.
Molecules 27 02625 g005
Figure 6. Equilibrium schemes involving various N sites in metal ion complexes of acyclic nucleoside phosphonates containing a purine moiety.
Figure 6. Equilibrium schemes involving various N sites in metal ion complexes of acyclic nucleoside phosphonates containing a purine moiety.
Molecules 27 02625 g006
Figure 7. Evidence for enhanced stabilities of some M(PME2A6DMAP) (●) complexes in comparison with those of the corresponding M(PMEA) (⊗, M(PMEDAP) (▽), and M(PME-R) (◇) species, based on the relationship between log K M ( R - PO 3 ) M and p K H ( R - PO 3 ) H   for M(R-PO3) complexes of some simple phosphate monoester and phosphonate ligands ( R - PO 3 2 ) (◯) (see legend to Figure 4). The points due to the equilibrium constants for the M2+/PME2A6DMAP systems (●) are based on the values listed in Tables 1 and 2 of Ref [61]; those for the M2+/PMEA systems (⊗) are from Ref [44], for the M2+/PMEDAP systems (▽) from Ref [62], and those for the M2+/PME-R systems (◇) are based on p K H ( PME - R ) H = 6.99 ± 0.04 (average of the pKa values for H(PME) (7.02) and H(PMEC) (6.95) [63]) and the stability enhancements listed in Ref [63]. The vertical broken lines emphasize the stability differences to the reference lines; they equal log ΔM/PE, as defined in Equation (12) for the M(PE) complexes in general. All the plotted equilibrium constants refer to aqueous solutions at 25 °C and I = 0.1 M (NaNO3). Reproduced with permission from our publication in the Canadian Journal of Chemistry [61]; copyright 2014, NRC Research Press.
Figure 7. Evidence for enhanced stabilities of some M(PME2A6DMAP) (●) complexes in comparison with those of the corresponding M(PMEA) (⊗, M(PMEDAP) (▽), and M(PME-R) (◇) species, based on the relationship between log K M ( R - PO 3 ) M and p K H ( R - PO 3 ) H   for M(R-PO3) complexes of some simple phosphate monoester and phosphonate ligands ( R - PO 3 2 ) (◯) (see legend to Figure 4). The points due to the equilibrium constants for the M2+/PME2A6DMAP systems (●) are based on the values listed in Tables 1 and 2 of Ref [61]; those for the M2+/PMEA systems (⊗) are from Ref [44], for the M2+/PMEDAP systems (▽) from Ref [62], and those for the M2+/PME-R systems (◇) are based on p K H ( PME - R ) H = 6.99 ± 0.04 (average of the pKa values for H(PME) (7.02) and H(PMEC) (6.95) [63]) and the stability enhancements listed in Ref [63]. The vertical broken lines emphasize the stability differences to the reference lines; they equal log ΔM/PE, as defined in Equation (12) for the M(PE) complexes in general. All the plotted equilibrium constants refer to aqueous solutions at 25 °C and I = 0.1 M (NaNO3). Reproduced with permission from our publication in the Canadian Journal of Chemistry [61]; copyright 2014, NRC Research Press.
Molecules 27 02625 g007
Figure 8. Dependence of the initial rate, v0 = d[PO4]/dt (M s−1), of the Cu2+-promoted dephosphorylation of ATP ([Cu2+]tot = [ATP]tot = 10−3 M) in aqueous solution on the addition of further divalent metal ions (empty circles, half-filled circles, full circles) or the Cu2+ 1:1 complex with diethylenetriamine (= dien = 1,4,7-triazaheptane (▽)) at pH0 5.50; I = 0.1 M, NaClO4; 50 °C. The broken line portions indicate uncertainty due to precipitation. The figure is reproduced with permission from our publication in Inorganica Chimica Acta [55]; copyright 1992, Elsevier Sequoia (see also [14]).
Figure 8. Dependence of the initial rate, v0 = d[PO4]/dt (M s−1), of the Cu2+-promoted dephosphorylation of ATP ([Cu2+]tot = [ATP]tot = 10−3 M) in aqueous solution on the addition of further divalent metal ions (empty circles, half-filled circles, full circles) or the Cu2+ 1:1 complex with diethylenetriamine (= dien = 1,4,7-triazaheptane (▽)) at pH0 5.50; I = 0.1 M, NaClO4; 50 °C. The broken line portions indicate uncertainty due to precipitation. The figure is reproduced with permission from our publication in Inorganica Chimica Acta [55]; copyright 1992, Elsevier Sequoia (see also [14]).
Molecules 27 02625 g008
Figure 9. Variation in the proportions of ATP present in the monomer (1), dimer (2), trimer (3), …, and octamer (8) in D2O solutions in dependence on the total concentration of ATP4− (K = 1.3 M−1), Mg(ATP)2− (K = 4.0 M−1), and Zn(ATP)2− ( K D * = 20 M−1 and Kst = 4 M−1) (see text in Section 6.2 and Equations (18)–(20)) (27 °C, I = 0.1 to ~2 M (NaNO3) in D2O). This figure is reproduced with permission from our publication in the Journal of the American Chemical Society [89]; copyright 1981, American Chemical Society (see also [53]).
Figure 9. Variation in the proportions of ATP present in the monomer (1), dimer (2), trimer (3), …, and octamer (8) in D2O solutions in dependence on the total concentration of ATP4− (K = 1.3 M−1), Mg(ATP)2− (K = 4.0 M−1), and Zn(ATP)2− ( K D * = 20 M−1 and Kst = 4 M−1) (see text in Section 6.2 and Equations (18)–(20)) (27 °C, I = 0.1 to ~2 M (NaNO3) in D2O). This figure is reproduced with permission from our publication in the Journal of the American Chemical Society [89]; copyright 1981, American Chemical Society (see also [53]).
Molecules 27 02625 g009
Figure 10. TOP: Proposed structure of the reactive [M2(ATP)]2(OH) dimer, which occurs in low concentration during the metal ion-promoted dephosphorylation of ATP. The intramolecular attack of OH is indicated on the right-hand side, while the left-hand side is ready to transfer also into the reactive state by deprotonation of the coordinated water molecule or to undergo an intramolecular water attack, corresponding to the dimeric [M2(ATP)]2 species. BOTTOM: Probable structure of the reactive Cu3(ATP)(AMP)(OH) species. The intramolecular attack of OH is indicated on the right-hand side, while the left-hand side shows the metal ion bridge stabilizing the purine stack by M2+ coordination to the phosphate group of AMP2− and to N7 of ATP4− [14,80,94]. The structures are adapted from Ref [80].
Figure 10. TOP: Proposed structure of the reactive [M2(ATP)]2(OH) dimer, which occurs in low concentration during the metal ion-promoted dephosphorylation of ATP. The intramolecular attack of OH is indicated on the right-hand side, while the left-hand side is ready to transfer also into the reactive state by deprotonation of the coordinated water molecule or to undergo an intramolecular water attack, corresponding to the dimeric [M2(ATP)]2 species. BOTTOM: Probable structure of the reactive Cu3(ATP)(AMP)(OH) species. The intramolecular attack of OH is indicated on the right-hand side, while the left-hand side shows the metal ion bridge stabilizing the purine stack by M2+ coordination to the phosphate group of AMP2− and to N7 of ATP4− [14,80,94]. The structures are adapted from Ref [80].
Molecules 27 02625 g010
Figure 11. Influence of AMP (◯) and PMEA (●) as well as of ε-AMP (Figure 12) on the initial rate v0 (M s−1) of the dephosphorylation of the Cu2+/ATP 2:1 system ([Cu2+]tot = 2 × 10−3 M and [ATP]tot = 10−3 M) in aqueous solution at pH 6.70 (I = 0.1 M, NaClO4; 50 °C). This figure is adapted with permission from Figure 3 of our publication in Chem. Commun. [94]; copyright 1998, The Royal Society of Chemistry.
Figure 11. Influence of AMP (◯) and PMEA (●) as well as of ε-AMP (Figure 12) on the initial rate v0 (M s−1) of the dephosphorylation of the Cu2+/ATP 2:1 system ([Cu2+]tot = 2 × 10−3 M and [ATP]tot = 10−3 M) in aqueous solution at pH 6.70 (I = 0.1 M, NaClO4; 50 °C). This figure is adapted with permission from Figure 3 of our publication in Chem. Commun. [94]; copyright 1998, The Royal Society of Chemistry.
Molecules 27 02625 g011
Figure 12. Chemical structures of AMP2− analogues: tubercidin 5′-monophosphate (TuMP2− = 7-deaza-AMP2−), 1,N6-ethenoadenosine 5′-monophosphate (ε-AMP2−), and adenosine 5′-monophosphate N(1)-oxide (AMP·NO2−).
Figure 12. Chemical structures of AMP2− analogues: tubercidin 5′-monophosphate (TuMP2− = 7-deaza-AMP2−), 1,N6-ethenoadenosine 5′-monophosphate (ε-AMP2−), and adenosine 5′-monophosphate N(1)-oxide (AMP·NO2−).
Molecules 27 02625 g012
Figure 13. Simplified picture of a matrix-assisted ATP stack indicating an electron and thus an information transfer through this stack from one end to the other. This figure is reproduced with permission from our publication in Pure and Applied Chemistry [21]; copyright 2004, International Union of Pure and Applied Chemistry (IUPAC).
Figure 13. Simplified picture of a matrix-assisted ATP stack indicating an electron and thus an information transfer through this stack from one end to the other. This figure is reproduced with permission from our publication in Pure and Applied Chemistry [21]; copyright 2004, International Union of Pure and Applied Chemistry (IUPAC).
Molecules 27 02625 g013
Figure 14. Probable schematic structure in solution of an intramolecular stack in the Cu(bpy)(ATP)2− complex [105]; adapted from similar structures for Cu(bpy)(AMP), as shown in Refs [124,125,126].
Figure 14. Probable schematic structure in solution of an intramolecular stack in the Cu(bpy)(ATP)2− complex [105]; adapted from similar structures for Cu(bpy)(AMP), as shown in Refs [124,125,126].
Molecules 27 02625 g014
Figure 15. Tentative and schematic structure of the stacked isomer in solution of an M(ATP)(Trp)3− complex. Adapted from Ref [101].
Figure 15. Tentative and schematic structure of the stacked isomer in solution of an M(ATP)(Trp)3− complex. Adapted from Ref [101].
Molecules 27 02625 g015
Figure 16. Evidence for an enhanced stability of the Cu(arm)(PME) and Cu(arm)(PMEA) complexes based on the relationship between log K Cu ( arm ) ( R - PO 3 ) Cu ( arm ) and p K H ( R - PO 3 ) H for the ternary. Cu(bpy)(PME) and Cu(bpy)(PMEA) (◯), as well as for the Cu(phen)(PME) and Cu(phen)(PMEA) (●) complexes in aqueous solution at I = 0.1 M (NaNO3) and 25 °C. The plotted data are from Table 6 in Section 8 of [141] (= Table 7 here). The two reference lines represent the log K versus pKa relationship for Cu(arm)(R-PO3) complexes; it should be emphasized that R-P O 3 2 symbolizes here phosphonates (or phosphate monoesters) with an R group unable to undergo any kind of hydrophobic, stacking, or other type of interaction; the broken line holds for arm = bpy and the solid line for arm = phen. Both straight lines are calculated with the parameters listed in Table 5 of Ref [141], and they represent the situation for ternary complexes without an intramolecular ligand–ligand interaction. Redrawn in a slightly altered way from Ref [141] with permission of the Royal Society of Chemistry and then reproduced with permission from our publication in Coordination Chemistry Reviews [45]; copyright 1995, Elsevier Science S.A., Lausanne, Switzerland.
Figure 16. Evidence for an enhanced stability of the Cu(arm)(PME) and Cu(arm)(PMEA) complexes based on the relationship between log K Cu ( arm ) ( R - PO 3 ) Cu ( arm ) and p K H ( R - PO 3 ) H for the ternary. Cu(bpy)(PME) and Cu(bpy)(PMEA) (◯), as well as for the Cu(phen)(PME) and Cu(phen)(PMEA) (●) complexes in aqueous solution at I = 0.1 M (NaNO3) and 25 °C. The plotted data are from Table 6 in Section 8 of [141] (= Table 7 here). The two reference lines represent the log K versus pKa relationship for Cu(arm)(R-PO3) complexes; it should be emphasized that R-P O 3 2 symbolizes here phosphonates (or phosphate monoesters) with an R group unable to undergo any kind of hydrophobic, stacking, or other type of interaction; the broken line holds for arm = bpy and the solid line for arm = phen. Both straight lines are calculated with the parameters listed in Table 5 of Ref [141], and they represent the situation for ternary complexes without an intramolecular ligand–ligand interaction. Redrawn in a slightly altered way from Ref [141] with permission of the Royal Society of Chemistry and then reproduced with permission from our publication in Coordination Chemistry Reviews [45]; copyright 1995, Elsevier Science S.A., Lausanne, Switzerland.
Molecules 27 02625 g016
Figure 17. Tentative and simplified structure of a species with an intramolecular stack for Cu(phen)(PMEA) in solution [45,141]. Adapted from Ref [45].
Figure 17. Tentative and simplified structure of a species with an intramolecular stack for Cu(phen)(PMEA) in solution [45,141]. Adapted from Ref [45].
Molecules 27 02625 g017
Figure 18. Evidence for an enhanced stability of the Cu(PME) (●) complex in 1,4-dioxane-water mixtures as solvents, based on the relationship between log K Cu ( R - PO 3 ) Cu and p K H ( R - PO 3 ) H for the Cu2+ 1:1 complexes of 4-nitrophenyl phosphate (1), phenyl phosphate (2), D-ribose 5-monophosphate (3), n-butyl phosphate (4),uridine 5′-monophosphate (5), thymidine [= 1-(2′-deoxy-β-D-ribofuranosyl)thymine] 5′-monophosphate (6), methanephosphonate (7), and ethanephosphonate (8) in water and in water containing 30% or 50% (v/v) 1,4-dioxane. The least-squares lines are drawn in each case through the data sets shown [44,150]; the equations for these reference lines are given in Ref [149]. The data points due to the methanephosphonate system in the mixed solvents (⊗,7) (see Ref [149]) are shown to prove that simple phosphonates fit within the experimental error limits on the reference lines established with phosphate monoester systems (see also Ref [149). The points due to the Cu2+ 1:1 complexes formed with PME2− (●) in the three mentioned solvents demonstrate the enhanced complex stabilities. The vertical broken lines emphasize the stability differences to the reference lines; these differences equal log ΔCu(R-PO3), as defined in Equation (12). All the plotted equilibrium constants refer to 25 °C and I = 0.1 M (NaNO3). This figure is adapted with permission from our publication in Inorganic Chemistry [149]; copyright 1993, American Chemical Society (see also [45]).
Figure 18. Evidence for an enhanced stability of the Cu(PME) (●) complex in 1,4-dioxane-water mixtures as solvents, based on the relationship between log K Cu ( R - PO 3 ) Cu and p K H ( R - PO 3 ) H for the Cu2+ 1:1 complexes of 4-nitrophenyl phosphate (1), phenyl phosphate (2), D-ribose 5-monophosphate (3), n-butyl phosphate (4),uridine 5′-monophosphate (5), thymidine [= 1-(2′-deoxy-β-D-ribofuranosyl)thymine] 5′-monophosphate (6), methanephosphonate (7), and ethanephosphonate (8) in water and in water containing 30% or 50% (v/v) 1,4-dioxane. The least-squares lines are drawn in each case through the data sets shown [44,150]; the equations for these reference lines are given in Ref [149]. The data points due to the methanephosphonate system in the mixed solvents (⊗,7) (see Ref [149]) are shown to prove that simple phosphonates fit within the experimental error limits on the reference lines established with phosphate monoester systems (see also Ref [149). The points due to the Cu2+ 1:1 complexes formed with PME2− (●) in the three mentioned solvents demonstrate the enhanced complex stabilities. The vertical broken lines emphasize the stability differences to the reference lines; these differences equal log ΔCu(R-PO3), as defined in Equation (12). All the plotted equilibrium constants refer to 25 °C and I = 0.1 M (NaNO3). This figure is adapted with permission from our publication in Inorganic Chemistry [149]; copyright 1993, American Chemical Society (see also [45]).
Molecules 27 02625 g018
Figure 19. Formation degrees of the chelated isomers, CuLcl, in the Cu(AMP) and Cu(PMEA) systems (Table 11) as well as for comparison in the Cu(DHAP) [47], Cu(G1P) [47], and Cu(AnP) [41,47,154] complex systems (Figure 20) as a function of the percentage of 1,4-dioxane added to an aqueous reagent mixture. All the plotted equilibrium constants refer to 25 °C and I = 0.1 M (NaNO3). This figure is reproduced with permission from our publication in Coordination Chemistry Reviews [47]; copyright 2000, Elsevier, Science S.A., Lausanne, Switzerland).
Figure 19. Formation degrees of the chelated isomers, CuLcl, in the Cu(AMP) and Cu(PMEA) systems (Table 11) as well as for comparison in the Cu(DHAP) [47], Cu(G1P) [47], and Cu(AnP) [41,47,154] complex systems (Figure 20) as a function of the percentage of 1,4-dioxane added to an aqueous reagent mixture. All the plotted equilibrium constants refer to 25 °C and I = 0.1 M (NaNO3). This figure is reproduced with permission from our publication in Coordination Chemistry Reviews [47]; copyright 2000, Elsevier, Science S.A., Lausanne, Switzerland).
Molecules 27 02625 g019
Figure 20. Chemical structures of dihydroxyacetone phosphate (DHAP2−) and glycerol 1-phosphate (G1P)2−. Acetonylphosphonate (AnP2−) is often employed as a model compound for acetyl phosphate (AcP2−), which is biologically very relevant but difficult to study because it is very hydrolysis-sensitive. It may be added that in the so-called [155] α-glycerophosphate shuttle, DHAP2− and G1P2− are interconverted into each other [155,156].
Figure 20. Chemical structures of dihydroxyacetone phosphate (DHAP2−) and glycerol 1-phosphate (G1P)2−. Acetonylphosphonate (AnP2−) is often employed as a model compound for acetyl phosphate (AcP2−), which is biologically very relevant but difficult to study because it is very hydrolysis-sensitive. It may be added that in the so-called [155] α-glycerophosphate shuttle, DHAP2− and G1P2− are interconverted into each other [155,156].
Molecules 27 02625 g020
Figure 21. Irving–Williams sequence-type plots [158] for the stability constants log K M ( N ) M (N = nucleotide) of the 1:1 complexes of Ba2+ through Cd2+ formed with mono- (R-MP2−), di- (R-DP3−), and triphosphate monoesters (R-TP4−) (◯), as well as of those for AMP2−, ADP3−, and ATP4− (●). The plotted data of the phosphate ligands are from Table 13 of [23]; they also represent the stability constants for the M2+ complexes of the pyrimidine-nucleoside 5′-mono-, di- or triphosphates (except for Cu(CTP)2− [74]). The log stability constants for the M2+ complexes of AMP2− and ADP3− are from Table 4 of [153]; those for ATP4− are from Table II of [74] (25 °C; I = 0.1 M, NaNO3) ([23,74,163,164] as well as pp. 364–365 in [165]).
Figure 21. Irving–Williams sequence-type plots [158] for the stability constants log K M ( N ) M (N = nucleotide) of the 1:1 complexes of Ba2+ through Cd2+ formed with mono- (R-MP2−), di- (R-DP3−), and triphosphate monoesters (R-TP4−) (◯), as well as of those for AMP2−, ADP3−, and ATP4− (●). The plotted data of the phosphate ligands are from Table 13 of [23]; they also represent the stability constants for the M2+ complexes of the pyrimidine-nucleoside 5′-mono-, di- or triphosphates (except for Cu(CTP)2− [74]). The log stability constants for the M2+ complexes of AMP2− and ADP3− are from Table 4 of [153]; those for ATP4− are from Table II of [74] (25 °C; I = 0.1 M, NaNO3) ([23,74,163,164] as well as pp. 364–365 in [165]).
Molecules 27 02625 g021
Table 1. Negative logarithm of the acidity constants (Equations (1) and (2)) for methyl phosphoric acid and methylphosphonic acid [H2(P)] in aqueous solution at 25 °C and I = 0.1 M (NaNO3) a.
Table 1. Negative logarithm of the acidity constants (Equations (1) and (2)) for methyl phosphoric acid and methylphosphonic acid [H2(P)] in aqueous solution at 25 °C and I = 0.1 M (NaNO3) a.
H2(P) p K H 2 ( P ) H p K H ( P ) H Ref b
CH3OPO(OH)21.1 ± 0.26.36 ± 0.01[40]
CH3PO(OH)22.10 ± 0.037.51 ± 0.01[41]
a So-called practical, mixed or Brønsted constants [39] are listed—i.e., for the H+ ion not the concentration but the activity (as measured by a glass electrode) is used in the calculation of an acidity constant (Equations (1) and (2)). The practical constant may be converted into the concentration constant by subtracting 0.02 from the listed pKa values (25 °C; I = 0.1 M, NaNO3).—The error limits correspond to three times the standard error of the mean value or the sum of the probable systematic errors, whichever is larger. b See also Refs [42,43].
Table 2. Stability enhancements (Equation (12)) for the M(PMEA) and M(PME) complexes and their comparison according to Equation (13). In addition, the extent of chelate formation according to Equilibrium (5) is given as expressed by the dimension-less equilibrium constant KI (Equations (6) and (10)) and the percentages of M(PME)cl (Equation (14)). (Aqueous solution; 25 °C; I = 0.1 M, NaNO3) a,b.
Table 2. Stability enhancements (Equation (12)) for the M(PMEA) and M(PME) complexes and their comparison according to Equation (13). In addition, the extent of chelate formation according to Equilibrium (5) is given as expressed by the dimension-less equilibrium constant KI (Equations (6) and (10)) and the percentages of M(PME)cl (Equation (14)). (Aqueous solution; 25 °C; I = 0.1 M, NaNO3) a,b.
M2+log ΔM(PMEA)log ΔM(PME)Δ log ΔKI/PME% M(PME)cl
Mg2+0.16 ± 0.050.22 ± 0.03–0.06 ± 0.060.66 ± 0.1240 ± 4
Ca2+0.11 ± 0.070.14 ± 0.05–0.03 ± 0.090.38 ± 0.1628 ± 9
Sr2+0.07 ± 0.050.07 ± 0.050.00 ± 0.070.17 ± 0.1415 ± 10
Ba2+0.08 ± 0.060.10 ± 0.05–0.02 ± 0.080.26 ± 0.1421 ± 9
Mn2+0.21 ± 0.080.27 ± 0.05–0.06 ± 0.090.86 ± 0.2346 ± 7
Co2+0.28 ± 0.070.29 ± 0.06–0.01 ± 0.090.95 ± 0.2849 ± 7
Ni2+0.30 ± 0.070.19 ± 0.050.11 ± 0.090.55 ± 0.1935 ± 8
Cu2+0.77 ± 0.070.48 ± 0.070.29 ± 0.102.02 ± 0.4767 ± 5
Zn2+0.30 ± 0.100.34 ± 0.06–0.04 ± 0.121.19 ± 0.3254 ± 7
Cd2+0.33 ± 0.050.30 ± 0.050.03 ± 0.081.00 ± 0.2550 ± 6
a The values are compiled from Tables 2 and 3 of [45]. For the error limits, see the terminal sentence of Footnote ‘a’ in Table 1; derived data were calculated according to the error propagation after Gauss. b Acidity constants: p K H 2 ( PMEA ) H = 4.16 ± 0.02 (deprotonation of the (N1)H+ site), p K H ( PMEA ) H = 6.90 ± 0.01 and p K H ( PME ) H = 7.02 ± 0.01 (deprotonation of the P(O)2(OH) residue) (from Table 2 in [45]).
Table 3. Formation degrees of the isomeric complexes, M(PME2AP)op, M(PME2AP)cl/O, and M(PME2AP)cl/N7 (Equilibrium Scheme (15) in Figure 6), in aqueous solution at 25 °C and I = 0.1 M (NaNO3) a.
Table 3. Formation degrees of the isomeric complexes, M(PME2AP)op, M(PME2AP)cl/O, and M(PME2AP)cl/N7 (Equilibrium Scheme (15) in Figure 6), in aqueous solution at 25 °C and I = 0.1 M (NaNO3) a.
M2+% M(PME2AP)op% M(PME2AP)cl/O% M(PME2AP)cl/N7
Mg2+63 ± 737 ± 7~0
Ca2+ 74 ± 14 26 ± 14~0
Mn2+59 ± 841 ± 8~0
Co2+24 ± 414 ± 662 ± 7
Ni2+11 ± 2 4 ± 385 ± 4
Cu2+9 ± 218 ± 673 ± 6
Zn2+13 ± 312 ± 575 ± 6
Cd2+35 ± 7 35 ± 1130 ± 13
a Regarding the error limits (3σ) see Footnote ‘a’ in Table 2. This table is based with permission on our publication in Dalton Transactions [43], copyright 2010, The Royal Society of Chemistry.
Table 4. Stability enhancements (Equation (12)) for the M(PMEDAP) and M(PME-R) complexes and their comparison according to Equation (13). In addition, the extent of chelate formation according to Equilibrium (5) is given as expressed by the dimension-less equilibrium constant KI (Equations (6) and (10)) and the percentages of M(PMEDAP)cl/O (Equation (14)). (Aqueous solution; 25 °C; I = 0.1 M, NaNO3) a.
Table 4. Stability enhancements (Equation (12)) for the M(PMEDAP) and M(PME-R) complexes and their comparison according to Equation (13). In addition, the extent of chelate formation according to Equilibrium (5) is given as expressed by the dimension-less equilibrium constant KI (Equations (6) and (10)) and the percentages of M(PMEDAP)cl/O (Equation (14)). (Aqueous solution; 25 °C; I = 0.1 M, NaNO3) a.
M2+log ΔM(PMEDAP)log ΔM(PME-R)Δ log ΔKI/PMEDAP%M(PMEDAP)cl/O
Mg2+0.17 ± 0.050.16 ± 0.040.01 ± 0.060.48 ± 0.1732 ± 8
Ca2+0.12 ± 0.060.12 ± 0.050.00 ± 0.080.32 ± 0.1924 ± 11
Sr2+0.08 ± 0.040.09 ± 0.05–0.01 ± 0.060.20 ± 0.1217 ± 9
Ba2+0.10 ± 0.060.11 ± 0.05–0.01 ± 0.080.26 ± 0.1921 ± 12
Mn2+0.18 ± 0.060.19 ± 0.06–0.01 ± 0.080.51 ± 0.2034 ± 9
Co2+0.33 ± 0.070.20 ± 0.06(0.13 ± 0.09) b(1.14 ± 0.36)c(53 ± 8) c
Ni2+0.48 ± 0.080.14 ± 0.07(0.34 ± 0.11) b(2.02 ± 0.54)c(67 ± 6) c
Cu2+0.73 ± 0.070.48 ± 0.07(0.25 ± 0.10) b(4.37 ± 0.89) c(81 ± 3) c
Zn2+0.40 ± 0.11d0.29 ± 0.070.11 ± 0.131.51 ± 0.6360 ± 10 d
Cd2+0.32 ± 0.130.30 ± 0.050.02 ± 0.141.09 ± 0.6352 ± 14
a The values are compiled from Tables 4 and 5 of [62]. For the error limits (3δ), see Footnote ‘a’ in Table 2. b These values are larger than zero, indicating that also a metal ion–nucleobase interaction occurs (see text). c Because Δ log Δ > 0, the KI and the % M(PMEDAP)cl/O values also contain a contribution of the M(PMEDAP)cl/N species. d The stability of the Zn(PMEDAP) complex contains probably also a contribution from a nucleobase–Zn2+ interaction (for details see Table 6 and Footnote ‘b’ of Table 6 in [62]).
Table 5. Stability enhancements (Equation (12)) for the M(PME2A6DMAP) and M(PME-R) complexes and their comparison according to Equation (13). In addition, the extent of chelate formation according to Equilibrium (5) is given as expressed by the dimension-less equilibrium constant KI (Equations (6) and (10)) and the percentages of M(PME2A6DMAP)cl/O (Equation (14)). The values for % M(PME-R)cl/O are given for comparison. (Aqueous solution; 25 °C; I = 0.1 M, NaNO3) a.
Table 5. Stability enhancements (Equation (12)) for the M(PME2A6DMAP) and M(PME-R) complexes and their comparison according to Equation (13). In addition, the extent of chelate formation according to Equilibrium (5) is given as expressed by the dimension-less equilibrium constant KI (Equations (6) and (10)) and the percentages of M(PME2A6DMAP)cl/O (Equation (14)). The values for % M(PME-R)cl/O are given for comparison. (Aqueous solution; 25 °C; I = 0.1 M, NaNO3) a.
PE2− = PME2A6DMAP2−
M2+Log
ΔM(PE)
log
ΔM(PME-R)
Δ log ΔKI/PE%
M(PE)cl/O
%
M(PME-R)cl/O
Mg2+0.16 ± 0.040.16 ± 0.040.00 ± 0.060.45 ± 0.1331 ± 631 ± 6
Ca2+0.10 ± 0.050.12 ± 0.05−0.02 ± 0.07 0.26 ± 0.1421 ± 924 ± 9
Sr2+0.11 ± 0.050.09 ± 0.050.02 ± 0.070.29 ± 0.1522 ± 919 ± 9
Ba2+0.08 ± 0.060.11 ± 0.05−0.03 ± 0.080.20 ± 0.1717 ± 1122 ± 9
Mn2+0.19 ± 0.060.19 ± 0.060.00 ± 0.080.55 ± 0.2135 ± 935 ± 9
Co2+0.21 ± 0.070.20 ± 0.060.01 ± 0.090.62 ± 0.2638 ± 1037 ± 9
Ni2+0.17 ± 0.070.14 ± 0.070.03 ± 0.100.48 ± 0.2432 ± 1128 ± 12
Cu2+0.56 ± 0.070.48 ± 0.070.08 ± 0.102.63 ± 0.5972 ± 467 ± 5
Zn2+0.28 ± 0.070.29 ± 0.07−0.01 ± 0.100.91 ± 0.3148 ± 849 ± 8
Cd2+0.25 ± 0.060.30 ± 0.05−0.05 ± 0.080.78 ± 0.2544 ± 850 ± 6
a The values are compiled from Tables 4 and 5 of [61]. For the error limits (3σ), see Footnote ‘a’ in Table 2.
Table 6. Equilibrium constants defined according to the isodesmic model for a non-cooperative self-association for the self-stacking (Equations (18) and (19)) of adenosine and its 5′-phosphates as well as of M(ATP)2− complexes, as determined by 1H NMR shift measurements [90,94] in D2O (27 °C; I = 0.1 − ~2 M NaNO3) (error limits: 2σ).
Table 6. Equilibrium constants defined according to the isodesmic model for a non-cooperative self-association for the self-stacking (Equations (18) and (19)) of adenosine and its 5′-phosphates as well as of M(ATP)2− complexes, as determined by 1H NMR shift measurements [90,94] in D2O (27 °C; I = 0.1 − ~2 M NaNO3) (error limits: 2σ).
No.SystemK (M−1)
1Ado15 ± 2
2AMP2−2.1 ± 0.4
3ADP3−1.8 ± 0.5
4ATP4−1.3 ± 0.2
5Mg(ATP)2−4.0 ± 0.5
6Zn(ATP)2−~11.1 ± 4.5
7Cd(ATP)2−~17
Reproduced with permission from our publication in Pure and Applied Chemistry [21]; copyright 2004, International Union of Pure and Applied Chemistry (IUPAC).
Table 8. Intramolecular equilibrium constants for the formation of the three differently structured Cu(arm)(PMEA) species shown in the Equilibrium Scheme (24), together with the percentages in which these species occur in aqueous solution at 25 °C and I = 0.1 M (NaNO3) a.
Table 8. Intramolecular equilibrium constants for the formation of the three differently structured Cu(arm)(PMEA) species shown in the Equilibrium Scheme (24), together with the percentages in which these species occur in aqueous solution at 25 °C and I = 0.1 M (NaNO3) a.
armKI = KI/tot
(Equation (29))
% Cu(arm)(PMEA)int/tot
(cf. Equation (29))
% Cu(arm)(PMEA)op
(Equation (24))
KI/O
(Equation (27))
KI/st
(Equations (27) and (29)) c
% Cu(arm)(PMEA)cl/O
(Equations (5) and (24)) b
% Cu(arm)(PMEA)st
(Equation (24)) c
bpy29.20 ± 4.8796.69 ± 0.533.31 ± 0.532.89 ± 0.6826.31 ± 4.9210 ± 387 ± 3 (21)
phen53.95 ± 8.8698.18 ± 0.291.82 ± 0.293.17 ± 0.6950.78 ± 8.896 ± 292 ± 2 (25)
a The values listed in the second column are from the fifth column in the lower part of Table 7. These values for KI = KI/tot follow from Equation (29), and % Cu(arm)(PMEA)int/tot is calculated analogously to Equation (14). The values given in the fourth column for % Cu(arm)(PMEA)op follow from 100 − % Cu(arm)(PMEA)int/tot. The constants KI/O of Column 5 are from the upper part in Column 5 of Table 7 (for the corresponding justification see text in Section 7.2); with Equation (26) and the now known values for KI and KI/tot that for KI/st may be calculated (Column 6). All error limits in this table correspond to three times the standard deviation (3σ); they were calculated according to the error propagation after Gauss. b These values were calculated via Equation (26) with KI/O and % Cu(arm)(PMEA)op. c The values for % Cu(arm)(PMEA)st follow from the difference % Cu(arm)(PMEA)int/tot − % Cu(arm)(PMEA)cl/O (cf. Equation (29)); % Cu(arm)(PMEA)st may also be calculated via Equation (27) with KI/st and % Cu(arm)(PMEA)op. The results are the same for both calculation methods (aside from the possible differences in the last digit due to differences in rounding), yet the error limits (which are given in parenthesis) are understandably larger for the second method. This table is reproduced (slight alterations were made) with permission from Table 9 in our publication in Coordination Chemistry Reviews [45]; copyright 1995, Elsevier Science S.A., Lausanne, Switzerland.
Table 9. Extent of chelate formation according to Equilibrium (5) in the Cu(PME) species as quantified by the dimension-less equilibrium constant KI (Equations (6) and (10)) and the percentage of Cu(PME)cl (Equation (14)) in aqueous solution and in water containing 30 or 50% (v/v) 1,4-dioxane together with properties of the mentioned solvents (25 °C; I = 0.1 M, NaNO3)a.
Table 9. Extent of chelate formation according to Equilibrium (5) in the Cu(PME) species as quantified by the dimension-less equilibrium constant KI (Equations (6) and (10)) and the percentage of Cu(PME)cl (Equation (14)) in aqueous solution and in water containing 30 or 50% (v/v) 1,4-dioxane together with properties of the mentioned solvents (25 °C; I = 0.1 M, NaNO3)a.
No.% Dioxane (v/v)mol
Fraction
ε blog ΔCu(PME)
Equation (12)
KI
Equations (6) and (10)
% Cu(PME)cl
Equation (14)
10078.50.48 ± 0.072.02 ± 0.4767 ± 5
2300.08352.70.39 ± 0.041.45 ± 0.2459 ± 4
3500.17535.20.41 ± 0.071.57 ± 0.4061 ± 6
a The values for log ΔCu(PME) (see also Figure 18), KI, and % Cu(PME)cl are from Table 5 of [149] (see also Table 6 of [45]). The information about the solvents is from Table 2 of [149] (see also Table 5 of [45]). For the error limits (3σ) see Footnote ‘a’ in Table 2. b The dielectric constants for the dioxane-water mixtures are interpolated from the data given in [147,148] (see also [45,149]).
Table 10. Effect of 1,4-dioxane (v/v) on the stability of stacks formed in various systems and various solvents a.
Table 10. Effect of 1,4-dioxane (v/v) on the stability of stacks formed in various systems and various solvents a.
Solvent
No.SystemWater30% Diox. a50% Diox. aComments
(1.a)phen (self.)31.1 ± 3.42.63 ± 0.440.63 ± 0.131H NMR; D2O b,c
(1.b)bpy (self.)7.4 ± 2.30.86 ± 0.210.38 ± 0.061H NMR; D2O b,c
(2.a)(phen)(ATP)4−26.8 ± 7.4 (38)4.1 ± 1.1 (4.8)1.6 ± 0.7 (1.8)1H NMR; D2O c,d
(2.b)(bpy)(ATP)4−13.6 ± 3.9 (16)3 ± 2 (3.5)0.3 ± 0.2 (0.4)1H NMR; D2O c,d
(3.a)Cu(phen)(ATP)2−0.56 ± 0.08–0.03 ± 0.05–0.41 ± 0.09Pot. e,f
(3.b)Cu(bpy)(ATP)2−0.33 ± 0.04–0.14 ± 0.05–0.46 ± 0.07Pot. e,f
(4.a)Cu(UTP)2−5.81 ± 0.066.16 ± 0.056.24 ± 0.03Pot. f,g
(4.b)Cu(ATP)2−6.32 ± 0.046.40 ± 0.056.34 ± 0.05Pot. f,g
a The values were collected from Tables 1, 2, and 4 of Ref [127]. For the properties of the dioxane-water mixtures, see Table 9. In the NMR experiments, deuterated dioxane was used. b  K arm self (self-association constant according to the isodesmic model; see Equations (18) and (19)). c The error limits correspond to twice the standard deviation (2σ). 27 °C; I = 0.1 M (NaNO3). d Stability constant K ( arm ) ( ATP ) ATP for the (arm)(ATP)4− adduct; the values in parentheses are corrected for the self-association of arm (= phen or bpy). e Listed are the values for Δ log KCu/arm/ATP (Equation (30)) which quantify Equilibrium (31). f Values based on potentiometric pH titrations [127]. The errors given are three times the standard error of the mean value (3σ); where needed the error propagation after Gauss was employed. 25 °C; I = 0.1 M, NaNO3. g Listed are the log stability constants for the binary Cu(NTP)2− complexes—that is, log K Cu ( NTP ) Cu [127].
Table 11. Extent of chelate formation involving Cu2+ and the adenine residue (AD) of AMP2− and PMEA2− (analogous to Equation (14)) in aqueous solution and in water containing varying percentages (v/v) of 1,4-dioxane at 25 °C and I = 0.1 M (NaNO3). The alterations are quantified via log ΔCu(AMP) and compared with log ΔCu(PMEA) (Equation (12)). These data are based on the stability of the Cu(AMP) and Cu(PMEA) complexes (analogous to Equation (8)) a,b.
Table 11. Extent of chelate formation involving Cu2+ and the adenine residue (AD) of AMP2− and PMEA2− (analogous to Equation (14)) in aqueous solution and in water containing varying percentages (v/v) of 1,4-dioxane at 25 °C and I = 0.1 M (NaNO3). The alterations are quantified via log ΔCu(AMP) and compared with log ΔCu(PMEA) (Equation (12)). These data are based on the stability of the Cu(AMP) and Cu(PMEA) complexes (analogous to Equation (8)) a,b.
ADc% (v/v)
Dioxane
log   K Cu ( AD ) Cu log   K Cu ( AD ) op Cu log ΔCu(AD)KI% Cu(AD)cl
AMP2−03.14 ± 0.012.87 ± 0.08 0.27 ± 0.080.86 ± 0.35 46 ± 10
203.56 ± 0.013.49 ± 0.040.07 ± 0.040.17 ± 0.1115 ± 8
303.86 ± 0.023.82 ± 0.030.04 ± 0.040.10 ± 0.09 9 ± 8
404.30 ± 0.024.19 ± 0.060.11 ± 0.060.29 ± 0.19 22 ± 11
504.73 ± 0.044.45 ± 0.020.28 ± 0.040.91 ± 0.2048 ± 5
PMEA2−03.96 ± 0.043.19 ± 0.060.77 ± 0.074.89 ± 0.9883 ± 3 (49 ± 10) c
304.70 ± 0.054.18 ± 0.030.52 ± 0.062.31 ± 0.4470 ± 4 (26 ± 10) c
505.54 ± 0.074.79 ± 0.030.75 ± 0.084.62 ± 0.9982 ± 3 (54 ± 9) c
a The data for Cu(AMP) are from the Table in [146] and those for Cu(PMEA) from Tables 5 and 6 of [45] (for Cu(AMP) see also [149]). For further details (error limits, etc.), see Footnote ‘f’ of Table 10. b Relevant definitions are given in Equation (6) through Equations (12) and (14), where PE2− = AMP2− or PMEA2−. c AMP2− forms in Cu(AMP) macrochelates with N7 (see, e.g., [21,46,50,53,95,153]). PMEA2− may form several chelates in M(PMEA) complexes (see Section 4.2); we consider here the total amounts of chelates formed, Cu(PMEA)cl/tot, but it may be assumed that Cu(PMEA)cl/O/N3 dominates (Section 3 and Section 4.2). For this reason, the calculated values for Cu(PMEA)cl/O/N3 listed in Table 7 of [45] are also given in parentheses; the small concentrations of Cu(PMEA)cl/N7 are ignored in the above numbers.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sigel, A.; Sigel, H.; Sigel, R.K.O. Coordination Chemistry of Nucleotides and Antivirally Active Acyclic Nucleoside Phosphonates, including Mechanistic Considerations. Molecules 2022, 27, 2625. https://doi.org/10.3390/molecules27092625

AMA Style

Sigel A, Sigel H, Sigel RKO. Coordination Chemistry of Nucleotides and Antivirally Active Acyclic Nucleoside Phosphonates, including Mechanistic Considerations. Molecules. 2022; 27(9):2625. https://doi.org/10.3390/molecules27092625

Chicago/Turabian Style

Sigel, Astrid, Helmut Sigel, and Roland K. O. Sigel. 2022. "Coordination Chemistry of Nucleotides and Antivirally Active Acyclic Nucleoside Phosphonates, including Mechanistic Considerations" Molecules 27, no. 9: 2625. https://doi.org/10.3390/molecules27092625

APA Style

Sigel, A., Sigel, H., & Sigel, R. K. O. (2022). Coordination Chemistry of Nucleotides and Antivirally Active Acyclic Nucleoside Phosphonates, including Mechanistic Considerations. Molecules, 27(9), 2625. https://doi.org/10.3390/molecules27092625

Article Metrics

Back to TopTop