Evidence for a high pKa of an aspartic acid residue in the active site of CALB by a fully atomistic multiscale approach.

Candida antarctica Lipase B (CALB) is a paradigm for the family of lipases. At pH 7, the optimal pH for catalysis, the protonation state of an aspartic acid of the active site (Asp134) could not be conclusively assigned. In fact, the pKa estimate provided by a widely used computational tool, namely PropKa, that predicts pKa values of ionizable groups in proteins based on the crystallographic structure, is only slightly above 7 (pKa = 7.25). This, along with the lack of an experimental evaluation, makes the assignment of its protonation state at neutral pH challenging. Here, we calculate the pKa of Asp134 by means of a fully atomistic multiscale computational approach based on classical molecular dynamics (MD) simulation and the perturbed matrix method (PMM), namely the MD-PMM approach. MD-PMM is able to take into account the dynamics of the system and, at the same time, to treat the deprotonation step at the quantum level. The calculations provide a pKa value of 8.9 ± 1.1, hence suggesting that Asp134 in CALB should be protonated at neutral, and even at slightly basic, pH.Communicated by Ramaswamy H. Sarma.

Reducing Properties of Superalkalis on Pyridinic Graphene Surfaces: a Computational Study.

The hyperlithiated species Li k + 1 F k (k=1, 2, 3, and 4) have been studied by quantum mechanical (QM) methods. Different structures have been localized for each molecule by the CBS-QB3 composite method: all the isomers show superalkali properties and strong tendency to donate an electron to carbon dioxide forming stable Li k + 1 F k · · · CO 2 complexes. With the aim to find molecular systems able to stabilize superalkalis, geometries of complexes between superalkalis and pyridine and superalkalis and graphene surfaces doped with a pyridinic vacancy were calculated. The pyridinic graphene sheets were modeled with two finite molecular systems C69 H21 N3 and C117 H27 N3 . The interaction with one pyridine molecule is quite weak and the superalkali maintains its structure and electron properties. The affinity for graphene sheets is instead stronger and the superalkalis tend to deform their geometry to better interact with the graphene surface. However, the superalkalis continue to show the tendency to transfer electrons to carbon dioxide reducing CO2 , as found in graphene absence.