Comparative structural analysis and substrate specificity engineering of the hyperthermostable beta-glucosidase CelB from Pyrococcus furiosus

The substrate specificity of the beta-glucosidase (CelB) from the hyperthermophilic archaeon Pyrococcus furiosus, a family 1 glycosyl hydrolase, has been studied at a molecular level. Following crystallization and X-ray diffraction of this enzyme, a 3.3 A resolution structural model has been obtained by molecular replacement. CelB shows a homo-tetramer configuration, with subunits having a typical (betaalpha)(8)-barrel fold. Its active site has been compared to the one of the previously determined 6-phospho-beta-glycosidase (LacG) from the mesophilic bacterium Lactococcus lactis. The overall design of the substrate binding pocket is very well conserved, with the exception of three residues that have been identified as a phosphate binding site in LacG. To verify the structural model and alter its substrate specificity, these three residues have been introduced at the corresponding positions in CelB (E417S, M424K, F426Y) in different combinations: single, double, and triple mutants. Characterization of the purified mutant CelB enzyme revealed that F426Y resulted in an increased affinity for galactosides, whereas M424K gave rise to a shifted pH optimum (from 5.0 to 6.0). Analysis of E417S revealed a 5-fold and a 3-fold increase of the efficiency of hydrolyzing o-nitrophenol-beta-D-galactopyranoside-6-phosphate, in the single and triple mutants, respectively. In contrast, their activity on nonphosphorylated sugars was largely reduced (30-300-fold). The residue at position E417 in CelB seems to be the determining factor for the difference in substrate specificity between the two types of family 1 glycosidases.     

Cargo sequences are important for Som1p-dependent signal peptide cleavage in yeast mitochondria

The inner membrane protease (IMP) has two catalytic subunits, Imp1p and Imp2p, that exhibit nonoverlapping substrate specificity in mitochondria of the yeast Saccharomyces cerevisiae. The IMP also has at least one noncatalytic subunit, Som1p, which is required to cleave signal peptides from a subset of Imp1p substrates. To understand how Som1p mediates Imp1p substrate specificity, we addressed the possibility that Som1p functions as a molecular chaperone, which binds to specific substrates and directs them to the catalytic site. Our results show that cargo sequences attached to the signal peptide are important for Som1p-dependent presequence cleavage; however, no specific cargo sequence is required. Indeed, we show that a substrate normally destined for Imp2p is cleaved in a Som1p-dependent manner when the substrate is directed to Imp1p. These results argue against the notion that Som1p is a molecular chaperone. Instead, we propose that the cargo of some Imp1p substrates can assume a conformation incompatible with presequence cleavage. Som1p could thus act through Imp1p to improve cleavage efficiency early during substrate maturation.     

Substitution of the heme binding module in hemoglobin alpha- and beta-subunits. Implication for different regulation mechanisms of the heme proximal structure between hemoglobin and myoglobin

In our previous work, we demonstrated that the replacement of the "heme binding module," a segment from F1 to G5 site, in myoglobin with that of hemoglobin alpha-subunit converted the heme proximal structure of myoglobin into the alpha-subunit type (Inaba, K., Ishimori, K. and Morishima, I. (1998) J. Mol. Biol. 283, 311-327). To further examine the structural regulation by the heme binding module in hemoglobin, we synthesized the betaalpha(HBM)-subunit, in which the heme binding module (HBM) of hemoglobin beta-subunit was replaced by that of hemoglobin alpha-subunit. Based on the gel chromatography, the betaalpha(HBM)-subunit was preferentially associated with the alpha-subunit to form a heterotetramer, alpha(2)[betaalpha(HBM)(2)], just as is native beta-subunit. Deoxy-alpha(2)[betaalpha(HBM)(2)] tetramer exhibited the hyperfine-shifted NMR resonance from the proximal histidyl N(delta)H proton and the resonance Raman band from the Fe-His vibrational mode at the same positions as native hemoglobin. Also, NMR spectra of carbonmonoxy and cyanomet alpha(2)[betaalpha(HBM)(2)] tetramer were quite similar to those of native hemoglobin. Consequently, the heme environmental structure of the betaalpha(HBM)-subunit in tetrameric alpha(2)[betaalpha(HBM)(2)] was similar to that of the beta-subunit in native tetrameric Hb A, and the structural conversion by the module substitution was not clear in the hemoglobin subunits. The contrastive structural effects of the module substitution on myoglobin and hemoglobin subunits strongly suggest different regulation mechanisms of the heme proximal structure between these two globins. Whereas the heme proximal structure of monomeric myoglobin is simply determined by the amino acid sequence of the heme binding module, that of tetrameric hemoglobin appears to be closely coupled to the subunit interactions.     

