Correlation among sites of limited proteolysis, enzyme accessibility and segmental mobility

The relationship among accessibility to an enzyme, flexibility, and limited proteolysis was explored. Regions accessible to large probes, comparable in size to proteolytic enzymes, were computed in the crystallographic structures of thermolysin, trypsinogen and ribonuclease. Positions of these accessible regions were compared with sites of autolytic/proteolytic attacks, and with locations of flexible backbone segments. All the proteolytic sites were found to be exceptionally accessible. Most of them were also flexible, but at least one prominent site in trypsinogen appeared to be rigid. Thus, surface exposure seems to be more essential to proteolysis than flexibility.     

Global and local motions in ribonuclease A: a molecular dynamics study

The understanding of protein dynamics is one of the major goals of structural biology. A direct link between protein dynamics and function has been provided by x-ray studies performed on ribonuclease A (RNase A) (B. F. Rasmussen et al., Nature, 1992, Vol. 357, pp. 423-424; L. Vitagliano et al., Proteins: Structure, Function, and Genetics, 2002, Vol. 46, pp. 97-104). Here we report a 3 ns molecular dynamics simulation of RNase A in water aimed at characterizing the dynamical behavior of the enzyme. The analysis of local and global motions provides interesting insight on the dynamics/function relationship of RNase A. In agreement with previous crystallographic reports, the present study confirms that the RNase A active site is constituted by rigid (His12, Asn44, Thr45) and flexible (Lys41, Asp83, His119, Asp121) residues. The analysis of the global motions, performed using essential dynamics, shows that the two beta-sheet regions of RNase A move coherently in opposite directions, thus modifying solvent accessibility of the active site, and that the mixed alpha/3(10)-helix (residues 50-60) behaves as a mechanical hinge during the breathing motion of the protein. These data demonstrate that this motion, essential for RNase A substrate binding and release, is an intrinsic dynamical property of the ligand-free enzyme.     

Molecular dynamics studies of alanine racemase: a structural model for drug design

Alanine racemase (AlaR) is a bacterial enzyme that catalyzes the interconversion of L- and D-alanine, which is an essential constituent of the peptidoglycan layer of the bacterial cell wall and requires pyridoxal 5'-phosphate (PLP) as a cofactor. The enzyme is universal to bacteria, including mycobacteria, making it an attractive target for drug design. To investigate the effects of flexibility on the binding modes of the substrate and an inhibitor and to analyze how the active site is affected by the presence of the substrate versus inhibitor, a molecular dynamics simulation on the full AlaR dimer from Bacillus stearothermophilus (pdb code: 1SFT) with a D-alanine molecule in one active site and the noncovalent inhibitor, propionate, in the second site has been carried out. Within the time scale of the simulation, we show that the active site becomes more stabilized in the presence of substrate versus inhibitor. The results of this simulation are in agreement with the proposed mechanism of alanine racemase reaction in which the substrate carboxyl group directly participates in the catalysis by acting cooperatively with Tyr 265' and Lys 39. A structural water molecule in contact with both substrate and inhibitor (i.e., in both active sites) and bridging residues in both active sites was identified. It shows a remarkably low mobility and does not exchange with bulk water. This water molecule can be taken into account for the design of specific AlaR inhibitors by either utilizing it as a bridging group or displacing it with an inhibitor atom. The results presented here provide insights into the dynamics of the alanine racemase in the presence of substrate/inhibitor, which will be used for the rational design of novel inhibitors.     

Structure and dynamics of Candida rugosa lipase: the role of organic solvent

The effect of organic solvent on the structure and dynamics of proteins was investigated by multiple molecular dynamics simulations (1 ns each) of Candida rugosa lipase in water and in carbon tetrachloride. The choice of solvent had only a minor structural effect. For both solvents the open and the closed conformation of the lipase were near to their experimental X-ray structures (C rms deviation 1–1.3 Å). However, the solvents had a highly specific effect on the flexibility of solvent-exposed side chains: polar side chains were more flexible in water, but less flexible in organic solvent. In contrast, hydrophobic residues were more flexible in organic solvent, but less flexible in water. As a major effect solvent changed the dynamics of the lid, a mobile element involved in activation of the lipase, which fluctuated as a rigid body about its average position. While in water the deviations were about 1.6 Å, organic solvent reduced flexibility to 0.9 Å. This increase rigidity was caused by two salt bridges (Lys85–Asp284, Lys75–Asp79) and a stable hydrogen bond (Lys75–Asn 292) in organic solvent. Thus, organic solvents stabilize the lid but render the side chains in the hydrophobic substrate-binding site more mobile. Figure Superimposition of open (black, PDB entry 1CRL) and closed (gray, PDB entry 1TRH) conformers of C. rugosa lipase. The mobile lid is indicatedThis revised version was published online in October 2004 with corrections to the Graphical Abstract.     

