Molecular simulation-based structural prediction of protein complexes in mass spectrometry: the human insulin dimer

Protonation State-Induced Unfolding of Protein Secondary Structure in the Gas Phase

The combination of infrared spectroscopy (IR) and ion mobility mass spectrometry (IM-MS) has revealed that protein secondary structures are retained upon transformation from aqueous solution to the gas phase under gentle conditions. Yet the details about where and how these structural elements are embedded in the gas phase remain elusive. In this study, we employ long time scale molecular dynamics (MD) simulations to examine the extent to which proteins retain their solution structures and the impact of protonation state on the stability of secondary structures in the gas phase. Our investigation focuses on two well-studied proteins, myoglobin and β-lactoglobulin, representing typical helical and β-sheet proteins, respectively. Our simulations accurately reproduce the experimental collision cross section (CCS) data measured by IM-MS. Based on accurately reproducing previous experimental collision cross section data and dominant secondary structural species obtained from IM-MS and IR, we confirm that both proteins largely retain their native secondary structural components upon passing from aqueous solution to the gas phase. However, we observe significant reductions in secondary structure contents (19.2 ± 1.2% for myoglobin and 7.3 ± 0.6% for β-lactoglobulin) in specific regions predominantly composed of ionizable residues. Further mechanistic analysis suggests that alterations in protonation states of these residues after phase transition induce changes in their local interaction networks and backbone dihedral angles, which potentially promote the unfolding of secondary structures in the gas phase. We anticipate that similar protonation state induced unfolding may be observed in other proteins possessing distinct secondary structures. Further studies on a broader array of proteins will be essential to refine our understanding of protein structural behavior during the transition to the gas phase.

High-Performance Molecular Dynamics Simulations for Native Mass Spectrometry of Large Protein Complexes with the Fast Multipole Method

Native mass spectrometry (MS) is widely employed to study the structures and assemblies of proteins ranging from small monomers to megadalton complexes. Molecular dynamics (MD) simulation is a useful complement as it provides the spatial detail that native MS cannot offer. However, MD simulations performed in the gas phase have suffered from rapidly increasing computational costs with the system size. The primary bottleneck is the calculation of electrostatic forces, which are effective over long distances and must be explicitly computed for each atom pair, precluding efficient use of methods traditionally used to accelerate condensed-phase simulations. As a result, MD simulations have been unable to match the capacity of MS in probing large multimeric protein complexes. Here, we apply the fast multipole method (FMM) for computing the electrostatic forces, recently implemented by Kohnke et al. (J. Chem. Theory Comput., 2020, 16, 6938-6949), showing that it significantly enhances the performance of gas-phase simulations of large proteins. We assess how to achieve adequate accuracy and optimal performance with FMM, finding that it expands the accessible size range and time scales dramatically. Additionally, we simulate a 460 kDa ferritin complex over microsecond time scales, alongside complementary ion mobility (IM)-MS experiments, uncovering conformational changes that are not apparent from the IM-MS data alone.

Identification and Tetramer Structure of Hemin-Binding Protein SPD_0310 Linked to Iron Homeostasis and Virulence of Streptococcus pneumoniae