Functional differences between two major ubiquitin receptors in the proteasome; S5a and hRpn13

It is well known that S5a and hRpn13 are two major ubiquitin (Ub) receptors in the proteasome but little is known about their functional difference in recruiting ubiquitinated substrates. In this study using siRNA-mediated knockdown of S5a or hRpn13, we found that two Ub receptors had different substrate specificity although similar level of accumulation of high molecular weight Ub-conjugates was observed. Interesting enough, depletion of S5a, but not hRpn13, resulted in the Ub-containing aggregates and induced ER chaperones such as Grp78 and Grp94. ERAD substrates such as α-TCR and α1-antitrypsin were also stabilized by the depletion of S5a but not hRpn13. Our results suggest that there is different substrate specificity between S5a and hRpn13 at the level of delivery and S5a may be the major docking site for ERAD substrates.     

Endopolygalacturonases reveal molecular features for processivity pattern and tolerance towards acetylated pectin

Endopolygalacturonases (EndoPGs) hydrolyse the 1-4 linkages between two alpha-d-galacturonic acids (GalA) of the smooth homogalacturonan regions of pectin. GalA may be methyl-esterified on the carboxylic group and acetyl-esterified on the hydroxylic groups. EndoPG activity most often decreases with such increasing degree of substitution. In this paper, we used bioinformatics and molecular modelling technics to explain the tolerance profile at the molecular scale and processivity scheme of three endoPGs with respect to acetylated pectin substrate; the first two enzymes originate from Aspergillus niger (AnPGI and AnPGII) and the third from Fusarium moniliforme (FmPG). Partly acetylated and methylated homogalacturonan fragments in complex with the three PGs were successively modelled in silico. The amino acid residues involved in substrate binding were identified for each enzyme. Similarly, the docking pattern of the differently decorated oligomers in the catalytic groove was individually characterized for each enzyme. This work shows full agreement with our previous extensive mass spectrometry analysis of the hydrolytic products that established distinct tolerance profiles for the three endoPGs and earlier work that ascertained processivity, specifically for AnPGI. In our previous work, AnPGI was shown to be the most powerful enzyme among the three enzymes with an enhanced tolerance towards O2- and O3-acetylated substrates. We report here amino acids of AnPGI that are unique in binding the pectin backbone and that are identified as possibly crucial for its specificity, namely S191(An)(PGI)/D240(An)(PGI). Similarly, topologically equivalent residues in AnPGII and FmPG were identified that could impede such binding; S234(An)(PGII)/S91(An)(PGII) and S245(Fm)(PG)/V89(Fm)(PG). In addition, we report here, from normal mode analysis computed on AnPG1, a shear bending motion of 15 A of amplitude that fully accredits the processive action pattern for this enzyme, with D240(An)(PGI) and R96(An)(PGI) working as crampons to favour the sliding of the substrate. Conversely, the same method clearly evidences a hinge binding motion for AnPGII and FmPG that should only authorize one hydrolytic event per enzyme/substrate encounter.     

Directed evolution of (betaalpha)(8)-barrel enzymes