RNA kink-turns as molecular elbows: hydration, cation binding, and large-scale dynamics

The presence of Kink-turns (Kt) at key functional sites in the ribosome (e.g., A-site finger and L7/L12 stalk) suggests that some Kink-turns can confer flexibility on RNA protuberances that regulate the traversal of tRNAs during translocation. Explicit solvent molecular dynamics demonstrates that Kink-turns can act as flexible molecular elbows. Kink-turns are associated with a unique network of long-residency static and dynamical hydration sites that is intimately involved in modulating their conformational dynamics. An implicit solvent conformational search confirms the flexibility of Kink-turns around their X-ray geometries and identifies a second low-energy region with open structures that could correspond to Kink-turn geometries seen in solution experiments. An extended simulation of Kt-42 with the factor binding site (helices 43 and 44) shows that the local Kt-42 elbow-like motion fully propagates beyond the Kink-turn, and that there is no other comparably flexible site in this rRNA region. Kink-turns could mediate large-scale adjustments of distant RNA segments.     

Local fluctuations vs. global unfolding of proteins investigated by limited proteolysis

The structure of a protein molecule consists of both rigid and flexible sections to satisfy the demands for stability and catalysis. Because the flexibility of a protein segment is indispensable for a proteolytic attack, limited proteolysis is a superb tool to analyse both confined local fluctuations and global unfolding events in proteins. While the identification of the primary cleavage products allows the assignment of the flexible regions to the primary structure, the kinetics of proteolytic degradation enables differentiation between local fluctuations in the native protein molecule and the global unfolding process during denaturation. Modifications of the amino acid sequence in the concerned regions can tune proteolytic susceptibility and alter protein stability. In the present paper, we summarise our results on native-state and unfolded-state proteolysis of ribonuclease A (RNase A) and the effect of mutations in the detected flexible regions on the stability and unfolding of the RNase A molecule.     

Local fluctuations vs. global unfolding of proteins investigated by limited proteolysis

The structure of a protein molecule consists of both rigid and flexible sections to satisfy the demands for stability and catalysis. Because the flexibility of a protein segment is indispensable for a proteolytic attack, limited proteolysis is a superb tool to analyse both confined local fluctuations and global unfolding events in proteins. While the identification of the primary cleavage products allows the assignment of the flexible regions to the primary structure, the kinetics of proteolytic degradation enables differentiation between local fluctuations in the native protein molecule and the global unfolding process during denaturation. Modifications of the amino acid sequence in the concerned regions can tune proteolytic susceptibility and alter protein stability. In the present paper, we summarise our results on native-state and unfolded-state proteolysis of ribonuclease A (RNase A) and the effect of mutations in the detected flexible regions on the stability and unfolding of the RNase A molecule.     

Nonspecific yet decisive: Ubiquitination can affect the native‐state dynamics of the modified protein

Ubiquitination is one of the most common post‐translational modifications of proteins, and mediates regulated protein degradation among other cellular processes. A fundamental question regarding the mechanism of protein ubiquitination is whether and how ubiquitin affects the biophysical nature of the modified protein. For some systems, it was shown that the position of ubiquitin within the attachment site is quite flexible and ubiquitin does not specifically interact with its substrate. Nevertheless, it was revealed that polyubiquitination can decrease the thermal stability of the modified protein in a site‐specific manner because of alterations of the thermodynamic properties of the folded and unfolded states. In this study, we used detailed atomistic simulations to focus on the molecular effects of ubiquitination on the native structure of the modified protein. As a model, we used Ubc7, which is an E2 enzyme whose in vivo ubiquitination process is well characterized and known to lead to degradation. We found that, despite the lack of specific direct interactions between the ubiquitin moiety and Ubc7, ubiquitination decreases the conformational flexibility of certain regions of the substrate Ubc7 protein, which reduces its entropy and thus destabilizes it. The strongest destabilizing effect was observed for systems in which Lys48‐linked tetra‐ubiquitin was attached to sites used for in vivo degradation. These results reveal how changes in the configurational entropy of the folded state may modulate the stability of the protein's native state. Overall, our results imply that ubiquitination can modify the biophysical properties of the attached protein in the folded state and that, in some proteins, different ubiquitination sites will lead to different biophysical outcomes. We propose that this destabilizing effect of polyubiquitin on the substrate is linked to the functions carried out by the modification, and in particular, regulatory control of protein half‐life through proteasomal degradation.     

