The previously developed particle mesh Ewald method is reformulated in terms of efficient B‐spline interpolation of the structure factors. This reformulation allows a natural extension of the method to potentials of the form 1/r^p with p≥1. Furthermore, efficient calculation of the virial tensor follows. Use of B‐splines in place of Lagrange interpolation leads to analytic gradients as well as a significant improvement in the accuracy. We demonstrate that arbitrary accuracy can be achieved, independent of system size N, at a cost that scales as N log(N). For biomolecular systems with many thousands of atoms this method permits the use of Ewald summation at a computational cost comparable to that of a simple truncation method of 10 Å or less.

Many‐body effects were examined in a series of molecular dynamics computer simulations on the ionic aqueous clusters Na⁺(H₂O)_n (n=4,5,6,14) and Cl⁻(H₂O)_n (n=4,5,6,7,8,14). Two potential models were used in the simulations. In one model (TIP4P) the potential was pairwise additive, while in the second model (SPCE/POL) the many body effects were explicitly included through a self‐consistent polarization routine. The two models produce equilibrium structures which are significantly different in energy and geometry. The SPCE/POL model consistently predicts energetically more stable products. In addition, for the anion cluster systems the SPCE/POL model places the Cl⁻ on the surface of the water cluster.

Transcription factor-dependent cellular reprogramming is integral to normal development and is central to production of induced pluripotent stem cells. This process typically requires pioneer transcription factors (TFs) to induce de novo formation of enhancers at previously closed chromatin. Mechanistic information on this process is currently sparse. Here we explore the mechanistic basis by which GATA3 functions as a pioneer TF in a cellular reprogramming event relevant to breast cancer, the mesenchymal to epithelial transition (MET). In some instances, GATA3 binds previously inaccessible chromatin, characterized by stable, positioned nucleosomes where it induces nucleosome eviction, alters local histone modifications, and remodels local chromatin architecture. At other loci, GATA3 binding induces nucleosome sliding without concomitant generation of accessible chromatin. Deletion of the transactivation domain retains the chromatin binding ability of GATA3 but cripples chromatin reprogramming ability, resulting in failure to induce MET. These data provide mechanistic insights into GATA3-mediated chromatin reprogramming during MET, and suggest unexpected complexity to TF pioneering. Successful reprogramming requires stable binding to a nucleosomal site; activation domain-dependent recruitment of co-factors including BRG1, the ATPase subunit of the SWI/SNF chromatin remodeling complex; and appropriate genomic context. The resulting model provides a new conceptual framework for de novo enhancer establishment by a pioneer TF.

Molecular dynamics simulations were performed using a modified amoeba force field to determine hydration and dynamical properties of the divalent cations Ca2+ and Mg2+. The extension of amoeba to divalent cations required the introduction of a cation specific parametrization. To accomplish this, the Tholé polarization damping model parametrization was modified based on the ab initio polarization energy computed by a constrained space orbital variation energy decomposition scheme. Excellent agreement has been found with condensed phase experimental results using parameters derived from gas phase ab initio calculations. Additionally, we have observed that the coordination of the calcium cation is influenced by the size of the periodic water box, a recurrent issue in first principles molecular dynamics studies.