Dominance of extreme statistics in a prototype many-body Brownian ratchet

Many forms of cell motility rely on Brownian ratchet mechanisms that involve multiple stochastic processes. We present a computational and theoretical study of the nonequilibrium statistical dynamics of such a many-body ratchet, in the specific form of a growing polymer gel that pushes a diffusing obstacle. We find that oft-neglected correlations among constituent filaments impact steady-state kinetics and significantly deplete the gel's density within molecular distances of its leading edge. These behaviors are captured quantitatively by a self-consistent theory for extreme fluctuations in filaments' spatial distribution.Comment: 5 pages with 3 figures + 20 pages of Supplementary Material with 2 figures. Updated to agree with published version; published as a Communication in J. Chem. Phy

Coexistence between fluid and crystalline phases of proteins in photosynthetic membranes

Photosystem II (PSII) and its associated light-harvesting complex II (LHCII) are highly concentrated in the stacked grana regions of photosynthetic thylakoid membranes. Within the membrane, PSII-LHCII supercomplexes can be arranged in disordered packings, ordered arrays, or mixtures thereof. The physical driving forces underlying array formation are unknown, complicating attempts to determine a possible functional role for arrays in regulating light harvesting or energy conversion efficiency. Here we introduce a coarse-grained model of protein interactions in coupled photosynthetic membranes, focusing on just two particle types that feature simple shapes and potential energies motivated by structural studies. Reporting on computer simulations of the model's equilibrium fluctuations, we demonstrate its success in reproducing diverse structural features observed in experiments, including extended PSII-LHCII arrays. Free energy calculations reveal that the appearance of arrays marks a phase transition from the disordered fluid state to a system-spanning crystal, which can easily be arrested by thermodynamic constraints or slow dynamics. The region of fluid-crystal coexistence is broad, encompassing much of the physiologically relevant parameter regime. Our results suggest that grana membranes lie at or near phase coexistence, conferring significant structural and functional flexibility to this densely packed membrane protein system.Comment: 11 pages, 5 figure

Precision shooting: Sampling long transition pathways

The kinetics of collective rearrangements in solution, such as protein folding and nanocrystal phase transitions, often involve free energy barriers that are both long and rough. Applying methods of transition path sampling to harvest simulated trajectories that exemplify such processes is typically made difficult by a very low acceptance rate for newly generated trajectories. We address this problem by introducing a new generation algorithm based on the linear short-time behavior of small disturbances in phase space. Using this ``precision shooting'' technique, arbitrarily small disturbances can be propagated in time, and any desired acceptance ratio of shooting moves can be obtained. We demonstrate the method for a simple but computationally problematic isomerization process in a dense liquid of soft spheres. We also discuss its applicability to barrier crossing events involving metastable intermediate states.Comment: 9 pages, 12 figures, submitted to J. Chem. Phy

Two-stage coarsening mechanism in a kinetically constrained model of an attractive colloid

We study an attractive version of the East model using the real-space renormalization group (RG) introduced by Stella et al. The former is a kinetically constrained model with an Ising-like interaction between excitations, and shows striking agreement with the phenomonology of attractive colloidal systems. We find that the RG predicts two nonuniversal dynamic exponents, which suggests that in the out-of-equilibrium regime the model coarsens via a two-stage mechanism. We explain this mechanism physically, and verify this prediction numerically. In addition, we find that the characteristic relaxation time of the model is a non-monotonic function of attraction strength, again in agreement with numerical results.Comment: 10 page