How do substrates enter and products exit the buried active site of cytochrome P450cam? 1. Random expulsion molecular dynamics investigation of ligand access channels and mechanisms.

Cytochrome P450s form a ubiquitous protein family with functions including the synthesis and degradation of many physiologically important compounds and the degradation of xenobiotics. Cytochrome P450cam from Pseudomonas putida has provided a paradigm for the structural understanding of cytochrome P450s. However, the mechanism by which camphor, the natural substrate of cytochrome P450cam, accesses the buried active site is a long-standing puzzle. While there is recent crystallographic and simulation evidence for opening of a substrate-access channel in cytochrome P450BM-3, for cytochrome P450cam, no such conformational changes have been observed either in different crystal structures or by standard molecular dynamics simulations. Here, a novel simulation method, random expulsion molecular dynamics, is presented, in which substrate-exit channels from the buried active site are found by imposing an artificial randomly oriented force on the substrate, in addition to the standard molecular dynamics force field. The random expulsion molecular dynamics method was tested in simulations of the substrate-bound structure of cytochrome P450BM-3, and then applied to complexes of cytochrome P450cam with different substrates and with product. Three pathways were identified, one of which corresponds to a channel proposed earlier on the basis of crystallographic and site-directed mutagenesis data. Exit via the water-filled channel, which was previously suggested to be a product exit channel, was not observed. The pathways obtained by the random expulsion molecular dynamics method match well with thermal motion pathways obtained by an analysis of crystallographic B-factors. In contrast to large backbone motions (up to 4 A) observed in cytochrome P450BM-3 for the exit of palmitoleic acid, passage of camphor through cytochrome P450cam only requires small backbone motions (less than 2.4 A) in conjunction with side-chain rotations. Concomitantly, in almost all the exit trajectories, salt-links that have been proposed to act as ionic tethers between secondary structure elements of the protein, are perturbed.

How do substrates enter and products exit the buried active site of cytochrome P450cam? 2. Steered molecular dynamics and adiabatic mapping of substrate pathways.

Three possible channels by which substrates and products can exit from the buried active site of cytochrome P450cam have been identified by means of random expulsion molecular dynamics simulations. In the investigation described here, we computed estimates of the relative probabilities of ligand passage through the three channels using steered molecular dynamics and adiabatic mapping. For comparison, the same techniques are also applied to investigate substrate egress from cytochrome P450-BM3. The channel in cytochrome P450cam, for which there is the most supporting evidence from experiments (which we name pathway 2a), is computed to be the most probable ligand exit channel. It has the smallest computed unbinding work and force. For this channel, the ligand exits between the F/G loop and the B' helix. Two mechanistically distinct, but energetically similar routes through this channel were observed, showing that multiple pathways along one channel are possible. The probability of ligand exit via the next most probable channel (pathway 3), which is located between the I helix and the F and G helices, is estimated to be less than 1/10 of the probability of exit along pathway 2a. Low-frequency modes of the protein extracted from an essential dynamics analysis of a 1 ns duration molecular dynamics simulation of cytochrome P450cam with camphor bound, support the opening of pathway 2a on a longer timescale. On longer timescales, it is therefore expected that this pathway becomes more dominant than estimated from the present computations.

Towards molecular dynamics simulation of large proteins with a hydration shell at constant pressure.

Molecular dynamics simulation of a large protein in explicit water with periodic boundary conditions is extremely demanding in terms of computation time. Consequently, we have sought approximations of the solvent environment that model its important features. Here, we describe our SAPHYR (Shell Approximation for Protein HYdRation) model in which the protein is surrounded by a shell of water molecules maintained at constant pressure. In addition to the usual pairwise interatomic interactions, these water molecules are subjected to forces approximating van der Waals and dipole-dipole interactions with the implicit surrounding bulk solvent. The SAPHYR model is tested for a system of one argon atom in water and for the protein ubiquitin, and then applied to cytochrome P450cam, a protein with over 400 residues. The results demonstrate that structural and dynamic properties of the simulated systems are improved by use of the SAPHYR model, and that this model provides a significant computational saving over simulations with periodic boundary conditions.

Electrostatic steering and ionic tethering in enzyme-ligand binding: insights from simulations.

To bind at an enzyme's active site, a ligand must diffuse or be transported to the enzyme's surface, and, if the binding site is buried, the ligand must diffuse through the protein to reach it. Although the driving force for ligand binding is often ascribed to the hydrophobic effect, electrostatic interactions also influence the binding process of both charged and nonpolar ligands. First, electrostatic steering of charged substrates into enzyme active sites is discussed. This is of particular relevance for diffusion-influenced enzymes. By comparing the results of Brownian dynamics simulations and electrostatic potential similarity analysis for triose-phosphate isomerases, superoxide dismutases, and beta-lactamases from different species, we identify the conserved features responsible for the electrostatic substrate-steering fields. The conserved potentials are localized at the active sites and are the primary determinants of the bimolecular association rates. Then we focus on a more subtle effect, which we will refer to as "ionic tethering." We explore, by means of molecular and Brownian dynamics simulations and electrostatic continuum calculations, how salt links can act as tethers between structural elements of an enzyme that undergo conformational change upon substrate binding, and thereby regulate or modulate substrate binding. This is illustrated for the lipase and cytochrome P450 enzymes. Ionic tethering can provide a control mechanism for substrate binding that is sensitive to the electrostatic properties of the enzyme's surroundings even when the substrate is nonpolar.