Modeling pH-Dependent Protein Dynamics by Integrating Coarse-Grained Molecular Simulation and the Deep Neural Network
Yanhang Liu1,2, Huaqi Peng1,2, Hengyan Huang1,2
1Department of Physics, National Laboratory of Solid State Microstructure, Nanjing University, Nanjing 210093, China.
Journal of Chemical Theory and Computation
|May 4, 2026
Summary
DeepCGpKa predicts protein pKa values using deep learning on coarse-grained structures, matching all-atom accuracy. This advances simulations of protein behavior and pH-dependent dynamics.
Area of Science:
- Computational Biology
- Biophysics
- Structural Biology
Background:
- Protein pKa values are crucial for understanding protein behavior, influencing folding, dynamics, and interactions.
- Current pKa prediction methods require detailed all-atom structures, limiting their use with widely adopted coarse-grained models.
- Estimating pKa in coarse-grained models is challenging due to the lack of explicit ionizable groups.
Purpose of the Study:
- To develop a deep-learning-based pKa predictor for coarse-grained protein structures.
- To enable accurate pKa prediction within coarse-grained simulation frameworks.
- To improve the modeling of electrostatic interactions in coarse-grained biomolecular simulations.
Main Methods:
- Introduction of DeepCGpKa, a novel deep-learning algorithm for pKa prediction.
- Application of DeepCGpKa to coarse-grained protein structures.
- Benchmarking against state-of-the-art all-atom pKa prediction methods.
- Integration with coarse-grained molecular dynamics simulations.
Main Results:
- DeepCGpKa achieves accuracy comparable to existing all-atom prediction methods.
- The predictor demonstrates robust performance on partially unfolded protein structures.
- Coupling DeepCGpKa with simulations effectively captures pH-dependent protein conformational changes.
Conclusions:
- DeepCGpKa offers a powerful tool for pKa prediction in coarse-grained protein models.
- This approach enhances the treatment of electrostatic interactions in coarse-grained simulations.
- The integration of data-driven prediction and physics-based simulation provides a practical solution for coarse-grained biomolecular modeling.
Related Concept Videos
Protein Dynamics in Living Cells
1.9K
Different fluorescence-based techniques are used to study the protein dynamics in living cells. These techniques include FRAP, FRET, and PET.
Fluorescent recovery after photobleaching (FRAP) is a fluorescent-protein-based detection technique used to quantify protein movement rates within the cell. This method exposes a small portion of the cell to an intense laser beam. The laser beam causes permanent photobleaching of the fluorophore-tagged proteins in the exposed region. As the bleached...
Fluorescent recovery after photobleaching (FRAP) is a fluorescent-protein-based detection technique used to quantify protein movement rates within the cell. This method exposes a small portion of the cell to an intense laser beam. The laser beam causes permanent photobleaching of the fluorophore-tagged proteins in the exposed region. As the bleached...
1.9K
Pharmacodynamic Models: Overview
159
Pharmacodynamic (PD) responses describe the interaction between a drug and its biological target, culminating in a physiological effect. These responses can be classified into different types: continuous variables, such as blood glucose levels; categorical outcomes, like survival rates; and time-to-event metrics, such as disease progression. Understanding and modeling PD responses are critical for optimizing drug efficacy and safety.PD models describe the relationship between drug concentration...
159
Physiological Pharmacokinetic Models: Assumption with Protein Binding
378
Physiological models with protein binding in pharmacokinetics offer a sophisticated approach to understanding drug disposition. These models consider drug-protein interactions, enabling them to effectively predict drug concentrations in different organs and tissues. This precision aids in accurate drug dosing, providing a significant advantage over conventional models. A key process within these models is equilibration, which ensures that drug concentrations achieve a steady state within the...
378
Modeling with Differential Equations
334
Population dynamics can be described mathematically by considering the population size P(t) as a function of time. The rate of change of the population is then represented by the derivative of P(t). A simple assumption is that the rate of growth is proportional to the size of the population itself. This leads to an exponential growth model, where the population increases rapidly without bound. While this is a useful first approximation, it does not reflect realistic long-term...
334
Model Approaches for Pharmacokinetic Data: Distributed Parameter Models
333
Pharmacokinetic models are mathematical constructs that represent and predict the time course of drug concentrations in the body, providing meaningful pharmacokinetic parameters. These models are categorized into compartment, physiological, and distributed parameter models.
The distributed parameter models are specifically designed to account for variations and differences in some drug classes. This model is particularly useful for assessing regional concentrations of anticancer or...
The distributed parameter models are specifically designed to account for variations and differences in some drug classes. This model is particularly useful for assessing regional concentrations of anticancer or...
333


