Molecular dynamics simulations of proteins: can the explicit water model be varied?

In molecular mechanics simulations of biological systems, the solvation water is typically represented by a default water model which is an integral part of the force field. Indeed, protein nonbonding parameters are chosen in order to obtain a balance between water-water and protein-water interactions and hence a reliable description of protein solvation. However, less attention has been paid to the question of whether the water model provides a reliable description of the water properties under the chosen simulation conditions, for which more accurate water models often exist. Here we consider the case of the CHARMM protein force field, which was parametrized for use with a modified TIP3P model. Using quantum mechanical and molecular mechanical calculations, we investigate whether the CHARMM force field can be used with other water models: TIP4P and TIP5P. Solvation properties of N-methylacetamide (NMA), other small solute molecules, and a small protein are examined. The results indicate differences in binding energies and minimum energy geometries, especially for TIP5P, but the overall description of solvation is found to be similar for all models tested. The results provide an indication that molecular mechanics simulations with the CHARMM force field can be performed with water models other than TIP3P, thus enabling an improved description of the solvent water properties.

Testing isotherm models and recovering empirical relationships for adsorption in microporous carbons using virtual carbon models and grand canonical Monte Carlo simulations

Using the plausible model of activated carbon proposed by Harris and co-workers and grand canonical Monte Carlo simulations, we study the applicability of standard methods for describing adsorption data on microporous carbons widely used in adsorption science. Two carbon structures are studied, one with a small distribution of micropores in the range up to 1 nm, and the other with micropores covering a wide range of porosity. For both structures, adsorption isotherms of noble gases (from Ne to Xe), carbon tetrachloride and benzene are simulated. The data obtained are considered in terms of Dubinin-Radushkevich plots. Moreover, for benzene and carbon tetrachloride the temperature invariance of the characteristic curve is also studied. We show that using simulated data some empirical relationships obtained from experiment can be successfully recovered. Next we test the applicability of Dubinin's related models including the Dubinin-Izotova, Dubinin-Radushkevich-Stoeckli, and Jaroniec-Choma equations. The results obtained demonstrate the limits and applications of the models studied in the field of carbon porosity characterization.

Molecular modelling studies of binding of DACD derivatives into G-Quadruplex DNA: comparison of force field and quantum polarized ligand docking methods

The DNA G-qadruplexes are one of the targets being actively explored for anti-cancer therapy by inhibiting them through small molecules. This computational study was conducted to predict the binding strengths and orientations of a set of novel dimethyl-amino-ethyl-acridine (DACA) analogues that are designed and synthesized in our laboratory, but did not diffract in Synchrotron light.Thecrystal structure of DNA G-Quadruplex(TGGGGT)4(PDB: 1O0K) was used as target for their binding properties in our studies.We used both the force field (FF) and QM/MM derived atomic charge schemes simultaneously for comparing the predictions of drug binding modes and their energetics. This study evaluates the comparative performance of fixed point charge based Glide XP docking and the quantum polarized ligand docking schemes. These results will provide insights on the effects of including or ignoring the drug-receptor interfacial polarization events in molecular docking simulations, which in turn, will aid the rational selection of computational methods at different levels of theory in future drug design programs. Plenty of molecular modelling tools and methods currently exist for modelling drug-receptor or protein-protein, or DNA-protein interactionssat different levels of complexities.Yet, the capasity of such tools to describevarious physico-chemical propertiesmore accuratelyis the next step ahead in currentresearch.Especially, the usage of most accurate methods in quantum mechanics(QM) is severely restricted by theirtedious nature. Though the usage of massively parallel super computing environments resulted in a tremendous improvement in molecular mechanics (MM) calculations like molecular dynamics,they are still capable of dealing with only a couple of tens to hundreds of atoms for QM methods. One such efficient strategy that utilizes thepowers of both MM and QM are the QM/MM hybrid methods. Lately, attempts have been directed towards the goal of deploying several different QM methods for betterment of force field based simulations, but with practical restrictions in place. One of such methods utilizes the inclusion of charge polarization events at the drug-receptor interface, that is not explicitly present in the MM FF.

Connecting mean field models of neural activity to EEG and fMRI data

