Calcium-Induced Structural Transitions of the Calmodulin−Melittin System Studied by Electrospray Mass Spectrometry: Conformational Subpopulations and Metal-Unsaturated Intermediates

A Conformation-Specific Approach to Native Top-down Mass Spectrometry

Native top-down mass spectrometry is a powerful approach for characterizing proteoforms and has recently been applied to provide similarly powerful insights into protein conformation. Current approaches, however, are limited such that structural insights can only be obtained for the entire conformational landscape in bulk or without any direct conformational measurement. We report a new ion-mobility-enabled method for performing native top-down MS in a conformation-specific manner. Our approach identified conformation-linked differences in backbone dissociation for the model protein calmodulin, which simultaneously informs upon proteoform variations and provides structural insights. We also illustrate that our method can be applied to protein-ligand complexes, either to identify components or to probe ligand-induced structural changes.

Inter-Domain Repulsion of Dumbbell-Shaped Calmodulin during Electrospray Ionization Revealed by Molecular Dynamics Simulations

The mechanisms whereby protein ions are released from nanodroplets at the liquid-gas interface have continued to be controversial since electrospray ionization (ESI) mass spectrometry was widely applied in biomolecular structure analysis in solution. Several viable pathways have been proposed and verified for single-domain proteins. However, the ESI mechanism of multi-domain proteins with more complicated and flexible structures remains unclear. Herein, dumbbell-shaped calmodulin was chosen as a multi-domain protein model to perform molecular dynamics simulations to investigate the structural evolution during the ESI process. For [Ca₄CAM], the protein followed the classical charge residue model. As the inter-domain electrostatic repulsion increased, the droplet was found to split into two sub-droplets, while stronger-repulsive apo-calmodulin unfolded during the early evaporation stage. We designated this novel ESI mechanism as the domain repulsion model, which provides new mechanistic insights into further exploration of proteins containing more domains. Our results suggest that greater attention should be paid to the effect of domain-domain interactions on structure retention during liquid-gas interface transfer when mass spectrometry is used as the developing technique in gas phase structural biology.

Calmodulin's Interdomain Linker Is Optimized for Dynamics Signal Transmission and Calcium Binding

Linkers are ubiquitous in multidomain proteins. These linkers are integral to protein functions, and accumulating evidence suggests that the linkers' versatile roles are encoded in their sequences. However, a molecular picture of how amino acid differences in the linker influence protein function is still lacking. By using extensive Gaussian-accelerated MD coupled with dynamic network analysis, we reveal the molecular bases underlying the linker's role in Calmodulin (CaM), a highly conserved Ca²⁺-signaling hub in eukaryotes. Three CaM constructs comprising a wild-type linker, a flexible linker (four glycines at position D78-S81), and a rigid linker (four prolines at position D78-S81) were simulated. We show that the flexible linker resembles the wild type in allowing CaM to sample a large ensemble of conformations while the rigid linker confines the sampling. Our simulations recapture experimental observations that target binding enhances the Ca²⁺ affinity to CaM's EF-hand sites at the N-domain. However, only the wild-type linker can both correctly capture the Ca²⁺ binding order and maintain the α-helical structure of the domain. The other two constructs either bind Ca²⁺ in an incorrect order or exhibit unfolding of an N-domain helix. We demonstrate that the wild-type linker achieves these outcomes by transmitting interdomain dynamics efficiently. This was evidenced by stronger (anti)correlations among the linker residues, decoupling of the hydrogen bonds between A1-A15 and V35-E45, and structuring of the N-domain for Ca²⁺ binding. This decoupling was not evident for the other two constructs. Lastly, we show that the wild-type linker's optimal transmission stems from its thermodynamically favorable strain and solvation relative to the other two constructs. Our results show how the linker sequence tunes CaM function, suggesting possible mechanisms for changes in linker properties such as mutations or post-translational modifications to modulate protein/substrate binding.

Subsecond Time-Resolved Mass Spectrometry in Dynamic Structural Biology

