Abstract:Molecular dynamics (MD) is an extremely powerful, highly effective, and widely used approach to understanding the nature of chemical processes in atomic details for proteins. The accuracy of results from MD simulations is highly dependent on force fields. Currently, molecular mechanical (MM) force fields are mainly utilized in MD simulations because of their low computational cost. Quantum mechanical (QM) calculation has high accuracy, but it is exceedingly time consuming for protein simulations. Machine learning (ML) provides the capability for generating accurate potential at the QM level without increasing much computational effort for specific systems that can be studied at the QM level. However, the construction of general machine learned force fields, needed for broad applications and large and complex systems, is still challenging. Here, general and transferable neural network (NN) force fields based on CHARMM force fields, named CHARMM-NN, are constructed for proteins by training NN models on 27 fragments partitioned from the residue-based systematic molecular fragmentation (rSMF) method. The NN for each fragment is based on atom types and uses new input features that are similar to MM inputs, including bonds, angles, dihedrals, and non-bonded terms, which enhance the compatibility of CHARMM-NN to MM MD and enable the implementation of CHARMM-NN force fields in different MD programs. While the main part of the energy of the protein is based on rSMF and NN, the nonbonded interactions between the fragments and with water are taken from the CHARMM force field through mechanical embedding. The validations of the method for dipeptides on geometric data, relative potential energies, and structural reorganization energies demonstrate that the CHARMM-NN local minima on the potential energy surface are very accurate approximations to QM, showing the success of CHARMM-NN for bonded interactions. However, the MD simulations on peptides and proteins indicate that more accurate methods to represent protein-water interactions in fragments and non-bonded interactions between fragments should be considered in the future improvement of CHARMM-NN, which can increase the accuracy of approximation beyond the current mechanical embedding QM/MM level.

Automatically Constructed Neural Network Potentials for Molecular Dynamics Simulation of Zinc Proteins

Automatically constructed neural network potentials for molecular dynamics simulation of zinc proteins

Molecular Dynamics Simulation of Zinc Ion in Water with an Ab Initio Based Neural Network Potential

Ab Initio Molecular Dynamics Simulation of Zinc Metalloproteins with Enhanced Self-Organizing Incremental High Dimensional Neural Network

Modeling Zinc Complexes Using Neural Networks

Artificial neural network potential for pure zinc

Interaction of divalent cations and amino acids in bulk water: Molecular simulations with neural network potentials

A New Quantum Calibrated Force Field for Zinc-Protein Complex

Development of machine learning interatomic potential for zinc

A Fragment Quantum Mechanical Method for Metalloproteins.

Toward a general neural network force field for protein simulations: Refining the intramolecular interaction in protein

A Transferable Non-bonded Pairwise Force Field to Model Zinc Interactions in Metalloproteins.

A Physics-Aware Neural Network for Protein-Ligand Interactions with Quantum Chemical Accuracy

Neural Network Potential with Multi-Resolution Approach Enables Accurate Prediction of Reaction Free Energies in Solution

Neural network potentials

Modeling of Metal Nanoparticles: Development of Neural-Network Interatomic Potential Inspired by Features of the Modified Embedded-Atom Method

A physics-aware neural network for protein–ligand interactions with quantum chemical accuracy

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Prediction of structures of zinc-binding proteins through explicit modeling of metal coordination geometry.