Gas-phase Conformations Research Articles

Structural Elucidation of Ubiquitin via Gas-Phase Ion/Ion Cross-Linking Reactions Using Sodium-Cationized Reagents Coupled with Infrared Multiphoton Dissociation.

Accurate structural determination of proteins is critical to understanding their biological functions and the impact of structural disruption on disease progression. Gas-phase cross-linking mass spectrometry (XL-MS) via ion/ion reactions between multiply charged protein cations and singly charged cross-linker anions has previously been developed to obtain low-resolution structural information on proteins. This method significantly shortens experimental time relative to conventional solution-phase XL-MS but has several technical limitations: (1) the singly deprotonated N-hydroxysulfosuccinimide (sulfo-NHS)-based cross-linker anions are restricted to attachment at neutral amine groups of basic amino acid residues and (2) analyzing terminal cross-linked fragment ions is insufficient to unambiguously localize sites of linker attachment. Herein, we demonstrate enhanced structural information for alcohol-denatured A-state ubiquitin obtained from an alternative gas-phase XL-MS approach. Briefly, singly sodiated ethylene glycol bis(sulfosuccinimidyl succinate) (sulfo-EGS) cross-linker anions enable covalent cross-linking at both ammonium and amine groups. Additionally, covalently modified internal fragment ions, along with terminal b-/y-type counterparts, improve the determination of linker attachment sites. Molecular dynamics simulations validate experimentally obtained gas-phase conformations of denatured ubiquitin. This method has identified four cross-linking sites across 8+ ubiquitin, including two new sites in the N-terminal region of the protein that were originally inaccessible in prior gas-phase XL approaches. The two N-terminal cross-linking sites suggest that the N-terminal half of ubiquitin is more compact in gas-phase conformations. By comparison, the two C-terminal linker sites indicate the signature transformation of this region of the protein from a native to a denatured conformation. Overall, the results suggest that the solution-phase secondary structures of the A-state ubiquitin are conserved in the gas phase. This method also provides sufficient sensitivity to differentiate between two gas-phase conformers of the same charge state with subtle structural variations.

Molecular Gas-Phase Conformational Ensembles.

Accurately determining the global minima of a molecular structure is important in diverse scientific fields, including drug design, materials science, and chemical synthesis. Conformational search engines serve as valuable tools for exploring the extensive conformational space of molecules and for identifying energetically favorable conformations. In this study, we present a comparison of Auto3D, CREST, Balloon, and ETKDG (from RDKit), which are freely available conformational search engines, to evaluate their effectiveness in locating global minima. These engines employ distinct methodologies, including machine learning (ML) potential-based, semiempirical, and force field-based approaches. To validate these methods, we propose the use of collisional cross-section (CCS) values obtained from ion mobility-mass spectrometry studies. We hypothesize that experimental gas-phase CCS values can provide experimental evidence that we likely have the global minimum for a given molecule. To facilitate this effort, we used our gas-phase conformation library (GPCL) which currently consists of the full ensembles of 20 small molecules and can be used by the community to validate any conformational search engine. Further members of the GPCL can be readily created for any molecule of interest using our standard workflow used to compute CCS values, expanding the ability of the GPCL in validation exercises. These innovative validation techniques enhance our understanding of the conformational landscape and provide valuable insights into the performance of conformational generation engines. Our findings shed light on the strengths and limitations of each search engine, enabling informed decisions for their utilization in various scientific fields, where accurate molecular structure determination is crucial for understanding biological activity and designing targeted interventions. By facilitating the identification of reliable conformations, this study significantly contributes to enhancing the efficiency and accuracy of molecular structure determination, with particular focus on metabolite structure elucidation. The findings of this research also provide valuable insights for developing effective workflows for predicting the structures of unknown compounds with high precision.

The role of solvation on the conformational landscape of α-synuclein.

