An extrapolation method for computing protein solvation energies based on density fragmentation of a graphical...

Area of Science:

• Computational chemistry
• Biomolecular modeling
• Quantum mechanics/molecular mechanics (QM/MM)

Background:

• Continuum polarizable models are crucial for simulating solvent effects in proteins.
• Current methods using surface tessellation face computational limitations (matrix inversion) for large systems.

Purpose of the Study:

• To develop a self-consistent approach for conductor-like screening charges in protein simulations.
• To improve the efficiency and scalability of modeling solvent effects in large biomolecular systems.

Main Methods:

• A density fragmentation approach applied to a graphical surface tessellation.
• Implementation with fixed charge and QM/MM (moving domain) protein models.
• Development of a parallelizable scheme for conductor-like models.

Main Results:

• The developed method provides accurate conductor-like screening charges for protein applications.
• The approach is suitable for large systems, overcoming limitations of traditional matrix inversion.
• Integration with QM/MM successfully captures protein self-polarization and solvent polarization effects.

Conclusions:

• The density fragmentation method offers a scalable and efficient way to model solvent effects in proteins.
• This approach facilitates more accurate simulations of chemical events within proteins.
• The method is amenable to parallelization, enhancing computational performance.