Probing site-specific calmodulin calcium and lanthanide affinity by grafting

Ca2+ binding is essential for the biological functions of calmodulin (CaM) as a trigger/sensor protein to regulate many biological processes in the Ca2+ -signaling cascade. A challenge in understanding the mechanism of Ca2+ signaling is to obtain site-specific information about the Ca2+ binding properties of individual Ca2+ -binding sites of EF-hand proteins, especially for CaM. In this paper, we report the first estimation of the intrinsic Ca2+ affinities of the four EF-hand loops of calmoduin (I-IV) by individually grafting into the domain 1 of CD2. Taking advantage of the Trp residues in the host protein, we first determined metal-binding affinities for Tb3+, Ca2+, and La3+ for all four grafted EF-loops using Tb3+ aromatic resonance energy transfer. EF-loop I exhibits the strongest binding affinity for Ca2+, La3+, and Tb3+, while EF-loop IV has the weakest metal-binding affinity. EF-loops I-IV of CaM have dissociation constants for Ca2+ of 34, 245, 185, and 814 microM, respectively, with the order I > III approximately equal to II > IV. These findings support a charge-ligand-balanced model in which both the number of negatively charged ligand residues and the balanced electrostatic dentate-dentate repulsion by the adjacent charged residues are two major determinants for the relative Ca2+ -binding affinities of EF-loops in CaM. Our grafting method provides a new strategy to obtain site-specific Ca2+ binding properties and a better estimation of the cooperativity and conformational change contributions of coupled EF-hand proteins.     

Spectroscopic investigation of calcium binding sites in the neurotoxin vipoxin and its components-relation with the X-ray structure

Vipoxin is a neurotoxin from the venom of Vipera ammodytes meridionalis, the most toxic snake in Europe. It is a unique complex of a toxic phospholipase A2 (PLA2) and a non-toxic PLA2-like protein inhibitor (Inh) which probably evolved from the enzyme and reduces its activity and toxicity. The enzymatic activity of Vipoxin is Ca2+-dependent and the interaction of this metal ion with the neurotoxic complex and its separated components was investigated using the fluorescent probe ANS. Vipoxin binds two calcium ions, one per each subunit. The X-ray model of the Ca2+-free neurotoxin shows that the potential metal-binding sites require minor structural changes to bind calcium. The dissociation constants K(2+)Ca of the calcium complexes of Vipoxin and its components, PLA2 and Inh, were determined to be 16, 10 and 9 mM, respectively. The affinity for calcium of Vipoxin is reduced in comparison to those of PLA2 and Inh. The X-ray model shows that the potential Ca2+-binding sites in the two components are partially 'shielded' in the complex. The affinity of the neurotoxin to Sr2+ and Ba2+ is lower and the respective K(2+)Ca are 20 and 30 mM. The saturation of Ca2+-binding sites increased the melting point Tm of Vipoxin by 11 degrees C and the activation energy for the thermal deactivation of the excited tryptophans Ea by 11 kJ mol(-1) x Ca2+ is important not only for the enzymatic activity of Vipoxin but also for its thermostability.     

The use of ESI-MS to probe the binding of divalent cations to calmodulin

Proteins have evolved with distinct sites for binding particular metal ions. This allows metalloproteins to perform a myriad of specialized tasks with conformations tailor-made by the combination of its primary sequence and the effect on this of the ligated metal ion. Here we investigate the selectivity of the calcium trigger protein calmodulin for divalent metal ions. This ubiquitous and highly abundant protein exists in equilibrium between its apo and its holo form wherein four calcium ions are bound. Amongst its many functions, calmodulin modulates the calcium concentration present in cells, but this functional property renders it a target for competition from other metal ions. We study the competition posed by four other divalent cations for the calcium binding sites in calmodulin using electrospray ionization mass spectrometry (ESI-MS). We have chosen two other group II cations Mg²⁺, Sr²⁺, and two heavy metals Cd²⁺, Pb²⁺. The ease with which each of these metals binds to apo and to holo CaM[4Ca] is described. We find that each metal ion has different properties with respect to calmodulin binding and competition with calcium. The order of affinity for apo CaM is Ca²⁺ ≫ Sr²⁺ ∼ Mg²⁺ > Pb²⁺ ∼ Cd²⁺. In the presence of calcium the affinity alters to Pb²⁺ > Ca²⁺ > Cd²⁺ > Sr²⁺ > Mg²⁺. Once complexes have been formed between the metal ions and protein (CaM:[xM]) we investigate whether the structural change which must accompanies calcium ligation to allow target binding takes place for a given CaM:[xM] system. We use a 20 residue target peptide, which forms the CaM binding site within the enzyme neuronal nitric-oxide synthase. Our earlier work (Shirran et al. 2005) [1] has demonstrated the particular selectivity of this system for CaM:4Ca²⁺. We find that along with Ca²⁺ only Pb²⁺ forms complexes of the form CaM:4M²⁺:nNOS. This work demonstrates the affinity for calcium above all other metals, but also warns about the ability of lead to replace calcium with apparent ease.     

