Crystal structure of LAAO from Calloselasma rhodostoma with an L-phenylalanine substrate: insights into structure and mechanism

L-Amino acid oxidase is a dimeric glycosylated flavoenzyme, a major constituent of the venom-from the snake Calloselasma rhodostoma. The enzyme exhibits apoptosis inducing effects as well as antibacterial and anti-HIV activities. The structure of l-amino acid oxidase with its substrate (L-phenylalanine) has been refined to a resolution of 1.8 A. The complex structure reveals the substrate bound to the reduced flavin (FADred). Alternative conformations for the key residues His223 and Arg322 are evident, suggesting a dynamic active site. Furthermore, conformational changes are apparent for the isoalloxazine ring; the three-ring system exhibits more bending around the N5-N10 axis compared to the oxidized flavin. The implications of the observed dynamics on the mechanism of catalysis are discussed. Inspection of buried surfaces in the enzyme reveals a Y-shaped channel system extending from the external surface of the protein to the active site. One portion of this channel may serve as the entry path for O2 during the oxidative half-reaction. The second region, separated from the proposed O2 channel by the N terminus (residues 8-16) of the protein, may play a role in H2O2 release. Interestingly, the latter portion of the channel would direct the H2O2 product to the exterior surface of the protein, near the glycan moiety, thought to anchor the enzyme to the host cell. This channel location may explain the ability of the enzyme to localize H2O2 to the targeted cell and thus induce the apoptotic effect.     

Structure of the DPS-like protein from Sulfolobus solfataricus reveals a bacterioferritin-like dimetal binding site within a DPS-like dodecameric assembly

The superfamily of ferritin-like proteins has recently expanded to include a phylogenetically distinct class of proteins termed DPS-like (DPSL) proteins. Despite their distinct genetic signatures, members of this subclass share considerable similarity to previously recognized DPS proteins. Like DPS, these proteins are expressed in response to oxidative stress, form dodecameric cage-like particles, preferentially utilize H(2)O(2) in the controlled oxidation of Fe(2+), and possess a short N-terminal extension implicated in stabilizing cellular DNA. Given these extensive similarities, the functional properties responsible for the preservation of the DPSL signature in the genomes of diverse prokaryotes have been unclear. Here, we describe the crystal structure of a DPSL protein from the thermoacidophilic archaeon Sulfolobus solfataricus. Although the overall fold of the polypeptide chain and the oligomeric state of this protein are indistinguishable from those of authentic DPS proteins, several important differences are observed. First, rather than a ferroxidase site at the subunit interface, as is observed in all other DPS proteins, the ferroxidase site in SsDPSL is buried within the four-helix bundle, similar to bacterioferritin. Second, the structure reveals a channel leading from the exterior surface of SsDPSL to the bacterioferritin-like dimetal binding site, possibly allowing divalent cations and/or H(2)O(2) to access the active site. Third, a pair of cysteine residues unique to DPSL proteins is found adjacent to the dimetal binding site juxtaposed between the exterior surface of the protein and the active site channel. The cysteine residues in this thioferritin motif may play a redox active role, possibly serving to recycle iron at the ferroxidase center.     

Fungal catalases: Function, phylogenetic origin and structure

