Helices and Sheets in vacuo

High-Performance Molecular Dynamics Simulations for Native Mass Spectrometry of Large Protein Complexes with the Fast Multipole Method

Native mass spectrometry (MS) is widely employed to study the structures and assemblies of proteins ranging from small monomers to megadalton complexes. Molecular dynamics (MD) simulation is a useful complement as it provides the spatial detail that native MS cannot offer. However, MD simulations performed in the gas phase have suffered from rapidly increasing computational costs with the system size. The primary bottleneck is the calculation of electrostatic forces, which are effective over long distances and must be explicitly computed for each atom pair, precluding efficient use of methods traditionally used to accelerate condensed-phase simulations. As a result, MD simulations have been unable to match the capacity of MS in probing large multimeric protein complexes. Here, we apply the fast multipole method (FMM) for computing the electrostatic forces, recently implemented by Kohnke et al. (J. Chem. Theory Comput., 2020, 16, 6938-6949), showing that it significantly enhances the performance of gas-phase simulations of large proteins. We assess how to achieve adequate accuracy and optimal performance with FMM, finding that it expands the accessible size range and time scales dramatically. Additionally, we simulate a 460 kDa ferritin complex over microsecond time scales, alongside complementary ion mobility (IM)-MS experiments, uncovering conformational changes that are not apparent from the IM-MS data alone.

Infrared Multiphoton Dissociation Enables Top-Down Characterization of Membrane Protein Complexes and G Protein-Coupled Receptors

Membrane proteins are challenging to analyze by native mass spectrometry (MS) as their hydrophobic nature typically requires stabilization in detergent micelles that are removed prior to analysis via collisional activation. There is however a practical limit to the amount of energy which can be applied, which often precludes subsequent characterization by top-down MS. To overcome this barrier, we have applied a modified Orbitrap Eclipse Tribrid mass spectrometer coupled to an infrared laser within a high-pressure linear ion trap. We show how tuning the intensity and time of incident photons enables liberation of membrane proteins from detergent micelles. Specifically, we relate the ease of micelle removal to the infrared absorption of detergents in both condensed and gas phases. Top-down MS via infrared multiphoton dissociation (IRMPD), results in good sequence coverage enabling unambiguous identification of membrane proteins and their complexes. By contrasting and comparing the fragmentation patterns of the ammonia channel with two class A GPCRs, we identify successive cleavage of adjacent amino acids within transmembrane domains. Using gas-phase molecular dynamics simulations, we show that areas prone to fragmentation maintain aspects of protein structure at increasing temperatures. Altogether, we propose a rationale to explain why and where in the protein fragment ions are generated.

Infrared Multiphoton Dissociation Enables Top-Down Characterization of Membrane Protein Complexes and G Protein-Coupled Receptors

Membrane proteins are challenging to analyze by native mass spectrometry (MS) as their hydrophobic nature typically requires stabilization in detergent micelles that are removed prior to analysis via collisional activation. There is however a practical limit to the amount of energy which can be applied, which often precludes subsequent characterization by top-down MS. To overcome this barrier, we have applied a modified Orbitrap Eclipse Tribrid mass spectrometer coupled to an infrared laser within a high-pressure linear ion trap. We show how tuning the intensity and time of incident photons enables liberation of membrane proteins from detergent micelles. Specifically, we relate the ease of micelle removal to the infrared absorption of detergents in both condensed and gas phases. Top-down MS via infrared multiphoton dissociation (IRMPD), results in good sequence coverage enabling unambiguous identification of membrane proteins and their complexes. By contrasting and comparing the fragmentation patterns of the ammonia channel with two class A GPCRs, we identify successive cleavage of adjacent amino acids within transmembrane domains. Using gas-phase molecular dynamics simulations, we show that areas prone to fragmentation maintain aspects of protein structure at increasing temperatures. Altogether, we propose a rationale to explain why and where in the protein fragment ions are generated.

D-Amino acid recognition of tripeptides studied by ultraviolet photodissociation spectroscopy of hydrogen-bonded clusters

To understand the roles of D-amino acids, evaluating their chemical properties in living organisms is essential. Herein, D-amino acid recognition of peptides was investigated using a tandem mass spectrometer equipped with an electrospray ionization source and a cold ion trap. Ultraviolet (UV) photodissociation spectroscopy and water adsorption of hydrogen-bonded protonated clusters of tryptophan (Trp) enantiomers and tripeptides (SAA, ASA, and AAS, where S and A denote L-serine and L-alanine, respectively) were carried out at 8 K in the gas phase. In the UV photodissociation spectrum of H⁺(D-Trp)ASA, the bandwidth of the S₁-S₀ transition, which corresponds to the ππ* state of the Trp indole ring, was narrower than those of the other five clusters, H⁺(D-Trp)SAA, H⁺(D-Trp)AAS, H⁺(L-Trp)SAA, H⁺(L-Trp)ASA, and H⁺(L-Trp)AAS. In the UV photoexcitation of H⁺(D-Trp)ASA(H₂O)_n, which were formed via water adsorption on gas-phase H⁺(D-Trp)ASA, the evaporation of water molecules was the main photodissociation pathway. An NH₂CHCOOH-eliminated ion and H⁺ASA were observed in the product ion spectrum. By contrast, water molecules adsorbed on the other five clusters remained on the product ions for NH₂CHCOOH elimination and Trp detachment after the UV photoexcitation. The results indicated that the indole ring of Trp was located on the surface of H⁺(D-Trp)ASA, and the amino and carboxyl groups of Trp formed hydrogen bonds in H⁺(D-Trp)ASA. For the other five clusters, the indole rings of Trp were hydrogen bonded in the clusters, and the amino and carboxyl groups of Trp were present on the cluster surfaces.

Diserinol Isophthalamide: A Novel Reagent for Complexation with Biomolecular Anions in Electrospray Ionization Mass Spectrometry