Natural molecular evolution supplies us with manifold examples of protein engineering. The imitation of these natural processes in the design of new enzymes has led to surprising and insightful results. Well-suited for design by evolutionary methods are enzymes with the common and versatile (betaalpha)(8)-barrel fold. Studies of enzyme stability, folding and design as well as the evolution of (betaalpha)(8)-barrel enzymes are discussed.     

Dynamic profiling of double-stranded RNA binding proteins

Double-stranded (ds) RNA is a key player in numerous biological activities in cells, including RNA interference, anti-viral immunity and mRNA transport. The class of proteins responsible for recognizing dsRNA is termed double-stranded RNA binding proteins (dsRBP). However, little is known about the molecular mechanisms underlying the interaction between dsRBPs and dsRNA. Here we examined four human dsRBPs, ADAD2, TRBP, Staufen 1 and ADAR1 on six dsRNA substrates that vary in length and secondary structure. We combined single molecule pull-down (SiMPull), single molecule protein-induced fluorescence enhancement (smPIFE) and molecular dynamics (MD) simulations to investigate the dsRNA-dsRBP interactions. Our results demonstrate that despite the highly conserved dsRNA binding domains, the dsRBPs exhibit diverse substrate specificities and dynamic properties when in contact with different RNA substrates. While TRBP and ADAR1 have a preference for binding simple duplex RNA, ADAD2 and Staufen1 display higher affinity to highly structured RNA substrates. Upon interaction with RNA substrates, TRBP and Staufen1 exhibit dynamic sliding whereas two deaminases ADAR1 and ADAD2 mostly remain immobile when bound. MD simulations provide a detailed atomic interaction map that is largely consistent with the affinity differences observed experimentally. Collectively, our study highlights the diverse nature of substrate specificity and mobility exhibited by dsRBPs that may be critical for their cellular function.     

Crystal structure of D-aminoacylase from Alcaligenes faecalis DA1. A novel subset of amidohydrolases and insights into the enzyme mechanism

D-Aminoacylase is an attractive candidate for commercial production of D-amino acids through its catalysis in the hydrolysis of N-acyl-D-amino acids. We report here the first D-aminoacylase crystal structure from A. faecalis at 1.5-A resolution. The protein comprises a small beta-barrel, and a catalytic (betaalpha)(8)-barrel with a 63-residue insertion. The enzyme structure shares significant similarity to the alpha/beta-barrel amidohydrolase superfamily, in which the beta-strands in both barrels superimpose well. Unexpectedly, the enzyme binds two zinc ions with widely different affinities, although only the tightly bound zinc ion is required for activity. One zinc ion is coordinated by Cys(96), His(220), and His(250), while the other is loosely chelated by His(67), His(69), and Cys(96). This is the first example of the metal ion coordination by a cysteine residue in the superfamily. Therefore, D-aminoacylase defines a novel subset and is a mononuclear zinc metalloenzyme but containing a binuclear active site. The preferred substrate was modeled into a hydrophobic pocket, revealing the substrate specificity and enzyme catalysis. The 63-residue insertion containing substrate-interacting residues may act as a gate controlling access to the active site, revealing that the substrate binding would induce a closed conformation to sequester the catalysis from solvent.     

Exceptional disfavor for proline at the P + 1 position among AGC and CAMK kinases establishes reciprocal specificity between them and the proline-directed kinases

To precisely regulate critical signaling pathways, two kinases that phosphorylate distinct sites on the same protein substrate must have mutually exclusive specificity. Evolution could assure this by designing families of kinase such as basophilic kinases and proline-directed kinase with distinct peptide specificity; their reciprocal peptide specificity would have to be very complete, since recruitment of substrate allows phosphorylation of even rather poor phosphorylation sites in a protein. Here we report a powerful evolutionary strategy that assures distinct substrates for basophilic kinases (PKA, PKG and PKC (AGC) and calmodulin-dependent protein kinase (CAMK)) and proline-directed kinase, namely by the presence or absence of proline at the P + 1 position in substrates. Analysis of degenerate and non-degenerate peptides by in vitro kinase assays reveals that proline at the P + 1 position in substrates functions as a "veto" residue in substrate recognition by AGC and CAMK kinases. Furthermore, analysis of reported substrates of two typical basophilic kinases, protein kinase C and protein kinase A, shows the lowest occurrence of proline at the P + 1 position. Analysis of crystal structures and sequence conservation provides a molecular basis for this disfavor and illustrate its generality.     

