Signatures in SARS-CoV-2 spike protein conferring escape to neutralizing antibodies.

Understanding SARS-CoV-2 evolution and host immunity is critical to control COVID-19 pandemics. At the core is an arms-race between SARS-CoV-2 antibody and angiotensin-converting enzyme 2 (ACE2) recognition, a function of the viral protein spike. Mutations in spike impacting antibody and/or ACE2 binding are appearing worldwide, imposing the need to monitor SARS-CoV2 evolution and dynamics in the population. Determining signatures in SARS-CoV-2 that render the virus resistant to neutralizing antibodies is critical. We engineered 25 spike-pseudotyped lentiviruses containing individual and combined mutations in the spike protein, including all defining mutations in the variants of concern, to identify the effect of single and synergic amino acid substitutions in promoting immune escape. We confirmed that E484K evades antibody neutralization elicited by infection or vaccination, a capacity augmented when complemented by K417N and N501Y mutations. In silico analysis provided an explanation for E484K immune evasion. E484 frequently engages in interactions with antibodies but not with ACE2. Importantly, we identified a novel amino acid of concern, S494, which shares a similar pattern. Using the already circulating mutation S494P, we found that it reduces antibody neutralization of convalescent and post-immunization sera, particularly when combined with E484K and with mutations able to increase binding to ACE2, such as N501Y. Our analysis of synergic mutations provides a signature for hotspots for immune evasion and for targets of therapies, vaccines and diagnostics.

F508del disturbs the dynamics of the nucleotide binding domains of CFTR before and after ATP hydrolysis.

The cystic fibrosis transmembrane conductance regulator (CFTR) channel is an ion channel responsible for chloride transport in epithelia and it belongs to the class of ABC transporters. The deletion of phenylalanine 508 (F508del) in CFTR is the most common mutation responsible for cystic fibrosis. Little is known about the effect of the mutation in the isolated nucleotide binding domains (NBDs), on dimer dynamics, ATP hydrolysis and even on nucleotide binding. Using molecular dynamics simulations of the human CFTR NBD dimer, we showed that F508del increases, in the prehydrolysis state, the inter-motif distance in both ATP binding sites (ABP) when ATP is bound. Additionally, a decrease in the number of catalytically competent conformations was observed in the presence of F508del. We used the subtraction technique to study the first 300 ps after ATP hydrolysis in the catalytic competent site and found that the F508del dimer evidences lower conformational changes than the wild type. Using longer simulation times, the magnitude of the conformational changes in both forms increases. Nonetheless, the F508del dimer shows lower C-α RMS values in comparison to the wild-type, on the F508del loop, on the residues surrounding the catalytic site and the portion of NBD2 adjacent to ABP1. These results provide evidence that F508del interferes with the NBD dynamics before and after ATP hydrolysis. These findings shed a new light on the effect of F508del on NBD dynamics and reveal a novel mechanism for the influence of F508del on CFTR.

The rice cold-responsive calcium-dependent protein kinase OsCPK17 is regulated by alternative splicing and post-translational modifications.

Plant calcium-dependent protein kinases (CDPKs) are key proteins implicated in calcium-mediated signaling pathways of a wide range of biological events in the organism. The action of each particular CDPK is strictly regulated by many mechanisms in order to ensure an accurate signal translation and the activation of the adequate response processes. In this work, we investigated the regulation of a CDPK involved in rice cold stress response, OsCPK17, to better understand its mode of action. We identified two new alternative splicing (AS) mRNA forms of OsCPK17 encoding truncated versions of the protein, missing the CDPK activation domain. We analyzed the expression patterns of all AS variants in rice tissues and examined their subcellular localization in onion epidermal cells. The results indicate that the AS of OsCPK17 putatively originates truncated forms of the protein with distinct functions, and different subcellular and tissue distributions. Additionally, we addressed the regulation of OsCPK17 by post-translational modifications in several in vitro experiments. Our analysis indicated that OsCPK17 activity depends on its structural rearrangement induced by calcium binding, and that the protein can be autophosphorylated. The identified phosphorylation sites mostly populate the OsCPK17 N-terminal domain. Exceptions are phosphosites T107 and S136 in the kinase domain and S558 in the C-terminal domain. These phosphosites seem conserved in CDPKs and may reflect a common regulatory mechanism for this protein family.

A LysM Domain Intervenes in Sequential Protein-Protein and Protein-Peptidoglycan Interactions Important for Spore Coat Assembly in Bacillus subtilis.