Most fungi have several monofunctional heme-catalases. Filamentous ascomycetes (Pezizomycotina) have two types of large-size subunit catalases (L1 and L2). L2-type are usually induced by different stressors and are extracellular enzymes; those from the L1-type are not inducible and accumulate in asexual spores. L2 catalases are important for growth and the start of cell differentiation, while L1 are required for spore germination. In addition, pezizomycetes have one to four small-size subunit catalases. Yeasts (Saccharomycotina) do not have large-subunit catalases and generally have one peroxisomal and one cytosolic small-subunit catalase. Small-subunit catalases are inhibited by substrate while large-subunit catalases are activated by H(2)O(2). Some small-subunit catalases bind NADPH preventing inhibition by substrate. We present a phylogenetic analysis revealing one or two events of horizontal gene transfers from Actinobacteria to a fungal ancestor before fungal diversification, as the origin of large-size subunit catalases. Other possible horizontal transfers of small- and large-subunit catalases genes were detected and one from bacteria to the fungus Malassezia globosa was analyzed in detail. All L2-type catalases analyzed presented a secretion signal peptide. Mucorales preserved only L2-type catalases, with one containing a secretion signal if two or more are present. Basidiomycetes have only L1-type catalases, all lacking signal peptide. Fungal small-size catalases are related to animal catalases and probably evolved from a common ancestor. However, there are several groups of small-size catalases. In particular, a conserved group of fungal sequences resemble plant catalases, whose phylogenetic origin was traced to a group of bacteria. This group probably has the heme orientation of plant catalases and could in principle bind NADPH. From almost a hundred small-subunit catalases only one fourth has a peroxisomal localization signal and in fact many fungi lack a peroxisomal catalase. Catalases have a deep buried active site and H(2)O(2) has to go through a long passage to reach it. In all known structures of catalases, the major channel has common features, particularly in the straight and narrow final section that is positioned perpendicular to the heme. Besides, other conserved channels are present in catalases whose function remains to be elucidated. One of these channels intercommunicates the major channels from the two R-related subunits. In three of the four known large-subunits catalase structures, the heme b is partially transformed into heme d. In Neurospora crassa, this occurs in vivo and is related to oxidative stress conditions in which singlet oxygen is produced. A pure source of singlet oxygen oxidizes catalases purified from different sources and singlet oxygen quenchers prevent oxidation. A second modification is observed in N. crassa catalase-1, in which the tyrosine that forms the fifth coordination bound to the heme iron makes a covalent bond with a vicinal cysteine, similarly to the tyrosine-histidine bonding found in Escherichia coli hydroperoxidase II. Molecular dynamics has been used to determine how H(2)O(2) reaches the enzyme active site and how products exit the protein. We found that the bottleneck of the major channel seems to disappear in water and is wide open in the presence of substrate. Amino acid residues exhibiting an increased residence time for H(2)O(2) are abundant at the protein surface and at the entrances to the major channel. The net effect of this is an increased H(2)O(2)/H(2)O ratio in the major channel. Once in the final section of this channel, H(2)O(2) is retained and tends to occupy specific sites while water molecules have a higher turnover rate and occupy different sites. Despite the intense study of catalases our knowledge of this enzyme is still limited and in need of new studies and different approaches.     

On the role of bicarbonate in peroxidations catalyzed by Cu,Zn superoxide dismutase

The interaction of Cu,ZnSOD with H2O2 generates an oxidant at the active site that can then cause either the inactivation of this enzyme or the oxidation of a variety of exogenous substrates. We show that the rate of inactivation, imposed by 10-mM H2O2 at 25 degrees C and pH 7.2, is not influenced by 10-mM HCO3-; whereas the oxidation of 2,2'-azino-bis-[3-ethylbenzothiazoline sulfonate] (ABTS=) is virtually completely dependent upon HCO3-. The reduction of the active site Cu(II) by H2O2, which precedes inactivation of the enzyme, occurred at the same rate in phosphate buffer with or without bicarbonate added. These results indicate that HCO3- does not play a role in facilitating the interaction of H2O2 with the active site copper, but they can be accommodated by the proposal that HCO3- is oxidized to HCO3*, which then diffuses from that site and causes the oxidation of substrates, such as ABTS=, that are too large to traverse the solvent access channel to the Cu(II).     

Pathways of H2 toward the active site of [NiFe]-hydrogenase

Hydrogenases catalyze the reversible oxidation of molecular hydrogen (H(2)), but little is known about the diffusion of H(2) toward the active site. Here we analyze pathways for H(2) permeation using molecular dynamics (MD) simulations in explicit solvent. Various MD simulation replicates were done, to improve the sampling of the system states. H(2) easily permeates hydrogenase in every simulation and it moves preferentially in channels. All H(2) molecules that reach the active site made their approach from the side of the Ni ion. H(2) is able to reach distances of <4 A from the active site, although after 6 A permeation is difficult. In this region we mutated Val-67 into alanine and perform new MD simulations. These simulations show an increase of H(2) inside the protein and at lower distances from the active site. This valine can be a control point in the H(2) access to the active center.     