Bioinformatic approaches for predicting substrates of proteases

Proteases have central roles in "life and death" processes due to their important ability to catalytically hydrolyze protein substrates, usually altering the function and/or activity of the target in the process. Knowledge of the substrate specificity of a protease should, in theory, dramatically improve the ability to predict target protein substrates. However, experimental identification and characterization of protease substrates is often difficult and time-consuming. Thus solving the "substrate identification" problem is fundamental to both understanding protease biology and the development of therapeutics that target specific protease-regulated pathways. In this context, bioinformatic prediction of protease substrates may provide useful and experimentally testable information about novel potential cleavage sites in candidate substrates. In this article, we provide an overview of recent advances in developing bioinformatic approaches for predicting protease substrate cleavage sites and identifying novel putative substrates. We discuss the advantages and drawbacks of the current methods and detail how more accurate models can be built by deriving multiple sequence and structural features of substrates. We also provide some suggestions about how future studies might further improve the accuracy of protease substrate specificity prediction.     

Substrate specificity of Sphingobium chlorophenolicum 2,6-dichlorohydroquinone 1,2-dioxygenase

PcpA is an aromatic ring-cleaving dioxygenase that is homologous to the well-characterized Fe(II)-dependent catechol extradiol dioxygenases. This enzyme catalyzes the oxidative cleavage of 2,6-dichlorohydroquinone in the catabolism of pentachlorophenol by Sphingobium chlorophenolicum ATCC 39723. (1)H NMR and steady-state kinetics were used to determine the regiospecificity of ring cleavage and the substrate specificity of the enzyme. PcpA exhibits a high degree of substrate specificity for 2,6-disubstituted hydroquinones, with halogens greatly preferred at those positions. Notably, the k(cat)(app)/K(mA)(app) of 2,6-dichlorohydroquinone is ~40-fold higher than that of 2,6-dimethylhydroquinone. The asymmetric substrate 2-chloro-6-methylhydroquinone yields a mixture of 1,2- and 1,6-cleavage products. These two modes of cleavage have different K(mO(2))(app) values (21 and 260 μM, respectively), consistent with a mechanism in which the substrate binds in two catalytically productive orientations. In contrast, monosubstituted hydroquinones show a limited amount of ring cleavage but rapidly inactivate the enzyme in an O(2)-dependent fashion, suggesting that oxidation of the Fe(II) may be the cause. Potent inhibitors of PcpA include ortho-disubstituted phenols and 3-bromocatechol. 2,6-Dibromophenol is the strongest competitive inhibitor, consistent with PcpA's substrate specificity. Several factors that could yield this specificity for halogen substituents are discussed. Interestingly, 3-bromocatechol also inactivates the enzyme, while 2,6-dihalophenols do not, indicating a requirement for two hydroxyl groups for ring cleavage and for enzyme inactivation. These results provide mechanistic insights into the hydroquinone dioxygenases.     

Structural variation governs substrate specificity for organic anion transporter (OAT) homologs. Potential remote sensing by OAT family members