Life at the molecular level is a dynamic world, where the key players-proteins, oligonucleotides, lipids, and carbohydrates-are in a perpetual state of structural flux, shifting rapidly between local minima on their conformational free energy landscapes. The techniques of classical structural biology, X-ray crystallography, structural NMR, and cryo-electron microscopy (cryo-EM), while capable of extraordinary structural resolution, are innately ill-suited to characterize biomolecules in their dynamically active states. Subsecond time-resolved mass spectrometry (MS) provides a unique window into the dynamic world of biological macromolecules, offering the capacity to directly monitor biochemical processes and conformational shifts with a structural dimension provided by the electrospray charge-state distribution, ion mobility, covalent labeling, or hydrogen-deuterium exchange. Over the past two decades, this suite of techniques has provided important insights into the inherently dynamic processes that drive function and pathogenesis in biological macromolecules, including (mis)folding, complexation, aggregation, ligand binding, and enzyme catalysis, among others. This Review provides a comprehensive account of subsecond time-resolved MS and the advances it has enabled in dynamic structural biology, with an emphasis on insights into the dynamic drivers of protein function.

Millisecond Time-Resolved Solid-State NMR Reveals a Two-Stage Molecular Mechanism for Formation of Complexes between Calmodulin and a Target Peptide from Myosin Light Chain Kinase

Calmodulin (CaM) mediates a wide range of biological responses to changes in intracellular Ca²⁺ concentrations through its calcium-dependent binding affinities to numerous target proteins. Binding of two Ca²⁺ ions to each of the two four-helix-bundle domains of CaM results in major conformational changes that create a potential binding site for the CaM binding domain of a target protein, which also undergoes major conformational changes to form the complex with CaM. Details of the molecular mechanism of complex formation are not well established, despite numerous structural, spectroscopic, thermodynamic, and kinetic studies. Here, we report a study of the process by which the 26-residue peptide M13, which represents the CaM binding domain of skeletal muscle myosin light chain kinase, forms a complex with CaM in the presence of excess Ca²⁺ on the millisecond time scale. Our experiments use a combination of selective ¹³C labeling of CaM and M13, rapid mixing of CaM solutions with M13/Ca²⁺ solutions, rapid freeze-quenching of the mixed solutions, and low-temperature solid state nuclear magnetic resonance (ssNMR) enhanced by dynamic nuclear polarization. From measurements of the dependence of 2D ¹³C-¹³C ssNMR spectra on the time between mixing and freezing, we find that the N-terminal portion of M13 converts from a conformationally disordered state to an α-helix and develops contacts with the C-terminal domain of CaM in about 2 ms. The C-terminal portion of M13 becomes α-helical and develops contacts with the N-terminal domain of CaM more slowly, in about 8 ms. The level of structural order in the CaM/M13/Ca²⁺ complexes, indicated by ¹³C ssNMR line widths, continues to increase beyond 27 ms.

Probing Changes in Ca2+-Induced Interaction Forces between Calmodulin and Melittin by Atomic Force Microscopy

Mechanobiology studies the means by which physical forces and mechanical properties change intra- or inter- biological macromolecules. Calmodulin (CaM) is involved in physiological activities and various metabolic processes in eukaryotic cells. Although the configuration changes in the interaction between calmodulin and melittin have been studied, the biomechanical relationship of their interaction has rarely been explored. Here, we measured the adhesion forces between calmodulin and melittin in solutions of gradient concentration of calcium ions using atomic force microscopy (AFM). We found that the specific (Fi) and nonspecific (F0) adhesion forces between single melittin and calmodulin in a PBS solution were 69.4 ± 5.0 and 29.3 ± 8.9 pN, respectively. In the presence of 10-7 to 10-3 M Ca2+ PBS solution, the Fi increased significantly to 93.8 ± 5.0, 139.9 ± 9.0, 140.4 ± 9.7, 171.5 ± 9.0, and 213.3 ± 17.8 pN, indicating that the unbinding force between melittin and calmodulin increased in the presence of Ca2+ in a concentration-dependent manner. These findings demonstrated that biomechanical studies based on AFM could help us better understand the melittin/calmodulin-binding processes in the presence of calcium and help us design and screen peptide drugs based on calmodulin.

Charge-dependent modulation of specific and nonspecific protein-metal ion interactions in nanoelectrospray ionization mass spectrometry

RATIONALE

Previous studies found that charge state could affect both specific and nonspecific binding of protein-metal ion interactions in nanoelectrospray ionization mass spectrometry (nESI-MS). However, the two kinds of interactions have been studied individually in spite of the problem that they often coexist in the same system. Thus, it is necessary to study the effects of charge state on specific and nonspecific protein-metal ion interactions in one system to reveal more accurate binding state.

