Flexibility and enzymatic cold-adaptation: a comparative molecular dynamics investigation of the elastase family

Molecular dynamics simulations of representative mesophilic and psycrophilic elastases have been carried out at different temperatures to explore the molecular basis of cold adaptation inside a specific enzymatic family. The molecular dynamics trajectories have been compared and analyzed in terms of secondary structure, molecular flexibility, intramolecular and protein-solvent interactions, unravelling molecular features relevant to rationalize the efficient catalytic activity of psychrophilic elastases at low temperature. The comparative molecular dynamics investigation reveals that modulation of the number of protein-solvent interactions is not the evolutionary strategy followed by the psycrophilic elastase to enhance catalytic activity at low temperature. In addition, flexibility and solvent accessibility of the residues forming the catalytic triad and the specificity pocket are comparable in the cold- and warm-adapted enzymes. Instead, loop regions with different amino acid composition in the two enzymes, and clustered around the active site or the specificity pocket, are characterized by enhanced flexibility in the cold-adapted enzyme. Remarkably, the psycrophilic elastase is characterized by reduced flexibility, when compared to the mesophilic counterpart, in some scattered regions distant from the functional sites, in agreement with hypothesis suggesting that local rigidity in regions far from functional sites can be beneficial for the catalytic activity of psychrophilic enzymes.     

Characterization of an anionic trypsin from the eel (Anguilla japonica)

An enzyme, isolated from the pancreas of the eel Anguilla japonica and designated as anionic trypsin 1, had a molecular weight of 26,000 and an isoelectric point of 5.5. The amino acid composition of the enzyme was similar to that of bovine cationic trypsin as well as anionic trypsins from other species of fish. The enzyme was stable at pH 6 to 9 in the presence of calcium ions. Km and kcat values of the enzyme for N-tosyl-L-arginine methyl ester and N-tosyl-L-lysine methyl ester were quite similar to those of catfish anionic and bovine cationic trypsins.     

Electrostatic effects play a central role in cold adaptation of trypsin.

Organisms that live in constantly cold environments have to adapt their metabolism to low temperatures, but mechanisms of enzymatic adaptation to cold environments are not fully understood. Cold active trypsin catalyses reactions more efficiently and binds ligands more strongly in comparison to warm active trypsin. We have addressed this issue by means of comparative free energy calculations studying the binding of positively charged ligands to two trypsin homologues. Stronger inhibition of the cold active trypsin by benzamidine and positively charged P1-variants of BPTI is caused by rather subtle electrostatic effects. The different affinity of benzamidine originates solely from long range interactions, while the increased binding of P1-Lys and -Arg variants of BPTI is attributed to both long and short range effects that are enhanced in the cold active trypsin compared to the warm active counterpart. Electrostatic interactions thus provide an efficient strategy for cold adaptation of trypsin.     

Immunoreactive trypsins in sera from dogs before and after induction of experimental pancreatitis

Radioimmunoassays for anionic and cationic dog trypsins are described. Characterization of the immunoreactivities in sera from fasting dogs demonstrated the presence of the two proenzymes only. Fasting sera from 10 dogs contained anionic and cationic trypsinogen in concentrations between 17-110 micrograms/l and 7-19 micrograms/l, respectively. Induction of experimental pancreatitis in dogs was accompanied by a large increase of immunoreactive anionic and cationic trypsins in the circulation. During the progress of the pancreatitis, immunoreactive trypsin with larger molecular weight than trypsinogen appeared. This high molecular weight immunoreactive trypsin was not seen in serum after intravenous injection of pancreatic juice in dogs. The high molecular weight immunoreactive trypsin probably consists of trypsin in complex with protease inhibitors. In vitro studies showed that the immunoreactivity of trypsin decreased considerably after binding to alpha 1-antitrypsin or alpha-macroglobulins.     

Comparative thermal unfolding study of psychrophilic and mesophilic subtilisin-like serine proteases by molecular dynamics simulations