Substrate flow in catalases deduced from the crystal structures of active site variants of HPII from Escherichia coli.

The active site of heme catalases is buried deep inside a structurally highly conserved homotetramer. Channels leading to the active site have been identified as potential routes for substrate flow and product release, although evidence in support of this model is limited. To investigate further the role of protein structure and molecular channels in catalysis, the crystal structures of four active site variants of catalase HPII from Escherichia coli (His128Ala, His128Asn, Asn201Ala, and Asn201His) have been determined at approximately 2.0-A resolution. The solvent organization shows major rearrangements with respect to native HPII, not only in the vicinity of the replaced residues but also in the main molecular channel leading to the heme distal pocket. In the two inactive His128 variants, continuous chains of hydrogen bonded water molecules extend from the molecular surface to the heme distal pocket filling the main channel. The differences in continuity of solvent molecules between the native and variant structures illustrate how sensitive the solvent matrix is to subtle changes in structure. It is hypothesized that the slightly larger H(2)O(2) passing through the channel of the native enzyme will promote the formation of a continuous chain of solvent and peroxide. The structure of the His128Asn variant complexed with hydrogen peroxide has also been determined at 2.3-A resolution, revealing the existence of hydrogen peroxide binding sites both in the heme distal pocket and in the main channel. Unexpectedly, the largest changes in protein structure resulting from peroxide binding are clustered on the heme proximal side and mainly involve residues in only two subunits, leading to a departure from the 222-point group symmetry of the native enzyme. An active role for channels in the selective flow of substrates through the catalase molecule is proposed as an integral feature of the catalytic mechanism. The Asn201His variant of HPII was found to contain unoxidized heme b in combination with the proximal side His-Tyr bond suggesting that the mechanistic pathways of the two reactions can be uncoupled.     

The ins and outs of cytochrome P450s

The active site of cytochromes P450 is situated deep inside the protein next to the heme cofactor. Consequently, enzyme specificity and kinetics can be influenced by how substrates pass through the protein to access the active site and how products egress from the active site. We previously analysed the channels between the active site and the protein surface in P450 crystal structures available in October 2003 [R.C. Wade, P.J. Winn, I. Schlichting, Sudarko, A survey of active site access channels in cytochromes P450, J. Inorg. Biochem. 98 (2004) 1175-1182]. Since then, 52 new P450 structures have been made available, including entries for ten isozymes for which structures were not previously available. We present an updated survey covering all P450 crystal structures available in March 2006. This survey shows channels not observed earlier in crystal structures, some of which were identified in previous molecular dynamics simulations. The crystal structures demonstrate how some of the channels can merge when the protein structure opens up resulting in a wide cleft to the active site, caused largely by movements of the F-G helix-loop-helix and the B-C loop. Significant differences were observed between the channels in the crystal structures of the mammalian and bacterial enzymes. The multiplicity of channels suggests possibilities for substrate channelling to and from the P450s.     

Role of protein and substrate dynamics in catalysis by Pseudomonas putida cytochrome P450cam

