Quantitative Structural Insight into Human Variegate Porphyria Disease*

Background: Defects in the human protoporphyrinogen oxidase (PPO) causes the variegate porphyria (VP) disease. Results: The activity of 44 clinically reported VP-causing mutants have been accurately predicted. Conclusion: A quantitative understanding into the molecular basis of VP by PPO mutation was established. Significance: The quantitative insight into VP should be helpful for the diagnosis, precaution, and treatment of this disease. Defects in the human protoporphyrinogen oxidase (hPPO) gene, resulting in ∼50% decreased activity of hPPO, is responsible for the dominantly inherited disorder variegate porphyria (VP). To understand the molecular mechanism of VP, we employed the site-directed mutagenesis, biochemical assays, structural biology, and molecular dynamics simulation studies to investigate VP-causing hPPO mutants. We report here the crystal structures of R59Q and R59G mutants in complex with acifluorfen at a resolution of 2.6 and 2.8 Å. The r.m.s.d. of the Cα atoms of the active site structure of R59G and R59Q with respect to the wild-type was 0.20 and 0.15 Å, respectively. However, these highly similar static crystal structures of mutants with the wild-type could not quantitatively explain the observed large differences in their enzymatic activity. To understand how the hPPO mutations affect their catalytic activities, we combined molecular dynamics simulation and statistical analysis to quantitatively understand the molecular mechanism of VP-causing mutants. We have found that the probability of the privileged conformations of hPPO can be correlated very well with the kcat/Km of PPO (correlation coefficient, R2 > 0.9), and the catalytic activity of 44 clinically reported VP-causing mutants can be accurately predicted. These results indicated that the VP-causing mutation affect the catalytic activity of hPPO by affecting the ability of hPPO to sample the privileged conformations. The current work, together with our previous crystal structure study on the wild-type hPPO, provided the quantitative structural insight into human variegate porphyria disease.

Inhibition of PPO in plants leads the accumulation of protogen, which exports to the cytoplasm and forms porphyrin through nonenzymatic oxidation. The photosensitizing porphyrin can lead to singlet oxygen, which can cause peroxidation of membrane lipids and cell death (18,19). In human, because the deficiency of PPO leads to the tissue accumulation and excessive excretion of porphyrin, it is hypothesized that the sensitivity of VP patients to light should be similar to the condition in plants (20,21).
Though the study of VP has been performed for more than half a century (8,17,21,22), the detailed molecular mechanism of VP is still unclear. To address this important issue, we have resolved the crystal structure of human PPO in a resolution of 1.9 Å and characterized 44 VP-causing mutations in vitro (22). Based on this high-resolution structure of hPPO and the kinetic data, we attempted to explain the defects of 44 missense mutants. The mapping of these 44 mutations on the hPPO crystal structure indicated these mutants distributed on all of the three structure domains of hPPO. However, our understanding of the molecular mechanism of VP-causing mutants was limited because our insights mostly depend on structural analysis and in silico docking (22); the detailed information of molecular motion on substrate binding was especially lacking. In this work, the site-directed mutagenesis, biochemical assays, structural biology, and molecular dynamics simulation studies have been carried out to understand molecular mechanism of VPcausing mutants. We report here the crystal structures of R59Q and R59G mutants in complex with acifluorfen at a resolution of 2.6 and 2.8 Å. However, these highly similar static crystal structures of mutants with the wild-type could not quantitatively explain the observed large differences in their enzymatic activity. We further combined molecular dynamics simulation and statistical analysis to quantitatively understand the molecular mechanism of VP-causing mutants. We have found that the VP-causing mutations affect the catalytic activity of hPPO by affecting the ability to sample privileged conformations. The effects of 44 clinically important mutations on the catalytic activity of hPPO have been accurately predicted with R 2 ϭ 0.94 in this work.