Molecular dynamics (MD) simulations of a subtilisin-like serine protease VPR from the psychrophilic marine bacterium Vibrio sp. PA-44 and its mesophilic homologue, proteinase K (PRK), have been performed for 20 ns at four different temperatures (300, 373, 473, and 573 K). The comparative analyses of MD trajectories reveal that at almost all temperatures, VPR exhibits greater structural fluctuations/deviations, more unstable regular secondary structural elements, and higher global flexibility than PRK. Although these two proteases follow similar unfolding pathways at high temperatures, VPR initiates unfolding at a lower temperature and unfolds faster at the same high temperatures than PRK. These observations collectively indicate that VPR is less stable and more heat-labile than PRK. Analyses of the structural/geometrical properties reveal that, when compared to PRK, VPR has larger radius of gyration (Rg), less intramolecular contacts and hydrogen bonds (HBs), more protein-solvent HBs, and smaller burial of nonpolar area and larger exposure of polar area. These suggest that the increased flexibility of VPR would be most likely caused by its reduced intramolecular interactions and more favourable protein-solvent interactions arising from the larger exposure of the polar area, whereas the enhanced stability of PRK could be ascribed to its increased intramolecular interactions arising from the better optimized hydrophobicity. The factors responsible for the significant differences in local flexibility between these two proteases were also analyzed and ascertained. This study provides insights into molecular basis of thermostability of homologous serine proteases adapted to different temperatures.     

Cold adaption of enzymes: Structural comparison between salmon and bovine trypsins

The crystal structure of an anionic form of salmon trypsin has been determined at 1.82 Å resolution. We report the first structure of a trypsin from a phoikilothermic organism in a detailed comparison to mammalian trypsins in order to look for structural rationalizations for the cold-adaption features of salmon trypsin. This form of salmon trypsin (T II) comprises 222 residues, and is homologous to bovine trypsin (BT) in about 65% of the primary structure. The tertiary structures are similar, with an overall displacement in main chain atomic positions between salmon trypsin and various crystal structures of bovine trypsin of about 0.8 Å. Intramolecular hydrogen bonds and hydrophobic interactions are compared and discussed in order to estimate possible differences in molecular flexibility which might explain the higher catalytic efficiency and lower thermostability of salmon trypsin compared to bovine trypsin. No overall differences in intramolecular interactions are detected between the two structures, but there are differences in certain regions of the structures which may explain some of the observed differences in physical properties. The distribution of charged residues is different in the two trypsins, and the impact this might have on substrate affinity has been discussed. © 1994 Wiley-Liss, Inc.     

Novel cationic triazine dyes for protein purification

The effectiveness of a new immobilized cationic triazine dye was investigated alongside two new amphoteric triazine dyes and two well known anionic triazine dyes, Procion Red H-3B and Procion Blue H-B, as chromatographic media for binding four familiar proteases-trypsin, chymotrypsin, thrombin and carboxypeptidase-B-as well as a typical oxidoreductase, lactate dehydrogenase, and human serum albumin. The new affinity adsorbent, CL-Sepharose-immobilized Cationic Dye, specifically binds trypsin-like proteases such as trypsin, thrombin, and carboxypeptidase-B, but none of the other proteins tested. In contrast, the amphoteric and anionic immobilized dyes bind all the other proteins tested in a similar fashion. The specificity of the cationic dye was exploited in the resolution of trypsin and chymotrypsin from a crude activated bovine pancreatic extract. The procedure described here affords trypsin with specific activity of 7400 units/mg with a 79% overall yield in a single step. The immobilized cationic dye, unlike previously reported adsorbents for trypsin, is inexpensive, readily synthesized, and displays a workable capacity of 4000 trypsin units or 0.55 mg protein/g moist weight gel (1.2 mumol dye/g moist weight gel) from a crude bovine pancreatic extract and, thus, is potentially amenable to process-scale operations.     