The role of protein structural flexibility and substrate dynamics in catalysis by cytochrome P450 enzymes is an area of current interest. We have addressed these in cytochrome P450(cam) (P450(cam)) and its Y96A mutant with camphor and its related compounds using fluorescence spectroscopy. Previously [Prasad et al. (2000) FEBS Lett. 477, 157-160], we provided experimental support to dynamic fluctuations in P450(cam), and substrate access into the active site region via the channel next to the flexible F-G helix-loop-helix segment. In the investigation described here, we show that the dynamic fluctuations in the enzyme are substrate dependent as reflected by tryptophan fluorescence quenching experiments. The orientation of tryptophan relative to heme (kappa(2)) for W42 obtained from time-resolved tryptophan fluorescence measurements show variation with type of substrate bound to P450(cam) suggesting regions distant from heme-binding site are affected by physicochemical and steric characteristics/protein-substrate interactions of P450(cam) active site. We monitored substrate dynamics in the active site region of P450(cam) by time-resolved substrate anisotropy measurements. The anisotropy decay of substrates bound to P450(cam) indicate that mobility of substrates is modulated by physicochemical and steric characteristics/protein-substrate interactions of local active site structure, and provides an understanding of factors controlling observed hydroxylated products for substrate bound P450(cam) complexes. The present study shows that P450(cam) local and peripheral structural flexibility and heterogeneity along with substrate mobility play an important role in regulating substrate binding orientation during catalysis and accommodating diverse range of substrates within P450(cam) heme pocket.     

Structural and dynamic properties of cytochrome P450 BM-3 in pure water and in a dimethylsulfoxide/water mixture

Solvent molecules play an important role for the structural and dynamical properties of proteins. A major focus of current protein engineering is the development of enzymes that are catalytically active in the presence of organic solvents. The monooxygenase P450 BM-3 is one of the best-studied enzymes and promising for industrial applications but with limited activity in the presence of organic solvents or cosolvents. To gain insights into the structural and dynamical properties of the heme domain of this enzyme in solution, molecular dynamics simulations in pure water and in a 14% DMSO/water mixture were performed. The results of the simulations show overall similar structural fluctuations in both solvent systems, with no indication of partial or global unfolding. In 14% DMSO, the regions comprising the helices E, F, and the EF loop (implicated in controlling the entry to the active site channel) undergo a large shift. Significant changes were also observed near the active site access channel at the residue R47. During the simulation, no DMSO molecule penetrated the active site. However, a significant accumulation of DMSO molecules close to the substrate-binding site and to the Flavin Mononucleotide (FMN) reductase domain interface was observed.     

Water as a structural element in a channel: gating in the Kcsa channel, and implications for voltage-gated ion channels

