Charge-Induced Conformational Changes of PEG-(Na+)n in a Vacuum and Aqueous Nanodroplets

Modeling a Continuously Operating Electrospray Ionization Emitter Using Molecular Dynamics Simulations: From Bulk Solution to Gaseous Ions

Numerous analytical workflows involve electrospray ionization (ESI), a process that converts solution species into gaseous ions for detection by mass spectrometry (MS). Upon exposure of a solution-filled emitter capillary to an electric field, ESI proceeds via a liquid jet that forms at the apex of a Taylor cone. Breakup of this jet generates charged droplets that ultimately release analyte ions into the gas phase. Many facets of these events remain incompletely understood. Molecular dynamics (MD) simulations of ESI droplets have become an important tool for ESI mechanistic investigations. Here, we extend these modeling efforts to entire ESI emitter capillaries. Such simulations face significant challenges. Chief among these is the fact that experimental ESI emitters operate in steady-state mode, where ejected solution is constantly replenished from an upstream reservoir. Previous MD simulations employed capillaries with a finite amount of solution, resulting in only short bursts of sprayer operation. The current work develops an MD algorithm that combines trajectory stitching with "solution recycling". This algorithm periodically removes ejected droplets, while simultaneously replenishing the emitter inlet with fresh solution. For the first time, this approach provides continuous MD-based ESI operation, while keeping the instantaneous number of atoms (and hence computational cost) manageable. For a 50 ns simulation, we illustrate a reduction of wall clock time from an estimated 5000 days to 20 days. The ESI emitter MD simulations performed here are the most detailed to date, providing insights into the properties of the Taylor cone and the behavior of nascent ESI droplets. Focusing on aqueous NaCl solution, this work provides the first comprehensive perspective of the entire ESI process, from analyte solution in the emitter to MS-detectable gaseous ions.

MD Simulations of Peptide-Containing Electrospray Droplets: Effects of Parameter Settings on the Predicted Mechanisms of Gas Phase Ion Formation

Electrospray ionization (ESI) mass spectrometry is widely used for interrogating peptides, proteins, and other biomolecular analytes. A growing number of laboratories use molecular dynamics (MD) simulations for uncovering ESI mechanisms by modeling the behavior of highly charged nanodroplets. The outcome of any MD simulation depends on certain assumptions and parameter settings, and it is desirable to optimize these factors by benchmarking computational data against experiments. Unfortunately, benchmarking of ESI simulations is difficult because experimentally generated gaseous ions do not generally retain any features that would reveal their formation pathway [e.g., the charged residue mechanism (CRM) or the ion evaporation mechanism (IEM)]. Here, we tackle this problem by examining the effects of various MD settings on the ESI behavior of the 9-residue peptide bradykinin in acidic aqueous droplets. Several parameters were found to significantly affect the kinetic competition between peptide IEM and CRM. By systematically probing the droplet behavior, we uncovered problems associated with certain settings, including peptide/solvent temperature imbalances, unexpected peptide deceleration during IEM, and a dependence of the ESI mechanism on the water model. We also noted different simulation outcomes for different force fields. On the basis of comprehensive tests, we propose a set of "best practice" parameter settings for MD simulations of ESI droplets. The strategies used here should be transferable to other types of droplet simulations, paving the way toward a more solid understanding of ESI mechanisms.

Gas-phase structure of polymer ions: Tying together theoretical approaches and ion mobility spectrometry

An increasing number of studies take advantage of ion mobility spectrometry (IMS) coupled to mass spectrometry (IMS-MS) to investigate the spatial structure of gaseous ions. Synthetic polymers occupy a unique place in the field of IMS-MS. Indeed, due to their intrinsic dispersity, they offer a broad range of homologous ions with different lengths. To help rationalize experimental data, various theoretical approaches have been described. First, the study of trend lines is proposed to derive physicochemical and structural parameters. However, the evaluation of data fitting reflects the overall behavior of the ions without reflecting specific information on their conformation. Atomistic simulations constitute another approach that provide accurate information about the ion shape. The overall scope of this review is dedicated to the synergy between IMS-MS and theoretical approaches, including computational chemistry, demonstrating the essential role they play to fully understand/interpret IMS-MS data.

