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.

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.

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/Å) 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.