Determination of the metal-binding cooperativity of wild-type and mutant calbindin D9K by electrospray ionization mass spectrometry

Since the initial reports showing the ability of electrospray ionization mass spectrometry (ESI-MS) to study intact noncovalent biomolecular complexes, an increasing number of uses for this technique in studying biochemical systems is emerging. We have investigated the ability of ESI-MS to characterize the metal-binding properties of calcium (Ca2+) binding proteins by studying the incorporation of Ca2+ and cadmium (Cd2+) into wild-type and mutant calbindin D9K. ESI-MS showed that wild-type calbindin D9K binds two Ca2+ ions with similar affinities while the binding of two Cd2+ ions is sequential, as is the binding of the two Ca2+ or Cd2+ ions to the N56A mutant of calbindin. The binding of Ca2+ to the wild-type protein was clearly seen to be cooperative. These results demonstrate the potential efficacy of ESI-MS to discriminate between cooperative and independent site metal binding to metalloproteins.     

钙调素与重金属Pb2+结合反应的方波极谱与循环伏安法研究

采用方波极谱法研究了重金属Pb2+与钙调素(CaM)的结合反应, 直接检测到Pb2+-CaM配合物的存在, 并进一步利用循环伏安法研究了Pb2+-CaM的电极反应. 在pH=6.5时, 用方波极谱法在Pb2+-CaM体系中检测出2个还原峰, 峰电位分别为-0.44~-0.47 V和-0.73~-0.77 V, 说明在Pb2+-CaM体系中铅有2种存在形式, -0.44~-0.47 V的还原峰对应于游离态Pb2+, 电位更负的还原峰对应于配合物[Pb2+-CaM]. 2个还原峰的峰电流均随着cPb2+/cCaM比值增大而增大; 至cPb2+/cCaM≥10后, 配合物[Pb2+-CaM]的峰电流基本不再变化, 而游离态Pb2+的峰电流则继续增大. 利用极谱滴定曲线的拐点可判断出Pb2+在CaM中有10个结合位点. 进一步的测量结果表明, 循环伏安曲线出现游离态Pb2+的氧化峰和还原峰, 而络合态的[Pb2+-CaM]只有其还原峰, 反向电压扫描时不出现阳极波, 即没有相对应的氧化峰出现.     

Binding studies of nNOS-active amphibian peptides and Ca2+ calmodulin, using negative ion electrospray ionisation mass spectrometry

Amphibian peptides which inhibit the formation of nitric oxide by neuronal nitric oxide synthase (nNOS) do so by binding to the protein cofactor, Ca2+calmodulin (Ca2+CaM). Complex formation between active peptides and Ca2+CaM has been demonstrated by negative ion electrospray ionisation mass spectrometry using an aqueous ammonium acetate buffer system. In all cases studied, the assemblies are formed with a 1:1:4 calmodulin/peptide/Ca2+ stoichiometry. In contrast, the complex involving the 20-residue binding domain of the plasma Ca2+ pump C20W (LRRGQILWFRGLNRIQTQIK-OH) with CaM has been shown by previous two-dimensional nuclear magnetic resonance (2D NMR) studies to involve complexation of the C-terminal end of CaM. Under identical conditions to those used for the amphibian peptide study, the ESI complex between C20W and CaM shows specific 1:1:2 stoichiometry. Since complex formation with the studied amphibian peptides requires Ca2+CaM to contain its full complement of four Ca2+ ions, this indicates that the amphibian peptides require both ends of the CaM to effect complex formation. Charge-state analysis and an H/D exchange experiment (with caerin 1.8) suggest that complexation involves Ca2+CaM undergoing a conformational change to a more compact structure.     

