Does Thermal Breathing Affect Collision Cross Sections of Gas-Phase Peptide Ions? An Ab Initio Molecular Dynamics Study

Second-Order Mass-Weighting Scheme for Atom-Centered Density Matrix Propagation Molecular Dynamics

The atom-centered density matrix propagation (ADMP) method is an extended Lagrangian approach to ab initio molecular dynamics, which includes the density matrix in an orthonormalized atom-centered Gaussian basis as additional, fictitious, electronic degrees of freedom, classically propagated along with the nuclear ones. A high adiabaticity between the nuclear and electronic subsystems is mandatory in order to keep the trajectory close to the Born-Oppenheimer (BO) surface. In this regard, the fictitious electronic mass μ, being a symmetric, nondiagonal matrix in its most general form, represents a free parameter, exploitable to optimize the propagation of the electronic density. Although mass-weighting schemes in ADMP exist, a systematic procedure to define an optimal value of the fictitious masses is not available yet. In this work, in order to rationally evaluate the electronic mass, fictitious electronic normal modes are defined through the diagonalization of the Hessian of the electronic density matrix. If the same frequency is imposed on all such modes (compatible with the chosen integration time step), then the corresponding μ matrix can be calculated and then employed for the following propagation. Analysis of several ADMP test simulations reveals that such Hessian-based mass-weighting approach is able to ensure, together with a 0.1/0.2 fs time steps, a high separation between the (real) nuclear and the (fictitious) electronic frequencies, which determines a high adiabaticity. This high, unprecedented, accuracy in the propagation leads, in turn, to low errors in the estimated nuclear vibrational frequencies, making the ADMP method totally comparable to a fully converged BO molecular dynamics simulation but more computationally efficient. This work, therefore, contributes to a further development of the ADMP ab initio molecular dynamics method, aimed at improving its accuracy through a more rational evaluation of the fictitious electronic mass parameter.

Effects of charge on protein ion structure: Lessons from cation-to-anion, proton-transfer reactions

Collision cross-section values, which can be determined using ion mobility experiments, are sensitive to the structures of protein ions and useful for applications to structural biology and biophysics. Protein ions with different charge states can exhibit very different collision cross-section values, but a comprehensive understanding of this relationship remains elusive. Here, we review cation-to-anion, proton-transfer reactions (CAPTR), a method for generating a series of charge-reduced protein cations by reacting quadrupole-selected cations with even-electron monoanions. The resulting CAPTR products are analyzed using a combination of ion mobility, mass spectrometry, and collisional activation. We compare CAPTR to other charge-manipulation strategies and review the results of various CAPTR-based experiments, exploring their contribution to a deeper understanding of the relationship between protein ion structure and charge state.

Imaging of femtosecond bond breaking and charge dynamics in ultracharged peptides

X-ray free-electrons lasers have revolutionized the method of imaging biological macromolecules such as proteins, viruses and cells by opening the door to structural determination of both single particles and crystals at room temperature. By utilizing high intensity X-ray pulses on femtosecond timescales, the effects of radiation damage can be reduced. Achieving high resolution structures will likely require knowledge of how radiation damage affects the structure on an atomic scale, since the experimentally obtained electron densities will be reconstructed in the presence of radiation damage. Detailed understanding of the expected damage scenarios provides further information, in addition to guiding possible corrections that may need to be made to obtain a damage free reconstruction. In this work, we have quantified the effects of ionizing photon-matter interactions using first principles molecular dynamics. We utilize density functional theory to calculate bond breaking and charge dynamics in three ultracharged molecules and two different structural conformations that are important to the structural integrity of biological macromolecules, comparing to our previous studies on amino acids. The effects of the ultracharged states and subsequent bond breaking in real space are studied in reciprocal space using coherent diffractive imaging of an ensemble of aligned biomolecules in the gas phase.

Probing d- and l-Adrenaline Binding to β2-Adrenoreceptor Peptide Motifs by Gas-Phase Photodissociation Cross-Linking and Ion Mobility Mass Spectrometry

