Structural stability of electrosprayed proteins: temperature and hydration effects

Collision induced unfolding and molecular dynamics simulations of norovirus capsid dimers reveal strain-specific stability profiles

Collision induced unfolding (CIU) is a method used with ion mobility mass spectrometry to examine protein structures and their stability. Such experiments yield information about higher order protein structures, yet are unable to provide details about the underlying processes. That information can however be provided using molecular dynamics simulations. Here, we investigate the gas-phase unfolding of norovirus capsid dimers from the Norwalk and Kawasaki strains by employing molecular dynamics simulations over a range of temperatures, representing different levels of activation, together with CIU experiments. The dimers have highly similar structures, but their CIU reveals different stability that can be explained by the different dynamics that arises in response to the activation seen in the simulations, including a part of the sequence with previously observed strain-specific dynamics in solution. Our findings show how similar protein variants can be examined using mass spectrometric techniques in conjunction with atomistic molecular dynamics simulations to reveal differences in stability as well as differences in how and where unfolding takes place upon activation.

Enhanced EMC-Advantages of partially known orientations in x-ray single particle imaging

Single particle imaging of proteins in the gas phase with x-ray free-electron lasers holds great potential to study fast protein dynamics, but is currently limited by weak and noisy data. A further challenge is to discover the proteins' orientation as each protein is randomly oriented when exposed to x-rays. Algorithms such as the expand, maximize, and compress (EMC) exist that can solve the orientation problem and reconstruct the three-dimensional diffraction intensity space, given sufficient measurements. If information about orientation were known, for example, by using an electric field to orient the particles, the reconstruction would benefit and potentially reach better results. We used simulated diffraction experiments to test how the reconstructions from EMC improve with particles' orientation to a preferred axis. Our reconstructions converged to correct maps of the three-dimensional diffraction space with fewer measurements if biased orientation information was considered. Even for a moderate bias, there was still significant improvement. Biased orientations also substantially improved the results in the case of missing central information, in particular in the case of small datasets. The effects were even more significant when adding a background with 50% the strength of the averaged diffraction signal photons to the diffraction patterns, sometimes reducing the data requirement for convergence by a factor of 10. This demonstrates the usefulness of having biased orientation information in single particle imaging experiments, even for a weaker bias than what was previously known. This could be a key component in overcoming the problems with background noise that currently plague these experiments.

Cryo-EM of soft-landed β-galactosidase: Gas-phase and native structures are remarkably similar

Native mass spectrometry (MS) has become widely accepted in structural biology, providing information on stoichiometry, interactions, homogeneity, and shape of protein complexes. Yet, the fundamental assumption that proteins inside the mass spectrometer retain a structure faithful to native proteins in solution remains a matter of intense debate. Here, we reveal the gas-phase structure of β-galactosidase using single-particle cryo-electron microscopy (cryo-EM) down to 2.6-Å resolution, enabled by soft landing of mass-selected protein complexes onto cold transmission electron microscopy (TEM) grids followed by in situ ice coating. We find that large parts of the secondary and tertiary structure are retained from the solution. Dehydration-driven subunit reorientation leads to consistent compaction in the gas phase. By providing a direct link between high-resolution imaging and the capability to handle and select protein complexes that behave problematically in conventional sample preparation, the approach has the potential to expand the scope of both native mass spectrometry and cryo-EM.

Observation of a single protein by ultrafast X-ray diffraction

The idea of using ultrashort X-ray pulses to obtain images of single proteins frozen in time has fascinated and inspired many. It was one of the arguments for building X-ray free-electron lasers. According to theory, the extremely intense pulses provide sufficient signal to dispense with using crystals as an amplifier, and the ultrashort pulse duration permits capturing the diffraction data before the sample inevitably explodes. This was first demonstrated on biological samples a decade ago on the giant mimivirus. Since then, a large collaboration has been pushing the limit of the smallest sample that can be imaged. The ability to capture snapshots on the timescale of atomic vibrations, while keeping the sample at room temperature, may allow probing the entire conformational phase space of macromolecules. Here we show the first observation of an X-ray diffraction pattern from a single protein, that of Escherichia coli GroEL which at 14 nm in diameter is the smallest biological sample ever imaged by X-rays, and demonstrate that the concept of diffraction before destruction extends to single proteins. From the pattern, it is possible to determine the approximate orientation of the protein. Our experiment demonstrates the feasibility of ultrafast imaging of single proteins, opening the way to single-molecule time-resolved studies on the femtosecond timescale.

