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

Dipole orientation of hydrated gas phase proteins

In the aerosolization of single proteins from solution, the proteins may be covered by a layer of water. This is relevant to consider in sample delivery for single particle imaging (SPI) with X-ray free-electron lasers. Previous studies suggest that the presence of a 3 Å water layer stabilizes the molecular structure and decreases structural heterogeneity which is important since it facilitates the structure determination in SPI. It has also been shown that SPI would benefit from the possibility of controlling the particle orientation in the interaction region. It has been proposed that such control would be possible by applying a DC electric field that interacts with the intrinsic dipole of the particle. This study investigates how SPI experiments, including dipole orientation, would be affected by the presence of a hydration layer covering the proteins. We investigated this by performing classical MD simulations of a globular protein in gas phase interacting with an external electric field. Two hydration levels were used: a fully desolvated molecule and one with a water layer corresponding to 3 Å covering the proteins surface. Our simulations show that a water layer enables the molecules to orient at lower field amplitudes, and on shorter time scales, as compared to the desolvated case. We also see a marginally larger stability of the molecular structure in the hydrated case at field strengths below 2 V nm-1. The presence of a water layer, in combination with an electric field, also tend to stabilize the dipole axis significantly within the molecular structure.

Laser-Induced Alignment of Nanoparticles and Macromolecules for Coherent-Diffractive-Imaging Applications

Laser-induced alignment of particles and molecules was long envisioned to support three-dimensional structure determination using "single-molecule diffraction" with X-ray free-electron lasers [PRL 92, 198102 (2004)]. However, the alignment of isolated macromolecules has not yet been demonstrated also because quantitative modeling is very expensive. We computationally demonstrated that the alignment of nanorods and proteins is possible with a standard laser technology. We performed a comprehensive analysis on the dependence of the degree of alignment on molecular properties and experimental details, e.g., particle temperature and laser-pulse energy. Considering the polarizability anisotropy of about 150,000 proteins, our analysis revealed that most of these proteins can be aligned using realistic experimental parameters.

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.

Electric Polarization-Dependent Absorption and Photocurrent Generation in Limnospira indica Immobilized on Boron-Doped Diamond

We present the change of light absorption of cyanobacteria in response to externally applied electrical polarization. Specifically, we studied the relation between electrical polarization and changes in light absorbance for a biophotoelectrode assembly comprising boron-doped diamond as semiconducting electrode and live Limnospira indicaPCC 8005 trichomes embedded in either polysaccharide (agar) or conductive conjugated polymer (PEDOT-PSS) matrices. Our study involves the monitoring of cyanobacterial absorbance and the measurement of photocurrents at varying wavelengths of illumination for switched electric fields, i.e., using the bioelectrode either as an anode or as cathode. We observed changes in the absorbance characteristics, indicating a direct causal relationship between electrical polarization and absorbing properties of L. indica. Our finding opens up a potential avenue for optimization of the performance of biophotovoltaic devices through controlled polarization. Furthermore, our results provide fundamental insights into the wavelength-dependent behavior of a bio photovoltaic system using live cyanobacteria.

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.

Full control of the orientation of non-symmetric molecules using weak and moderate electric fields

We investigate the full control over the orientation of a planar non-symmetric molecule by using moderate and weak electric fields. Quantum optimal control techniques allow us to orient any axis of 6-chloropyridazine-3-carbonitrile, which is taken as prototype example here, along the electric field direction. We perform a detailed analysis by exploring the impact on the molecular orientation of the time scale and strength of the control field. The underlying physical phenomena allowing for the control of the orientation are interpreted in terms of the frequencies contributing to the field-dressed dynamics and to the driving field by a spectral analysis.

Electrochemistry in sensing of molecular interactions of proteins and their behavior in an electric field

Electrochemical methods can be used not only for the sensitive analysis of proteins but also for deeper research into their structure, transport functions (transfer of electrons and protons), and sensing their interactions with soft and solid surfaces. Last but not least, electrochemical tools are useful for investigating the effect of an electric field on protein structure, the direct application of electrochemical methods for controlling protein function, or the micromanipulation of supramolecular protein structures. There are many experimental arrangements (modalities), from the classic configuration that works with an electrochemical cell to miniaturized electrochemical sensors and microchip platforms. The support of computational chemistry methods which appropriately complement the interpretation framework of experimental results is also important. This text describes recent directions in electrochemical methods for the determination of proteins and briefly summarizes available methodologies for the selective labeling of proteins using redox-active probes. Attention is also paid to the theoretical aspects of electron transport and the effect of an external electric field on the structure of selected proteins. Instead of providing a comprehensive overview, we aim to highlight areas of interest that have not been summarized recently, but, at the same time, represent current trends in the field.

Conformational Fluctuations and Induced Orientation of a Protein, Its Solvation Shell, and Bulk Water in Weak Non-Unfolding External Electric Fields

