Investigation by MD simulation of the key residues related to substrate-binding and heme-release in human ferrochelatase

Molecular dynamics (MD) simulations of three models based on the crystal structure of the E343K variant of human ferrochelatase were performed in this study. The “open” and “closed” conformations of the enzyme obtained by simulations are in agreement with the corresponding crystal structures. The snapshots and the structure analysis indicate that alterations of the hydrogen bonds and the positions of E347 and E351 lead to a conformational change in the π-helix. The hydrogen bonded form of residue R164 could be regarded as a signal indicating alteration of the active site conformation. When R164 forms a hydrogen bond with D95, the active site is closed, and when a hydrogen bond is formed with E171, the active site is open. Interestingly, the protoporphyrin with Fe²⁺ is observed to move noticeably out of the enzyme while the protoporphyrin lacking Fe²⁺ remains almost fixed. Alterations of the hydrogen bonds between the propionate of the heme and R115, K118 and S303 trigger movement of the heme out of the active site. Residues E347 and E351, which are located on the π-helix and form an acidic path leading to a salt bridge interaction with the propionate of the heme, accelerate the release process.     

Aromatic stacking as a determinant of the thermal stability of CYP119 from Sulfolobus solfataricus

Two notable features of the thermophilic CYP119, an Arg154-Glu212 salt bridge between the F-G loop and the I helix and an extended aromatic cluster, were studied to determine their contributions to the thermal stability of the enzyme. Site-specific mutants of the salt bridge (Arg154, Glu212) and aromatic cluster (Tyr2, Trp4, Trp231, Tyr250, Trp281) were expressed and purified. The substrate-binding and kinetic constants for lauric acid hydroxylation are little affected in most mutants, but the E212D mutant is inactive and the R154Q mutant has higher K(s),K(m), and k(cat) values. The salt bridge mutants, like wild-type CYP119, melt at 91+/-1 degrees C, whereas mutation of individual residues in the extended aromatic cluster lowers the T(m) by 10-15 degrees C even though no change is observed on mutation of an unrelated aromatic residue. The extended aromatic cluster, but not the Arg154-Glu212 salt bridge, contributes to the thermal stability of CYP119.     

A shift in the salt bridge interaction of residues D620 and E621 mediates the constitutive activation of Jak2-H538Q/K539L

Jak2 mutations in the exon 14 and exon 12 regions that cause constitutive activation have been associated with myeloproliferative neoplasms. We have previously shown that a pi stacking interaction between F617 and F595 is important for the constitutive activation of Jak2-V617F (Gnanasambandan et al., Biochemistry 49:9972-9984, 2010). Here, using a combination of molecular dynamics (MD) simulations and in vitro mutagenesis, we studied the molecular mechanism for the constitutive activation of the Jak2 exon 12 mutation, H538Q/K539L. The activation levels of Jak2-H538Q/K539L were found to be similar to that of Jak2-V617F, and Jak2-H538Q/K539L/V617F. Data from MD simulations indicated a shift in the salt bridge interactions of D620 and E621 with K539 in Jak2-WT to R541 in Jak2-H538Q/K539L. When compared to Jak2-WT, K539A mutation resulted in increased activation, while K539D or K539E mutations diminished Jak2 activation by 50 %. In the context of Jak2-H538Q/K539L, R541A mutation reduced its activation by 50 %, while R541D and R541E mutations returned its activation levels to that of Jak2-WT. Collectively, these results indicate that a shift in the salt bridge interaction of D620 and E621 with K539 in Jak2-WT to R541 in Jak2-H538Q/K539L is critical for constitutive activation of this Jak2 exon 12 mutant.     

Glu88 in the non-catalytic domain of acylpeptide hydrolase plays dual roles: charge neutralization for enzymatic activity and formation of salt bridge for thermodynamic stability