At a late stage in spore development in Bacillus subtilis, the mother cell directs synthesis of a layer of peptidoglycan known as the cortex between the two forespore membranes, as well as the assembly of a protective protein coat at the surface of the forespore outer membrane. SafA, the key determinant of inner coat assembly, is first recruited to the surface of the developing spore and then encases the spore under the control of the morphogenetic protein SpoVID. SafA has a LysM peptidoglycan-binding domain, SafALysM, and localizes to the cortex-coat interface in mature spores. SafALysM is followed by a region, A, required for an interaction with SpoVID and encasement. We now show that residues D10 and N30 in SafALysM, while involved in the interaction with peptidoglycan, are also required for the interaction with SpoVID and encasement. We further show that single alanine substitutions on residues S11, L12, and I39 of SafALysM that strongly impair binding to purified cortex peptidoglycan affect a later stage in the localization of SafA that is also dependent on the activity of SpoVE, a transglycosylase required for cortex formation. The assembly of SafA thus involves sequential protein-protein and protein-peptidoglycan interactions, mediated by the LysM domain, which are required first for encasement then for the final localization of the protein in mature spores.IMPORTANCEBacillus subtilis spores are encased in a multiprotein coat that surrounds an underlying peptidoglycan layer, the cortex. How the connection between the two layers is enforced is not well established. Here, we elucidate the role of the peptidoglycan-binding LysM domain, present in two proteins, SafA and SpoVID, that govern the localization of additional proteins to the coat. We found that SafALysM is a protein-protein interaction module during the early stages of coat assembly and a cortex-binding module at late stages in morphogenesis, with the cortex-binding function promoting a tight connection between the cortex and the coat. In contrast, SpoVIDLysM functions only as a protein-protein interaction domain that targets SpoVID to the spore surface at the onset of coat assembly.

Unexplored Nucleotide Binding Modes for the ABC Exporter MsbA.

The ATP-binding cassette (ABC) transporter MsbA is an ATP-driven lipid-A flippase. It belongs to the ABC protein superfamily whose members are characterized by conserved motifs in their nucleotide binding domains (NBDs), which are responsible for ATP hydrolysis. Recently, it was found that MsbA could catalyze a reverse adenylate kinase (rAK)-like reaction in addition to ATP hydrolysis. Both reactions are connected and mediated by the same conserved NBD domains. Here, the structural foundations underlying the nucleotide binding to MsbA were therefore explored using a concerted approach based on conventional- and DNP-enhanced solid-state NMR, pulsed-EPR, and MD simulations. MsbA reconstituted into lipid bilayers was trapped in various catalytic states corresponding to intermediates of the coupled ATPase-rAK mechanism. The analysis of nucleotide-binding dependent chemical shift changes, and the detection of through-space contacts between bound nucleotides and MsbA within these states provides evidence for an additional nucleotide-binding site in close proximity to the Q-loop and the His-Switch. By replacing Mg2+ with Mn2+ and employing pulsed EPR spectroscopy, evidence is provided that this newly found nucleotide binding site does not interfere with the coordination of the required metal ion. Molecular dynamic (MD) simulations of nucleotide and metal binding required for the coupled ATPase-rAK mechanism have been used to corroborate these experimental findings and provide additional insight into nucleotide location, orientation, and possible binding modes.

Molecular structure of FoxE, the putative iron oxidase of Rhodobacter ferrooxidans SW2.

The ancient metabolism of photoferrotrophy is likely to have played a key role in the biogeochemical cycle of iron on Early Earth leading to the deposition of Banded Iron Formations prior to the emergence of oxygenic photosynthesis. Extant organisms still performing this metabolism provide a convenient window to peer into its molecular mechanisms. Here we report the molecular structure of FoxE, the putative terminal iron oxidase of Rhodobacter ferrooxidans SW2. This protein is organized as a trimer with two hemes and a disulfide bridge per monomer. The distance between hemes, their solvent exposure and the surface electrostatics ensure a controlled electron transfer rate. They also guarantee segregation between electron capture from ferrous iron and electron release to downstream acceptors, which do not favor the precipitation of ferric iron. Combined with the functional characterization of this protein, the structure reveals how iron oxidation can be performed in the periplasmic space of this Gram-negative bacterium at circumneutral pH, while minimizing the risk of mineral precipitation and cell encrustation.

The key role of glutamate 172 in the mechanism of type II NADH:quinone oxidoreductase of Staphylococcus aureus.