Transferring biomolecules from solution to vacuum facilitates a detailed analysis of molecular structure and dynamics by isolating molecules of interest from a complex environment. However, inherent in the ion desolvation process is the loss of solvent hydrogen bonding partners, which are critical for the stability of a condensed-phase structure. Thus, transfer of ions to vacuum can favor structural rearrangement, especially near solvent-accessible charge sites, which tend to adopt intramolecular hydrogen bonding motifs in the absence of solvent. Complexation of monoalkylammonium moieties (e.g., lysine side chains) with crown ethers such as 18-crown-6 can disfavor structural rearrangement of protonated sites, but no equivalent ligand has been investigated for deprotonated groups. Herein we describe diserinol isophthalamide (DIP), a novel reagent for the gas-phase complexation of anionic moieties within biomolecules. Complexation is observed to the C-terminus or side chains of the small model peptides GD, GE, GG, DF-OMe, VYV, YGGFL, and EYMPME in electrospray ionization mass spectrometry (ESI-MS) studies. In addition, complexation is observed with the phosphate and carboxylate moieities of phosphoserine and phosphotyrosine. DIP performs favorably in comparison to an existing anion recognition reagent, 1,1'-(1,2-phenylene)bis(3-phenylurea), that exhibits moderate carboxylate binding in organic solvent. This improved performance in ESI-MS experiments is attributed to reduced steric constraints to complexation with carboxylate groups of larger molecules. Overall, diserinol isophthalamide is an effective complexation reagent that can be applied in future work to study retention of solution-phase structure, investigate intrinsic molecular properties, and examine solvation effects.

Chemical properties of inner and surficial regions of hydrogen-bonded clusters of biological molecules: ultraviolet photodissociation and water adsorption analyses in the gas phase

The chemical and physical properties of cold, gas-phase hydrogen-bonded clusters of L-alanine (L-Ala), L-trialanine (L-Ala₃), L-tetraalanine (L-Ala₄), and tryptophan (Trp) enantiomers were investigated using tandem mass spectrometry with an electrospray ionization source and cold ion trap. From the ultraviolet (UV) photodissociation spectra at 265-290 nm, the electronic structures of homochiral H⁺(L-Trp)(L-Ala) at 8 K were found to be different from those of heterochiral H⁺(D-Trp)(L-Ala) and protonated Trp. The number of water molecules adsorbed on the surface of gas-phase H⁺(D-Trp)(L-Ala) was larger than that of H⁺(L-Trp)(L-Ala), indicating stronger intermolecular interactions of L-Ala with H⁺(L-Trp) than those with H⁺(D-Trp). The product ion spectrum obtained by 265 nm photoexcitation of H⁺(L-Trp)(L-Ala₃)(H₂O)_n formed via gas-phase water adsorption on H⁺(L-Trp)(L-Ala₃) showed that the evaporation of water molecules was the main photodissociation process. In the case of H⁺(L-Trp)(L-Ala₄)(H₂O)_n, signals of H⁺(L-Ala₄) (H₂O)_n formed via L-Trp evaporation were observed in the product ion spectra, and the cross-section for UV photoinduced L-Trp evaporation became larger as the number of adsorbed water molecules increased. This observation indicates that water molecules were selectively adsorbed on the H⁺(L-Ala₄) side of H⁺(L-Trp)(L-Ala₄) and weakened the intermolecular interactions between L-Trp and H⁺(L-Ala₄) in the hydrogen-bonded cluster.

Forming a Double-Helix Phase of Single Polymer Chains by the Cooperation between Local Structure and Nonlocal Attraction

Double-helix structures, such as DNA, are formed in nature to realize many unique functions. Inspired by this, researchers are pursuing strategies to design such structures from polymers. A key question is whether the double helix can be formed from the self-folding of a single polymer chain without specific interactions. Here, using Langevin dynamics simulation and theoretical analysis, we find that a stable double-helix phase can be achieved by the self-folding of single semiflexible polymers as a result of the cooperation between local structure and nonlocal attraction. The critical temperature of double-helix formation approximately follows T^{cri}∼ln(k_{θ}) and T^{cri}∼ln(k_{τ}), where k_{θ} and k_{τ} are the polymer bending and torsion stiffness, respectively. Furthermore, the double helix can exhibit major and minor grooves due to symmetric break for better packing. Our results provide a novel guide to the experimental design of the double helix.

Viral Prefusion Targeting Using Entry Inhibitor Peptides: The Case of SARS-CoV-2 and Influenza A virus

In this study, peptide entry inhibitors against the fusion processes of severe acute respiratory syndrome coronavirus-2 (SCV2) and influenza A virus (IAV) were designed and evaluated. Fusion inhibitor peptides targeting the conformational shift of the viral fusion protein were designed based on the relatively conserved sequence of HR2 from SCV2 spike protein and the conserved fusion peptide from hemagglutinin (HA) of IAV. Helical HR2 peptides bind more efficiently to HR1 trimer, while helical amphipathic anti-IAV peptides have higher cell penetration and endosomal uptake. The initial sequences were mutated by increasing the amphipathicity, using helix favoring residues, and residues likely to form salt- and disulfide-bridges. After docking against their targets, all anti-SCV2 designed peptides bonded with the HR1 3-helical bundle's hydrophobic crevice, while AntiSCV2P1, AntiSCV2P3, AntiSCV2P7, and AntiSCV2P8 expected to form coiled coils with at least one of the HR1 strands. Four of the designed anti-IAV peptides were cell-penetrating (AntiIAVP2, AntiIAVP3, AntiIAVP4, AntiIAVP7). All of them interacted with the fusion peptide of HA and some of the residues in the conserved hydrophobic pocket of HA2 in H1N1, H3N1, and H5N1 subtypes of IAV. AntiIAVP3 and AntiIAVP4 peptides had the best binding to HA2 conserved hydrophobic pocket, while, AntiIAVP2 and AntiIAVP6 showed the best binding to the fusion peptide region. According to analyses for in-vivo administration, AntiSCV2P1, AntiSCV2P7, AntiIAVP2, and AntiIAVP7 were the best candidates. AntiSCV2 and AntiIAV peptides were also conjugated using an in vivo cleavable linker sensitive to TMPRSS2 applicable as a single therapeutic in coinfections or uncertain diagnosis.