Organic anion transporters (OATs, SLC22) interact with a remarkably diverse array of endogenous and exogenous organic anions. However, little is known about the structural features that determine their substrate selectivity. We examined the substrate binding preferences and transport function of olfactory organic anion transporter, Oat6, in comparison with the more broadly expressed transporter, Oat1 (first identified as NKT). In analyzing interactions of both transporters with over 40 structurally diverse organic anions, we find a correlation between organic anion potency (pKi) and hydrophobicity (logP) suggesting a hydrophobicity-driven association with transporter-binding sites, which appears particularly prominent for Oat6. On the other hand, organic anion binding selectivity between Oat6 and Oat1 is influenced by the anion mass and net charge. Smaller mono-anions manifest greater potency for Oat6 and di-anions for Oat1. Comparative molecular field analysis confirms these mechanistic insights and provides a model for predicting new OAT substrates. By comparative molecular field analysis, both hydrophobic and charged interactions contribute to Oat1 binding, although it is predominantly the former that contributes to Oat6 binding. Together, the data suggest that, although the three-dimensional structures of these two transporters may be very similar, the binding pockets exhibit crucial differences. Furthermore, for six radiolabeled substrates, we assessed transport efficacy (Vmax) for Oat6 and Oat1. Binding potency and transport efficacy had little correlation, suggesting that different molecular interactions are involved in substrate binding to the transporter and translocation across the membrane. Substrate specificity for a particular transporter may enable design of drugs for targeting to specific tissues (e.g. olfactory mucosa). We also discuss how these data suggest a possible mechanism for remote sensing between OATs in different tissue compartments (e.g. kidney, olfactory mucosa) via organic anions.     

Specificity Profiling of Dual Specificity Phosphatase Vaccinia VH1-related (VHR) Reveals Two Distinct Substrate Binding Modes

Vaccinia VH1-related (VHR) is a dual specificity phosphatase that consists of only a single catalytic domain. Although several protein substrates have been identified for VHR, the elements that control the in vivo substrate specificity of this enzyme remain unclear. In this work, the in vitro substrate specificity of VHR was systematically profiled by screening combinatorial peptide libraries. VHR exhibits more stringent substrate specificity than classical protein-tyrosine phosphatases and recognizes two distinct classes of Tyr(P) peptides. The class I substrates are similar to the Tyr(P) motifs derived from the VHR protein substrates, having sequences of (D/E/φ)(D/S/N/T/E)(P/I/M/S/A/V)pY(G/A/S/Q) or (D/E/φ)(T/S)(D/E)pY(G/A/S/Q) (where φ is a hydrophobic amino acid and pY is phosphotyrosine). The class II substrates have the consensus sequence of (V/A)P(I/L/M/V/F)X_1–6pY (where X is any amino acid) with V/A preferably at the N terminus of the peptide. Site-directed mutagenesis and molecular modeling studies suggest that the class II peptides bind to VHR in an opposite orientation relative to the canonical binding mode of the class I substrates. In this alternative binding mode, the Tyr(P) side chain binds to the active site pocket, but the N terminus of the peptide interacts with the carboxylate side chain of Asp¹⁶⁴, which normally interacts with the Tyr(P) + 3 residue of a class I substrate. Proteins containing the class II motifs are efficient VHR substrates in vitro, suggesting that VHR may act on a novel class of yet unidentified Tyr(P) proteins in vivo.     

Molecular dynamics of substrate complexes with hamster cytochrome P450c17 (CYP17): mechanistic approach to understanding substrate binding and activities

The cytochrome P450c17 isoforms from various animal species have different substrate selectivity, especially for 17,20-lyase activity. In particular, the human P450c17 selectively produces dehydroepiandrosterone with little androstenedione (AD). Hamster P450c17, on the other hand, produces both of these steroids at comparable rates. We thus investigated if computational analysis could explain the difference in activity profiles. Therefore, we inserted the four P450c17 substrates-pregnenolone, progesterone, and their 17alpha-hydroxylated forms-inside our hamster P450c17 model, which we derived from our human P450c17 model based on the crystal structure of P450BMP. We performed molecular dynamics (MD) simulations on the complexes and analyzed the resultant trajectories to identify amino acids that interact with substrates. Starting with substrates in two different orientations, we obtained two sets of binding trajectories in each case. The first set of trajectories reveal structural rearrangements that occur during binding, whereas the second set of trajectories reflects substrate orientations during catalysis. Our modeling suggests that three distinct steps are required for substrate selectivity and binding to the hamster P450c17: (1) recognition of the substrate at the putative substrate entrance, characterized by a pocket at the surface of the hamster P450c17 containing charged residues R96 and D116; (2) entry of the substrate into the active site, in an intermediate position directed by possible hydrogen bonding of the substrates with the heme D-ring propionate group, R96, R440, and T306; followed by (3) 90 degrees counterclockwise rotation of the substrates, positioning them in optimal position for reactivity, a process that may be directed by hydrogen bonding to the 110-112 region of the hamster P450c17. With some substrates, we obtained trajectories which suggest that major distortions in the I-helix and opening of the H-I loop occur during substrate binding. In conclusion, these modeling exercises provide insight to possible structural reorganizations that occur during substrate binding and suggest that amino acids that participate in three distinct steps of this process may all contribute to substrate binding and activity.     

