Mechanisms underlying dioxygen reduction in laccases. Structural and modelling studies focusing on proton transfer
© Bento et al; licensee BioMed Central Ltd. 2010
Received: 26 November 2009
Accepted: 7 September 2010
Published: 7 September 2010
Laccases are enzymes that couple the oxidation of substrates with the reduction of dioxygen to water. They are the simplest members of the multi-copper oxidases and contain at least two types of copper centres; a mononuclear T1 and a trinuclear that includes two T3 and one T2 copper ions. Substrate oxidation takes place at the mononuclear centre whereas reduction of oxygen to water occurs at the trinuclear centre.
In this study, the CotA laccase from Bacillus subtilis was used as a model to understand the mechanisms taking place at the molecular level, with a focus in the trinuclear centre. The structures of the holo-protein and of the oxidised form of the apo-protein, which has previously been reconstituted in vitro with Cu(I), have been determined. The former has a dioxygen moiety between the T3 coppers, while the latter has a monoatomic oxygen, here interpreted as a hydroxyl ion. The UV/visible spectra of these two forms have been analysed in the crystals and compared with the data obtained in solution. Theoretical calculations on these and other structures of CotA were used to identify groups that may be responsible for channelling the protons that are needed for reduction of dioxygen to water.
These results present evidence that Glu 498 is the only proton-active group in the vicinity of the trinuclear centre. This strongly suggests that this residue may be responsible for channelling the protons needed for the reduction. These results are compared with other data available for these enzymes, highlighting similarities and differences within laccases and multicopper oxidases.
Multicopper oxidases (MCOs) are a group of enzymes that are able to couple the oxidation of a variety of different substrates concomitantly with dioxygen reduction to water [1–3]. This protein family comprises laccases, (p-diphenol: dioxygen oxidoreductase, EC 22.214.171.124), metalloxidases, ascorbate oxidase and ceruloplasmin [2, 4, 5]. They are found in prokaryotes, eukaryotes as well as in archea. In plants, they have been associated with lignin formation; in fungi, with pigment formation, lignin degradation and detoxification processes; in yeast and mammals with iron metabolism; in bacteria with copper homeostasis and manganese oxidation ([6, 7] and references therein). Overall, they are capable of oxidise substrates that vary from organic compounds, such as ascorbic acid, phenolate siderophores and organic amines, to inorganic ones, such as metal ions like Fe(II), Cu(I) and Mn(II) ([6–8]and references therein). In addition, they are able to transfer their electrons to molecular oxygen with the concomitant production of water, a mechanism that can also constitute a key factor in the management of oxygen in aerobic organisms . These characteristics made them important targets for structure-function studies. Furthermore, laccases have been show to have diverse biotechnological applications constituting a model for this type of studies .
MCOs catalyze one-electron oxidation processes, and four molecules of substrate are oxidised in order to reduce a dioxygen molecule to two waters molecules. Substrate oxidation occurs at the mononuclear T1 centre and then the electrons are shuttled, along a T1 coordinating cisteine, to the two histidines that are coordinating the T3 coppers of the trinuclear centre (Figure 1b), where reduction of dioxygen occurs [2, 4, 11]. This constitutes a HCH conserved motif characteristically found in MCOs. Four electrons, as well as four protons, are used to reduce a molecule of dioxygen with the concomitant formation of two water molecules [1, 2, 12]. Electrons needed for this process are obtained through the oxidation of a variety of substrates, but not much is known about the mechanism of proton transfer during this process. Recently, site directed mutagenesis studies suggest that Asp112 of CueO (Asp 116 in CotA; Figure 1c) [13, 14], a conserved residue that is located in the exit channel in close proximity to the T2 copper ion, plays an important role in the protonation process. A similar role has been suggested for Asp94 in Fet3p [15–17], which is also the equivalent residue of CotA Asp116. This internalised acidic group participates in a network of hydrogen bonds involving ligands of the trinuclear centre (see figure 1c), being connected with Cu4 through a chain of water molecules. Site directed mutagenesis studies, on the prokaryotic multicopper oxidase CueO  and on the eukaryotic multicopper oxidase Fet3p [16, 17], provided evidence for the importance of acidic residues located in the entrance channel in the reduction of dioxygen by these enzymes. These are Glu506 in CueO  (equivalent to Glu498 in CotA (see figure 1c)), and Glu487 in Fet3P [16, 17] (located in the other side of the channel). Despite the fact that these two residues are not conserved in sequence in all MCOs, the authors provide evidence for their participation in proton transfer to the reduction of dioxygen. In many structures, this acidic group can form a hydrogen bond with a water molecule that interacts with the oxygen species between the T3 coppers (see figure 1c).
In order to study the mechanism of dioxygen reduction to water by MCOs, CotA, a thermostable and thermoactive laccase, from Bacillus subtilis, was used as a model system . The X-ray structure determination of CotA-laccase, in different redox states and in complex with hydrogen peroxide and azide, together with the crystals structures from other MCOs already deposited in the PDB http://www.pdb.org, allowed the proposal of a mechanism for such an enzymatic process . However, many questions remain to be addressed, including the mechanism of protonation, the nature of the resting state of these enzymes and the identification of key residues that are involved in the reduction of dioxygen. The first MCO X-ray crystal structure to be reported was that of ascorbate oxidase , in which a hydroxyl moiety bridged the two T3 copper ions; this structure has frequently been taken as the template for the native state. Most spectroscopic measurements have also indicated that MCOs, in their resting state, have a hydroxyl group as a bridging ligand. However, nearly all crystal structures of CotA that have been determined [20, 22, 23], and in particular the one where full copper content was achieved by soaking the crystals with CuCl2 (from now on defined as holoCu(II) , PDB code: 1W6L), indicate that the structure of the oxidised state of this enzyme has a dioxygen moiety bridging the two T3 copper ions  (Figures 1b and 1c). It is therefore pertinent to ask whether the CotA-laccase has the same resting state as others MCOs. In order to try to answer this question the X-ray structures of two new preparations of CotA laccase  were determined. The first one corresponds to an in vivo fully copper loaded CotA protein (from now on called holoCotA), isolated from cells grown in copper supplemented media under microaerobic conditions. The second one corresponds to an in vitro reconstituted (with Cu(I)) fully copper loaded protein (from now on called apoCu(I)). This apo protein was obtained from cells grown aerobically in unsupplemented copper media, and despite having been soaked with Cu(I), its spectra shows that it contains Cu(II) , meaning that the protein became oxidised when exposed to the oxygen. The structural comparison between these structures and the previous determined holoCu(II)  will be performed in order to study the conformational effects of the different ways for incorporating copper in CotA. Moreover, the new structures determined here, together with the structures of other states, obtained for the same protein, were used in equilibrium protonation simulations in order to locate the groups likely to be involved in proton transfer. This combined structural biology approach, where X-ray structures of a system are obtained in controlled conditions by the same team (crystallization conditions kept as similar as possible, as well as refinement procedures), and their conformational differences clearly identified in the electron density maps, is fundamental for the success of this type of calculations and mechanistic analysis, which are sensitive to small details. The application of careful procedures (see below in Methods) was tested before with success in the analysis of multihaem cytochromes [25, 26] as well as in CotA mutants . This type of proton transfer studies on laccases may also be helpful for understanding other dioxygen processing enzymes, where proton transfer plays a key role. The best known examples are the haem-copper oxidases , where electron transfer and dioxygen reduction, at the core of the protein, are associated with proton transfer of "chemical" protons, necessary for the reaction, as well as with vectorial proton pumping across the membrane. Similar types of approaches as the ones applied here for CotA laccase, based on continuum electrostatics/Monte Carlo calculations (CE/MC) done on structural data, were used in haem-copper oxidases in order to understand their proton transfer determinants [28–32]. Despite different molecular architectures, the physical principles behind proton transfer mechanisms in multicopper oxidases should be similar, at least in part, with the ones present in haem-copper, especially with what concerns "chemical" protons.
Results and Discussion
Structural characterization of holoCotA and of apoCu(I)
Occupancies of Copper sites.
Refinement and Quality of Refined Models
No. of protein atoms
No. of solvent atoms
No. of hetero atoms
Final free R-factor
Mean B values (Å2)
Estimated overall coordinate uncertainty (Å) ‡
Distance deviations †
Bond distances (Å)
Bond angles (°)
Planar groups (Å)
Chiral volume deviation (Å3)
Ramachandran analysis %
Mechanism of dioxygen reduction to water: insights into the characterization of the resting state of MCOs
Mechanism of dioxygen reduction to water: groups involved in the proton transfer mechanism
The availability of three-dimensional structures for several states of the trinuclear centre, with different oxygen species bound, allows the investigation of the protonation of ionisable groups involved in the enzymatic mechanism. A working hypothesis for the mechanism in this type of enzymes is presented on figure 4, only considering the water molecules that are consistently observed in all deposited three-dimensional structures (this is a revised version of the one presented on ). CE/MC calculations were done on these structures, not taking into account any other water molecules in an explicit way. However, mobile water molecules that should exist inside protein structures are modelled implicitly by this procedure, since it considers empty cavities inside the protein as high dielectric zones containing solvent. Despite the fact that several studies (see for instance [37–40]) document the existence of hydrophobic cavities inside proteins, which are devoid of water and are, therefore, not correctly treated by this procedure, no such cavities exist inside CotA near the two active sites. The closest cavity in CotA, which can potentially display such characteristics, is almost 20 Å away from the dioxygen binding site and more than 22 Å from the T1 copper. Not considering water molecules here, and in other proton transfer processes, is a simplification, and clearly we cannot describe the whole proton transfer process. However, this is a necessary limitation of the present methodology, given the high energetic cost of deprotonating water molecules. For this reason, their deprotonation is transient in nature, and they may act through a Grotthuss like mechanism (see ) for a review), being present in "proton wires". Nevertheless, in long distance proton transfer, one usually finds protonatable groups that can form protonation intermediates in the "proton wire", allowing protons to be transferred in smaller steps. The aim here is to identify these protonatable groups that were likely to be involved in the process of proton transfer, necessary for the reduction of molecular oxygen to water. These groups are probably located close to the trinuclear centre in order to transfer protons to the intermediate oxygen species formed along the mechanism (see figure 4). Analysis of the protonation equilibrium simulations evidenced that the only group actively titrating around pH 7, near the zone of interest, is Glu 498, an acidic group present in the dioxygen entrance channel in CotA laccase. The pH titration of Glu 498 obtained for all the available structurally characterised states of the trinuclear centre of CotA is presented in figure 5 and will be discussed in the following paragraph. In view of the limitations of this type of simulations (which can have root mean square deviations below 1 pK units ), the rationalisation of the results should consider pH ranges, rather than exact values of pH. Additionally, what we aim at analysing is the difference between the different redox and ligand binding states of the protein, which should be less prone to the errors described above; therefore we will concentrate on analysing differential behaviour mostly.
Starting from the fully functional enzyme with dioxygen between the T3 copper ions ((see Figure 4a); holoCotA), the simulations show that the protonation of Glu498 is negligible around pH 7. This situation changes dramatically for the peroxide adduct, which corresponds to a 2 electron reduced species of the previous state (Figure 4b; CotA-H2O2): Glu 498 increases its protonation by a large degree, showing a considerable protonated fraction around pH 7. The next two states in the mechanism in figure 4 (4c and 4d) have not been structurally characterised in CotA, and therefore we could not perform calculations on them. The state in figure 4c is observed if more electrons and protons (possibly coming from Glu 498) are delivered to the trinuclear centre, resulting in the native intermediate that contains two hydroxyl groups; this is observed in the laccase structure from T. hirsuta (Figure 4c) . Upon further structural rearrangements another configuration for these two hydroxyl groups is reached; this is the state observed in the ascorbate oxidase structure (Figure 4d) . The next state in the mechanism could be structurally characterised in CotA and it is a state that results from a release, upon protonation, of one of the hydroxyl groups as water; this is the structure of apoCu(I) (Figure 4e) determined in this work. Calculations using this structure show that the pH titration curve of Glu 498 does not evidence a very marked protonation around pH 7, but a higher affinity for protons is observed instead, when compared with the dioxygen bound structure. Upon further protonation of this hydroxyl moiety (hypothetically by protons coming from Glu 498), its release as a water molecule, and the entrance of another dioxygen molecule, the holoCotA structure is restored, closing the cycle. The fully reduced structure ((Figure 4f); CotA-red) , has Glu 498 predominantly protonated.
One could argue that the results presented here, by having been done with one particular enzyme of the bacterial multicopper oxidase family - CotA laccase - may be specific for this particular enzyme. Namely, one could raise the possibility that the residues equivalent to Glu 498 and Asp 116 in other multicopper oxidases, which have been subjected to site-directly mutagenesis, may behave differently. In order to investigate this, we set up test calculations using the structure of the multicopper oxidase CueO from E. coli. Given that the only structurally characterised state in this multicopper oxidase is the oxidised, hydroxyl bridged T3 state (corresponding to apoCu(I) state of CotA), we had to setup mimics for the two hydroxyl state (another hydroxyl bound to T2) and for the fully reduced state. These studies evidence (results not shown) that Glu 506, similar to its equivalent Glu 498 in CotA, is very sensitive to the redox and liganted states of the trinuclear centre, being mostly deprotonated in the oxidised, hydroxyl bridged T3 state (at around pH 7), but becoming mostly protonated in the fully reduced state (similar to the behaviour of Glu 498 in Figure 5). On the other hand, the behaviour of Asp 112 is also similar to its equivalent Asp 116 in CotA, being mostly deprotonated in the whole of the pH range investigated, in all three states studied (results not shown). Therefore, these studies strongly suggest that the reported behaviour that we found in CotA laccase is probably general for other bacterial multicopper oxidases.
The present work suggests the importance of Glu 498 in the proton transfer mechanisms operating at the trinuclear centre and its vicinity. However, it should be remembered that this acidic group is not conserved in all laccases, although it is found, to our knowledge, in the same position in all the prokaryotic laccases that have been structurally characterised. In structurally characterised eukaryotic laccases, no acidic residue is found in this exact position, but an acid is always observed in the oxygen access channel (some examples can be found in the figure 4 containing the mechanism), suggesting a similar role for this group. As stated before, mutations were made  on one of these acidic groups on the eukaryotic multicopper oxidase Fet3p (Glu 487- equivalent to Glu 510 in ascorbate oxidase; figure 4d)), showing its importance in the mechanism; additionally, the authors suggested its role in proton transfer. It is foreseeable that equivalent acidic groups, or other acidic groups in the oxygen access channel of other multicopper oxidases, play a similar role. Nevertheless, the available structural characterisation of multicopper oxidases evidences (results not shown) that the distance between these acidic groups and the trinuclear centre is larger in eukaryotic than in prokaryotic enzymes. This difference may have some consequences in the details of the proton transfer mechanism, but the overall reported importance of these acid groups corroborates the hypothesis that they are responsible for the proton transfer mechanisms in multicopper oxidases.
The identification of Glu 498 as involved in proton transfer, should not rule out other residues in the vicinity of the trinuclear centre that may be important. One should not forget that our calculations address the thermodynamics of proton transfer, and their redox coupling, not considering the other component of this process, which is proton transfer kinetics that is much more complex to model (see for example [48, 49]). Therefore, there may be other important groups in this kinetic process that have not been identified yet, despite the fact that no clear-cut experimental evidence has been provided so far. Nevertheless, the influence on the mechanism evidenced by mutational and other studies of Asp 116 (CotA numbering) in Fet3p [16, 45] and in CueO [13, 14] may suggest such a role.
The differences in structurally characterised resting states observed for different laccases may result from the molecular details in the vicinity of trinuclear site. In addition to being observed in the CotA laccase , the dioxygen entity is only observed in the eukaryotic laccase from Melanocarpus albomices. However, in this latter case, the structure has a Cl- anion bound to the T2 copper ion, instead of a water or hydroxyl group, as observed in other laccases. It is tempting to speculate that the proximity of Glu 498, or equivalent groups, to the trinuclear centre in prokaryotic laccases (and to the dioxygen molecule itself) can be a factor for the stability of the dioxygen state, but other factors may exist. One example of the sensitivity of the trinuclear centre to the effects of its surroundings is the influence of the acidic group equivalent to Asp 116 in CotA (Figure 6) on the calculated differential stability of different conformations of the peroxyl intermediate . In fact, the trinuclear centre in multicopper oxidases seems to be fine tuned to achieve proper catalysis, as evidenced by a recent study focused on T3 copper sites . Further studies are clearly needed to further characterise this phenomenon.
In summary, these studies clearly indicate that in the CotA laccase, Glu498 in the entrance channel for the dioxygen moiety plays a role in channelling protons during the reduction process. They also show that the observation of a shoulder at 330 nm in the absorption spectrum, due to the antiferromagnetic coupling of the T3 copper ions, is present when the bridging species is either a hydroxyl or a dioxygen moiety. The results also strongly support the mechanism of dioxygen reduction as indicated in Figure 4. However, one key stage that requires further study is the mechanism by which the hydroxyl or water molecules produced by the reaction migrate past the T2 copper ion, before leaving the enzyme through the exit channel.
Recently, after the submission of this work, we made mutants of Glu498 in CotA, which show severe catalytic impairment when this group is substituted by threonine or leucine  Substitution by an aspartate renders an enzyme with 10% catalytic activity. These results are totally consistent with the results presented here, evidencing the need for an acid residue at this structural position.
Crystallization and Structure Solution
Data Collection statistics
Beam line at ESRF
Detector Distance (mm)
Cell parameters (Å), a
Oscillation range (°)
Oscillation angle (°)
No. of unique hkl
CE/MC calculations (described below) require special care concerning the consistency of the structural models to be compared (different redox situations and different liganted states of active site are examples). The main reason for this is the fact that CE/MC calculations use rigid structures (or semi-rigid, if on considers the inclusion of tautomers), making the simulation of the titration of protonatable groups quite susceptible to conformational changes occurring in polar and charged groups in their vicinity. This effect is especially significant in protein interiors. The inclusion of flexibility in the methodology would alleviate this susceptibility, but this is quite difficult nowadays. Therefore, to overcome these limitations of the methodology, besides using protein crystals obtained in conditions as similar as possible for the different situations, a special crystallographic refinement procedure was applied here; when electron density was not enough to clearly identify the conformation of residue side-chains (evidencing natural flexibility), these were placed in a conformation as close as possible to the previously determined holoCu(II) structure . This is very important, especially for charged side-chains, given that they can substantially influence the proton and electron binding calculations described next. Using this methodology, which we developed and used before in the analysis of multihaem cytochromes [25, 26] and CotA mutants , we can distinguish clear effects from noise coming from undetermined structural details.
UV/Visible absorption spectra for holoCotA and on apoCu(I) crystals were collected under cryo conditions (110 K) using an offline microspectrophotometer at the Cryobench Laboratory, ESRF, Grenoble, France .
Simulated pH and redox titrations were performed using methodologies, developed in our laboratory, for studying the binding equilibrium of protons and electrons [44, 47, 61]. These methodologies are based on continuum electrostatic (CE) methods and Monte-Carlo (MC) sampling of binding states. The CE methods (see  for a review), can be considered semi-microscopic and treat a protein like a low dielectric model, immersed in a high dielectric corresponding to the solvent (normally water), with charges placed at atomic positions. Electrostatic desolvation of charges and interactions between them are, therefore, treated by these methods, which normally solve the linear form of the Poisson-Boltzmann equation. The CE methods are used to compute the individual (pKint) and interaction terms of the binding free energy of the protons and electrons. These free energy terms are then used in a MC method (see details in [44, 47, 61]) that samples the binding states of protons and electrons. Additionally, and given that the protein is considered rigid in CE methods, this MC method treats protonatable groups in a multiconfigurational way that considers alternative positions for the protons (tautomers) . By the same limitation, the method also considers alternative configurations of the proton of alcohol groups (Ser and Thr), considering these configurations (three) as tautomers. However, given the high pKa of alcohol groups, here these are not considerable titrable. The procedure used allows us to deal with protonation-deprotonation and tautomeric exchange in a physically sound manner. The CE calculations were made using the package MEAD (version 2.2.0) [63, 64]. The sets of atomic radii and partial charges were taken from the GROMOS96 43A1 force field [65, 66], except for the metal centres, where quantum chemical calculations (see below) were used to derive charges. The dielectric constants used were 80 for the solvent and 10 for the protein, which are values within the range where pKa prediction is optimised . The solvent probe radius was 1.4 Å, the ion exclusion layer 2.0 Å, the ionic strength 0.1 M and the temperature 27°C. The program PETIT [47, 61] was used for the MC sampling of proton and electron binding states. Site pairs were selected for double moves when at least one pairwise term was greater than 2 pK units. Averages were computed using 105 MC steps. In all simulations, only the T1 centre was considered titrable, while the trinuclear centre was considered fixed in each state considered. Contrary to simulated pH titrations, which are done against a model compound reference (reflected in the pKa mod) and can be directly compared with experimental values, simulated redox titrations are usually made against an unknown reference, due to the unavailability of a corresponding Emod, i.e. the redox potential of an adequate model compound in water. This is due to the difficulty in finding such a model compound for protein metal co-factors, as it is the case for the T1 copper centre. The only option is to fit the value of Emod to experimental data on redox titrations of the same group in different proteins, but this is only usually exact if these proteins have only one type of redox centre (see for instance  for the case of bis-histidinyl coordinated tetrahaem cytochromes). However, in the case of multicopper oxidases, more than one redox centres is present, and further approximations need to be used. Here, and in a similar manner to previous work , the experimental value obtained for the holoCotA enzyme was used to fit the calculations of the enzyme structure with a dioxygen species in the trinuclear centre. This is a crude approximation, since the trinuclear centres are not allowed to titrate in the simulations, and, in the real system, the titrations of all redox centres of the protein may be coupled (and therefore exert an influence on the experimentally observed redox potential of the T1 copper centre). Nevertheless, in this work, we are interested, not on the exact redox potential, but on the relative behaviour of the T1 centre, and on the effects that other alterations of the protein (for instance, the redox state of the trinuclear centre) may bring to its redox affinity. For consistency, we have used the conformational state found in holoCu(II)  to model the holoCotA protein form. Since the present holoCotA structure contains an extra copper centre that has not been observed before, and we want to make comparisons with other CotA structures that do not contain it, we decided to use the equivalent holoCu(II) structure.
Partial charges for the T1 copper and the trinuclear centre were calculated considering model compounds for the different redox and liganted states analysed. These model compounds were created using coordinates of the centres as present in the structures with dioxygen, peroxide and in the fully reduced state as obtained previously  and the new structure of CotA containing hydroxyl obtained here. The ligands of the metals were included down to the C-beta carbon (C-alpha in the case of the cysteine residue of the T1 copper). A water molecule bound to the T2 copper was also included. Single point calculations were performed using Gaussian 98 , using B3LYP and the 6-31G(d) basis set for all atoms, with the exception of copper atoms where the 6-31G(2df) basis set was used. These calculations were used to derive electrostatic potentials, which were then fitted using RESP  to calculate the partial charges. These charges were used directly with no scaling applied. While the GROMOS 43A1 charges are not derived using the RESP procedure, we believe that, from a physical standpoint, this is the most convenient way of deriving partial charges for metal centres. For the T1 copper, partial charges were calculated for the oxidised and reduced states, since these two situations are needed to simulate the redox titration of this group.
António M. Baptista and Miguel Teixeira are gratefully acknowledged for helpful discussions. Maria Arménia Carrondo is gratefully acknowledged for support. This work was supported by project grants from Fundação para a Ciência e Tecnlogia, Portugal (POCI/BIO/57083/2004) and from European Union (BIORENEW - FP6-2004-NMP-NI-4/026456). C. Silva holds a Ph.D (SFRH/BD/40586/2007) fellowship and Z. Chen holds a Post-doc fellowship (SFRH/BBD/27104/2006) from Fundação para a Ciência e a Tecnologia, Portugal. The provision of synchrotron radiation facilities and the assistance of the macromolecular crystallography staff are sincerely acknowledged at the European Synchrotron Radiation Facility in Grenoble, France.
- Messerschmidt A: Multi-copper Oxidases. Singapore: World Science Press; 1997.Google Scholar
- Lindley PF: Multi-copper oxidases. Basel, New York: Marcel Dekker, Inc; 2001.Google Scholar
- Bento I, Carrondo MA, Lindley PF: Reduction of dioxygen by enzymes containing copper. J Biol Inorg Chem 2006, 11(5):539–547. 10.1007/s00775-006-0114-9View ArticlePubMedGoogle Scholar
- Stoj CS, Kosman DJ: Copper Oxidases. New York 2005.Google Scholar
- Hoegger PJ, Kilaru S, James TY, Thacker JR, Kues U: Phylogenetic comparison and classification of laccase and related multicopper oxidase protein sequences. FEBS J 2006, 273(10):2308–2326. 10.1111/j.1742-4658.2006.05247.xView ArticlePubMedGoogle Scholar
- Nakamura K, Go N: Function and molecular evolution of multicopper blue proteins. Cell Mol Life Sci 2005, 62(18):2050–2066. 10.1007/s00018-004-5076-xView ArticlePubMedGoogle Scholar
- Giardina P, Faraco V, Pezzella C, Piscitelli A, Vanhulle S, Sannia G: Laccases: a never-ending story. Cell Mol Life Sci 2010, 67(3):369–385. 10.1007/s00018-009-0169-1View ArticlePubMedGoogle Scholar
- Kosman DJ: Multicopper oxidases: a workshop on copper coordination chemistry, electron transfer, and metallophysiology. J Biol Inorg Chem 2010, 15(1):15–28. 10.1007/s00775-009-0590-9View ArticlePubMedGoogle Scholar
- Malmström BG: Enzymology of oxygen. Ann Rev Biochem 1982, 51: 21–59. 10.1146/annurev.bi.51.070182.000321View ArticlePubMedGoogle Scholar
- Solomon EI, Sundaram UM, Machonkin TE: Multicopper oxidases and oxygenases. Chem Rev 1996, 96: 2563–2605. 10.1021/cr950046oView ArticlePubMedGoogle Scholar
- Solomon EI, Chen P, Metz M, Lee SK, Palmer AE: Oxygen Binding, Activation, and Reduction to Water by Copper Proteins. Angew Chem Int Ed Engl 2001, 40(24):4570–4590. 10.1002/1521-3773(20011217)40:24<4570::AID-ANIE4570>3.0.CO;2-4View ArticlePubMedGoogle Scholar
- Sakurai T, Kataoka K: Structure and function of type I copper in multicopper oxidases. Cell Mol Life Sci 2007, 64(19–20):2642–2656. 10.1007/s00018-007-7183-yView ArticlePubMedGoogle Scholar
- Ueki Y, Inoue M, Kurose S, Kataoka K, Sakurai T: Mutations at Asp112 adjacent to the trinuclear Cu center in CueO as the proton donor in the four-electron reduction of dioxygen. FEBS Lett 2006, 580(17):4069–4072. 10.1016/j.febslet.2006.06.049View ArticlePubMedGoogle Scholar
- Kataoka K, Sugiyama R, Hirota S, Inoue M, Urata K, Minagawa Y, Seo D, Sakurai T: Four-electron reduction of dioxygen by a multicopper oxidase, CueO, and roles of Asp112 and Glu506 located adjacent to the trinuclear copper center. J Biol Chem 2009, 284(21):14405–14413. 10.1074/jbc.M808468200PubMed CentralView ArticlePubMedGoogle Scholar
- Quintanar L, Stoj C, Wang TP, Kosman DJ, Solomon EI: Role of aspartate 94 in the decay of the peroxide intermediate in the multicopper oxidase Fet3p. Biochemistry 2005, 44(16):6081–6091. 10.1021/bi047379cView ArticlePubMedGoogle Scholar
- Augustine AJ, Quintanar L, Stoj CS, Kosman DJ, Solomon EI: Spectroscopic and kinetic studies of perturbed trinuclear copper clusters: the role of protons in reductive cleavage of the O-O bond in the multicopper oxidase Fet3p. J Am Chem Soc 2007, 129(43):13118–13126. 10.1021/ja073905mPubMed CentralView ArticlePubMedGoogle Scholar
- Solomon EI, Augustine AJ, Yoon J: O2 reduction to H2O by the multicopper oxidases. Dalton Trans 2008, (30):3921–3932. 10.1039/b800799cGoogle Scholar
- Martins LO, Soares CM, Pereira MM, Teixeira M, Costa T, Jones GH, Henriques AO: Molecular and biochemical characterization of a highly stable bacterial laccase that occurs as a structural component of the Bacillus subtilis endospore coat. J Biol Chem 2002, 277(21):18849–18859. 10.1074/jbc.M200827200View ArticlePubMedGoogle Scholar
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res 2000, 28(1):235–242. 10.1093/nar/28.1.235PubMed CentralView ArticlePubMedGoogle Scholar
- Bento I, Martins LO, Gato Lopes G, Armenia Carrondo M, Lindley PF: Dioxygen reduction by multi-copper oxidases; a structural perspective. Dalton Trans 2005(21):3507–3513.Google Scholar
- Messerschmidt A, Ladenstein R, Huber R, Bolognesi M, Avigliano L, Petruzzelli R, Rossi A, Finazzi-Agro A: Refined crystal structure of ascorbate oxidase at 1.9 A resolution. J Mol Biol 1992, 224(1):179–205. 10.1016/0022-2836(92)90583-6View ArticlePubMedGoogle Scholar
- Durao P, Bento I, Fernandes AT, Melo EP, Lindley PF, Martins LO: Perturbations of the T1 copper site in the CotA laccase from Bacillus subtilis: structural, biochemical, enzymatic and stability studies. J Biol Inorg Chem 2006, 11(4):514–526. 10.1007/s00775-006-0102-0View ArticlePubMedGoogle Scholar
- Durao P, Chen Z, Silva CS, Soares CM, Pereira MM, Todorovic S, Hildebrandt P, Bento I, Lindley PF, Martins LO: Proximal mutations at the type 1 copper site of CotA laccase: spectroscopic, redox, kinetic and structural characterization of I494A and L386A mutants. Biochem J 2008, 412(2):339–346. 10.1042/BJ20080166View ArticlePubMedGoogle Scholar
- Durao P, Chen Z, Fernandes AT, Hildebrandt P, Murgida DH, Todorovic S, Pereira MM, Melo EP, Martins LO: Copper incorporation into recombinant CotA laccase from Bacillus subtilis: characterization of fully copper loaded enzymes. J Biol Inorg Chem 2008, 13(2):183–193. 10.1007/s00775-007-0312-0View ArticlePubMedGoogle Scholar
- Louro RO, Bento I, Matias PM, Catarino T, Baptista AM, Soares CM, Carrondo MA, Turner DL, Xavier AV: Conformational Component in the Coupled Transfer of Multiple Electrons and Protons in a Monomeric Tetraheme Cytochrome*. JBiolChem 2001, 276: 44044–44051.Google Scholar
- Bento I, Matias PM, Baptista AM, da Costa PN, van Dongen WM, Saraiva LM, Schneider TR, Soares CM, Carrondo MA: Molecular basis for redox-Bohr and cooperative effects in cytochrome c3 from Desulfovibrio desulfuricans ATCC 27774: crystallographic and modeling studies of oxidized and reduced high-resolution structures at pH 7.6. Proteins 2004, 54(1):135–152. 10.1002/prot.10431View ArticlePubMedGoogle Scholar
- Pereira MM, Sousa FL, Verissimo AF, Teixeira M: Looking for the minimum common denominator in haem-copper oxygen reductases: towards a unified catalytic mechanism. Biochim Biophys Acta 2008, 1777(7–8):929–934. 10.1016/j.bbabio.2008.05.441View ArticlePubMedGoogle Scholar
- Fadda E, Yu CH, Pomes R: Electrostatic control of proton pumping in cytochrome c oxidase. Biochim Biophys Acta 2008, 1777(3):277–284. 10.1016/j.bbabio.2007.11.010View ArticlePubMedGoogle Scholar
- Kannt A, Lancaster CR, Michel H: The coupling of electron transfer and proton translocation: electrostatic calculations on Paracoccus denitrificans cytochrome c oxidase. Biophys J 1998, 74(2 Pt 1):708–721. 10.1016/S0006-3495(98)73996-7PubMed CentralView ArticlePubMedGoogle Scholar
- Popovic DM, Stuchebrukhov AA: Electrostatic study of the proton pumping mechanism in bovine heart cytochrome C oxidase. J Am Chem Soc 2004, 126(6):1858–1871. 10.1021/ja038267wView ArticlePubMedGoogle Scholar
- Popovic DM, Stuchebrukhov AA: Proton exit channels in bovine cytochrome c oxidase. J Phys Chem B 2005, 109(5):1999–2006. 10.1021/jp0464371View ArticlePubMedGoogle Scholar
- Soares CM, Baptista AM, Pereira MM, Teixeira M: Investigation of protonatable residues in Rhodothermus marinus caa3 haem-copper oxygen reductase: comparison with Paracoccus denitrificans aa3 haem-copper oxygen reductase. J Biol Inorg Chem 2004, 9(2):124–134. 10.1007/s00775-003-0509-9View ArticlePubMedGoogle Scholar
- Enguita FJ, Martins LO, Henriques AO, Carrondo MA: Crystal Structure of a Bacterial Endospore Coat Component: A LACCASE WITH ENHANCED THERMOSTABILITY PROPERTIES. J Biol Chem 2003, 278(21):19416–19425. 10.1074/jbc.M301251200View ArticlePubMedGoogle Scholar
- Hakulinen N, Kruus K, Koivula A, Rouvinen J: A crystallographic and spectroscopic study on the effect of X-ray radiation on the crystal structure of Melanocarpus albomyces laccase. Biochem Biophys Res Commun 2006, 350(4):929–934. 10.1016/j.bbrc.2006.09.144View ArticlePubMedGoogle Scholar
- Shin W, Sundaram UM, Cole JL, Zhang HH, Hedman B, Hodgson KO, Solomon EI: Chemical and Spectroscopic Definition of the Peroxide-Level Intermediate in the Multicopper Oxidases: Relevance to the Catalytic Mechanism of Dioxygen Reduction to Water. J Am Chem Soc 1996, 118: 3202–3215. 10.1021/ja953621eView ArticleGoogle Scholar
- Palmer AE, Lee SK, Solomon EI: Decay of the peroxide intermediate in laccase: reductive cleavage of the O-O bond. J Am Chem Soc 2001, 123(27):6591–6599. 10.1021/ja010365zView ArticlePubMedGoogle Scholar
- Young T, Abel R, Kim B, Berne BJ, Friesner RA: Motifs for molecular recognition exploiting hydrophobic enclosure in protein-ligand binding. Proc Natl Acad Sci USA 2007, 104(3):808–813. 10.1073/pnas.0610202104PubMed CentralView ArticlePubMedGoogle Scholar
- Qvist J, Davidovic M, Hamelberg D, Halle B: A dry ligand-binding cavity in a solvated protein. Proc Natl Acad Sci USA 2008, 105(17):6296–6301. 10.1073/pnas.0709844105PubMed CentralView ArticlePubMedGoogle Scholar
- Berne BJ, Weeks JD, Zhou R: Dewetting and hydrophobic interaction in physical and biological systems. Annu Rev Phys Chem 2009, 60: 85–103. 10.1146/annurev.physchem.58.032806.104445PubMed CentralView ArticlePubMedGoogle Scholar
- Young T, Hua L, Huang X, Abel R, Friesner R, Berne BJ: Dewetting transitions in protein cavities. Proteins 2010, 78(8):1856–1869.PubMed CentralPubMedGoogle Scholar
- Cukierman S: Et tu, Grotthuss! and other unfinished stories. Biochim Biophys Acta 2006, 1757(8):876–885. 10.1016/j.bbabio.2005.12.001View ArticlePubMedGoogle Scholar
- Teixeira VH, Cunha CA, Machuqueiro M, Oliveira AS, Victor BL, Soares CM, Baptista AM: On the use of different dielectric constants for computing individual and pairwise terms in poisson-boltzmann studies of protein ionization equilibrium. J Phys Chem B 2005, 109(30):14691–14706. 10.1021/jp052259fView ArticlePubMedGoogle Scholar
- Polyakov KM, Fedorova TV, Stepanova EV, Cherkashin EA, Kurzeev SA, Strokopytov BV, Lamzin VS, Koroleva OV: Structure of native laccase from Trametes hirsuta at 1.8 A resolution. Acta Cryst D 2009, 65: 611–617. 10.1107/S0907444909011950View ArticleGoogle Scholar
- Baptista AM, Martel PJ, Soares CM: Simulation of electron-proton coupling with a Monte Carlo method: application to cytochrome c 3 using continuum electrostatics. BiophysJ 1999, 76: 2978–2998. 10.1016/S0006-3495(99)77452-7View ArticleGoogle Scholar
- Yoon J, Liboiron BD, Sarangi R, Hodgson KO, Hedman B, Solomon EI: The two oxidized forms of the trinuclear Cu cluster in the multicopper oxidases and mechanism for the decay of the native intermediate. Proc Natl Acad Sci USA 2007, 104(34):13609–13614. 10.1073/pnas.0705137104PubMed CentralView ArticlePubMedGoogle Scholar
- Roberts SA, Weichsel A, Grass G, Thakali K, Hazzard JT, Tollin G, Rensing C, Montfort WR: Crystal structure and electron transfer kinetics of CueO, a multicopper oxidase required for copper homeostasis in Escherichia coli. Proc Natl Acad Sci USA 2002, 99(5):2766–2771. 10.1073/pnas.052710499PubMed CentralView ArticlePubMedGoogle Scholar
- Teixeira VH, Soares CM, Baptista AM: Studies of the reduction and protonation behaviour of tetrahaem cytochromes using atomic detail. J BioI Inorg Chem 2002, 7: 200–216. 10.1007/s007750100287View ArticleGoogle Scholar
- Olsson MH, Siegbahn PE, Blomberg MR, Warshel A: Exploring pathways and barriers for coupled ET/PT in cytochrome c oxidase: a general framework for examining energetics and mechanistic alternatives. Biochim Biophys Acta 2007, 1767(3):244–260. 10.1016/j.bbabio.2007.01.015PubMed CentralView ArticlePubMedGoogle Scholar
- Olsson MH, Warshel A: Monte Carlo simulations of proton pumps: on the working principles of the biological valve that controls proton pumping in cytochrome c oxidase. Proc Natl Acad Sci USA 2006, 103(17):6500–6505. 10.1073/pnas.0510860103PubMed CentralView ArticlePubMedGoogle Scholar
- Hakulinen N, Kiiskinen LL, Kruus K, Saloheimo M, Paananen A, Koivula A, Rouvinen J: Crystal structure of a laccase from Melanocarpus albomyces with an intact trinuclear copper site. Nat Struct Biol 2002, 9(8):601–605.PubMedGoogle Scholar
- Yoon J, Fujii S, Solomon EI: Geometric and electronic structure differences between the type 3 copper sites of the multicopper oxidases and hemocyanin/tyrosinase. Proc Natl Acad Sci USA 2009, 106(16):6585–6590. 10.1073/pnas.0902127106PubMed CentralView ArticlePubMedGoogle Scholar
- Chen Z, Durao P, Silva CS, Pereira MM, Todorovic S, Hildebrandt P, Bento I, Lindley PF, Martins LO: The role of Glu(498) in the dioxygen reactivity of CotA-laccase from Bacillus subtilis. Dalton Trans 2010, 39(11):2875–2882. 10.1039/b922734bView ArticlePubMedGoogle Scholar
- Leslie A: Joint CCP4, ESF-EAMCB. Newsletter on Protein Crystallography 1992., 26:Google Scholar
- Leslie AG: The integration of macromolecular diffraction data. Acta Crystallogr D 2006, 62(Pt 1):48–57.View ArticlePubMedGoogle Scholar
- CCP4 CCP, Number 4: The CCP4 suite: programs for protein crystallography. Acta Cryst 1994, D50: 760–763.Google Scholar
- Vagin A, Teplyakov A: MOLREP: an Automated Program for Molecular Replacement. J Appl Cryst 1997, 30: 1022–1025. 10.1107/S0021889897006766View ArticleGoogle Scholar
- Murshudov GN, Vagin AA, Lebedev A, Wilson KS, Dodson EJ: Efficient anisotropic refinement of macromolecular structures using FFT. Acta Crystallogr D 1999, 55(Pt 1):247–255. 10.1107/S090744499801405XView ArticlePubMedGoogle Scholar
- Emsley P, Cowtan K: Coot: model-building tools for molecular graphics. Acta Crystallogr D Biol Crystallogr 2004, 60(Pt 12 Pt 1):2126–2132. 10.1107/S0907444904019158View ArticlePubMedGoogle Scholar
- Colloc'h N, Gabison L, Monard G, Altarsha M, Chiadmi M, Marassio G, Sopkova-de Oliveira Santos J, El Hajji M, Castro B, Abraini JH, et al.: Oxygen pressurized X-ray crystallography: probing the dioxygen binding site in cofactorless urate oxidase and implications for its catalytic mechanism. Biophys J 2008, 95(5):2415–2422. 10.1529/biophysj.107.122184PubMed CentralView ArticlePubMedGoogle Scholar
- Bourgeois D, Vernède X, Adam V, Fioravanti E, Ursby T: A microspectrophotometer for UV-visible and fluorescence studies of protein crystals. Journal of Applied Crystallography 2002, 35: 319–326. 10.1107/S0021889802003837View ArticleGoogle Scholar
- Baptista AM, Soares CM: Some theoretical and computational aspects of the inclusion of proton isomerism in the protonation equilibrium of proteins. JPhysChem B 2001, 105: 293–309.Google Scholar
- Bashford D: Macroscopic electrostatic models for protonation states in proteins. Front Biosci 2004, 9: 1082–1099. 10.2741/1187View ArticlePubMedGoogle Scholar
- Bashford D, Gerwert K: Electrostatic calculations of the pKa values of ionizable groups in bacteriorhodopsin. JMolBiol 1992, 224: 473–486.Google Scholar
- Bashford D: An Object-Oriented Programming Suite for Electrostatic Effects in Biological Molecules. In Scientific Computing in Object-Oriented Parallel Environments. Volume 1343. Edited by: Ishikawa Y, Oldehoeft RR, Reynders JVW, Tholburn M. Berlin: ISCOPE97, Springer; 1997:233–240. full_textView ArticleGoogle Scholar
- van Gunsteren WF, Billeter SR, Eising AA, Hunenberger PH, Kruger P, Mark AE, Scott WRP, Tironi IG: Biomolecular simulation: The GROMOS96 manual and user guide. Zurich, Groninger: BIOMOS b.v.; 1996.Google Scholar
- Scott WRP, Hünenberger PH, Tironi IG, Mark AE, Billeter SR, Fennen J, Torda AE, Huber T, Krüger P, van Gunsteren WF: The GROMOS biomolecular simulation program package. JPhysChem 1999, 103: 3596–3607.Google Scholar
- Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Zakrzewski VG, Montgomery J, J A, Stratmann RE, Burant JC, et al.: Gaussian 98, Revision A.7. Pittsburgh PA: Gaussian, Inc; 1998.Google Scholar
- Bayly CI, Cieplak P, Cornell WD, Kollman PA: A well-behaved electrostatic potential based method using charge restraints for deriving atomic charges: the RESP model. JPhysChem 1993, 97: 10269–10280.Google Scholar
- DeLano WL: The PyMOL Molecular Graphics System. San Carlos, CA, USA: DeLano Scientific; 2002.Google Scholar
- Ramakrishnan C, Ramachandran GN: Stereochemical criteria for polypeptide and protein chain conformations. II. Allowed conformations for a pair of peptide units. Biophys J 1965, 5(6):909–933. 10.1016/S0006-3495(65)86759-5PubMed CentralView ArticlePubMedGoogle Scholar
- Lovell SC, Davis IW, Arendall WB, de Bakker PI, Word JM, Prisant MG, Richardson JS, Richardson DC: Structure validation by Calpha geometry: phi,psi and Cbeta deviation. Proteins 2003, 50(3):437–450. 10.1002/prot.10286View ArticlePubMedGoogle Scholar