Limitations of Atomistic Molecular Dynamics to Reveal Ejection of Proteins from Charged Nanodroplets

Atomistic molecular dynamics (MD) is frequently used to unravel the mechanisms of macroion release from electrosprayed droplets. However, atomistic MD is currently feasible for only the smallest window of droplet sizes appearing at the end steps of a droplet's lifetime. The relevance of the observations made to the actual droplet evolution, which is much longer than the simulated sizes, has not been addressed yet in the literature. Here, we perform a systematic study of the desolvation mechanisms of poly(ethylene glycol) (PEG), protonated peptides of different compositions, and proteins, to (a) obtain insight into the charging mechanism of macromolecules in larger droplets than those that are currently amenable to atomistic MD and (b) examine whether currently used atomistic MD modeling can establish the extrusion mechanism of proteins from droplets. To mimic larger droplets that are not amenable to MD modeling, we scale down the systems, by simulating a large droplet size relative to the macromolecule. MD of PEG charging reveals that, above a critical droplet size, ions are available near the backbone of the macromolecule, but charging occurs only transiently by transfer of ions from the solvent to the macroion, while below the critical size, the capture of the ion from PEG has a lifetime sufficiently long for the extrusion of a charged PEG from the aqueous droplet. This is the first report of the role of droplet curvature in the relation between macroion conformation and charging. Simulations of protonated peptides with a high degree of hydrophobicity show that partial extrusion of a peptide from the droplet surface is rare relative to desolvation by drying-out. Different from what has been presented in the literature, we argue that atomistic MD simulations have not sufficiently established the extrusion mechanism of proteins from droplets and their charging mechanism. We also argue that release of highly charged proteins can occur at an earlier stage of a droplet's lifetime than predicted by atomistic MD. In this earlier stage, we emphasize the key role of jets emanating from a droplet at the point of charge-induced instability in the release of proteins.

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.

Atomistic Modeling of Jet Formation in Charged Droplets

The first atomistic simulations that reveal the mechanism of Rayleigh fission are presented. It is demonstrated that simple ion or macroion ejection takes place through droplet deformation from a spherical into a distinct "tear" shape that contains a conical protrusion. We assert that the latter state is a free-energy minimum along an order parameter that measures the degree of droplet asphericity. The charged droplet's long-time evolution proceeds by alternating between the two minima above and below the critical value that are reached through solvent evaporation and ion ejection, respectively. For the first time, this mechanism allows one to explain the nature of the progeny droplets and the percentage of charge lost during fission. The cone half angle is estimated and found to be in good agreement with the value predicted from the solution of the electrostatic equation for the dielectric liquid. It is found that the conical deformation is independent of the effect of electrohydrodynamic forces reported in experiments. Contrary to the experimental observations of two diametrically opposite jets for droplets suspended in the electric field, we find that a single jet is formed at the Rayleigh limit. The study provides insight into the mechanism of capture of a macroion in jets appearing in electrospray ionization mass spectrometry (ESI-MS) experiments and may explain the tolerance of the ESI-MS spectrum to salt contamination of the sample.

Formation of Gaseous Peptide Ions from Electrospray Droplets: Competition between the Ion Evaporation Mechanism and Charged Residue Mechanism

