Efficient global biopolymer sampling with end-transfer configurational bias Monte Carlo

Area of Science:

• Computational Biology
• Biophysics
• Statistical Mechanics

Background:

• Efficient thermodynamic sampling of complex biopolymers is crucial for understanding their behavior.
• Traditional Monte Carlo methods can struggle with sampling large, complex systems like chromatin.

Purpose of the Study:

• To develop and evaluate a novel "end-transfer configurational bias Monte Carlo" method.
• To assess its efficiency for sampling oligonucleosome models under varying salt conditions.

Main Methods:

• The end-transfer method involves deleting a motif from one end and regrowing it at the other.
• Performance was compared against local moves, pivot rotations, and standard configurational bias.
• Sampling efficiency was measured across translational, rotational, and internal degrees of freedom.

Main Results:

• The end-transfer method showed superior sampling efficiency for all degrees of freedom at high salt concentrations (weak electrostatics).
• It was less effective than pivot rotations for internal and rotational sampling at low-to-moderate salt concentrations (strong electrostatics).
• The end-transfer method was orders of magnitude more efficient than standard configurational bias, with sampling time scaling quadratically vs. exponentially with length.

Conclusions:

• The end-transfer configurational bias Monte Carlo method offers significant performance improvements for global biomolecular simulations.
• It is particularly effective in condensed systems with weak nonbonded interactions.
• The method can be combined with local enhancements for improved overall sampling efficiency.