Solidification and oral delivery of biologics to the colon- A review

The oral delivery of biologics such as therapeutic proteins, peptides and oligonucleotides for the treatment of colon related diseases has been the focus of increasing attention over the last years. However, the major disadvantage of these macromolecules is their degradation propensity in liquid state which can lead to the undesirable and complete loss of function. Therefore, to increase the stability of the biologic and reduce their degradation propensity, formulation techniques such as solidification can be performed to obtain a stable solid dosage form for oral administration. Due to their fragility, stress exerted on the biologic during solidification has to be reduced with the incorporation of stabilizing excipients into the formulation. This review focuses on the state-of-the-art solidification techniques required to obtain a solid dosage form for the oral delivery of biologics to the colon and the use of suitable excipients for adequate stabilization upon solidification. The solidifying processes discussed within this review are spray drying, freeze drying, bead coating and also other techniques such as spray freeze drying, electro spraying, vacuum- and supercritical fluid drying. Further, the colon as site of absorption in both healthy and diseased state is critically reviewed and possible oral delivery systems for biologics are discussed.

Evaluation for Ion Heating of H2A-H2B Dimer in Ion Mobility Spectrometry-Mass Spectrometry

Ion mobility spectrometry-mass spectrometry (IMS-MS) provides m/z values and collision cross sections (CCSs) of gas-phase ions. In our previous study, an intrinsically disordered protein, the H2A-H2B dimer, was analyzed using IMS-MS, resulting in two conformational populations of CCS. Based on experimental and theoretical approaches, this resulted from a structural diversity of intrinsically disordered regions. We predicted that this phenomenon is related to ion heating in the IMS-MS instrument. In this study, to reveal the effect of ion heating from parameters in the IMS-MS instrument on the conformational population of the H2A-H2B dimer, we investigated the arrival time distributions of the H2A-H2B dimer by changing values of three instrumental parameters, namely, cone voltage located in the first vacuum chamber, trap collision energy (trap CE) for tandem mass spectrometry, and trap bias voltage for the entrance of IMS. These results revealed that the two populations observed for the H2A-H2B dimer were due to the trap bias voltage. Furthermore, to evaluate the internal energies of the analyte ions with respect to each parameter, benzylpyridinium derivatives were used as temperature-sensitive probes. The results showed that the trap CE voltage imparts greater internal energy to the ions than the trap bias voltage. In addition, this slight change in the internal energy caused by the trap bias voltage resulted in the structural diversity of the H2A-H2B dimer. Therefore, the trap bias voltage should be set with attention to the properties of the analytes, even if the effect of the trap bias voltage on the internal energy is negligible.

Water layer and radiation damage effects on the orientation recovery of proteins in single-particle imaging at an X-ray free-electron laser

The noise caused by sample heterogeneity (including sample solvent) has been identified as one of the determinant factors for a successful X-ray single-particle imaging experiment. It influences both the radiation damage process that occurs during illumination as well as the scattering patterns captured by the detector. Here, we investigate the impact of water layer thickness and radiation damage on orientation recovery from diffraction patterns of the nitrogenase iron protein. Orientation recovery is a critical step for single-particle imaging. It enables to sort a set of diffraction patterns scattered by identical particles placed at unknown orientations and assemble them into a 3D reciprocal space volume. The recovery quality is characterized by a "disconcurrence" metric. Our results show that while a water layer mitigates protein damage, the noise generated by the scattering from it can introduce challenges for orientation recovery and is anticipated to cause problems in the phase retrieval process to extract the desired protein structure. Compared to these disadvantageous effects due to the thick water layer, the effects of radiation damage on the orientation recovery are relatively small. Therefore, minimizing the amount of residual sample solvent should be considered a crucial step in improving the fidelity and resolution of X-ray single-particle imaging experiments.

Rehydration Post-orientation: Investigating Field-Induced Structural Changes via Computational Rehydration

