We develop a numerical method to simulate mechanical objects in a viscous medium at a scale where inertia is negligible. Fibers, spheres and other voluminous objects are represented with points. Different types of connections are used to link the points together and in this way create composite mechanical structures. The motion of such structures in a Brownian environment is described by a first-order multivariate Langevin equation. We propose a computationally efficient method to integrate the equation, and illustrate the applicability of the method to cytoskeletal modeling with several examples.

Based on model structures with atomic resolution, a coarse-grained model for the nucleosome geometry was implemented. The dependence of the chromatin fiber conformation on the spatial orientation of nucleosomes and the path and length of the linker DNA was systematically explored by Monte Carlo simulations. Two fiber types were analyzed in detail that represent nucleosome chains without and with linker histones, respectively: two-start helices with crossed-linker DNA (CL conformation) and interdigitated one-start helices (ID conformation) with different nucleosome tilt angles. The CL conformation was derived from a tetranucleosome crystal structure that was extended into a fiber. At thermal equilibrium, the fiber shape persisted but relaxed into a structure with a somewhat lower linear mass density of 3.1 +/- 0.1 nucleosomes/11 nm fiber. Stable ID fibers required local nucleosome tilt angles between 40 degrees and 60 degrees. For these configurations, much higher mass densities of up to 7.9 +/- 0.2 nucleosomes/11 nm fiber were obtained. A model is proposed, in which the transition between a CL and ID fiber is mediated by relatively small changes of the local nucleosome geometry. These were found to be in very good agreement with changes induced by linker histone H1 binding as predicted from the high resolution model structures.

The cytoskeleton is essential for the maintenance of cell morphology in eukaryotes. In fission yeast, for example, polarized growth sites are organized by actin, whereas microtubules (MTs) acting upstream control where growth occurs. Growth is limited to the cell poles when MTs undergo catastrophes there and not elsewhere on the cortex. Here, we report that the modulation of MT dynamics by forces as observed in vitro can quantitatively explain the localization of MT catastrophes in Schizosaccharomyces pombe. However, we found that it is necessary to add length-dependent catastrophe rates to make the model fully consistent with other previously measured traits of MTs. We explain the measured statistical distribution of MT–cortex contact times and re-examine the curling behavior of MTs in unbranched straight tea1Δ cells. Importantly, the model demonstrates that MTs together with associated proteins such as depolymerizing kinesins are, in principle, sufficient to mark the cell poles.