Modulating the Stability of Collagen Triple Helices by Terminal Charged Residues

Cationic or anionic residues are frequently located at the termini of proteins because their charged side chain can form electrostatic interactions with a terminal carboxylate or ammonium group to stabilize the structure under physiological conditions. Here, we used collagen-mimetic peptides (CMPs) to examine how the terminal charge-charge interactions affect the collagen triple helix stability. We designed a series of CMPs with either a Lys or Glu incorporated into the terminus and measured their pH-dependent stability. The results showed that the terminal electrostatic attractions stabilized the triple helix, while the terminal electrostatic repulsions destabilized the trimer. The data also revealed that the repulsions imposed a greater effect than did the attractions on the triple helix. An amino acid with a shorter side chain, such as aspartate and ornithine, was also installed to investigate the length effect on electrostatic interactions, which was found to be insignificant. Meanwhile, simultaneously incorporating cationic and anionic residues into the termini showed slight additive stabilization effects but pronounced additive destabilization consequences. We have demonstrated that the collagen triple helix stability can be modulated by introducing a cationic or anionic residue into the terminus of a peptide, giving useful information for the design of collagen-associated materials.

Investigations into pesticide charge site isomers using conventional IM and cIM systems

A growing number of pesticides are being used around the world necessitating strict regulatory policies to guarantee consumer safety. Liquid Chromatography - Mass Spectrometry (LC-MS) is a highly sensitive method for pesticide screening, which provides retention time, mass/charge ratios and the relative abundances of characteristic product ions. Variability in the latter necessitates relatively large tolerances (±30%, SANCO/12682/2019, current EU regulation). One cause of this variability may stem from the presence of different charge-site isomers (charge carrier being a proton, sodium cation, potassium cation and alike); each yielding a set of different product ions, of which the relative ratios are influenced by solution and ion source conditions. Consequently, varying relative abundances may be observed for analyte ions produced from calibration standards, chemical residues in food matrices and across different instruments. Ion Mobility Spectrometry (IMS) is a fast, gas phase separation technique which can resolve charge-site isomers based on differences in their collisional cross sections (CCSs). We previously used the IM device embedded in LC-IM-MS geometry to generate a pesticide CCS database and subsequently focussed upon identification of pesticides which form charge-site isomers. Latterly, we applied this approach to screen food commodities for pesticide residues. In some instances, isomer separation was clear, however sometimes broad, unresolved distributions were observed. Using a high-resolution cyclic IM device (cIM) we resolved and determined CCS values of species of indoxacarb, spinosad, fenpyroximate, epoxiconazole, metaflumizone and avermectin. Furthermore, utilising novel cIM functionalities (tandem-IM) we discovered that two spinosyn sodimers can interconvert in the gas phase.

Sequence-Specific Model for Predicting Peptide Collision Cross Section Values in Proteomic Ion Mobility Spectrometry

The contribution of peptide amino acid sequence to collision cross section values (CCS) has been investigated using a dataset of ∼134 000 peptides of four different charge states (1+ to 4+). The migration data were acquired using a two-dimensional liquid chromatography (LC)/trapped ion mobility spectrometry/quadrupole/time-of-flight mass spectrometry (MS) analysis of HeLa cell digests created using seven different proteases and was converted to CCS values. Following the previously reported modeling approaches using intrinsic size parameters (ISP), we extended this methodology to encode the position of individual residues within a peptide sequence. A generalized prediction model was built by dividing the dataset into eight groups (four charges for both tryptic/nontryptic peptides). Position-dependent ISPs were independently optimized for the eight subsets of peptides, resulting in prediction accuracy of ∼0.981 for the entire population of peptides. We find that ion mobility is strongly affected by the peptide's ability to solvate the positively charged sites. Internal positioning of polar residues and proline leads to decreased CCS values as they improve charge solvation; conversely, this ability decreases with increasing peptide charge due to electrostatic repulsion. Furthermore, higher helical propensity and peptide hydrophobicity result in a preferential formation of extended structures with higher than predicted CCS values. Finally, acidic/basic residues exhibit position-dependent ISP behavior consistent with electrostatic interaction with the peptide macrodipole, which affects the peptide helicity. The MS raw data files have been deposited with the ProteomeXchange Consortium via the jPOST partner repository (http://jpostdb.org) with the dataset identifiers PXD021440/JPST000959, PXD022800/JPST001017, and PXD026087/ JPST001176.

Exploring Conformational Landscapes Using Trap and Release Tandem Ion Mobility Spectrometry

The dynamics and thermodynamics of structural changes in isolated glu-fibrinopeptide B (GluFib) were investigated by tandem ion mobility spectrometry (IMS). Doubly protonated GluFib²⁺ ions were first selected by IMS and then stored for a controlled duration in a thermalized ion trap. Temperature-induced conformational changes were finally monitored by IMS as a function of trapping time. Based on this procedure, isomerization rates and equilibrium populations of the different conformers were determined as a function of temperature. We demonstrate that the measured thermodynamic quantities can be directly compared to simulated observables from ensemble molecular modeling based on appropriate order parameters. We obtained good qualitative agreement with replica-exchange molecular dynamics simulations based on the AMOEBA force field and processed using the weighted histogram analysis method. This suggests that the balance between Coulomb repulsion and optimal charge solvation is the main source of the observed conformational bistability. Our results emphasize the differences between the kinetically driven quasi-equilibrium distributions obtained after collisional activation and the thermodynamically driven distributions from the present equilibrium experiments due to entropic effects. As a consequence, our measurements not only allow straightforward determination of Arrhenius activation energies but also yield the relative enthalpy and entropy changes associated to a structural transition.