Proteins can be oriented in the gas phase using strong electric fields, which brings advantages for structure determination using X-ray free electron lasers. Both the vacuum conditions and the electric-field exposure risk damaging the protein structures. Here, we employ molecular dynamics simulations to rehydrate and relax vacuum and electric-field exposed proteins in aqueous solution, which simulates a refinement of structure models derived from oriented gas-phase proteins. We find that the impact of the strong electric fields on the protein structures is of minor importance after rehydration, compared to that of vacuum exposure and ionization in electrospraying. The structures did not fully relax back to their native structure in solution on the simulated timescales of 200 ns, but they recover several features, including native-like intra-protein contacts, which suggests that the structures remain in a state from which the fully native structure is accessible. Our findings imply that the electric fields used in native mass spectrometry are well below a destructive level, and suggest that structures inferred from X-ray diffraction from gas-phase proteins are relevant for solution and in vivo conditions, at least after in silico rehydration.

GPU-Enhanced DFTB Metadynamics for Efficiently Predicting Free Energies of Biochemical Systems

Metadynamics calculations of large chemical systems with ab initio methods are computationally prohibitive due to the extensive sampling required to simulate the large degrees of freedom in these systems. To address this computational bottleneck, we utilized a GPU-enhanced density functional tight binding (DFTB) approach on a massively parallelized cloud computing platform to efficiently calculate the thermodynamics and metadynamics of biochemical systems. To first validate our approach, we calculated the free-energy surfaces of alanine dipeptide and showed that our GPU-enhanced DFTB calculations qualitatively agree with computationally-intensive hybrid DFT benchmarks, whereas classical force fields give significant errors. Most importantly, we show that our GPU-accelerated DFTB calculations are significantly faster than previous approaches by up to two orders of magnitude. To further extend our GPU-enhanced DFTB approach, we also carried out a 10 ns metadynamics simulation of remdesivir, which is prohibitively out of reach for routine DFT-based metadynamics calculations. We find that the free-energy surfaces of remdesivir obtained from DFTB and classical force fields differ significantly, where the latter overestimates the internal energy contribution of high free-energy states. Taken together, our benchmark tests, analyses, and extensions to large biochemical systems highlight the use of GPU-enhanced DFTB simulations for efficiently predicting the free-energy surfaces/thermodynamics of large biochemical systems.

Cryo-EM samples of gas-phase purified protein assemblies using native electrospray ion-beam deposition

An increasing number of studies on biomolecular function indirectly combine mass spectrometry (MS) with imaging techniques such as cryo electron microscopy (cryo-EM). This approach allows information on the homogeneity, stoichiometry, shape, and interactions of native protein complexes to be obtained, complementary to high-resolution protein structures. We have recently demonstrated TEM sample preparation via native electrospray ion-beam deposition (ES-IBD) as a direct link between native MS and cryo-EM. This workflow forms a potential new route to the reliable preparation of homogeneous cryo-EM samples and a better understanding of the relation between native solution-phase and native-like gas-phase structures. However, many aspects of the workflow need to be understood and optimized to obtain performance comparable to that of state-of-the-art cryo-EM. Here, we expand on the previous discussion of key factors by probing the effects of substrate type and deposition energy. We present and discuss micrographs from native ES-IBD samples with amorphous carbon, graphene, and graphene oxide, as well as landing energies in the range between 2 and 150 eV per charge.

Stability and conformational memory of electrosprayed and rehydrated bacteriophage MS2 virus coat proteins

Proteins are innately dynamic, which is important for their functions, but which also poses significant challenges when studying their structures. Gas-phase techniques can utilise separation and a range of sample manipulations to transcend some of the limitations of conventional techniques for structural biology in crystalline or solution phase, and isolate different states for separate interrogation. However, the transfer from solution to the gas phase risks affecting the structures, and it is unclear to what extent different conformations remain distinct in the gas phase, and if resolution in silico can recover the native conformations and their differences. Here, we use extensive molecular dynamics simulations to study the two distinct conformations of dimeric capsid protein of the MS2 bacteriophage. The protein undergoes notable restructuring of its peripheral parts in the gas phase, but subsequent simulation in solvent largely recovers the native structure. Our results suggest that despite some structural loss due to the experimental conditions, gas-phase structural biology techniques provide meaningful data that inform not only about the structures but also conformational dynamics of proteins.

Mass-selective and ice-free electron cryomicroscopy protein sample preparation via native electrospray ion-beam deposition

