Abstract:Tunnels are structural conduits in biomolecules responsible for transporting chemical compounds and solvent molecules from the active site. They have been shown to be present in a wide variety of enzymes across all functional and structural classes. However, the study of such pathways is experimentally challenging, because they are typically transient. Computational methods, such as molecular dynamics (MD) simulations, have been successfully proposed to explore tunnels. Conventional MD (cMD) provides structural details to characterize tunnels but suffers from sampling limitations to capture rare tunnel openings on longer time scales. Therefore, in this study, we explored the potential of Gaussian accelerated MD (GaMD) simulations to improve the exploration of complex tunnel networks in enzymes. We used the haloalkane dehalogenase LinB and its two variants with engineered transport pathways, which are not only well-known for their application potential but have also been extensively studied experimentally and computationally regarding their tunnel networks and their importance in multistep catalytic reactions. Our study demonstrates that GaMD efficiently improves tunnel sampling and allows the identification of all known tunnels for LinB and its two mutants. Furthermore, the improved sampling provided insight into a previously unknown transient side tunnel (ST). The extensive conformational landscape explored by GaMD simulations allowed us to investigate in detail the mechanism of ST opening. We determined variant-specific dynamic properties of ST opening, which were previously inaccessible due to limited sampling of cMD. Our comprehensive analysis supports multiple indicators of the functional relevance of the ST, emphasizing its potential significance beyond structural considerations. In conclusion, our research proves that the GaMD method can overcome the sampling limitations of cMD for the effective study of tunnels in enzymes, providing further means for identifying rare tunnels in enzymes with the potential for drug development, precision medicine, and rational protein engineering.

Memetic Algorithms for Ligand Expulsion from Protein Cavities

Ligand diffusion in proteins via enhanced sampling in molecular dynamics

Incorporating Prior Knowledge in the Seeds of Adaptive Sampling Molecular Dynamics Simulations of Ligand Transport in Enzymes with Buried Active Sites

Molecular Recognition of Metabotropic Glutamate Receptor Type 1 (mglur1): Synergistic Understanding with Free Energy Perturbation and Linear Response Modeling

Mapping transiently formed and sparsely populated conformations on a complex energy landscape

A Workflow for Exploring Ligand Dissociation from a Macromolecule: Efficient Random Acceleration Molecular Dynamics Simulation and Interaction Fingerprints Analysis of Ligand Trajectories

Identifying ligand binding sites and poses using GPU-accelerated Hamiltonian replica exchange molecular dynamics

Reinforcing Tunnel Network Exploration in Proteins Using Gaussian Accelerated Molecular Dynamics

Guided macro-mutation in a graded energy based genetic algorithm for protein structure prediction

Calculating Protein-Ligand Residence Times Through State Predictive Information Bottleneck based Enhanced Sampling

Computational Study of Protein-Ligand Unbinding for Enzyme Engineering

Simulation of Ligand Transport in Receptors Using CaverDock

The free energy landscape of GABA binding to a pentameric ligand-gated ion channel and its disruption by mutations

Data Driven Classification of Ligand Unbinding Pathways

Rationalizing diverse binding mechanisms to the same protein fold: insights for ligand recognition and biosensor design

Coevolution-Driven Method for Efficiently Simulating Conformational Changes in Proteins Reveals Molecular Details of Ligand Effects in the β2AR Receptor

Enhanced molecular dynamics sampling of drug target conformations

Understanding Protein Dynamics with L1-Regularized Reversible Hidden Markov Models

LARMD: Integration of Bioinformatic Resources to Profile Ligand-Driven Protein Dynamics with a Case on the Activation of Estrogen Receptor.

Sampling of conformational ensemble for virtual screening using molecular dynamics simulations and normal mode analysis