Molecular dynamics simulations of human L-asparaginase1: Insights into structural determinants of enzymatic activity

The L-asparaginase enzyme is used in cancer therapy, mainly acute lymphoid leukemia (ALL). Commercial enzymes (EcASNase2) cause adverse reactions during treatment, such as immunogenicity. A human enzyme could be a non-immunogenic substitute. However, no candidate was found showing efficient kinetic properties. HASNase1 is an L-asparaginase that comes from the N-terminal domain of a protein called 60 kDa-lysophospholipase and its 3D structure has not been resolved. HASNase1 is homologous to EcASNase1 and gpASNase1, and this last one has shown efficient kinetic properties. Homology modeling was used to find the 3D structure of hASNase1, so one could submit it to Molecular Dynamics (MD), in order to understand structural differences that lead to different catalytic efficiency compared to EcASNase2 and gpASNase1. The interaction potential between L-Asn and active site residues showed that the substrate can rotate in the site when Region1 is open. Region1 residues sequence favors deformations and movements as shown in MD. Region2-A is linear in gpASNase1, and it features a helix portion in hASNase1, which leaves the Tyr308 position projected to the active site ratifying its role in catalytic efficiency. Analysis of Lys188 orientation and movement showed the effect of positive cooperativity in hASNase1. It was found that the presence of Asn at the allosteric site helps, not only in Region1 stabilization, but also in Lys188 stabilization for the maintenance of the triad. Despite structural similarities in hASNase1, gpASNase1, and EcASNase2, there are differences in structural determinants that, in addition to allosterism, may explain the different kinetic properties.