EXPERIMENTAL PROCEDURES
Plasmid Construction of hPPO Mutants-The recombinant plasmid (pHPPO-X) was a generous gift from Dr. H. A. Dailey (University of Georgia). All mutations of hPPO described here were generated from the recombinant plasmid (pHPPO-X) using DpnI-mediated site-directed mutagenesis methods and confirmed by DNA sequencing.
Expression and Purification of Mutant hPPO-Plasmids of wild-type and various hPPO mutants were initially transformed into BL21(DE3)pLysS Escherichia coli cells. The single colony of resulting cells was grown in 10 ml of 2ϫYT medium with 100 g/ml ampicillin at 37°C overnight. Total culture were then inoculated into one liter of 2ϫYT medium containing 100 g/ml ampicillin and grown with shaking at 37°C to A 600 of 0.6. The expression of the recombinant hPPO wild-type enzyme and its mutants were induced by adding isopropyl 1-thio-␤-Dgalactopyranoside to a final concentration of 1 mM. Cells were grown for an additional 4 h at 25°C.
Enzyme Assay for Kinetic Analysis-PPO activity was assayed by measuring the constant velocity of formation of protopor-phyrin IX from protoporphyrinogen IX using a 96-well plate with the continuous fluorometric methods at 25°C (23)(24)(25). The product has a maximum excitation wavelength at 410 nm and a maximum emission wavelength at 630 nm. The total volume of the reaction mixture was 200 l consisting of 0.1 M potassium phosphate buffer (pH 7.4), 5 M FAD, 5 mM DTT, 1 mM EDTA, 0.2 M imidazole, and 0.03% Tween 80 (v/v). The enzymatic reaction was started by the addition of substrate. The autoxidation rates were determined concomitantly and were subsequently subtracted. Kinetic parameters, including the Michaelis-Menten constant (K m ), the maximal velocity (V max ), and the catalytic constant (k cat ), were determined by a Lineweaver-Burk plot.
Crystallization and Data Collection-Crystals of the mutant hPPOs⅐acifluorfen complexes were grown using the same method as the wild-type protein described previously (22). Diffraction data were collected on beam station BL17U1 of the Shanghai Synchrotron Radiation Facility. The data of R59Q and R59G were processed and scaled to 2.6 and 2.8 Å, respectively, using the HKL2000 software package (26). The crystals of the complex of mutant hPPO with acifluorfen belong to the space group R3.
Structure Determination and Refinement-The initial phases were obtained by molecular replacement using the structure of PPO enzyme from Homo sapiens (Protein Data Bank code 3NKS (22)) as a template. The program PHASER was able to locate one molecule in the asymmetric unit (27). The models were built manually using the program COOT (28) and refined using the program CNS (29) and Phenix. The final structure of R59Q mutant had an R cryst value of 17.9% and an R free value of 24.3%. The Ramachandran plot calculated by the program PROCHECK (30) showed that 90.3% of the residues were in their most favored regions, 9.7% of the residues were in additionally allowed regions, and no residues were in generously allowed regions and disallowed regions. The structure of R59G mutant had an R cryst value of 17.5% and an R free value of 24.4%. The Ramachandran plot calculated by the program PRO-CHECK (30) showed that 89.2% of the residues were in their most favored regions, 10.5% of the residues were in additionally allowed regions, 0.3% of the residues were in generously allowed regions, and no residues disallowed regions. Detailed data collection and refinement statistics are summarized in supplemental Table S1.
Molecular Dynamics Simulations-MD simulations for wildtype and mutant hPPO in complexes with substrate protogen were carried out to generate the equilibrium conformations for each of the studied systems. The crystal structures of the wildtype hPPO (Protein Data Bank code 3NKS (22)) and R59G and