Diazirine-tagged d- and l-adrenaline derivatives formed abundant noncovalent gas-phase ion complexes with peptides N-Ac-SSIVSFY-NH₂ (peptide S) and N-Ac-VYILLNWIGY-NH₂ (peptide V) upon electrospray ionization. These peptide sequences represent the binding motifs in the β₂-adrenoreceptor. The structures of the gas-phase complexes were investigated by selective laser photodissociation of the diazirine chromophore at 354 nm, which resulted in a loss of N₂ and formation of a transient carbene intermediate in the adrenaline ligand without causing its expulsion. The photolyzed complexes were analyzed by collision-induced dissociation (CID-MS³ and CID-MS⁴) in an attempt to detect cross-links and establish the binding sites. However, no cross-linking was detected in the complexes regardless of the peptide and d- or l-configuration in adrenaline. Cyclic ion mobility measurements were used to obtain collision cross sections (CCS) in N₂ for the peptide S complexes. These showed identical values, 334 ± 0.9 Ų, for complexes of the l- and d-adrenaline derivatives, respectively. Identical CCS were also obtained for peptide S complexes with natural l- and d-adrenaline, 317 ± 1.2 Ų, respectively. Born-Oppenheimer molecular dynamics (BOMD) in combination with full geometry optimization by density functional theory calculations provided structures for the complexes that were used to calculate theoretical CCS with the ion trajectory method. A close match (337 Ų) was found for a single low Gibbs energy structure that displayed a binding pocket with Ser 2 and Ser 5 residues forming hydrogen bonds to the adrenaline catechol hydroxyls. Analysis of the BOMD trajectories revealed a small number of contacts between the incipient carbene carbon atom in the ligand and X-H bonds in the peptide, which was consistent with the lack of cross-linking. Temperature dependence of the internal dynamics of peptide S-adrenaline complexes as well as the specifics of the adrenaline carbene reactions are discussed. In particular, peptide amide hydrogen transfer to the carbene carbon atom was calculated to require crossing a potential energy barrier, which may hamper cross-linking in competition with carbene internal rearrangements.

A Not Obvious Correlation Between the Structure of Green Fluorescent Protein Chromophore Pocket and Hydrogen Bond Dynamics: A Choreography From ab initio Molecular Dynamics

The Green Fluorescent Protein (GFP) is a widely studied chemical system both for its large amount of applications and the complexity of the excited state proton transfer responsible of the change in the protonation state of the chromophore. A detailed investigation on the structure of the chromophore environment and the influence of chromophore form (either neutral or anionic) on it is of crucial importance to understand how these factors could potentially influence the protein function. In this study, we perform a detailed computational investigation based on the analysis of ab-initio molecular dynamics simulations, to disentangle the main structural quantities determining the fine balance in the chromophore environment. We found that specific hydrogen bonds interactions directly involving the chromophore (or not), are correlated to quantities, such as the volume of the cavity in which the chromophore is embedded and that it is importantly affected by the chromophore protonation state. The cross-correlation analysis performed on some of these hydrogen bonds and the cavity volume, demonstrates a direct correlation among them and we also identified the ones specifically involved in this correlation. We also found that specific interactions among residues far in the space are correlated, demonstrating the complexity of the chromophore environment and that many structural quantities have to be taken into account to properly describe and understand the main factors tuning the active site of the protein. From an overall evaluation of the results obtained in this work, it is shown that the residues which a priori are perceived to be spectators play instead an important role in both influencing the chromophore environment (cavity volume) and its dynamics (cross-correlations among spatially distant residues).

Multiresolution continuous wavelet transform for studying coupled solute-solvent vibrations via ab initio molecular dynamics