Characterizing Native and Hydrocarbon-Stapled Enfuvirtide Conformations with Ion Mobility Mass Spectrometry and Hydrogen-Deuterium Exchange

The number of approved peptide therapeutics, as well as those in development, has been increasing in recent years. Frequently, the biological activity of such peptides is elicited through the adoption of secondary structural elements upon interaction with their cellular target. However, many therapeutic peptides are unstructured in solution and accordingly exhibit a poor bioavailability due to rapid proteolysis in vivo. To combat this degradation, numerous naturally occurring peptides with therapeutic properties contain stabilizing features, such as N-to-C cyclization or disulfide bonds. Recently, hydrocarbon stapling via non-native amino acid substitution followed by ring-closing metathesis has been shown to induce a dramatic stabilization of α-helical peptides. Identifying the ideal staple location along the peptide backbone is a critical developmental step, and methods to streamline this optimization are needed. Mass spectrometry-based methods such as ion mobility (IM) and hydrogen-deuterium exchange (HDX) can detect multiple discrete peptide conformations, a significant advantage over bulk spectroscopic techniques. In this study we use IM-MS and HDX-MS to demonstrate that the native 36-residue enfuvirtide peptide is highly dynamic in solution and the conformational ensemble populated by stabilized constructs depends heavily on the staple location. Further, our measurements yielded results that correlate well with the average α-helical content measured by circular dichroism. The MS-based approaches described herein represent sensitive and potentially high-throughput methods for characterizing and identifying optimally stapled peptides.

Deep learning the collisional cross sections of the peptide universe from a million experimental values

The size and shape of peptide ions in the gas phase are an under-explored dimension for mass spectrometry-based proteomics. To investigate the nature and utility of the peptide collisional cross section (CCS) space, we measure more than a million data points from whole-proteome digests of five organisms with trapped ion mobility spectrometry (TIMS) and parallel accumulation-serial fragmentation (PASEF). The scale and precision (CV < 1%) of our data is sufficient to train a deep recurrent neural network that accurately predicts CCS values solely based on the peptide sequence. Cross section predictions for the synthetic ProteomeTools peptides validate the model within a 1.4% median relative error (R > 0.99). Hydrophobicity, proportion of prolines and position of histidines are main determinants of the cross sections in addition to sequence-specific interactions. CCS values can now be predicted for any peptide and organism, forming a basis for advanced proteomics workflows that make full use of the additional information.

Proteomics has been advanced by algorithms that can predict different peptide features, but predicting peptide collisional cross sections (CCS) has remained challenging. Here, the authors measure over one million CCS values of tryptic peptides and develop a deep learning model for peptide CCS prediction.

Mechanistic Insights into an Unusual Side-Chain-Mediated N-Cα Bond Cleavage under Collision-Induced Dissociation Conditions in the Disulfide-Containing Peptide Conopressin

Conopressin, a nonapeptide disulfide CFIRNCPKG amide present in cone snail venom, undergoes a facile cleavage at the Cys6-Pro7 peptide bond to yield a disulfide bridged b₆ ion. Analysis of the mass spectral fragmentation pattern reveals the presence of a major fragment ion, which is unambiguously assigned as the tripeptide sequence IRN amide. The sequence dependence of this unusual fragmentation process has been investigated by comparing it with the fragmentation patterns of related peptides, oxytocin (CYIQNCPLG amide), Lys-vasopressin (CYFQNCPKG amide), and a series of synthetic analogues. The results establish the role of the Arg4 residue in facilitating the unusual N-C^α bond cleavage at Cys6. Structures are proposed for a modified disulfide bridged fragment containing the Cys1 and Cys6 residues. Gas-phase molecular dynamics simulations provide evidence for the occurrence of conformational states that permit close approach of the Arg4 side chain to the Cys6 C^β methylene protons.

Secondary structure effects on internal proton transfer in poly-peptides

A pump-probe approach was designed to determine the internal proton transfer (PT) rate in a series of poly-peptide radical cations containing both histidine and tryptophan. The proton transfer is driven by the gas-phase basicity difference between residues. The fragmentation scheme indicates that the gas-phase basicity of histidine is lower than that of radical tryptophan so that histidine is always pulling the proton away from tryptophan. However, the proton transfer requires the two basic sites to be in close proximity, which is rate limited by the peptide conformational dynamics. PT rate measurements were used to probe and explore the peptide conformational dynamics in several poly-glycines/prolines/alanines. For small and unstructured peptides, the PT rate decreases with the size, as expected from a statistical point of view in a flat conformational space. Conversely, if structured conformations are accessible, the structural flexibility of the peptide is decreased. This slows down the occurrence of conformations favorable to proton transfer. A dramatic decrease in the PT rates was observed for peptides HAnW, when n changes from 5 to 6. This is attributed to the onset of a stable helix for n = 6. No such discontinuity is observed for poly-glycines or poly-prolines. In HAnW, the gas-phase basicity and helix propensity compete for the position of the charge. Interestingly, in this competition between PT and helix formation in HA6W, the energy gain associated with helix formation is large enough to slow down the PT beyond experimental time but does not ultimately prevail over the proton preference for histidine.

Automated exploration of the low-energy chemical space with fast quantum chemical methods