Acylpeptide hydrolase of Aeropyrum pernix K1 is composed of a catalytic alpha/beta hydrolase domain and a non-catalytic beta-propeller domain. The Glu88 residue of the propeller domain is highly conserved in the prolyl oligopeptidase family and forms an inter-domain salt bridge with Arg526, a key residue for substrate binding. We have dissected the functions of Glu88 using site-directed mutagenesis, steady-state kinetics analyses, and molecular dynamics simulations. In E88A and E88A/R526K mutants, with a broken inter-domain salt bridge and a positive charge at position 526, catalytic activities for both a peptidase substrate and an esterase substrate were almost abolished. Analysis of the pH dependence of the mutants' reaction kinetics indicates that these mutations lead to changes in the electrostatic environment of the active site, which can be modulated by chloride ions. These findings indicate that the neutralization at position 526 is favorable for the activity of the enzyme, which is also verified by the catalytic behavior of E88A/R526V mutant. All mutants have lower thermodynamic stability than the wild-type. Therefore, Glu88 plays two major roles in the function of the enzyme: neutralizing the positive charge of Arg526, thereby increasing the enzymatic activity, and forming the Glu88-Arg526 salt bridge, thereby stabilizing the protein.     

Latch and trigger role for R445 in DAT transport explains molecular basis of DTDS

A recent study reports on five different mutations as sources of dopamine transporter (DAT) deficiency syndrome (DTDS). One of these mutations, R445C, is believed to be located on the intracellular side of DAT distal to the primary (S1) or secondary (S2) sites to which substrate binding is understood to occur. Thus, the molecular mechanism by which the R445C mutation results in DAT transport deficiency has eluded explanation. However, the recently reported X-ray structures of the endogenous amine transporters for dDAT and hSERT revealed the presence of a putative salt bridge between R445 and E428 suggesting a possible mechanism. To evaluate whether the R445C effect is a result of a salt bridge interaction, the mutants R445E, E428R, and the double mutant E428R/R445E were generated. The single mutants R445E and E428R displayed loss of binding and transport properties of the substrate [³H]DA and inhibitor [³H]CFT at the cell surface while the double mutant E428R/R445E, although nonfunctional, restored [³H]DA and [³H]CFT binding affinity to that of WT. Structure based analyses of these results led to a model wherein R445 plays a dual role in normal DAT function. R445 acts as a component of a latch in its formation of a salt bridge with E428 which holds the primary substrate binding site (S1) in place and helps enforce the inward closed protein state. When this salt bridge is broken, R445 acts as a trigger which disrupts a local polar network and leads to the release of the N-terminus from its position inducing the inward closed state to one allowing the inward open state. In this manner, both the loss of binding and transport properties of the R445C variant are explained.     

Molecular dynamics simulations of lignin peroxidase in solution

The dynamical and structural properties of lignin peroxidase and its Trp171Ala mutant have been investigated in aqueous solution using molecular dynamics (MD) simulations. In both cases, the enzyme retained its overall backbone structure and all its noncovalent interactions in the course of the MD simulations. Very interestingly, the analysis of the MD trajectories showed the presence of large fluctuations in correspondence of the residues forming the heme access channel; these movements enlarge the opening and facilitate the access of substrates to the enzyme active site. Moreover, steered molecular dynamics docking simulations have shown that lignin peroxidase natural substrate (veratryl alcohol) can easily approach the heme edge through the access channel.     

Inter-subunit disulfide locking of the human P2X3 receptor elucidates ectodomain movements associated with channel gating

P2X3 receptors (P2X3R) are trimeric ATP-gated cation channels involved in sensory neurotransmission and inflammatory pain. We used homology modeling and molecular dynamic simulations of the hP2X3R to identify inter-subunit interactions of residues that are instrumental to elucidate conformational changes associated with gating of the hPX3R. We identified an ionic interaction between E112 and R198 of the head domain and dorsal fin domain, respectively, and E57 and T263 of the lower body domains of adjacent subunits and detected a marked rearrangement of these domains during gating of the hP3X3R. Double-mutant cycle analysis of the inter-subunit residue pairs E112/R198 and E57/T263 revealed significant interaction-free energies. Disulfide locking of the hP2X3R E112C/R198C or the E57C/T263C double cysteine mutants markedly reduced the ATP-induced current responses. The decreased current amplitude following inter-subunit disulfide cross-linking indicates that disulfide locking of the head and dorsal fin domains or at the level of the lower body domains of the hP2X3R prevents the gating-induced conformational rearrangement of the subunits with respect to each other. The distinct reorganization of the subunit interfaces during gating of the hP2X3R is generally consistent with the gating mechanism of other P2XRs. Charge-reversal mutagenesis and methanethiosulfonate (MTS)-modification of substituted cysteines demonstrated that E112 and R198 interact electrostatically. Both disulfide locking and salt bridge breaking of the E112/R198 interaction reduced the hP2X3R function. We conclude that the inter-subunit salt bridge between E112 and R198 of the head and dorsal fin domains, respectively, serves to control the mobility of these domains during agonist-activation of the hP2X3R.     

