Abstract:Metalloproteins widely exist in biology and play important roles in various processes. To accurately simulate metalloprotein systems, modeling polarization and charge transfer effects is vital. The fluctuating charge (FQ) model can efficiently generate atomic charges and simulate the charge transfer effect; it has been developed for a wide range of applications, but few models have been specifically tailored for metalloproteins. In this study, we present a fluctuating charge model specifically for zinc-containing metalloproteins based on the extended charge equilibration (EQeq) scheme. Our model was parametrized to reproduce CM5 charges instead of RESP/CHELPG charges because the former is less dependent on the conformation or basis set, does not suffer from unphysical charges for buried atoms, and is still being able to well reproduce the molecular dipoles. During our study, we found that adding the Pauling-bond-order-like term (referred to as the "+C term" in a previous study) between the zinc ion and ligating atoms significantly improves the model's performance. Although our model was trained for four-coordinated zinc sites only, our results indicated it can well describe the atomic charges in diverse zinc sites. Morever, our model was used to generate partial charges for the metal sites in three different zinc-containing metalloproteins (with four-, five-, and six-coordinated metal sites, respectively). These charges exhibited performance comparable to that of the RESP charges in molecular dynamics (MD) simulations. Additional tests indicated our model could also well reproduce the CM5 charges when geometric changes were involved. Those results indicate that our model can efficiently calculate the atomic charges for metal sites and well simulate the charge transfer effect, which marks an important step toward developing versatile polarizable force fields for metalloproteins.

DP/MM: A Hybrid Model for Zinc-Protein Interactions in Molecular Dynamics

A Discontinuous Potential Model for Protein-Protein Interactions.

A New Quantum Calibrated Force Field for Zinc-Protein Complex

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

Structure and dynamics of a dizinc metalloprotein: Effect of charge transfer and polarization

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

Development of a Fluctuating Charge Model for Zinc-Containing Metalloproteins

A Force Balanced Fragmentation Method for Ab Initio Molecular Dynamic Simulation of Protein

QDπ: A Quantum Deep Potential Interaction Model for Drug Discovery

Systematic Derivation of AMBER Force Field Parameters Applicable to Zinc-Containing Systems.

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

Quantum Fragment Based ab Initio Molecular Dynamics for Proteins.

Toward On-The-Fly Quantum Mechanical/Molecular Mechanical (QM/MM) Docking: Development and Benchmark of a Scoring Function

Structural Survey of Zinc Containing Proteins and the Development of the Zinc AMBER Force Field (ZAFF).

A Fragment Quantum Mechanical Method for Metalloproteins.

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

Comparison of Empirical Zn2+ Models in Protein–DNA Complexes

Parameterization of a Fluctuating Charge Model for Complexes Containing 3d Transition Metals

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

A Comparison of Bonded and Nonbonded Zinc(II) Force Fields with NMR Data

Completely Multipolar Model for Many-Body Water-Ion and Ion-Ion Interactions