We propose and discuss an efficient scheme for the in silico sampling for parts of the molecular chemical space by semiempirical tight-binding methods combined with a meta-dynamics driven search algorithm. The focus of this work is set on the generation of proper thermodynamic ensembles at a quantum chemical level for conformers, but similar procedures for protonation states, tautomerism and non-covalent complex geometries are also discussed. The conformational ensembles consisting of all significantly populated minimum energy structures normally form the basis of further, mostly DFT computational work, such as the calculation of spectra or macroscopic properties. By using basic quantum chemical methods, electronic effects or possible bond breaking/formation are accounted for and a very reasonable initial energetic ranking of the candidate structures is obtained. Due to the huge computational speedup gained by the fast low-cost quantum chemical methods, overall short computation times even for systems with hundreds of atoms (typically drug-sized molecules) are achieved. Furthermore, specialized applications, such as sampling with implicit solvation models or constrained conformational sampling for transition-states, metal-, surface-, or noncovalently bound complexes are discussed, opening many possible applications in modern computational chemistry and drug discovery. The procedures have been implemented in a freely available computer code called CREST, that makes use of the fast and reliable GFNn-xTB methods.

Hydrogen-Deuterium Exchange and Electron Capture Dissociation to Interrogate the Conformation of Gaseous Melittin Ions

There is a need in the field of biological mass spectrometry for structural tools which can report on regional, rather than solely global, structure of gaseous protein ions. Site-specific hydrogen-deuterium (H/D) exchange has shown promise in fulfilling this need, but requires additional method development to prove its utility. In this study, we use H/D exchange and electron capture dissociation (ECD) to probe the gaseous structure of two peptides which are α-helical in solution and which differ by a single point mutation. Global H/D exchange levels, ECD fragmentation profiles, and region specific H/D exchange profiles are compared between wild type (WT) melittin, which adopts a hinged helix conformation in solution, and a mutant P14A melittin which folds into a single helix in solution. High protection from H/D exchange by both peptides is consistent with retention of secondary structure in the gas phase (or refolding into some other compact structure). The P14A mutant melittin exhibits lower ECD fragmentation efficiency than WT melittin, suggesting that it contains more secondary structure in the gas phase, which may indicate that these peptides retain some memory of their solution-phase structures. Examination of the isotopic distributions of fragment ions derived from H/D exchange with subsequent ECD reveals that the C-terminus of these peptides adopts multiple conformations. The results reported here offer insight into the stability of alpha helices in the gas phase, and also highlight the value of combining gas-phase H/D exchange with electron capture dissociation to interrogate gaseous peptide conformation.

Solvent Mediation of Peptide Conformations: Polyproline Structures in Water, Methanol, Ethanol, and 1-Propanol as Determined by Ion Mobility Spectrometry-Mass Spectrometry

Ion mobility spectrometry and circular dichroism spectroscopy are used to examine the populations of the small model peptide, polyproline-13 in water, methanol, ethanol, and 1-propanol over a range of solution temperatures (from 288 to 318 K). At low temperatures, the less-polar solvents (1-propanol and ethanol) favor the all-cis polyproline I helix (PPI); as the temperature is increased, the trans-configured polyproline II helix (PPII) is formed. In polar solvents (methanol and water), PPII is favored at all temperatures. From the experimental data, we determine the relative stabilities of the eight structures in methanol, ethanol, and 1-propanol, as well as four in water, all with respect to PPII. Although these conformers show relatively small differences in free energies, substantial variability is observed in the enthalpies and entropies across the structures and solvents. This requires that enthalpies and entropies be highly correlated: in 1-propanol, cis-configured PPI conformations are energetically favorable but entropically disfavored. In more polar solvents, PPI is enthalpically less favorable and entropy favors trans-configured forms. While either ΔH0 or ΔS0 can favor different structures, no conformation in any solvent is simultaneously energetically and entropically stabilized. These data present a rare opportunity to examine the origin of conformational stability. Graphical Abstract ᅟ.

The axial/equatorial conformational landscape and intramolecular dispersion: new insights from the rotational spectra of monoterpenoids

Using rotational spectroscopy and quantum chemistry calculations, we show that intramolecular dispersion stabilises the axial conformers of monoterpenoids, and that an accurate account of these interactions is challenging for theoretical methods.

Monitoring Metallo-Macromolecular Assembly Equilibria by Ion Mobility-Mass Spectrometry

Ion mobility-mass spectrometry (IM-MS) allows the separation of isomeric and isobaric species on the basis of their size, shape, and charge. The fast separation timescale (ms) and high sensitivity of these measurements make IM-MS an ideally suitable method for monitoring changes in macromolecular structure, such as those occurring in interconverting terpyridine-based metallosupramolecular self-assemblies. IM-MS is used to verify the elemental composition (size) and architecture (shape) of the self-assembled products. Additionally, this article demonstrates its applicability to the elucidation of concentration-driven association-dissociation (fusion-fission) equilibria between isobaric structures. IM-MS enables both quantitative separation and identification of the interconverting complexes as well as derivation of the corresponding equilibrium constants (i.e., thermodynamic information) from extracted IM-MS abundance data.

Conformer Selection by Matter-Wave Interference

We establish that matter-wave diffraction at near-resonant ultraviolet optical gratings can be used to spatially separate individual conformers of complex molecules. Our calculations show that the conformational purity of the prepared beam can be close to 100% and that all molecules remain in their electronic ground state. The proposed technique is independent of the dipole moment and the spin of the molecule and thus paves the way for structure-sensitive experiments with hydrocarbons and biomolecules, such as neurotransmitters and hormones, which have evaded conformer-pure isolation so far.

Isolated Collagen Mimetic Peptide Assemblies Have Stable Triple-Helix Structures

