Can enzymes adopt a self-inhibited form? Results of x-ray crystallographic studies of chymosin

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.

Allosteric-Activation Mechanism of Bovine Chymosin Revealed by Bias-Exchange Metadynamics and Molecular Dynamics Simulations

The aspartic protease, bovine chymosin, catalyzes the proteolysis of κ-casein proteins in milk. The bovine chymosin-κ-casein complex is of industrial interest as the enzyme is used extensively in the manufacturing of processed dairy products. The apo form of the enzyme adopts a self-inhibited conformation in which the side chain of Tyr77 occludes the binding site. On the basis of kinetic, mutagenesis, and crystallographic data, it has been widely reported that a HPHPH sequence in the P8-P4 residues of the natural substrate κ-casein acts as the allosteric activator, but the mechanism by which this occurs has not previously been elucidated due to the challenges associated with studying this process by experimental methods. Here we have employed two computational techniques, molecular dynamics and bias-exchange metadynamics simulations, to study the mechanism of allosteric activation and to compute the free energy surface for the process. The simulations reveal that allosteric activation is initiated by interactions between the HPHPH sequence of κ-casein and a small α-helical region of chymosin (residues 112-116). A small conformational change in the α-helix causes the side chain of Phe114 to vacate a pocket that may then be occupied by the side chain of Tyr77. The free energy surface for the self-inhibited to open transition is significantly altered by the presence of the HPHPH sequence of κ-casein.

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.

Hot-spot mapping of the interactions between chymosin and bovine κ-casein

Chymosin is a commercially important enzyme in the manufacturing of cheese. Chymosin cleaves the milk protein κ-casein, which initiates the clotting process. Recently, it has been shown that camel chymosin has superior enzymatic properties toward cow's milk, compared to bovine chymosin. The two enzymes possess a high degree of homology. There are only minor differences in the binding cleft; hence, these must be important for binding the substrate. Models for the binding of a 16 amino acid fragment, consisting of the chymosin-sensitive region of bovine κ-casein (97-112), to both enzymes have previously been presented. Computational alanine scanning for mutating 39 residues in the substrate and the bovine enzyme are presented herein, and warm- (ΔΔG > 1 kcal/mol) and hot-spot (ΔΔG > 2 kcal/mol) residues in the bovine enzyme are identified. These residues are relevant for site-directed mutagenesis, with the aim of modifying the binding affinity and in turn affecting the catalytic efficacy of the enzyme.

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.

Camel and bovine chymosin: the relationship between their structures and cheese-making properties

Bovine and camel chymosin are aspartic peptidases that are used industrially in cheese production. They cleave the Phe105-Met106 bond of the milk protein κ-casein, releasing its predominantly negatively charged C-terminus, which leads to the separation of the milk into curds and whey. Despite having 85% sequence identity, camel chymosin shows a 70% higher milk-clotting activity than bovine chymosin towards bovine milk. The activities, structures, thermal stabilities and glycosylation patterns of bovine and camel chymosin obtained by fermentation in Aspergillus niger have been examined. Different variants of the enzymes were isolated by hydrophobic interaction chromatography and showed variations in their glycosylation, N-terminal sequences and activities. Glycosylation at Asn291 and the loss of the first three residues of camel chymosin significantly decreased its activity. Thermal differential scanning calorimetry revealed a slightly higher thermal stability of camel chymosin compared with bovine chymosin. The crystal structure of a doubly glycosylated variant of camel chymosin was determined at a resolution of 1.6 Å and the crystal structure of unglycosylated bovine chymosin was redetermined at a slightly higher resolution (1.8 Å) than previously determined structures. Camel and bovine chymosin share the same overall fold, except for the antiparallel central β-sheet that connects the N-terminal and C-terminal domains. In bovine chymosin the N-terminus forms one of the strands which is lacking in camel chymosin. This difference leads to an increase in the flexibility of the relative orientation of the two domains in the camel enzyme. Variations in the amino acids delineating the substrate-binding cleft suggest a greater flexibility in the ability to accommodate the substrate in camel chymosin. Both enzymes possess local positively charged patches on their surface that can play a role in interactions with the overall negatively charged C-terminus of κ-casein. Camel chymosin contains two additional positive patches that favour interaction with the substrate. The improved electrostatic interactions arising from variation in the surface charges and the greater malleability both in domain movements and substrate binding contribute to the better milk-clotting activity of camel chymosin towards bovine milk.

Initial stage of cheese production: a molecular modeling study of bovine and camel chymosin complexed with peptides from the chymosin-sensitive region of κ-casein