Native ion mobility mass spectrometry has been used extensively to characterize ensembles of intrinsically disordered protein (IDP) conformers, but the extent to which the gaseous measurements provide realistic pictures of the solution conformations for such flexible proteins remains unclear. Therefore, we systematically studied the relationship between the solution and gaseous structural ensembles by measuring electrospray charge state and collision cross section (CCS) distributions for cationic and anionic forms of α-synuclein (αSN), an anionic protein in solution, as well as directly probed gas phase residue to residue distances via ion/ion reactions between gaseous α-synuclein cations and disulfonic acid linkers that form strong electrostatic bonds. We also combined results from in-solution protein crosslinking identified from native tandem mass spectrometry (MS/MS) with an initial αSN ensemble generated computationally by IDPConformerGenerator to generate an experimentally restrained solution ensemble of αSN. CCS distributions were directly calculated for the solution ensembles determined by NMR and compared to predicted gaseous conformers. While charge state and collision cross section distributions are useful for qualitatively describing the relative structural dynamics of proteins and major conformational changes induced by changes to solution states, the predicted and measured gas phase conformers include subpopulations that are significantly different than those expected from completely "freezing" solution conformations and preserving them in the gas phase. However, insights were gained on the various roles of solvent in stabilizing various conformers for extremely dynamic proteins like α-synuclein.

Gas Phase Conformation of Trisaccharides and Core Pentasaccharide: A Three-Step Tree-Based Sampling and Quantum Mechanical Computational Approach.

As an important component of N-linked glycoproteins, the core pentasaccharide is highly crucial to the potential application prospect of glycoprotein. However, the gas phase conformation study is a challenging one due to the size and complexity of the molecule, together with the necessity to rely on quantum chemistry modeling for relevant energetics and structures. In this paper, the structures of the trisaccharides and core pentasaccharides in N-linked glycans in the gas phase were constructed by a three-step tree-based (TSTB) sampling. Since single point energies of all the conformers are calculated at the temperature of zero, it is necessary to evaluate the stability at a high temperature. We calculate the Gibbs free energies using the standard thermochemistry model (T = 298.15 K). For trimannose, the energetic ordering at 298.15 K can be strongly changed compared to 0 K. Moreover, two structures of trimannose with high energies at 0 K are considered to provide a much better match of IR vibration signatures with the low Gibbs free energies. On this basis, the core pentasaccharide was constructed in three ways. The building configurations of core pentasaccharide were optimized to obtain reasonable low-energy stable conformers. Fortunately, the lowest-energy structure of core pentasaccharide is eventually the minimum at 0 K and 298.15 K. Furthermore, spectrum analysis of core pentasaccharide was carried out. Although poorly resolved, its contour from the experiment was in qualitative correspondence with the computed IR spectrum associated with its minimum free energy structure. A large number of strongly and weakly hydrogen-bonded hydroxyl and acetylamino groups contribute to a highly congested set of overlapping bands. Compared with traditional conformation generators, the TSTB sampling is employed to efficiently and comprehensively obtain preferred conformers of larger saccharides with lower energy.

An NMR-based approach for the design of crystallization processes for rigid molecules and effects of molecular rigidity and impurities on crystallization pathways

A methodology based on NMR analysis to enable development of crystallization procedures for relatively rigid molecules is presented. The number of relatively stable conformers in solution and their relative proportions are determined from peak splitting of equivalent protons and their relative areas as measured by NMR. Experimental data indicated that slow conformational interconversion was limiting crystallization at low temperatures. Variable temperature NMR analysis was utilized to determine rotamers’ peak coalescence and identify a suitable crystallization solvent. This NMR-based approach led to the development of a manufacturing crystallization recipe of a relatively rigid molecule that was impossible to crystallize via a scalable process. Furthermore, prediction of relatively stable conformer in gas phase, determined by Molecular Dynamics (MD) calculations were in acceptable agreement with data obtained from NMR analyses. Finally, the effects of molecular rigidity and impurities on crystallization and oiling-out pathways are described and an approach to avoid oiling-out is presented.

Systematic Comparison of Experimental Crystallographic Geometries and Gas-Phase Computed Conformers for Torsion Preferences.

