Ondrej
Gutten
and
Lubomír
Rulíšek
*
Institute of Organic Chemistry and Biochemistry, Gilead Sciences Research Center & IOCB, Academy of Sciences of the Czech Republic, Flemingovo nám. 2, 166 10 Praha 6, Czech Republic. E-mail: rulisek@uochb.cas.cz; Fax: +420-220-183-578; Tel: +420-220-183-263
First published on 17th February 2015
The metal-ion selectivity in biomolecules represents one of the most important phenomena in bioinorganic chemistry. The open question to what extent is the selectivity in the complex bioinorganic structures such as metallopeptides determined by the first-shell ligands of the metal ion is answered herein using six model peptides complexed with the set of divalent metal ions (Mn2+, Fe2+, Co2+, Ni2+, Cu2+, Zn2+, Cd2+, and Hg2+) and their various first-shell representations. By calculating the differences among the free energies of complexation of metal ions in these peptides and their model (truncated) systems it is quantitatively shown that the definition of the first shell is paramount to this discussion and revolves around the chemical nature of the binding site. Despite the vast conceivable diversity of peptidic structures, that suggest certain fluidity of this definition, major contributing factors are identified and assessed based on their importance for capturing metal-ion selectivity. These factors include soft/hard character of ligands and various non-covalent interactions in the vicinity of the binding site. The relative importance of these factors is considered and specific suggestions for effective construction of the models are made. The relationship of first-shell models and their corresponding parent peptides is discussed thoroughly, both with respect to their chemical similarity and potential disparity introduced by generally “non-alignable” conformational flexibility of the two systems. It is concluded that, in special cases, this disparity can be negligible and that heeding the chemical factors contributing to selectivity during construction of the model can successfully result in models that retain the affinity profile for various metal ions with high fidelity.
Most of the experimental and computational findings in the area of metal-ion selectivity have been very recently reviewed by Dudev and Lim.6 As highlighted in their review, there are two ‘external’ factors (i.e., independent of the constitution of particular metal-binding site) one always needs to take into account in any considerations of the metal-ion selectivity: (i) average concentrations of metal ions in intracellular and extracellular fluids,8 and (ii) the inherent ‘binding properties’ of metal ions. Concerning the former, it can be reminded that the concentrations of unbound metal ions in cytosol range from millimolar (Na+, K+, Mg2+) through micro- (Mn2+, Fe2+, Ca2+) nano- (Co2+, Ni2+) to femto- (Zn2+) and attomolar (Cu+/Cu2+).9 The latter external factor is mostly exemplified by the Irving–Williams series10 which qualitatively ranks the stabilities of complexes formed by divalent metal ions. This series has its physico-chemical origin in the second ionization enthalpies of the metals11 and predicts the following order:12
Mg2+ < Mn2+ < Fe2+ < Co2+ < Ni2+ < Cu2+ > Zn2+. |
Thus, copper(II) in general forms the most stable complexes with a ‘generic’ set of ligands, followed by Zn2+/Ni2+ whereas complexes of Mg2+ are expected to have the lowest stability constants (β). In tuning the metal-ion selectivity, one has to consider the Irving–Williams series as a ‘baseline’ with respect to which the modifications in the computed or measured metal-ion stabilities are evaluated.
Last, but not least authors of the title review6 clearly summarize most of the areas where metal-ion selectivity is of key importance, notably in function of sodium/potassium/calcium ion channels13–15 and in structure–function of almost all metalloproteins (both regulatory and enzymatic).16–18
Over the last two decades, many studies were reported that addressed the problem of metal-ion selectivity from a computational and quantum chemical perspective.19–28 These most often involved quantum chemical calculations of the small model complexes, both in vacuo and polarized dielectric continuum (to address the effects of the environment, such as solution or protein). Quite often, the results were qualitatively correlated with the experimentally determined stability constants21 or phenomenological information obtained from the abundance of metal ions in the sites of metalloproteins.29–31 In a recent study,32 an attempt was made to conceive a robust and accurate computational protocol that would yield stability constants (β) of selected metal ions (Mn2+, Fe2+, Co2+, Ni2+, Cu2+, Zn2+, Cd2+, and Hg2+) in small model complexes. The protocol was benchmarked on the series of complexes with the known experimental values of β. It was concluded that current computational approaches are likely to suffer from both metal-dependent and ligand-dependent systematic shifts and the straightforward ‘ab initio’ predictions of the ‘absolute’ values of these thermodynamic properties are likely still beyond their grasp. At the same time, it was demonstrated that a relatively easy procedure can be followed that partially accounts for these systematic shifts and the metal-ion selectivity for a particular model site can be in many cases predicted to 1–2 kcal mol−1 accuracy.32
Thus, under an optimistic assumption that computational issues pertinent to the binding of metal ions in small complexes (e.g. metal-binding sites represented by the first-shell ligands) were at least partially solved, a question emanates to what extent are the second-shell effects (which may also include water mediated binding of certain ions to metal binding sites or ion channels) important for the metal-ion selectivity of the site. Some of these issues were addressed in computational studies of ion channels33,34 where second-shell residues are expected to have stabilizing effect on certain type of coordination geometries which, in turn, favor the binding of specific metal ion (and not of its counterpart). In an analogous way, the computations together with statistical survey in Protein Data Bank were used to correlate structure and composition of the outer coordination sphere of metal sites in metalloproteins with those of the inner sphere.35–37 Despite these achievements, the comprehensive study that may yield the robust computational protocol to treat these effects rigorously and quantitatively, is lacking.
This immediately leads to the central question addressed in this work: “How much of the metal-ion selectivity is captured by the first-shell ligand residues of a particular site in metalloprotein?” In answering such a question, one needs to quantitatively assess the stability constants of metal ions in peptidic scaffolds from the first principles. To accomplish such a computational task involves exploring the limits of contemporary quantum chemistry and computational modeling, notably accurate calculations of solvation energies of charged systems.
While the ultimate (and perhaps still distant) goal of this and related studies is the quantitative prediction of stability constants of metal ions in the peptide structures with an accuracy of 1–2 pK (logβ) units, we expect that several important questions can (and will) be raised and answered on the way towards its accomplishment. As mentioned above, these include (i) the suitability of the first-shell-only representation (that is of the model complex [M(Yi)n]c+ where Yi are truncated metal binding amino acids) of the whole metal-binding peptide in calculations of complexation energies and stability constants with special attention paid to their relative values (which, in turn, determine the metal-ion selectivity) and (ii) justifiability of various constraints during the geometry optimization to preserve the original coordination environment in the whole [metal + peptide] complex. The positive answer to the first question would not only provide us with an exciting opportunity to quickly scan various metal sites (both catalytic and functional) for their inherent metal-binding properties and shed further insight into the role of metal ions in biomolecules, but also open new avenues in bottom-up approach for design of novel specific metal-binding sites.
Six model systems studied in this work, schematically depicted in Fig. 1, and described in more details in Section 2.2., are considered. They were carefully selected to represent the considerable part of the ‘spectrum’ of the experimentally observed binding modes. Two of them represent in silico designed metal-binders that were tested experimentally both in gas-phase and solution.39 Unfortunately, no experimental structural information exists on the nature of the binding of metal ions by these peptides. The other four systems represent the continuous metal-binding amino acid sequences in the cores of selected metalloproteins for which the crystal structure was available (cf., Section 2.2.).
![]() | ||
Fig. 1 The six studied peptides. The peptides are abbreviated by the amino-acid residues participating in binding of a metal ion. The larger peptides (CHCC, DHHD, DNDO, DDSOEE) are based on structures found in MESPEUS database.38 The smaller two peptides (CC, MM) come from previous research in our laboratory.39 Vid Section 2.2. for details. |
The series of metal ions with the potential of metal-binding by the peptides includes Mn2+, Fe2+, Co2+, Ni2+, Cu2+, Zn2+, Cd2+, and Hg2+ ions as they represent the most common divalent ions and coincides with the selection in our previous work.32
All metal-ions (Mn, Fe, Co, Ni, Cu, Zn, Cd, Hg) were considered in their +2 oxidation state. For open-shell metal ions in the series (Mn–Cu), high-spin states were considered, in line with our previous work.26,32 These are assumed to be ground electronic states for all of the studied complexes containing mostly low-field ligands. Stuttgart/Dresden pseudopotentials for Cd and Hg ions were applied.47 We consider it as a plausible approximation to relativistic effects, since non-scalar relativistic effects (such as spin–orbit coupling) are estimated to play negligible role for the reaction energetics. The use of ECPs is also in line with our previous benchmark studies.26,32
Free energies of interaction with solvent of all studied species were calculated using the COSMO-RS method48,49 (conductor-like screening model for realistic solvation) as implemented in the COSMOtherm program,50 using the “BP_TZVP_C30_1201.ctd” parametrization file.
The selected PDB structures were required to possess a binding site for a metal ion with all of the binding partners within a relatively short sequence of amino acids. The structures were identified using the information contained in the MESPEUS database.38 The whole metalloprotein (PDB) structure was reduced (truncated) to the minimalistic metal-binding continuous sequence and metal-bound water molecules; terminated with acetyl group at the N-terminus and N-methyl at the C-terminus. Namely, the four peptide sequences considered were CNHEPGTVCPIC (PDB code: 1G71; referred to as CHCC according to the metal-binding residues), DQDKSGFIEEDE (PDB code: 2PAL, referred to as DDSOEE, O standing for backbone carbonyl), DKNGDGE (PDB code: 1IGV; includes 2 water molecules; referred to as DNDO), DHDDVQQHVD (PDB code: 1B9M; includes 1 water molecule; referred to as DHHD).
The remaining two peptides were CGSC (referred to as CC) and MINM (referred to as MM). The binding modes of these two peptides are merely putative. Nevertheless, these two peptides were previously synthesized and binding of ions in the gas phase experimentally determined by mass spectrometry.39 Unfortunately, neither the structural information, nor the solution thermodynamics was obtained for these two systems.
Metal-binding cysteine side chains of CHCC and CC peptides (which were the only Cys residues on both peptides), as well as all aspartate and glutamate side-chains of CHCC, DHHD, DNDO and DDSOEE peptides, were considered in deprotonated form. The studied systems are depicted in Fig. 1. Each system was optimized with each of the eight metal ions, resulting in 6 (peptides) × 8 (metal ions) = 48 structures.
![]() | ||
Scheme 1 Ligands and their representation in individual models. {B} signifies attachment to the peptide backbone, i.e.![]() |
In general, only the substituted hydrogens were optimized in these first-shell models, whereas the remaining atoms were kept frozen in the Cartesian coordinates of their parent system.
The models, by definition, lack the interaction with distant groups and their conformational freedom is unrelated to that of the parent system (see Section 2.4.2). These features are essential for the models to serve their primary purpose – a comprehensive tool for isolating and studying metal-specific interactions of the binding site and the metal ion, and identifying their features and behavior in various settings. The models can be utilized in other ways too, although one has to keep in mind the limitations these features imply.
Firstly, the models can be used in the ‘top-down’ approach, where they can serve as small and cheap models for the selectivity of a specific conformation of the parent system. If more than a single conformation of the binding site is relevant, a single calculation of the parent system for each conformation can be complemented with the model systems for each of the metal ions. This is cheaper than full scale calculation for each metal ion in each conformation.
Alternatively, in the ‘bottom-up’ approach, the models can be used as templates for constructing (e.g. peptidic) scaffolds that can support the binding geometry. Since the model is constructed in a way that it contains almost all selective interactions, the resulting system (e.g. polypeptide chain with the metal-binding site) may retain the selectivity properties identical to those of the model. Rigid scaffolds can also eliminate the need for conformational sampling.
GM,L = Eel + GIS | (1) |
The definition deliberately ignores zero-point vibrational energy and thermal corrections to vibrational, translational and rotational partition function. As discussed thoroughly in previous work,32 these two terms present a non-trivial technical challenge – stemming from finding minima for large (and ideally also solvated) systems, and differences between solvent and gas-phase structures. As these are not expected to be heavily metal-dependent, we preferred to use eqn (1) as the practical (albeit not theoretically pure) approach and plausible approximation to the free energy of the complex. In our previous work, this simple approach somewhat surprisingly yielded better relative stability constants (differences in experimental measure of binding/complexation Gibbs free energies) for a diverse set of model complexes (8 metal ions + 6 sets of ligands) when compared to experimental data.32
In discussing metal-ion selectivity, it is usually advantageous to consider relative values of quantities, rather than their absolute values. These can be brought into spotlight by shifting the quantity equally for all metal ions – an operation that does not change the differences among the values of a quantity. Where applicable, we use quantities shifted in this fashion and indicate the fact by “REL” superscript:
XRELM,L = XM,L − max![]() | (2) |
This realization is in no conflict with the intended purpose of the study – quite the contrary, it allows isolating factors which can influence selectivity but are not part of the model. Consider two cases of a peptide binding a different metal-ion in each case. These two systems differ not only in a metal ion bound, but also in geometry of the peptide/ligands. Comparing such two systems can be viewed as a 2-step transition process: (I) changing a metal ion and (II) relaxing a structure (Fig. 2).
The most important difference between a peptide and its model is the presence of the 2nd shell (i.e. the part of a peptide not included in the model) in the former system. Thus, the overall difference in selectivity of a peptide and its model will be determined by the interaction of the 2nd shell with a metal ion (and the 1st shell) during both of the substeps – giving rise to two contributions: a metal-induced selectivity (step I in Fig. 2) and a structure-induced selectivity (step II in Fig. 2). It is instructive to examine these two contributions separately.
Rather than using e.g. hexahydrated metal ion as reference and examining the parent system and the model separately i.e. (GM,L,peptide − GM,(H2O)6)REL and (GM,L,model − GM,(H2O)6)REL, it is beneficial to consider the difference of these two, i.e.
G2ndM,L,model = (GM,L,peptide − GM,L,model)REL | (3) |
Using the parent system as a reference state, this quantity directly describes the interaction of a metal ion with the 2nd shell. In an ideal case, the model would capture all of the selectivity of a peptide, resulting in this quantity being invariantly zero for all metal ions (M) and systems (L).
In reality, this will not be the case, and the variation of this quantity, examined as average absolute deviation and maximum absolute deviation (eqn (4)), will determine the faithfulness of a model. As will turn out later it might be advantageous that the average is taken over all metal ions (M) but not over different systems (L), as the nature of the systems and even specific structural details will prove to be relevant for the discussion.
![]() | (4) |
Using quantities defined in eqn (3) and (4) and shown in Table 1, it can now be immediately seen that the model works relatively well for CC, MM and CHCC systems (incidentally, these have the largest range of (GM,L − GM,(H2O)6)REL, i.e. the most selective, cf.Table 1). On the other hand, the average absolute deviation is ∼3 kcal mol−1 (maximum absolute deviation >6 kcal mol−1) for DNDO and DDSOEE systems, which is even more significant in light of their low selectivity (range of (GM,L − GM,(H2O)6)REL ∼ 15 kcal mol−1).
System | M2+ | |||||||||
---|---|---|---|---|---|---|---|---|---|---|
Quantity | Mn2+ | Fe2+ | Co2+ | Ni2+ | Cu2+ | Zn2+ | Cd2+ | Hg2+ |
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
|
a All values are in kcal mol−1. Metal-dependent shifts, ΔGcorr, pertaining to the protocol and reference states used and shown to lead to the best computational estimates to the experimental β in our previous work (ref. 32) were applied. ΔGcorr = 0.2, −0.9, 1.9, −1.1, 2.5, −1.6, 0.4, −1.4 kcal mol−1 for Mn2+, Fe2+, Co2+, Ni2+, Cu2+, Zn2+, Cd2+, and Hg2+, respectively. b Lower values of AADL,ALPHA and MADL,ALPHA indicate high retention of metal-ion selectivity in the model. See eqn (4) for definition. c These large values can be traced to a different position of a 2nd-shell charged group, resulting in significantly different interaction with a metal-ion. | ||||||||||
CC | (GM,L,peptide − GM,(H2O)6)REL | 0.0 | −4.7 | −4.7 | −5.2 | −25.7 | −11.2 | −20.6 | −57.8 | |
(GM,L,APLHA − GM,(H2O)6)REL | 0.0 | −5.7 | −4.1 | −3.9 | −23.9 | −11.0 | −18.7 | −54.6 | ||
G2ndM,L,APLHA | −0.9 | 0.0 | −1.5 | −2.2 | −2.7 | −1.1 | −2.8 | −4.2 |
![]() ![]() ![]() |
|
MM | (GM,L,peptide − GM,(H2O)6)REL | 0.0 | −2.5 | −6.5 | −8.6 | N/A | −10.1 | −21.5 | −59.8 | |
(GM,L,ALPHA − GM,(H2O)6)REL | 0.0 | −3.3 | −7.5 | −9.6 | −10.8 | −21.4 | −57.8 | |||
G2ndM,L,APLHA | −1.0 | −0.1 | 0.0 | 0.0 | −1.0 | −0.2 | −1.1 | −2.9 |
![]() ![]() ![]() |
|
DHHD | (GM,L,peptide − GM,(H2O)6)REL | −8.2 | −7.3 | −9.8 | 0.0 | −6.1 | −10.7 | −7.1 | −18.8 | |
(GM,L,ALPHA − GM,(H2O)6)REL | −11.0 | −11.8 | −14.3 | 0.0 | −10.6 | −14.2 | −6.2 | −18.3 | ||
G2ndM,L,APLHA | −1.8 | 0.0 | −0.1 | −4.6 | −0.1 | −1.0 | −5.4 | −5.0 |
![]() ![]() ![]() |
|
DNDO | (GM,L,peptide − GM,(H2O)6)REL | −9.1 | −6.1 | −7.1 | −5.7 | 0.0 | −6.9 | −9.6 | −15.1 | |
(GM,L,ALPHA − GM,(H2O)6)REL | −5.9 | 0.0 | −9.3 | −6.6 | −6.7 | −8.6 | −13.0 | −20.5 | ||
G2ndM,L,APLHA | −10.0 | −12.8 | −4.5 | −5.9 | 0.0 | −5.1 | −3.4 | −1.4 |
![]() ![]() ![]() |
|
CHCC | (GM,L,peptide − GM,(H2O)6)REL | 0.0 | −10.1 | −12.8 | −7.2 | −22.5 | −9.1 | −11.0 | −40.7 | |
(GM,L,ALPHA − GM,(H2O)6)REL | 0.0 | −8.7 | −11.2 | −5.0 | −20.4 | −8.9 | −14.2 | −40.5 | ||
G2ndM,L,APLHA | −3.2 | −4.6 | −4.8 | −5.5 | −5.3 | −3.4 | 0.0 | −3.5 |
![]() ![]() ![]() |
|
DDSOEE | (GM,L,peptide − GM,(H2O)6)REL | −15.3 | −12.4 | 0.0 | −0.7 | −1.8 | −2.9 | −6.1 | −16.6 | |
(GM,L,ALPHA − GM,(H2O)6)REL | −11.8 | −9.7 | −2.6 | 0.0 | −1.2 | −5.2 | −13.4 | −21.9 | ||
G2ndM,L,APLHA | −10.9 | −10.1 | −4.8 | −8.1 | −8.0 | −5.1 | 0.0 | −2.0 |
![]() ![]() ![]() |
G 2ndM,L describes the information about the metal selectivity retained by a model much more clearly than the original quantity, (GM,L − GM,(H2O)6)REL. The set of eight values (one for each metal ion) can be comprehensively reduced to AADL,model and MADL,model, quantities that will be the cornerstone of the ensuing discussion.
The values of G2nd,REL,hybridM,L are collected in Tables S1–S24 (ESI†). These values exhibit a strong trend present in virtually all structures and systems, decreasing gradually across the row to reach a minimum (the most negative value; implying the strongest interaction with the 2nd shell) at Ni2+ or Cu2+, while the maximum (the weakest influence) is usually exhibited by Cd2+ systems. Thus, the models skew the relative free energy between e.g. Cu2+ and Cd2+ in favor of the latter. While this trend follows, at least in partial, the Irving–Williams series, its physico-chemical origin remains unknown to us. Bearing this trend in mind, we choose to present only the AADhybridL,model (and MADhybridL,model).
These values for one of the structures (selected as a representative example), of all six peptide systems and for all four representations are presented in Table 2. Despite the diversity of the systems, most of the AADhybridL,model values are below 1 kcal mol−1 even for the SMALL model, which translate into MADhybridL,model values below ∼2 kcal mol−1 (videTable 2). Examining relationship between these values and the nature of the systems provides insight into the importance of the 2nd shell for selectivity and basis for constructing a robust and well-balanced model.
Model | System | |||||
---|---|---|---|---|---|---|
CC | MM | DHHD | DNDO | CHCC | DDSOEE | |
a All values in kcal mol−1. | ||||||
TINY |
![]() ![]() ![]() |
![]() ![]() ![]() |
![]() ![]() ![]() |
![]() ![]() ![]() |
![]() ![]() ![]() |
![]() ![]() ![]() |
SMALL |
![]() ![]() ![]() |
![]() ![]() ![]() |
![]() ![]() ![]() |
![]() ![]() ![]() |
![]() ![]() ![]() |
![]() ![]() ![]() |
ALPHA |
![]() ![]() ![]() |
![]() ![]() ![]() |
![]() ![]() ![]() |
![]() ![]() ![]() |
![]() ![]() ![]() |
![]() ![]() ![]() |
FULL_AA |
![]() ![]() ![]() |
![]() ![]() ![]() |
![]() ![]() ![]() |
![]() ![]() ![]() |
![]() ![]() ![]() |
![]() ![]() ![]() |
In general, the AADhybridL,model values decrease as the size of the model increases (Tables S1–S24, ESI†):
AADhybridL,FULL_AA < AADhybridL,ALPHA < AADhybridL,SMALL < AADhybridL,TINY |
There are a few remarkable exceptions to this progression, but the magnitude of this discrepancy is diminutive (<few tenths of kcal mol−1) and they are not highlighted in further discussions.
The influence is most visibly exhibited by the TINY model. While systems with predominantly soft residues (CC, MM, CHCC) have AADhybridL,model values of ∼2 kcal mol−1, those of harder ligands are 1 kcal mol−1 or less.
The hardness characteristic can be thought of as the distance at which the ligand can be polarized by its environment. Thus, while all systems converge to the original peptide as we increase the size of the model, the soft ligands will be more sensitive to this change than the hard ligands. Moving from TINY to SMALL model, the AADhybridL,model values drop below 1 kcal mol−1 in all cases, but the decrease continues upon moving further to ALPHA model only for soft ligands; while being almost nonexistent for hard ligands. In general, a similar level of selectivity is achieved by larger model (ALPHA) for sulfur based ligands, while smaller model (SMALL or even TINY) is sufficient for harder ligands.
It should be noted that only hydrogen bonds to atoms directly ligating the central metal ion are considered, as these are the ones that can be expected to chemically (and possibly selectively) influence the binding of a metal ion. Other hydrogen bonds and interactions are unlikely to depend on the identity of the metal ion and can influence the binding only indirectly through structural change (not addressed in the “hybrid” systems studied here).
The following analysis is performed for each relevant peptide (bearing at least some hydrogen bonds to binding atoms) on one of the structures, selected as a representative example. It comes at no surprise (Table S25, ESI†) that there are almost none hydrogen bonds in the smaller representations of the full peptide (up to FULL_AA model). Thus, the effect of a hydrogen bond, if any, is present in the full peptide and FULL_AA model but usually not in the smaller models – which can contribute to different selectivity of a peptide and its model. The “ester” systems, on the other hand, do not possess these hydrogen bonds in either a model or its parent, which results in models being a more faithful representation of its parent system. The FULL_AA model (of a peptide system) usually contains these hydrogen bonds, thus we expect little differences in values of AADhybridL,FULL_AA based on peptides and those based on their “ester” analogue.
The AADhybridL,model values (Table 3) are almost identical in all but few cases, the differences ranging from −0.2 to 0.5 kcal mol−1. The smallest differences are observed in the DHHD, DNDO and DDSOEE systems all of which include hydrogen bonds to hard ligands (carboxyl groups/water molecule/serine alcohol group). On the other hand, AADhybridL,model differ by 0.2–0.5 kcal mol−1 in case of CC and CHCC systems, where the hydrogen bond acceptor is amide bond carbonyl and cysteine thiolate groups, respectively. This amount constitutes 20–70% of the total AADhybridL,model. Moreover, in CHCC system the hydrogen bonds are far from the ideal orientation and in case of strong hydrogen bond, the influence can be expected to be much higher. We thus conclude that the effect of hydrogen bond is potentially significant for strong hydrogen bonds to softer ligands.
Model | System | |||||
---|---|---|---|---|---|---|
CC | DHHD | DNDO | CHCC | DDSOEE | ||
a All values in kcal mol−1. | ||||||
TINY | Peptide | 2.3 | 0.4 | 0.7 | 2.1 | 0.8 |
Ester | 1.9 | 0.3 | 0.6 | 1.9 | 1.0 | |
SMALL | Peptide | 1.0 | 0.3 | 0.6 | 0.9 | 0.6 |
Ester | 0.4 | 0.3 | 0.6 | 0.7 | 0.6 | |
ALPHA | Peptide | 0.7 | 0.3 | 0.4 | 0.9 | 0.6 |
Ester | 0.2 | 0.4 | 0.4 | 0.7 | 0.4 | |
FULL_AA | Peptide | 0.2 | 0.3 | 0.3 | 0.3 | 0.1 |
Ester | 0.2 | 0.2 | 0.2 | 0.4 | 0.2 |
A ligand can also act as a hydrogen bond donor, as exemplified by terminal glutamate of DNDO peptide, which can interact with one of the metal-bound water molecules. The validity of attributing the differences in selectivity to the single carboxyl moiety is confirmed by comparing AADhybridL,model values of original systems to a “mutated” peptide (glutamate → norvaline), presented in Table 4. However, a water molecule is the smallest of the studied ligands, which allows for relative proximity of the charged moiety to a metal ion, which is probably the cause of the observed differences, rather than the existence of the hydrogen bond.
Structure | Glutamate leaning ina | Glutamate leaning outb | ||
---|---|---|---|---|
a Fe2+-optimized structure, see ESI for coordinates. b Co2+-optimized structure, see ESI for coordinates. c Distance between a metal-ion and proximal oxygen of the terminal glutamate side-chain carboxyl group in Ångströms. d Values in kcal mol−1. | ||||
C-terminal residue | Glutamate | Norvaline | Glutamate | Norvaline |
d(M–OOC)c | 4.22 | — | 5.29 | — |
AADhybridL,modeld | ||||
---|---|---|---|---|
TINY | 1.2 | 0.9 | 0.6 | 0.7 |
SMALL | 1.1 | 0.8 | 0.6 | 0.5 |
ALPHA | 0.9 | 0.6 | 0.4 | 0.4 |
FULL_AA | 0.5 | 0.2 | 0.2 | 0.2 |
As in the previous section, we strive to compare systems that differ significantly in this respect (i.e. in the screening of metal-ion), while being as similar as possible in other respects. The most straightforward procedure is to change the number of water molecules in the systems, as these constitute ligands that are not covalently bound to the rest of the system and, thus, no other changes are necessary. Representative structures of peptides CC and DNDO that demonstrate the range of this influence vividly were chosen for discussion. In case of CC peptide, 2 water molecules were added. In case of DNDO, 2 water molecules were removed. In both cases the structures were subsequently optimized, however, in the latter case both the non-optimized and optimized structures were used for selectivity analysis, i.e. for construction and calculation of “hybrid” systems. The results are presented in Table 5.
AADhybridL,modeld | ||||||||
---|---|---|---|---|---|---|---|---|
TINY | 1.7 | 1.5 | 0.6 | 0.9 | 1.1 | 0.8 | 1.3 | 3.1 |
SMALL | 0.8 | 0.8 | 0.6 | 0.7 | 0.9 | 0.6 | 1.1 | 2.3 |
ALPHA | 0.6 | 0.9 | 0.4 | 0.5 | 0.4 | 0.7 | 0.4 | 1.5 |
FULL_AA | 0.2 | 0.1 | 0.2 | 0.4 | 0.3 | 0.3 | 0.4 | 0.9 |
The CC system mode of binding is a distorted plane with the ligands forming a triangle. The two water molecules are added in a line approximately perpendicular to this plane, resulting in a square pyramidal formation (Fig. 4A). There is no part of the peptide being screened by these water molecules. Correspondingly, the values of AADhybridL,model for the two systems are almost identical (see Table 5).
The DNDO system, on the other hand, has an octahedral binding mode with 2 water molecules in cis-positions that “eclipse” two charged groups – a glutamate (already discussed in Section 3.2.2) and a lysine side-chain (Fig. 4B). The removal of the water molecules exposes these charged moieties – which are not part of any of the models – to the metal ion. In all studied cases, this leads to a significant increase in 2nd-shell selectivity. The situation is largely remedied upon optimizing the structures, which distort to almost tetrahedral geometry. The screening is nevertheless deteriorated, compared to the original octahedral structure, and the values of AADhybridL,model are thus higher.
Results for two DNDO structures are shown in Table 5. The original “2-water” systems and even optimized “0-water” systems show similar 2nd-shell selectivity, but there are immense differences in the selectivity of the “0-water” non-optimized structures. The fact that both of these structures were equally submitted to the selectivity analysis (as performed in the previous cases) shows that the difference must pertain to the structure of the peptide. This is alarming, as the two structures are virtually identical. An extensive investigation suggests that the difference is not due to coordination geometry, ligand–metal distances, nor solvent cavity construction. Despite all our efforts, we have not been able to pinpoint the source of this discrepancy. Although the cause of this phenomenon can be artificial, its disappearance in optimized structure does manifest the importance of a properly screened binding site; the absence of which will severely hamper the accuracy regardless of the model used.
A different point of view on the addition/removal of water molecules is to think of it as an exchange of an explicit water molecule for an implicit one or vice versa (it can be reminded that despite advanced treatment of solvation in the COSMO-RS method, the method still belongs to a class of polarized continuum models, PCM). While in the case of the CC system, where the water molecules do not screen a part of the peptide, this seems to leave the selectivity undisturbed, the values of GM,L should still be compared only among systems of identical 1st shell composition.
Table 6 compares the selectivity of one of the hybrid systems, a representative example of metal-induced selectivity, with the fully relaxed structures, representative of the overall selectivity.
Model | System | ||||||
---|---|---|---|---|---|---|---|
CC | MM | DHHD | DNDO | CHCC | DDSOEE | ||
TINY | Hybrid | 2.2 | 3.8 | 0.6 | 0.8 | 2.1 | 1.0 |
Fully relaxed | 2.3 | 3.7 | 1.1 | 2.7 | 2.6 | 2.6 | |
SMALL | Hybrid | 0.9 | 1.5 | 0.3 | 0.5 | 0.8 | 0.7 |
Fully relaxed | 0.9 | 1.0 | 1.8 | 2.9 | 1.4 | 2.8 | |
ALPHA | Hybrid | 0.4 | 0.4 | 0.3 | 0.4 | 0.7 | 0.7 |
Fully relaxed | 1.0 | 0.7 | 2.1 | 3.1 | 1.3 | 3.2 | |
FULL_AA | Hybrid | 0.2 | 0.1 | 0.1 | 0.2 | 0.3 | 0.2 |
Fully relaxed | 0.4 | 0.3 | 1.7 | 3.2 | 0.7 | 2.1 |
In some of the studied systems (DHHD, DNDO, DDSOEE), the disparity is destructive for the utility of the models. The most contributing structural changes are flexible charged groups which cause significantly altered interaction with a metal-ion. In the other cases (CC, MM, CHCC), the structure and selectivity of 2nd shell is virtually unchanged, or even lower after full relaxation. This results from interplay of metal-induced and structure-induced selectivity. While these are comprehensible and practical concepts for discussion, they are not independent but tend to partially cancel each other in these scenarios of diminutive changes of 2nd-shell structure.
As no conformational sampling has been performed, this merely shows that there exists a local minimum where the effect of structure-induced selectivity is minimal. A proper sampling is likely to reveal energetically more relevant structures that do not possess this property. However, this does demonstrate that the assumption of low structure-induced selectivity on the overall selectivity is not inherently incorrect, signaling a green light to the start of search for systems with this characteristics.
This section, together with previous discussion, points out that in cases where the 2nd-shell structure is not sensitive to the identity of bound metal ion, the selectivity is to a large extent confined to the 1st shell. The rigidity condition mentioned in Section 2.4.2 has its origins in this discussion.
For these rigid systems, using a single 1st-shell model structure for each metal ion is sufficient for reproducing selectivity of the parent systems with reasonable accuracy. Looking back at Tables 1 and 2, the satisfactory performance of the models in these cases should now be fully comprehensible.
For ‘bottom-up’ design procedures, the structures of 1st-shell models have to be obtained without prior knowledge of the parent structure. This can be aided by constrained optimization reflecting structural support of the 2nd shell and should consider specifics of the system.
The factors relevant for selectivity identified in this study are of ‘chemical’ and structural nature. Their containment into the described models, that we believe has been sufficiently demonstrated herein, renders them largely independent of the nature of the 2nd shell. As long as the structural and conformational aspects are treated properly, the utility of the models can be expanded outside of the peptide context examined here. Thus, systems such as polymeric scaffolds, cluster models, zeolites, etc. can all be expected to be equally well described by these models. From the methodological perspective, both ‘top-down’ and ‘bottom-up’ approaches to studying and designing metal-binding systems can benefit from the models. The protocol and the models could conceivably be used for e.g. predicting differences in redox potentials of bound metal ions. Caution and further testing is advised, however, as the higher oxidation states can significantly affect not only the quality of treatment but also e.g. the required size of the model.
The metal-induced selectivity of a system is a complex function of multiple factors, that define measure of its localization to a binding site, i.e. the possibility of capturing it using a 1st shell model of the system. It is essential that this model includes all of the ligands bound to a metal-ion. This requirement alone usually ensures that the binding is described reasonably well by the model and that the metal-ion is properly “shielded” from parts of a system further away, bringing the average absolute error to values no higher than few units kcal mol−1. However, some binding geometries can prove to be problematic due to insufficient shielding, e.g. linear or square planar geometry.
Achieving more satisfactory accuracy requires consideration of the nature of ligands. Softer ligands, represented in the realm of metalloproteins by cysteine and methionine residues, require bulkier representation, e.g. the whole side-chain of a residue. Hard ligands, like carboxyl groups, alcohols, amines, can be described at a similar level of accuracy more economically by the functional group plus single methyl group. Imidazole is sufficient to represent histidine side-chain, N-methylmethylamide to represent backbone amide group, and methylamide to represent glutamine and asparagine ligands. This level of description can be expected to result in average absolute deviations below 1 kcal mol−1, which translates to maximum relative error (among studied metal ions) below 3 kcal mol−1. Moreover, most of this relative error derives from bias that appears consistently in the model and can be corrected for heuristically, although this may not be desirable, due to insufficient understanding of this bias.
Non-covalent partners of ligands usually influence retention of selectivity in a negative way. The presence of charged groups can be influential, especially in case of insufficient screening or proximity to metal ion due to small 1st shell ligand, e.g. water molecule. In case of hydrogen bonding the effect is much less significant, influencing average absolute deviation by tenths of kcal mol−1. The strength of the hydrogen bond as well as softness of the acceptor can increase this value.
The overall error in reproducing metal-ion selectivity of a system will stem from inaccuracy of the protocol (∼2 kcal mol−1 in relative free energies32), approximation introduced by the model (<1 kcal mol−1) and assumption of low structure-induced selectivity (the error can be arbitrarily large but, as argued in Section 3.3, conceivably negligible). While some cancellation of errors can be expected, the average errors in the order of units of kcal mol−1 (translating to individual relative errors of roughly up to 5 kcal mol−1) must be expected. As the ranges of free energies of binding are in the order of several tens of kcal mol−1, the protocol holds a promise as a tool for fast estimation of selectivity.
Footnote |
† Electronic supplementary information (ESI) available: The equilibrium geometries of all the molecules studied and Tables S1–S35. See DOI: 10.1039/c4cp04876h |
This journal is © the Owner Societies 2015 |