An Integrated Mass Spectrometry and Molecular Dynamics Simulations Approach Reveals the Spatial Organization Impact of Metal-Binding Sites on the Stability of Metal-Depleted Metallothionein-2 Species

Mammalian metallothioneins (MTs) are a group of cysteine-rich proteins that bind metal ions in two α- and β-domains and represent a major cellular Zn(II)/Cu(I) buffering system in the cell. At cellular free Zn(II) concentrations (10-11-10-9 M), MTs do not exist in fully loaded forms with seven Zn(II)-bound ions (Zn7MTs). Instead, MTs exist as partially metal-depleted species (Zn4-6MT) because their Zn(II) binding affinities are on the nano- to picomolar range comparable to the concentrations of cellular Zn(II). The mode of action of MTs remains poorly understood, and thus, the aim of this study is to characterize the mechanism of Zn(II) (un)binding to MTs, the thermodynamic properties of the Zn1-6MT2 species, and their mechanostability properties. To this end, native mass spectrometry (MS) and label-free quantitative bottom-up and top-down MS in combination with steered molecular dynamics simulations, well-tempered metadynamics (WT-MetaD), and parallel-bias WT-MetaD (amounting to 3.5 μs) were integrated to unravel the chemical coordination of Zn(II) in all Zn1-6MT2 species and to explain the differences in binding affinities of Zn(II) ions to MTs. Differences are found to be the result of the degree of water participation in MT (un)folding and the hyper-reactive character of Cys21 and Cys29 residues. The thermodynamics properties of Zn(II) (un)binding to MT2 are found to differ from those of Cd(II), justifying their distinctive roles. The potential of this integrated strategy in the investigation of numerous unexplored metalloproteins is attested by the results highlighted in the present study.

Metal binding properties of zinc fingers with a naturally altered metal binding site

Zinc fingers (ZFs) are among the most abundant motifs found in proteins, and are commonly known for their structural role. Classical ZFs (CCHH) are part of the transcription factors that participate in DNA binding. Although biochemical studies of classical ZFs have a long history, there is limited knowledge about the sequential and structural diversity of ZFs. We have found that classical ZFs, with metal binding sites consisting of amino acids other than conserved Cys or His residues, are frequently encoded in the human genome, and we refer to these peptides as ZFs with a naturally altered metal binding site. The biological role of the altered ZFs remains undiscovered. In this study, we characterized nine natural XCHH, CXHH, CCXH and CCHX ZFs in terms of their Zn(ii) and Co(ii) binding properties, such as complex stoichiometry, spectroscopic properties and metal-to-peptide affinity. We revealed that XCHH and CXHH ZFs form ML complexes that are 4-5 orders of magnitude weaker in comparison to CCHH ZFs. Nevertheless, spectroscopic studies demonstrate that, depending on the altered position, they may adopt an open coordination geometry with one or two water molecules bound to a central metal ion, which has not been demonstrated in natural ZFs before. Stability data show that both CCXH and CCHX peptides have high Zn(ii) affinity (with a K_d of 10^-9 to 10^-11 M), suggesting their potential biological function. This study is a comprehensive overview of the relationship between the sequence, structure, and stability of ZFs.

Metal ions induced secondary structure rearrangements: mechanically interlocked lassovs.unthreaded branched-cyclic topoisomers

Metal ions can play a significant role in a variety of important functions in protein systems including cofactor for catalysis, protein folding, assembly, structural stability and conformational change.

Nickel(ii)-promoted specific hydrolysis of zinc finger proteins

The (S/T)XH sequence in Cys₂His₂zinc fingers can be hydrolytically cleaved by Ni(ii) ions. This reaction can be applied for purification, inhibition or activation of designed zinc finger fusion proteins.

Exploring the boundaries of direct detection and characterization of labile isomers - a case study of copper(ii)-dipeptide systems

The investigation of the linkage isomers of biologically essential and kinetically labile metal complexes in aqueous solutions poses a challenge, as these microspecies cannot be separately studied. Therefore, derivatives are commonly used to initially determine the stability or spectral characteristics of at least one of the isomers. Here we directly detect the isomers, describe the metal ion coordination sphere, speciation and thermodynamic parameters by a synergistic application of temperature dependent EPR and CD spectroscopic measurements in copper(ii)-dipeptide systems including His-Gly and His-Ala ligands. The ΔH = (-23 ± 4) kJ mol^-1 value of the standard enthalpy change corresponding to the peptide-type to histamine-type isomerisation equilibrium of the [CuL]⁺ complex was corroborated by several techniques. The preferential coordination of the side-chains was observed at lower temperatures, whereas, metal-binding of the backbone atoms became favourable upon increasing temperature. This study exemplifies the necessity of using temperature dependent multiple methodologies for a reliable description of similar systems for upstream applications.

The multiple conformational charge states of zinc(II) coordination by 2His-2Cys oligopeptide investigated by ion mobility-mass spectrometry, density functional theory and theoretical collision cross sections