Despite tremendous advances in sample preparation and classification algorithms for electron cryomicroscopy (cryo-EM) and single-particle analysis (SPA), sample heterogeneity remains a major challenge and can prevent access to high-resolution structures. In addition, optimization of preparation conditions for a given sample can be time-consuming. In the current work, it is demonstrated that native electrospray ion-beam deposition (native ES-IBD) is an alternative, reliable approach for the preparation of extremely high-purity samples, based on mass selection in vacuum. Folded protein ions are generated by native electrospray ionization, separated from other proteins, contaminants, aggregates, and fragments, gently deposited on cryo-EM grids, frozen in liquid nitrogen, and subsequently imaged by cryo-EM. We demonstrate homogeneous coverage of ice-free cryo-EM grids with mass-selected protein complexes. SPA reveals that the complexes remain folded and assembled, but variations in secondary and tertiary structures are currently limiting information in 2D classes and 3D EM density maps. We identify and discuss challenges that need to be addressed to obtain a resolution comparable to that of the established cryo-EM workflow. Our results show the potential of native ES-IBD to increase the scope and throughput of cryo-EM for protein structure determination and provide an essential link between gas-phase and solution-phase protein structures.

Protonation of Serine in Gas and Condensed and Microsolvated States in Aqueous Solution

Identification of molecules and elucidation of their chemical structure are ubiquitous problems in chemistry. Mass spectrometry (MS) can be used due to its sensitivity and versatility. For detection to occur, analytes must be ionized and transferred to the gas phase. Soft ionization processes such as electrospray ionization are popular; however, resulting microsolvated phases can alter the chemistry of analytes and therefore detection and identification. To understand these processes, we use computational methods to probe the ionization propensity of serine in the gas phase, aqueous microsolvated clusters, and aqueous solution. We show that the tautomeric form of serine is altered by the presence of water, as five water molecules can stabilize the zwitterionic tautomer. Inclusion of cosolutes such as ions can stabilize the zwitterion with as few as one or two water molecules present. We demonstrate that ionization propensity, as measured by gas phase bacisity, can increase by over 100 kJ/mol when placed in a small water-serine cluster, showing the sensitivity of the chemistry of microsolvated analytes. Finally, detailed analysis reveals that small droplets (less than seven water molecules) are extremely sensitive to addition of further water molecules. Beyond this limit, structural and electronic properties change little with droplet size.

In a flash of light: X-ray free electron lasers meet native mass spectrometry

During the last years, X-ray free electron lasers (XFELs) have emerged as X-ray sources of unparalleled brightness, delivering extreme amounts of photons in femtosecond pulses. As such, they have opened up completely new possibilities in drug discovery and structural biology, including studying high resolution biomolecular structures and their functioning in a time resolved manner, and diffractive imaging of single particles without the need for their crystallization. In this perspective, we briefly review the operation of XFELs, their immediate uses for drug discovery and focus on the potentially revolutionary single particle diffractive imaging technique and the challenges which remain to be overcome to fully realize its potential to provide high resolution structures without the need for crystallization, freezing or the need to keep proteins stable at extreme concentrations for long periods of time. As the issues have been to a large extent sample delivery related, we outline a way for native mass spectrometry to overcome these and enable so far impossible research with a potentially huge impact on structural biology and drug discovery, such as studying structures of transient intermediate species in viral life cycles or during functioning of molecular machines.

Protein orientation in time-dependent electric fields: orientation before destruction

Proteins often have nonzero electric dipole moments, making them interact with external electric fields and offering a means for controlling their orientation. One application that is known to benefit from orientation control is single-particle imaging with x-ray free-electron lasers, in which diffraction is recorded from proteins in the gas phase to determine their structures. To this point, theoretical investigations into this phenomenon have assumed that the field experienced by the proteins is constant or a perfect step function, whereas any real-world pulse will be smooth. Here, we explore the possibility of orienting gas-phase proteins using time-dependent electric fields. We performed ab initio simulations to estimate the field strength required to break protein bonds, with 45 V/nm as a breaking point value. We then simulated ubiquitin in time-dependent electric fields using classical molecular dynamics. The minimal field strength required for orientation within 10 ns was on the order of 0.5 V/nm. Although high fields can be destructive for the structure, the structures in our simulations were preserved until orientation was achieved regardless of field strength, a principle we denote “orientation before destruction.”

Reproducibility in the unfolding process of protein induced by an external electric field

