“…The starting configurations of the enzyme were based on the 4PEP crystallographic structure of porcine pepsin at 1.8 Å resolution available from the protein data bank [21]. The total electrostatic charge of pepsin, −6|e|, was determined previously by calculating the pKa of the protein side chains at pH = 3.6 [26], where the protein loading in SBA-15 is maximal. In particular, His53, Lys319, Arg307 and Arg315 were modelled as positively charged residues, whereas 10 Asp residues (11,87,118,138,171,195,215,242,290,314) were unprotonated (negatively charged).…”