We performed exhaustive torsion sampling on more than 3 million compounds using the GFN2-xTB method and performed a comparison of experimental crystallographic and gas-phase conformers. Many conformer sampling methods derive torsional angle distributions from experimental crystallographic data, limiting the torsion preferences to molecules that must be stable, synthetically accessible, and able to be crystallized. In this work, we evaluate the differences in torsional preferences of experimental crystallographic geometries and gas-phase computed conformers from a broad selection of compounds to determine whether torsional angle distributions obtained from semiempirical methods are suitable priors for conformer sampling. We find that differences in torsion preferences can be mostly attributed to a lack of available experimental crystallographic data with small deviations derived from gas-phase geometry differences. GFN2 demonstrates the ability to provide accurate and reliable torsional preferences that can provide a basis for new methods free from the limitations of experimental data collection. We provide Gaussian-based fits and sampling distributions suitable for torsion sampling and propose an alternative to the widely used "experimental torsion and knowledge distance geometry" (ETKDG) method using quantum torsion-derived distance geometry (QTDG) methods.

Fingerprinting of Underivatized Monosaccharide Stereoisomers Using High-Resolution Ion Mobility Spectrometry and Its Implications for Carbohydrate Sequencing.

Although carbohydrates are the most abundant biopolymers on Earth, there is currently no streamlined method to elucidate their complete sequence. Mass spectrometry (MS) alone is blind to many cases of isomerism and thus gives incomplete information for carbohydrates. Notably, the coexistence of numerous stereoisomeric monosaccharide subunits is of special concern. Over the last 10 years, the coupling of ion mobility spectrometry (IMS) with MS has kept gaining momentum─especially with the advent of high-resolution (HR) IMS devices such as cyclic IMS (cIMS). In fact, IMS is sensitive to the gas-phase conformations of molecules and, thus, to stereoisomerisms. In this article, we present innovative ion mobility methods on a cIMS instrument that allowed us to build a database of HR-IMS fingerprints for various underivatized monosaccharide stereoisomers. The conditions were fully compatible with MS/MS fragmentation approaches. We further verify that these fingerprints afford the identification of monosaccharidic fragments released upon collisional fragmentation of oligosaccharides. Overall, these results pave the way toward direct sequencing of carbohydrates at the monosaccharide level using HR-IMS.

Gas-Phase Behavior of Galactofuranosides upon Collisional Fragmentation: A Multistage High-Resolution Ion Mobility Study

Carbohydrates are ubiquitous in nature but are among the least conserved biomolecules in life. These biopolymers pose a particular challenge to analytical chemists because of their high diversity and structural heterogeneity. In addition, they contain many isomerisms that complicate their structural characterization, notably by mass spectrometry. The tautomerism of the constitutive subunits is of particular interest. A given cyclized monosaccharide unit can take two forms: a most common 6-membered ring (pyranose, p) and a more flexible 5-membered ring (furanose, f). The tautomers impact the biological properties of polysaccharides, resulting in interesting properties of the derived oligosaccharides. From an analytical point of view, the impact of tautomerism on the gas-phase behavior of ions has scarcely been described in the literature. In this work, we study the behavior of Galf-containing oligosaccharides, ionized as [M+Li]+ species, under collisional dissociation (CID) conditions using high-resolution and multistage ion mobility (IMS) on a Cyclic IMS platform. In the first part of this work, we studied whether disaccharidic fragments released from Galf-containing (Gal)1(Man)2 trisaccharides (and their Galp counterpart) would match the corresponding disaccharide standards, and─despite the fragments generally being a good match─we showed the possibility of Galf migrations and other unidentified alterations in the IMS profile. Next, we expanded on these unknown features using multistage IMS and molecular dynamics, unveiling the contributions of additional gas-phase conformers in the profile of fragments from a Galf-containing trisaccharide compared with the corresponding disaccharides.

The Role of the Transient Atropisomerism and Chirality of Flurbiprofen Unveiled by Laser-Ablation Rotational Spectroscopy.

The combination of atropisomerism and chirality in flurbiprofen is shown to be relevant concerning its pharmacological activity. The two most stable conformers of a total of eight theoretically predicted for each R- or S- flurbiprofen enantiomers have been isolated in the cooling conditions of a supersonic jet and structurally characterized by laser ablation Fourier transform microwave spectroscopy. The detected conformers, whose structure is mainly defined by three dihedral angles, only differ in the sign of the phenyl torsion angle giving rise to Sa and Ra atropisomers. A comparison with the structures available for the R- and S- enantiomers complexed to COX isoforms reveals that the enzymes select only the Sa atropisomers, resulting in a diastereoisomer-specific recognition. The most stable gas phase conformer is exclusively selected when using the S- enantiomer while the second is recognized only for the R- enantiomer. These experimental results highlight the importance of atropisomerism in drug design.

