Water self-diffusion and the dipolar response of the selectivity filter within human aquaporin 4 have been studied using molecular dynamics (MD) simulations in the absence and presence of pulses of external static and alternating electric fields. The pulses were approximately 50 and 100 ns in duration and 0.0065 V/Å in (r.m.s.) intensity and were either static or else 2.45 or 100 GHz in frequency and applied both along and perpendicular to the channels. In addition, the relaxation of the aquaporin, water self-diffusion and gating dynamics following cessation of the impulses was studied. In previous work it was determined that switches in the dihedral angle of the selectivity filter led to boosting of water permeation events within the channels, in the presence of identical external static and alternating electric fields, although applied continuously. Here the application of field impulses (and subsequently, upon removal) has shown that it is the dipolar orientation of the histidine-201 residue in the selectivity filter which governs the dihedral angle, and hence influences water self-diffusion; this constitutes an appropriate order parameter. The dipolar response of this residue to the applied field leads to the adoption of four distinct states, which we modelled as time-homogeneous Markov jump processes, and may be distinguished in the potential of mean force (PMF) as a function of the dipolar orientation of histidine-201. The observations of enhanced "dipolar flipping" of H201 serve to explain increased levels of water self-diffusion within aquaporin channels during, and immediately following, field impulses, although the level of statistical certainty here is lower. Given the appreciable size of the energy barriers evident in PMFs computed directly from deterministic MD (whether in the absence or presence of external fields), metadynamics calculations were undertaken to explore the free-energy landscape of histidine-201 orientation with greater accuracy and precision. These indicate that electric fields do alter the free-energy profile of the H201 side-chain orientation, wherein a perturbation of the symmetric bimodal state evident in the zero-field case is observed. These effects are dependent on the field intensities.

Protein folding and unfolding under the effect of exogenous perturbations remains a topic of great interest, further enhanced by recent technological developments in the field of signal generation that allow the use of intense ultrashort electric pulses to directly interact at microscopic level with biological matter. In this paper, we show results from molecular dynamics (MD) simulations of a single myoglobin molecule in water exposed to pulsed and static electric fields, ranging from 10(8) to 10(9) V/m, compared to data with unexposed conditions. We have found that the highest intensity (10(9) V/m) produced a fast transition (occurring within a few hundreds of picoseconds) between folded and unfolded states, as inferred by secondary structures and geometrical analysis. Fields of 10(8) V/m, on the contrary, produced no significant denaturation, although a relevant effect on the protein dipolar behavior was detected.

The understanding of possible nonthermal bio-effects has been an open question during the last five decades. In this paper, the authors present a critical literature review of the models of the interaction mechanisms, together with an overview of all the publications finding positive results for in vitro and in vivo studies. The systematic approach consisted of pooling together the positive studies on the basis of the endpoints and the biological systems, to identify specific plausible targets of the action of the electromagnetic fields and the related pathways. Such a classification opens the way to the discussion of some hypotheses of interaction mechanisms considered as first transduction step. The authors conclude that only through a multiscale methodology it is possible to perform a comprehensive study of the nonthermal effects, based on affordable and realistic in silico models

Intense pulsed electric fields are known to act at the cell membrane level and are already being exploited in biomedical and biotechnological applications. However, it is not clear if electric pulses within biomedically-attainable parameters could directly influence intra-cellular components such as cytoskeletal proteins. If so, a molecular mechanism of action could be uncovered for therapeutic applications of such electric fields. To help clarify this question, we first identified that a tubulin heterodimer is a natural biological target for intense electric fields due to its exceptional electric properties and crucial roles played in cell division. Using molecular dynamics simulations, we then demonstrated that an intense - yet experimentally attainable - electric field of nanosecond duration can affect the bβ-tubulin’s C-terminus conformations and also influence local electrostatic properties at the GTPase as well as the binding sites of major tubulin drugs site. Our results suggest that intense nanosecond electric pulses could be used for physical modulation of microtubule dynamics. Since a nanosecond pulsed electric field can penetrate the tissues and cellular membranes due to its broadband spectrum, our results are also potentially significant for the development of new therapeutic protocols.

Non-equilibrium molecular dynamics simulations of liquid water have been performed at 298 K in the presence of external time-varying electric fields, approximating a square wave, of varying peak intensity (0.005-0.1 V/A) in the microwave to far-infrared frequency range (20-500 GHz). Significant non-thermal field effects were noted in terms of dipolar response and acceleration of hydrogen-bond kinetics. The coupling between the total dipole moment and the external field has been investigated and autocorrelation functions (ACFs) of both the total dipole moment and the average of the individual molecular dipole moment along the laboratory axis of the applied fields exhibited coupling, with the former showing a stronger coupling and the latter showing coupling to lower magnitude fields. The maximum alignment achieved has been computed as a function of field intensities and frequencies: the lower frequencies show a greater maximum alignment as the system had more time within each field cycle to respond. The normalised probability distribution and the hydrogen-bond ACFs have been computed: the ACF showed a clear effect over shortening the hydrogen-bond relaxation time. The field effects over the molecules' transitions from four to five hydrogen bonds have been computed. There was an enhancement of fewer molecules undergoing transitions and a dampening for a larger proportion of molecules, depending on the external fields' periods

In this paper we present a computational approach, based on NVT molecular dynamics trajectories, that allows the direct evaluation of the protein partial molecular volume. The results obtained for five different globular proteins demonstrate the accuracy of this computational procedure in reproducing protein partial molecular volumes, providing quantitative characterization of the hydration shell in terms of the protein excluded volume, hydration shell ellipsoidal volume and related solvent density. Remarkably, our data indicate for the hydration shell a ≈10% solvent density increase with respect to the liquid water bulk density, in excellent agreement with the available experimental data.

No abstract

Water has many intriguing and anomalous physical properties that have puzzled and titillated the scientific community for centuries, perhaps none more so than the proposition that water may retain some (permanent) “memory” of conditions (e.g., dilution) or electric fields to which it has been subject. Here, we have performed non-equilibrium molecular dynamics simulations of liquid water in external electric-field nanosecond pulses, at 260–310 K, and gauged significant non-thermal field effects in terms of dipolar response. Response of both system- and individual-dipoles has been investigated, and autocorrelation functions of both show more significant effects in stronger fields, with more sluggish relaxation. Crucially, we show that once the field is removed, the dipoles relax, exhibiting no memory or permanent dipolar alignment. We also quantify the time scales for system dynamical-dipolar properties to revert to zero-field equilibrium behaviour.