Bovine chymosin has long been the preferred enzyme used to coagulate cow's milk, in the initial stage of cheese production, during which it cleaves a specific bond in the milk protein κ-casein. Recently, camel chymosin has been shown to have a 70% higher clotting activity toward cow's milk and, moreover, to cleave κ-casein more selectively. Bovine chymosin, on the other hand, is a poor clotting agent toward camel's milk. This paper reports a molecular modeling study aimed at understanding this disparity, based on homology modeling and molecular dynamics simulations using peptide fragments of κ-casein from cow and camel in both bovine and camel chymosin. The results show that the complex between bovine chymosin and the fragment of camel κ-casein is indeed less stable in the binding pocket. The results also indicate that this in part may be due to charge repulsion between a lysine residue in bovine chymosin and an arginine residue in the P4 position of camel κ-casein.

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.

Bovine chymosin: a computational study of recognition and binding of bovine kappa-casein

Bovine chymosin is an aspartic protease that selectively cleaves the milk protein kappa-casein. The enzyme is widely used to promote milk clotting in cheese manufacturing. We have developed models of residues 97-112 of bovine kappa-casein complexed with bovine chymosin, using ligand docking, conformational search algorithms, and molecular dynamics simulations. In agreement with limited experimental evidence, the model suggests that the substrate binds in an extended conformation with charged residues on either side of the scissile bond playing an important role in stabilizing the binding pose. Lys111 and Lys112 are observed to bind to the N-terminal domain of chymosin displacing a conserved water molecule. A cluster of histidine and proline residues (His98-Pro99-His100-Pro101-His102) in kappa-casein binds to the C-terminal domain of the protein, where a neighboring conserved arginine residue (Arg97) is found to be important for stabilizing the binding pose. The catalytic site (including the catalytic water molecule) is stable in the starting conformation of the previously proposed general acid/base catalytic mechanism for 18 ns of molecular dynamics simulations.

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.

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.

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.

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

The structures of the native Saccharomyces cerevisiae proteinase A have been solved by molecular replacement in the monoclinic and trigonal crystal forms and refined at 2.6-2.7A resolution. These structures agree overall with those of other uninhibited aspartic proteinases. However, an unusual orientation for the side chain of Tyr75, a conserved residue on the flexible "flap" that covers the active site and is important for the activity of these enzymes, was found in the trigonal crystals. A similar conformation of Tyr75 occupying the S1 substrate-binding pocket was previously reported only for chymosin (where it was interpreted as representing a "self-inhibited" state of the enzyme), but for no other aspartic proteinases. Since this orientation of Tyr75 has now been seen in the structures of two members of the family of aspartic proteinases, it might indicate that the placement of that residue in the S1 substrate-binding pocket might have some functional significance, analogous to what was seen for self-inhibited structures of serine proteinases.

Post X-ray crystallographic studies of chymosin: the existence of two structural forms and the regulation of activity by the interaction with the histidine-proline cluster of kappa-casein

Calf chymosin molecules exist in the two alternative structural forms: the first one has S1 and S3 binding pockets occluded by its own Tyr77 residue (the self-inhibited form); the second has these pockets free for a substrate binding (the active form). The preliminary incubation of the enzyme with a pentapeptide corresponding to the histidine-proline cluster of the specific substrate kappa-casein results in a 200-fold increase of the hydrolysis rate for the enzyme 'slow substrate'. The result suggests that the cluster is an allosteric effector that promotes the conversion of the enzyme into the active form. These data provide the experimental ground for the explanation of chymosin specificity towards kappa-casein.

Restrained molecular dynamics study of the interaction between bovine kappa-casein peptide 98-111 and bovine chymosin and porcine pepsin

The cleavage of bovine kappa-casein at the Phe105-Met106 bond by chymosin or pepsin is the first stage in casein micelle coagulation and casein digestion. The nature of the interaction of the peptide His98-Pro-His-Pro-His-Leu-Ser-Phe105-Met-Ala-Ile-Pro-Pro- Lys111 with chymosin and porcine pepsin was investigated using molecular modelling and energy minimization techniques. This study verified and extended a proposed model that electrostatic binding (involving His98, His100, His102 and Lys111 or Lys112) at either end of the active site cleft of chymosin is important for the positioning of residues 103-108 in the cleft. The peptide conformation remained unchanged in going from solution to binding into the active site cleft, with the exception that optimum binding of substrate to chymosin required the isomerization of the His98-Pro99 peptide bond from the trans to the cis conformation. The study also identified an acidic region in porcine pepsin that is in a position to form strong electrostatic interactions with the histidines at the N-terminus of the peptide.