Accurate free energy calculation along optimized paths.

Area of Science:

• Computational Chemistry
• Molecular Dynamics
• Biophysics

Background:

• Path-based free energy calculation methods (e.g., thermodynamic integration, free energy perturbation) are theoretically simple but practically challenging.
• Difficulty arises from the lack of smooth transition paths, particularly for large molecules.

Purpose of the Study:

• To present a novel method for constructing smooth and short transition paths for peptide free energy calculations.
• To validate the efficiency and accuracy of this new method compared to conventional approaches.

Main Methods:

• Utilizes harmonic potentials to restrain nonhydrogen atom dihedrals in initial and final states.
• Employs a series of geometrical optimization steps to build the transition path.
• Applies the method to a 10-ALA peptide and the beta-hairpin trpzip2.

Main Results:

• Successfully constructs smooth and short transition paths for peptide systems.
• Calculated free energy changes for helix-helix and helix-hairpin transitions were self-convergent and cross-convergent.
• Demonstrates higher efficiency than conventional molecular dynamics for accurate free energy calculations.

Conclusions:

• The novel path construction method offers an efficient and accurate approach for free energy calculations in molecular systems.
• This method overcomes limitations of traditional path-based techniques, especially for complex molecular structures.