Whether traveling wave ion mobility-mass spectrometry (IM-MS), B3LYP/LanL2DZ density functional theory, and ion size scaled Lennard-Jones (LJ) collision cross sections (CCS) from the B3LYP optimized structures could be used to determine the type of Zn(II) coordination by the oligopeptide acetyl-His₁ -Cys₂ -Gly₃ -Pro₄ -Tyr₅ -His₆ -Cys₇ (amb₅ ) was investigated. The IM-MS analyses of a pH titration of molar equivalents of Zn(II):amb₅ showed that both negatively and positively charged complexes formed and coordination of Zn(II) increased as the His and Cys deprotonated near their pKa values. The B3LYP method was used to generate a series of alternative coordination structures to compare with the experimental results. The method predicted that the single negatively charged complex coordinated Zn(II) in a distorted tetrahedral geometry via the 2His-2Cys substituent groups, whereas, the double negatively charged and positively charged complexes coordinated Zn(II) via His, carbonyl oxygens and the C-terminus. The CCS of the B3LYP complexes were calculated using the LJ method and compared with those measured by IM-MS for the various charge state complexes. The LJ method provided CCS that agreed with five of the alternative distorted tetrahedral and trigonal bipyramidal coordinations for the doubly charged complexes, but provided CCS that were 15 to 31 Ų larger than those measured by IM-MS for the singly charged complexes. Collision-induced dissociation of the Zn(II) complexes and a further pH titration study of amb_5B , which included amidation of the C-terminus, suggested that the 2His-2Cys coordination was more significant than coordinations that included the C-terminus. Copyright © 2016 John Wiley & Sons, Ltd.

The use of ion mobility mass spectrometry to assist protein design: a case study on zinc finger fold versus coiled coil interactions

The dramatic conformational change in zinc fingers on binding metal ions for DNA recognition makes their structure-function behaviour an attractive target to mimic in de novo designed peptides. Mass spectrometry, with its high throughput and low sample consumption provides insight into how primary amino acid sequence can encode stable tertiary fold. We present here the use of ion mobility mass spectrometry (IM-MS) coupled with molecular dynamics (MD) simulations as a rapid analytical platform to inform de novo design efforts for peptide-metal and peptide-peptide interactions. A dual peptide-based synthetic system, ZiCop based on a zinc finger peptide motif, and a coiled coil partner peptide Pp, have been investigated. Titration mass spectrometry determines the relative binding affinities of different divalent metal ions as Zn(2+) > Co(2+) ≫ Ca(2+). With collision induced dissociation (CID), we probe complex stability, and establish that peptide-metal interactions are stronger and more 'specific' than those of peptide-peptide complexes, and the anticipated hetero-dimeric complex is more stable than the two homo-dimers. Collision cross-sections (CCS) measurements by IM-MS reveal increased stability with respect to unfolding of the metal-bound peptide over its apo-form, and further, larger collision cross sections for the hetero-dimeric forms suggest that dimeric species formed in the absence of metal are coiled coil like. MD supports these structural assignments, backed up by data from visible light absorbance measurements.

Analysis of triacetone triperoxide complexes with alkali metal ions by electrospray and extractive electrospray ionisation combined with ion mobility spectrometry and mass spectrometry

The complexation of triacetone triperoxide (TATP) with a range of alkali metals has been studied by electrospray ionisation-mass spectrometry yield [M+Cat](+) ions for all of the alkali metals. The formation of [2TATP+Li+LiX](+) (X = Br, Cl) sandwich complexes was also observed. Collision cross- sections for the lithium-containing complexes of TATP were measured by travelling wave ion mobility spectrometry mass spectrometry, and compared well with computationally determined structures. Extractive electrospray ionisation (EESI) using a lithium doped electrospray is demonstrated for the detection of TATP vapours desorbed from a metal surface. The limit of detection for EESI was shown to be 20 ng using the [TATP+Li](+) ion.

Characterization of the Zn(II) binding properties of the human Wilms' tumor suppressor protein C-terminal zinc finger peptide