Pancreatic anionic trypsin: evidence for the existence of a 30 kDa form.

1. An anionic form of trypsin has been isolated from pancreas of various species (rat, pig, dog and cow). 2. The purification procedure included affinity chromatography on STI-Sepharose 4B and ion-exchange chromatography on DEAE-Sephadex A-50. 3. The preparation was homogenous as checked by SDS-polyacrylamide slab gel electrophoresis, resulting in an estimated molecular weight of 30 kilodaltons (kDa) for this anionic form. 4. Antibodies against the anionic form from rat pancreas cross-reacted towards the anionic enzyme from porcine pancreas but not with the dog or bovine enzyme, nor with all the studied cationic forms. 5. Limited proteolysis of tubulin, a cytoskeletal protein, with an anionic or cationic form of trypsin showed striking differences in the size of produced peptides.     

Activities of anionic and cationic trypsins in the temperature range from 5 to 37 degrees C. Mutant anionic trypsins as a model of cold-adapted (psychrophilic) enzymes

Temperature dependences of kinetic constants (k _cat and K _m) were studied for enzymatic hydrolysis of N-succinyl-L-alanyl-L-alanyl-L-prolyl-L-arginine-p-nitroanilide and N-succinyl-L-alanyl-L-alanyl-L-prolyl-L-lysine-p-nitroanilide by bovine cationic and rat anionic (wild-type and mutant) trypsins. The findings were compared with the corresponding literature data for hydrolysis of N-benzoyl-DL-arginine-p-nitroanilide by bovine cationic trypsin and natural trypsins of coldadapted fishes. The anionic and cationic trypsins were found to differ in organization of the S₁ -substrate-binding pocket. The difference in the binding of lysine and arginine residues to this site (S₁) was also displayed by opposite temperature dependences of hydrolysis constants for the corresponding substrates by the anionic and cationic trypsins. The data suggest that the effect of any factor on the binding of substrates (the K _m value) to the anionic and cationic trypsins and on the catalytic activity k _cat should be compared only with the corresponding data for the natural enzyme of the same type. Mutants of rat anionic trypsin at residues K188 or Y228 were prepared by site-directed mutagenesis as approximate models of natural psychrophilic trypsins. Substitution of the charged lysine residue in position 188 by hydrophobic phenylalanine residue shifted the pH optimum of the resulting mutant trypsin K188F from 8.0 to 9.0-10.0, similarly to the case of some natural psychrophilic trypsins, and also 1.5-fold increased its catalytic activity at low temperatures as compared to the wild-type enzyme.     

Leap-dynamics: efficient sampling of conformational space of proteins and peptides in solution

A molecular simulation scheme, called Leap-dynamics, that provides efficient sampling of protein conformational space in solution is presented. The scheme is a combined approach using a fast sampling method, imposing conformational 'leaps' to force the system over energy barriers, and molecular dynamics (MD) for refinement. The presence of solvent is approximated by a potential of mean force depending on the solvent accessible surface area. The method has been successfully applied to N-acetyl-L-alanine-N-methylamide (alanine dipeptide), sampling experimentally observed conformations inaccessible to MD alone under the chosen conditions. The method predicts correctly the increased partial flexibility of the mutant Y35G compared to native bovine pancreatic trypsin inhibitor. In particular, the improvement over MD consists of the detection of conformational flexibility that corresponds closely to slow motions identified by nuclear magnetic resonance techniques.     

Change in protein flexibility upon complex formation: analysis of Ras-Raf using molecular dynamics and a molecular framework approach