Elastic Electron scattering by thermal mixture of glycine conformers in gas phase

A theoretical study of electron scattering by a thermal mixture of glycine molecules in the energy range from 1 to 10 eV is performed using the UK-RMol codes which are based on the R-Matrix method. The six lowest relative Gibbs free energies glycine conformers considered, Ip, IIp, IIn, IIIp, IIIn and IVn, are significantly populated in thermal mixtures. All these conformers present similar resonance structures in the eigenphase sums: a lower-energy resonance state near 1.8 eV and another at higher-energy above 7 eV. For the six conformers the lowest resonance lies between 1.75 eV and 2.21 eV. The very large dipole moments of 6.32 D and 5.67 D for IIp and IIn, respectively, makes the magnitude of their cross sections significantly larger than other conformers, which increases the average cross sections in thermal mixtures compared with the cross sections of the lowest energy Ip conformer. Three conformer population sets are used to calculate the averaged differential and integral cross sections: two theoretical sets based on the relative Gibbs free energies and another set that aims to mimic experiment based on the observed populations. The averaged cross sections are similar for all population sets, but differ from the Ip conformer cross section. This suggests that, for large and flexible molecules, the computed average cross sections should be used when comparing with experimental data.

Infrared action spectroscopy of the deprotonated formic acid trimer, trapped in helium nanodroplets.

Hydrogen bonding interactions are essential in the structural stabilization and physicochemical properties of complex molecular systems, and carboxylic acid functional groups are common participants in these motifs. Consequently, the neutral formic acid (FA) dimer has been extensively investigated in the past, as it represents a useful model system to investigate proton donor-acceptor interactions. The analogous deprotonated dimers, in which two carboxylate groups are bound by a single proton, have also served as informative model systems. In these complexes, the position of the shared proton is mainly determined by the proton affinity of the carboxylate units. However, very little is known about the nature of the hydrogen bonding interactions in systems containing more than two carboxylate units. Here we report a study on the deprotonated (anionic) FA trimer. IR spectra are recorded in the 400-2000 cm-1 spectral range by means of vibrational action spectroscopy of FA trimer ions embedded in helium nanodroplets. Characterization of the gas-phase conformer and assignment of the vibrational features is achieved by comparing the experimental results with electronic structure calculations. To assist in the assignments, the 2H and 18O FA trimer anion isotopologues are also measured under the same experimental conditions. Comparison between the experimental and computed spectra, especially the observed shifts in spectral line positions upon isotopic substitution of the exchangeable protons, suggests that the prevalent conformer, under the experimental conditions, exhibits a planar structure that resembles the crystalline structure of formic acid.

Native Electrospray Ionization of Multi-Domain Proteins via a Bead Ejection Mechanism.

Native ion mobility mass spectrometry is potentially useful for the biophysical characterization of proteins, as the electrospray charge state distribution and the collision cross section distribution depend on their solution conformation. We examine here the charging and gas-phase conformation of multi-domain therapeutic proteins comprising globular domains tethered by disordered linkers. The charge and collision cross section distributions are multimodal, suggesting several conformations in solution, as confirmed by solution hydrogen/deuterium exchange. The most intriguing question is the ionization mechanism of these structures: a fraction of the population does not follow the charged residue mechanism but cannot ionize by pure chain ejection either. We deduce that a hybrid mechanism is possible, wherein globular domains are ejected one at a time from a parent droplet. The charge vs solvent accessible surface area correlations of denatured and intrinsically disordered proteins are also compatible with this "bead ejection mechanism", which we propose as a general tenet of biomolecule electrospray.

Structure and Conformation of Novel BODIPY Ugi Adducts.

