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.

Mass Spectrometry of Nucleic Acid Noncovalent Complexes

Nucleic acids have been among the first targets for antitumor drugs and antibiotics. With the unveiling of new biological roles in regulation of gene expression, specific DNA and RNA structures have become very attractive targets, especially when the corresponding proteins are undruggable. Biophysical assays to assess target structure as well as ligand binding stoichiometry, affinity, specificity, and binding modes are part of the drug development process. Mass spectrometry offers unique advantages as a biophysical method owing to its ability to distinguish each stoichiometry present in a mixture. In addition, advanced mass spectrometry approaches (reactive probing, fragmentation techniques, ion mobility spectrometry, ion spectroscopy) provide more detailed information on the complexes. Here, we review the fundamentals of mass spectrometry and all its particularities when studying noncovalent nucleic acid structures, and then review what has been learned thanks to mass spectrometry on nucleic acid structures, self-assemblies (e.g., duplexes or G-quadruplexes), and their complexes with ligands.

Aptamer-ligand recognition studied by native ion mobility-mass spectrometry

The range of applications for aptamers, small oligonucleotide-based receptors binding to their targets with high specificity and affinity, has been steadily expanding. Our understanding of the mechanisms governing aptamer-ligand recognition and binding is however lagging, stymieing the progress in the rational design of new aptamers and optimization of the known ones. Here we demonstrate the capabilities and limitations of native ion mobility-mass spectrometry for the analysis of their higher-order structure and non-covalent interactions. A set of related cocaine-binding aptamers, displaying a range of folding properties and ligand binding affinities, was used as a case study in both positive and negative electrospray ionization modes. Using carefully controlled experimental conditions, we probed their conformational behavior and interactions with the high-affinity ligand quinine as a surrogate for cocaine. The ratios of bound and unbound aptamers in the mass spectra were used to rank them according to their apparent quinine-binding affinity, qualitatively matching the published ranking order. The arrival time differences between the free aptamer and aptamer-quinine complexes were consistent with a small ligand-induced conformational change, and found to inversely correlate with the affinity of binding. This mass spectrometry-based approach provides a fast and convenient way to study the molecular basis of aptamer-ligand recognition.

Stoichiometry of Nucleotide Binding to Proteasome AAA+ ATPase Hexamer Established by Native Mass Spectrometry

AAA+ ATPases constitute a large family of proteins that are involved in a plethora of cellular processes including DNA disassembly, protein degradation and protein complex disassembly. They typically form a hexametric ring-shaped structure with six subunits in a (pseudo) 6-fold symmetry. In a subset of AAA+ ATPases that facilitate protein unfolding and degradation, six subunits cooperate to translocate protein substrates through a central pore in the ring. The number and type of nucleotides in an AAA+ ATPase hexamer is inherently linked to the mechanism that underlies cooperation among subunits and couples ATP hydrolysis with substrate translocation. We conducted a native MS study of a monodispersed form of PAN, an archaeal proteasome AAA+ ATPase, to determine the number of nucleotides bound to each hexamer of the WT protein. We utilized ADP and its analogs (TNP-ADP and mant-ADP), and a nonhydrolyzable ATP analog (AMP-PNP) to study nucleotide site occupancy within the PAN hexamer in ADP- and ATP-binding states, respectively. Throughout all experiments we used a Walker A mutant (PANK217A) that is impaired in nucleotide binding as an internal standard to mitigate the effects of residual solvation on mass measurement accuracy and to serve as a reference protein to control for nonspecific nucleotide binding. This approach led to the unambiguous finding that a WT PAN hexamer carried - from expression host - six tightly bound ADP molecules that could be exchanged for ADP and ATP analogs. Although the Walker A mutant did not bind ADP analogs, it did bind AMP-PNP, albeit at multiple stoichiometries. We observed variable levels of hexamer dissociation and an appearance of multimeric species with the over-charged molecular ion distributions across repeated experiments. We posit that these phenomena originated during ESI process at the final stages of ESI droplet evolution.

Ion Mobility Mass Spectrometry Uncovers the Impact of the Patterning of Oppositely Charged Residues on the Conformational Distributions of Intrinsically Disordered Proteins

The global dimensions and amplitudes of conformational fluctuations of intrinsically disordered proteins are governed, in part, by the linear segregation versus clustering of oppositely charged residues within the primary sequence. Ion mobility-mass spectrometry (IM-MS) affords unique advantages for probing the conformational consequences of the linear patterning of oppositely charged residues because it measures and separates proteins electrosprayed from solution on the basis of charge and shape. Here, we use IM-MS to measure the conformational consequences of charge patterning on the C-terminal intrinsically disordered region (p27 IDR) of the cell cycle inhibitory protein p27^Kip1. We report the range of charge states and accompanying collisional cross section distributions for wild-type p27 IDR and two variants with identical amino acid compositions, κ14 and κ56, distinguished by the extent of linear mixing versus segregation of oppositely charged residues. Wild-type p27 IDR (κ31) and κ14, where the oppositely charged residues are more evenly distributed, exhibit a broad distribution of charge states. This is concordant with high degrees of conformational heterogeneity in solution. By contrast, κ56 with linear segregation of oppositely charged residues leads to limited conformational heterogeneity and a narrow distribution of charged states. Gas-phase molecular dynamics simulations demonstrate that the interplay between chain solvation and intrachain interactions (self-solvation) leads to conformational distributions that are modulated by salt concentration, with the wild-type sequence showing the most sensitivity to changes in salt concentration. These results suggest that the charge patterning within the wild-type p27 IDR may be optimized to sample both highly solvated and self-solvated conformational states.

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 ᅟ.