Quantitative Insight into Variegate Porphyria
R59Q mutants were taken as the starting point for MD simulations of the complexes. All missing residues of the structure were added using InsightII (Accelrys, Inc., San Diego, CA). The hydrogen atoms were added to hPPO structure using leap module of AMBER (version 9) (31). The structures of other hPPO mutants were generated from the structure by mutation of corresponding residues using InsightII (Accelrys, Inc.). A 3000step minimization was performed on the mutated residues and the adjacent residues within 4.0 Å from the mutated residue to remove the steric clash. The optimization were performed with SANDER module of AMBER (version 9) (31) and AMBER (version 03) force fields (32,33).
The coordinates of protogen were taken from the previous docking studies where protogen was docked into the crystal structure of hPPO (22). The force field parameters of cofactor FAD and protogen were generated with the Gaussian 03 (34) and AMBER (version 9).
Each of the complex systems was immersed in a truncated octahedral box of explicit TIP3P water molecules and was neutralized by adding counter ions (Cl Ϫ or Na ϩ ). All MD simulations were carried out with the SANDER module of AMBER (version 9) (31). Prior to MD simulations, the molecular systems were subjected to a series of energy minimizations to relax the system. In MD simulations, particle mesh Ewald (35) and periodic boundary condition were employed to deal with the long range electrostatic interactions. The SHAKE method (36) was used to constrain bonds involving hydrogen atoms to tolerate 2-fs time steps. The system was gradually heated from 0 to 300 K over 50 ps. Equilibrating calculation was executed at 1 atm and at 300 K for 50 ps. Then, 5-ns MD simulations were performed under 300 K and 1 atm. The snapshot of the system was taken every 1 ps. Trajectories of reactants protogen and FAD after the initial 2 ns of equilibration were considered for analysis; thus, 3000 MD structures, collected during the time period of 2001-5000 ps, were generated for each of the systems using the ptraj module of AMBER (version 9). To test the consistency of data from MD simulations, the simulation of wildtype hPPO system was extended to 44 ns.
Description of the Conformational Space of the Reactants-Based on the proposed PPO catalytic mechanism (4, 37) and our previous study (22), the distance between C5 atom of protogen and N5 atom of FAD (D C5-N5 ) and the angle between C5, N5 and N10 (A C5-N5-N10 ) (shown in Fig. 2) as geometry param-eters were used to describe the reactant conformational space. Each conformation of the reactants was represented by a pair of D C5-N5 -A C5-N5-N10 . For each of the studied systems, conformational ensemble of the reactants was represented by a set of D C5-N5 -A C5-N5-N10 pairs, corresponding to a set of discrete data points in a D C5-N5 -A C5-N5-N10 plane.
Construction of Conformation Probability Density Function (CPDF)-To calculate the probability of the reactant conformation in any given interval of D C5-N5 and A C5-N5-N10 , for each of the studied systems, we converted the sets of discrete data points of D C5-N5 -A C5-N5-N10 pairs to the CPDF, denoted as Z(A,D). The conversion was performed in the following steps: (i) dividing the D C5-N5 -A C5-N5-N10 plane into square bins with 0.1 Å width and 5°length; (ii) taking count of the number of data points for each bin; (iii) calculating the density of each bin by the number of D C5-N5 -A C5-N5-N10 data points fall in the corresponding bin being divided by the total number of the angledistance data points; (iv) fitting the data with the two-dimensional Gaussian function to construct Z(A,D) by the non-linear surface fitting methods implemented in Origin Lab (version 8.0, Northampton, MA). The general form of the CPDF can be presented as follows in Equation 1, where Z denotes the conformation probability density in D C5-N5 and A C5-N5-N10 for the observing systems. The factor C in this expression ensures that the total area under the curve Z(A,D) is equal to one. The parameter is called the mean of the normal distribution, and it determines the coordinate of highest point of the Gaussian surface. The parameter is called the S.D. and describes how concentrated the distribution is around its mean.
Calculation of the Probability of PC-The probability of any conformation of the reactants at equilibrium was calculated by taking the integral of the CPDF over the given integration intervals. Based on the x-ray structures of the flavoenzyme complexes with substrate analogues (38), we defined the conformations with D C5-N5 within 3.0 -3.8 Å and A C5-N5-N10 within 96 -117°as the privileged conformations (PC). For each studied system, the probability of PC, denoted P PC , was calculated by the bivariate integration of corresponding CPDF over the interval (96, 117) for A C5-N5-N10 , and (3.0, 3.8) for D C5-N5 (Equation 2) using Matlab software (The Math Works, Natick, MA). The general integral form was expressed as follows, Correlation of the P PC with the k cat /K m of PPO Mutants-Nonlinear regression analyses were performed in Origin Lab (version 8.0) to explore whether the P PC for the studied systems could be correlated with corresponding k cat /K m values. Best fitting was assessed by correlation coefficient, square root of mean square error, and F-test comparisons. Accession Codes-Coordinates and associated structure factors of solved structures for human protoporphyrinogen IX oxidase R59Q and R59G mutants in complex with acifluorfen have been deposited in the Protein Data Bank (codes 4IVO and 4IVM, respectively).