The dynamics of proteins are crucial for their function. However, commonly used techniques for studying protein structures are limited in monitoring time-resolved dynamics at high resolution. Combining electric fields with existing techniques to study gas-phase proteins, such as single particle imaging using free-electron lasers and gas-phase small angle X-ray scattering, has the potential to open up a new era in time-resolved studies of gas-phase protein dynamics. Using molecular dynamics simulations, we identify well-defined unfolding pathways of a protein, induced by experimentally achievable external electric fields. Our simulations show that strong electric fields in conjunction with short-pulsed X-ray sources such as free-electron lasers can be a new path for imaging dynamics of gas-phase proteins at high spatial and temporal resolution.

Structural Heterogeneity in Single Particle Imaging Using X-ray Lasers

One of the challenges facing single particle imaging with ultrafast X-ray pulses is the structural heterogeneity of the sample to be imaged. For the method to succeed with weakly scattering samples, the diffracted images from a large number of individual proteins need to be averaged. The more the individual proteins differ in structure, the lower the achievable resolution in the final reconstructed image. We use molecular dynamics to simulate two globular proteins in vacuum, fully desolvated as well as with two different solvation layers, at various temperatures. We calculate the diffraction patterns based on the simulations and evaluate the noise in the averaged patterns arising from the structural differences and the surrounding water. Our simulations show that the presence of a minimal water coverage with an average 3 Å thickness will stabilize the protein, reducing the noise associated with structural heterogeneity, whereas additional water will generate more background noise.

Perspectives on single particle imaging with x rays at the advent of high repetition rate x-ray free electron laser sources

X-ray free electron lasers (XFELs) now routinely produce millijoule level pulses of x-ray photons with tens of femtoseconds duration. Such x-ray intensities gave rise to the idea that weakly scattering particles-perhaps single biomolecules or viruses-could be investigated free of radiation damage. Here, we examine elements from the past decade of so-called single particle imaging with hard XFELs. We look at the progress made to date and identify some future possible directions for the field. In particular, we summarize the presently achieved resolutions as well as identifying the bottlenecks and enabling technologies to future resolution improvement, which in turn enables application to samples of scientific interest.

Is radiation damage the limiting factor in high-resolution single particle imaging with X-ray free-electron lasers?

The prospect of single particle imaging with atomic resolution is one of the scientific drivers for the development of X-ray free-electron lasers. The assumption since the beginning has been that damage to the sample caused by intense X-ray pulses is one of the limiting factors for achieving subnanometer X-ray imaging of single particles and that X-ray pulses need to be as short as possible. Based on the molecular dynamics simulations of proteins in X-ray fields of various durations (5 fs, 25 fs, and 50 fs), we show that the noise in the diffracted signal caused by radiation damage is less than what can be expected from other sources, such as sample inhomogeneity and X-ray shot-to-shot variations. These findings show a different aspect of the feasibility of high-resolution single particle imaging using free-electron lasers, where employing X-ray pulses of longer durations could still provide a useful diffraction signal above the noise due to the Coulomb explosion.

Reproducibility of single protein explosions induced by X-ray lasers

Single particle imaging (SPI) using X-ray pulses has become increasingly attainable with the advent of high-intensity free electron lasers. Eliminating the need for crystallized samples enables structural studies of molecules previously inaccessible by conventional crystallography. While this emerging technique already demonstrates substantial promise, some obstacles need to be overcome before SPI can reach its full potential. One such problem is determining the spatial orientation of the sample at the time of X-ray interaction. Existing solutions rely on diffraction data and are computationally demanding and sensitive to noise. In this in silico study, we explore the possibility of aiding these methods by mapping the ion distribution as the sample undergoes a Coulomb explosion following the intense ionization. By detecting the ions ejected from the fragmented sample, the orientation of the original sample should be possible to determine. Knowledge of the orientation has been shown earlier to be of substantial advantage in the reconstruction of the original structure. 150 explosions of each of twelve separate systems - four polypeptides with different amounts of surface bound water - were simulated with molecular dynamics (MD) and the average angular distribution of carbon and sulfur ions was investigated independently. The results show that the explosion maps are reproducible in both cases, supporting the idea that orientation information is preserved. Additional water seems to restrict the carbon ion trajectories further through a shielding mechanism, making the maps more distinct. For sulfurs, water has no significant impact on the trajectories, likely due to their higher mass and greater ionization cross section, indicating that they could be of particular interest. Based on these findings, we conclude that explosion data can aid spatial orientation in SPI experiments and could substantially improve the capabilities of the novel technique.