Water is becoming understood as a structural element in proteins. Here we are concerned with one particular type of protein, ion channels. The S. Lividans KcsA K(+) channel, the X-ray structure of which is known, is gated by protons (i.e, by a drop in pH). Ab initio calculations suggest that an H(5)O(2) group, partially charged, connects the E118 residues in the gating region, when the four residues have a -2 net charge, but that the hydrogen bonding is not strong enough to do this when the charge becomes -1. The H(5)O(2) group would block the channel, in the -2 state, and prevent motion of the four transmembrane (TM) segments of the protein, by binding them. With the weaker bond in the -1 state, the TM segments would be able to separate (as they have been found to do experimentally, opening the channel. Voltage gated channels have four additional TM segments for each of the four domains of the channel protein. These appear to allow motion of protons; in fact there is evidence that the initial step in gating must be the transfer of a proton. We have earlier shown that the transfer of a single proton between two methylamines under the influence of a field is possible, as proton tunneling. Subsequent steps are hypothesized to result from four proton transfer cascades of about three protons each, triggered by the initial proton transfer. We suggest that the extra 4 TM segments of the voltage gated channel act as a voltage to proton-current transducer. Water, held by hydrogen bonds, is also suggested as the source of the accessibility data found with MTS reagents, based largely on simulations, our earlier Monte Carlo simulations as well as molecular dynamics studies reported by others. These waters may also play a structural role in the protein.     

Mechanisms of hydrogen peroxide-induced calcium dysregulation in PC12 cells

The mechanisms of H(2)O(2)-induced elevated calcium baselines in PC12 cells were investigated in the present study by using fura-2-fluorescent image analysis. The results showed that the calcium comes from both intracellular and extracellular sources. Although the major intracellular source was mitochondria, only the extracellular calcium influx was responsible for the sustained post-H(2)O(2)-exposure increases. This calcium influx was partially blocked by calcium channel antagonists [verapamil (L-type) or mibefradil (nonselective)] and was more effectively blocked by the sodium channel antagonist, tetrodotoxin (TTX). Membrane depolarization following H(2)O(2) exposure contributed to the opening of the ion channels. The H(2)O(2)-induced calcium influx was blocked by TTX even in a sodium-free buffer, indicating that calcium directly fluxed through sodium channels. Sodium-calcium exchangers (NCX) on the plasma membrane did not play a role, because use of a specific reverse mode NCX inhibitor, No. 7943, was ineffective in blocking the influx. The H(2)O(2)-induced calcium influx was mimicked by using a thiol-selective oxidizing reagent, 2', 2'-dithiodipyridine, and in both situations, the calcium levels were completely reversed by a thiol-selective reducing reagent, dithiothreitol. Our results indicated that mechanisms of oxidant-induced elevated calcium baselines in PC12 cells involved calcium influx through sodium and calcium channels that may be directly or indirectly attributed to thiol oxidation.     

Structural dynamics of ligand diffusion in the protein matrix: A study on a new myoglobin mutant Y(B10) Q(E7) R(E10)

A triple mutant of sperm whale myoglobin (Mb) [Leu(B10) --> Tyr, His(E7) --> Gln, and Thr(E10) --> Arg, called Mb-YQR], investigated by stopped-flow, laser photolysis, crystallography, and molecular dynamics (MD) simulations, proved to be quite unusual. Rebinding of photodissociated NO, O2, and CO from within the protein (in a "geminate" mode) allows us to reach general conclusions about dynamics and cavities in proteins. The 3D structure of oxy Mb-YQR shows that bound O2 makes two H-bonds with Tyr(B10)29 and Gln(E7)64; on deoxygenation, these two residues move toward the space occupied by O2. The bimolecular rate constant for NO binding is the same as for wild-type, but those for CO and O2 binding are reduced 10-fold. While there is no geminate recombination with O2 and CO, geminate rebinding of NO displays an unusually large and very slow component, which is pretty much abolished in the presence of xenon. These results and MD simulations suggest that the ligand migrates in the protein matrix to a major "secondary site," located beneath Tyr(B10)29 and accessible via the motion of Ile(G8)107; this site is different from the "primary site" identified by others who investigated the photolyzed state of wild-type Mb by crystallography. Our hypothesis may rationalize the O2 binding properties of Mb-YQR, and more generally to propose a mechanism of control of ligand binding and dissociation in hemeproteins based on the dynamics of side chains that may (or may not) allow access to and direct temporary sequestration of the dissociated ligand in a docking site within the protein. This interpretation suggests that very fast (picosecond) fluctuations of amino acid side chains may play a crucial role in controlling O2 delivery to tissue at a rate compatible with physiology.     

Anatomy of enzyme channels

Background

Enzyme active sites can be connected to the exterior environment by one or more channels passing through the protein. Despite our current knowledge of enzyme structure and function, surprisingly little is known about how often channels are present or about any structural features such channels may have in common.

Results

Here, we analyze the long channels (i.e. >15 Å) leading to the active sites of 4,306 enzyme structures. We find that over 64% of enzymes contain two or more long channels, their typical length being 28 Å. We show that amino acid compositions of the channel significantly differ both to the composition of the active site, surface and interior of the protein.

Conclusions

The majority of enzymes have buried active sites accessible via a network of access channels. This indicates that enzymes tend to have buried active sites, with channels controlling access to, and egress from, them, and that suggests channels may play a key role in helping determine enzyme substrate.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0379-x) contains supplementary material, which is available to authorized users.     

Ca(2+)-activated K+ channel inhibition by reactive oxygen species

