An unusual orientation for Tyr75 in the active site of the aspartic proteinase from Saccharomyces cerevisiae

Exploring the pH-Dependent Structure-Dynamics-Function Relationship of Human Renin

Renin is a pepsin-like aspartyl protease and an important drug target for the treatment of hypertension; despite three decades' research, its pH-dependent structure-function relationship remains poorly understood. Here, we employed continuous constant pH molecular dynamics (CpHMD) simulations to decipher the acid/base roles of renin's catalytic dyad and the conformational dynamics of the flap, which is a common structural feature among aspartyl proteases. The calculated pK_a's suggest that catalytic Asp38 and Asp226 serve as the general base and acid, respectively, in agreement with experiment and supporting the hypothesis that renin's neutral optimum pH is due to the substrate-induced pK_a shifts of the aspartic dyad. The CpHMD data confirmed our previous hypothesis that hydrogen bond formation is the major determinant of the dyad pK_a order. Additionally, our simulations showed that renin's flap remains open regardless of pH, although a Tyr-inhibited state is occasionally formed above pH 5. These findings are discussed in comparison to the related aspartyl proteases, including β-secretases 1 and 2, cathepsin D, and plasmepsin II. Our work represents a first step toward a systematic understanding of the pH-dependent structure-dynamics-function relationships of pepsin-like aspartyl proteases that play important roles in biology and human disease states.

Heterologous expression and characterization of the aspartic endoprotease Pep4um from Ustilago maydis, a homolog of the human Chatepsin D, an important breast cancer therapeutic target

The pep4um gene (um04926) of Ustilago maydis encodes a protein related to either vacuolar or lysosomal aspartic proteases. Bioinformatic analysis of the Pep4um protein revealed that it is a soluble protein with a signal peptide suggesting that it likely passes through the secretory pathway, and it has two probable self-activation sites, which are similar to those in Saccharomyces cerevisiae PrA. Moreover, the active site of the Pep4um has the two characteristic aspartic acid residues of aspartyl proteases. The pep4um gene was cloned, expressed in Pichia pastoris and a 54 kDa recombinant protein was observed. Pep4um-rec was confirmed to be an aspartic protease by specifically inhibiting its enzymatic activity with pepstatin A. Pep4um-rec enzymatic activity on acidic hemoglobin was optimal at pH 4.0 and at 40 °C. To the best of our knowledge this is the first report about the heterologous expression of an aspartic protease from a basidiomycete. An in-depth in silico analysis suggests that Pep4um is homolog of the human cathepsin D protein. Thus, the Pep4um-rec protein may be used to test inhibitors of human cathepsin D, an important breast cancer therapeutic target.

Activation of Mst11 and Feedback Inhibition of Germ Tube Growth in Magnaporthe oryzae

Appressorium formation and invasive growth are two important steps in the infection cycle of Magnaporthe oryzae that are regulated by the Mst11-Mst7-Pmk1 mitogen-activated protein kinase (MAPK) pathway. However, the molecular mechanism involved in the activation of Mst11 MAPK kinase kinase is not clear in the rice blast fungus. In this study, we functionally characterized the regulatory region of Mst11 and its self-inhibitory binding. Deletion of the middle region of Mst11, which contains the Ras-association (RA) domain and two conserved phosphorylation sites (S453 and S458), blocked Pmk1 activation and appressorium formation. However, the MST11(ΔRA) transformant MRD-2 still formed appressoria, although it was reduced in virulence. Interestingly, over 50% of its germ tubes branched and formed two appressoria by 48 h, which was suppressed by treatments with exogenous cAMP. The G18V dominant active mutation enhanced the interaction of Ras2 with Mst11, suggesting that Mst11 has stronger interactions with the activated Ras2. Furthermore, deletion and site-directed mutagenesis analyses indicated that phosphorylation at S453 and S458 of Mst11 is important for appressorium formation and required for the activation of Pmk1. We also showed that the N-terminal region of Mst11 directly interacted with its kinase domain, and the S789G mutation reduced their interactions. Expression of the MST11(S789G) allele rescued the defect of the mst11 mutant in plant infection and resulted in the formation of appressoria on hydrophilic surfaces, suggesting the gain-of-function effect of the S789G mutation. Overall, our results indicate that the interaction of Mst11 with activated Ras2 and phosphorylation of S453 and S458 play regulatory roles in Mst11 activation and infection-related morphogenesis, possibly by relieving its self-inhibitory interaction between its N-terminal region and the C-terminal kinase domain. In addition, binding of Mst11 to Ras2 may be involved in the feedback inhibition of cAMP signaling and further differentiation of germ tubes after appressorium formation.