Iron and iron-containing compounds are essential for bacterial virulence and host infection. Hemin is an important supplement compound for bacterial survival in an iron-deficient environment. Despite strong interest in hemin metabolism, the detailed mechanism of hemin transportation in Gram-positive bacteria is yet to be reported. The results of our study revealed that the homologous proteins of SPD_0310 were significantly conservative in Gram-positive bacteria (P < 0.001), and these proteins were identified as belonging to an uncharacterized protein family (UPF0371). The results of thermodynamic and kinetic studies have shown that SPD_0310 has a high hemin-binding affinity. Interestingly, we found that the crystal structure of SPD_0310 presented a homotetramer conformation, which is required for hemin binding. SPD_0310 can interact with many hemin-binding proteins (SPD_0090, SPD_1609, and GAPDH) located on the cell surface, which contributes to hemin transfer to the cytoplasm. It also has a high affinity with other iron transporters in the cytoplasm (SPD_0226 and SPD_0227), which facilitates iron redistribution in cells. More importantly, the knockout of the spd_0310 gene (Δspd_0310) resulted in a decrease in the iron content and protein expression levels of many bacterial adhesion factors. Moreover, the animal model showed that the Δspd_0310 strain has a lower virulence than the wild type. Based on the crystallographic and biochemical studies, we inferred that SPD_0310 is a hemin intermediate transporter which contributes to iron homeostasis and further affects the virulence of Streptococcus pneumoniae in the host. Our study provides not only an important theoretical basis for the in-depth elucidation of the hemin transport mechanism in bacteria but also an important candidate target for the development of novel antimicrobial agents based on metal transport systems. IMPORTANCE Iron is an essential element for bacterial virulence and infection of the host. The detailed hemin metabolism in Gram-positive bacteria has rarely been studied. SPD_0310 belongs to the UPF0371 family of proteins, and results of homology analysis and evolutionary tree analysis suggested that it was widely distributed and highly conserved in Gram-positive bacteria. However, the function of the UPF0371 family remains unknown. We successfully determined the crystal structure of apo-SPD_0310, which is a homotetramer. We found that cytoplasmic protein SPD_0310 with a special tetramer structure has a strong hemin-binding ability and interacts with many iron transporters, which facilitates hemin transfer from the extracellular space to the cytoplasm. The results of detailed functional analyses indicated that SPD_0310 may function as a hemin transporter similar to hemoglobin in animals and contributes to bacterial iron homeostasis and virulence. This study provides a novel target for the development of antimicrobial drugs against pathogenic Gram-positive bacteria.

Mass Spectrometric and Bio-Computational Binding Strength Analysis of Multiply Charged RNAse S Gas-Phase Complexes Obtained by Electrospray Ionization from Varying In-Solution Equilibrium Conditions

We investigated the influence of a solvent’s composition on the stability of desorbed and multiply charged RNAse S ions by analyzing the non-covalent complex’s gas-phase dissociation processes. RNAse S was dissolved in electrospray ionization-compatible buffers with either increasing organic co-solvent content or different pHs. The direct transition of all the ions and the evaporation of the solvent from all the in-solution components of RNAse S under the respective in-solution conditions by electrospray ionization was followed by a collision-induced dissociation of the surviving non-covalent RNAse S complex ions. Both types of changes of solvent conditions yielded in mass spectrometrically observable differences of the in-solution complexation equilibria. Through quantitative analysis of the dissociation products, i.e., from normalized ion abundances of RNAse S, S-protein, and S-peptide, the apparent kinetic and apparent thermodynamic gas-phase complex properties were deduced. From the experimental data, it is concluded that the stability of RNAse S in the gas phase is independent of its in-solution equilibrium but is sensitive to the complexes’ gas-phase charge states. Bio-computational in-silico studies showed that after desolvation and ionization by electrospray, the remaining binding forces kept the S-peptide and S-protein together in the gas phase predominantly by polar interactions, which indirectly stabilized the in-bulk solution predominating non-polar intermolecular interactions. As polar interactions are sensitive to in-solution protonation, bio-computational results provide an explanation of quantitative experimental data with single amino acid residue resolution.

SUMOylation of spastin promotes the internalization of GluA1 and regulates dendritic spine morphology by targeting microtubule dynamics

Dendritic spines are specialized structures involved in neuronal processes on which excitatory synaptic contact occurs. The microtubule cytoskeleton is vital for maintaining spine morphology and mature synapses. Spastin is related to microtubule-severing proteases and is involved in synaptic bouton formation. However, it is not yet known if spastin can be modified by Small Ubiquitin-like Modifier (SUMO) or how this modification regulates dendritic spines. Spastin was shown to be SUMOylated at K427, and its deSUMOylation promoted microtubule stability. In addition, SUMOylation of spastin was shown to affect signalling pathways associated with long term synaptic depression. SUMOylated spastin promoted the development of dendrites and dendritic spines. Moreover, SUMOylated spastin regulated endocytosis and affected the transport of the AMPA receptor, GluA1. Our findings suggest that SUMOylation of spastin promotes GluA1 internalization and regulates dendritic spine morphology through targeting of microtubule dynamics.