Progress in functional neuroimaging of the brain increasingly relies on the integration of data from complementary imaging modalities in order to improve spatiotemporal resolution and interpretability. However, the usefulness of merely statistical combinations is limited, since neural signal sources differ between modalities and are related non-trivially. We demonstrate here that a mean field model of brain activity can simultaneously predict EEG and fMRI BOLD with proper signal generation and expression. Simulations are shown using a realistic head model based on structural MRI, which includes both dense short-range background connectivity and long-range specific connectivity between brain regions. The distribution of modeled neural masses is comparable to the spatial resolution of fMRI BOLD, and the temporal resolution of the modeled dynamics, importantly including activity conduction, matches the fastest known EEG phenomena. The creation of a cortical mean field model with anatomically sound geometry, extensive connectivity, and proper signal expression is an important first step towards the model-based integration of multimodal neuroimages.

Modeling electrocortical activity through improved local approximations of integral neural field equations

Neural field models of firing rate activity typically take the form of integral equations with space-dependent axonal delays. Under natural assumptions on the synaptic connectivity we show how one can derive an equivalent partial differential equation (PDE) model that properly treats the axonal delay terms of the integral formulation. Our analysis avoids the so-called long-wavelength approximation that has previously been used to formulate PDE models for neural activity in two spatial dimensions. Direct numerical simulations of this PDE model show instabilities of the homogeneous steady state that are in full agreement with a Turing instability analysis of the original integral model. We discuss the benefits of such a local model and its usefulness in modeling electrocortical activity. In particular, we are able to treat “patchy” connections, whereby a homogeneous and isotropic system is modulated in a spatially periodic fashion. In this case the emergence of a “lattice-directed” traveling wave predicted by a linear instability analysis is confirmed by the numerical simulation of an appropriate set of coupled PDEs.

Temporal evolution and electric potential structure of the auroral acceleration region from multispacecraft measurements

Bright aurorae can be excited by the acceleration of electrons into the atmosphere in violation of ideal magnetohydrodynamics. Modelling studies predict that the accelerating electric potential consists of electric double layers at the boundaries of an acceleration region but observations suggest that particle acceleration occurs throughout this region. Using multi-spacecraft observations from Cluster we have examined two upward current regions on 14 December 2009. Our observations show that the potential difference below C4 and C3 changed by up to 1.7 kV between their respective crossings, which were separated by 150 s. The field-aligned current density observed by C3 was also larger than that observed by C4. The potential drop above C3 and C4 was approximately the same in both crossings. Using a novel technique of quantitatively comparing the electron spectra measured by Cluster 1 and 3, which were separated in altitude, we determine when these spacecraft made effectively magnetically conjugate observations and use these conjugate observations to determine the instantaneous distribution of the potential drop in the AAR. Our observations show that an average of 15% of the potential drop in the AAR was located between C1 at 6235 km and C3 at 4685 km altitude, with a maximum potential drop between the spacecraft of 500~V and that the majority of the potential drop was below C3. By assuming a spatial invariance along the length of the upward current region, we discuss these observations in terms of temporal changes and the vertical structure of the electrostatic potential drop and in the context of existing models and previous observations single- and multi-spacecraft observations.

The effect of the magnetic field of the secondary star in dwarf novae

We investigate the effect of a secondary star magnetic field on the accretion disc dynamics of dwarf novae. Simulations have been carried out with a particle code and a dipolar magnetic field structure. The magnetic field acts to remove angular momentum from the disc material, increasing the inward mass flow. This makes the accretion disc more centrally condensed, causing a reduction in the recurrence time for dwarf nova outbursts. We have produced Doppler tomograms and light curves which may be compared with observations. These tomograms are significantly different from those produced in the absence of a magnetic field on the secondary. We derive an upper limit to the magnetic moment of the secondary star in UGem of mu_2<2x10^32 A m^2. The magnetic truncation of the accretion disc produces resonance phenomena similar to those seen in the superoutbursts of SUUMa systems. While these have not been observed for systems like UGem, observations of the SUUMa systems provide us with a useful diagnostic of the disc-field interaction. We are able to place an upper limit on the magnetic moment of the secondary in ZCha of mu_2<1x10^30 A m^2.

Flow structure and near-field dispersion in arrays of building-like obstacles

