Strengths and Weaknesses of Molecular Simulations of Electrosprayed Droplets

Molecular Dynamics Simulations of Native Protein Charging via Proton Transfer during Electrospray Ionization with Grotthuss Diffuse H3O. Anal Chem 2024;96:4146-4153. [PMID: 38427846 PMCID: PMC11337394 DOI: 10.1021/acs.analchem.3c05089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Unraveling the mechanism by which native proteins are charged through electrospray ionization (ESI) has been the focus of considerable research because observable charge states can be correlated to biophysical characteristics, such as protein folding and, thus, solution conformation. Difficulties in characterizing electrosprayed droplets have catalyzed the use of molecular dynamics (MD) to provide insights into the mechanisms by which proteins are charged and transferred to the gas phase. However, prior MD studies have utilized metal ions, primarily Na+, as charge carriers, even though proteins are primarily detected as protonated ions in the mass spectra. Here, we propose a modified MD protocol for simulating discrete Grotthuss diffuse H3O+ that is capable of dynamically altering amino-acid protonation states to model electrospray charging and gaseous ion formation of model proteins, ubiquitin, and myoglobin. Application of the protocol to the evaporation of acidic droplets enables a molecular perspective of H3O+ coordination and proton transfer to/from proteins, which is unfeasible with the metal charge carriers used in previous MD studies of ESI. Our protocol recreates experimentally observed charge-state distributions and supports the charge residue model (CRM) as the dominant mechanism of native protein ionization during ESI. Additionally, our results suggest that protonation is highly specific to individual residues and is correlated to the formation of localized hydrated regions on the protein surface as droplets desolvate. Considering the use of discrete H3O+ instead of Na+, the developed protocol is a necessary step toward developing a more comprehensive model of protein ionization during ESI.

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.

Release of nanodiscs from charged nano-droplets in the electrospray ionization revealed by molecular dynamics simulations

Electrospray ionization (ESI) is essential for application of mass spectrometry in biological systems, as it prevents the analyte being split into fragments. However, due to lack of a clear understanding of the mechanism of ESI, the interpretation of mass spectra is often ambiguous. This is a particular challenge for complex biological systems. Here, we focus on systems that include nanodiscs as membrane environment, which are essential for membrane proteins. We performed microsecond atomistic molecular dynamics simulations to study the release of nanodiscs from highly charged nano-droplets into the gas phase, the late stage of ESI. We observed two distinct major scenarios, highlighting the diversity of morphologies of gaseous product ions. Our simulations are in reasonable agreement with experimental results. Our work provides a detailed atomistic view of the ESI process of a heterogeneous system (lipid nanodisc), which may give insights into the interpretation of mass spectra of all lipid-protein systems.

Salt Enrichment and Dynamics in the Interface of Supercooled Aqueous Droplets

The interconversion reaction of NaCl between the contact-ion pair (CIP) and the solvent-separated ion pair (SSIP) as well as the free-ion state in cold droplets has not yet been investigated. We report direct computational evidence that the lower is the temperature, the closer to the surface the ion interconversion reaction takes place. In supercooled droplets the enrichment of the subsurface in salt becomes more evident. The stability of the SSIP relative to the CIP increases as the ion-pairing is transferred toward the droplet's outer layers. In the free-ion state, where the ions diffuse independently in the solution, the number density of Cl^- shows a broad maximum in the interior in addition to the well-known maximum in the surface. In the study of the reaction dynamics, we find a weak coupling between the interionic NaCl distance reaction coordinate and the solvent degrees of freedom, which contrasts with the diffusive crossing of the free energy barrier found in bulk solution modeling. The H₂O self-diffusion coefficient is found to be at least an order of magnitude larger than that in the bulk solution. We propose to exploit the enhanced surface ion concentration at low temperature to eliminate salts from droplets in native mass spectrometry ionization methods.

Grotthuss Molecular Dynamics Simulations for Modeling Proton Hopping in Electrosprayed Water Droplets

Excess protons in water exhibit unique transport properties because they can rapidly hop along H-bonded water wires. Considerable progress has been made in unraveling this Grotthuss diffusion mechanism using quantum mechanical-based computational techniques. Unfortunately, high computational cost tends to restrict those techniques to small systems and short times. Molecular dynamics (MD) simulations can be applied to much larger systems and longer time windows. However, standard MD methods do not permit the dissociation/formation of covalent bonds, such that Grotthuss diffusion cannot be captured. Here, we bridge this gap by combining atomistic MD simulations (using Gromacs and TIP4P/2005 water) with proton hopping. Excess protons are modeled as hydronium ions that undergo H₃O⁺ + H₂O → H₂O + H₃O⁺ transitions. In accordance with ab initio MD data, these Grotthuss hopping events are executed in "bursts" with quasi-instantaneous hopping across one or more waters. The bursts are separated by regular MD periods during which H₃O⁺ ions undergo Brownian diffusion. The resulting proton diffusion coefficient agrees with the literature value. We apply this Grotthuss MD technique to highly charged water droplets that are in a size regime encountered during electrospray ionization (5 nm radius, ∼17,000 H₂O). The droplets undergo rapid solvent evaporation and occasional H₃O⁺ ejection, keeping them at ca. 81% of the Rayleigh limit. The simulated behavior is consistent with phase Doppler anemometry data. The Grotthuss MD technique developed here should be useful for modeling the behavior of various proton-containing systems that are too large for high-level computational approaches. In particular, we envision future applications related to electrospray processes, where earlier simulations used metal cations while in reality excess protons dominate.

Sulfolane-Induced Supercharging of Electrosprayed Salt Clusters: An Experimental/Computational Perspective

It is well-known that supercharging agents (SCAs) such as sulfolane enhance the electrospray ionization (ESI) charge states of proteins, although the mechanistic origins of this effect remain contentious. Only very few studies have explored SCA effects on analytes other than proteins or peptides. This work examines how sulfolane affects electrosprayed NaI salt clusters. Such alkali metal halide clusters have played a key role for earlier ESI mechanistic studies, making them interesting targets for supercharging investigations. ESI of aqueous NaI solutions predominantly generated singly charged [Na_nI_(n-1)]⁺ clusters. The addition of sulfolane resulted in abundant doubly charged [Na_nI_(n-2)Sulfolane_s]²⁺ species. These experimental data for the first time demonstrate that electrosprayed salt clusters can undergo supercharging. Molecular dynamics (MD) simulations of aqueous ESI nanodroplets containing Na⁺/I^- with and without sulfolane were conducted to obtain atomistic insights into the supercharging mechanism. The simulations produced [Na_nI_i]^z+ and [Na_nI_iSulfolane_s]^z+ clusters similar to those observed experimentally. The MD trajectories demonstrated that these clusters were released into the gas phase upon droplet evaporation to dryness, in line with the charged residue model. Sulfolane was found to evaporate much more slowly than water. This slow evaporation, in conjunction with the large dipole moment of sulfolane, resulted in electrostatic stabilization of the shrinking ESI droplets and the final clusters. Hence, charge-dipole stabilization causes the sulfolane-containing droplets and clusters to retain more charge, thereby providing the mechanistic foundation of salt cluster supercharging.

Molecular Characterization of the Surface Excess Charge Layer in Droplets

The surface excess charge layer (SECL) in droplets has often been associated with distinct chemistry. We examine the effect of the nature of ions in the composition and structure of SECL by using molecular dynamics. We find that in the presence of simple ions the thickness of SECL is invariant not only with respect to droplet size but also with respect to the nature of the ions. In the presence of simple ions, this layer has a thickness of ∼1.5-1.7 nm but in the presence of macroions it may extend to ∼2.0 nm. The proportion of ions contained in SECL depends on the nature of the ions and the droplet size. For the same droplet size, I^- and model H₃O⁺ ions show considerably higher concentration than Na⁺ and Cl^- ions. We identify the maximum ion concentration region, which, in nanodrops, may partially overlap with SECL. As the relative shape fluctuations decrease when microdrop size is approached, the overlap between SECL and maximum ion concentration region increases. We suggest the extension of the bilayer droplet structure assumed in the equilibrium partitioning model of Enke to include the maximum ion concentration region that may not coincide with SECL in nanodrops. We compute the ion concentrations in SECL, which are those that should enter the kinetic equation in the ion-evaporation mechanism, instead of the overall drop ion concentration that has been used.

An overview of biological applications and fundamentals of new inlet and vacuum ionization technologies

RATIONALE

The developments of new ionization technologies based on processes previously unknown to mass spectrometry (MS) have gained significant momentum. Herein we address the importance of understanding these unique ionization processes, demonstrate the new capabilities currently unmet by other methods, and outline their considerable analytical potential.

METHODS

The inlet and vacuum ionization methods of solvent-assisted ionization (SAI), matrix-assisted ionization (MAI), and laserspray ionization can be used with commercial and dedicated ion sources producing ions from atmospheric or vacuum conditions for analyses of a variety of materials including drugs, lipids, and proteins introduced from well plates, pipet tips and plate surfaces with and without a laser using solid or solvent matrices. Mass spectrometers from various vendors are employed.

RESULTS

Results are presented highlighting strengths relative to ionization methods of electrospray ionization (ESI) and matrix-assisted laser desorption/ionization. We demonstrate the utility of multi-ionization platforms encompassing MAI, SAI, and ESI and enabling detection of what otherwise is missed, especially when directly analyzing mixtures. Unmatched robustness is achieved with dedicated vacuum MAI sources with mechanical introduction of the sample to the sub-atmospheric pressure (vacuum MAI). Simplicity and use of a wide array of matrices are attained using a conduit (inlet ionization), preferably heated, with sample introduction from atmospheric pressure. Tissue, whole blood, urine (including mouse, chicken, and human origin), bacteria strains and chemical on-probe reactions are analyzed directly and, especially in the case of vacuum ionization, without concern of carryover or instrument contamination.

CONCLUSIONS

Examples are provided highlighting the exceptional analytical capabilities associated with the novel ionization processes in MS that reduce operational complexity while increasing speed and robustness, achieving mass spectra with low background for improved sensitivity, suggesting the potential of this simple ionization technology to drive MS into areas currently underserved, such as clinical and medical applications.

Ion Formation from Rapidly Heated Aqueous Droplets by Droplet-Assisted Ionization

When aqueous droplets travel through a temperature-controlled capillary from atmospheric pressure into a vacuum, they undergo aerodynamic and/or thermal breakup to give charged progeny droplets that subsequently produce gas-phase molecular ions from solutes that were in the original droplets. This phenomenon is the basis of droplet-assisted ionization, a method that was recently developed for online characterization of aerosols by mass spectrometry. The conditions allowing initial droplets to break up into progeny droplets were studied by computational fluid dynamics (CFD) with a droplet evaporation model. The CFD results were then used to interpret experimental measurements of ion current vs capillary wall temperature. For capillary wall temperatures below about 150 °C, the abilities of droplets to undergo either aerodynamic or thermal breakup are strongly temperature dependent. Above this temperature, the mode of initial droplet breakup becomes temperature independent, and the temperature dependence of the ion signal intensity can be explained in relation to ion formation from charged progeny droplets. Activation energies for ion formation fall into two main categories: ∼41 kJ mol^-1 for droplets containing predominantly nonionic solutes, which matches the enthalpy of vaporization for water and suggests a charge residue process for ion formation, and ∼24 kJ mol^-1 for droplets containing salts, which suggests an ion evaporation process where the ion is ejected from the droplet surface within a cluster of solvent molecules.

Enhancing Protein Electrospray Charge States by Multivalent Metal Ions: Mechanistic Insights from MD Simulations and Mass Spectrometry Experiments

The structure and reactivity of electrosprayed protein ions is governed by their net charge. Native proteins in non-denaturing aqueous solutions produce low charge states. More highly charged ions are formed when electrospraying proteins that are unfolded and/or exposed to organic supercharging agents. Numerous studies have explored the electrospray process under these various conditions. One phenomenon that has received surprisingly little attention is the charge enhancement caused by multivalent metal ions such as La³⁺ when electrospraying proteins out of non-denaturing solutions. Here, we conducted mass spectrometry and ion mobility spectrometry experiments, in combination with molecular dynamics (MD) simulations, to uncover the mechanistic basis of this charge enhancement. MD simulations of aqueous ESI droplets reproduced the experimental observation that La³⁺ boosts protein charge states relative to monovalent metals (e.g., Na⁺). The simulations showed that gaseous proteins were released by solvent evaporation to dryness, consistent with the charged residue model. Metal ion ejection kept the shrinking droplets close to the Rayleigh limit until ∼99% of the solvent had left. For droplets charged with Na⁺, metal adduction during the final stage of solvent evaporation produced low protein charge states. Droplets containing La³⁺ showed a very different behavior. The trivalent nature of La³⁺ favored adduction to the protein at a very early stage, when most of the solvent had not evaporated yet. This irreversible binding via multidentate contacts suppressed La³⁺ ejection from the vanishing droplets, such that the resulting gaseous proteins carried significantly more charge. Our results illustrate that MD simulations are suitable for uncovering intricate aspects of electrospray mechanisms, paving the way toward an atomistic understanding of mass spectrometry based analytical workflows.

Computational Insights into Compaction of Gas-Phase Protein and Protein Complex Ions in Native Ion Mobility-Mass Spectrometry

Native ion mobility-mass spectrometry (IM-MS) is a rapidly growing field for studying the composition and structure of biomolecules and biomolecular complexes using gas-phase methods. Typically, ions are formed in native IM-MS using gentle nanoelectrospray ionization conditions, which in many cases can preserve condensed-phase stoichiometry. Although much evidence shows that large-scale condensed-phase structure, such as quaternary structure and topology, can also be preserved, it is less clear to what extent smaller-scale structure is preserved in native IM-MS. This review surveys computational and experimental efforts aimed at characterizing compaction and structural rearrangements of protein and protein complex ions upon transfer to the gas phase. A brief summary of gas-phase compaction results from molecular dynamics simulations using multiple common force fields and a wide variety of protein ions is presented and compared to literature IM-MS data.