The origin of the triple-helix structure and high stability of collagen has been debated for many years. As models of the triple helix and building blocks for new biomaterials, collagen mimetic peptide (CMP) assemblies have been deeply studied in the condensed phase. In particular, it was found that hydroxylation of proline, an abundant post-translational modification in collagen, increases its stability. Two main hypotheses emerged to account for this behavior: 1) intra-helix stereoelectronic effects, and 2) the role of water molecules H-bound to hydroxyproline side-chains. However, in condensed-phase investigations, the influence of water cannot be fully removed. Therefore, we employed a combination of tandem ion mobility and mass spectrometries to assess the structure and stability of CMP assemblies in the gas phase. These results show a conservation of the structure and stability properties of triple helix models in the absence of solvent, supporting an important role of stereoelectronic effects. Moreover, evidence that small triple helix assemblies with controlled stoichiometry can be studied in the gas phase is given, which opens new perspectives in the understanding of the first steps of collagen fiber growth.

Assessing the stability of alanine-based helices by conformer-selective IR spectroscopy

Polyalanine based peptides that carry a lysine at the C-terminus ([Ac-AlanLys + H](+)) are known to form α-helices in the gas phase. Three factors contribute to the stability of these helices: (i) the interaction between the helix macro dipole and the charge, (ii) the capping of dangling C[double bond, length as m-dash]O groups by lysine and (iii) the cooperative hydrogen bond network. In previous studies, the influence of the interaction between the helix dipole and the charge as well as the impact of the capping was studied intensively. Here, we complement these findings by systematically assessing the third parameter, the H-bond network. In order to selectively remove one H-bond along the backbone, we use amide-to-ester substitutions. The resulting depsi peptides were analyzed by ion-mobility and m/z-selective infrared spectroscopy as well as theoretical calculations. Our results indicate that peptides which contain only one ester bond still maintain the helical conformation. We conclude that the interaction between the charge and the helix macro-dipole is most crucial for the formation of the α-helical conformation and a single backbone H-bond has only little influence on the overall stability.

Combining Structural Probes in the Gas Phase - Ion Mobility-Resolved Action-FRET

In the context of native mass spectrometry, the development of gas-phase structural probes sensitive to the different levels of structuration of biomolecular assemblies is necessary to push forward conformational studies. In this paper, we provide the first example of the combination of ion mobility (IM) and Förster resonance energy transfer (FRET) measurements within the same experimental setup. The possibility to obtain mass- and mobility-resolved FRET measurements is demonstrated on a model peptide and applied to monitor the collision-induced unfolding of ubiquitin. Graphical Abstract ᅟ.

Protein folding, misfolding and aggregation: The importance of two-electron stabilizing interactions

Proteins associated with neurodegenerative diseases are highly pleiomorphic and may adopt an all-α-helical fold in one environment, assemble into all-β-sheet or collapse into a coil in another, and rapidly polymerize in yet another one via divergent aggregation pathways that yield broad diversity of aggregates’ morphology. A thorough understanding of this behaviour may be necessary to develop a treatment for Alzheimer’s and related disorders. Unfortunately, our present comprehension of folding and misfolding is limited for want of a physicochemical theory of protein secondary and tertiary structure. Here we demonstrate that electronic configuration and hyperconjugation of the peptide amide bonds ought to be taken into account to advance such a theory. To capture the effect of polarization of peptide linkages on conformational and H-bonding propensity of the polypeptide backbone, we introduce a function of shielding tensors of the C^α atoms. Carrying no information about side chain-side chain interactions, this function nonetheless identifies basic features of the secondary and tertiary structure, establishes sequence correlates of the metamorphic and pH-driven equilibria, relates binding affinities and folding rate constants to secondary structure preferences, and manifests common patterns of backbone density distribution in amyloidogenic regions of Alzheimer’s amyloid β and tau, Parkinson’s α-synuclein and prions. Based on those findings, a split-intein like mechanism of molecular recognition is proposed to underlie dimerization of Aβ, tau, αS and PrP^C, and divergent pathways for subsequent association of dimers are outlined; a related mechanism is proposed to underlie formation of PrP^Sc fibrils. The model does account for: (i) structural features of paranuclei, off-pathway oligomers, non-fibrillar aggregates and fibrils; (ii) effects of incubation conditions, point mutations, isoform lengths, small-molecule assembly modulators and chirality of solid-liquid interface on the rate and morphology of aggregation; (iii) fibril-surface catalysis of secondary nucleation; and (iv) self-propagation of infectious strains of mammalian prions.

Enantioselective Collision-Activated Dissociation of Gas-Phase Tryptophan Induced by Chiral Recognition of Protonated L-Alanine Peptides

Enantioselective dissociation in the gas phase is important for enantiomeric enrichment and chiral transmission processes in molecular clouds regarding the origin of homochirality in biomolecules. Enantioselective collision-activated dissociation (CAD) of tryptophan (Trp) and the chiral recognition ability of L-alanine peptides (L-Ala n ; n = 2-4) were examined using a linear ion trap mass spectrometer. CAD spectra of gas-phase heterochiral H+(D-Trp)(L-Ala n ) and homochiral H+(L-Trp)(L-Ala n ) noncovalent complexes were obtained as a function of the peptide size n. The H2O-elimination product was observed in CAD spectra of both heterochiral and homochiral complexes for n = 2 and 4, and in homochiral H+(L-Trp)(L-Ala3), indicating that the proton is attached to the L-alanine peptide, and H2O loss occurs from H+(L-Ala n ) in the noncovalent complexes. H2O loss did not occur in heterochiral H+(D-Trp)(L-Ala3), where NH3 loss and (H2O + CO) loss were the primary dissociation pathways. In heterochiral H+(D-Trp)(L-Ala3), the protonation site is the amino group of D-Trp, and NH3 loss and (H2O + CO) loss occur from H+(D-Trp). L-Ala peptides recognize D-Trp through protonation of the amino group for peptide size n = 3. NH3 loss and (H2O + CO) loss from H+(D-Trp) proceeds via enantioselective CAD in gas-phase heterochiral H+(D-Trp)(L-Ala3) at room temperature, whereas L-Trp dissociation was not observed in homochiral H+(L-Trp)(L-Ala3). These results suggest that enantioselective dissociation induced by chiral recognition of L-Ala peptides through protonation could play an important role in enantiomeric enrichment and chiral transmission processes of amino acids.