METHODS

The HIV-1 nucleocapsid protein (NCp7(31-55)) which can bind specifically and nonspecifically to Zn²⁺ served as the model to show the charge-dependent protein-metal ion interactions. Hydrogen/deuterium exchange (HDX) and photodissociation (PD) were used to demonstrate that specific binding state was correlated with protein structure. In addition to NCp7(31-55), three other model proteins were used to investigate the reason for the charge-dependent nonspecific binding.

RESULTS

For specific binding, we proposed that protein ions with different charge states had different conformations. The HDX results showed that labile protons in the NCp7(31-55)-Zn complex were exchanged in a charge-state-dependent way. The PD experiments revealed differential fragment yields for different charge states. For nonspecific binding, higher charge states had more Zn²⁺ additions, but less SO₄ ^2- additions. The effects of charge states on nonspecific binding levels were entirely the opposite for Zn²⁺ and SO₄ ^2- . These results could reveal that the nonspecific binding was caused by electrostatic interaction.

CONCLUSIONS

For specific binding, NCp7(31-55) with lower charge states have folding and undenatured structures. The binding states of lower charge states can better reflect more native binding states. For nonspecific binding, when multiple metal ions adduct to proteins, the proteins have more net positive charges, which tend to generate higher charge ions during electrospray.

Combining the biased and unbiased sampling strategy into one convenient free energy calculation method

Constructing a free energy landscape for a large molecule is difficult. One has to use either a high temperature or a strong driving force to enhance the sampling on the free energy barriers. In this work, we propose a mixed method that combines these two kinds of acceleration strategies into one simulation. First, it applies an adaptive biasing potential to some replicas of the molecule. These replicas are particularly accelerated in a collective variable space. Second, it places some unbiased and exchangeable replicas at various temperature levels. These replicas generate unbiased sampling data in the canonical ensemble. To improve the sampling efficiency, biased replicas transfer their state variables to the unbiased replicas after equilibrium by Monte Carlo trial moves. In comparison to previous integrated methods, it is more convenient for users. It does not need an initial reference biasing potential to guide the sampling of the molecule. And it is also unnecessary to insert many replicas for the requirement of passing the free energy barriers. The free energy calculation is accomplished in a single stage. It samples the data as fast as a biased simulation and it processes the data as simple as an unbiased simulation. The method provides a minimalist approach to the construction of the free energy landscape. © 2019 Wiley Periodicals, Inc.

Reliable Tracking In-Solution Protein Unfolding via Ultrafast Thermal Unfolding/Ion Mobility-Mass Spectrometry

Sequential unfolding of monomeric proteins is important for the global understanding of local conformational elements (e.g., secondary structures and domain connections) within those protein assemblies. Ion mobility-mass spectrometry (IM-MS) is an emerging and promising technique for probing gradual protein structural perturbations in the gas phase. However, it is still challenging to track sequential unfolding in the solution phase. Here, we extended IM-MS to track in-solution sequential unfolding of monomeric proteins having single and/or multidomains. The present method combines ultrafast local heating effect (LHE)-driven sequential unfolding with IM-MS identification. Protein sequential unfolding in solution is demonstrated by the rapid and controllable IM-MS data switch between native and gradually unfolded states. Our results show that LHE induces gradual protein conformational transitions associated with biological functions, where IM-MS tracks the sequential unfolding of monomeric proteins.

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.

A Database of Transition-Metal-Coordinated Peptide Cross-Sections: Selective Interaction with Specific Amino Acid Residues

Ion mobility mass spectrometry (IMS-MS) techniques were used to generate a database of 2288 collision cross sections of transition-metal-coordinated tryptic peptide ions. This database consists of cross sections for 1253 [Pep + X]²⁺ and 1035 [Pep + X + H]³⁺, where X²⁺ corresponds to Mn²⁺, Co²⁺, Ni²⁺, Cu²⁺, or Zn²⁺. This number of measurements enables the extraction of structural trends for transition-metal-coordinated peptide ions. The range of structures and changes in collision cross sections for X²⁺-coordinated species (compared with protonated species of the same charge state) is similar to Mg²⁺-coordinated species. This suggests that the structures are largely determined by similarities in cation size with differences among the cross section distributions presumably caused by X²⁺ interactions with specific functional groups offered by the residue R-groups or the peptide backbone. Cross section contributions for individual residues upon X²⁺ solvation are assessed with the derivation of intrinsic size parameters (ISPs). The comparison of the [Pep + X]²⁺ ISPs with those previously reported for [Pep + Mg]²⁺ ions displays a lower contribution to the cross section for His, carboxyamidomethylated Cys, and Met, and is consistent with specific metal-residue interactions identified within protein X-ray crystallography databases. Graphical Abstract ᅟ.