What factors determine the stability of a weak protein-protein interaction in a charged aqueous droplet?

Maintaining the interface of a weak transient protein complex transferred from bulk solution to the gaseous state via evaporating droplets is a critical question in the detection of the complex association (dissociation) constant by using electrospray ionization mass spectrometry (ESI-MS). Here we explore the factors that may affect the stability of a protein-protein interaction (PPI) using atomistic molecular dynamics (MD) modelling of a complex of ubiquitin (Ub) and the ubiquitin-associated domain (UbA) (RCSB PDB code ) and a non-covalent complex of diubiquitin (RCSB PDB code ) in aqueous droplets. A general method is presented to determine the protonation states of the complexes we investigate in particular, and that of a protein in general, under various pH conditions that an evaporating droplet acquires due to its change in size. We find that the combination of high temperature and high charge states of the protein complexes may destabilize the interface by creating new interfaces instead of a direct rupture of the initial stable interface. We provide evidence that highly charged protein complexes are found in droplets that form conical extrusions of the solvent on the surface due to charge-induced instability. This distinct droplet morphology leads to a higher solvent evaporation rate that assists in transferring the complex in the gaseous state without dissociation. The conical solvent protrusions expose on the droplet surface certain amino acids that otherwise would be solvated in a droplet with the protein complex of low charge states. The new vapor-protein interface does not have a direct effect on the stability of the PPI. A common way in experiments to stabilize the protein complexes in droplets is to reduce the protonation state of the proteins. Here we find that weakly bound protein complexes even at high protonation states can be stabilized by the presence of a small number of counterions, without affecting the protonation state of the protein. Our findings may provide guiding principles in ESI-MS experiments to stabilize weak transient PPIs.

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 ᅟ.

The Evolution of Electrospray Generated Droplets is Not Affected by Ionization Mode

Ionization efficiency and mechanism in ESI is strongly affected by the properties of mobile phase. The use of mobile-phase properties to accurately describe droplets in ESI source is convenient but may be inadequate as the composition of the droplets is changing in the plume due to electrochemical reactions occurring in the needle tip as well as continuous drying and fission of droplets. Presently, there is paucity of research on the effect of the polarity of the ESI mode on mobile phase composition in the droplets. In this paper, the change in the organic solvent content, pH, and droplet size are studied in the ESI plume in both ESI+ and ESI- ionization mode. We introduce a rigorous way - the absolute pH (pH_abs^H₂^O) - to describe pH change in the plume that takes into account organic solvent content in the mobile phase. pH_abs^H₂^O enables comparing acidities of ESI droplets with different organic solvent contents. The results are surprisingly similar for both ionization modes, indicating that the dynamics of the change of mobile-phase properties is independent from the ESI mode used. This allows us to conclude that the evolution of ESI droplets first of all proceeds via the evaporation of the organic modifier and to a lesser extent via fission of smaller droplets from parent droplets. Secondly, our study shows that qualitative findings related to the ESI process obtained on the ESI+ mode can almost directly be applied also in the ESI- mode. Graphical Abstract ᅟ.

Compaction of Duplex Nucleic Acids upon Native Electrospray Mass Spectrometry

We report on the fate of nucleic acids conformation in the gas phase as sampled using native mass spectrometry coupled to ion mobility spectrometry. On the basis of several successful reports for proteins and their complexes, the technique has become popular in structural biology, and the conformation survival becomes more and more taken for granted. Surprisingly, we found that DNA and RNA duplexes, at the electrospray charge states naturally obtained from native solution conditions (≥100 mM aqueous NH₄OAc), are significantly more compact in the gas phase compared to the canonical solution structures. The compaction is observed for all duplex sizes (gas-phase structures are more compact than canonical B-helices by ∼20% for 12-bp, and by up to ∼30% for 36-bp duplexes), and for DNA and RNA alike. Molecular modeling (density functional calculations on small helices, semiempirical calculations on up to 12-bp, and molecular dynamics on up to 36-bp duplexes) demonstrates that the compaction is due to phosphate group self-solvation prevailing over Coulomb repulsion. Molecular dynamics simulations starting from solution structures do not reproduce the experimental compaction. To be experimentally relevant, molecular dynamics sampling should reflect the progressive structural rearrangements occurring during desolvation. For nucleic acid duplexes, the compaction observed for low charge states results from novel phosphate-phosphate hydrogen bonds formed across both grooves at the very late stages of electrospray.

Stability of a Transient Protein Complex in a Charged Aqueous Droplet with Variable pH

Electrospray ionization mass spectrometry (ESI-MS) has the potential to become a high-throughput robust experimental method for the detection of protein-protein equilibrium constants. Poorly understood processes that affect the stability of weak noncovalent protein complexes in the intervening droplet environment are a significant factor that precludes the advancement of the method. We use molecular dynamics to study the stability of a ubiquitin and ubiquitin-associated domain complex (RCSB PDB code 2MRO ) in an aqueous droplet with changing size and charge concentration. We present evidence that a weak protein complex changes conformation and may dissociate in shrinking droplets. Then, the droplets containing these dissociated proteins divide. Our findings suggest that in some cases ESI-MS does not measure the correct association constants. The study intends to stimulate research for systematic development of experimental protocols that stabilize weakly bound protein interfaces in droplets.

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.