Monte Carlo simulation of a statistical mechanical model of multiple protein sequence alignment

Area of Science:

• Computational biology
• Bioinformatics
• Statistical mechanics

Background:

• Multiple protein sequence alignment is crucial for understanding protein function and evolution.
• Existing models often struggle to incorporate long-range interactions and variable-length insertions effectively.
• The lattice gas model (LGM) offers a framework for sequence analysis, but requires robust simulation methods.

Purpose of the Study:

• To present a grand canonical Monte Carlo (MC) algorithm for the lattice gas model (LGM) of multiple protein sequence alignment.
• To combine long-range interactions and variable-length insertions within a unified computational framework.
• To explore sequence subspaces and assess model parameters using MC simulations.

Main Methods:

• Development and application of a grand canonical Monte Carlo (MC) algorithm.
• Parameter optimization and production runs using MC simulations.
• Comparison of MC simulation results with the mean-field approximation.

Main Results:

• The MC algorithm successfully integrates long-range interactions and variable-length insertions in the LGM.
• Preliminary MC simulations reveal a potential two-state transition in the sequence space for the SH3 domain family.
• The mean-field approximation was found to be inappropriate for the LGM in this context.

Conclusions:

• The developed MC algorithm provides a powerful tool for studying complex protein sequence alignments.
• The findings suggest a phase transition in protein sequence space, offering new insights into evolutionary dynamics.
• The study underscores the limitations of simplified approximations in modeling intricate biological systems.