Effective energy function for proteins in lipid membranes.

Area of Science:

• Biophysics
• Computational Chemistry
• Molecular Modeling

Background:

• Accurate modeling of biomolecules requires considering heterogeneous membrane-aqueous environments.
• Existing energy functions often lack sufficient representation for these complex systems.

Purpose of the Study:

• To extend the EEF1 energy function for improved simulations in heterogeneous membrane-aqueous media.
• To develop a robust computational model for predicting protein behavior at membrane interfaces.

Main Methods:

• Developed solvation parameters for nonpolar phases using experimental amino acid transfer data.
• Introduced a coordinate-dependent solvation free energy and modified dielectric screening for membrane environments.
• Adjusted the electrostatic model to match experimental data for peptides and protein dimers.
• Performed 1-ns molecular dynamics simulations with the new IMM1 energy function.

Main Results:

• The IMM1 function accurately reproduces the preference of tryptophan (Trp) and tyrosine (Tyr) for membrane interfaces.
• It provides reasonable energies for protein insertion and adsorption onto membranes.
• Stable molecular dynamics simulations of the glycophorin A dimer were achieved.
• Melittin's lowest-energy orientation in bilayers was found to be dependent on hydrocarbon layer thickness.

Conclusions:

• The IMM1 energy function offers a significant improvement for modeling systems in membrane-aqueous environments.
• This model facilitates more accurate predictions of protein-membrane interactions and conformational dynamics.
• The findings highlight the importance of membrane thickness in determining peptide orientation.