A simple clustering algorithm can be accurate enough for use in calculations of pKs in macromolecules

Flavodoxin: a compromise between efficiency and versatility in the electron transfer from Photosystem I to Ferredoxin-NADP(+) reductase

Under iron-deficient conditions Flavodoxin (Fld) replaces Ferredoxin in Anabaena as electron carrier from Photosystem I (PSI) to Ferredoxin-NADP(+) reductase (FNR). Several residues modulate the Fld interaction with FNR and PSI, but no one appears as specifically critical for efficient electron transfer (ET). Fld shows a strong dipole moment, with its negative end directed towards the flavin ring. The role of this dipole moment in the processes of interaction and ET with positively charged surfaces exhibited by PSI and FNR has been analysed by introducing single and multiple charge reversal mutations on the Fld surface. Our data confirm that in this system interactions do not rely on a precise complementary surface of the reacting molecules. In fact, they indicate that the initial orientation driven by the alignment of dipole moment of the Fld molecule with that of the partner contributes to the formation of a bunch of alternative binding modes competent for the efficient ET reaction. Additionally, the fact that Fld uses different interaction surfaces to dock to PSI and to FNR is confirmed.

Analysis of basic clustering algorithms for numerical estimation of statistical averages in biomolecules

In statistical mechanics, the equilibrium properties of a physical system of particles can be calculated as the statistical average over accessible microstates of the system. In general, these calculations are computationally intractable since they involve summations over an exponentially large number of microstates. Clustering algorithms are one of the methods used to numerically approximate these sums. The most basic clustering algorithms first sub-divide the system into a set of smaller subsets (clusters). Then, interactions between particles within each cluster are treated exactly, while all interactions between different clusters are ignored. These smaller clusters have far fewer microstates, making the summation over these microstates, tractable. These algorithms have been previously used for biomolecular computations, but remain relatively unexplored in this context. Presented here, is a theoretical analysis of the error and computational complexity for the two most basic clustering algorithms that were previously applied in the context of biomolecular electrostatics. We derive a tight, computationally inexpensive, error bound for the equilibrium state of a particle computed via these clustering algorithms. For some practical applications, it is the root mean square error, which can be significantly lower than the error bound, that may be more important. We how that there is a strong empirical relationship between error bound and root mean square error, suggesting that the error bound could be used as a computationally inexpensive metric for predicting the accuracy of clustering algorithms for practical applications. An example of error analysis for such an application-computation of average charge of ionizable amino-acids in proteins-is given, demonstrating that the clustering algorithm can be accurate enough for practical purposes.

Evaluation of the ‘side door’ in carboxylesterase-mediated catalysis and inhibition

Structures of mammalian carboxylesterases (CEs) reveal the presence of a ‘side door’ that is proposed to act as an alternative pore for the trafficking of substrates and products. p-Nitrobenzyl esterase (pnb CE) from Bacillus subtilis exhibits close structural homology and a similar side-door domain as mammalian CEs. We investigated the role of a specific ‘gate’ residue at the side door (i.e., Leu 362) during pnb CE-catalyzed hydrolysis of model esters, pesticides, and lipids. Recombinant pnb CE proteins containing mutations at position 362 demonstrated markedly lower k cat and k cat/K m values. The mutation with the most significant impact on catalysis was the L362R mutant (k cat/K m was 22-fold lower). Moreover, the ability of the L362R mutant to be inhibited by organophosphates (OP) was also lower. Investigation into the altered catalytic proficiency using pH-activity studies indicated that the catalytic triad of the mutant enzyme was preserved. Furthermore, viscosity variation and carbamate inhibition experiments indicated that rates of substrate association and acylation/deacylation were lower. Finally, recombinant CEs were found to possess lipolytic activity toward cholesteryl oleate and 2-arachidonylglycerol. In summary, the L362R mutant CE markedly slowed the rate of ester hydrolysis and was less sensitive to OP inhibition. The apparent causes of the diminished catalysis are discussed.