Evolution of flavone synthase I from parsley flavanone 3beta-hydroxylase by site-directed mutagenesis

Flavanone 3beta-hydroxylase (FHT) and flavone synthase I (FNS I) are 2-oxoglutarate-dependent dioxygenases with 80% sequence identity, which catalyze distinct reactions in flavonoid biosynthesis. However, FNS I has been reported exclusively from a few Apiaceae species, whereas FHTs are more abundant. Domain-swapping experiments joining the N terminus of parsley (Petroselinum crispum) FHT with the C terminus of parsley FNS I and vice versa revealed that the C-terminal portion is not essential for FNS I activity. Sequence alignments identified 26 amino acid substitutions conserved in FHT versus FNS I genes. Homology modeling, based on the related anthocyanidin synthase structure, assigned seven of these amino acids (FHT/FNS I, M106T, I115T, V116I, I131F, D195E, V200I, L215V, and K216R) to the active site. Accordingly, FHT was modified by site-directed mutagenesis, creating mutants encoding from one to seven substitutions, which were expressed in yeast (Saccharomyces cerevisiae) for FNS I and FHT assays. The exchange I131F in combination with either M106T and D195E or L215V and K216R replacements was sufficient to confer some FNS I side activity. Introduction of all seven FNS I substitutions into the FHT sequence, however, caused a nearly complete change in enzyme activity from FHT to FNS I. Both FHT and FNS I were proposed to initially withdraw the beta-face-configured hydrogen from carbon-3 of the naringenin substrate. Our results suggest that the 7-fold substitution affects the orientation of the substrate in the active-site pocket such that this is followed by syn-elimination of hydrogen from carbon-2 (FNS I reaction) rather than the rebound hydroxylation of carbon-3 (FHT reaction).     

Surface salt bridges contribute to the extreme thermal stability of an FN3-like domain from a thermophilic bacterium

This study uses differential scanning calorimetry, X-ray crystallography, and molecular dynamics simulations to investigate the structural basis for the high thermal stability (melting temperature 97.5°C) of a FN3-like protein domain from thermophilic bacteria Thermoanaerobacter tengcongensis (FN3tt). FN3tt adopts a typical FN3 fold with a three-stranded beta sheet packing against a four-stranded beta sheet. We identified three solvent exposed arginine residues (R23, R25, and R72), which stabilize the protein through salt bridge interactions with glutamic acid residues on adjacent strands. Alanine mutation of the three arginine residues reduced melting temperature by up to 22°C. Crystal structures of the wild type (WT) and a thermally destabilized (?Tm ?19.7°C) triple mutant (R23L/R25T/R72I) were found to be nearly identical, suggesting that the destabilization is due to interactions of the arginine residues. Molecular dynamics simulations showed that the salt bridge interactions in the WT were stable and provided a dynamical explanation for the cooperativity observed between R23 and R25 based on calorimetry measurements. In addition, folding free energy changes computed using free energy perturbation molecular dynamics simulations showed high correlation with melting temperature changes. This work is another example of surface salt bridges contributing to the enhanced thermal stability of thermophilic proteins. The molecular dynamics simulation methods employed in this study may be broadly useful for in silico surface charge engineering of proteins.     

Engineering out motion: introduction of a de novo disulfide bond and a salt bridge designed to close a dynamic cleft on the surface of cytochrome b5