Dispersion in the near-field region of localised releases in urban areas is difficult to predict because of the strong influence of individual buildings. Effects include upstream dispersion, trapping of material into building wakes and enhanced concentration fluctuations. As a result, concentration patterns are highly variable in time and mean profiles in the near field are strongly non-Gaussian. These aspects of near-field dispersion are documented by analysing data from direct numerical simulations in arrays of building-like obstacles and are related to the underlying flow structure. The mean flow structure around the buildings is found to exert a strong influence over the dispersion of material in the near field. Diverging streamlines around buildings enhance lateral dispersion. Entrainment of material into building wakes in the very near field gives rise to secondary sources, which then affect the subsequent dispersion pattern. High levels of concentration fluctuations are also found in this very near field; the fluctuation intensity is of order 2 to 5.

Self-consistent wave-particle interactions in dispersive scale long-period field-line-resonances

Using 1D Vlasov drift-kinetic computer simulations, it is shown that electron trapping in long period standing shear Alfven waves (SAWs) provides an efficient energy sink for wave energy that is much more effective than Landau damping. It is also suggested that the plasma environment of low altitude auroral-zone geomagnetic field lines is more suited to electron acceleration by inertial or kinetic scale Alfven waves. This is due to the self-consistent response of the electron distribution function to SAWs, which must accommodate the low altitude large-scale current system in standing waves. We characterize these effects in terms of the relative magnitude of the wave phase and electron thermal velocities. While particle trapping is shown to be significant across a wide range of plasma temperatures and wave frequencies, we find that electron beam formation in long period waves is more effective in relatively cold plasma.

Creep and recovery of magnetorheological fluids: experiments and simulations

A direct comparative study on the creep-recovery behavior of conventional MR fluids is carried out using magnetorheometry and particle-level simulations. Two particle concentrations are investigated (ϕ=0.05 and 0.30) at two different magnetic field strengths (53 kA•m-1 and 173 kA•m-1) in order to match the yield stresses developed in both systems for easier comparison. Simulations are mostly started with random initial structures with some additional tests of using preassembled single chains in the low concentration case. Experimental and simulation data are in good qualitative agreement. The results demonstrate three regions in the creep curves: i) In the initial viscoelastic region, the chain-like (at ϕ=0.05) or percolated three-dimensional network (at ϕ=0.30) structures fill up the gap and the average cluster size remains constant; ii) Above a critical strain of 10 %, in the retardation region, these structures begin to break and rearrange under shear. At large enough imposed stress values, they transform into thin sheet-like or thick lamellar structures, depending on the particle concentration; iii) Finally in the case of larger strain values either the viscosity diverges (at low stress values) or reaches a constant low value (at high stress values), showing a clear bifurcation behavior. For stresses below the bifurcation point the MR fluid is capable to recover the strain by a certain fraction. However, no recovery is observed for large stress values.

Computer simulations of the structure of colloidal ferrofluids

The structure of a ferrofluid under the influence of an external magnetic field is expected to become anisotropic due to the alignment of the dipoles into the direction of the external field, and subsequently to the formation of particle chains due to the attractive head to tail orientations of the ferrofluid particles. Knowledge about the structure of a colloidal ferrofluid can be inferred from scattering data via the measurement of structure factors. We have used molecular-dynamics simulations to investigate the structure of both monodispersed and polydispersed ferrofluids. The results for the isotropic structure factor for monodispersed samples are similar to previous data by Camp and Patey that were obtained using an alternative Monte Carlo simulation technique, but in a different parameter region. Here we look in addition at bidispersed samples and compute the anisotropic structure factor by projecting the q vector onto the XY and XZ planes separately, when the magnetic field was applied along the z axis. We observe that the XY- plane structure factor as well as the pair distribution functions are quite different from those obtained for the XZ plane. Further, the two- dimensional structure factor patterns are investigated for both monodispersed and bidispersed samples under different conditions. In addition, we look at the scaling exponents of structure factors. Our results should be of value to interpret scattering data on ferrofluids obtained under the influence of an external field.

On a resurgence of management simulations and games

