19 research outputs found
Expressions for forces and torques in molecular simulations using rigid bodies
Expressions for intermolecular forces and torques, derived from pair
potentials between rigid non-spherical units, are presented. The aim is to give
compact and clear expressions, which are easily generalised, and which minimise
the risk of error in writing molecular dynamics simulation programs. It is
anticipated that these expressions will be useful in the simulation of liquid
crystalline systems, and in coarse-grained modelling of macromolecules
Scaling in Complex Systems: Analytical Theory of Charged Pores
In this paper we find an analytical solution of the equilibrium ion
distribution for a toroidal model of a ionic channel, using the Perfect
Screening Theorem (PST). The ions are charged hard spheres, and are treated
using a variational Mean Spherical Approximation (VMSA) .
Understanding ion channels is still a very open problem, because of the many
exquisite tuning details of real life channels. It is clear that the electric
field plays a major role in the channel behaviour, and for that reason there
has been a lot of work on simple models that are able to provide workable
theories. Recently a number of interesting papers have appeared that discuss
models in which the effect of the geometry, excluded volume and non-linear
behaviour is considered.
We present here a 3D model of ionic channels which consists of a charged,
deformable torus with a circular or elliptical cross section, which can be flat
or vertical (close to a cylinder). Extensive comparisons to MC simulations were
performed.
The new solution opens new possibilities, such as studying flexible pores,
and water phase transformations inside the pores using an approach similar to
that used on flat crystal surfaces
Average and extreme multi-atom Van der Waals interactions: Strong coupling of multi-atom Van der Waals interactions with covalent bonding
<p>Abstract</p> <p>Background</p> <p>The prediction of ligand binding or protein structure requires very accurate force field potentials – even small errors in force field potentials can make a 'wrong' structure (from the billions possible) more stable than the single, 'correct' one. However, despite huge efforts to optimize them, currently-used all-atom force fields are still not able, in a vast majority of cases, even to keep a protein molecule in its native conformation in the course of molecular dynamics simulations or to bring an approximate, homology-based model of protein structure closer to its native conformation.</p> <p>Results</p> <p>A strict analysis shows that a specific coupling of multi-atom Van der Waals interactions with covalent bonding can, in extreme cases, increase (or decrease) the interaction energy by about 20–40% at certain angles between the direction of interaction and the covalent bond. It is also shown that on average multi-body effects decrease the total Van der Waals energy in proportion to the square root of the electronic component of dielectric permittivity corresponding to dipole-dipole interactions at small distances, where Van der Waals interactions take place.</p> <p>Conclusion</p> <p>The study shows that currently-ignored multi-atom Van der Waals interactions can, in certain instances, lead to significant energy effects, comparable to those caused by the replacement of atoms (for instance, C by N) in conventional pairwise Van der Waals interactions.</p