pH-Dependent Population Shift Regulates BACE1 Activity and Inhibition

BACE1, a major therapeutic target for treatment of Alzheimer's disease, functions within a narrow pH range. Despite tremendous effort and progress in the development of BACE1 inhibitors, details of the underlying pH-dependent regulatory mechanism remain unclear. Here we elucidate the pH-dependent conformational mechanism that regulates BACE1 activity using continuous constant-pH molecular dynamics (MD). The simulations reveal that BACE1 mainly occupies three conformational states and that the relative populations of the states shift according to pH. At intermediate pH, when the catalytic dyad is monoprotonated, a binding-competent state is highly populated, while at low and high pH a Tyr-inhibited state is dominant. Our data provide strong evidence supporting conformational selection as a major mechanism for substrate and peptide-inhibitor binding. These new insights, while consistent with experiment, greatly extend the knowledge of BACE1 and have implications for further optimization of inhibitors and understanding potential side effects of targeting BACE1. Finally, the work highlights the importance of properly modeling protonation states in MD simulations.

Crystal structure of a putative aspartic proteinase domain of the Mycobacterium tuberculosis cell surface antigen PE_PGRS16

We report the crystal structure of the first prokaryotic aspartic proteinase-like domain identified in the genome of Mycobacterium tuberculosis. A search in the genomes of Mycobacterium species showed that the C-terminal domains of some of the PE family proteins contain two classic DT/SG motifs of aspartic proteinases with a low overall sequence similarity to HIV proteinase. The three-dimensional structure of one of them, Rv0977 (PE_PGRS16) of M. tuberculosis revealed the characteristic pepsin-fold and catalytic site architecture. However, the active site was completely blocked by the N-terminal His-tag. Surprisingly, the enzyme was found to be inactive even after the removal of the N-terminal His-tag. A comparison of the structure with pepsins showed significant differences in the critical substrate binding residues and in the flap tyrosine conformation that could contribute to the lack of proteolytic activity of Rv0977.

The role of tyrosine 71 in modulating the flap conformations of BACE1

β-Site amyloid precursor protein cleaving enzyme 1 (BACE1) is a potential target for treating Alzheimer's disease. BACE1's binding site is partially covered by a flexible loop on its N-terminal domain, known as the "flap," which has been found in several conformations in crystal structures of BACE1 and other aspartyl proteases. The side chain of the invariant residue Tyr71 on the flap adopts several rotameric orientations, leading to our hypothesis that the orientation of this residue dictates the movement and conformations available to the flap. We investigated this hypothesis by performing 220 ns of molecular dynamics simulations of bound and unbound wild-type BACE1 as well as the unbound Y71A mutant. Our findings indicate that the flap exhibits various degrees of mobility and adopts different conformations depending on the Tyr71 orientation. Surprisingly, the "self-inhibited" form is stable in our simulations, making it a reasonable target for drug design. The alanine mutant, lacking a large side chain at position 71, displays significant differences in flap dynamics from wild type, freely sampling very open and closed conformations. Our simulations show that Tyr71, in addition to its previously determined functions in catalysis and substrate binding, has the important role of modulating flap conformations in BACE1.