The past few years have seen a significant resurgence of interest in ‘management games’ and ‘management flight simulators’, one particularly active source of such work being the system dynamics community. After proposing a distinction between games and simulations, this paper provides some background to these developments by briefly describing the historical roots of the field and the fundamental ideas of the system dynamics community, which are now giving rise to ‘microworlds’. The training advantages of management simulations and games are then discussed. The paper closes with a note on the research and findings of the system dynamics field and by offering some words of warning on the perils of simulation and game use. Two scenarios for how the use of simulations and games as management education devices might develop in the future are proposed. An Appendix describes five examples of very different types of management simulations and games.

The evolution of the sun’s open magnetic flux: II. full solar cycle simulations

In this paper the origin and evolution of the Sun’s open magnetic flux is considered by conducting magnetic flux transport simulations over many solar cycles. The simulations include the effects of differential rotation, meridional flow and supergranular diffusion on the radial magnetic field at the surface of the Sun as new magnetic bipoles emerge and are transported poleward. In each cycle the emergence of roughly 2100 bipoles is considered. The net open flux produced by the surface distribution is calculated by constructing potential coronal fields with a source surface from the surface distribution at regular intervals. In the simulations the net open magnetic flux closely follows the total dipole component at the source surface and evolves independently from the surface flux. The behaviour of the open flux is highly dependent on meridional flow and many observed features are reproduced by the model. However, when meridional flow is present at observed values the maximum value of the open flux occurs at cycle minimum when the polar caps it helps produce are the strongest. This is inconsistent with observations by Lockwood, Stamper and Wild (1999) and Wang, Sheeley, and Lean (2000) who find the open flux peaking 1–2 years after cycle maximum. Only in unrealistic simulations where meridional flow is much smaller than diffusion does a maximum in open flux consistent with observations occur. It is therefore deduced that there is no realistic parameter range of the flux transport variables that can produce the correct magnitude variation in open flux under the present approximations. As a result the present standard model does not contain the correct physics to describe the evolution of the Sun’s open magnetic flux over an entire solar cycle. Future possible improvements in modeling are suggested.

Large plasma velocities along the magnetic field line in the auroral zone

In the auroral zone, ionospheric plasma often moves horizontally at more than 1 km s−1, driven by magnetospheric electric fields, but it has usually been assumed that vertical velocities are much smaller. On occasions, however, plasma has been seen to move upwards along the magnetic field line at several hundred m s−1. These upward velocities are associated with electric fields strong enough to heat the ion population and drive it into a non-thermal state1,2. Here we report observations of substantial upwards acceleration of plasma, to velocities as high as 500 m s−1. An initial upthrust was provided by a combined upwelling of the neutral atmosphere and ionosphere but the continued acceleration at greater heights is explained by a combination of enhanced plasma pressure and the 'hydrodynamic mirror force'3. This acceleration marks an important stage in the transport of plasma from the ionosphere into the magnetosphere.

Emergence of spatially heterogeneous burst suppression in a neural field model of electrocortical activity

Burst suppression in the electroencephalogram (EEG) is a well-described phenomenon that occurs during deep anesthesia, as well as in a variety of congenital and acquired brain insults. Classically it is thought of as spatially synchronous, quasi-periodic bursts of high amplitude EEG separated by low amplitude activity. However, its characterization as a “global brain state” has been challenged by recent results obtained with intracranial electrocortigraphy. Not only does it appear that burst suppression activity is highly asynchronous across cortex, but also that it may occur in isolated regions of circumscribed spatial extent. Here we outline a realistic neural field model for burst suppression by adding a slow process of synaptic resource depletion and recovery, which is able to reproduce qualitatively the empirically observed features during general anesthesia at the whole cortex level. Simulations reveal heterogeneous bursting over the model cortex and complex spatiotemporal dynamics during simulated anesthetic action, and provide forward predictions of neuroimaging signals for subsequent empirical comparisons and more detailed characterization. Because burst suppression corresponds to a dynamical end-point of brain activity, theoretically accounting for its spatiotemporal emergence will vitally contribute to efforts aimed at clarifying whether a common physiological trajectory is induced by the actions of general anesthetic agents. We have taken a first step in this direction by showing that a neural field model can qualitatively match recent experimental data that indicate spatial differentiation of burst suppression activity across cortex.