Interaction of different metal ions with carboxylic acid group: a quantitative study

The binding strength of the carboxylic acid group (-COOH) with different divalent metal ions displays considerable variation in arachidic acid (AA) thin films. It is considered that in AA thin films the metal ions straddle the hydrophilic regions of the stacked bilayers of AA molecules via formation of carboxylates. In this study first the uptake of different divalent cations in films of AA is estimated by atomic absorption spectroscopy (AAS). Through the amount of cation uptake, it is found that the strength of binding of different cations varies as Ca2+>Co2+>Pb2+>Cd2+. Variation in the binding strength of different ions is also manifested in experiments where AA thin films are exposed to metal ion mixtures. The higher binding strength of AA with certain metal ions when exposed individually, as well as the preference over the other metal ions when exposed to mixtures, reveal some interesting deviation from the expected behavior based on considerations of ionic radii. For example, Pb2+ is always found to bind to AA much more strongly than Cd2+ even though the latter has smaller ionic radius, indicating that other factors also play an important role in governing the binding strength trends apart from the effects of ionic radii. Then, to get a more meaningful knowledge regarding the binding capability, first-principles calculations based on density functional theory have been applied to study the interaction of different cations with the simplest carboxylic acid, acetic acid, that can result in formation of metal diacetates. Their electronic and molecular structures, cohesive energies, and stiffness of the local potential energy well at the cation (M) site are determined and attempts are made to understand the diversity in geometry and the properties of binding of different metal ions with -COOH group. We find that the calculated M-O bond energies depend sensitively on the chemistry of M atom and follow the experimentally observed trends quite accurately. The trends in M-O bond energies and hence the total M-acetate binding energy trends can actually be related to their molecular structures that fall into different categories: Ca and Cd have tetrahedral coordination; Fe, Ni, and Co exhibit planar 4-fold coordination; and Pb is off-centered from the planar structure (forming pyramidal structure) due to its stereochemically active lone pair of electrons.     

Complexes of bivalent metal cations in electrospray mass spectra of common organic compounds

In the standard electrospray ionization mass spectra of many common, low molecular mass organic compounds dissolved in methanol, peaks corresponding to ions with formula [3M + Met](2+) (M = organic molecule, Met = bivalent metal cation) are observed, sometimes with significant abundances. The most common are ions containing Mg(2+), Ca(2+) and Fe(2+). Their presence can be easily rationalized on the basis of typical organic reaction work-up procedures. The formation of [3M + Met](2+) ions has been studied using N-FMOC-proline methyl ester as a model organic ligand and Mg(2+), Ca(2+), Sr(2+), Ba(2+), Fe(2+), Ni(2+), Mn(2+), Co(2+) and Zn(2+) chlorides or acetates as the sources of bivalent cation. It was found that all ions studied form [3M + Met](2+) complexes with N-FMOC-proline methyl ester, some of them at very low concentrations. Transition metal cations generally show higher complexation activity in comparison with alkaline earth metal cations. They are also more specific in the formation of [3M + Met](2+) complexes. In the case of alkaline earth metal cations [2M + Met](2+) and [4M + Met](2+) complex ions are also observed. It has been found that [3M + Met](2+) complex ions undergo specific fragmentation at relatively low energy, yielding fluorenylmethyl cation as a major product. [M + Na](+) ions are much more stable and their fragmentation is not as specific.     

Metal ion binding to parvalbumin. A proton NMR study

The 1H NMR spectra of carp parvalbumin saturated with Ca2+, Cd2+, La3+ and Lu3+ were compared, using 2D 1H NMR techniques as well as conventional 1H NMR spectra. The Ca2+ and Cd2+ saturated parvalbumin (with both high affinity Ca2+-binding sites occupied) gave rise to very similar spectra. This shows that these two species have almost identical protein conformations. The 1H NMR spectrum from the Ln3+ saturated parvalbumins deviated from the other two and it was therefore concluded that Cd2+ is a better probe for Ca2+ than Ln3+ in parvalbumin and probably also for related calcium binding proteins. The addition of excess of divalent metal ions, such as Mg2+ or Ca2+, causes small changes in the chemical shift of some methyl resonances. This is presumably caused by binding of these metal ions to a third site close to the CD site which is made up of the carboxylic groups from Glu 60 and Asp 61.     