Calculation of the Local Free Energy Landscape in the Restricted Region by the Modified Tomographic Method

The free energy landscape is the most important information in the study of the reaction mechanisms of the molecules. However, it is difficult to calculate. In a large collective variable space, a molecule must take a long time to obtain the sufficient sampling during the simulation. To save the calculation quantity, decreasing the sampling region and constructing the local free energy landscape is required in practice. However, the restricted region in the collective variable space may have an irregular shape. Simply restricting one or more collective variables of the molecule cannot satisfy the requirement. In this paper, we propose a modified tomographic method to perform the simulation. First, it divides the restricted region by some hyperplanes and connects the centers of hyperplanes together by a curve. Second, it forces the molecule to sample on the curve and the hyperplanes in the simulation and calculates the free energy data on them. Finally, all the free energy data are combined together to form the local free energy landscape. Without consideration of the area outside the restricted region, this free energy calculation can be more efficient. By this method, one can further optimize the path quickly in the collective variable space.

Structural analysis of calmodulin binding by nNOS inhibitory amphibian peptides

Calmodulin (CaM) is a ubiquitous protein in nature and plays a regulatory role in numerous biological processes, including the upregulation of nitric oxide (NO) synthesis in vivo. Several peptides that prevent NO production by interacting with CaM have been isolated in the cutaneous secretions of Australian amphibians, and are thought to serve as a defense mechanism against predators. In this work, we probe the mechanism by which three of these peptides, namely, caerin 1.8, dahlein 5.6, and a synthetic modification of citropin 1.1, interact with CaM to inhibit NO signaling. Isothermal titration calorimetry was used to determine thermodynamic parameters of the binding interactions and revealed that all the peptides bind to CaM in a similar fashion, with the peptide encapsulated between the two lobes of CaM. Ion mobility-mass spectrometry was used to investigate the changes in collision cross section that occur as a result of complexation, providing additional evidence for this binding mode. Finally, nuclear magnetic resonance spectroscopy was used to track chemical shift changes upon binding. The results obtained confirm that these complexes adopt canonical collapsed structures and demonstrate the strength of the interaction between the peptides and CaM. An understanding of these molecular recognition events provides insights into the underlying mechanism of the amphibian host-defense system.

Effects of metal ion adduction on the gas-phase conformations of protein ions

Changes in protein ion conformation as a result of nonspecific adduction of metal ions to the protein during electrospray ionization (ESI) from aqueous solutions were investigated using traveling wave ion mobility spectrometry (TWIMS). For all proteins examined, protein cations (and in most cases anions) with nonspecific metal ion adducts are more compact than the fully protonated (or deprotonated) ions with the same charge state. Compaction of protein cations upon nonspecific metal ion binding is most significant for intermediate charge state ions, and there is a greater reduction in collisional cross section with increasing number of metal ion adducts and increasing ion valency, consistent with an electrostatic interaction between the ions and the protein. Protein cations with the greatest number of adducted metal ions are no more compact than the lowest protonated ions formed from aqueous solutions. These results show that smaller collisional cross sections for metal-attached protein ions are not a good indicator of a specific metal-protein interaction in solution because nonspecific metal ion adduction also results in smaller gaseous protein cation cross sections. In contrast, the collisional cross section of α-lactalbumin, which specifically binds one Ca(2+), is larger for the holo-form compared with the apo-form, in agreement with solution-phase measurements. Because compaction of protein cations occurs when metal ion adduction is nonspecific, elongation of a protein cation may be a more reliable indicator that a specific metal ion-protein interaction occurs in solution.

A tale of a tail: Structural insights into the conformational properties of the polyglutamine protein ataxin-3