Type II NADH:quinone oxidoreductases (NDH-2s) are membrane bound enzymes that deliver electrons to the respiratory chain by oxidation of NADH and reduction of quinones. In this way, these enzymes also contribute to the regeneration of NAD+, allowing several metabolic pathways to proceed. As for the other members of the two-Dinucleotide Binding Domains Flavoprotein (tDBDF) superfamily, the enzymatic mechanism of NDH-2s is still little explored and elusive. In this work we addressed the role of the conserved glutamate 172 (E172) residue in the enzymatic mechanism of NDH-2 from Staphylococcus aureus. We aimed to test our earlier hypothesis that E172 plays a key role in proton transfer to allow the protonation of the quinone. For this we performed a complete biochemical characterization of the enzyme's variants E172A, E172Q and E172S. Our steady state kinetic measurements show a clear decrease in the overall reaction rate, and our substrate interaction studies indicate the binding of the two substrates is also affected by these mutations. Interestingly our fast kinetic results show quinone reduction is more affected than NADH oxidation. We have also determined the X-ray crystal structure of the E172S mutant (2.55Ǻ) and compared it with the structure of the wild type (2.32Ǻ). Together these results support our hypothesis for E172 being of central importance in the catalytic mechanism of NDH-2, which may be extended to other members of the tDBDF superfamily.

Coupling between protonation and conformation in cytochrome c oxidase: Insights from constant-pH MD simulations.