Extreme external electric fields have been reported to disrupt the tertiary structure of stably folded proteins; however, the effects of weaker electric fields on many biomolecules, especially net-uncharged proteins, and on the surrounding aqueous environment have been rarely discussed. To explore these effects at the atomic level, here, we have used molecular dynamics simulations to estimate rotational motion and induced structural fluctuations in the model protein ubiquitin and its hydration layer due to applied non-unfolding electrostatic fields. When exposed to weak electric fields of up to 0.2 V nm-1, ubiquitin displayed competition between internal structure-maintaining molecular interactions and the external orienting force, which disrupted the local structure in certain regions of the protein. Moreover, relative to hydration water, bulk water showed a greater tendency to align with the electric field, indicating that the presence of protein caused hydration water to acquire rotational mobility different from that in a pure-water system. The differential influence of the applied electric field on the hydration and bulk water surrounding ubiquitin will be common to almost all (nonmembrane) biomacromolecules. Our findings highlight the importance of local dipoles and their electric polarizability even in net-uncharged biomolecules.

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.

Coherent diffractive imaging of proteins and viral capsids: simulating MS SPIDOC

MS SPIDOC is a novel sample delivery system designed for single (isolated) particle imaging at X-ray Free-Electron Lasers that is adaptable towards most large-scale facility beamlines. Biological samples can range from small proteins to MDa particles. Following nano-electrospray ionization, ionic samples can be m/z-filtered and structurally separated before being oriented at the interaction zone. Here, we present the simulation package developed alongside this prototype. The first part describes how the front-to-end ion trajectory simulations have been conducted. Highlighted is a quadrant lens; a simple but efficient device that steers the ion beam within the vicinity of the strong DC orientation field in the interaction zone to ensure spatial overlap with the X-rays. The second part focuses on protein orientation and discusses its potential with respect to diffractive imaging methods. Last, coherent diffractive imaging of prototypical T = 1 and T = 3 norovirus capsids is shown. We use realistic experimental parameters from the SPB/SFX instrument at the European XFEL to demonstrate that low-resolution diffractive imaging data (q < 0.3 nm^-1) can be collected with only a few X-ray pulses. Such low-resolution data are sufficient to distinguish between both symmetries of the capsids, allowing to probe low abundant species in a beam if MS SPIDOC is used as sample delivery.

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.

Electric Field Effect on Inhibiting the Co-fibrillation of Amyloid Peptides by Modulating the Aggregation Pathway

With the revelation of the close link between Alzheimer's disease (AD) and type II diabetes (T2D) and the possible assembly of multiple amyloid peptides therein, it is critical to understand and regulate the co-fibrillation pathway between related amyloid peptides. Here, we show experimentally and theoretically that electric field (EF) inhibited hybrid amyloid fibrillation of β-amyloid peptide (Aβ) and human islet amyloid peptide (hIAPP) by modulating the hetero-aggregation pathway. Experimental results confirm that the β-sheet secondary structure of amyloid peptides would be disrupted under small static EF and accompanied by transforming fibril aggregates into amorphous particles in vitro. Molecular dynamics simulations further demonstrate that even with the transformation of the secondary structure from β-sheet to random coil, the strong interaction between Aβ and hIAPP peptides would remain largely unaffected under the small static EF, leading to the formation of amorphous nanoparticles observed in the experiments. This inhibitory effect of EF on the co-fibrillation of multiple amyloid peptides might contribute to reducing the mutual deterioration of different degenerative diseases and show great potential for the noninvasive treatment of amyloid-related diseases.

Electromagnetic bioeffects: a multiscale molecular simulation perspective

Electromagnetic bioeffects remain an enigma from both the experimental and theoretical perspectives despite the ubiquitous presence of related technologies in contemporary life. Multiscale computational modelling can provide valuable insights into biochemical systems and predict how they will be perturbed by external stimuli. At a microscopic level, it can be used to determine what (sub)molecular scale reactions various stimuli might induce; at a macroscopic level, it can be used to examine how these changes affect dynamic behaviour of essential molecules within the crowded biomolecular milieu in living tissues. In this review, we summarise and evaluate recent computational studies that examined the impact of externally applied electric and electromagnetic fields on biologically relevant molecular systems. First, we briefly outline the various methodological approaches that have been employed to study static and oscillating field effects across different time and length scales. The practical value of such modelling is then illustrated through representative case-studies that showcase the diverse effects of electric and electromagnetic field on the main physiological solvent - water, and the essential biomolecules - DNA, proteins, lipids, as well as some novel biomedically relevant nanomaterials. The implications and relevance of the theoretical multiscale modelling to practical applications in therapeutic medicine are also discussed. Finally, we summarise ongoing challenges and potential opportunities for theoretical modelling to advance the current understanding of electromagnetic bioeffects for their modulation and/or beneficial exploitation in biomedicine and industry.