A combined experimental and theoretical study of divalent metal ion selectivity and function in proteins: application to E. coli ribonuclease H1

Structural and thermodynamic aspects of alkaline earth metal dication (Mg(2+), Ca(2+), Sr(2+), Ba(2+)) binding to E. coli ribonuclease H1 (RNase H1) have been investigated using both experimental and theoretical methods. The various metal-binding modes of the enzyme were explored using classical molecular dynamics simulations, and relative binding free energies were subsequently evaluated by free energy simulations. The trends in the free energies of model systems based on the simulation structures were subsequently verified using a combination of density functional theory and continuum dielectric methods. The calculations provide a physical basis for the experimental results and suggest plausible role(s) for the metal cation and the catalytically important acidic residues in protein function. Magnesium ion indirectly activates water attack of the phosphorus atom by freeing one of the active site carboxylate residues, D70, to act as a general base through its four first-shell water molecules, which prevent D70 from binding directly to Mg(2+). Calcium ion, on the other hand, inhibits enzyme activity by preventing D70 from acting as a general base through bidentate interactions with both carboxylate oxygen atoms of D70. These additional interactions to D70, in addition to the D10 and E48 monodentate interactions found for Mg(2+), enable Ca(2+) to bind tighter than the other divalent ions. However, a bare Mg(2+) ion with two or less water molecules in the first shell could bind directly to the three active-site carboxylates, in particular D70, thus inhibiting enzymatic activity. The present analyses and results could be generalized to other members of the RNase H family that possess the same structural fold and show similar metal-binding site and Mg(2+)-dependent activity.     

The significance of the alkene size and the nature of the metal ion in metal-alkene complexes: a theoretical study

Cation interactions with π-systems are a problem of outstanding contemporary interest and the nature of these interactions seems to be quite different for transition and main group metal ions. In this paper, we have systematically analyzed the contrast in the bonding of Cu(+) and main group metal ions. The molecular structures and energetics of the complexes formed by various alkenes (A = C(n)H(2n), n = 2-6; C(n)H(2n- 2), n = 3-8 and C(n)H(2n + 2), n = 5-10) and metal ions (M = Li(+), Na(+), K(+), Ca(2+), Mg(2+), Cu(+) and Zn(2+)) are investigated by employing ab initio post Hartree-Fock (MP2/6-311++G**) calculations and are reported in the current study. The study, which also aims to evaluate the effect of the size of the alkyl portion attached to the π-system on the complexation energy, indicates a linear relationship between the two. The decreasing order of complexation energy with various metal ion-alkene complexes follows the order Zn(2+)-A > Mg(2+)-A > Ca(2+)-A > Cu(+)-A > Li(+)-A > Na(+)-A > K(+)-A. The increased charge transfer and the electron density at (3,-1) intermolecular bond critical point corroborates well with the size of the π-system and the complexation energy. The observed deviation from the linear dependency of the Cu(+)-A complexes is attributed to the dπ→π* back bonding interaction. An energy decomposition analysis via the reduced variational space (RVS) procedure was also carried out to analyze which component among polarization, charge transfer, coulomb and exchange repulsion contributes to the increase in the complexation energy. The RVS results suggest that the polarization component significantly contributes to the increase in the complexation energy when the alkene size increases.     

Separation and evaluation of soybean protein hydrolysates prepared by immobilized metal ion affinity chromatography with different metal ions

Metal ion affinity chromatography is widely used to purify peptides on the basis of the dissimilarities of their amino acids. However, researchers are interested in the separation differences between different metal ions in this method. In our study, four kinds of commonly used metal ions are compared by the amount of immobilized metal ion on iminodiacetic acid-Sepharose and binding amount of soybean peptide to immobilized iminodiacetic acid-Mn(+) adsorbents and evaluated by high-performance liquid chromatography (HPLC) profiles. The results show that due to the different adsorption behaviors of metal ions, the binding ability order of soybean protein peptide on the column should be Fe(3+) > Cu(2+) > Zn(2+) > Ca(2+). The HPLC profiles show that peptides adsorbed by four kinds of metal ions display similar strong hydrophobic characteristics.     

Electrospray ionization mass spectrometry studies of noncovalent myosin VI complexes reveal a new specific calmodulin binding site