Slow Dynamics of Tryptophan-Water Networks in Proteins

Water has a profound effect on the dynamics of biomolecules and governs many biological processes, leading to the concept that function is slaved to solvent dynamics within and surrounding the biomolecule. Protein conformational changes on μs-ms time scales are frequently associated with protein function, but little is known about the behavior of protein-bound water on these time scales. Here we have used NMR relaxation dispersion measurements to probe the tryptophan indoles in the enzyme dihydrofolate reductase (DHFR). We find that during structural changes on the μs-ms time scale, large chemical shift changes are often observed for the NH proton on the indole ring, while relatively smaller chemical shift changes are observed for the ring nitrogen atom. Comparison with experimental chemical shifts and density functional theory-based chemical shift predictions show that during the structural change the tryptophan indole NHs remain bound to water, but the geometry of the protein-bound water networks changes. These results establish that relaxation dispersion measurements can indirectly probe water dynamics and indicate that water can influence, or be influenced by, protein conformational changes on the μs-ms time scale. Our data show that structurally conserved bound water molecules can play a critical role in transmitting information between functionally important regions of the protein and provide evidence that internal protein motions can be coupled through the mediation of hydrogen-bonded water bound in the protein structure.

Controlling Protein Orientation in Vacuum Using Electric Fields

Single-particle imaging using X-ray free-electron lasers is an emerging technique that could provide high-resolution structures of macromolecules in the gas phase. One of the largest difficulties in realizing this goal is the unknown orientation of the individual sample molecules at the time of exposure. Preorientation of the molecules has been identified as a possible solution to this problem. Using molecular dynamics simulations, we identify a range of electric field strengths where proteins become oriented without losing their structure. For a number of experimentally relevant cases we show that structure determination is possible only when orientation information is included in the orientation-recovery process. We conclude that nondestructive field orientation of intact proteins is feasible and that it enables a range of new structural investigations with single-particle imaging.

Terahertz Frequency Spectroscopy to Determine Cold Shock Protein Stability upon Solvation and Evaporation - A Molecular Dynamics Study

Infrared (IR) and Terahertz (THz) spectroscopy simulations were carried out using CHARMM35b2 to determine protein stability. The stabilities of three bacterial cold shock proteins (Csps) originating from mesophiles, thermophiles and hyper- thermophiles respectively were investigated in this study. The three different Csps were investigated by Normal-Mode analysis and Molecular Dynamics simulation of THz spectra using the Hessian matrix for solvated systems, interpreted in the harmonic approximation at optimum near-melting temperatures of each homologue, by incorporating differences in the hydrous and anhydrous states of the Csps. The results show slight variations in the large scale protein motion. However, the IR spectra of Csps observed at the low frequency saddle surface region, clearly distinguishes the thermophilic and mesophilic proteins based on their stability. Further studies on protein stability employing low-frequency collective modes have the potential to reveal functionally important conformational changes that are biologically significant.

Integrating mass spectrometry with MD simulations reveals the role of lipids in Na+/H+ antiporters

Na+/H+ antiporters are found in all kingdoms of life and exhibit catalysis rates that are among the fastest of all known secondary-active transporters. Here we combine ion mobility mass spectrometry and molecular dynamics simulations to study the conformational stability and lipid-binding properties of the Na+/H+ exchanger NapA from Thermus thermophilus and compare this to the prototypical antiporter NhaA from Escherichia coli and the human homologue NHA2. We find that NapA and NHA2, but not NhaA, form stable dimers and do not selectively retain membrane lipids. By comparing wild-type NapA with engineered variants, we show that the unfolding of the protein in the gas phase involves the disruption of inter-domain contacts. Lipids around the domain interface protect the native fold in the gas phase by mediating contacts between the mobile protein segments. We speculate that elevator-type antiporters such as NapA, and likely NHA2, use a subset of annular lipids as structural support to facilitate large-scale conformational changes within the membrane.

Water Determines the Structure and Dynamics of Proteins

