The electrostatically driven Monte Carlo method: application to conformational analysis of decaglycine

Area of Science:

• Computational chemistry
• Biophysics
• Molecular modeling

Background:

• Peptide conformation is crucial for biological function.
• Understanding conformational preferences requires accurate energy models and simulation methods.
• The Empirical Conformational Energy Program for Peptides (ECEPP/2) is a widely used model for peptide energy calculations.

Purpose of the Study:

• To investigate the conformational behavior of a glycine decamer using the Electrostatically Driven Monte Carlo (EDMC) method.
• To determine the influence of end-group modifications on the stability of alpha-helical conformations.
• To assess the effectiveness of the EDMC method as a global energy optimizer for peptide systems.

Main Methods:

• Application of the Electrostatically Driven Monte Carlo (EDMC) method.
• Utilizing the Empirical Conformational Energy Program for Peptides (ECEPP/2) potential energy model.
• Simulation of a decamer of glycine with both free and blocked N- and C-termini.

Main Results:

• Unblocked glycine decamers showed lower potential energies in non-alpha-helical conformations.
• Blocking N- and C-termini with acetyl and methyl amide groups reduced unsatisfied hydrogen bonds at helix termini.
• Blocked termini favored alpha-helical conformations as the lowest energy state due to enhanced hydrogen bonding.

Conclusions:

• End-group modification significantly impacts peptide conformational energy landscapes.
• The EDMC method effectively identifies low-energy conformations and avoids local minima.
• Optimizing hydrogen bonding through terminus blocking can stabilize alpha-helical structures in peptides.