Endosomal proteolysis of insulin-like growth factor-I at its C-terminal D-domain by cathepsin B

IGF-I is degraded within the endosomal apparatus as a consequence of receptor-mediated endocytosis. However, the nature of the responsible protease and the position of the cleavage sites in the IGF-I molecule remain undefined. In vitro proteolysis of IGF-I using an endosomal lysate required an acidic pH and was sensitive to CA074, an inhibitor of the cathepsin B enzyme. By nondenaturing immunoprecipitation, the acidic IGF-I-degrading activity was attributed to the luminal species of endosomal cathepsin B with apparent molecular masses of 32- and 28-kDa. The cathepsin B precursor, procathepsin B, was processed in vitro within isolated endosomes at pH 5 or at 7 in the presence of ATP, the substrate of the vacuolar H(+)-ATPase. The rate of IGF-I hydrolysis using an endosomal lysate or pure cathepsin B was found to be optimal at pH 5-6 and moderate at pH 4 and 7. Competition studies revealed that EGF and IGF-I share a common binding site on the cathepsin B enzyme, with native IGF-I displaying the lowest affinity for the protease (IC50 approximately 1.5 microM). Hydrolysates of IGF-I generated at low pH by endosomal IGF-I-degrading activity and analyzed by reverse-phase HPLC and mass spectrometry revealed cleavage sites at Lys68-Ser69, Ala67-Lys68, Pro66-Ala67 and Lys65-Pro66 within the C-terminal D-domain of IGF-I. Treatment of human HepG2 hepatoma cells with the cathepsin B proinhibitor CA074-Me reduced, in vivo, the intracellular degradation of internalized [125I]IGF-I and, in vitro, the degradation of exogenous [125I]IGF-I incubated with the cell-lysates at pH 5. Inhibitors of cathepsin B and pro-cathepsin B processing, which abolish endosomal proteolysis of IGF-I and alter tumor cell growth and IGF-I receptor signalling, merit investigation as antimetastatic drugs.     

Mass spectrometric identification of the trypsin cleavage pathway in lysyl-proline containing oligotuftsin peptides.

Trypsin cleaves specifically peptide bonds at the C-terminal side of lysine and arginine residues, except for -Arg-Pro- and -Lys-Pro- bonds which are normally resistant to proteolysis. Here we report evidence for a -Lys-Pro- tryptic cleavage in modified oligotuftsin derivatives, Ac-[TKPKG]4-NH2) (1), using high-resolution mass spectrometry and HPLC as primary methods for analysis of proteolytic reactions. The proteolytic susceptibility of -Lys-Pro- bonds was strongly dependent on flanking residues, and the flexibility of the peptide backbone might be a prerequisite for this unusual cleavage. While -Lys-Gly- bonds in 1 were rapidly cleaved, the modification of these Lys residues by the attachment of a ss-amyloid(4-10) epitope to yield -Lys(X)-Gly derivatives prevented cleavage of this bond, and provided trypsin cleavage of -Lys-Pro- bonds, the pathway of this degradation being independent on the type of Lys-N(epsilon)-side chains (acetyl group, amino acid, peptide). Substitution of the Lys residues by Ala at the P'2 positions decreased the tryptic cleavage, while replacement of the bulky side chain of Thr at the P2 positions strongly increased the cleavage of -Lys-Pro- bonds. Circular dichroism (CD) data of the modified oligotuftsin derivatives are in accord with enhanced flexibility of the peptide backbone, as a prerequisite for increased susceptibility to cleavage of -Lys-Pro- bonds. These results obtained of oligotuftsin derivatives might have implications for the proteolytic degradation of target peptides that require specific conformational preconditions.     

