Benchmarking Rapid TLES Simulations of Gas Diffusion in Proteins: Mapping O2 Migration and Escape in Myoglobin as a Case Study.

Standard molecular dynamics (MD) simulations of gas diffusion consume considerable computational time and resources even for small proteins. To combat this, temperature-controlled locally enhanced sampling (TLES) examines multiple diffusion trajectories per simulation by accommodating multiple noninteracting copies of a gas molecule that diffuse independently, while the protein and water molecules experience an average interaction from all copies. Furthermore, gas migration within a protein matrix can be accelerated without altering protein dynamics by increasing the effective temperature of the TLES copies. These features of TLES enable rapid simulations of gas diffusion within a protein matrix at significantly reduced (∼98%) computational cost. However, the results of TLES and standard MD simulations have not been systematically compared, which limits the adoption of the TLES approach. We address this drawback here by benchmarking TLES against standard MD in the simulation of O2 diffusion in myoglobin (Mb) as a case study since this model system has been extensively characterized. We find that 2 ns TLES and 108 ns standard simulations map the same network of diffusion tunnels in Mb and uncover the same docking sites, barriers, and escape portals. We further discuss the influence of simulation time as well as the number of independent simulations on the O2 population density within the diffusion tunnels and on the sampling of Mb's conformational space as revealed by principal component analysis. Overall, our comprehensive benchmarking reveals that TLES is an appropriate and robust tool for the rapid mapping of gas diffusion in proteins when the kinetic data provided by standard MD are not required. Furthermore, TLES provides explicit ligand diffusion pathways, unlike most rapid methods.

Quaternary-Linked Changes in Structure and Dynamics That Modulate O2 Migration within Hemoglobin's Gas Diffusion Tunnels.

Atomistic molecular dynamics simulations of diffusion of O2 from the hemes to the external solvent in the α- and ß-subunits of the human hemoglobin (HbA) tetramer reveal transient gas tunnels that are not seen in crystal structures. We find here that the tunnel topology, which encompasses the reported experimental Xe binding cavities, is identical in HbA's T, R, and R2 quaternary states. However, the O2 population in the cavities and the preferred O2 escape portals vary significantly with quaternary structure. For example, most O2 molecules escape from the T ß-subunit via the cavity at the center of the tetramer, but direct exit from the distal heme pocket dominates in the R2 ß-subunit. To understand what triggers the quaternary-linked redistribution of O2 within its tunnels, we examined how the simulated tertiary structure and dynamics of each subunit differs among T, R, and R2 and report that minor adjustments in α-chain dynamics and ß-heme position modulate O2 distribution and escape in HbA. Coupled to the ß-heme position, residue ßF71 undergoes quaternary-linked conformations that strongly regulate O2 migration between the ß-subunit and HbA's central cavity. Remarkably, the distal histidine (HisE7) remains in a closed conformation near the α- and ß-hemes in all states, but this does not prevent an average of 23, 31, and 46% of O2 escapes from the distal heme pockets of T, R, and R2, respectively, via several distal portals, with the balance of escapes occurring via the interior tunnels. Furthermore, preventing or restricting the access of O2 to selected cavities by mutating HisE7 and other heme pocket residues to tryptophan reveals how O2 migration adjusts to the bulky indole ring and sheds light on the experimental ligand binding kinetics of these variants. Overall, our simulations underscore the high gas porosity of HbA in its T, R, and R2 quaternary states and provide new mechanistic insights into why undergoing transitions among these states likely ensures effective O2 delivery by this tetrameric protein.

O2 and Water Migration Pathways between the Solvent and Heme Pockets of Hemoglobin with Open and Closed Conformations of the Distal HisE7.