Probing the Gaseous Structure of a β-Hairpin Peptide with H/D Exchange and Electron Capture Dissociation

An improved understanding of the extent to which native protein structure is retained upon transfer to the gas phase promises to enhance biological mass spectrometry, potentially streamlining workflows and providing fundamental insights into hydration effects. Here, we investigate the gaseous conformation of a model β-hairpin peptide using gas-phase hydrogen-deuterium (H/D) exchange with subsequent electron capture dissociation (ECD). Global gas-phase H/D exchange levels, and residue-specific exchange levels derived from ECD data, are compared among the wild type 16-residue peptide GB1p and several variants. High protection from H/D exchange observed for GB1p, but not for a truncated version, is consistent with the retention of secondary structure of GB1p in the gas phase or its refolding into some other compact structure. Four alanine mutants that destabilize the hairpin in solution show levels of protection similar to that of GB1p, suggesting collapse or (re)folding of these peptides upon transfer to the gas phase. These results offer a starting point from which to understand how a key secondary structural element, the β-hairpin, is affected by transfer to the gas phase. This work also demonstrates the utility of a much-needed addition to the tool set that is currently available for the investigation of the gaseous conformation of biomolecules, which can be employed in the future to better characterize gaseous proteins and protein complexes. Graphical Abstract ᅟ.

The fate of proteins in outer space

It is well established that any properly conducted biophysical studies of proteins must take appropriate account of solvent. For water-soluble proteins it has been an article of faith that water is largely responsible for stabilizing the fold, a notion that has recently come under increasing scrutiny. Further, there are some instances when proteins are studied experimentally in the absence of solvent, as in matrix-assisted laser desorption/ionization or electrospray mass spectrometry, for example, or in organic solvents for protein engineering purposes. Apart from these considerations, there is considerable speculation as to whether there is life on planets other than Earth, where conditions including the presence of water (both in liquid or vapor form and indeed ice), temperature and pressure may be vastly different from those prevailing on Earth. Mars, for example, has only 0.6% of Earth's mean atmospheric pressure which presents profound problems to protein structures, as this paper and a large corpus of experimental work demonstrate. Similar objections will most likely apply in the case of most exoplanets and other bodies such as comets whose chemistry and climate are still largely unknown. This poses the question, how do proteins survive in these different environments? In order to cast some light on these issues we have conducted a series of molecular dynamics simulations on protein dehydration under a variety of conditions. We find that, while proteins undergoing dehydration can retain their integrity for a short duration they ultimately become disordered, and we further show that the disordering can be retarded if superficial water is kept in place on the surface. These findings are compared with other published results on protein solvation in an astrobiological and astrochemical setting. Inter alia, our results suggest that there are limits as to what to expect in terms of the existence of possible extraterrestrial forms as well to what can be achieved in experimental investigations on living systems despatched from Earth. This finding may appear to undermine currently held hopes that life will be found on nearby planets, but it is important to be aware that the presence of ice and water are by themselves not sufficient; there has to be an atmosphere which includes water vapor at a sufficiently high partial pressure for proteins to be active. A possible scenario in which there has been a history of adequate water vapor pressure which allowed organisms to prepare for a future dessicated state by forming suitable protective capsules cannot of course be ruled out.

Mass Spectrometry and Ion Mobility Characterization of Bioactive Peptide-Synthetic Polymer Conjugates

The bioconjugate BMP2-(PEO-HA)₂, composed of a dendron with two monodisperse poly(ethylene oxide) (PEO) branches terminated by a hydroxyapatite binding peptide (HA), and a focal point substituted with a bone growth stimulating peptide (BMP2), has been comprehensively characterized by mass spectrometry (MS) methods, encompassing matrix-assisted laser desorption ionization (MALDI), electrospray ionization (ESI), tandem mass spectrometry (MS²), and ion mobility mass spectrometry (IM-MS). MS² experiments using different ion activation techniques validated the sequences of the synthetic, bioactive peptides HA and BMP2, which contained highly basic amino acid residues either at the N-terminus (BMP2) or C-terminus (HA). Application of MALDI-MS, ESI-MS, and IM-MS to the polymer-peptide biomaterial confirmed its composition. Collision cross-section measurements and molecular modeling indicated that BMP2-(PEO-HA)₂ exists in several folded and extended conformations, depending on the degree of protonation. Protonation of all basic sites of the hybrid material nearly doubles its conformational space and accessible surface area.

Long-Lived Intermediates in a Cooperative Two-State Folding Transition

Biomolecular folding often occurs through a cooperative two-state reactant ↔ product transition; the term cooperative does not convey that intermediate structures are nonexistent but rather that these states are not observable by existing experimental techniques. Because of this, few intermediates have been studied and characterized. Recently, ion mobility spectrometry (IMS) measurements revealed that the oligomer polyproline-13 (Pro13, which in propanol (PrOH) favors the right-handed helical PPI structure having adjacent pyrrolidine rings in a cis configuration) folds through six sequential long-lived intermediates as it converts to the all-trans-configured PPII structure that is favored in aqueous solutions. Here, we examine the PPI_PrOH → PPII_aq folding transition for a HisPro13 sequence, i.e., Pro13 having a single histidine residue added to the N-terminus. Remarkably, the IMS measurements show that, upon addition of histidine, all of the IMS peaks associated with intermediate structures disappear. Instead, HisPro13 folds via a cooperative two-state transition, delayed by a significant induction period. The induction period is temperature dependent-shifting the transition to longer times at lower temperatures. Equilibrium studies show that the HisPro13 PPI_PrOH → PPII_aq transition is endothermic but favored entropically. From these clues, we propose a sequential folding mechanism and develop a model that suggests that ∼13-17 long-lived intermediates are likely responsible for the induction period. In this model, intermediates are separated by average individual activation barriers of ∼90 kJ·mol^-1, and are entropically favorable.