Vibrational analysis in solution and the theoretical determination of infrared and Raman spectra are of key importance in many fields of chemical interest. Vibrational band dynamics of molecules and their sensitivity to the environment can also be captured by these spectroscopies in their time dependent version. However, it is often difficult to provide an interpretation of the experimental data at the molecular scale, such as molecular mechanisms or the processes hidden behind them. In this work, we present a theoretical-computational protocol based on ab initio molecular dynamics simulations and a combination of normal-like (generalized) mode analysis of solute-solvent clusters with a wavelet transform, for the first time. The case study is the vibrational dynamics of N-methyl-acetamide (NMA) in water solution, a well-known model of hydration of peptides and proteins. Amide modes are typical bands of peptide and protein backbone, and their couplings with the environment are very challenging in terms of the accurate prediction of solvent induced intensity and frequency shifts. The contribution of water molecules surrounding NMA to the composition of generalized and time resolved modes is introduced in our vibrational analysis, showing unequivocally its influence on the amide mode spectra. It is also shown that such mode compositions need the inclusion of the first shell solvent molecules to be accurately described. The wavelet analysis is proven to be strongly recommended to follow the time evolution of the spectra, and to capture vibrational band couplings and frequency shifts over time, preserving at the same time a well-balanced time-frequency resolution. This peculiar feature also allows one to perform a combined structural-vibrational analysis, where the different strengths of hydrogen bond interactions can quantitatively affect the amide bands over time at finite temperature. The proposed method allows for the direct connection between vibrational modes and local structural changes, providing a link from the spectroscopic observable to the structure, in this case the peptide backbone, and its hydration layouts.

Ab-initio molecular dynamics and hybrid explicit-implicit solvation model for aqueous and nonaqueous solvents: GFP chromophore in water and methanol solution as case study

Solute-solvent interactions are proxies for understanding how the electronic density of a chromophore interacts with the environment in a more exhaustive way. The subtle balance between polarization, electrostatic, and non-bonded interactions need to be accurately described to obtain good agreement between simulations and experiments. First principles approaches providing accurate configurational sampling through molecular dynamics may be a suitable choice to describe solvent effects on solute chemical-physical properties and spectroscopic features, such as optical absorption of dyes. In this context, accurate energy potentials, obtained by hybrid implicit/explicit solvation methods along with employing nonperiodic boundary conditions, are required to represent bulk solvent around a large solute-solvent cluster. In this work, a novel strategy to simulate methanol solutions is proposed combining ab initio molecular dynamics, a hybrid implicit/explicit flexible solvent model, nonperiodic boundary conditions, and time dependent density functional theory. As case study, the robustness of the proposed protocol has been gauged by investigating the microsolvation and electronic absorption of the anionic green fluorescent protein chromophore in methanol and aqueous solution. Satisfactory results are obtained, reproducing the microsolvation layout of the chromophore and, as a consequence, the experimental trends shown by the optical absorption in different solvents.

Untangling Hydrogen Bond Networks with Ion Mobility Spectrometry and Quantum Chemical Calculations: A Case Study on H+XPGG

Ion mobility spectrometry-mass spectrometry and quantum chemical calculations are used to determine the structures and stabilities of singly protonated XaaProGlyGly peptides: H⁺DPGG, H⁺NPGG, H⁺EPGG, and H⁺QPGG. The IMS distributions are similar, suggesting the peptides adopt closely related structures in the gas phase. Quantum chemical calculations show that all conformers seen in the experimental spectrum correspond to the cis configuration about the Xaa-Pro peptide bond, significantly different from the behavior seen previously for H⁺GPGG. Density functional theory and quantum theory of atoms in molecules (QTAIM) investigations uncover a silent drama as a minor conformer not observed in the H⁺DPGG spectrum becomes the preferred conformer in H⁺QPGG, with both conformers being coincident in collision cross section. Investigation of the highly coupled hydrogen bond network, replete with CH···O interactions and bifurcated hydrogen bonds, reveals the cause of this effect as well as the absence of trans conformers from the spectra. A series of generalized observations are provided to aid in enzyme and ligand design using these coupled hydrogen bond motifs.

Ion Mobility-Mass Spectrometry of Iodine Pentoxide-Iodic Acid Hybrid Cluster Anions in Dry and Humidified Atmospheres