Zinc finger proteins that bind Zn(II) using a Cys2His2 coordination motif within a ββα protein fold are the most abundant DNA binding transcription factor domains in eukaryotic systems. These classic zinc fingers are typically unfolded in the apo state and spontaneously fold into their functional ββα folds upon incorporation of Zn(II). These metal-induced protein folding events obscure the free energy cost of protein folding by coupling the protein folding and metal-ion binding thermodynamics. Herein, we determine the formation constant of a Cys2His2/ββα zinc finger domain, the C-terminal finger of the Wilms' tumor suppressor protein (WT1-4), for the purposes of determining its free energy cost of protein folding. Measurements of individual conditional dissociation constants, Kd values, at pH values from 5 to 9 were determined using fluorescence spectroscopy by direct or competition titration. Potentiometric titrations of apo-WT1-4 followed by NMR spectroscopy provided the intrinsic pKa values of the Cys2His2 residues, and corresponding potentiometric titrations of Zn(II)-WT1-4 followed by fluorescence spectroscopy yielded the effective pKa(eff) values of the Cys2His2 ligands bound to Zn(II). The Kd, pKa, and pKa(eff) values were combined in a minimal, complete equilibrium model to yield the pH-independent formation constant value for Zn(II)-WT1-4, Kf(ML) value of 7.5 × 10(12) M(-1), with a limiting Kd value of 133 fM. This shows that Zn(II) binding to the Cys2His2 site in WT1-4 provides at least -17.6 kcal/mol in driving force to fold the protein scaffold. A comparison of the conditional dissociation constants of Zn(II)-WT1-4 to those from the model peptide Zn(II)-GGG-Cys2His2 over the pH range 5.0 to 9.0 and a comparison of their pH-independent Kf(ML) values demonstrates that the free energy cost of protein folding in WT1-4 is less than +2.1 kcal/mol. These results validate our GGG model system for determining the cost of protein folding in natural zinc finger proteins and support the conclusion that the cost of protein folding in most zinc finger proteins is ≤+4.2 kcal/mol, a value that pales in comparison to the free energy contribution of Zn(II) binding, -17.6 kcal/mol.

A database of alkaline-earth-coordinated peptide cross sections: insight into general aspects of structure

A database of 1470 collision cross sections (666 doubly- and 804 triply-charged) of alkaline-earth-coordinated tryptic peptide ions [where the cation (M(2+)) correspond to Mg(2+), Ca(2+), or Ba(2+)] is presented. The utility of such an extensive set of measurements is illustrated by extraction of general properties of M(2+)-coordinated peptide structures. Specifically, we derive sets of intrinsic size parameters (ISPs) for individual amino acid residues for M(2+)-coordinated peptides. Comparison of these parameters with existing ISPs for protonated peptides suggests that M(2+) binding occurs primarily through interactions with specific polar aliphatic residues (Asp, Ser, and Thr) and the peptide backbone. A comparison of binding interactions for these alkaline-earth metals with interactions reported previously for alkali metals is provided. Finally, we describe a new analysis in which ISPs are used as probes for assessing peptide structure based on amino acid composition.

Small-molecule inhibition of c-MYC:MAX leucine zipper formation is revealed by ion mobility mass spectrometry

The leucine zipper interaction between MAX and c-MYC has been studied using mass spectrometry and drift time ion mobility mass spectrometry (DT IM-MS) in addition to circular dichroism spectroscopy. Peptides comprising the leucine zipper sequence with (c-MYC-Zip residues 402-434) and without a postulated small-molecule binding region (c-MYC-ZipΔDT residues 406-434) have been synthesized, along with the corresponding MAX leucine zipper (MAX-Zip residues 74-102). c-MYC-Zip:MAX-Zip complexes are observed both in the absence and in the presence of the reported small-molecule inhibitor 10058-F4 for both forms of c-MYC-Zip. DT IM-MS, in combination with molecular dynamics (MD), shows that the c-MYC-Zip:MAX-Zip complex [M+5H](5+) exists in two conformations, one extended with a collision cross section (CCS) of 1164 ± 9.3 Å(2) and one compact with a CCS of 982 ± 6.6 Å(2); similar values are observed for the two forms of c-MYC-ZipΔDT:MAX-Zip. Candidate geometries for the complexes have been evaluated with MD simulations. The helical leucine zipper structure previously determined from NMR measurements (Lavigne, P.; et al. J. Mol. Biol. 1998, 281, 165), altered to include the DT region and subjected to a gas-phase minimization, yields a CCS of 1247 Å(2), which agrees with the extended conformation we observe experimentally. More extensive MD simulations provide compact complexes which are found to be highly disordered, with CCSs that correspond to the compact form from experiment. In the presence of the ligand, the leucine zipper conformation is completely inhibited and only the more disordered species is observed, providing a novel method to study the effect of interactions of disordered systems and subsequent inhibition of the formation of an ordered helical complex.

Ion mobility mass spectrometry for peptide analysis

The use of ion mobility mass spectrometry has grown rapidly over the last two decades. This powerful analytical platform now forms an attractive prospect for comprehensive analysis of many different molecular species, including chemically complex biological molecules. This paper describes the application of IM-MS to the study of peptides. We focus on three different ion mobility devices that are most frequently found in tandem with mass spectrometers. These are instruments using linear drift tubes (LDT), those using travelling wave ion guides (TWIGS) and those employing high field asymmetric ion mobility spectrometry (FAIMS). Each technique is described. Examples are given on the use of IM-MS for the determination of peptide structure, the study of peptides that form amyloid fibrils, and the study of complex peptide mixtures in proteomic investigations. We describe and comment on the methodologies used and the outlook for this developing analytical technique.