Hemoglobin transports O2 by binding the gas at its four hemes. Hydrogen bonding between the distal histidine (HisE7) and heme-bound O2 significantly increases the affinity of human hemoglobin (HbA) for this ligand. HisE7 is also proposed to regulate the release of O2 to the solvent via a transient E7 channel. To reveal the O2 escape routes controlled by HisE7 and to evaluate its role in gating heme access, we compare simulations of O2 diffusion from the distal heme pockets of the T and R states of HbA performed with HisE7 in its open (protonated) and closed (neutral) conformations. Irrespective of HisE7's conformation, we observe the same four or five escape routes leading directly from the α- or ß-distal heme pockets to the solvent. Only 21-53% of O2 escapes occur via these routes, with the remainder escaping through routes that encompass multiple internal cavities in HbA. The conformation of the distal HisE7 controls the escape of O2 from the heme by altering the distal pocket architecture in a pH-dependent manner, not by gating the E7 channel. Removal of the HisE7 side chain in the GlyE7 variant exposes the distal pockets to the solvent, and the percentage of O2 escapes to the solvent directly from the α- or ß-distal pockets of the mutant increases to 70-88%. In contrast to O2, the dominant water route from the bulk solvent is gated by HisE7 because protonation and opening of this residue dramatically increase the rate of influx of water into the empty distal heme pockets. The occupancy of the distal heme site by a water molecule, which functions as an additional nonprotein barrier to binding of the ligand to the heme, is also controlled by HisE7. Overall, analysis of gas and water diffusion routes in the subunits of HbA and its GlyE7 variant sheds light on the contribution of distal HisE7 in controlling polar and nonpolar ligand movement between the solvent and the hemes.

Effective simulations of gas diffusion through kinetically accessible tunnels in multisubunit proteins: O2 pathways and escape routes in T-state deoxyhemoglobin.

The diffusion of small gases to special binding sites within polypeptide matrices pivotally defines the biochemical specificity and reactivity of proteins. We investigate here explicit O(2) diffusion in adult human hemoglobin (HbA) as a case study employing the recently developed temperature-controlled locally enhanced sampling (TLES) method and vary the parameters to greatly increase the simulation efficiency. The method is carefully validated against standard molecular dynamics (MD) simulations and available experimental structural and kinetic data on ligand diffusion in T-state deoxyHbA. The methodology provides a viable alternative approach to traditional MD simulations and/or potential of mean force calculations for: (i) characterizing kinetically accessible diffusion tunnels and escape routes for light ligands in porous proteins; (ii) very large systems when realistic simulations require the inclusion of multiple subunits of a protein; and (iii) proteins that access short-lived conformations relative to the simulation time. In the case of T-state deoxyHbA, we find distinct ligand diffusion tunnels consistent with the experimentally observed disparate Xe cavities in the α- and ß-subunits. We identify two distal barriers including the distal histidine (E7) that control access to the heme. The multiple escape routes uncovered by our simulations call for a review of the current popular hypothesis on ligand escape from hemoglobin. Larger deviations from the crystal structure during simulated diffusion in isolated α- and ß-subunits highlight the dampening effects of subunit interactions and the importance of including all subunits of multisubunit proteins to map realistic kinetically accessible diffusion tunnels and escape routes.

Mechanisms of guanosine triphosphate hydrolysis by Ras and Ras-GAP proteins as rationalized by ab initio QM/MM simulations.

The hydrolysis reaction of guanosine triphosphate (GTP) by p21(ras) (Ras) has been modeled by using the ab initio type quantum mechanical-molecular mechanical simulations. Initial geometry configurations have been prompted by atomic coordinates of the crystal structure (PDBID: 1QRA) corresponding to the prehydrolysis state of Ras in complex with GTP. Multiple searches of minimum energy geometry configurations consistent with the hydrogen bond networks have been performed, resulting in a series of stationary points on the potential energy surface for reaction intermediates and transition states. It is shown that the minimum energy reaction path is consistent with an assumption of a two-step mechanism of GTP hydrolysis. At the first stage, a unified action of the nearest residues of Ras and the nearest water molecules results in a substantial spatial separation of the gamma-phosphate group of GTP from the rest of the molecule (GDP). This phase of hydrolysis process proceeds through the low barrier (16.7 kcal/mol) transition state TS1. At the second stage, the inorganic phosphate is formed in consequence of proton transfers mediated by two water molecules and assisted by the Gln61 residue from Ras. The highest transition state at this segment, TS3, is estimated to have an energy 7.5 kcal/mol above the enzyme-substrate complex. The results of simulations are compared to the previous findings for the GTP hydrolysis in the Ras-GAP (p21(ras)-p120(GAP)) protein complex. Conclusions of the modeling lead to a better understanding of the anticatalytic effect of cancer causing mutation of Gln61 from Ras, which has been debated in recent years.