Structures of Metalloporphyrin-Oligomer Multianions: Cofacial versus Coplanar Motifs as Resolved by Ion Mobility Spectrometry

We have combined ion mobility mass spectrometry with quantum chemical calculations to investigate the gas-phase structures of multiply negatively charged oligomers of meso-tetra(4-sulfonatophenyl)metalloporphyrins comprising the divalent metal centers Zn^II, Cu^II, and Pd^II. Sets of candidate structures were obtained by geometry optimizations based on calculations at both the semiempirical PM7 and density functional theory (DFT) levels. The corresponding theoretical cross sections were calculated with the projection approximation and also with the trajectory method. By comparing these collision cross sections with the respective experimental values we were able to assign oligomer structures up to the tetramer. In most cases the cross sections of the lowest energy isomers predicted by theory were found to agree with the measurements to within the experimental uncertainty (2%). Specifically, we find that for a given oligomer size the structures are independent of the metal center but depend strongly on the charge state. Oligomers in low charge states with a correspondingly larger number of sodium counterions tend to form stacked, cofacial structures reminiscent of H-aggregate motifs observed in solution. By contrast, in higher charge states, the stack opens to form coplanar structures.

Collision cross section calibrants for negative ion mode traveling wave ion mobility-mass spectrometry

Unlike traditional drift-tube ion mobility-mass spectrometry, traveling-wave ion mobility-mass spectrometry typically requires calibration in order to generate collision cross section (CCS) values. Although this has received a significant amount of attention for positive-ion mode analysis, little attention has been paid for CCS calibration in negative ion mode. Here, we provide drift-tube CCS values for [M - H](-) ions of two calibrant series, polyalanine and polymalic acid, and evaluate both types of calibrants in terms of the accuracy and precision of the traveling-wave ion mobility CCS values that they produce.

Molecular Dynamics Simulation for the Dynamics and Kinetics of Folding Peptides in the Gas Phase

The conformations of flexible molecular species, such as oligomers and oligopeptides, and their interconversion in the gas phase have been probed by ion mobility spectrometry measurements. The ion motion is interpreted through the calculation of effective cross sections in the case of stable conformations of the macromolecules. However, when the molecular structures transform to each other as the ions collide with gas atoms during their flight through the drift tube, the introduction of an average cross section is required. To provide a direct way for the reproduction of the ion motion, we employ a nonequilibrium molecular dynamics simulation method and consider a molecular model that consists of two connected stiff cylindrical bodies interacting through an intramolecular model potential. With this procedure we have calculated the ion mobility as a function of temperature for a prototype peptide that converts between a helical and an extended globular form. The results are in good agreement with ion mobility spectrometry data confirming that an angular vibration coordinate can be used for the interpretation of the shifting of the drift-time distributions at high temperatures. The approach produces mean kinetic energies as well as various combined distributions of the ion degrees of freedom. It is easily applied to flexible macromolecular ions and can be extended to include additional degrees of freedom.

Going clean: structure and dynamics of peptides in the gas phase and paths to solvation

The gas phase is an artificial environment for biomolecules that has gained much attention both experimentally and theoretically due to its unique characteristic of providing a clean room environment for the comparison between theory and experiment. In this review we give an overview mainly on first-principles simulations of isolated peptides and the initial steps of their interactions with ions and solvent molecules: a bottom up approach to the complexity of biological environments. We focus on the accuracy of different methods to explore the conformational space, the connections between theory and experiment regarding collision cross section evaluations and (anharmonic) vibrational spectra, and the challenges faced in this field.

Extending native mass spectrometry approaches to integral membrane proteins

Recent developments in native mass spectrometry and ion mobility have made it possible to analyze the composition and structure of membrane protein complexes in the gas-phase. In this short review we discuss the experimental strategies that allow to elucidate aspects of the dynamic structure of these important drug targets, such as the structural effects of lipid binding or detection of co-populated conformational and assembly states during gating on an ion channel. As native mass spectrometry relies on nano-electrospray of natively reconstituted proteins, a number of commonly used lipid- and detergent-based reconstitution systems have been evaluated for their compatibility with this approach, and parameters for the release of intact, native-like folded membrane proteins studied in the gas-phase. The strategy thus developed can be employed for the investigation of the subunit composition and stoichiometry, oligomeric state, conformational changes, and lipid and drug binding of integral membrane proteins.

Capping Motif for Peptide Helix Formation

It is known that a C-terminal lysine stabilizes helix formation in polyalanine peptides that have seven or more residues. Using a combination of cold ion spectroscopy and DFT calculations, we demonstrate that even a three-residue peptide, Ac-Phe-Ala-LysH(+), adopts a structure in which the lysine side chain forms three hydrogen bonds with backbone carbonyls, reproducing the capping motif of larger polyalanine helices. This is confirmed by comparison with Ac-Phe-(Ala)5-LysH(+), which forms a 310 helix containing the same structural feature. In both molecules, we identified the vibrational bands of the N- and C-terminal amide NH stretches, which lack strong hydrogen bonds with carbonyls and consequently appear in a characteristic region above 3400 cm(-1). A similar pattern is also present in the even longer peptide Ac-Phe-(Ala)10-LysH(+), illustrating the generality of this capping motif. The two longer peptides contain additional, characteristic amide NH stretch bands below 3400 cm(-1), which form the core of the helix.

Native like helices in a specially designed β peptide in the gas phase

First principles simulations and gas phase spectroscopy suggest equilibrium of helices for an oligomer of open chain β amino acids.