Two novel BODIPY‐Ugi (boron dipyrromethene) adducts exhibit peculiar room temperature (T=20 °C) H‐1 NMR spectra in that several protons located at the aromatic aniline‐type ring are lost in the baseline. This observation revealed the existence of a dynamic conformational process where rotation around the C−N bond is hindered. Variable‐temperature H‐1 and C‐13 NMR spectroscopic analysis confirmed this conclusion; that is, low‐temperature spectra show distinct signals for all four aromatic protons below coalescence, whereas average signals are recorded above coalescence (T=+120 °C). Particularly interesting was the rather large difference in chemical shifts for the ortho protons below coalescence, Δδ=1.45 ppm, which was explained based on DFT computational analysis. Indeed, the calculated lowest‐energy gas‐phase conformation of the BODIPY Ugi adducts locates one half of the aniline‐type ring in the shielding anisotropic cone of the bridge phenyl ring in the BODIPY segment. This is in contrast to the solid‐state conformation established by X‐ray diffraction analysis that shows a nearly parallel arrangement of the aromatic rings, probably induced by crystal packing forces.

Ion Mobility-Mass Spectrometry Coupled to Droplet Microfluidics for Rapid Protein Structure Analysis and Drug Discovery.

Native mass spectrometry coupled to ion mobility (IM-MS) has become an important tool for the investigation of protein structure and dynamics upon ligand binding. Additionally, collisional activation or collision induced unfolding (CIU) can further probe conformational changes induced by ligand binding; however, larger scale screens have not been implemented due to limitations associated with throughput and sample introduction. In this work we explore the high-throughput capabilities of CIU fingerprinting. Fingerprint collection times were reduced 10-fold over traditional data collections through the use of improved smoothing and interpolation algorithms. Fast-CIU was then coupled to a droplet sample introduction approach using 40 nL droplet sample volumes and 2 s dwell times at each collision voltage. This workflow, which increased throughput by ∼16-fold over conventional nanospray CIU methods, was applied to a 96-compound screen against Sirtuin-5, a protein target of clinical interest. Over 20 novel Sirtuin-5 binders were identified, and it was found that Sirtuin-5 inhibitors will stabilize specific Sirtuin-5 gas-phase conformations. This work demonstrates that droplet-CIU can be implemented as a high-throughput biophysical characterization approach. Future work will focus on improving the throughput of this workflow and on automating data acquisition and analysis.

Broad chemical transferability in structure-based coarse-graining.

Compared to top-down coarse-grained (CG) models, bottom-up approaches are capable of offering higher structural fidelity. This fidelity results from the tight link to a higher resolution reference, making the CG model chemically specific. Unfortunately, chemical specificity can be at odds with compound-screening strategies, which call for transferable parameterizations. Here, we present an approach to reconcile bottom-up, structure-preserving CG models with chemical transferability. We consider the bottom-up CG parameterization of 3441 C7O2 small-molecule isomers. Our approach combines atomic representations, unsupervised learning, and a large-scale extended-ensemble force-matching parameterization. We first identify a subset of 19 representative molecules, which maximally encode the local environment of all gas-phase conformers. Reference interactions between the 19 representative molecules were obtained from both homogeneous bulk liquids and various binary mixtures. An extended-ensemble parameterization over all 703 state points leads to a CG model that is both structure-based and chemically transferable. Remarkably, the resulting force field is on average more structurally accurate than single-state-point equivalents. Averaging over the extended ensemble acts as a mean-force regularizer, smoothing out both force and structural correlations that are overly specific to a single-state point. Our approach aims at transferability through a set of CG bead types that can be used to easily construct new molecules while retaining the benefits of a structure-based parameterization.

A round dance of acetaldehyde molecular ensembles on Rh(111) surface; formation and decomposition of various paraldehyde conformers

The adsorption, desorption and polymerization of acetaldehyde (AA) on Rh(111) surface were studied by TDS, HREELS, work function measurements and theoretical calculations. Polymerization, as a new feature after adsorption of AA on Rh(111), makes the interpretation of the surface processes more exciting. The single molecule itself has at least two possibilities for adsorption on Rh(111) surface, namely η1−(O)−CH3CHOa and η2−(O,C)−CH3CHOa adsorption forms are considered. Moreover, η1−(O)−CH3CHOa aldehyde molecules can react with each other forming oligomers or even longer polymers, as well. Paraldehyde, the triple molecular mass species of AA, shows some surprising kinetic and thermodynamic features on clean Rh(111) surface. However, the observed surface formation of the oligomer and its thermodynamically most stable gas phase conformer seems to be in contradiction. The interpretation of this phenomenon is the goal of the present study.