Changes in flexibility upon protein-protein complex formation of H-Ras and the Ras-binding domain of C-Raf1 have been investigated using the molecular framework approach FIRST (Floppy Inclusion and Rigid Substructure Topology) and molecular dynamics simulations (MD) of in total approximately 35 ns length. In a computational time of about one second, FIRST identifies flexible and rigid regions in a single, static three-dimensional molecular framework, whose vertices represent protein atoms and whose edges represent covalent and non-covalent (hydrogen bond and hydrophobic) constraints and fixed bond angles within the protein. The two methods show a very good agreement with respect to the identification of changes in flexibility in both binding partners on a local scale. This implies that flexibility can be successfully predicted by identifying which bonds limit motion within a molecule and how they are coupled. In particular, as identified by MD, the beta-sheet in Raf shows considerably more pronounced orientational correlations in the bound state compared to the unbound state. Similarly, FIRST assigns the beta-sheet to the largest rigid cluster of the complex. Interestingly, FIRST allows us to identify that interactions across the interface (but not conformational changes upon complex formation) result in the observed rigidification. Since regions of the beta-sheet of Raf that do not interact directly with Ras become rigidified, this also demonstrates the long-range aspect to rigidity percolation. Possible implications of the change of flexibility of the Ras-binding domain of Raf on the activation of Raf upon complex formation are discussed. Finally, the sensitivity of FIRST results with respect to the representation of non-covalent interactions used as constraints is probed.     

Dynamic properties of extremophilic subtilisin-like serine-proteases

The investigation of the structural determinants of enzymatic temperature adaptation is a crucial pre-requisite both in terms of fundamental research and industrial applications to develop new biocatalysts active at different temperature ranges. In several cases, the differences related to cold- or warm-adaptation are related to subtle structural and aminoacidic differences at the molecular level, often hard to detect. In this context, we present a comparative study of psychrophilic, mesophilic and thermophilic subtilisin-like serine proteases by all-atom molecular dynamics (MD) simulations in explicit solvent using a multiple-replica approach. Our results strongly enforce the current view on localized flexibility in crucial functional regions for cold-adapted serine proteases and point out a different optimization and usage of salt-bridge interactions and networks in cold- and warm-adapted enzymes. The analyses allow to identify a subset of structural and dynamic features strictly associated to cold adaptation and which change from cold- to heat-active subtilisins. In particular, the thermophilic subtilisin presents a high affinity calcium binding site which is not structurally conserved in the mesophilic and psychrophilic counterparts, which, as it turns out from the MD analyses, at the same position show a stable salt bridge network and no stabilizing intra-molecular interactions, respectively. These aspects, along with differential flexibility in regions close to the active site or substrate binding pocket, can be an indication of evolution at this protein site toward a lower stability moving from high to low temperature conditions.     

Isolation and nucleotide sequence of cDNA clone for bovine pancreatic anionic trypsinogen. Structural identity within the trypsin family.

A cDNA clone encoding an anionic form of bovine trypsinogen was isolated from a pancreatic cDNA library. The corresponding 855-nucleotide mRNA contains a short 5' noncoding region of 8 nucleotides and a long 3' noncoding region of 56 nucleotides in addition to a poly(A) tail of at least 50 nucleotides. The deduced amino acid sequence for the anionic pretrypsinogen (247 residues) includes the N-terminal 15-amino-acid signal peptide followed by an 8-amino-acid activation peptide. The zymogen (232 residues) contains an additional C-terminal serine, compared with the amino acid sequence of bovine cationic trypsinogen. The identity between the anionic and cationic forms of bovine trypsinogen (65%) is lower than that existing between the anionic protein and other mammalian anionic trypsinogens (73-85%), suggesting that trypsin gene duplication in mammals occurred prior to the evolutionary events responsible for the species divergence. Bovine pancreatic anionic trypsin possesses all the key amino acids characteristic of the serine protease family.     

Properties of lactate dehydrogenase from the isopod, Saduria entomon

