Mark D. Paulsen

Predicting the product specificity and coupling of cytochrome P450cam

SummaryWe present an analysis of several molecular dynamics trajectories of substrate-bound cytochrome P450cam. Trajectories were calculated for the native substrate, camphor, as well as for the alternative substrates, norcamphor and thiocamphor. The system modeled consisted of the crystallographically resolved amino acids, the heme group with a single oxygen atom as the distal ligand, the bound substrate, and the crystallographic waters. These trajectories of the presumptive ferryl oxygen intermediate were used to predict regiospecificity of hydroxylation and coupling between NADH consumption and product formation. Simple geometric criteria in combination with electronic considerations were used to calculate the probability of hydroxylation at specific sites on the substrate. We found that for all the cases examined, the predicted product ratios were in good agreement with the experimentally observed values. We also determined that these simple geometric criteria can be used to predict the degree of coupling between NADH consumption and product formation for a given substrate, which was in good agreement with the experimental values.

Ethylbenzene hydroxylation by cytochrome P450cam

The metabolism of ethylbenzene by cytochrome P450cam was analyzed by experimental and theoretical methods. The present experiments indicate that ethylbenzene is hydroxylated almost exclusively at the secondary ethyl carbon with about a 2:1 ratio of R:S product. Several molecular dynamics trajectories were performed with different starting conformations of ethylbenzene in the active site of P450cam. The stereochemistry of hydroxylation predicted from the molecular dynamics simulations was found to be in good agreement with the observed products.

Analysis of active site motions from a 175 picosecond molecular dynamics simulation of camphor-bound cytochrome P450cam.

The structure and internal motions of the active site residues of camphor-bound cytochrome P450cam have been evaluated on the basis of a 175 psec molecular dynamics simulation. The active site residues generally show very small deviations away from their starting crystal positions. These residues also generally show much smaller fluctuations than for the enzyme as a whole. Phe 87 is dynamically very unusual and is suggested to play a role in substrate movement into and/or out of the active site. The average distance between the heme iron and atoms C5, C6, and C3 of camphor is 5.3, 6.0, and 7.0 A, respectively. This trend is consistent with the experimentally observed stereospecificity of the hydroxylation reaction. On the basis of distance and angle criteria, both 5-exo and 5-endo hydrogen abstraction are predicted to occur during the hydroxylation reaction; although the 5-exo pathway is expected to be 3-fold more likely.

3 Nsec Molecular Dynamics Simulation of the Protein Ubiquitin and Comparison with X-ray Crystal and Solution NMR Structures

Mainly due to computational limitations, past protein molecular dynamics simulations have rarely been extended to 300 psec; we are not aware of any published results beyond 350 psec. The present work compares a 3000 psec simulation of the protein ubiquitin with the available x-ray crystallographic and solution NMR structures. Aside from experimental structure availability, ubiquitin was studied because of its relatively small size (76 amino acids) and lack of disulfide bridges. An implicit solvent model was used except for explicit treatment of waters of crystallization. We found that the simulated average structure retains most of the character of the starting x-ray crystal structure. In two highly surface accessible regions, the simulation was not in agreement with the x-ray structure. In addition, there are six backbone-backbone hydrogen bonds that are in conflict between the solution NMR and x-ray crystallographic structures; two are bonds that the NMR does not locate, and four are ones that the two methods disagree upon the donor. Concerning these six backbone-backbone hydrogen bonds, the present simulation agrees with the solution NMR structure in five out-of-the six cases, in that if a hydrogen bond is present in the x-ray structure and not in the NMR structure, the bond breaks within 700 psec. Of the two hydrogen bonds that are found in the NMR structure and not in the x-ray structure, one forms at 1400 psec and the other forms rarely. The present results suggest that relatively long molecular dynamics simulations, that use protein x-ray crystal coordinates for the starting structure and a computationally efficient solvent representation, may be used to gain an understanding of conformational and dynamic differences between the solid-crystal and dilute-solution states.

[38] Using molecular modeling and molecular dynamics simulation to predict P450 oxidation products

Publisher Summary This chapter focuses on the use of molecular modeling and molecular dynamics simulation to predict P450 oxidation products. The chapter overviews the application of computational methods to P450s and discusses the use of molecular dynamics simulations for the accurate assessment of oxidation specificity by P450 systems. It is observed that the substrates for the cytochromes P450I and P450II subfamilies cluster into two groups; substrates for the P450I subfamily are more planar in shape and exhibit smaller energy differences between the highest occupied and lowest unoccupied molecular orbitals (HOMO and LUMO, respectively) than the substrates for the P450II subfamily. Studies indicate that if due attention is given to ensure proper sampling of the important substrate conformations, accurate product profiles can be predicted for a number of reactions of P450, including hydroxylations, epoxidations, and sulfoxidations. Studies suggest that under appropriate conditions, many time-saving approximations such as the inclusion of little or no added explicit solvent or simulating the dynamics of only the substrate binding pocket can still lead to reliable results.

Active-site mobility inhibits reductive dehalogenation of 1,1,1-trichloroethane by cytochrome P450cam

SummaryRecent studies by Wackett and co-workers have shown that cytochrome P450cam is capable of reductively dehalogenating hexachloroethane at a significant rate, but that no appreciable dehalogenation of 1,1,1-trichloroethane is observed. A growing body of evidence indicates that differences in intrinsic reactivity can not completely explain this observation. We therefore explored the possible role of differences in preferred binding orientation and in active-site mobility. A detailed analysis of molecular dynamics trajectories with each of these substrates bound at the active site of P450cam is presented. While the dynamics and overall time-average structure calculated for the protein are similar in the two trajectories, the two substrates behave quite differently. The smaller substrate, 1,1,1-trichloroethane, is significantly more mobile than hexachloroethane and has a preferred orientation in which the substituted carbon is generally far from the heme iron. In contrast, for hexachloroethane, one of the chlorine atoms is nearly always in van der Waals contact with the heme iron, which should favor the initial electron transfer step.

Calculation of atom-centered partial charges for heme

Atom-centered partial charges are calculated for the Fe-heme in cytochrome P450cam for use in molecular dynamics simulations of polar substrates bound in the active site of the enzyme. Charges are fit to the electrostatic potential produced by ab initio UHF wavefunctions for an Fe-porphine model. Basis set dependence of these charges is observed using the LANL1DZ, LANL2DZ and augmented 6–31G levels of theory. Upon geometry optimization of the enzyme, these charge sets cause varying degrees of distortion of the porphyrin from its crystallographically observed conformation. Scaling the charges calculated from the augmented 6–31G basis by 75% reduces the heme distortion while preserving reasonable interactions with a polar substrate. A comparison of the calculated charges with other published values is presented.