A previous molecular dynamics (MD) simulation of cytochrome b5 (cyt b5) at 25 degrees C displayed localized dynamics on the surface of the protein giving rise to the periodic formation of a cleft that provides access to the heme through a protected hydrophobic channel [Storch and Daggett (1995) Biochemistry 34, 9682]. Here we describe the production and testing of mutants designed to prevent the cleft from opening using a combination of experimental and theoretical techniques. Two mutants have been designed to close the surface cleft: S18D to introduce a salt bridge and S18C:R47C to incorporate a disulfide bond. The putative cleft forms between two separate cores of the protein: one is structural in nature and can be monitored through the fluorescence of Trp 22, and the other binds the heme prosthetic group and can be tracked via heme absorbance. An increase in motion localized to the cleft region was observed for each protein, except for the disulfide-containing variant, in MD simulations at 50 degrees C compared to simulations at 25 degrees C. For the disulfide-containing variant, the cleft remained closed. Both urea and temperature denaturation curves were nearly identical for wild-type and mutant proteins when heme absorbance was monitored. In contrast, fluorescence studies revealed oxidized S18C:R47C to be considerably more stable based on the midpoints of the denaturation transitions, Tm and U1/2. Moreover, the fluorescence changes for each protein were complete at approximately 50 degrees C and a urea concentration of approximately 3.9 M, significantly below the temperature and urea concentration (62 degrees C, 5 M urea) required to observe heme release. In addition, solvent accessibility based on acrylamide quenching of Trp 22 was lower in the S18C:R47C mutant, particularly at 50 degrees C, before heme release [presented in the accompanying paper (58)]. The results suggest that a constraining disulfide bond can be designed to inhibit dynamic cleft formation on the surface of cyt b5. Located near the heme, the native dynamics of the cleft may be functionally important for protein-protein recognition and/or complex stabilization.     

Molecular basis for the substrate stereoselectivity in tryptophan dioxygenase

Tryptophan dioxygenase (TDO) and indoleamine 2,3-dioxygenase (IDO) are the only two heme proteins that catalyze the oxidation reaction of tryptophan (Trp) to N-formylkynurenine. While human IDO is able to oxidize both L- and D-Trp, human TDO (hTDO) displays major specificity for L-Trp. In this work, we aim to interrogate the molecular basis for the substrate stereoselectivity of hTDO. Our previous molecular dynamics simulation studies of Xanthomonas campestris TDO (xcTDO) showed that a hydrogen bond between T254 (T342 in hTDO) and the ammonium group of the substrate is present in the L-Trp-bound enzyme, but not in the D-Trp-bound enzyme. The fact that this is the only notable structural alteration induced by the change in the stereo structure of the substrate prompted us to produce and characterize the T342A mutant of hTDO to evaluate the structural role of T342 in controlling the substrate stereoselectivity of the enzyme. The experimental results indicate that the mutation only slightly perturbs the global structural properties of the enzyme but totally abolishes the substrate stereoselectivity. Molecular dynamics simulations of xcTDO show that T254 controls the substrate stereoselectivity of the enzyme by (i) modulating the hydrogen bonding interaction between the NH(3)(+) group and epoxide oxygen of the ferryl-indole 2,3-epoxide intermediate of the enzyme and (ii) regulating the dynamics of two active site loops, loop(250-260) and loop(117-130), critical for substrate binding.     

An elementary reaction step of the proton pump is revealed by mutation of tryptophan-164 to phenylalanine in cytochrome c oxidase from Paracoccus denitrificans

Cytochrome c oxidase couples reduction of dioxygen to water to translocation of protons over the inner mitochondrial or bacterial membrane. A likely proton acceptor for pumped protons is the Delta-propionate of heme a(3), which may receive the proton via water molecules from a conserved glutamic acid (E278 in subunit I of the Paracoccus denitrificans enzyme) and which receives a hydrogen bond from a conserved tryptophan, W164. Here, W164 was mutated to phenylalanine (W164F) to further explore the role of the heme a(3) Delta-propionate in proton translocation. FTIR spectroscopy showed changes in vibrations possibly attributable to heme propionates, and the midpoint redox potential of heme a(3) decreased by approximately 50 mV. The reaction of the oxidized W164F enzyme with hydrogen peroxide yielded substantial amounts of the intermediate F' even at high pH, which suggests that the mutation rearranges the local electric field in the binuclear center that controls the peroxide reaction. The steady-state proton translocation stoichiometry of the W164F enzyme dropped to approximately 0.5 H(+)/e(-) in cells and reconstituted proteoliposomes. Time-resolved electrometric measurements showed that when the fully reduced W164F enzyme reacted with O(2), the membrane potential generated in the fast phase of this reaction was far too small to account either for full proton pumping or uptake of a substrate proton from the inside of the proteoliposomes. Time-resolved optical spectroscopy showed that this fast electrometric phase occurred with kinetics corresponding to the transition from state A to P(R), whereas the subsequent transition to the F state was strongly delayed. This is due to a delay of reprotonation of E278 via the D-pathway, which was confirmed by observation of a slowed rate of Cu(A) oxidation and which explains the small amplitude of the fast charge transfer phase. Surprisingly, the W164F mutation thus mimics a weak block of the D-pathway, which is interpreted as an effect on the side chain isomerization of E278. The fast charge translocation may be due to transfer of a proton from E278 to a "pump site" above the heme groups and is likely to occur also in wild-type enzyme, though not distinguished earlier due to the high-amplitude membrane potential formation during the P(R)--> F transition.     