Multi-scaled explorations of binding-induced folding of intrinsically disordered protein inhibitor IA3 to its target enzyme

Biomolecular function is realized by recognition, and increasing evidence shows that recognition is determined not only by structure but also by flexibility and dynamics. We explored a biomolecular recognition process that involves a major conformational change – protein folding. In particular, we explore the binding-induced folding of IA3, an intrinsically disordered protein that blocks the active site cleft of the yeast aspartic proteinase saccharopepsin (YPrA) by folding its own N-terminal residues into an amphipathic alpha helix. We developed a multi-scaled approach that explores the underlying mechanism by combining structure-based molecular dynamics simulations at the residue level with a stochastic path method at the atomic level. Both the free energy profile and the associated kinetic paths reveal a common scheme whereby IA3 binds to its target enzyme prior to folding itself into a helix. This theoretical result is consistent with recent time-resolved experiments. Furthermore, exploration of the detailed trajectories reveals the important roles of non-native interactions in the initial binding that occurs prior to IA3 folding. In contrast to the common view that non-native interactions contribute only to the roughness of landscapes and impede binding, the non-native interactions here facilitate binding by reducing significantly the entropic search space in the landscape. The information gained from multi-scaled simulations of the folding of this intrinsically disordered protein in the presence of its binding target may prove useful in the design of novel inhibitors of aspartic proteinases.

The intrinsically disordered peptide IA3 is the endogenous inhibitor for the enzyme named yeast aspartic proteinase saccharopepsin (YPrA). In the presence of YPrA, IA3 folds itself into an amphipathic helix that blocks the active site cleft of the enzyme. We developed a multi-scaled approach to explore the underlying mechanism of this binding-induced ordering transition. Our approach combines a structure-based molecular dynamics model at the residue level with a stochastic path method at the atomic level. Our simulations suggest that IA3 inhibits YPrA through an induced-fit mechanism where the enzyme (YPrA) induces conformational change of its inhibitor (IA3). This expands the definition of an induced-fit model from its original meaning that the binding of substrate (IA3) drives conformational change in the protein (YPrA). Our result is consistent with recent kinetic experiments and provides a microscopic explanation for the underlying mechanism. We also discuss the important roles of non-native interactions and backtracking. These results enrich our understanding of the enzyme-inhibition mechanism and may have value in the design of drugs.

Crystal structure of an active form of BACE1, an enzyme responsible for amyloid beta protein production

BACE1 (beta-secretase) is a transmembrane aspartic protease that cleaves the beta-amyloid precursor protein and generates the amyloid beta peptide (Abeta). BACE1 cycles between the cell surface and the endosomal system many times and becomes activated interconvertibly during its cellular trafficking, leading to the production of Abeta. Here we report the crystal structure of the catalytically active form of BACE1. The active form has novel structural features involving the conformation of the flap and subsites that promote substrate binding. The functionally essential residues and water molecules are well defined and play a key role in the iterative activation of BACE1. We further describe the crystal structure of the dehydrated form of BACE1, showing that BACE1 activity is dependent on the dynamics of a catalytically required Asp-bound water molecule, which directly affects its catalytic properties. These findings provide insight into a novel regulation of BACE1 activity and elucidate how BACE1 modulates its activity during cellular trafficking.

The structure and function of Saccharomyces cerevisiae proteinase A

Saccharomyces cerevisiae proteinase A (saccharopepsin; EC 3.4.23.25) is a member of the aspartic proteinase superfamily (InterPro IPR001969), which are proteolytic enzymes distributed among a variety of organisms. Targeted to the vacuole as a zymogen, its activation at acidic pH can occur by two different pathways, a one-step process to release mature proteinase A, involving the intervention of proteinase B, or a step-wise pathway via the autoactivation product known as pseudo-proteinase A. Once active, S. cerevisiae proteinase A is essential to the activities of other yeast vacuolar hydrolases, including proteinase B and carboxypeptidase Y. The mature enzyme is bilobal, with each lobe providing one of the two catalytically essential aspartic acid residues in the active site. The crystal structure of free proteinase A reveals that the flap loop assumes an atypical position, pointing directly into the S(1) pocket of the enzyme. With regard to hydrolysis, proteinase A has a preference for hydrophobic residues with Phe, Leu or Glu at the P1 position and Phe, Ile, Leu or Ala at P1', and is inhibited by IA(3), a natural and highly specific inhibitor produced by S. cerevisiae. This review is the first comprehensive review of S. cerevisiae PrA.