A thermally sensitive loop in clostridial glutamate dehydrogenase detected by limited proteolysis

The structural flexibility and thermostability of glutamate dehydrogenase (GDH) from Clostridium symbiosum were examined by limited proteolysis using three proteinases with different specificities, trypsin, chymotrypsin, and endoproteinase Glu-C. Clostridial GDH resisted proteolysis by any of these enzymes at 25 degrees C. Above 30 degrees C, however, GDH became cleavable by chymotrypsin, apparently at a single site. SDS-PAGE indicated the formation of one large fragment with a molecular mass of approximately 44 kDa and one small one of <10 kDa. Proteolysis was accompanied by the loss of enzyme activity, which outran peptide cleavage, suggesting a cooperative conformational change. Proteolysis was prevented by either of the substrates 2-oxoglutarate or l-glutamate but not by the coenzymes NAD(+) or NADH. Circular dichroism spectroscopy indicated that the protective effects of these ligands resulted from fixation of flexible regions of the native structure of the enzyme. Size-exclusion chromatography and SDS-PAGE studies of chymotrypsin-treated GDH showed that the enzyme retained its hexameric structure and all of its proteolytic fragments. However, circular dichroism spectroscopy and analytical ultracentrifugation showed global conformational changes affecting the overall compactness of the protein structure. Chymotrypsin-catalyzed cleavage also diminished the thermostability of GDH and the cooperativity of the transition between its native and denatured states. N-terminal amino acid sequencing and mass spectrometry showed that heat-induced sensitivity to chymotrypsin emerged in the loop formed by residues 390-393 that lies between helices alpha(15) and alpha(16) in the folded structure of the enzyme.     

On the effect of mutations in bovine or camel chymosin on the thermodynamics of binding κ‐caseins

Bovine and camel chymosins are aspartic proteases that are used in dairy food manufacturing. Both enzymes catalyze proteolysis of a milk protein, κ‐casein, which helps to initiate milk coagulation. Surprisingly, camel chymosin shows a 70% higher clotting activity than bovine chymosin for bovine milk, while exhibiting only 20% of the unspecific proteolytic activity. By contrast, bovine chymosin is a poor coagulant for camel milk. Although both enzymes are marketed commercially, the disparity in their catalytic activity is not yet well understood at a molecular level, due in part to a lack of atomistic resolution data about the chymosin—κ‐casein complexes. Here, we report computational alanine scanning calculations of all four chymosin—κ‐casein complexes, allowing us to elucidate the influence that individual residues have on binding thermodynamics. Of the 12 sequence differences in the binding sites of bovine and camel chymosin, eight are shown to be particularly important for understanding differences in the binding thermodynamics (Asp112Glu, Lys221Val, Gln242Arg, Gln278Lys. Glu290Asp, His292Asn, Gln294Glu, and Lys295Leu. Residue in bovine chymosin written first). The relative binding free energies of single‐point mutants of chymosin are calculated using the molecular mechanics three dimensional reference interaction site model (MM‐3DRISM). Visualization of the solvent density functions calculated by 3DRISM reveals the difference in solvation of the binding sites of chymosin mutants.     

Characterization of guanosine diphospho-D-mannose dehydrogenase from Pseudomonas aeruginosa. Structural analysis by limited proteolysis.