Water is an essential participant in the stability, structure, dynamics, and function of proteins and other biomolecules. Thermodynamically, changes in the aqueous environment affect the stability of biomolecules. Structurally, water participates chemically in the catalytic function of proteins and nucleic acids and physically in the collapse of the protein chain during folding through hydrophobic collapse and mediates binding through the hydrogen bond in complex formation. Water is a partner that slaves the dynamics of proteins, and water interaction with proteins affect their dynamics. Here we provide a review of the experimental and computational advances over the past decade in understanding the role of water in the dynamics, structure, and function of proteins. We focus on the combination of X-ray and neutron crystallography, NMR, terahertz spectroscopy, mass spectroscopy, thermodynamics, and computer simulations to reveal how water assist proteins in their function. The recent advances in computer simulations and the enhanced sensitivity of experimental tools promise major advances in the understanding of protein dynamics, and water surely will be a protagonist.

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.

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

Protein electrospray ionization (ESI) mass spectrometry (MS)-based techniques are widely used to provide insight into structural proteomics under the assumption that non-covalent protein complexes being transferred into the gas phase preserve basically the same intermolecular interactions as in solution. Here we investigate the applicability of this assumption by extending our previous structural prediction protocol for single proteins in ESI-MS to protein complexes. We apply our protocol to the human insulin dimer (hIns2) as a test case. Our calculations reproduce the main charge and the collision cross section (CCS) measured in ESI-MS experiments. Molecular dynamics simulations for 0.075 ms show that the complex maximizes intermolecular non-bonded interactions relative to the structure in water, without affecting the cross section. The overall gas-phase structure of hIns2 does exhibit differences with the one in aqueous solution, not inferable from a comparison with calculated CCS. Hence, care should be exerted when interpreting ESI-MS proteomics data based solely on NMR and/or X-ray structural information.

Dynamics of lysozyme and its hydration water under an electric field

The effects of a static electric field on the dynamics of lysozyme and its hydration water are investigated by means of incoherent quasi-elastic neutron scattering (QENS). Measurements were performed on lysozyme samples, hydrated respectively with heavy water (D2O) to capture the protein dynamics and with light water (H2O), to probe the dynamics of the hydration shell, in the temperature range from 210 < T < 260 K. The hydration fraction in both cases was about ∼ 0.38 gram of water per gram of dry protein. The field strengths investigated were respectively 0 kV/mm and 2 kV/mm (~2 × 10(6) V/m) for the protein hydrated with D2O and 0 kV and 1 kV/mm for the H2O-hydrated counterpart. While the overall internal protons dynamics of the protein appears to be unaffected by the application of an electric field up to 2 kV/mm, likely due to the stronger intra-molecular interactions, there is also no appreciable quantitative enhancement of the diffusive dynamics of the hydration water, as would be anticipated based on our recent observations in water confined in silica pores under field values of 2.5 kV/mm. This may be due to the difference in surface interactions between water and the two adsorption hosts (silica and protein), or to the existence of a critical threshold field value Ec ~2-3 kV/mm for increased molecular diffusion, for which electrical breakdown is a limitation for our sample.

Dodecyl maltoside protects membrane proteins in vacuo

Molecular dynamics simulations have been used to characterize the effects of transfer from aqueous solution to a vacuum to inform our understanding of mass spectrometry of membrane-protein-detergent complexes. We compared two membrane protein architectures (an α-helical bundle versus a β-barrel) and two different detergent types (phosphocholines versus an alkyl sugar) with respect to protein stability and detergent packing. The β-barrel membrane protein remained stable as a protein-detergent complex in vacuum. Zwitterionic detergents formed conformationally destabilizing interactions with an α-helical membrane protein after detergent micelle inversion driven by dehydration in vacuum. In contrast, a nonionic alkyl sugar detergent resisted micelle inversion, maintaining the solution-phase conformation of the protein. This helps to explain the relative stability of membrane proteins in the presence of alkyl sugar detergents such as dodecyl maltoside.

Do charge state signatures guarantee protein conformations?

