Folding Lennard-Jones proteins by a contact potential

Area of Science:

• Computational chemistry
• Biophysics
• Polymer physics

Background:

• Lennard-Jones potentials are fundamental in molecular simulations.
• Accurate modeling of protein folding is crucial for understanding biological function.
• Simplifying complex interactions can accelerate computational studies.

Purpose of the Study:

• To investigate the feasibility of approximating Lennard-Jones interactions with pairwise contact potentials for protein-like heteropolymers.
• To determine if a contact potential can accurately represent the native states of heteropolymer sequences designed with Lennard-Jones potentials.
• To assess the effectiveness of contact potential-based energy minimization for predicting protein structures.

Main Methods:

• Designing off-lattice heteropolymer sequences using a Lennard-Jones potential.
• Identifying native conformations via molecular dynamics simulations.
• Developing and testing pairwise contact potentials to approximate Lennard-Jones interactions.
• Employing energy minimization in the space of contact maps.

Main Results:

• Exact approximation of Lennard-Jones interactions by a pairwise contact potential is not possible.
• Optimized contact energy parameters allow for energy minimization yielding structures close to native ones.
• Using these structures as starting points for molecular dynamics recovers native folds with high probability.

Conclusions:

• Pairwise contact potentials offer a viable, albeit approximate, approach to modeling protein folding.
• This method can efficiently guide simulations towards correct protein structures.
• The findings have implications for coarse-grained modeling and protein design.