Nanometer-scale clusters form from vapor-phase precursors and can subsequently grow into nanoparticles during atmospheric nucleation events. A particularly interesting set of clusters relevant to nucleation is hybrid iodine pentoxide-iodic acid clusters of the form (I₂O₅) _x(HIO₃) _y as these clusters have been observed in coastal region nucleation events in anomalously high concentrations. To better understand their properties, we utilized ion mobility-mass spectrometry to probe the structures of cluster anions of the form (I₂O₅) _x(HIO₃) _y(IO_α)^- ( x = 0-7, y = 0-1, α = 1-3), similar to those observed in coastal nucleation events. We show that (I₂O₅) _x(HIO₃) _y(IO_α)^- clusters are relatively stable against dissociation during mass spectrometric measurement, as compared to other clusters observed in nucleation events over continental sites, and that at atmospherically relevant relative humidity levels (65% and less) clusters can become sufficiently hydrated to facilitate complete conversion of iodine pentoxide to iodic acid but that water sorption beyond this level is limited, indicating that the clusters do not persist as nanometer-scale droplets in the ambient.

A parallelized molecular collision cross section package with optimized accuracy and efficiency

A new parallelized calculation package predicts collision cross sections with high accuracy and efficiency.

Replacing H+ by Na+ or K+ in phosphopeptide anions and cations prevents electron capture dissociation

By successively replacing H+ by Na+ or K+ in phosphopeptide anions and cations, we show that the efficiency of fragmentation into c and z˙ or c˙ and z fragments from N-Cα backbone bond cleavage by negative ion electron capture dissociation (niECD) and electron capture dissociation (ECD) substantially decreases with increasing number of alkali ions attached. In proton-deficient phosphopeptide ions with a net charge of 2-, we observed an exponential decrease in electron capture efficiency with increasing number of Na+ or K+ ions attached, suggesting that electrons are preferentially captured at protonated sites. In proton-abundant phosphopeptide ions with a net charge of 3+, the electron capture efficiency was not affected by replacing up to four H+ ions with Na+ or K+ ions, but the yield of c, z˙ and c˙, z fragments from N-Cα backbone bond cleavage generally decreased next to Na+ or K+ binding sites. We interpret the site-specific decrease in fragmentation efficiency as Na+ or K+ binding to backbone amide oxygen in competition with interactions of protonated sites that would otherwise lead to backbone cleavage into c, z˙ or c˙, z fragments. Our findings seriously challenge the hypothesis that the positive charge responsible for ECD into c, z˙ or c˙, z fragments can generally be a sodium or other metal ion instead of a proton.

Stereospecific control of peptide gas-phase ion chemistry with cis and trans cyclo ornithine residues

We report non-chiral amino acid residues cis- and trans-1,4-diaminocyclohexane-1-carboxylic acid (cyclo-ornithine, cO) that exhibit unprecedented stereospecific control of backbone dissociations of singly charged peptide cations and hydrogen-rich cation radicals produced by electron-transfer dissociation. Upon collision-induced dissociation (CID) in the slow heating regime, peptide cations containing trans-cO residues undergo facile backbone cleavages of amide bonds C-terminal to trans-cO. By contrast, peptides with cis-cO residues undergo dissociations at several amide bonds along the peptide ion backbone. Diastereoisomeric cO-containing peptides thus provide remarkably distinct tandem mass spectra. The stereospecific effect in CID of the trans-cO residue is explained by syn-facially directed proton transfer from the 4-ammonium group at cO to the C-terminal amide followed by neighboring group participation in the cleavage of the CO-NH bond, analogous to the aspartic acid and ornithine effects. Backbone dissociations of diastereoisomeric cO-containing peptide ions generate distinct [b_n ]⁺ -type fragment ions that were characterized by CID-MS³ spectra. Stereospecific control is also reported for electron-transfer dissociation of cis- and trans-cO containing doubly charged peptide ions. The stereospecific effect upon electron transfer is related to the different conformations of doubly charged peptide ions that affect the electron attachment sites and ensuing N-C_α bond dissociations.