Mandal_2026_J.Chem.Theory.Comput__

Transport tunnels in enzymes with buried active sites are critical gatekeepers of enzymatic function, controlling substrate access, product release, and catalytic efficiency. Despite their importance, the transient nature of these tunnels makes them difficult to study using conventional simulation methods. In this study, we systematically evaluate three coarse-grained (CG) molecular dynamics approachesMartini with Elastic network restraints, Martini with Go-model restraints, and SIRAHfor their ability to characterize tunnel structure and dynamics across diverse enzyme classes. Using haloalkane dehalogenase LinB and its engineered variants as model systems, we show that CG methods accurately reproduce the geometry of tunnel ensembles observed in all-atom (AA) simulations while providing notable computational speedups. The Martini-Go model performed particularly well, capturing subtle mutation-induced changes in tunnel dynamics, such as the closure of a main tunnel and the de novo opening of a transient auxiliary tunnel in LinB variants. In contrast, Martini with Elastic network restraints was limited in capturing tunnel dynamics due to the structural bias introduced by the restraints. We further validated these findings across nine enzymes from the oxidoreductase, transferase, and hydrolase classes with diverse structural folds. Although all CG methods reliably identified functionally relevant tunnels and provided fairly accurate estimates of their ensemble geometry and key bottleneck residues, they differed in their ability to replicate tunnel dynamics, with tunnel occurrences and ranking showing moderate to good correspondence with AA results. This comprehensive evaluation highlights the strengths and weaknesses of CG simulations, establishing them as powerful tools for high-throughput analysis of enzyme tunnels, which enables more efficient enzyme engineering and drug design efforts targeting these critical structural features.