An anionic trypsin from pyloric caeca of chum salmon (Oncorhynchus keta) was purified by ammonium sulfate and acetone fractionation followed by affinity chromatography, gel-filtration, and DEAE-anion exchange chromatography. The apparent molecular mass was about 24 kDa as determined by SDS-PAGE. The anionic chum salmon trypsin was moderately active toward esterase substrates such as tosyl-L-arginine methyl ester and tosyl-L-lysine methyl ester. Its amidase activity for benzoyl-L-arginine p-nitroanilide was comparative to those of bovine and Streptomyces griseus trypsins. Kinetic characteristics of anionic chum salmon, bovine, and Streptomyces griseus trypsins toward inverse substrate (p-amidinophenyl ester) were compared. Inverse substrate behaved as a specific substrate for anionic chum salmon trypsin with specific binding, efficient acylation, and relatively slow deacylation.     

Interactions in vitro and in vivo between anionic and cationic porcine trypsin and porcine plasma proteinase inhibitors

A method for purifying porcine anionic and cationic trypsin is presented. Reaction mixtures with increasing amounts of the two porcine trypsins and porcine serum were studied in vitro to evaluate the relative importance of alpha 1-macroglobulin and alpha 2-macroglobulin as well as alpha 1-proteinase inhibitor in the rapid binding of porcine anionic and cationic trypsin. Porcine cationic trypsin was preferentially bound to alpha 1-macroglobulin, while anionic trypsin exhibited equal binding to both alpha-macroglobulins. Both trypsins were also bound by the alpha 1-proteinase inhibitor but not until alpha 1-macroglobulin approached saturation. Trypsin-alpha-macroglobulin complexes were cleared from plasma with a half-life of 6 min. For trypsin-alpha 1-proteinase inhibitor-complexes the half-life was 120 min. These findings are in accordance with results for other mammalian species, including man.     

A comparative study of trypsin specificity based on QM/MM molecular dynamics simulation and QM/MM GBSA calculation

Hydrogen bonding and polar interactions play a key role in identification of protein-inhibitor binding specificity. Quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) simulations combined with DFT and semi-empirical Hamiltonian (AM1d, RM1, PM3, and PM6) methods were performed to study the hydrogen bonding and polar interactions of two inhibitors BEN and BEN1 with trypsin. The results show that the accuracy of treating the hydrogen bonding and polar interactions using QM/MM MD simulation of PM6 can reach the one obtained by the DFT QM/MM MD simulation. Quantum mechanics/molecular mechanics generalized Born surface area (QM/MM-GBSA) method was applied to calculate binding affinities of inhibitors to trypsin and the results suggest that the accuracy of binding affinity prediction can be significantly affected by the accurate treatment of the hydrogen bonding and polar interactions. In addition, the calculated results also reveal the binding specificity of trypsin: (1) the amidinium groups of two inhibitors generate favorable salt bridge interaction with Asp189 and form hydrogen bonding interactions with Ser190 and Gly214, (2) the phenyl of inhibitors can produce favorable van der Waals interactions with the residues His58, Cys191, Gln192, Trp211, Gly212, and Cys215. This systematic and comparative study can provide guidance for the choice of QM/MM MD methods and the designs of new potent inhibitors targeting trypsin.     

Tyrosine Sulfation of Human Trypsin Steers S2’ Subsite Selectivity towards Basic Amino Acids

Human cationic and anionic trypsins are sulfated on Tyr154, a residue which helps to shape the prime side substrate-binding subsites. Here, we used phage display technology to assess the significance of tyrosine sulfation for the specificity of human trypsins. The prime side residues P1′–P4′ in the binding loop of bovine pancreatic trypsin inhibitor (BPTI) were fully randomized and tight binding inhibitor phages were selected against non-sulfated and sulfated human cationic trypsin. The selection pattern for the two targets differed mostly at the P2′ position, where variants selected against non-sulfated trypsin contained primarily aliphatic residues (Leu, Ile, Met), while variants selected against sulfated trypsin were enriched also for Arg. BPTI variants carrying Arg, Lys, Ile, Leu or Ala at the P2′ position of the binding loop were purified and equilibrium dissociation constants were determined against non-sulfated and sulfated cationic and anionic human trypsins. BPTI variants harboring apolar residues at P2′ exhibited 3–12-fold lower affinity to sulfated trypsin relative to the non-sulfated enzyme, whereas BPTI variants containing basic residues at P2′ had comparable affinity to both trypsin forms. Taken together, the observations demonstrate that the tyrosyl sulfate in human trypsins interacts with the P2′ position of the substrate-like inhibitor and this modification increases P2′ selectivity towards basic side chains.     