A case of primary selective hypoaldosteronism carrying three mutations in the aldosterone synthase (Cyp11b2) gene

An infant with a clinical phenotype of early onset hypoaldosteronism has been screened for mutation analysis of the Cyp11b2 gene encoding aldosterone synthase enzyme. We have described a novel nonsense mutation in exon 3 (c.508C>T) that gave rise to a shorter protein (Q170X) and two known concurrent missense mutations (c.594A>C in exon 3 and c.1157T>C in exon 7) that led to substitution of glutamic acid for aspartic acid at amino acid position 198 (E198D) and of valine for alanine at amino acid position 386 (V386A). The father, who carried E198D plus V386A mutations, showed a fractional sodium excretion of 1.25% that was unmodified by dietary salt restriction, suggesting a mild haploinsufficiency. We examined by in silico analysis the effect of the mutations on the secondary and tertiary structures of aldosterone synthase to explain the inefficient enzymatic activity. The Q170X mutation produced a truncated protein, which was consequently associated with a loss of catalytic activity. As predicted by JPred web system and Dock 6.3 software, the concurrent expression of E198D and V386A mutations induced a significant secondary structure rearrangement and a shift of the heme group and the 18-hydroxycorticosterone substrate from their optimal placement.     

Evaluation of the mechanism of aromatase cytochrome P450. A site-directed mutagenesis study.