De Novo Crystal Structure Determination from Machine Learned Chemical Shifts.

Determination of the three-dimensional atomic-level structure of powdered solids is one of the key goals in current chemistry. Solid-state NMR chemical shifts can be used to solve this problem, but they are limited by the high computational cost associated with crystal structure prediction methods and density functional theory chemical shift calculations. Here, we successfully determine the crystal structures of ampicillin, piroxicam, cocaine, and two polymorphs of the drug molecule AZD8329 using on-the-fly generated machine-learned isotropic chemical shifts to directly guide a Monte Carlo-based structure determination process starting from a random gas-phase conformation.

On the Conformation of Anionic Peptoids in the Gas Phase.

Although N-(S)-phenylethyl peptoids are known to adopt helical structures in solutions, the corresponding positively charged ions lose their helical structure during the transfer from the solution to the gas phase due to the so-called charge solvation effect. We, here, considered negatively charged peptoids to investigate by ion mobility spectrometry-mass spectrometry whether the structural changes described in the positive ionization mode can be circumvented in the negative mode by a fine-tuning of the peptoid sequence, that is, by positioning the negative charge at the positive side of the helical peptoid macrodipole. N-(S)-(1-carboxy-2-phenylethyl) (Nscp) and N-(S)-phenylethyl (Nspe) were selected as the negative charge carrier and as the helix inductor, respectively. We, here, report the results of a joint theoretical and experimental study demonstrating that the structures adopted by the NspenNscp anions remain compactly folded in the gas phase for chains containing up to 10 residues, whereas no evidence of the presence of a helical structure was obtained, even if, for selected sequences and lengths, different gas phase conformations are detected.

Infrared multiple-photon dissociation spectroscopy of cationized glycine: effects of alkali metal cation size on gas-phase conformation.

The gas-phase structures of cationized glycine (Gly), including complexes with Li+, Na+, K+, Rb+, and Cs+, are examined using infrared multiple-photon dissociation (IRMPD) spectroscopy utilizing light generated by a free electron laser, in conjunction with ab initio calculations. To identify the structures present in the experimental studies, measured IRMPD spectra are compared to spectra calculated at B3LYP/6-311+G(d,p) for the Li+, Na+, and K+ complexes and at B3LYP/def2TZVP for the Rb+ and Cs+ complexes. Single-point energy calculations were carried out at the B3LYP, B3P86, and MP2(full) levels using the 6-311+G(2d,2p) basis set for Li+, Na+, K+ and the def2TZVPP basis set for Rb+ and Cs+. The Li+ and Na+ complexes are identified as metal cation coordination to the amino nitrogen and carbonyl oxygen, [N,CO]-tt, although Na+(Gly) may have contributions from additional structures. The heavier metal cations coordinate to either the carbonyl oxygen, [CO]-cc, or the carbonyl oxygen and hydroxy oxygen, [CO,OH]-cc, with the former apparently preferred for Rb+ and Cs+ and the latter for K+. These two structures reside in a double-well potential and different levels of theory predict very different relative stabilities. Some experimental evidence is provided that MP2(full) theory provides the most accurate relative energies.

Charge Engineering Reveals the Roles of Ionizable Side Chains in Electrospray Ionization Mass Spectrometry.

In solution, the charge of a protein is intricately linked to its stability, but electrospray ionization distorts this connection, potentially limiting the ability of native mass spectrometry to inform about protein structure and dynamics. How the behavior of intact proteins in the gas phase depends on the presence and distribution of ionizable surface residues has been difficult to answer because multiple chargeable sites are present in virtually all proteins. Turning to protein engineering, we show that ionizable side chains are completely dispensable for charging under native conditions, but if present, they are preferential protonation sites. The absence of ionizable side chains results in identical charge state distributions under native-like and denaturing conditions, while coexisting conformers can be distinguished using ion mobility separation. An excess of ionizable side chains, on the other hand, effectively modulates protein ion stability. In fact, moving a single ionizable group can dramatically alter the gas-phase conformation of a protein ion. We conclude that although the sum of the charges is governed solely by Coulombic terms, their locations affect the stability of the protein in the gas phase.