Among the myosin superfamily, myosin VI differs from all others by a reverse directionality and a particular motility. Little structural information is available for myosin VI. It is known that it binds one calmodulin (CaM) by means of a single "IQ motif" and that myosin VI contains a specific insert located at the junction between the motor domain (MD) and the lever arm, likely to play a critical role for the unusual motility previously observed. Electrospray ionization mass spectrometry (MS) was used to determine the CaM and Ca2+ stoichiometries in several myosin VI constructs. In particular, the experimental conditions required for the observation of multiprotein/Ca2+ noncovalent assemblies are detailed for two truncated MD constructs (less than 20 kDa) and for three full MD constructs (more than 90 KDa). The specificity of the detected stoichiometries is discussed for each construct and the resolving power of Time of Flight mass spectrometry is stressed, in particular for the detection of metal ions binding to high molecular weight complexes. MS reveals a new CaM binding site for myosin VI and highlights a different behavior for the five myosin VI constructs versus Ca2+ binding. In addition to these stoichiometry based experiments, gas-phase dissociation analyses on intact complexes are described. They reveal that Ca2+ transfer between protein partners occurs during the dissociation process for one construct with a full MD. Charge-transfer and dissociation behavior has allowed to draw structural assumptions for the interaction of the MD with the CaM N-terminal lobe.     

Relative calcium-binding strengths of amino acids determined using the kinetic method

We have measured the relative calcium-binding energies of amino acids using tandem mass spectrometry of Ca(2+)-bound trimeric amino acids. Although calcium-bound dimeric amino acid complexes coordinated too strongly to allow observation of the two competing dissociation products (calcium-bound monomeric ions) required for analysis of their metal binding affinities using the conventional kinetic method, the Ca(2+)-bound trimeric cluster ions dissociated readily to form dimeric cluster ions through simple ligand losses. The calcium-binding energies were obtained by comparing the ratio of the [Ca(2+)(A(1))(2) - H(+)](+) and [Ca(2+) (A(1))(A(2)) - H(+)](+) ions that dissociated from the [Ca(2+) (A(1))(2)(A(2)) - H(+)](+) ion and the ratio of the [Ca(2+)(A(2))(2) - H(+)](+) and [Ca(2+)(A(1)) (A(2)) - H(+)](+) ions that dissociated from the [Ca(2+)(A(1))(A(2))(2) - H(+)](+) ion, where A(1) and A(2) represent two amino acids. The energies deduced from this analysis represent the relative average binding energies of complexes having the form [Ca(2+)(A(1))(2) - H(+)](+). The relative Ca(2+)-binding strengths of the alpha-amino acid complexes follow the order Cys < Ser < Thr < Ile < Leu < Val < Gly < Ala < Pro < Phe < Met < Tyr < Asn < His < Gln < Trp < Lys < Arg. To our knowledge, this report provides the first example of using kinetic methods to determine the relative binding strengths of divalent metal-amino acid complexes.     

Interactions of the hormone oxytocin with divalent metal ions

The interaction of the cyclic nonapeptide oxytocin (OT) with a number of alkaline earth and divalent transition metal ions (X(2+)) was examined employing mass spectrometry (MS) and ion mobility spectrometry (IMS) techniques in combination with molecular dynamics (MD) and density functional theory (DFT) calculations. Under acidic conditions it was found that OT exhibits an exceptionally strong affinity for all divalent metal ions resulting in strong [OT + X](2+) peaks in the mass spectrum. Under basic conditions only Cu(2+) and Ni(2+)-OT complexes were detected and these were singly, doubly, triply, or quadruply deprotonated. Collision-induced dissociation of the [OT - 3H + Cu](-) complex yielded exclusively C-terminal Cu(2+)-containing fragments (Cu(2+)fragment(3-)), suggesting that the Cu(2+) ligation site includes deprotonated C-terminal backbone amide nitrogen atoms and the N-terminal amino nitrogen atom in [OT - 3H + Cu](-). MD and DFT calculations indicate a square-planar complex is consistent with these observations and with experimental collision cross sections. MD and DFT calculations also indicate either an octahedral or trigonal-bipyramidal complex between Zn(2+) and OT is lowest in energy with carbonyl oxygens being the primary ligation sites. Both complexes yield cross sections in agreement with experiment. The biological impact of the structural changes induced in OT by divalent metal ion coodination is discussed.     

