Double strand breaks (DSB) in the DNA backbone are the most lethal type of defect induced in the cell nucleus by chemical and radiation treatments of cancer. However, little is known about the outcomes of damage in nucleosomal DNA, and on its effects on damage repair. We performed microsecond-long molecular dynamics computer simulations of nucleosomes including a DSB at various sites, to characterize the early stages of the evolution of this DNA lesion. The damaged structures are studied by the essential dynamics of DNA and histones, and compared to the intact nucleosome, thus exposing key features of the interactions. All DSB configurations tend to remain compact, with only the terminal bases interacting with histone proteins. Umbrella sampling calculations show that broken DNA ends at the DSB must overcome a free-energy barrier to detach from the nucleosome core. Finally, by calculating the covariant mechanical stress, we demonstrate that the coupled bending and torsional stress can force the DSB free ends to open up straight, thus making it accessible to damage signalling proteins.

Maturation of vertebrate oocytes into haploid gametes relies on two consecutive meioses without intervening DNA replication. The temporal sequence of cellular transitions driving eggs from G2 arrest to meiosis I (MI) and then to meiosis II (MII) is controlled by the interplay between cyclin-dependent and mitogen-activated protein kinases. In this paper, we propose a dynamical model of the molecular network that orchestrates maturation of Xenopus laevis oocytes. Our model reproduces the core features of maturation progression, including the characteristic non-monotonous time course of cyclin-Cdks, and unveils the network design principles underlying a precise sequence of meiotic decisions, as captured by bifurcation and sensitivity analyses. Firstly, a coherent and sharp meiotic resumption is triggered by the concerted action of positive feedback loops post-translationally activating cyclin-Cdks. Secondly, meiotic transition is driven by the dynamic antagonism between positive and negative feedback loops controlling cyclin turnover. Our findings reveal a highly modular network in which the coordination of distinct regulatory schemes ensures both reliable and flexible cell-cycle decisions.

Electrostatic interactions play a crucial role in many biomolecular processes, including molecular recognition and binding. Biomolecular electrostatics is modulated to a large extent by the water surrounding the molecules. Here, we present a novel approach to the computation of electrostatic potentials which allows the inclusion of water structure into the classical theory of continuum electrostatics. Based on our recent purely differential formulation of nonlocal electrostatics [Hildebrandt, et al. (2004) Phys. Rev. Lett., 93, 108104] we have developed a new algorithm for its efficient numerical solution. The key component of this algorithm is a boundary element solver, having the same computational complexity as established boundary element methods for local continuum electrostatics. This allows, for the first time, the computation of electrostatic potentials and interactions of large biomolecular systems immersed in water including effects of the solvent's structure in a continuum description. We illustrate the applicability of our approach with two examples, the enzymes trypsin and acetylcholinesterase. The approach is applicable to all problems requiring precise prediction of electrostatic interactions in water, such as protein-ligand and protein-protein docking, folding and chromatin regulation. Initial results indicate that this approach may shed new light on biomolecular electrostatics and on aspects of molecular recognition that classical local electrostatics cannot reveal.