We studied the effect of H(2)O(2) on the gating behavior of large-conductance Ca(2+)-sensitive voltage-dependent K(+) (K(V,Ca)) channels. We recorded potassium currents from single skeletal muscle channels incorporated into bilayers or using macropatches of Xenopus laevis oocytes membranes expressing the human Slowpoke (hSlo) alpha-subunit. Exposure of the intracellular side of K(V,Ca) channels to H(2)O(2) (4-23 mM) leads to a time-dependent decrease of the open probability (P(o)) without affecting the unitary conductance. H(2)O(2) did not affect channel activity when added to the extracellular side. These results provide evidence for an intracellular site(s) of H(2)O(2) action. Desferrioxamine (60 microM) and cysteine (1 mM) completely inhibited the effect of H(2)O(2), indicating that the decrease in P(o) was mediated by hydroxyl radicals. The reducing agent dithiothreitol (DTT) could not fully reverse the effect of H(2)O(2). However, DTT did completely reverse the decrease in P(o) induced by the oxidizing agent 5,5'-dithio-bis-(2-nitrobenzoic acid). The incomplete recovery of K(V,Ca) channel activity promoted by DTT suggests that H(2)O(2) treatment must be modifying other amino acid residues, e.g., as methionine or tryptophan, besides cysteine. Noise analysis of macroscopic currents in Xenopus oocytes expressing hSlo channels showed that H(2)O(2) induced a decrease in current mediated by a decrease both in the number of active channels and P(o).     

Coordination of membrane excitability through a GIRK1 signaling complex in the atria

Control of heart rate is a complex process that integrates the function of multiple G protein-coupled receptors and ion channels. Among them, the G protein-regulated inwardly rectifying K+ (GIRK or KACh) channels of sinoatrial node and atria play a major role in beat-to-beat regulation of the heart rate. The atrial KACh channels are heterotetrameric proteins that consist of two pore-forming subunits, GIRK1 and GIRK4. Following m2-muscarinic acetylcholine receptor (M2R) stimulation, KACh channel activation is conferred by the direct binding of G protein betagamma subunits (Gbetagamma) to the channel. Here we show that atrial KACh channels are assembled in a signaling complex with Gbetagamma, G protein-coupled receptor kinase, cyclic adenosine monophosphate-dependent protein kinase, two protein phosphatases, PP1 and PP2A, receptor for activated C kinase 1, and actin. This complex would enable the KACh channels to rapidly integrate beta-adrenergic and M2R signaling in the membrane, and it provides insight into general principles governing spatial integration of different transduction pathways. Furthermore, the same complex might recruit protein kinase C (PKC) to the KACh channel following alpha-adrenergic receptor stimulation. Our electro-physiological recordings from single atrial KACh channels revealed a potent inhibition of Gbetagamma-induced channel activity by PKC, thus validating the physiological significance of the observed complex as interconnecting site where signaling molecules congregate to execute a coordinated control of membrane excitability.     

A survey of active site access channels in cytochromes P450

In cytochrome P450s, the active site is situated deep inside the protein next to the heme cofactor, and is often completely isolated from the surrounding solvent. To identify routes by which substrates may enter into and products exit from the active site, random expulsion molecular dynamics simulations were performed for three cytochrome P450s: CYP101, CYP102A1 and CYP107A1 [J. Mol. Biol. 303 (2000) 797; Proc. Natl. Acad. Sci. USA 99 (2002) 5361]. Amongst the different pathways identified, one pathway was found to be common to all three cytochrome P450s although the mechanism of ligand passage along it was different in each case and apparently adapted to the substrate specificity of the enzyme. Recently, a number of new crystal structures of cytochrome P450s have been solved. Here, we analyse the open channels leading to the active site that these structures reveal. We find that in addition to showing the common pathway, they provide experimental evidence for the existence of three additional channels that were identified by simulation. We also discuss how the location of xenon binding sites in CYP101 suggests a role for one of the pathways identified by molecular dynamics simulations as a route for gaseous species, such as oxygen, to access the active site.     

The dynamics of camphor in the cytochrome P450 CYP101D2