The transfer of peptide ions from solution into the gas phase by electrospray ionization (ESI) is an integral component of mass spectrometry (MS)-based proteomics. The mechanisms whereby gaseous peptide ions are released from charged ESI nanodroplets remain unclear. This is in contrast to intact protein ESI, which has been the focus of detailed investigations using molecular dynamics (MD) simulations and other methods. Under acidic liquid chromatography/MS conditions, many peptides carry a solution charge of 3+ or 2+. Because of this pre-existing charge and their relatively small size, prevailing views suggest that peptides follow the ion evaporation mechanism (IEM). The IEM entails analyte ejection from ESI droplets, driven by electrostatic repulsion between the analyte and droplet. Surprisingly, recent peptide MD investigations reported a different behavior, that is, the release of peptide ions via droplet evaporation to dryness which represents the hallmark of the charged residue mechanism (CRM). Here, we resolved this conundrum by performing MD simulations on a common model peptide (bradykinin) in Rayleigh-charged aqueous droplets. The primary focus was on pH 2 conditions (bradykinin solution charge = 3+), but we also verified that our MD strategy captured pH-dependent charge state shifts seen in ESI-MS experiments. In agreement with earlier simulations, we found that droplets with initial radii of 1.5-3 nm predominantly release peptide ions via the CRM. In contrast, somewhat larger radii (4-5 nm) favor IEM behavior. It appears that these are the first MD data to unequivocally demonstrate the viability of peptide IEM events. Electrostatic arguments can account for the observed droplet size dependence. In summary, both CRM and IEM can be operative in peptide ESI-MS. The prevalence of one over the other mechanism depends on the droplet size distribution in the ESI plume.

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.

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.

Effects of electrospray mechanisms and structural relaxation on polylactide ion conformations in the gas phase: insights from ion mobility spectrometry and molecular dynamics simulations

Gas-phase polymer ions may retain structural features associated with their electrospray formation mechanisms.

Gas-phase conformations of intrinsically disordered proteins and their complexes with ligands: Kinetically trapped states during transfer from solution to the gas phase

Flexible structures of intrinsically disordered proteins (IDPs) are crucial for versatile functions in living organisms, which involve interaction with diverse partners. Electrospray ionization ion mobility mass spectrometry (ESI-IM-MS) has been widely applied for structural characterization of apo-state and ligand-associated IDPs via two-dimensional separation in the gas phase. Gas-phase IDP structures have been regarded as kinetically trapped states originated from conformational features in solution. However, an implication of the states remains elusive in the structural characterization of IDPs, because it is unclear what structural property of IDPs is preserved. Recent studies have indicated that the conformational features of IDPs in solution are not fully reproduced in the gas phase. Nevertheless, the molecular interactions captured in the gas phase amplify the structural differences between IDP conformers. Therefore, an IDP conformational change that is not observed in solution is observable in the gas-phase structures obtained by ESI-IM-MS. Herein, we have presented up-to-date researches on the key implications of kinetically trapped states in the gas phase with a brief summary of the structural dynamics of IDPs in ESI-IM-MS.

Determination of Gas-Phase Ion Structures of Locally Polar Homopolymers Through High-Resolution Ion Mobility Spectrometry-Mass Spectrometry

The strong synergy arising from coupling two orthogonal analytical techniques such as ion mobility and mass spectrometry can be used to separate complex mixtures and determine structural information of analytes in the gas phase. A tandem study is performed using two systems with different gases and pressures to ascertain gas-phase conformations of homopolymer ions. Aside from spherical and stretched configurations, intermediate configurations formed by a multiply charged globule and a "bead-on-a-string" appendix are confirmed for polyethylene-glycol (PEG), polycaprolactone (PCL), and polydimethylsiloxane (PDMS). These intermediate configurations are shown to be ubiquitous for all charge states and masses present. For each charge state, configurations evolve in two distinctive patterns: an inverse evolution which occurs as an elementary charge attached to the polymer leaves the larger globule and incorporates itself into the appendage, and a forward evolution which reduces the globule without relinquishing a charge while leaving the appendix relatively constant. Forward evolutions are confirmed to form self-similar family shapes that transcend charge states for all polymers. Identical structural changes occur at the same mass over charge regardless of the system, gas or pressure strongly suggesting that conformations are only contingent on number of charges and chain length, and start arranging once the ion is at least partially ejected from the droplet, supporting a charge extrusion mechanism. Configurational changes are smoother for PDMS which is attributed to the larger steric hindrance caused by protruding pendant groups. This study has implications in the study of the configurational space of more complex homopolymers and heteropolymers. Graphical Abstract.

Charging and supercharging of proteins for mass spectrometry: recent insights into the mechanisms of electrospray ionization

