Molecular Dynamics and Ion Mobility Spectrometry Study of Model β-Hairpin Peptide, Trpzip1

Resolving Isomers of Star-Branched Poly(Ethylene Glycols) by IMS-MS Using Multiply Charged Ions

Ion mobility spectrometry (IMS) mass spectrometry (MS) centers on the ability to separate gaseous structures by size, charge, shape, and followed by mass-to-charge (m/z). For oligomeric structures, improved separation is hypothesized to be related to the ability to extend structures through repulsive forces between cations electrostatically bonded to the oligomers. Here we show the ability to separate differently branched multiply charged ions of star-branched poly(ethylene glycol) oligomers (up to 2000 Da) regardless of whether formed by electrospray ionization (ESI) charged solution droplets or from charged solid particles produced directly from a surface by matrix-assisted ionization. Detailed structural characterization of isomers of the star-branched compositions was first established using a home-built high-resolution ESI IMS-MS instrument. The doubly charged ions have well-resolved drift times, achieving separation of isomers and also allowing differentiation of star-branched versus linear oligomers. An IMS-MS "snapshot" approach allows visualization of architectural dispersity and (im)purity of samples in a straightforward manner. Analyses capabilities are shown for different cations and ionization methods using commercially available traveling wave IMS-MS instruments. Analyses directly from surfaces using the new ionization processes are, because of the multiply charging, not only associated with the benefits of improved gas-phase separations, relative to that of ions produced by matrix-assisted laser desorption/ionization, but also provide the potential for spatially resolved measurements relative to ESI and other ionization methods.

Conformational investigation of the structure-activity relationship of GdFFD and its analogues on an achatin-like neuropeptide receptor of Aplysia californica involved in the feeding circuit

Proteins and peptides in nature are almost exclusively made from l-amino acids, and this is even more absolute in the metazoan. With the advent of modern bioanalytical techniques, however, previously unappreciated roles for d-amino acids in biological processes have been revealed. Over 30 d-amino acid containing peptides (DAACPs) have been discovered in animals where at least one l-residue has been isomerized to the d-form via an enzyme-catalyzed process. In Aplysia californica, GdFFD and GdYFD (the lower-case letter "d" indicates a d-amino acid residue) modulate the feeding behavior by activating the Aplysia achatin-like neuropeptide receptor (apALNR). However, little is known about how the three-dimensional conformation of DAACPs influences activity at the receptor, and the role that d-residues play in these peptide conformations. Here, we use a combination of computational modeling, drift-tube ion-mobility mass spectrometry, and receptor activation assays to create a simple model that predicts bioactivities for a series of GdFFD analogs. Our results suggest that the active conformations of GdFFD and GdYFD are similar to their lowest energy conformations in solution. Our model helps connect the predicted structures of GdFFD analogs to their activities, and highlights a steric effect on peptide activity at position 1 on the GdFFD receptor apALNR. Overall, these methods allow us to understand ligand-receptor interactions in the absence of high-resolution structural data.

Towards the Analysis of High Molecular Weight Proteins and Protein complexes using TIMS-MS

In the present work, we demonstrate the potential and versatility of TIMS for the analysis of proteins, DNA-protein complexes and protein-protein complexes in their native and denatured states. In addition, we show that accurate CCS measurement are possible and in good agreement with previously reported CCS values using other IMS analyzers (<5% difference). The main challenges for the analysis of high mass proteins and protein complexes in the mobility and m/z domain are described. That is, the analysis of high molecular weight systems in their native state may require the use of higher electric fields or a compromise in the TIMS mobility resolution by reducing the bath gas velocity in order to effectively trap at lower electric fields. This is the first report of CCS measurements of high molecular weight biomolecules and biomolecular complexes (~ 150 kDa) using TIMS-MS.

Conformational Shift of a β-Hairpin Peptide upon Complex Formation with an Oligo-proline Peptide Studied by Mass Spectrometry

So-called super-secondary structures such as the β-hairpin, studied here, form an intermediate hierarchy between secondary and tertiary structures of proteins. Their sequence-derived 'pure' peptide backbone conformation is combined with 'remote' interstrand or interresidue contacts reminiscent of the 3D-structure of full-length proteins. This renders them ideally suited for studying potential nucleation sites of protein folding reactions as well as intermolecular interactions. But β-hairpins do not merely serve as model systems; their unique structure characteristics warrant a central role in structural studies on their own. In this study we applied photo cross-linking in combination with high-resolution mass spectrometry and computational modeling as well as with ion mobility-mass spectrometry to elucidate these structural properties. Using variants of a known β-hairpin representative, the so-called trpzip peptide and its ligands, we found evidence for a conformational transition of the β-hairpin and its impact on ligand binding.

