Combined Enhanced Conformational Sampling and Protein-Ligand Dockings on CYP3A4 Improves in Silico Prediction of Drug Metabolism
4 affiliations
- 1Insitute of Biotechnology and Biomedicine & Departament de Química
- 2Universitat Autònoma de Barcelona
- 3Jealott's Hill International Research Centre
- 4Syngenta Crop Protection, Schaffhauserstrasse 101, Stein, Switzerland, 4332
Abstract
Cytochrome P450 3A4 (CYP3A4) is the primary drug-metabolizing enzyme in humans. This is largely due to its highly flexible structure, which enables it to bind a broad spectrum of xenobiotics with a wide variety of shapes and physicochemical properties. This plasticity, however, presents challenges for in silico substrate prediction methods, in particular those that consider only a limited conformational space of the receptor or that rely on static crystal structures. In this study, we aim to explore the conformational landscape of CYP3A4 using Gaussian accelerated Molecular Dynamics (GaMD) and then to assess the viability of using this method to identify pro-reactive conformations through protein-ligand docking. GaMD simulations were initiated from distinct X-ray structures in triplicate 500 ns runs. The results demonstrate that the simulations cover a very large conformational space, including all reported X-ray structures as well as conformations with significantly expanded active-site geometries unseen in the experimental data. Protein-ligand docking of known CYP3A4 substrates onto the main cluster representatives from GaMD simulations allows the identification of metabolically active poses of the substrates. This is true even in cases where substrate-bound X-ray structures or docking against a single structure fails to yield productive results, as in the case of CYP3A4–Erythromycin. Furthermore, a classifier trained on ensemble docking scores discriminates CYP3A4 substrates from non-substrates better than any single conformation (ROC-AUC 0.735 vs 0.697), with performance increasing with ensemble size. Overall, the study shows that enhanced conformational sampling is crucial for accurately modelling highly flexible enzymes like CYP3A4 and can significantly improve metabolic predictions.
Related
- PLUMED-JAX: GPU-native enhanced sampling for atomistic simulations
- Enhanced surface hydrophobicity in the combined treatment of chaplins and rodlins from Streptomyces griseus
- METALICA: METAdynamics and repLICA exchange for enhanced diffusion sampling
- Reactive calcium carbonate precipitation from an atomic cluster expansion potential and enhanced sampling
- Enhanced sampling enables binding site water reorganization in polarizable absolute binding free energy calculations
BibTeX
@misc{boulaamane2026combined,
title = {{Combined Enhanced Conformational Sampling and Protein-Ligand Dockings on CYP3A4 Improves in Silico Prediction of Drug Metabolism}},
author = {Yassir Boulaamane and Christopher M. Baker and Lea Talmann and Jean-Didier Marechal},
howpublished = {ChemRxiv},
year = {2026},
month = oct,
doi = {10.26434/chemrxiv.15010077/v1},
url = {https://chemrxiv.org/doi/full/10.26434/chemrxiv.15010077/v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Combined Enhanced Conformational Sampling and Protein-Ligand Dockings on CYP3A4 Improves in Silico Prediction of Drug Metabolism
AU - Yassir Boulaamane
AU - Christopher M. Baker
AU - Lea Talmann
AU - Jean-Didier Marechal
PY - 2026
DA - 2026/10/08
PB - ChemRxiv
DO - 10.26434/chemrxiv.15010077/v1
UR - https://chemrxiv.org/doi/full/10.26434/chemrxiv.15010077/v1
AB - Cytochrome P450 3A4 (CYP3A4) is the primary drug-metabolizing enzyme in humans. This is largely due to its highly flexible structure, which enables it to bind a broad spectrum of xenobiotics with a wide variety of shapes and physicochemical properties. This plasticity, however, presents challenges for in silico substrate prediction methods, in particular those that consider only a limited conformational space of the receptor or that rely on static crystal structures. In this study, we aim to explore the conformational landscape of CYP3A4 using Gaussian accelerated Molecular Dynamics (GaMD) and then to assess the viability of using this method to identify pro-reactive conformations through protein-ligand docking. GaMD simulations were initiated from distinct X-ray structures in triplicate 500 ns runs. The results demonstrate that the simulations cover a very large conformational space, including all reported X-ray structures as well as conformations with significantly expanded active-site geometries unseen in the experimental data. Protein-ligand docking of known CYP3A4 substrates onto the main cluster representatives from GaMD simulations allows the identification of metabolically active poses of the substrates. This is true even in cases where substrate-bound X-ray structures or docking against a single structure fails to yield productive results, as in the case of CYP3A4–Erythromycin. Furthermore, a classifier trained on ensemble docking scores discriminates CYP3A4 substrates from non-substrates better than any single conformation (ROC-AUC 0.735 vs 0.697), with performance increasing with ensemble size. Overall, the study shows that enhanced conformational sampling is crucial for accurately modelling highly flexible enzymes like CYP3A4 and can significantly improve metabolic predictions.
ER -