Alginate is believed to be a major virulence factor in the pathogenicity of Pseudomonas aeruginosa in the lungs of patients suffering from cystic fibrosis. Guanosine diphospho-D-mannose dehydrogenase (GDPmannose dehydrogenase, EC 1.1.1.132) is a key enzyme in the alginate biosynthetic pathway which catalyzes the oxidation of guanosine diphospho-D-mannose (GDP-D-mannose) to GDP-D-mannuronic acid. In this paper, we report the structural analysis of GMD by limited proteolysis using three different proteases, trypsin, submaxillary Arg-C protease, and chymotrypsin. Treatment of GMD with these proteases indicated that the amino-terminal part of this enzyme may fold into a structural domain with an apparent molecular mass of 25-26 kDa. Multiple proteolytic cleavage sites existed at the carboxyl-terminal end of this domain, indicating that this segment may represent an exposed region of the protein. Initial proteolysis also generated a carboxyl-terminal fragment with an apparent molecular mass of 16-17 kDa which was further digested into smaller fragments by trypsin and chymotrypsin. The proteolytic cleavage sites were localized by partial amino-terminal sequencing of the peptide fragments. Arg-295 was identified as the initial cleavage site for trypsin and Tyr-278 for chymotrypsin. Catalytic activity of GMD was totally abolished by the initial cleavage. However, binding of the substrate, GDP-D-mannose, increased stability toward proteolysis and inhibited the loss of enzyme activity. GMP and GDP (guanosine 5'-mono- and diphosphates) also blocked the initial cleavage, but NAD and mannose showed no effect. These results suggest that binding of the guanosine moiety at the catalytic site of GMD may induce a conformational change that reduces the accessibility of the cleavage sites to proteases. Binding of [14C]GDP-D-mannose to the amino-terminal domain was not affected by the removal of the carboxyl-terminal 16-kDa fragment. Furthermore, photoaffinity labeling of GMD with [32P]arylazido-beta-alanine-NAD followed by proteolysis demonstrated that the radioactive NAD was covalently linked to the amino-terminal domain. These observations imply that the amino-terminal domain (25-26 kDa) contains both the substrate and cofactor binding sites. However, the carboxyl-terminal fragment (16-17 kDa) may possess amino acid residues essential for catalysis. Thus, proteolysis had little effect on substrate binding, but totally eliminated catalysis. These biochemical data are in complete agreement with amino acid sequence analysis for the existence of substrate and cofactor sites of GMD. A linear peptide map of GMD was constructed for future structure/functional studies.     

Comparison of SARS and NL63 papain-like protease binding sites and binding site dynamics: inhibitor design implications

The human severe acute respiratory syndrome coronavirus (SARS-CoV) and the NL63 coronaviruses are human respiratory pathogens for which no effective antiviral treatment exists. The papain-like cysteine proteases encoded by the coronavirus (SARS-CoV: PLpro; NL63: PLP1 and PLP2) represent potential targets for antiviral drug development. Three recent inhibitor-bound PLpro structures highlight the role of an extremely flexible six-residue loop in inhibitor binding. The high binding site plasticity is a major challenge in computational drug discovery/design efforts. From conventional molecular dynamics and accelerated molecular dynamics (aMD) simulations, we find that with conventional molecular dynamics simulation, PLpro translationally samples the open and closed conformation of BL2 loop on a picosecond-nanosecond timescale but does not reproduce the peptide bond inversion between loop residues Tyr269 and Gln270 that is observed on inhibitor GRL0617 binding. Only aMD simulation, starting from the closed loop conformation, reproduced the 180° ?-ψ dihedral rotation back to the open loop state. The Tyr-Gln peptide bond inversion appears to involve a progressive conformational change of the full loop, starting at one side, and progressing to the other. We used the SARS-CoV apo X-ray structure to develop a model of the NL63-PLP2 catalytic site. Superimposition of the PLP2 model on the PLpro X-ray structure identifies binding site residues in PLP2 that contribute to the distinct substrate cleavage site specificities between the two proteases. The topological and electrostatic differences between the two protease binding sites also help explain the selectivity of non-covalent PLpro inhibitors.     

Effect of oxidative sulfitolysis of disulfide bonds of bovine serum albumin on its structural properties: A physicochemical study

The effect of S-S bond cleavage of bovine serum albumin (BSA) on some of its structural properties, including solubility, viscosity, and conformation, were investigated. Cleavage of S-S bonds decreased the solubility of serum albumin and also shifted its isoelectric point to lowerpH values. S-S bond cleavage resulted in changes in shape and hydrodynamic volume of the protein, increasing the specific viscosity, with cleavage of up to 14 S-S bonds, followed by a decrease with further cleavage. Both UV difference and fluorescence spectral measurements indicated that conformational flexibility increases with S-S bond cleavage. Secondary structure estimations by far UV-CD suggested a gradual decrease in -helical content of the protein with progressive cleavage of its S-S bonds. However, fully S-S bond cleaved protein maintained some -helical structure. Sulfitolysis of the protein also decreased its I, 8-anilino-naphthalene sulfonate-binding ability.     

Increased proteolytic resistance of ribonuclease A by protein engineering.

Although highly stable toward unfolding, native ribonuclease A is known to be cleaved by unspecific proteases in the flexible loop region near Ala20. With the aim to create a protease-resistant ribonuclease A, Ala20 was substituted for Pro by site-directed mutagenesis. The resulting mutant enzyme was nearly identical to the wild-type enzyme in the near-UV and far-UV circular dichroism spectra, in its activity to 2',3'-cCMP and in its thermodynamic stability. However, the proteolytic resistance to proteinase K and subtilisin Carlsberg was extremely increased. Pseudo-first-order rate constants of proteolysis, determined by densitometric analysis of the bands of intact protein in SDS-PAGE, decreased by two orders of magnitude. In contrast, the rate constant of proteolysis with elastase was similar to that of the wild-type enzyme. These differences can be explained by the analysis of the fragments occurring in proteolysis with elastase. Ser21-Ser22 was identified as the main primary cleavage site in the degradation of the mutant enzyme by elastase. Obviously, this bond is not cleavable by proteinase K or subtilisin Carlsberg. The results demonstrate the high potential of a single mutation in protein stabilization to proteolytic degradation.     

Limited proteolysis of Escherichia coli cytidine 5'-triphosphate synthase. Identification of residues required for CTP formation and GTP-dependent activation of glutamine hydrolysis.

Cytidine 5'-triphosphate synthase catalyses the ATP-dependent formation of CTP from UTP using either ammonia or l-glutamine as the source of nitrogen. When glutamine is the substrate, GTP is required as an allosteric effector to promote catalysis. Limited trypsin-catalysed proteolysis, Edman degradation, and site-directed mutagenesis were used to identify peptide bonds C-terminal to three basic residues (Lys187, Arg429, and Lys432) of Escherichia coli CTP synthase that were highly susceptible to proteolysis. Lys187 is located at the CTP/UTP-binding site within the synthase domain, and cleavage at this site destroyed all synthase activity. Nucleotides protected the enzyme against proteolysis at Lys187 (CTP > ATP > UTP > GTP). The K187A mutant was resistant to proteolysis at this site, could not catalyse CTP formation, and exhibited low glutaminase activity that was enhanced slightly by GTP. K187A was able to form tetramers in the presence of UTP and ATP. Arg429 and Lys432 appear to reside in an exposed loop in the glutamine amide transfer (GAT) domain. Trypsin-catalyzed proteolysis occurred at Arg429 and Lys432 with a ratio of 2.6 : 1, and nucleotides did not protect these sites from cleavage. The R429A and R429A/K432A mutants exhibited reduced rates of trypsin-catalyzed proteolysis in the GAT domain and wild-type ability to catalyse NH3-dependent CTP formation. For these mutants, the values of kcat/Km and kcat for glutamine-dependent CTP formation were reduced approximately 20-fold and approximately 10-fold, respectively, relative to wild-type enzyme; however, the value of Km for glutamine was not significantly altered. Activation of the glutaminase activity of R429A by GTP was reduced 6-fold at saturating concentrations of GTP and the GTP binding affinity was reduced 10-fold. This suggests that Arg429 plays a role in both GTP-dependent activation and GTP binding.     

Molecular dynamics study of a hyperthermophilic and a mesophilic rubredoxin

In recent years, increased interest in the origin of protein thermal stability has gained attention both for its possible role in understanding the forces governing the folding of a protein and for the design of new highly stable engineered biocatalysts. To study the origin of thermostability, we have performed molecular dynamics simulations of two rubredoxins, from the mesophile Clostridium pasteurianum and from the hyperthermophile Pyrococcus furiosus. The simulations were carried out at two temperatures, 300 and 373 K, for each molecule. The length of the simulations was within the range of 6-7.2 ns. The rubredoxin from the hyperthermophilic organism was more flexible than its mesophilic counterpart at both temperatures; however, the overall flexibility of both molecules at their optimal growth temperature was the same, despite 59% sequence homology. The conformational space sampled by both molecules was larger at 300 K than at 373 K. The essential dynamics analysis showed that the principal overall motions of the two molecules are significantly different. On the contrary, each molecule showed similar directions of motion at both temperatures.     

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.