Engineering the substrate specificity of Staphylococcus aureus Sortase A. The beta6/beta7 loop from SrtB confers NPQTN recognition to SrtA

The Staphylococcus aureus transpeptidase Sortase A (SrtA) anchors virulence and colonization-associated surface proteins to the cell wall. SrtA selectively recognizes a C-terminal LPXTG motif, whereas the related transpeptidase Sortase B (SrtB) recognizes a C-terminal NPQTN motif. In both enzymes, cleavage occurs after the conserved threonine, followed by amide bond formation between threonine and the pentaglycine cross-bridge of cell wall peptidoglycan. Genetic and biochemical studies strongly suggest that SrtA and SrtB exhibit exquisite specificity for their recognition motifs. To better understand the origins of substrate specificity within these two isoforms, we used sequence and structural analysis to predict residues and domains likely to be involved in conferring substrate specificity. Mutational analyses and domain swapping experiments were conducted to test their function in substrate recognition and specificity. Marked changes in the specificity profile of SrtA were obtained by replacing the beta6/beta7 loop in SrtA with the corresponding domain from SrtB. The chimeric beta6/beta7 loop swap enzyme (SrtLS) conferred the ability to acylate NPQTN-containing substrates, with a k(cat)/K(m)(app) of 0.0062 +/- 0.003 m(-1) s(-1). This enzyme was unable to perform the transpeptidation stage of the reaction, suggesting that additional domains are required for transpeptidation to occur. The overall catalytic specificity profile (k(cat)/K(m)(app)(NPQTN)/k(cat)/K(m)(app)(LPETG)) of SrtLS was altered 700,000-fold from SrtA. These results indicate that the beta6/beta7 loop is an important site for substrate recognition in sortases.     

In silico analyses of substrate interactions with human serum paraoxonase 1

Human paraoxonase (HuPON1) is a serum enzyme that exhibits a broad spectrum of hydrolytic activities, including the hydrolysis of various organophosphates, esters, and recently identified lactone substrates. Despite intensive site-directed mutagenesis and other biological studies, the structural basis for the specificity of substrate interactions of HuPON1 remains elusive. In this study, we apply homology modeling, docking, and molecular dynamic (MD) simulations to probe the binding interactions of HuPON1 with representative substrates. The results suggest that the active site of HuPON1 is characterized by two distinct binding regions: the hydrophobic binding site for arylesters/lactones, and the paraoxon binding site for phosphotriesters. The unique binding modes proposed for each type of substrate reveal a number of key residues governing substrate specificity. The polymorphic residue R/Q192 interacts with the leaving group of paraoxon, suggesting it plays an important role in the proper positioning of this substrate in the active site. MD simulations of the optimal binding complexes show that residue Y71 undergoes an "open-closed" conformational change upon ligand binding, and forms strong interactions with substrates. Further binding free energy calculations and residual decomposition give a more refined molecular view of the energetics and origin of HuPON1/substrate interactions. These studies provide a theoretical model of substrate binding and specificity associated with wild type and mutant forms of HuPON1, which can be applied in the rational design of HuPON1 variants as bioscavengers with enhanced catalytic activity.     

Structural, kinetic, and thermodynamic studies of specificity designed HIV-1 protease