RESULTS AND DISCUSSION
Crystal Structure of R59Q and R59G Mutant-The R59W mutation is by far the most prevalent VP-causing mutation identified, and Ͼ94% of the patients in South Africa carry this mutation (16,17). Biochemical studied showed that the R59W mutation dramatically decreased the catalytic constant of hPPO in vitro (17,21,41,42). Herein, we have determined the crystal structures of R59Q and R59G mutants in complex with acifluorfen at a resolution of 2.6 and 2.8 Å, respectively, with 1 molecule/asymmetric unit. Data collection and refinement statistics appear in supplemental Table S1.
Overall, the structures of R59Q, R59G, and wild-type hPPO were highly similar, with r.m.s.d. in C␣ positions of 0.26 and 0.29 Å, respectively (Fig. 3, A and B). The mutation caused minor structural change in the active site residues. The r.m.s.d. of the active site residues (residues within 4.0 Å distance with acifluorfen) of R59Q and R59G mutants with wild-type was 0.20 and 0.15 Å, respectively. In addition, after being subjected to the molecular dynamics simulation with bound substrate, the r.m.s.d. of the active site residues of R59Q and R59G mutants with wild-type was 0.27 and 0.39 Å, respectively. These data suggested that the static active site architecture was not significantly affected by the mutations (Fig. 3, A and B). These highly similar static crystal structures of mutants with the wildtype could not quantitatively explain the observed large differences in their enzymatic activity (Fig. 3C). To answer the question how the hPPO mutations affect their catalytic activities, the dynamics motion of the structure of hPPO mutants should be considered.
Conformational Distribution of the Reactants-Dynamics simulation methods are widely used to obtain information on the time evolution of conformations of proteins and other biological macromolecules, which present a powerful complement to static structures of proteins and other biomolecules. Convergence of results from simulation with experimental data provides a foundation for using simulations to enhance experimental observations and study features of biomolecules that cannot be seen directly in experiments. Because it was difficult to study the effect of the mutants only from the comparison of the crystal structures of wild-type and mutant hPPO, we have therefore chosen to study the molecular mechanism of how VP-causing mutants affect the activity of hPPO through molecular dynamics simulations of the hPPO⅐substrate complex. First, we chose wild-type hPPO and 16 hPPO mutants (including the R59G and R59Q mutants) as the studied systems (Fig. 4A), because all of these mutations are involved in the binding of substrate or FAD and clinically found in the VP patients (7, 9 -11). MD simulations for the 17 enzymes in complexes with substrate protogen were carried using AMBER (version 9) (31). For each of the studied system, 3000 conformers extracted from the MD trajectory during the time period of 2001-5000 ps were used as the reduced representation of the conformational ensemble at equilibrium state. For each studied system, the fluctuations of temperature, density, and total energy during the MD simulation were examined, and the r.m.s.d. of the backbone heavy atoms with respect to the starting structures were reported as a function of time. These physical quantities from the wild-type hPPO system were shown in supplemental Fig. S1. Other detailed data about the MD trajectories of studied systems were provided in supplemental Figs. S2 and S3 and supplemental Table S2).
It was essential to describe conformational space of the reactants using a set of proper parameters for capturing the salient features of conformational distribution of the reactants at equilibrium system. The oxidation of protogen to protoporphyrin IX catalyzed by PPO involves the rupture of C-H bond of methylene of macrocycle in protogen, coupled to the transfer of two electrons to FAD (4). The carbon atoms involved in these C-H bonds were considered as the sites of oxidative attack, and N5 of the tricyclic isoalloxazine ring system of FAD was believed to take part directly in substrate dehydrogenation (37,(43)(44)(45). In the previous studies, we also explored the binding mode of protogen at the active site of hPPO based on the high-resolution crystal structure of hPPO by the molecular docking, MD simulation, and quantum mechanical calculations. Our studies indicated that the methylene bridge C5 atom of protogen (see Fig. 2) was close to the N5 atom of FAD (22). Thus, the distance between C5 and N5 atoms (D C5-N5 ) and the angle between C5, N5, and N10 (A C5-N5-N10 ) (shown in Fig. 2) as a set of geometry parameters were used to describe the conformational space in the reaction process. In conformational ensembles, each conformation of the reactants (protogen and FAD) was represented by a pair of D C5-N5 -A C5-N5-N10 , whereas each conforma- . Conformation probability analysis of the PPO mutants. A, A total of 17 mutants listed at left column were used to construct the prediction model, and four mutants shown at the right column were the mutants constructed and measured later as the initial test set. The mutated residues, the cofactor FAD, and the substrate protogen were shown in stick representation with gray, cyan, and green carbon atoms, respectively. B, D C5-N5 versus A C5-N5-N10 scatter plots. The distribution of the pairs of D C5-N5 -A C5-N5-N10 for wild-type hPPO from MD simulation was shown in blue, and the distributions for the mutants were shown in red. C, probability distribution of D C5-N5 observed for wild-type hPPO. Horizontal coordinates indicate D C5-N5 value; vertical coordinates indicate probability density. D, probability distribution of A C5-N5-N10 observed for wild-type hPPO. Horizontal coordinates indicate A C5-N5-N10 value; vertical coordinates indicate probability density. E, conformation probability density function surface for wild-type hPPO.
tional ensemble of the reactants was represented by a set of D C5-N5 -A C5-N5-N10 pairs, corresponding to a set of discrete data points in a D C5-N5 -A C5-N5-N10 plane.
In Fig. 4B, the D C5-N5 -A C5-N5-N10 pairs as data points were plotted in a D C5-N5 -A C5-N5-N10 plane, in which each of the scatter plots for the mutants was superposed on the one for wildtype system to compare the conformational spaces of the reactants of wild-type hPPO and mutants. The D C5-N5 versus A C5-N5-N10 scatter plots for the mutants were provided in supplemental Fig. S4. The ranges of D C5-N5 and A C5-N5-N10 and the average of the two values for each conformational ensemble were listed in supplemental Table S3.
From the pairs of D C5-N5 -A C5-N5-N10 scatter plots, we obtained the qualitative understanding of the conformational distribution for each of the studied system. As seen in Fig. 4B and supplemental Table S3, despite the fact that the wild-type, V347A, and F331A system had similar regions where the data points appeared, the data pattern was different for other studied systems. The studied systems with higher catalytic efficiency, for example, wild-type hPPO and F331A, showed more condensed data distribution than that of the systems with lower catalytic efficiency. The results suggested that the catalytic efficiency may be related to the conformational distribution of reactants for our studied systems.
As an estimate of the probability distribution of a continuous variable, the histograms of probability of D C5-N5 and A C5-N5-N10 were investigated for each of the studied systems. Of these, the histograms for the wild-type hPPO, as an example, were displayed in Fig. 4, C and D, in which a bell-shaped probability distributions were observed for both of D C5-N5 and A C5-N5-N10 . The same behaviors were also observed in all of the mutants. For all of the studied systems, the probability distributions of D C5-N5 and A C5-N5-N10 were well approximated by Gaussian distributions, with R 2 of 0.95-0.99 for D C5-N5 , and 0.99 for A C5-N5-N10 , respectively (supplemental Table S4 and supplemental Fig. S5).
Considering that both of distribution of D C5-N5 and A C5-N5-N10 for each studied system followed Gaussian distribution, we fitted the binned data of D C5-N5 -A C5-N5-N10 to a two-dimensional Gaussian distribution using a nonlinear regression implemented in Origin Lab (version 8.0). For all of the 17 studied systems, they were well approximated by Gaussian distributions (Equation 1), and the best estimates of the mean and standard deviation obtained from the fitting procedures for each of studied systems were reported in Table 1. The CPDF surface for the wild-type and mutant hPPO were displayed in Fig. 4E and supplemental Fig. S6.
The Probability of the PC-Based on the x-ray crystal structures of the flavoenzyme complexes with substrate analogues, Fraaije and Mattevi (38) have studied the catalytic features of the flavoenzymes that perform a dehydrogenation reaction. The studies indicated that the geometry parameters between substrate and FAD corresponding to the distance D C5-N5 and the angle A C5-N5-N10 of hPPO were generally conserved for the studied flavoenzymes, within 3.0 -3.8 Å for D C5-N5 , 96 to 117°f or A C5-N5-N10 . Therefore, we defined the conformations with both of D C5-N5 within 3.0 -3.8 Å and A C5-N5-N10 within 96 -117°as the PC. In fact, there is widely accepted proposition that the reactant particles must collide with appropriate orientations to achieve an enzymatic reaction, and enzymes bind to their substrates for the formation of a substrate⅐enzyme complex prior to the reaction step. Here, we quantify the appropriate orientation collision probability between the reactants as the P PC , which was obtained by the bivariate integration of the CPDF over the reactable geometry ranges. Such a quantitation may correlate with the catalytic activity of enzymes.
For each of 17 studied systems, the probability of reactant conformation in any given interval of D C5-N5 and A C5-N5-N10 can be calculated through the corresponding CPDF. Thus, P PC was calculated by the bivariate integration of the corresponding CPDF over the interval (96, 117) for A C5-N5-N10 , and (3.0, 3.8) for D C5-N5 using Matlab software (supplemental Fig. S7). P PC for each of the studied systems were listed in Table 1. The kinetic parameters of the wild-type hPPO and 16 mutants were also listed in Table 1. It can be found that P PC of each studied system goes the same trend with the k cat /K m of enzyme

Quantitative Insight into Variegate Porphyria
mutants, which is in keeping with our first hypothesis that the probability of the appropriate orientated collision between the reactants should correlate with the catalytic activity of enzymes.
Correlation of the Probability of PC with the Catalytic Efficiency of the Enzymes-The k cat /K m values of 17 studied systems were plotted against the corresponding P PC values. Upon inspection of Fig. 5A, we found that the P PC increased with a nonlinear increasing of the k cat /K m . Thus, the paired data were fitted to BoxLucas1 functions, and the best fit was assessed by square root of mean square error values and F-test comparisons in Origin Lab (version 8.0). Finally, the correlation equation with R 2 ϭ 0.90, square root of mean square error ϭ 0.22, and F ϭ 140.82 was expressed as follows in Equation 3.
The predicted values of the k cat /K m from this regression were listed in Table 1. The experimental k cat /K m values are plotted against the predicted k cat /K m in Fig. 5B with correlation R 2 ϭ 0.91. The relatively high R 2 between the experimental and predicted k cat /K m indicated that a high proportion of the data can be explained by Equation 3. Thus, these mutants affect the catalytic activity of hPPO by affecting the P PC. Four more mutants such as F331T, L334T, R59C, and V347T have thus been constructed, and kinetics have been measured to test our model. Their experimental k cat /K m values again correlated very well with the calculated values ( Fig. 5C and supplemental Table S5), which further indicated that the probability of the privileged conformations in the PPO-catalyzed reaction process made the major contribution. Encouraged by these results, we decided to use this computational approach to predict the catalytic activities of all clinically reported VP-causing mutants.
The Molecular Mechanism of the VP-causing Mutants-The mapping of these 44 VP-causing mutants on the hPPO crystal structure indicated these mutants distributed on all of the three structure domains of hPPO (Fig. 6A). We categorized these mutants into five potential categories, i.e. affecting the substrate binding, FAD binding, hydrophobic core, secondary structure, and surface of hPPO. However, our understanding of the molecular mechanism of VP-causing mutants was limited because our insights mostly depend on the static structural analysis and in silico docking. In this work, based on molecular dynamics simulation and statistical analysis, we found the positive correlation of the P PC with k cat /K m of the 17 PPO variants. We set out to further test whether the VP-causing mutants affect the catalytic activity of hPPO by affecting the distribution of PC. Among the 44 mutants, 36 mutants display detectable activity, whereas eight mutants were inactive. (The catalytic activity of the inactive mutants was designated as zero.) The k cat /K m values of these mutations were then predicted using the Equation 3, and the data were listed in supplemental Table S4. Fig. 6B depicts a correlation plot of the predicted and experimental k cat /K m values, in which the correlation coefficient (R 2 ) is 0.94.
Because the quantitative correlation between experimental k cat /K m and the computational P PC could be established, we can therefore assess the molecular mechanism of the VP disease by VP-causing mutants from dynamics simulation. Based on our molecular dynamics simulation and statistical analysis, we found that the VP-causing mutations, even those far from the active site and FAD could significantly impact the P PC of hPPO. These mutations lead to the changes of structural and conformational motions of the entire hPPO that alter the ability of sampling privileged conformations, thereby affecting the catalytic activity of hPPO. Now we can quantitatively understand how the catalytic activities were affected by hPPO mutations, a result cannot come out of the static crystal structure of hPPO mutants.
From dynamics simulation, it could be seen that in the wildtype, the side chain of Pro-58 is stacking with the isoalloxazine ring of FAD, and the average distance of all simulated ternary conformations between the side chain center of Pro-58 and the center of ring B of isoalloxazine ring of FAD was kept at ϳ3.98 Å during the simulation. While in the R59Q mutant, the aver-age distance increased to 5.43 Å (Fig. 6, C and D), which indicated that the stacking interaction between Pro-58 and FAD is disrupted in the reaction process by the mutation. In the dynamics simulation of R59G mutants, the average r.m.s.d. of 1.24 Å for the loop of Leu-56 -Ile-61 of R59G respect to wildtype (higher than the overall structure r.m.s.d., 0.78 Å) indicated that the loop of Leu-56 -Ile-61 undergo a large structural motion compared with the wild-type in the reaction process (Fig. 6, C and E). Because this loop locates above the isoalloxazine ring of FAD, the structural motion of this loop should affect the binding of the isoalloxazine ring of FAD and therefore affect the activity of PPO. The detailed information from molecular dynamics study about the effects of the clinically reported mutations on enzymatic activity were provided in supplemental Table S7. These data indicated that the VP-causing mutations affect the catalytic activity of hPPO by affecting the ability to sample privileged conformations, which is in accord-FIGURE 6. Prediction of the catalytic activity and molecular mechanism of VP-causing mutants. A, the mapping of the VP-causing mutations on the crystal structure of hPPO. Mutants are indicated by colored spheres, substrate and FAD are shown in stick representation. B, plot of the experimental and predicted k cat /K m for the clinically reported or potential VP-causing mutants from the Equation 3. The data points are colored by their effects on the activity of hPPO. The diagonal red line corresponding to (k cat /K m ) exp ϭ (k cat /K m ) pred , the correlation showed R 2 ϭ 0.94. C, the average conformation of the active site from dynamics simulation for the wild-type. D, the average conformation of the active site from dynamics simulation for R59Q. E, the average conformation of the active site from dynamics simulation for R59G.

Quantitative Insight into Variegate Porphyria
ance with the widely accepted proposition that the reactant particles must collide with appropriate orientations to achieve an enzymatic reaction.
Computational Approach for Predicting the Catalytic Efficiency of Enzyme Mutants-Many reports have shown that an enzyme activity depends on the conformational distribution of reactants in enzyme active site (39, 40, 46 -53). In this study, we have shown that the probability of the reactable conformations of reactants could be employed to quantitatively predict the catalytic efficiency of the hPPO mutants. This protocol for predicting catalytic efficiency of mutant hPPO by combining MD simulation/statistical analysis was called Prenzyme (Fig. 7). The protocol was based on the notion that a reaction can occur only through the collisions by the reactants with the privileged conformations, and this favorable collision can be evaluated by the probability of the privileged conformations of ternary catalytic complex. The conformation probability density function can be obtained from statistical analysis of the conformational distribution of the ternary catalytic complex from MD and then gives the probability value by bivariate integration over the intervals of the corresponding geometric parameters of the privileged conformations. The obtained probability should be correlated with k cat /K m of enzyme mutants.