Cytochrome c oxidases (CcOs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes reduce dioxygen (O(2)) to water and, simultaneously, generate a transmembrane electrochemical proton gradient. Despite their importance in the aerobic metabolism and the large amount of structural and biochemical data available for the A1-type CcO family, there is still no consensually accepted description of the molecular mechanisms operating in this protein. A substantial number of questions about the CcO's working mechanism remain to be answered, including how the protonation behavior of some key residues is modulated during a reduction cycle and how is the conformation of the protein affected by protonation. The main objective of this work was to study the protonation-conformation coupling in CcOs and identify the molecular factors that control the protonation state of some key residues. In order to directly capture the interplay between protonation and conformational effects, we have performed constant-pH MD simulations of an A1-type CcO inserted into a lipid bilayer in two redox states (oxidized and reduced) at physiological pH. From the simulations, we were able to identify several groups with unusual titration behavior that are highly dependent on the protein redox state, including the A-propionate from heme a and the D-propionate from heme a3, two key groups possibly involved in proton pumping. The protonation state of these two groups is heavily influenced by subtle conformational changes in the protein (notably of R481(I) and R482(I)) and by small changes in the hydrogen bond network.

Effect of a pH Gradient on the Protonation States of Cytochrome c Oxidase: A Continuum Electrostatics Study.

Cytochrome c oxidase (CcO) couples the reduction of dioxygen to water with transmembrane proton pumping, which leads to the generation of an electrochemical gradient. In this study we analyze how one of the components of the electrochemical gradient, the difference in pH across the membrane, or ΔpH, influences the protonation states of residues in CcO. We modified our continuum electrostatics/Monte Carlo (CE/MC) method in order to include the ΔpH and applied it to the study of CcO, in what is, to our best knowledge, the first CE/MC study of CcO in the presence of a pH gradient. The inclusion of a transmembrane pH gradient allows for the identification of residues whose titration behavior depends on the pH on both sides of the membrane. Among the several residues with unusual titration profiles, three are well-known key residues in the proton transfer process of CcO: E286I, Y288I, and K362I. All three residues have been previously identified as being critical for the catalytic or proton pumping functions of CcO. Our results suggest that when the pH gradient increases, these residues may be part of a regulatory mechanism to stem the proton flow.

Structural and Functional Characterization of an Ancient Bacterial Transglutaminase Sheds Light on the Minimal Requirements for Protein Cross-Linking.

Transglutaminases are best known for their ability to catalyze protein cross-linking reactions that impart chemical and physical resilience to cellular structures. Here, we report the crystal structure and characterization of Tgl, a transglutaminase from the bacterium Bacillus subtilis. Tgl is produced during sporulation and cross-links the surface of the highly resilient spore. Tgl-like proteins are found only in spore-forming bacteria of the Bacillus and Clostridia classes, indicating an ancient origin. Tgl is a single-domain protein, produced in active form, and the smallest transglutaminase characterized to date. We show that Tgl is structurally similar to bacterial cell wall endopeptidases and has an NlpC/P60 catalytic core, thought to represent the ancestral unit of the cysteine protease fold. We show that Tgl functions through a unique partially redundant catalytic dyad formed by Cys116 and Glu187 or Glu115. Strikingly, the catalytic Cys is insulated within a hydrophobic tunnel that traverses the molecule from side to side. The lack of similarity of Tgl to other transglutaminases together with its small size suggests that an NlpC/P60 catalytic core and insulation of the active site during catalysis may be essential requirements for protein cross-linking.

Exploring O2 diffusion in A-type cytochrome c oxidases: molecular dynamics simulations uncover two alternative channels towards the binuclear site.

Cytochrome c oxidases (Ccoxs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes couple dioxygen (O2) reduction to the generation of a transmembrane electrochemical proton gradient. Despite decades of research and the availability of a large amount of structural and biochemical data available for the A-type Ccox family, little is known about the channel(s) used by O2 to travel from the solvent/membrane to the heme a3-CuB binuclear center (BNC). Moreover, the identification of all possible O2 channels as well as the atomic details of O2 diffusion is essential for the understanding of the working mechanisms of the A-type Ccox. In this work, we determined the O2 distribution within Ccox from Rhodobacter sphaeroides, in the fully reduced state, in order to identify and characterize all the putative O2 channels leading towards the BNC. For that, we use an integrated strategy combining atomistic molecular dynamics (MD) simulations (with and without explicit O2 molecules) and implicit ligand sampling (ILS) calculations. Based on the 3D free energy map for O2 inside Ccox, three channels were identified, all starting in the membrane hydrophobic region and connecting the surface of the protein to the BNC. One of these channels corresponds to the pathway inferred from the X-ray data available, whereas the other two are alternative routes for O2 to reach the BNC. Both alternative O2 channels start in the membrane spanning region and terminate close to Y288I. These channels are a combination of multiple transiently interconnected hydrophobic cavities, whose opening and closure is regulated by the thermal fluctuations of the lining residues. Furthermore, our results show that, in this Ccox, the most likely (energetically preferred) routes for O2 to reach the BNC are the alternative channels, rather than the X-ray inferred pathway.

Self-assembly molecular dynamics simulations shed light into the interaction of the influenza fusion Peptide with a membrane bilayer.

Influenza virus is one of the most devastating human pathogens. In order to infect host cells, this virus fuses its membrane with the host membrane in a process mediated by the glycoprotein hemagglutinin. During fusion, the N-terminal region of hemagglutinin, which is known as the fusion peptide (FP), inserts into the host membrane, promoting lipid mixing between the viral and host membranes. Therefore, this peptide plays a key role in the fusion process, but the exact mechanism by which it promotes lipid mixing is still unclear. To shed light into this matter, we performed molecular dynamics (MD) simulations of the influenza FP in different environments (water, dodecylphosphocholine (DPC) micelles, and a dimyristoylphosphatidylcholine (DMPC) membrane). While in pure water the peptide lost its initial secondary structure, in simulations performed in the presence of DPC micelles it remained stable, in agreement with previous experimental observations. In simulations performed in the presence of a preassembled DMPC bilayer, the peptide became unstructured and was unable to insert into the membrane as a result of technical limitations of the method used. To overcome this problem, we used a self-assembly strategy, assembling the membrane together with the peptide. These simulations revealed that the peptide can adopt a membrane-spanning conformation, which had not been predicted by previous MD simulation studies. The peptide insertion had a strong effect on the membrane, lowering the bilayer thickness, disordering nearby lipids, and promoting lipid tail protrusion. These results contribute to a better understanding of the role of the FP in the fusion process.

Mind the gap: cytochrome interactions reveal electron pathways across the periplasm of Shewanella oneidensis MR-1.

Extracellular electron transfer is the key metabolic trait that enables some bacteria to play a significant role in the biogeochemical cycling of metals and in bioelectrochemical devices such as microbial fuel cells. In Shewanella oneidensis MR-1, electrons generated in the cytoplasm by catabolic processes must cross the periplasmic space to reach terminal oxidoreductases found at the cell surface. Lack of knowledge on how these electrons flow across the periplasmic space is one of the unresolved issues related with extracellular electron transfer. Using NMR to probe protein-protein interactions, kinetic measurements of electron transfer and electrostatic calculations, we were able to identify protein partners and their docking sites, and determine the dissociation constants. The results showed that both STC (small tetrahaem cytochrome c) and FccA (flavocytochrome c) interact with their redox partners, CymA and MtrA, through a single haem, avoiding the establishment of stable redox complexes capable of spanning the periplasmic space. Furthermore, we verified that the most abundant periplasmic cytochromes STC, FccA and ScyA (monohaem cytochrome c5) do not interact with each other and this is likely to be the consequence of negative surface charges in these proteins. This reveals the co-existence of two non-mixing redox pathways that lead to extracellular electron transfer in S. oneidensis MR-1 established through transient protein interactions.

Computer-aided design and implementation of efficient biosynthetic pathways to produce high added-value products derived from tyrosine in Escherichia coli.

Developing efficient bioprocesses requires selecting the best biosynthetic pathways, which can be challenging and time-consuming due to the vast amount of data available in databases and literature. The extension of the shikimate pathway for the biosynthesis of commercially attractive molecules often involves promiscuous enzymes or lacks well-established routes. To address these challenges, we developed a computational workflow integrating enumeration/retrosynthesis algorithms, a toolbox for pathway analysis, enzyme selection tools, and a gene discovery pipeline, supported by manual curation and literature review. Our focus has been on implementing biosynthetic pathways for tyrosine-derived compounds, specifically L-3,4-dihydroxyphenylalanine (L-DOPA) and dopamine, with significant applications in health and nutrition. We selected one pathway to produce L-DOPA and two different pathways for dopamine-one already described in the literature and a novel pathway. Our goal was either to identify the most suitable gene candidates for expression in Escherichia coli for the known pathways or to discover innovative pathways. Although not all implemented pathways resulted in the accumulation of target compounds, in our shake-flask experiments we achieved a maximum L-DOPA titer of 0.71 g/L and dopamine titers of 0.29 and 0.21 g/L for known and novel pathways, respectively. In the case of L-DOPA, we utilized, for the first time, a mutant version of tyrosinase from Ralstonia solanacearum. Production of dopamine via the known biosynthesis route was accomplished by coupling the L-DOPA pathway with the expression of DOPA decarboxylase from Pseudomonas putida, resulting in a unique biosynthetic pathway never reported in literature before. In the context of the novel pathway, dopamine was produced using tyramine as the intermediate compound. To achieve this, tyrosine was initially converted into tyramine by expressing TDC from Levilactobacillus brevis, which, in turn, was converted into dopamine through the action of the enzyme encoded by ppoMP from Mucuna pruriens. This marks the first time that an alternative biosynthetic pathway for dopamine has been validated in microbes. These findings underscore the effectiveness of our computational workflow in facilitating pathway enumeration and selection, offering the potential to uncover novel biosynthetic routes, thus paving the way for other target compounds of biotechnological interest.

Putative dioxygen-binding sites and recognition of tigecycline and minocycline in the tetracycline-degrading monooxygenase TetX.

Expression of the aromatic hydroxylase TetX under aerobic conditions confers bacterial resistance against tetracycline antibiotics. Hydroxylation inactivates and degrades tetracyclines, preventing inhibition of the prokaryotic ribosome. X-ray crystal structure analyses of TetX in complex with the second-generation and third-generation tetracyclines minocycline and tigecycline at 2.18 and 2.30 Šresolution, respectively, explain why both clinically potent antibiotics are suitable substrates. Both tetracyclines bind in a large tunnel-shaped active site in close contact to the cofactor FAD, pre-oriented for regioselective hydroxylation to 11a-hydroxytetracyclines. The characteristic bulky 9-tert-butylglycylamido substituent of tigecycline is solvent-exposed and does not interfere with TetX binding. In the TetX-minocycline complex a second binding site for a minocycline dimer is observed close to the active-site entrance. The pocket is formed by the crystal packing arrangement on the surface of two neighbouring TetX monomers. Crystal structure analysis at 2.73 Šresolution of xenon-pressurized TetX identified two adjacent Xe-binding sites. These putative dioxygen-binding cavities are located in the substrate-binding domain next to the active site. Molecular-dynamics simulations were performed in order to characterize dioxygen-diffusion pathways to FADH2 at the active site.

Influence of the protein structure surrounding the active site on the catalytic activity of [NiFeSe] hydrogenases.

A combined experimental and theoretical study of the catalytic activity of a [NiFeSe] hydrogenase has been performed by H/D exchange mass spectrometry and molecular dynamics simulations. Hydrogenases are enzymes that catalyze the heterolytic cleavage or production of H2. The [NiFeSe] hydrogenases belong to a subgroup of the [NiFe] enzymes in which a selenocysteine is a ligand of the nickel atom in the active site instead of cysteine. The aim of this research is to determine how much the specific catalytic properties of this hydrogenase are influenced by the replacement of a sulfur by selenium in the coordination of the bimetallic active site versus the changes in the protein structure surrounding the active site. The pH dependence of the D2/H(+) exchange activity and the high isotope effect observed in the Michaelis constant for the dihydrogen substrate and in the single exchange/double exchange ratio suggest that a "cage effect" due to the protein structure surrounding the active site is modulating the enzymatic catalysis. This "cage effect" is supported by molecular dynamics simulations of the diffusion of H2 and D2 from the outside to the inside of the protein, which show different accumulation of these substrates in a cavity next to the active site.

A molecular perspective on nonaqueous biocatalysis: contributions from simulation studies.

The technological value of nonaqueous enzymology has been recognized for more than thirty years. A detailed understanding of the molecular determinants of enzyme behaviour in nonaqueous media is essential to explore their potential. Computer simulations have provided valuable contributions to this field, having elucidated how the solvent affects the structural and dynamic properties of enzymes, as well as their activity and enantioselectivity. They have also helped to shed light on the effect of hydration and the role of counterions. In this perspective, we describe the major challenges and achievements of molecular simulations of enzymes in different types of nonaqueous solvents, including organic solvents, ionic liquids and supercritical fluids.

Microbial production of the plant flavanone hesperetin from caffeic acid.

OBJECTIVE: Hesperetin is an important O-methylated flavonoid produced by citrus fruits and of potential pharmaceutical relevance. The microbial biosynthesis of hesperetin could be a viable alternative to plant extraction, as plant extracts often yield complex mixtures of different flavonoids making it challenging to isolate pure compounds. In this study, hesperetin was produced from caffeic acid in the microbial host Escherichia coli. We combined a previously optimised pathway for the biosynthesis of the intermediate flavanone eriodictyol with a combinatorial library of plasmids expressing three candidate flavonoid O-methyltransferases. Moreover, we endeavoured to improve the position specificity of CCoAOMT7, a flavonoid O-methyltransferase from Arabidopsis thaliana that has been demonstrated to O-methylate eriodictyol in both the para- and meta-position, thus leading to a mixture of hesperetin and homoeriodictyol. RESULTS: The best performing flavonoid O-methyltransferase in our screen was found to be CCoAOMT7, which could produce up to 14.6 mg/L hesperetin and 3.8 mg/L homoeriodictyol from 3 mM caffeic acid in E. coli 5-alpha. Using a platform for enzyme engineering that scans the mutational space of selected key positions, predicting their structures using homology modelling and inferring their potential catalytic improvement using docking simulations, we were able to identify a CCoAOMT7 mutant with a two-fold higher position specificity for hesperetin. The mutant's catalytic activity, however, was considerably diminished. Our findings suggest that hesperetin can be created from central carbon metabolism in E. coli following the introduction of a caffeic acid biosynthesis pathway.

Molecular modeling of hair keratin/peptide complex: Using MM-PBSA calculations to describe experimental binding results.

Molecular dynamics simulations of a keratin/peptide complex have been conducted to predict the binding affinity of four different peptides toward human hair. Free energy calculations on the peptides' interaction with the keratin model demonstrated that electrostatic interactions are believed to be the main driving force stabilizing the complex. The molecular mechanics-Poisson-Boltzmann surface area methodology used for the free energy calculations demonstrated that the dielectric constant in the protein's interior plays a major role in the free energy calculations, and the only way to obtain accordance between the free energy calculations and the experimental binding results was to use the average dielectric constant.

The role of Asp116 in the reductive cleavage of dioxygen to water in CotA laccase: assistance during the proton-transfer mechanism.

Multi-copper oxidases constitute a family of proteins that are capable of coupling the one-electron oxidation of four substrate equivalents to the four-electron reduction of dioxygen to two molecules of water. The main catalytic stages occurring during the process have already been identified, but several questions remain, including the nature of the protonation events that take place during the reductive cleavage of dioxygen to water. The presence of a structurally conserved acidic residue (Glu498 in CotA laccase from Bacillus subtilis) at the dioxygen-entrance channel has been reported to play a decisive role in the protonation mechanisms, channelling protons during the reduction process and stabilizing the site as a whole. A second acidic residue that is sequentially conserved in multi-copper oxidases and sited within the exit channel (Asp116 in CotA) has also been identified as being important in the protonation process. In this study, CotA laccase has been used as a model system to assess the role of Asp116 in the reduction process of dioxygen to water. The crystal structures of three distinct mutants, D116E, D116N and D116A, produced by site-saturation mutagenesis have been determined. In addition, theoretical calculations have provided further support for a role of this residue in the protonation events.