Molecular dynamics and binding specificity analysis of the bovine immunodeficiency virus BIV Tat-TAR complex

We have performed molecular dynamics (MD) simulations, with particle-mesh Ewald, explicit waters, and counterions, and binding specificity analyses using combined molecular mechanics and continuum solvent (MM-PBSA) on the bovine immunodeficiency virus (BIV) Tat peptide-TAR RNA complex. The solution structure for the complex was solved independently by Patel and co-workers and Puglisi and co-workers. We investigated the differences in both structures and trajectories, particularly in the formation of the U-A-U base triple, the dynamic flexibility of the Tat peptide, and the interactions at the binding interface. We observed a decrease in RMSD in comparing the final average RNA structures and initial RNA structures of both trajectories, which suggests the convergence of the RNA structures to a MD equilibrated RNA structure. We also calculated the relative binding of different Tat peptide mutants to TAR RNA and found qualitative agreement with experimental studies.     

Comparison of anionic and cationic trypsinogens: the anionic activation domain is more flexible in solution and differs in its mode of BPTI binding in the crystal structure

Unlike bovine cationic trypsin, rat anionic trypsin retains activity at high pH. This alkaline stability has been attributed to stabilization of the salt bridge between the N-terminal Ile16 and Asp194 by the surface negative charge (Soman K, Yang A-S, Honig B, Fletterick R., 1989, Biochemistry 28:9918-9926). The formation of this salt bridge controls the conformation of the activation domain in trypsin. In this work we probe the structure of rat trypsinogen to determine the effects of the surface negative charge on the activation domain in the absence of the Ile16-Asp194 salt bridge. We determined the crystal structures of the rat trypsin-BPTI complex and the rat trypsinogen-BPTI complex at 1.8 and 2.2 A, respectively. The BPTI complex of rat trypsinogen resembles that of rat trypsin. Surprisingly, the side chain of Ile16 is found in a similar position in both the rat trypsin and trypsinogen complexes, although it is not the N-terminal residue and cannot form the salt bridge in trypsinogen. The resulting position of the activation peptide alters the conformation of the adjacent autolysis loop (residues 142-153). While bovine trypsinogen and trypsin have similar CD spectra, the CD spectrum of rat trypsinogen has only 60% of the intensity of rat trypsin. This lower intensity most likely results from increased flexibility around two conserved tryptophans, which are adjacent to the activation domain. The NMR spectrum of rat trypsinogen contains high field methyl signals as observed in bovine trypsinogen. It is concluded that the activation domain of rat trypsinogen is more flexible than that of bovine trypsinogen, but does not extend further into the protein core.     

Molecular dynamics simulation of the P2Y14 receptor. Ligand docking and identification of a putative binding site of the distal hexose moiety

A rhodopsin-based homology model of the P2Y14 receptor was inserted into a phospholipid bilayer and refined by molecular dynamics (MD) simulation. The binding modes of several known agonists, namely UDP-glucose and its analogues, were proposed using automatic molecular docking combined with Monte Carlo Multiple Minimum calculations. Compared to other P2Y receptors, the P2Y14 receptor has an atypical binding mode of the nucleobase, ribose, and phosphate moieties. The diphosphate moiety interacts with only one cationic residue, namely Lys171 of EL2, while in other P2Y receptor subtypes three Arg or Lys residues interact with the phosphate chain. Two other conserved cationic residues, namely Arg253 (6.55) and Lys277 (7.35) of the P2Y14 receptor together with two anionic residues (Glu166 and Glu174, located in EL2), are likely involved in interactions with the distal hexose moiety.