Selectivity principle of the ligand escape process from a two-gate tunnel in myoglobin: molecular dynamics simulation

We proposed a selectivity principle for the ligand escape process from two fluctuating bottlenecks in a cavity with a multigate inside a myoglobin pocket. Our previous analytical theory proposed a fluctuating bottleneck model for a Brownian particle passing through two gates on a cavity surface of an enzyme protein and has determined the escape rate in terms of the time-dependent gate function and the competition effect. It illustrated that with two (or more than two) gates on a cavity surface the gate modulation, which is controlled by protein fluctuation, dominates the ligand escape pathway. We have performed a molecular dynamics simulation to investigate the selectivity principle of the ligand escape process from two-gate tunnel in myoglobin. The simulation results confirm our theoretical conjecture. It indicates that the escape process is actually entropy driven, and the ligand escape pathway is chosen via the gate modulation. This suggests an interesting intrinsic property, that is, the oxymyoglobin tertiary structure is favorable to the departure of the ligand from one direction rather than through a biased random walk.