Functional Plasticity in the Substrate Binding Site of β-Secretase

The aspartic protease beta-secretase (BACE) cleaves the amyloid precursor protein into a 42 residue beta-peptide, which is the principal biochemical marker of Alzheimer's disease. Multiple explicit-water molecular dynamics simulations of the apo and inhibitor bound structures of BACE indicate that both open- and closed-flap conformations are accessible at room temperature and should be taken into account for inhibitor design. Correlated motion is observed within each of the two lobes of BACE, as well as for the interfacial region. A self-inhibited conformation with the side chain of Tyr71 occupying the S(1) pocket is present in some of the unbound simulations. The reversible loss of the side chain hydrogen bond between the catalytic Asp32 and Ser35, due to the concomitant reorientation of the Ser35 hydroxyl group and a water molecule conserved in pepsin-like enzymes, provides further evidence for the suggestion that Ser35 assists in proton acceptance and release by Asp32 during catalysis.

X-ray Snapshots of Peptide Processing in Mutants of Tricorn-interacting Factor F1 from Thermoplasma acidophilum

The tricorn-interacting factor F1 of the archaeon Thermoplasma acidophilum cleaves small hydrophobic peptide products of the proteasome and tricorn protease. F1 mutants of the active site residues that are involved in substrate recognition and catalysis displayed distinct activity patterns toward fluorogenic test substrates. Crystal structures of the mutant proteins complexed with peptides Phe-Leu, Pro-Pro, or Pro-Leu-Gly-Gly showed interaction of glutamates 213 and 245 with the N termini of the peptides and defined the S1 and S1' sites and the role of the catalytic residues. Evidence was found for processive peptide cleavage in the N-to-C direction, whereby the P1' product is translocated into the S1 site. A functional interaction of F1 with the tricorn protease was observed with the inactive F1 mutant G37A. Moreover, small angle x-ray scattering measurements for tricorn and inhibited F1 have been interpreted as formation of transient and substrate-induced complexes.

Crystal structure of cockroach allergen Bla g 2, an unusual zinc binding aspartic protease with a novel mode of self-inhibition

The crystal structure of Bla g 2 was solved in order to investigate the structural basis for the allergenic properties of this unusual protein. This is the first structure of an aspartic protease in which conserved glycine residues, in two canonical DTG triads, are substituted by different amino acid residues. Another unprecedented feature revealed by the structure is the single phenylalanine residue insertion on the tip of the flap, with the side-chain occupying the S1 binding pocket. This and other important amino acid substitutions in the active site region of Bla g 2 modify the interactions in the vicinity of the catalytic aspartate residues, increasing the distance between them to approximately 4A and establishing unique direct contacts between the flap and the catalytic residues. We attribute the absence of substantial catalytic activity in Bla g 2 to these unusual features of the active site. Five disulfide bridges and a Zn-binding site confer stability to the protein, which may contribute to sensitization at lower levels of exposure than other allergens.

Apo and inhibitor complex structures of BACE (beta-secretase)

Human BACE, also known as beta-secretase, shows promise as a potential therapeutic target for Alzheimer's disease. We determined the apo structure of BACE to 1.75 A, and a structure of a hydroxyethylamine inhibitor complex derived by soaking. These show significant active-site movements compared to previously described BACE structures. Additionally, the structures reveal two pockets that could be targeted by structure-based drug design.