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

Effects of Nano-Electrospray Ionization Emitter Position on Unintentional In-Source Activation of Peptide and Protein Ions

Native ion mobility-mass spectrometry (IM-MS) typically introduces protein ions into the gas phase through nano-electrospray ionization (nESI). Many nESI setups have mobile stages for tuning the ion signal and extent of co-solute and salt adduction. However, tuning the position of the emitter capillary in nESI can have unintended downstream consequences for collision-induced unfolding or collision-induced dissociation (CIU/D) experiments. Here, we show that relatively small variations in the nESI emitter position can shift the midpoint (commonly called the "CID50" or "CIU50") potential of CID breakdown curves and CIU transitions by as much as 8 V on commercial instruments. A spatial "map" of the shift in CID50 for the loss of heme from holomyoglobin onto the emitter position on a Waters Synapt G2-Si mass spectrometer shows that emitter positions closer to the instrument inlet can result in significantly greater in-source activation, whereas different effects are found on an Agilent 6545XT instrument for the ions studied. A similar effect is observed for CID of the singly protonated leucine enkephalin peptide and Shiga toxin 1 subunit B homopentamer on the Waters Synapt G2-Si instrument. In-source activation effects on a Waters Synapt G2-Si are also investigated by examining the RMSD between CIU fingerprints acquired at different emitter positions and the shifts in CIU50 for structural transitions of bovine serum albumin and NIST monoclonal antibody.

Atomistic Insights into the Formation of Nonspecific Protein Complexes during Electrospray Ionization

Native electrospray ionization (ESI)-mass spectrometry (MS) is widely used for the detection and characterization of multi-protein complexes. A well-known problem with this approach is the possible occurrence of nonspecific protein clustering in the ESI plume. This effect can distort the results of binding affinity measurements, and it can even generate gas-phase complexes from proteins that are strictly monomeric in bulk solution. By combining experiments and molecular dynamics (MD) simulations, the current work for the first time provides detailed insights into the ESI clustering of proteins. Using ubiquitin as a model system, we demonstrate how the entrapment of more than one protein molecule in an ESI droplet can generate nonspecific clusters (e.g., dimers or trimers) via solvent evaporation to dryness. These events are in line with earlier proposals, according to which protein clustering is associated with the charged residue model (CRM). MD simulations on cytochrome c (which carries a large intrinsic positive charge) confirmed the viability of this CRM avenue. In addition, the cytochrome c data uncovered an alternative mechanism where protein-protein contacts were formed early within ESI droplets, followed by cluster ejection from the droplet surface. This second pathway is consistent with the ion evaporation model (IEM). The observation of these IEM events for large protein clusters is unexpected because the IEM has been thought to be associated primarily with low-molecular-weight analytes. In all cases, our MD simulations produced protein clusters that were stabilized by intermolecular salt bridges. The MD-generated charge states agreed with experiments. Overall, this work reveals that ESI-induced protein clustering does not follow a tightly orchestrated pathway but can proceed along different avenues.

Computationally-obtained structural insights into the molecular interactions between Pidilizumab and binding partners DLL1 and PD-1

Pidilizumab is a monoclonal antibody tested against several types of malignancies, such as lymphoma and metastatic melanoma, showing promising results. In 2016, the FDA put Pidilizumab's clinical studies on partial hold due to emerging evidence pointing to the antibody target uncertainty. Although initial studies indicated an interaction with the PD-1 checkpoint receptor, recent updates assert that Pidilizumab binds primarily to Notch ligand DLL1. However, a detailed description of which interactions coordinate antibody-antigen complex formation is lacking. Therefore, this study uses computational tools to identify molecular interactions between Pidilizumab and its reported targets PD-1 and DLL1. A docking methodology was validated and applied to determine the binding modes between modeled Pidilizumab scFvs and the two antigens. We used Molecular Dynamics (MD) simulations to verify the complexes' stability and submitted the resulting trajectory files to MM/PBSA and Principal Component Analysis. A set of different prediction tools determined scFv interface hot-spots. Whereas docking and MD simulations revealed that the antibody fragments do not interact straightforwardly with PD-1, ten scFv hot-spots, including Met93 and Leu112, mediated the interaction with the DLL1 C2 domain. The interaction triggered a conformational selection-like effect on DLL1, allowing new hydrogen bonds on the β3-β4 interface loop. The unprecedented structural data on Pidilizumab's interactions provided novel evidence that its legitimate target is the DLL1 protein and offered structural insight on how these molecules interact, shedding light on the pathways that could be affected by the use of this essential immunobiological.Communicated by Ramaswamy H. Sarma.

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.

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

Addressing a Common Misconception: Ammonium Acetate as Neutral pH "Buffer" for Native Electrospray Mass Spectrometry

Native ESI-MS involves the transfer of intact proteins and biomolecular complexes from solution into the gas phase. One potential pitfall is the occurrence of pH-induced changes that can affect the analyte while it is still surrounded by solvent. Most native ESI-MS studies employ neutral aqueous ammonium acetate solutions. It is a widely perpetuated misconception that ammonium acetate buffers the analyte solution at neutral pH. By definition, a buffer consists of a weak acid and its conjugate weak base. The buffering range covers the weak acid pK_a ± 1 pH unit. NH₄⁺ and CH₃-COO^- are not a conjugate acid/base pair, which means that they do not constitute a buffer at pH 7. Dissolution of ammonium acetate salt in water results in pH 7, but this pH is highly labile. Ammonium acetate does provide buffering around pH 4.75 (the pK_a of acetic acid) and around pH 9.25 (the pK_a of ammonium). This implies that neutral ammonium acetate solutions electrosprayed in positive ion mode will likely undergo acidification down to pH 4.75 ± 1 in the ESI plume. Ammonium acetate nonetheless remains a useful additive for native ESI-MS. It is a volatile electrolyte that can mimic the solvation properties experienced by proteins under physiological conditions. Also, a drop from pH 7 to around pH 4.75 is less dramatic than the acidification that would take place in pure water. It is hoped that the habit of referring to pH 7 solutions as ammonium acetate "buffer" will disappear from the literature. Ammonium acetate "solution" should be used instead. Graphical Abstract ᅟ.