Ataxin-3 is the protein responsible for the neurodegenerative polyglutamine disease Spinocerebellar ataxia type 3. Full structural characterisation of ataxin-3 is required to aid in understanding the mechanism of disease. Despite extensive study, little is known about the conformational properties of the full-length protein, in either its non-expanded healthy or expanded pathogenic forms, particularly since its polyglutamine-containing region has denied structural elucidation. In this work, travelling-wave ion mobility spectrometry-mass spectrometry and limited proteolysis have been used to compare the conformational properties of full-length non-expanded ataxin-3 (14Q) and its isolated N-terminal Josephin domain (JD). Limited proteolysis experiments have confirmed that the JD is stable, being extremely resistant to trypsin digestion, with the exception of the α2/α3 hairpin which is flexible and exposed to protease cleavage in solution. The C-terminal region of ataxin-3 which contains the glutamine-rich sequences is largely unstructured, showing little resistance to limited proteolysis. Using ion mobility spectrometry-mass spectrometry we show that ataxin-3 (14Q) adopts a wide range of conformational states in vitro conferred by the flexibility of its C-terminal tail and the α2/α3 hairpin of the N-terminal JD. This study highlights how the power of MS-based approaches to protein structural characterisation can be particularly useful when the target protein is aggregation-prone and has intrinsically unordered regions.

Binding free-energy calculation of an ion-peptide complex by constrained dynamics

Binding free energy is the most important physical parameter that describes the binding affinity of a receptor-ligand complex. Conventionally, it was obtained based on the thermodynamic cycle or alchemical reaction. These strategies have been widely used, but they would be problematic if the receptors and/or ligands have large conformational changes during the binding processes. In this paper, we present a way to calculate the binding free energy: constrained dynamics along a fragmental and high-dimensional transition path. This method directly considers unbound states in the simulation. The application to the calmodulin loop-calcium complexes shows that it is practical and the calculated relative binding affinities are in good agreement with experimental results.

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.

Dissociation kinetics of the streptavidin-biotin interaction measured using direct electrospray ionization mass spectrometry analysis

Dissociation rate constants (k (off)) for the model high affinity interaction between biotin (B) and the homotetramer of natural core streptavidin (S(4)) were measured at pH 7 and temperatures ranging from 15 to 45 °C using electrospray ionization mass spectrometry (ESI-MS). Two different approaches to data analysis were employed, one based on the initial rate of dissociation of the (S(4) + 4B) complex, the other involving nonlinear fitting of the time-dependent relative abundances of the (S(4) + iB) species. The two methods were found to yield k (off) values that are in good agreement, within a factor of two. The Arrhenius parameters for the dissociation of the biotin-streptavidin interaction in solution were established from the k (off) values determined by ESI-MS and compared with values measured using a radiolabeled biotin assay. Importantly, the dissociation activation energies determined by ESI-MS agree, within 1 kcal mol(-1), with the reported value. In addition to providing a quantitative measure of k (off), the results of the ESI-MS measurements revealed that the apparent cooperative distribution of (S(4) + iB) species observed at short reaction times is of kinetic origin and that sequential binding of B to S(4) occurs in a noncooperative fashion with the four ligand binding sites being kinetically and thermodynamically equivalent and independent.

Distinct mechanisms of calmodulin binding and regulation of adenylyl cyclases 1 and 8

Calmodulin (CaM), by mediating the stimulation of the activity of two adenylyl cyclases (ACs), plays a key role in integrating the cAMP and Ca(2+) signaling systems. These ACs, AC1 and AC8, by decoding discrete Ca(2+) signals can contribute to fine-tuning intracellular cAMP dynamics, particularly in neurons where they predominate. CaM comprises an α-helical linker separating two globular regions at the N-terminus and the C-terminus that each bind two Ca(2+) ions. These two lobes have differing affinities for Ca(2+), and they can interact with target proteins independently. This study explores previous indications that the two lobes of CaM can regulate AC1 and AC8 differently and thereby yield different responses to cellular transitions in [Ca(2+)](i). We first compared by glutathione S-transferase pull-down assays and offline nanoelectrospray ionization mass spectrometry the interaction of CaM and Ca(2+)-binding deficient mutants of CaM with the internal CaM binding domain (CaMBD) of AC1 and the two terminal CaMBDs of AC8. We then examined the influence of these three CaMBDs on Ca(2+) binding by native and mutated CaM in stopped-flow experiments to quantify their interactions. The three CaMBDs show quite distinct interactions with the two lobes of CaM. These findings establish the critical kinetic differences between the mechanisms of Ca(2+)-CaM activation of AC1 and AC8, which may underpin their different physiological roles.