From thermodynamics to kinetics: enhanced sampling of rare events

Despite great advances in molecular dynamics simulations, there remain large gaps between the simulations and experimental observations in terms of the time and length scales that can be approached. Developing fast and accurate algorithms and methods is of ultimate importance to bridge these gaps. In this Account, we briefly summarize recent efforts in such directions. In particular, we focus on integrated tempering sampling. The efficiency of this sampling method has been demonstrated by applications to a range of chemical and biological problems: protein folding, molecular cluster structure searches, and chemical reactions. The combination of integrated tempering sampling and a trajectory sampling method allows the calculation of rate constants and reaction pathways without predefined collective coordinates.

The pH Dependence of Product Ion Spectra Obtained from Precursor Ions with the Same Charge Number in ESI of Carbonic Anhydrase 2

The effect of solvent conditions, including pH, on product ion spectra obtained from precursor ions produced by electrospray ionization (ESI) was examined. Bovine carbonic anhydrase 2 was used as a model protein and the product ions generated by collision induced dissociation of the whole protein were measured under several different solvent conditions (pH 5.0, 3.7, and 0.1% HCOOH (pH 2.6)/MeCN (1/1)). The product ion spectra from precursor ions with the same charge number, the observed m/z values and the relative intensities of the product ions were similar. It therefore appears that the solvent conditions used have no effect on the product ion that is generated. On the other hand, different profiles of the product ion were obtained from precursor ions having different charge numbers. This indicates that the charge number of the precursor ion appears to be a major determinant of the product ion species and its relative intensity in product ion spectra of proteins.

Gas-phase ion isomer analysis reveals the mechanism of peptide sequence scrambling

Peptide sequence scrambling during mass spectrometry-based gas-phase fragmentation analysis causes misidentification of peptides and proteins. Thus, there is a need to develop an efficient approach to probing the gas-phase fragment ion isomers related to sequence scrambling and the underlying fragmentation mechanism, which will facilitate the development of bioinformatics algorithm for proteomics research. Herein, we report on the first use of electron transfer dissociation (ETD)-produced diagnostic fragment ions to probe the components of gas-phase peptide fragment ion isomers. In combination with ion mobility spectrometry (IMS) and formaldehyde labeling, this novel strategy enables qualitative and quantitative analysis of b-type fragment ion isomers. ETD fragmentation produced diagnostic fragment ions indicative of the precursor ion isomer components, and subsequent IMS analysis of b ion isomers provided their quantitative and structural information. The isomer components of three representative b ions (b9, b10, and b33 from three different peptides) were accurately profiled by this method. IMS analysis of the b9 ion isomers exhibited dynamic conversion among these structures. Furthermore, molecular dynamics simulation predicted theoretical drift time values, which were in good agreement with experimentally measured values. Our results strongly support the mechanism of peptide sequence scrambling via b ion cyclization, and provide the first experimental evidence to support that the conversion from molecular precursor ion to cyclic b ion (M → (c)b) pathway is less energetically (or kinetically) favored.

Contribution of the empirical dispersion correction on the conformation of short alanine peptides obtained by gas-phase QM calculations

In this work we analyze the effect of the inclusion of an empirical dispersion term to standard DFT (DFT-D) in the prediction of the conformational energy of the alanine dipeptide (Ala2) and in assessing the relative stabilities of short polyala-nine peptides in helical conformations, i.e., α and 310 helices, from Ala4 to Ala16. The Ala2 conformational energies obtained with the dispersion-corrected GGA functional B97-D are compared to previously published high level MP2 data. Meanwhile, the B97-D performance on larger polyalanine peptides is compared to MP2, B3LYP and RHF calculations obtained at a lower level of theory. Our results show that electron correlation affects the conformational energies of short peptides with a weight that increases with the peptide length. Indeed, while the contribution of vdW forces is significant for larger peptides, in the case of Ala2 it is negligible when compared to solvent effects. Even for short peptides, the inclusion of an empirical dispersion term greatly improves accuracy of DFT methods, providing results that correlate very well with the MP2 reference at no additional computational cost.

IMSPeptider: A computational peptide collision cross-section area calculator based on a novel molecular dynamics simulation protocol