Binding of Ca2+, Mg2+, and heparin by human serum amyloid P component in affinity capillary electrophoresis

Human serum amyloid P component (SAP) is a glycoprotein circulating in the blood and found in association with all types of amyloid (malfolded potein aggregates) examined so far. Despite uncertainties regarding the precise function of SAP in vivo, the lectin-like properties of this Ca(2+)-activated protein with affinity for anionic saccharides and malfolded proteins are well known. The propensity to form homomeric penta- or decamers in solution and the selfaggregation in the presence of Ca(2+) as well as the tendency of SAP to attach to uncoated fused silica have precluded the analysis of SAP by microelectrophoretic methods. We now work out conditions to characterize the binding of Ca(2+) and Mg(2+) and the binding of heparin to SAP in the presence of divalent metal ions by ACE. The results show a strong binding of heparin (sub-muM apparent dissociation constants) even in the abscence of Ca(2+) at low ionic strength, pH 8.2. Also, a selective interaction with Ca(2+) compared with Mg(2+) is demonstrated. The approach will further the use of microelectrophoretic methods to examine the interactions of SAP with ligands of putative pathophysiological relevance such as lipopolysaccharides and misfolded proteins.     

OFF-OFF-ON switching of fluorescence and electron transfer depending on stepwise complex formation of a host ligand with guest metal ions

Stepwise complex formation is observed between 2,3,5,6-tetrakis(2-pyridyl)pyrazine (TPPZ) and a series of metal ions (M(n+) = Sc3+, Y3+, Ho3+, Eu3+, Lu3+, Nd3+, Zn2+, Mg2+, Ca2+, Ba2+, Sr2+, Li+), where TPPZ forms a 2:1 complex [(TPPZ)2-M(n+)] and a 1:1 complex [TPPZ-M(n+)] with Mn+ at low and high concentrations of metal ions, respectively. The fluorescence intensity of TPPZ begins to increase at high concentrations of metal ions, when the 2:1 (TPPZ)2-M(n+) complex is converted to the fluorescent 1:1 TPPZ-M(n+) complex. This is regarded as an "OFF-OFF-ON" fluorescence sensor for metal ions depending on the stepwise complex formation between TPPZ and metal ions. The fluorescence quantum yields of the TPPZ-M(n+) complex vary depending on the metal valence state, in which the fluorescence quantum yields of the divalent metal complexes (TPPZ-M2+) are much larger than those of the trivalent metal complexes (TPPZ-M3+). On the other hand, the binding constants of (TPPZ)2-M(n+) (K1) and TPPZ-M(n+) (K2) vary depending on the Lewis acidity of metal ions (i.e., both K1 and K2 values increase with increasing Lewis acidity of metal ions). Sc3+, which acts as the strongest Lewis acid, forms the (TPPZ)2-Sc3+ and TPPZ-Sc3+ complexes stoichiometrically with TPPZ. In such a case, "OFF-OFF-ON" switching of electron transfer from cobalt(II) tetraphenylporphyrin (CoTPP) to O2 is observed in the presence of Sc3+ and TPPZ depending on the ratio of Sc3+ to TPPZ. Electron transfer from CoTPP to O2 occurs at Sc3+ concentrations above the 1:2 ratio ([Sc3+]/[TPPZ]0 > 0.5), when the (TPPZ)2-Sc3+ complex is converted to the TPPZ-Sc3+ complex and TPPZ-(Sc3+)2, which act as promoters of electron transfer (ON) by the strong binding of O2*- with Sc3+. In sharp contrast, no electron transfer occurs without metal ion (OFF) or in the presence at Sc3+ concentrations below the 1:2 ratio (OFF), when the (TPPZ)2-Sc3+ complex has no binding site available for O2*-.     

Factors governing the substitution of La3+ for Ca2+ and Mg2+ in metalloproteins: a DFT/CDM study