Role of hydrophobic residues for the gaseous formation of helical motifs

The secondary structure content of proteins and their complexes may change significantly on passing from aqueous solution to the gas phase (as in mass spectrometry experiments). In this work, we investigate the impact of hydrophobic residues on the formation of the secondary structure of a real protein complex in the gas phase. We focus on a well-studied protein complex, the amyloid-β (1-40) dimer (2Aβ). Molecular dynamics simulations reproduce the results of ion mobility-mass spectrometry experiments. In addition, a helix (not present in the solution) is identified involving 19FFAED23, consistent with infrared spectroscopy data on an Aβ segment. Our simulations further point to the role of hydrophobic residues in the formation of helical motifs - hydrophobic sidechains "shield" helices from being approached by residues that carry hydrogen bond sites. In particular, two hydrophobic phenylalanine residues, F19 and F20, play an important role for the helix, which is induced in the gas phase in spite of the presence of two carboxyl-containing residues.

Two zinc-binding domains in the transporter AdcA from facilitate high-affinity binding and fast transport of zinc

Zinc is an essential metal in bacteria. One important bacterial zinc transporter is AdcA, and most bacteria possess AdcA homologs that are single-domain small proteins due to better efficiency of protein biogenesis. However, a double-domain AdcA with two zinc-binding sites is significantly overrepresented in Streptococcus species, many of which are major human pathogens. Using molecular simulation and experimental validations of AdcA from Streptococcus pyogenes, we found here that the two AdcA domains sequentially stabilize the structure upon zinc binding, indicating an organization required for both increased zinc affinity and transfer speed. This structural organization appears to endow Streptococcus species with distinct advantages in zinc-depleted environments, which would not be achieved by each single AdcA domain alone. This enhanced zinc transport mechanism sheds light on the significance of the evolution of the AdcA domain fusion, provides new insights into double-domain transporter proteins with two binding sites for the same ion, and indicates a potential target of antimicrobial drugs against pathogenic Streptococcus species.

Proton Dynamics in Protein Mass Spectrometry

Native electrospray ionization/ion mobility-mass spectrometry (ESI/IM-MS) allows an accurate determination of low-resolution structural features of proteins. Yet, the presence of proton dynamics, observed already by us for DNA in the gas phase, and its impact on protein structural determinants, have not been investigated so far. Here, we address this issue by a multistep simulation strategy on a pharmacologically relevant peptide, the N-terminal residues of amyloid-β peptide (Aβ(1-16)). Our calculations reproduce the experimental maximum charge state from ESI-MS and are also in fair agreement with collision cross section (CCS) data measured here by ESI/IM-MS. Although the main structural features are preserved, subtle conformational changes do take place in the first ∼0.1 ms of dynamics. In addition, intramolecular proton dynamics processes occur on the picosecond-time scale in the gas phase as emerging from quantum mechanics/molecular mechanics (QM/MM) simulations at the B3LYP level of theory. We conclude that proton transfer phenomena do occur frequently during fly time in ESI-MS experiments (typically on the millisecond time scale). However, the structural changes associated with the process do not significantly affect the structural determinants.

Conformational effects in protein electrospray-ionization mass spectrometry

Electrospray-ionization mass spectrometry (ESI-MS) is a key tool of structural biology, complementing the information delivered by conventional biochemical and biophysical methods. Yet, the mechanism behind the conformational effects in protein ESI-MS is an object of debate. Two parameters-solvent-accessible surface area (As) and apparent gas-phase basicity (GBapp)-are thought to play a role in controlling the extent of protein ionization during ESI-MS experiments. This review focuses on recent experimental and theoretical investigations concerning the influence of these parameters on ESI-MS results and the structural information that can be derived. The available evidence supports a unified model for the ionization mechanism of folded and unfolded proteins. These data indicate that charge-state distribution (CSD) analysis can provide valuable structural information on normally folded, as well as disordered structures.