HIV-1 protease recognizes and cleaves more than 12 different substrates leading to viral maturation. While these substrates share no conserved motif, they are specifically selected for and cleaved by protease during viral life cycle. Drug resistant mutations evolve within the protease that compromise inhibitor binding but allow the continued recognition of all these substrates. While the substrate envelope defines a general shape for substrate recognition, successfully predicting the determinants of substrate binding specificity would provide additional insights into the mechanism of altered molecular recognition in resistant proteases. We designed a variant of HIV protease with altered specificity using positive computational design methods and validated the design using X-ray crystallography and enzyme biochemistry. The engineered variant, Pr3 (A28S/D30F/G48R), was designed to preferentially bind to one out of three of HIV protease's natural substrates; RT-RH over p2-NC and CA-p2. In kinetic assays, RT-RH binding specificity for Pr3 increased threefold compared to the wild-type (WT), which was further confirmed by isothermal titration calorimetry. Crystal structures of WT protease and the designed variant in complex with RT-RH, CA-p2, and p2-NC were determined. Structural analysis of the designed complexes revealed that one of the engineered substitutions (G48R) potentially stabilized heterogeneous flap conformations, thereby facilitating alternate modes of substrate binding. Our results demonstrate that while substrate specificity could be engineered in HIV protease, the structural pliability of protease restricted the propagation of interactions as predicted. These results offer new insights into the plasticity and structural determinants of substrate binding specificity of the HIV-1 protease.     

Substrate discrimination among mitogen-activated protein kinases through distinct docking sequence motifs

Mitogen-activated protein kinases (MAPKs) mediate cellular responses to a wide variety of extracellular stimuli. MAPK signal transduction cascades are tightly regulated, and individual MAPKs display exquisite specificity in recognition of their target substrates. All MAPK family members share a common phosphorylation site motif, raising questions as to how substrate specificity is achieved. Here we describe a peptide library screen to identify sequence requirements of the DEF site (docking site for ERK FXF), a docking motif separate from the phosphorylation site. We show that MAPK isoforms recognize DEF sites with unique sequences and identify two key residues on the MAPK that largely dictate sequence specificity. Based on these observations and computational docking studies, we propose a revised model for MAPK interaction with substrates containing DEF sites. Variations in DEF site sequence requirements provide one possible mechanism for encoding complex target specificity among MAPK isoforms.     

Enzyme specificity under dynamic control II: Principal component analysis of α-lytic protease using global and local solvent boundary conditions

The contributions of conformational dynamics to substrate specificity have been examined by the application of principal component analysis to molecular dynamics trajectories of alpha-lytic protease. The wild-type alpha-lytic protease is highly specific for substrates with small hydrophobic side chains at the specificity pocket, while the Met190-->Ala binding pocket mutant has a much broader specificity, actively hydrolyzing substrates ranging from Ala to Phe. Based on a combination of multiconformation analysis of cryo-X-ray crystallographic data, solution nuclear magnetic resonance (NMR), and normal mode calculations, we had hypothesized that the large alteration in specificity of the mutant enzyme is mainly attributable to changes in the dynamic movement of the two walls of the specificity pocket. To test this hypothesis, we performed a principal component analysis using 1-nanosecond molecular dynamics simulations using either a global or local solvent boundary condition. The results of this analysis strongly support our hypothesis and verify the results previously obtained by in vacuo normal mode analysis. We found that the walls of the wild-type substrate binding pocket move in tandem with one another, causing the pocket size to remain fixed so that only small substrates are recognized. In contrast, the M190A mutant shows uncoupled movement of the binding pocket walls, allowing the pocket to sample both smaller and larger sizes, which appears to be the cause of the observed broad specificity. The results suggest that the protein dynamics of alpha-lytic protease may play a significant role in defining the patterns of substrate specificity. As shown here, concerted local movements within proteins can be efficiently analyzed through a combination of principal component analysis and molecular dynamics trajectories using a local solvent boundary condition to reduce computational time and matrix size.