Trivalent lanthanide cations are extensively being used in biochemical experiments to probe various dication-binding sites in proteins; however, the factors governing the binding specificity of lanthanide cations for these binding sites remain unclear. Hence, we have performed systematic studies to evaluate the interactions between La3+ and model Ca2+ - and Mg2+ -binding sites using density functional theory combined with continuum dielectric methods. The calculations reveal the key factors and corresponding physical bases favoring the substitution of trivalent lanthanides for divalent Ca2+ and Mg2+ in holoproteins. Replacing Ca2+ or Mg2+ with La3+ is facilitated by (1) minimizing the solvent exposure and the flexibility of the metal-binding cavity, (2) freeing both carboxylate oxygen atoms of Asp/Glu side chains in the metal-binding site so that they could bind bidentately to La3+, (3) maximizing the number of metal-bound carboxylate groups in buried sites, but minimizing the number of metal-bound carboxylate groups in solvent-exposed sites, and (4) including an Asn/Gln side chain for sites lined with four Asp/Glu side chains. In proteins bound to both Mg2+ and Ca2+, La3+ would prefer to replace Ca2+, as compared to Mg2+. A second Mg2+-binding site with a net positive charge would hamper the Mg2+ --> La3+ exchange, as compared to the respective mononuclear site, although the La3+ substitution of the first native metal is more favorable than the second one. The findings of this work are in accord with available experimental data.     

Modelling peptide-metal dication interactions: formamide-Ca(2+) reactions in the gas phase

The collision induced dissociation of formamide-Ca(2+) complexes produced in the gas phase through nanoelectrospray ionization yields as main products ions [CaOH](+), [HCNH](+), [Ca(NH(2))](+), HCO(+) and [Ca(NH(3))](2+) and possibly [Ca(H(2)O)](2+) and [C,O,Ca](2+), the latter being rather minor. The mechanisms behind these fragmentation processes have been established by analyzing the topology of the potential energy surface by means of B3LYP calculations carried out with a core-correlated cc-pWCVTZ basis set. The Ca(2+) complexes formed by formamide itself and formimidic acid play a fundamental role. The former undergoes a charge separation reaction yielding [Ca(NH(2))](+) + HCO(+), and the latter undergoes the most favorable Coulomb explosion yielding [Ca-OH](+) + [HCNH](+) and is the origin of a multistep mechanism which accounts for the observed loss of water and HCN. Conversely, the other isomer of formamide, amino(hydroxyl)carbene, does not play any significant role in the unimolecular reactivity of the doubly charged molecular cation.     

Effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) and water coordination on the structure of glycine and zwitterionic glycine

Interactions between metal ions and amino acids are common both in solution and in the gas phase. Here, the effect of metal ions and water on the structure of glycine is examined. The effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) and water on structures of Gly.Mn+(H2O)m and GlyZwitt.Mn+(H2O)m (m = 0, 2, 5) complexes have been determined theoretically by employing the hybrid B3LYP exchange-correlation functional and using extended basis sets. Selected calculations were carried out also by means of CBS-QB3 model chemistry. The interaction enthalpies, entropies, and Gibbs energies of eight complexes Gly.Mn+ (Mn+ = Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) were determined at the B3LYP density functional level of theory. The computed Gibbs energies DeltaG degrees are negative and span a rather broad energy interval (from -90 to -1100 kJ mol(-1)), meaning that the ions studied form strong complexes. The largest interaction Gibbs energy (-1076 kJ mol(-1)) was computed for the NiGly2+ complex. Calculations of the molecular structure and relative stability of the Gly.Mn+(H2O)m and GlyZwitt.Mn+(H2O)m (Mn+ = Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+; m = 0, 2, and 5) systems indicate that in the complexes with monovalent metal cations the most stable species are the NO coordinated metal cations in non-zwitterionic glycine. Divalent cations Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+ prefer coordination via the OO bifurcated bonds of the zwitterionic glycine. Stepwise addition of two and five water molecules leads to considerable changes in the relative stability of the hydrated species. Addition of two water molecules at the metal ion in both Gly.Mn+ and GlyZwitt.Mn+ complexes reduces the relative stability of metallic complexes of glycine. For Mn+ = Li+ or Na+, the addition of five water molecules does not change the relative order of stability. In the Gly.K+ complex, the solvation shell of water molecules around K+ ion has, because of the larger size of the potassium cation, a different structure with a reduced number of hydrogen-bonded contacts. This results in a net preference (by 10.3 kJ mol(-1)) of the GlyZwitt.K+H2O5 system. Addition of five water molecules to the glycine complexes containing divalent cations Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+ results in a net preference for non-zwitterionic glycine species. The computed relative Gibbs energies are quite high (-10 to -38 kJ mol(-1)), and the NO coordination is preferred in the Gly.Mn+(H2O)5 (Mn+ = Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) complexes over the OO coordination.