Aromatase (CYP19) catalyzes three consecutive hydroxylation reactions converting C19 androgens to aromatic C18 estrogenic steroids. In this study, five human aromatase mutants (E302D, S478A, S478T, H480K, and H480Q) were prepared using a mammalian cell expression system. These mutants were evaluated by enzyme kinetic analysis, inhibitory profile studies, and reaction intermediate measurements. Three steroidal inhibitors [4-hydroxyandrostenedione (4-OHA), 7alpha-(4'-amino)phenylthio-1,4-androstandiene-3,17-dione (7alpha-APTADD), and bridge (2,19-methyleneoxy) androstene-3,17-dione (MDL 101003)], and four nonsteroidal inhibitors [aminoglutethimide (AG), CGS 20267, ICI D1033, and vorozole (R83842)] were used in the inhibitory profile studies. Our computer model of aromatase suggests that Glu302 is situated in the conserved I-helix region and located near the C-19 position of the steroid substrate. The model was supported by significant changes in kinetic parameters and a sevenfold increase in the Ki value of MDL 101,003 for the mutant E302D. As S478A was found to have kinetic properties similar to the wild-type enzyme and a much higher activity than S478T, Ser478 is thought to be situated in a rather restricted environment. There was a 10-fold increase in the Ki value of 7alpha-APTADD for S478T over that for the wild-type enzyme, suggesting that Ser478 might be near the C-7 position of the substrate. The reaction intermediate analysis revealed that significantly more 19-ol intermediate was generated by both S478A and S478T than the wild-type enzyme. These results would support a hypothesis that Ser478 plays a role in the first and second hydroxylation reactions. A positive charged amino acid is preferred at position 480 as shown by the fact that H480K has a significantly higher activity than H480Q. The Ki value of 4-OHA for H480Q was found to be three times that of the wild-type enzyme. In addition, significantly more 19-ol and 19-al intermediates were detected for both mutants H480K and H480Q than for the wild-type enzyme. Evaluation of the two mutations at His480 allows us to propose that this residue may participate in the aromatization reaction (the third step) by acting as a hydrogen bond donor for the C-3 keto group of the substrate. Furthermore, new products were generated when the enzyme was mutated at Ser478 and His480. Thus, these two residues must play an important role in the catalysis and are likely closer to the substrate binding site than previously predicted.     

Effect of single-site charge-reversal mutations on the catalytic properties of yeast cytochrome c peroxidase: evidence for a single, catalytically active, cytochrome c binding domain

Forty-six charge-reversal mutants of yeast cytochrome c peroxidase (CcP) have been constructed in order to determine the effect of localized charge on the catalytic properties of the enzyme. The mutants include the conversion of all 20 glutamate residues and 24 of the 25 aspartate residues in CcP, one at a time, to lysine residues. In addition, two positive-to-negative charge-reversal mutants, R31E and K149D, are included in the study. The mutants have been characterized by absorption spectroscopy and hydrogen peroxide reactivity at pH 6.0 and 7.5 and by steady-state kinetic studies using recombinant yeast iso-1 ferrocytochrome c (C102T) as substrate at pH 7.5. Many of the charge-reversal mutations cause detectable changes in the absorption spectrum of the enzyme reflecting increased amounts of hexacoordinate heme compared to wild-type CcP. The increase in hexacoordinate heme in the mutant enzymes correlates with an increase in H 2O 2-inactive enzyme. The maximum velocity of the mutants decreases with increasing hexacoordination of the heme group. Steady-state velocity studies indicate that 5 of the 46 mutations (R31E, D34K, D37K, E118K, and E290K) cause large increases in the Michaelis constant indicating a reduced affinity for cytochrome c. Four of the mutations occur within the cytochrome c binding site identified in the crystal structure of the 1:1 complex of yeast cytochrome c and CcP [Pelletier, H., and Kraut, J. (1992) Science 258, 1748-1755] while the fifth mutation site lies outside, but near, the crystallographic site. These data support the hypothesis that the CcP has a single, catalytically active cytochrome c binding domain, that observed in the crystal structures of the cytochrome c/CcP complex.     

Distal mutation modulates the heme sliding in mouse neuroglobin investigated by molecular dynamics simulation

Neuroglobin (Ngb), a hexa‐coordinated hemoprotein primarily expressed in the brain and retina, is thought to be involved in neuroprotection and signal transduction. Ngb can reversibly bind small ligands such as O₂ and CO to the heme iron by replacing the distal histidine which is bound to the iron as the endogenous ligand. In this work, molecular dynamics (MD) simulations were performed to investigate the functionally related structural properties and dynamical characteristics in carboxy mouse neuroglobin and three distal mutants including single mutants H64V, K67T and double mutant H64V/K67T. MD simulations suggest that the heme sliding motion induced by the binding of exogenous ligand is affected by the distal mutation obviously. Accompanying changes in loop flexibility and internal cavities imply the structural rearrangement of Ngb. Moreover, the solvent accessibility of heme and some crucial residues are influenced revealing an interactive network on the distal side. The work elucidates that the key residues K67 at E10 and H64 at E7 are significant in modulating the heme sliding and hence the structural and physiological function of Ngb. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Effects of noncovalent and covalent FAD binding on the redox and catalytic properties of p-cresol methylhydroxylase

Each flavoprotein subunit (alpha or PchF) of the alpha(2)beta(2) flavocytochrome p-cresol methylhydroxylase (PCMH) from Pseudomonas putida contains FAD covalently attached to Tyr384. PCMH oxidizes p-cresol to 4-hydroxybenzyl alcohol, which is oxidized subsequently by PCMH to 4-hydroxybenzaldehyde. The Y384F mutant form of PchF (apo-PchF[Y384F]) displayed stoichiometric noncovalent FAD binding. PchF[Y384F]FAD associated with the cytochrome subunit (beta or PchC) (producing PCMH[Y384F]), although not as avidly as with wild-type PchF containing covalently bound FAD (PchF(C)). Dramatic increases in the two-electron E(m,7) (NHE) values for FAD were observed when it bound noncovalently to either apo-PchF or apo-PchF[Y384F], and the two-electron E(m,7) value for FAD was increased further by about 75 mV upon covalent binding to PchF, i.e., PchF(C). The E(m,7) values increased by approximately 20 and 45 mV, respectively, when PchF(C) and PchF[Y384F]FAD associated with PchC. The two-electron E(m,7) for covalently bound FAD in PCMH is 84 mV, the highest measured for a flavoprotein. The values for the one-electron redox potentials (E(m,7), NHE) for FAD were measured also for various forms of PchF. Under anaerobiosis, the reduction of PchF[Y384F]FAD by substrates was similar to that observed previously for PchF containing noncovalently bound FAD. Stopped-flow kinetic studies indicated a rapid substrate reduction of the FAD and heme in PCMH[Y384F] which produced PchF[Y384F]FAD(rad) x PchC, the mutant enzyme containing the flavin radical and reduced heme. These experiments also revealed a slow reduction of unassociated PchC(ox) by PchF[Y384F]FAD(rad) x PchC. Steady-state kinetic studies of the reaction of PCMH[Y384F] with p-cresol indicated that the K(m) for this substrate was unchanged relative to that of PCMH, but that the k(cat) was diminished by an order of magnitude. The data indicate that the covalent attachment of FAD to PchF assists catalysis by raising the E(m,7) of the flavin. Contributions to this effect likely result from conformational changes.     

Interdomain salt‐bridges in the Ebola virus protein VP40 and their role in domain association and plasma membrane localization

The Ebola virus protein VP40 is a transformer protein that possesses an extraordinary ability to accomplish multiple functions by transforming into various oligomeric conformations. The disengagement of the C‐terminal domain (CTD) from the N‐terminal domain (NTD) is a crucial step in the conformational transformations of VP40 from the dimeric form to the hexameric form or octameric ring structure. Here, we use various molecular dynamics (MD) simulations to investigate the dynamics of the VP40 protein and the roles of interdomain interactions that are important for the domain–domain association and dissociation, and report on experimental results of the behavior of mutant variants of VP40. The MD studies find that various salt‐bridge interactions modulate the VP40 domain dynamics by providing conformational specificity through interdomain interactions. The MD simulations reveal a novel salt‐bridge between D45‐K326 when the CTD participates in a latch‐like interaction with the NTD. The D45‐K326 salt‐bridge interaction is proposed to help domain–domain association, whereas the E76‐K291 interaction is important for stabilizing the closed‐form structure. The effects of the removal of important VP40 salt‐bridges on plasma membrane (PM) localization, VP40 oligomerization, and virus like particle (VLP) budding assays were investigated experimentally by live cell imaging using an EGFP‐tagged VP40 system. It is found that the mutations K291E and D45K show enhanced PM localization but D45K significantly reduced VLP formation.     

Active site substitution A82W improves the regioselectivity of steroid hydroxylation by cytochrome P450 BM3 mutants as rationalized by spin relaxation nuclear magnetic resonance studies

Cytochrome P450 BM3 from Bacillus megaterium is a monooxygenase with great potential for biotechnological applications. In this paper, we present engineered drug-metabolizing P450 BM3 mutants as a novel tool for regioselective hydroxylation of steroids at position 16β. In particular, we show that by replacing alanine at position 82 with a tryptophan in P450 BM3 mutants M01 and M11, the selectivity toward 16β-hydroxylation for both testosterone and norethisterone was strongly increased. The A82W mutation led to a ≤42-fold increase in V(max) for 16β-hydroxylation of these steroids. Moreover, this mutation improves the coupling efficiency of the enzyme, which might be explained by a more efficient exclusion of water from the active site. The substrate affinity for testosterone increased at least 9-fold in M11 with tryptophan at position 82. A change in the orientation of testosterone in the M11 A82W mutant as compared to the orientation in M11 was observed by T(1) paramagnetic relaxation nuclear magnetic resonance. Testosterone is oriented in M11 with both the A- and D-ring protons closest to the heme iron. Substituting alanine at position 82 with tryptophan results in increased A-ring proton-iron distances, consistent with the relative decrease in the level of A-ring hydroxylation at position 2β.