Introduction of ion mobility mass spectrometry (IMS/MS) into the proteomic workflow provides an orthogonal separation to the widely used LC-MS platforms. IMS also provides structural information that could facilitate peptide identification. However, the lack of tools capable of predictive power in a high-throughput fashion makes peptide global profiling quite challenging. To target this issue, a computational workflow was developed based on biophysical principles to predict the collision cross-section area (CCS) of peptides as measured from IMS/MS experiments. Hosted on a web server, it allows the user to input a primary sequence (query) and retrieve information on peptide structure, sequence, and corresponding CCS. The current version is designed to identify peptide sequences up to 23 residues in length, in its higher charge state, based on a match of the molecule m/z and CCS. The protocol was validated against a 128-sequences-dataset and CCS predicted within 2.8% average error.

Do charge state signatures guarantee protein conformations?

The extent to which proteins in the gas phase retain their condensed-phase structure is a hotly debated issue. Closely related to this is the degree to which the observed charge state reflects protein conformation. Evidence from electron capture dissociation, hydrogen/deuterium exchange, ion mobility, and molecular dynamics shows clearly that there is often a strong correlation between the degree of folding and charge state, with the most compact conformations observed for the lowest charge states. In this article, we address recent controversies surrounding the relationship between charge states and folding, focussing also on the manipulation of charge in solution and its effect on conformation. 'Supercharging' reagents that have been used to effect change in charge state can promote unfolding in the electrospray droplet. However for several protein complexes, supercharging does not appear to perturb the structure in that unfolding is not detected. Consequently, a higher charge state does not necessarily imply unfolding. Whilst the effect of charge manipulation on conformation remains controversial, there is strong evidence that a folded, compact state of a protein can survive in the gas phase, at least on a millisecond timescale. The exact nature of the side-chain packing and secondary structural elements in these compact states, however, remains elusive and prompts further research.

Statistical analysis of ion mobility spectrometry. II. Adaptively biased methods and shape correlations

Following a recent effort [J. Am. Soc. Mass Spectrom. 23, 386-396 (2012)], we continue to explore computational methodologies for generating molecular conformations to support collisional cross sections suggested by ion mobility measurements. Here, adaptively biased molecular dynamics (ABMD) simulations are used to sample the configuration space and to achieve flat-histogram sampling along the reaction coordinates of the first two moments of the gyration tensor. The method is tested and compared with replica-exchange simulations on triply-protonated bradykinin and on a larger 25-residue peptide. It is found to have a much higher efficiency for producing large sets of conformations in a broad range of diffusion cross-sections, whereas it does not compete with conventional replica-exchange molecular dynamics in locating the lowest-energy structure. Nevertheless, the broad sampling obtained from the ABMD method allows to quantitatively correlate the diffusion cross-section Ω with other geometric order parameters that have simpler interpretation. The strong correlations found between the diffusion cross-section and the radius of gyration, the surface area and the volume of the convex hull suggest an optimal template for accurately mimicking the variations of Ω in a broad range of conformations, using only geometrical information and doing so at a very moderate computational cost. The existence of such a correlation is confirmed on the much larger protein α-lactalbumin.

Damping factor links periodic focusing and uniform field ion mobility for accurate determination of collision cross sections

The methodology for obtaining accurate ion-neutral collision cross section (Ω) values for peptides and proteins using periodic focusing ion mobility spectrometry (PF IMS) is presented. A mobility dampening factor (represented by the term α) is introduced to account for the relative increase in ion-neutral collisions in PF IMS compared to uniform field ion mobility spectrometry (UF IMS) for equivalent operating conditions. The results show that α may be easily quantified both theoretically and empirically for a specific PF IMS design operating at a given pressure based upon the charge state of the analyte. By simply incorporating an α term into traditional UF IMS expressions, accurate Ω values were obtained with excellent agreement (≤4% difference) compared to UF IMS measurements found in the current literature.

Number of solution states of bradykinin from ion mobility and mass spectrometry measurements

Ion mobility and mass spectrometry measurements have been used to examine the populations of different solution structures of the nonapeptide bradykinin. Over the range of solution compositions studied, from 0:100 to 100:0 methanol:water and 0:100 to 90:10 dioxane:water, evidence for 10 independent populations of bradykinin structures in solution is found. In some solutions as many as eight structures may coexist. The solution populations are substantially different than the gas-phase equilibrium distribution of ions, which exhibits only three distinct states. Such a large number of coexisting structures explains the inability of traditional methods of characterization such as nuclear magnetic resonance spectroscopy and crystallography to determine detailed structural features for some regions of this peptide.