Learned mappings for targeted free energy perturbation between peptide conformations
Soohaeng Yoo Willow1, Lulu Kang2, David D L Minh3
1Department of Chemistry, Illinois Institute of Technology, Chicago, Illinois 60616, USA.
The Journal of Chemical Physics
|December 21, 2023
Summary
Machine learning enhances free energy calculations by training neural networks to map between thermodynamic states. This approach accurately estimates deca-alanine
Area of Science:
- Computational chemistry
- Statistical mechanics
- Machine learning applications
Background:
- Targeted free energy perturbation methods rely on invertible mappings for accurate free energy estimates.
- Developing effective mappings for complex systems remains a significant challenge in computational chemistry.
- Prior work demonstrated machine learning's potential for mapping between thermodynamic states.
Purpose of the Study:
- To adapt a machine learning approach for calculating free energy differences of a flexible molecule.
- To investigate the efficacy of neural network mappings with harmonic biases and varying spring centers.
- To assess the accuracy and limitations of this method for different thermodynamic state separations.
Main Methods:
- Utilized a deep neural network approach to train mappings between Boltzmann distributions of different thermodynamic states.
- Applied harmonic biases with adjustable spring centers to the deca-alanine molecule.
- Employed "early stopping" criteria during neural network training based on test set loss.
- Calculated free energy differences and compared results with reference methods.
Main Results:
- Accurate free energy differences were obtained for thermodynamic states with spring centers separated by up to 2 Å.
- The neural network mapping successfully promoted configuration space overlap for closely spaced states.
- For more distant states, the mapping failed to generate representative structures, limiting accuracy.
Conclusions:
- Machine learning-based mappings offer a promising avenue for improving free energy calculations in molecular systems.
- "Early stopping" is a crucial technique for preventing overfitting and ensuring reliable results.
- The method's applicability is dependent on the separation between thermodynamic states, with limitations for highly divergent states.
Related Concept Videos
Conserved Binding Sites
4.2K
Many proteins’ biological role depends on their interactions with their ligands, small molecules that bind to specific locations on the protein known as ligand-binding sites. Ligand-binding sites are often conserved among homologous proteins as these sites are critical for protein function.
Binding sites are often located in large pockets, and if their location on a protein’s surface is unknown, it can be predicted using various approaches. The energetic method computationally...
Binding sites are often located in large pockets, and if their location on a protein’s surface is unknown, it can be predicted using various approaches. The energetic method computationally...
4.2K
Protein Folding
8.0K
Proteins are chains of amino acids linked together by peptide bonds. Upon synthesis, a protein folds into a three-dimensional conformation, critical to its biological function. Interactions between its constituent amino acids guide protein folding, and hence the protein structure is primarily dependent on its amino acid sequence.
Protein Structure Is Critical to Its Biological Function
Proteins perform a wide range of biological functions such as catalyzing chemical reactions, providing...
Protein Structure Is Critical to Its Biological Function
Proteins perform a wide range of biological functions such as catalyzing chemical reactions, providing...
8.0K
Protein-protein Interfaces
12.5K
Many proteins form complexes to carry out their functions, making protein-protein interactions (PPIs) essential for an organism's survival. Most PPIs are stabilized by numerous weak noncovalent chemical forces. The physical shape of the interfaces determines the way two proteins interact. Many globular proteins have closely-matching shapes on their surfaces, which form a large number of weak bonds. Additionally, many PPIs occur between two helices or between a surface cleft and a...
12.5K
Cooperative Allosteric Transitions
2.3K
2.3K
Noncovalent Attractions in Biomolecules
51.1K
Noncovalent attractions are associations within and between molecules that influence the shape and structural stability of complexes. These interactions differ from covalent bonding in that they do not involve sharing of electrons.
Four types of noncovalent interactions are hydrogen bonds, van der Waals forces, ionic bonds, and hydrophobic interactions.
Hydrogen bonding results from the electrostatic attraction of a hydrogen atom covalently bonded to a strong-electronegative atom like oxygen,...
Four types of noncovalent interactions are hydrogen bonds, van der Waals forces, ionic bonds, and hydrophobic interactions.
Hydrogen bonding results from the electrostatic attraction of a hydrogen atom covalently bonded to a strong-electronegative atom like oxygen,...
51.1K
Protein Organization
6.5K
Proteins are polymers of amino acid residues. They are versatile and responsible for different cellular functions, including DNA replication, molecular transport, catalysis, and structural support. Proteins have a hierarchical structure comprising at least three levels of organization: primary, secondary, and tertiary structure. Some large proteins have a quaternary structure where individual protein subunits are linked together.
The primary structure of a protein is its amino acid sequence....
The primary structure of a protein is its amino acid sequence....
6.5K