The recent crystal structures of CYP101D2, a cytochrome P450 protein from the oligotrophic bacterium Novosphingobium aromaticivorans DSM12444 revealed that both the native (substrate‐free) and camphor‐soaked forms have open conformations. Furthermore, two other potential camphor‐binding sites were also identified from electron densities in the camphor‐soaked structure, one being located in the access channel and the other in a cavity on the surface near the F‐helix side of the F‐G loop termed the substrate recognition site. These latter sites may be key intermediate positions on the pathway for substrate access to or product egress from the active site. Here, we show via the use of unbiased atomistic molecular dynamics simulations that despite the open conformation of the native and camphor‐bound crystal structures, the underlying dynamics of CYP101D2 appear to be very similar to other CYP proteins. Simulations of the native structure demonstrated that the protein is capable of sampling many different conformational substates. At the same time, simulations with the camphor positioned at various locations within the access channel or recognition site show that movement towards the active site or towards bulk solvent can readily occur on a short timescale, thus confirming many previously reported in silico studies using steered molecular dynamics. The simulations also demonstrate how the fluctuations of an aromatic gate appear to control access to the active site. Finally, comparison of camphor‐bound simulations with the native simulations suggests that the fluctuations can be of similar level and thus are more representative of the conformational selection model rather than induced fit.     

Structure and in silico simulations of a cold-active esterase reveals its prime cold-adaptation mechanism

Here we determined the structure of a cold active family IV esterase (EstN7) cloned from Bacillus cohnii strain N1. EstN7 is a dimer with a classical α/β hydrolase fold. It has an acidic surface that is thought to play a role in cold-adaption by retaining solvation under changed water solvent entropy at lower temperatures. The conformation of the functionally important cap region is significantly different to EstN7''s closest relatives, forming a bridge-like structure with reduced helical content providing greater access to the active site through more than one substrate access tunnel. However, dynamics do not appear to play a major role in cold adaption. Molecular dynamics at different temperatures, rigidity analysis, normal mode analysis and geometric simulations of motion confirm the flexibility of the cap region but suggest that the rest of the protein is largely rigid. Rigidity analysis indicates the distribution of hydrophobic tethers is appropriate to colder conditions, where the hydrophobic effect is weaker than in mesophilic conditions due to reduced water entropy. Thus, it is likely that increased substrate accessibility and tolerance to changes in water entropy are important for of EstN7''s cold adaptation rather than changes in dynamics.     

Dynamic mechanism of proton transfer in mannitol 2-dehydrogenase from Pseudomonas fluorescens: mobile GLU292 controls proton relay through a water channel that connects the active site with bulk solvent

The active site of mannitol 2-dehydrogenase from Pseudomonas fluorescens (PfM2DH) is connected with bulk solvent through a narrow protein channel that shows structural resemblance to proton channels utilized by redox-driven proton pumps. A key element of the PfM2DH channel is the "mobile" Glu(292), which was seen crystallographically to adopt distinct positions up and down the channel. It was suggested that the "down → up" conformational change of Glu(292) could play a proton relay function in enzymatic catalysis, through direct proton shuttling by the Glu or because the channel is opened for water molecules forming a chain along which the protons flow. We report evidence from site-directed mutagenesis (Glu(292) → Ala) substantiated by data from molecular dynamics simulations that support a role for Glu(292) as a gate in a water chain (von Grotthuss-type) mechanism of proton translocation. Occupancy of the up and down position of Glu(292) is influenced by the bonding and charge state of the catalytic acid base Lys(295), suggesting that channel opening/closing motions of the Glu are synchronized to the reaction progress. Removal of gatekeeper control in the E292A mutant resulted in a selective, up to 120-fold slowing down of microscopic steps immediately preceding catalytic oxidation of mannitol, consistent with the notion that formation of the productive enzyme-NAD(+)-mannitol complex is promoted by a corresponding position change of Glu(292), which at physiological pH is associated with obligatory deprotonation of Lys(295) to solvent. These results underscore the important role of conformational dynamics in the proton transfer steps of alcohol dehydrogenase catalysis.