The Bacillus subtilis endospore coat protein CotA shows laccase activity. By using comparative modeling techniques, we were able to derive a model for CotA based on the known x-ray structures of zucchini ascorbate oxidase and Cuprinus cereneus laccase. This model of CotA contains all the structural features of a laccase, including the reactive surface-exposed copper center (T1) and two buried copper centers (T2 and T3). Single amino acid substitutions in the CotA T1 copper center (H497A, or M502L) did not prevent assembly of the mutant proteins into the coat and did not alter the pattern of extractable coat polypeptides. However, in contrast to a wild type strain, both mutants produced unpigmented colonies and spores unable to oxidize syringaldazine (SGZ) and 2'2-azino-bis-(3-ethylbenzothiazoline-6-sulfonic acid) (ABTS). The CotA protein was purified to homogeneity from an overproducing Escherichia coli strain. The purified CotA shows an absorbance and a EPR spectra typical of blue multicopper oxidases. Optimal enzymatic activity was found at < or =pH 3.0 and at pH 7.0 for ABTS or SGZ oxidation, respectively. The apparent K(m) values for ABTS and SGZ at 37 degrees C were of 106 +/- 11 and 26 +/- 2 microm, respectively, with corresponding k(cat) values of 16.8 +/- 0.8 and 3.7 +/- 0.1 s(-1). Maximal enzyme activity was observed at 75 degrees C with ABTS as substrate. Remarkably, the coat-associated or the purified enzyme showed a half-life of inactivation at 80 degrees C of about 4 and 2 h, respectively, indicating that CotA is intrinsically highly thermostable.

A new method is proposed for performing constant-pH molecular dynamics (MD) simulations, that is, MD simulations where pH is one of the external thermodynamic parameters, like the temperature or the pressure. The protonation state of each titrable site in the solute is allowed to change during a molecular mechanics (MM) MD simulation, the new states being obtained from a combination of continuum electrostatics (CE) calculations and Monte Carlo (MC) simulation of protonation equilibrium. The coupling between the MM/MD and CE/MC algorithms is done in a way that ensures a proper Markov chain, sampling from the intended semigrand canonical distribution. This stochastic titration method is applied to succinic acid, aimed at illustrating the method and examining the choice of its adjustable parameters. The complete titration of succinic acid, using constant-pH MD simulations at different pH values, gives a clear picture of the coupling between the trans/gauche isomerization and the protonation process, making it possible to reconcile some apparently contradictory results of previous studies. The present constant-pH MD method is shown to require a moderate increase of computational cost when compared to the usual MD method.