Molecular dynamics simulations have uncovered mechanistic details of the protein ESI process under various experimental conditions.

Strengths and Weaknesses of Molecular Simulations of Electrosprayed Droplets

The origin and the magnitude of the charge in a macroion are critical questions in mass spectrometry analysis coupled to electrospray and other ionization techniques that transfer analytes from the bulk solution into the gaseous phase via droplets. In many circumstances, it is the later stages of the existence of a macroion in the containing solvent drop before the detection that determines the final charge state. Experimental characterization of small (with linear dimensions of several nanometers) and short-lived droplets is quite challenging. Molecular simulations in principle may provide insight exactly in this challenging for experiments regime. We discuss the strengths and weaknesses of the molecular modeling of electrosprayed droplets using molecular dynamics. We illustrate the limitations of the molecular modeling in the analysis of large macroions and specifically proteins away from their native states. Graphical Abstract ᅟ.

When droplets become stars: charged dielectric droplets beyond the Rayleigh limit

When a nano-drop comprising a single spherical central ion and dielectric solvent is charged above a well-defined threshold, it acquires a stable star morphology. In contrast, conducting droplets, will undergo fission. Here we report combined atomistic molecular dynamics and continuum modelling study of star formation of droplets that contain a highly charged ion. We assume that in the continuum model the dielectric response is linear. In this linear continuum model, which is an extension of Rayleigh model, the energy of the drop is comprised of terms analogous to those in Rayleigh model, which are surface energy and electrostatic energy of dielectric droplet charged by a central point charge. We present the stability analysis of the continuum model to determine the threshold of instability. Indeed we find that the model accounts well for the onset of the instabilities. Molecular dynamics show that the number of points of the star-shaped nano-drops depends only on the surface tension, dielectric constant and size of the droplet, and on the magnitude of the charge of the central ion, but not on its sign. Intuitively, it is expected that when a spherical dielectric drop becomes unstable it would transform into a non-spherical finite shape of the same volume as the initial spherical shape with the point charge located in the drop interior. To test whether the extended Rayleigh model can account for the observed droplet shapes, we performed numerical simulations of the linear continuum model. Contrary to the expectations, the simulations of the extended Rayleigh model does not reproduce the stable star shapes found in the atomistic simulations, not even when we account for the bending rigidity and spontaneous curvature of the surface. We argue that the assumption that the dielectric response is linear breaks down if the droplet surface approaches the central macro-ion, where the electric field strength is such that dielectric saturation sets in. We envisage that for certain solvents, these stars could be made permanent by cross-linking, opening the way to the production of a novel class of highly-non-convex colloids.

Charging and Release Mechanisms of Flexible Macromolecules in Droplets

We study systematically the charging and release mechanisms of a flexible macromolecule, modeled by poly(ethylene glycol) (PEG), in a droplet by using molecular dynamics simulations. We compare how PEG is solvated and charged by sodium Na⁺ ions in a droplet of water (H₂O), acetonitrile (MeCN), and their mixtures. Initially, we examine the location and the conformation of the macromolecule in a droplet bearing no net charge. It is revealed that the presence of charge carriers do not affect the location of PEG in aqueous and MeCN droplets compared with that in the neutral droplets, but the location of the macromolecule and the droplet size do affect the PEG conformation. PEG is charged on the surface of a sodiated aqueous droplet that is found close to the Rayleigh limit. Its charging is coupled to the extrusion mechanism, where PEG segments leave the droplet once they coordinate a Na⁺ ion or in a correlated motion with Na⁺ ions. In contrast, as PEG resides in the interior of a MeCN droplet, it is sodiated inside the droplet. The compact macro-ion transitions through partially unwound states to an extended conformation, a process occurring during the final stage of desolvation and in the presence of only a handful of MeCN molecules. For charged H₂O/MeCN droplets, the sodiation of PEG is determined by the H₂O component, reflecting its slower evaporation and preference over MeCN for solvating Na⁺ ions. We use the simulation data to construct an analytical model that suggests that the droplet surface electric field may play a role in the macro-ion-droplet interactions that lead to the extrusion of the macro-ion. This study provides the first evidence of the effect of the surface electric field by using atomistic simulations. Graphical Abstract ᅟ.