The extent to which proteins in the gas phase retain their condensed-phase structure is a hotly debated issue. Closely related to this is the degree to which the observed charge state reflects protein conformation. Evidence from electron capture dissociation, hydrogen/deuterium exchange, ion mobility, and molecular dynamics shows clearly that there is often a strong correlation between the degree of folding and charge state, with the most compact conformations observed for the lowest charge states. In this article, we address recent controversies surrounding the relationship between charge states and folding, focussing also on the manipulation of charge in solution and its effect on conformation. 'Supercharging' reagents that have been used to effect change in charge state can promote unfolding in the electrospray droplet. However for several protein complexes, supercharging does not appear to perturb the structure in that unfolding is not detected. Consequently, a higher charge state does not necessarily imply unfolding. Whilst the effect of charge manipulation on conformation remains controversial, there is strong evidence that a folded, compact state of a protein can survive in the gas phase, at least on a millisecond timescale. The exact nature of the side-chain packing and secondary structural elements in these compact states, however, remains elusive and prompts further research.

Role of cohesive energy in droplet fragmentation

Using molecular-dynamics simulation, we investigate the fragmentation behavior of droplets after collision with a wall. We demonstrate that the ratio of the impact to the cohesive energy E(coh) of the droplet is the key quantity characterizing the droplet fragmentation process. To show this both van der Waals-bonded Ar and N(2) droplets and polar H(2)O droplets are studied. If the impact energy per molecule E<(0.35-0.4)E(coh), the droplet is reflected without fragmenting. Beyond that impact energy fragmentation of the droplet abruptly starts. At E=E(coh), the fragmentation process already results in a fine dispersal of the droplet into daughter droplets; the maximum fragment contains only less than 4% of the initial droplet mass and around one-third of the droplet has been shattered into isolated molecules. The disintegration process continuously increases with collision energy. These findings are relevant for the process of droplet fragmentation as used in the method of impact desolvation of electrosprayed microdroplets mass spectrometry.

Molecular insight into conformational transition of amyloid β-peptide 42 inhibited by (-)-epigallocatechin-3-gallate probed by molecular simulations

Considerable experimental evidence indicates that (-)-epigallocatechin-3-gallate (EGCG) inhibits the fibrillogenesis of Aβ(42) and alleviates its associated cytotoxicity. However, the molecular mechanism of the inhibition effect of EGCG on the conformational transition of Aβ(42) remains unclear due to the limitations of current experimental techniques. In this work, molecular dynamics simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) analysis were coupled to better understand the issue. It was found that the direct interactions between EGCG and the peptide are the origin of its inhibition effects. Specifically, EGCG molecules expel water from the surface of the Aβ(42), cluster with each other, and interact directly with the peptide. The results of free energy decomposition calculated by MM-PBSA indicate that the nonpolar term contributes more than 71% to the binding free energy of the EGCG-Aβ(42) complex, while polar interactions (i.e., hydrogen bonding) play a minor role. It was identified that there are 12 important residues of Aβ(42) that strongly interact with EGCG (Phe4, Arg5, Phe19, Phe20, Glu22, Lys28, Gly29, Leu34-Gly37, and Ile41), while nonpolar interactions are mainly provided by the side chains of some hydrophobic residues (Phe, Met and Ile) and the main chains of some nonhydrophobic residues (Lys28 and Gly29). On the contrary, polar interactions are mainly formed by the main chain of Aβ(42), of which the main chains of Gly29 and Gly37 contribute greatly. The work has thus elucidated the molecular mechanism of the inhibition effect of EGCG on the conformational transition of Aβ(42), and the findings are considered critical for exploring more effective agents for the inhibition of Aβ(42) fibrillogenesis.

Coherent diffraction of a single virus particle: the impact of a water layer on the available orientational information

Coherent diffractive imaging using x-ray free-electron lasers (XFELs) may provide a unique opportunity for high-resolution structural analysis of single particles sprayed from an aqueous solution into the laser beam. As a result, diffraction images are measured from randomly oriented objects covered by a water layer. We analyze theoretically how the thickness of the covering water layer influences the structural and orientational information contained in the recorded diffraction images. This study has implications for planned experiments on single-particle imaging with XFELs.

Water evaporation and conformational changes from partially solvated ubiquitin

Using molecular dynamics simulation, we study the evaporation of water molecules off partially solvated ubiquitin. The evaporation and cooling rates are determined for a molecule at the initial temperature of 300 K. The cooling rate is found to be around 3 K/ns, and decreases with water temperature in the course of the evaporation. The conformation changes are monitored by studying a variety of intermediate partially solvated ubiquitin structures. We find that ubiquitin shrinks with decreasing hydration shell and exposes more of its hydrophilic surface area to the surrounding.