Protein flexibility is key to cisplatin crosslinking in calmodulin

Chemical crosslinking in combination with Fourier transform ion cyclotron resonance mass spectrometry (FTICR MS) has significant potential for studying protein structures and protein-protein interactions. Previously, cisplatin has been shown to be a crosslinker and crosslinks multiple methionine (Met) residues in apo-calmodulin (apo-CaM). However, the inter-residue distances obtained from nuclear magnetic resonance structures are inconsistent with the measured distance constraints by crosslinking. Met residues lie too far apart to be crosslinked by cisplatin. Here, by combining FTICR MS with a novel computational flexibility analysis, the flexible nature of the CaM structure is found to be key to cisplatin crosslinking in CaM. It is found that the side chains of Met residues can be brought together by flexible motions in both apo-CaM and calcium-bound CaM (Ca₄-CaM). The possibility of cisplatin crosslinking Ca₄-CaM is then confirmed by MS data. Therefore, flexibility analysis as a fast and low-cost computational method can be a useful tool for predicting crosslinking pairs in protein crosslinking analysis and facilitating MS data analysis. Finally, flexibility analysis also indicates that the crosslinking of platinum to pairs of Met residues will effectively close the nonpolar groove and thus will likely interfere with the binding of CaM to its protein targets, as was proved by comparing assays for cisplatin-modified/unmodified CaM binding to melittin. Collectively, these results suggest that cisplatin crosslinking of apo-CaM or Ca₄-CaM can inhibit the ability of CaM to recognize its target proteins, which may have important implications for understanding the mechanism of tumor resistance to platinum anticancer drugs.

Characterisation of Calmodulin Structural Transitions by Ion Mobility Mass Spectrometry

This study demonstrates the ability of travelling wave ion mobility-mass spectrometry to measure collision cross-sections of ions in the negative mode, using a calibration based approach. Here, negative mode ion mobility-mass spectrometry was utilised to understand structural transitions of calmodulin upon Ca2+ binding and complexation with model peptides melittin and the plasma membrane Ca2+ pump C20W peptide. Coexisting calmodulin conformers were distinguished on the basis of their mass and cross-section, and identified as relatively folded and unfolded populations, with good agreement in collision cross-section to known calmodulin geometries. Titration of calcium tartrate to physiologically relevant Ca2+ levels provided evidence for intermediately metalated species during the transition from apo- to holo-calmodulin, with collision cross-section measurements indicating that higher Ca2+ occupancy is correlated with more compact structures. The binding of two representative peptides which exemplify canonical compact (melittin) and extended (C20W) peptide-calmodulin binding models has also been interrogated by ion mobility mass spectrometry. Peptide binding to calmodulin involves intermediates with metalation states from 1–4 Ca2+, which demonstrate relatively collapsed structures, suggesting neither the existence of holo-calmodulin or a pre-folded calmodulin conformation is a prerequisite for binding target peptides or proteins. The biological importance of the different metal unsaturated calmodulin complexes, if any, is yet to be understood.

Charge-surface correlation in electrospray ionization of folded and unfolded proteins

Electrospray-ionization mass spectrometry (ESI-MS) is widely used for protein studies. It has been shown that the extent of protein ionization under nondenaturing conditions correlates well with the solvent-accessible surface area of the tridimensional structure, for either folded monomers or multimeric complexes. The goal of this study was to test whether this relation holds for unfolded proteins as well. In order to overcome the paucity of structural data, the server ProtSA was used to model the conformational ensembles of proteins in the unfolded state and generate estimates of the average solvent accessibility. The results are analyzed along with literature data or original measurements by ESI-MS. It is found that the charge-to-surface relation holds for proteins in the unfolded state, free from solvent effects. A double-log plot is derived, in close agreement with published data for folded proteins. These results suggest that the solvent-accessible surface area is a key factor determining the extent of protein ionization by electrospray, independent of the conformational state. This conclusion helps rationalizing conformational effects in protein ESI-MS. The here reported relation can be used to predict the average solvent accessibility and, hence, the state of folding of unknown proteins from ESI-MS data.