Characterization of "Star" Droplet Morphologies Induced by Charged Macromolecules

"Star" morphologies of charged liquid droplets are distinct droplet conformations that, for a certain charge squared to volume ratio, have lower energy than their spherically shaped analogues. For these shapes to appear, the charge should be carried by a single ionic species. A typical example of a charge carrier that we employ in this study is a fully charged double-stranded oligodeoxynucleotide (dsDNA) in an aqueous and an acetonitrile droplet. We characterize the structure and dynamics of the star-shaped droplets. We find that by increasing the charge squared to volume ratio, the droplet evolves from spherical to "spiky" shapes, by first passing from droplet sizes that undergo enhanced shape fluctuations relative to those of the larger spherical droplets. These fluctuations mark the onset of the instability. We also find that in the spiky droplet, the orientation of the solvent molecules in the first shell about the dsDNA is very close to that in the bulk solution. However, this orientation is substantially different farther away from the dsDNA. With regards to dynamics, the motion of the spikes is reflected in the autocorrelation functions of rotationally invariant order parameters that show a damped oscillator form of decay, indicative of the elastic motion of the spikes. We compare the formation of spikes with that of the ferrofluids and the dielectric materials in an electric field, and we conclude that they represent a different entity that deserves its own characterization. The study provides insight into the manner in which the charge distribution may give rise to well-controlled droplet morphologies and calls for experiments in this direction.

Charge-induced instabilities of droplets containing macromolecular complexes

Solvated macromolecular complexes are ubiquitous in nature, notably in biological systems containing proteins and nucleic acids. Studies of the interactions within a macromolecular complex and between the complex and the solvent in droplet environments are critical for understanding the stability of macromolecular complexes in electrospray ionization (ESI) and nanofluidic experiments. In this study, two distinct cases of macromolecular complexes in aqueous nanodrops are examined by using molecular dynamics simulations: (i) a pair of sodiated poly(ethylene) glycol (PEG) macroions and (ii) a double-stranded DNA (dsDNA). PEG represents a case in which the surface energy of the aqueous droplet is larger than the solvent–macromolecule energy. Conversely, in a droplet solvating dsDNA, the solvent–macromolecule interaction energy overcomes the solvent interaction energy. We report that charge-induced instabilities previously identified for single macroions also appear in the case of complexes, but with a higher level of complexity. In the case of a pair of PEG macroions, we found that their conformations on the surface of a droplet “sense” each other. The charged PEGs are each released from a droplet at different times through contiguous extrusion or drying-out mechanisms. In the case of the DNA, the charge-induced instability manifests as a spine droplet morphology. Narrow regions of the spines promote break down of the hydrogen bonds that hold the dsDNA together. The dsDNA separates into two single strands as it is increasingly exposed to vacuum. These findings elucidate charge-induced instabilities of macromolecular complexes in droplets, which are critical intermediates in ESI and nanofluidic experiments.

How do non-covalent complexes dissociate in droplets? A case study of the desolvation of dsDNA from a charged aqueous nanodrop

We present the desolvation mechanism of a double-stranded oligodeoxynucleotide (dsDNA) from an aqueous nanodrop studied by using atomistic molecular dynamics methods.

Manifestations of charge induced instability in droplets effected by charged macromolecules

Ion-release processes from droplets that contain excess charge are of central importance in determining the charge-state distributions of macromolecules in electrospray ionization methods. We develop an analytical theory to describe the mechanism of contiguous extrusion of a charged macromolecule from a droplet. We find that the universal parameter determining the system behavior is the ratio of solvation energy per unit length to the square of the ion charge density per unit length. Systems with the same value of the ratio will follow the same path in the course of droplet evaporation. The analytical model is compared with molecular simulations of charged polyethylene glycol macroion in aqueous droplets, and the results are in excellent agreement.