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

Characterization of a Novel Milk-Clotting Aspartic Protease from Penicillium sp. and Structural Explanation for its High Milk-Clotting Index

A novel milk-clotting enzyme isolated from Penicillium sp. ACCC 39790 (PsMCE) was prepared by heterologous expression. The recombinant PsMCE had an apparent molecular mass of 45 kDa and exhibited maximum casein hydrolysis activity at pH 4.0 and 50 °C. The PsMCE activity was enhanced by calcium ions and strongly inhibited by pepstatin A. Through hydrolysis pattern and cleavage site analyses, the milk-clotting activity of PsMCE was related to its specific hydrolysis between Phe₁₀₅ and Met₁₀₆ in the κ-casein proteins. The structural basis of PsMCE was characterized using homology modeling, molecular docking, and interactional analysis. The P1' region of PsMCE is critical for its selective binding to the hydrolytic site in κ-casein, and the hydrophobic forces play a decisive role in the specific cleavage of Phe₁₀₅ and Met₁₀₆. These interactional analyses between PsMCE and the ligand peptide clarified the fundamentals of its high milk-clotting index (MCI). PsMCE could be applied in cheese making due to its thermolability and high MCI value as a potential milk-clotting enzyme.

Can Recombinant Tree Shrew (Tupaia belangeri chinensis) Chymosin Coagulate Cow (Bos taurus) Milk?

Abstract

Genetically engineered chymosin from the tree shrew (Tupaia belangeri chinensis) has been obtained and partially characterized for the first time. The target enzyme was produced in Escherichia coli, strain BL21(DE3). It was shown that tree shrew recombinant chymosin coagulates cow milk (Bos taurus). The total and specific milk-clotting activity of the obtained enzyme was 0.7–5.3 IMCU/mL and 8.8–16.6 IMCU/mg. The nonspecific proteolytic activity of tree shrew recombinant chymosin in relation to total bovine casein was 30 and 117% higher than that of recombinant chymosin of cow and of single-humped camel respectively. It was found that in comparison with most of the known genetically engineered chymosins, the tree shrew enzyme showed exceptionally low thermal stability. After heating at 45°C, the coagulation ability of tree shrew recombinant chymosin decreased by more than 40%, and at 50°C the enzyme lost more than 90% of the initial milk-clotting activity. The Michaelis constant (Km), enzyme turnover number (kcat), and catalytic efficiency (kcat/Km) for genetically engineered chymosin from the tree shrew were 6.3 ± 0.1 µM, 11 927 ± 3169 s–1 and 1968 ± 620 µM–1 s–1, respectively. Comparative analysis showed that the primary structure of the chymosin-sensitive site of cow kappa-casein and the supposed similar sequence of tree shrew kappa-casein differed by 75%. The ability of tree shrew recombinant chymosin to coagulate cow’s milk, along with a low thermal stability and high catalytic efficiency with respect to the substrate, imitating the chymosin-sensitive site of cow kappa-casein, suggests that this enzyme is of potential interest for cheese making.

Assessing the performance of the MM/PBSA and MM/GBSA methods. 6. Capability to predict protein-protein binding free energies and re-rank binding poses generated by protein-protein docking

Understanding protein-protein interactions (PPIs) is quite important to elucidate crucial biological processes and even design compounds that interfere with PPIs with pharmaceutical significance. Protein-protein docking can afford the atomic structural details of protein-protein complexes, but the accurate prediction of the three-dimensional structures for protein-protein systems is still notoriously difficult due in part to the lack of an ideal scoring function for protein-protein docking. Compared with most scoring functions used in protein-protein docking, the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) and Molecular Mechanics/Poisson Boltzmann Surface Area (MM/PBSA) methodologies are more theoretically rigorous, but their overall performance for the predictions of binding affinities and binding poses for protein-protein systems has not been systematically evaluated. In this study, we first evaluated the performance of MM/PBSA and MM/GBSA to predict the binding affinities for 46 protein-protein complexes. On the whole, different force fields, solvation models, and interior dielectric constants have obvious impacts on the prediction accuracy of MM/GBSA and MM/PBSA. The MM/GBSA calculations based on the ff02 force field, the GB model developed by Onufriev et al. and a low interior dielectric constant (εin = 1) yield the best correlation between the predicted binding affinities and the experimental data (rp = -0.647), which is better than MM/PBSA (rp = -0.523) and a number of empirical scoring functions used in protein-protein docking (rp = -0.141 to -0.529). Then, we examined the capability of MM/GBSA to identify the possible near-native binding structures from the decoys generated by ZDOCK for 43 protein-protein systems. The results illustrate that the MM/GBSA rescoring has better capability to distinguish the correct binding structures from the decoys than the ZDOCK scoring. Besides, the optimal interior dielectric constant of MM/GBSA for re-ranking docking poses may be determined by analyzing the characteristics of protein-protein binding interfaces. Considering the relatively high prediction accuracy and low computational cost, MM/GBSA may be a good choice for predicting the binding affinities and identifying correct binding structures for protein-protein systems.

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.

Mutations and seeding of amylin fibril-like oligomers

Seeding a protein solution with preformed fibrils can dramatically enhance the growth rate of amyloids. As the seeds do not need to be of the same protein, seeding may account for the observed correlations between amyloid diseases. In an effort to understand better the molecular mechanisms behind cross seeding we have studied in silico the effect of mutations on the seeding of amylin fibrils. Our investigations of the structural stability of decamers of wild type amylin peptides, of Y37L mutants, and of heteroassemblies of wild-type and mutant amylin molecules show that the experimentally observed efficient cross seeding can be explained based on similarity in fibril structure of components. We find that amyloids with similar side chains packing at the β-sheet interface are structurally compatible, acting as a good template for the congruent incorporation of homologues peptides. In the Y37L mutants, lack of tyrosine-specific interactions causes significant higher flexibility of the C terminal than observed in the wild-type fibril. This effects elongation of the mutant fibril leading to the longer lag times during aggregation that are observed in experiments. Our study gives guidelines for the design of ligands that could stabilize amylin fibrils.