Protein side-chain motion and hydration in proton-transfer pathways. Results for cytochrome p450cam

In the footsteps of an inhibitor unbinding from the active site of human carbonic anhydrase II

The crystal structure of human carbonic anhydrase (HCA) II bound to an inhibitor molecule, 6-hydroxy-2-thioxocoumarin (FC5), shows FC5 to be located in a hydrophobic pocket at the active site. The present work employs classical molecular dynamics (MD) simulation to follow the FC5 molecule for 1 μs as it unbinds from its binding location, adopts the path of substrate/product diffusion (path 1) to leave the active site at around 75 ns. It is then found to undergo repeated binding and unbinding at different locations on the surface of the enzyme in water. Several transient excursions through different regions of the enzyme are also observed prior to its exit from the active site. These transient paths are combined with functionally relevant cavities/channels to enlist five additional pathways (path 2-6). Pathways 1-6 are subsequently explored using steered MD and umbrella sampling simulations. A free energy barrier of 0.969 kcal mol^-1 is encountered along path 1, while barriers in the range of 0.57-2.84 kcal mol^-1 are obtained along paths 2, 4 and 5. We also analyze in detail the interaction between FC5 and the enzyme along each path as the former leaves the active site of HCA II. Our results indicate path 1 to be the major exit pathway for FC5, although competing contributions may also come from the paths 2, 4 and 5.Communicated by Ramaswamy H. Sarma.

C-Graphs Tool with Graphical User Interface to Dissect Conserved Hydrogen-Bond Networks: Applications to Visual Rhodopsins

Dynamic hydrogen-bond networks provide proteins with structural plasticity required to translate signals such as ligand binding into a cellular response or to transport ions and larger solutes across membranes and, thus, are of central interest to understand protein reaction mechanisms. Here, we present C-Graphs, an efficient tool with graphical user interface that analyzes data sets of static protein structures or of independent numerical simulations to identify conserved, vs unique, hydrogen bonds and hydrogen-bond networks. For static structures, which may belong to the same protein or to proteins with different sequences, C-Graphs uses a clustering algorithm to identify sites of the hydrogen-bond network where waters are conserved among the structures. Using C-Graphs, we identify an internal protein-water hydrogen-bond network common to static structures of visual rhodopsins and adenosine A2A G protein-coupled receptors (GPCRs). Molecular dynamics simulations of a visual rhodopsin indicate that the conserved hydrogen-bond network from static structure can recruit dynamic hydrogen bonds and extend throughout most of the receptor. We release with this work the code for C-Graphs and its graphical user interface.

QM/MM Simulations of Enzymatic Hydrolysis of Cellulose: Probing the Viability of an Endocyclic Mechanism for an Inverting Cellulase

Glycoside hydrolases (GH) cleave carbohydrate glycosidic bonds and play pivotal roles in living organisms and in many industrial processes. Unlike acid-catalyzed hydrolysis of carbohydrates in solution, which can occur either via cyclic or acyclic oxocarbenium-like transition states, it is widely accepted that GH-catalyzed hydrolysis proceeds via a general acid mechanism involving a cyclic oxocarbenium-like transition state with protonation of the glycosidic oxygen. The GH45 subfamily C inverting endoglucanase from Phanerochaete chrysosporium (PcCel45A) defies the classical inverting mechanism as its crystal structure conspicuously lacks a general Asp or Glu base residue. Instead, PcCel45A has an Asn residue, a notoriously weak base in solution, as one of its catalytic residues at position 92. Moreover, unlike other inverting GHs, the relative position of the catalytic residues in PcCel45A impairs the proton abstraction from the nucleophilic water that attacks the anomeric carbon, a key step in the classical mechanism. Here, we investigate the viability of an endocyclic mechanism for PcCel45A using hybrid quantum mechanics/molecular mechanics (QM/MM) simulations, with the QM region treated with the self-consistent-charge density-functional tight-binding level of theory. In this mechanism, an acyclic oxocarbenium-like transition state is stabilized leading to the opening of the glucopyranose ring and formation of an unstable acyclic hemiacetal that can be readily decomposed into hydrolysis product. In silico characterization of the Michaelis complex shows that PcCel45A significantly restrains the sugar ring to the ⁴C₁ chair conformation at the −1 subsite of the substrate binding cleft, in contrast to the classical exocyclic mechanism in which ring puckering is critical. We also show that PcCel45A provides an environment where the catalytic Asn92 residue in its standard amide form participates in a cooperative hydrogen bond network resulting in its increased nucleophilicity due to an increased negative charge on the oxygen atom. Our results for PcCel45A suggest that carbohydrate hydrolysis catalyzed by GHs may take an alternative route from the classical mechanism.

Thermodynamic and Kinetic Transitions of Liquids in Nanoconfinement

Core principles of chemistry are ubiquitously invoked to shed light on the nature of molecular level interactions in nanoconfined fluids, which play a pivotal role in a wide range of processes in geochemistry, biology, and engineering. A detailed understanding of the physicochemical processes involved in the flow, structural transitions, and freezing or melting behavior of fluids confined within nanometer-sized pores of solid materials is thus of enormous importance for both basic research and technological applications.This Account provides a perspective on new insights into the thermodynamic and kinetic transitions of nanoconfined fluids in their stable and metastable forms. After briefly introducing the unique properties of mesoporous silicas from the SBA, MCM, and FDU families that serve as the confinement matrices, combining highly ordered single and bimodal mesopore architectures with tunable pore sizes in the ∼2-15 nm range and narrow size distributions, recent studies on melting/freezing behavior of water confined in these host matrices are reviewed. While differential scanning calorimetry (DSC) reveals a linear relationship between melting point depression and pore size (independent of the pore shape), as predicted by the Gibbs-Thomson relation, variable temperature ²H wide-line nuclear magnetic resonance (NMR) spectroscopy studies confirm the core-shell model of water and give evidence for a layer-by-layer freezing mechanism, which gives rise to an apparent fragile-to-strong transition in the solidification dynamics.In contrast to the freezing/melting behavior of water, the effect of nanoconfinement on the glass transition of supercooled liquids is nonuniversal and the glass transition temperature T_g can either increase or decrease with the dimensionality and extent of confinement. This nonuniversal behavior is exemplified by the two glass-forming molecular liquids, glycerol and ortho-terphenyl (OTP). While glycerol shows an increase in T_g and a pronounced slowdown of the rotational dynamics of the constituent molecules due to a change in the molecular packing between the bulk and the confined liquid, OTP displays a linear and confining-media-dependent depression of T_g with increased confinement that is strongly influenced by the pore-liquid interface characteristics.This Account concludes with a focus on recent experimental evidence of extreme spatial and dynamical heterogeneity in both freezing and glass transition processes. This discovery was enabled by the unique mesoporous structures of SBA-16 and FDU-5, possessing bimodal architectures with two interconnected pore types of different size and shape (spherical and cylindrical). For the very first time, two melting points for water and two glass transitions for supercooled OTP, corresponding to a specific pore type, were observed. Collectively, these observations strongly suggest a close mechanistic connection between the local fluctuations in the structure and dynamics of nanoconfined liquids. While the findings reviewed in this Account provide new insights into thermodynamic and kinetic transitions of fluids, there remain many unanswered questions regarding the effects of nanoconfinement on the fundamental properties of fluids, which offer exciting future opportunities in chemical research.

From an organic ligand to a metal–organic coordination polymer, and to a metal–organic coordination polymer–cocrystal composite: a continuous promotion of the proton conductivity of crystalline materials

A new strategy was proposed to increase proton conductivities in metal–organic coordination polymers (MOCPs) commencing from organic ligands, i.e. coordination inducement and MOCP–cocrystal composite formation.

What Are We Missing by Not Measuring the Net Charge of Proteins?

The net electrostatic charge (Z) of a folded protein in solution represents a bird's eye view of its surface potentials-including contributions from tightly bound metal, solvent, buffer, and cosolvent ions-and remains one of its most enigmatic properties. Few tools are available to the average biochemist to rapidly and accurately measure Z at pH≠pI. Tools that have been developed more recently seem to go unnoticed. Most scientists are content with this void and estimate the net charge of a protein from its amino acid sequence, using textbook values of pK_a . Thus, Z remains unmeasured for nearly all folded proteins at pH≠pI. When marveling at all that has been learned from accurately measuring the other fundamental property of a protein-its mass-one wonders: what are we missing by not measuring the net charge of folded, solvated proteins? A few big questions immediately emerge in bioinorganic chemistry. When a single electron is transferred to a metalloprotein, does the net charge of the protein change by approximately one elementary unit of charge or does charge regulation dominate, that is, do the pK_a values of most ionizable residues (or just a few residues) adjust in response to (or in concert with) electron transfer? Would the free energy of charge regulation (ΔΔG_z ) account for most of the outer sphere reorganization energy associated with electron transfer? Or would ΔΔG_z contribute more to the redox potential? And what about metal binding itself? When an apo-metalloprotein, bearing minimal net negative charge (e.g., Z=-2.0) binds one or more metal cations, is the net charge abolished or inverted to positive? Or do metalloproteins regulate net charge when coordinating metal ions? The author's group has recently dusted off a relatively obscure tool-the "protein charge ladder"-and used it to begin to answer these basic questions.

Hydrogen-Bonded Network and Water Dynamics in the D-channel of Cytochrome c Oxidase

Proton transfer in cytochrome c oxidase (CcO) from the cellular inside to the binuclear redox centre as well as proton pumping through the membrane takes place through proton entrance via two distinct pathways, the D- and K-channel. Both channels show a dependence of their hydration level on the protonation states of their key residues, K362 for the K-channel, and E286 or D132 for the D-channel. In the oxidative half of CcO's catalytic cycle the D-channel is the proton-conducting path. For this channel, an interplay of protonation state of the D-channel residues with the water and hydrogen-bond dynamics has been observed in molecular dynamics simulations of the CcO protein, embedded in a lipid bi-layer, modelled in different protonation states. Protonation of residue E286 at the end of the D-channel results in a hydrogen-bonded network pointing from E286 to N139, that is against proton transport, and favouring N139 conformations which correspond to a closed asparagine gate (formed by residues N121 and N139). Consequently, the hydration level is lower than with unprotonated E286. In those models, the Asn gate is predominantly open, allowing water molecules to pass and thus increase the hydration level. The hydrogen-bonded network in these states exhibits longer life times of the Asn residues with water than other models and shows the D-channel to be traversable from the entrance, D132, to exit, E286. The D-channel can thus be regarded as auto-regulated with respect to proton transport, allowing proton passage only when required, that is the proton is located at the lower part of the D-channel (D132 to Asn gate) and not at the exit (E286).

Distal Proton Shuttle Mechanism of Ribosome Catalysed Peptide Bond Formation-A Theoretical Study

In this work, we have investigated a novel distal proton shuttle mechanism of ribosome catalyzed peptide bond formation reaction. The reaction was found to follow a two-step mechanism. A distal water molecule located about 6.1 Å away from the attacking amine plays as a proton acceptor and results in a charge-separated intermediate that is stabilized by the N terminus of L27 and the A-site A76 5′-phosphate. The ribose A2451 bridges the proton shuttle pathway, thus plays critical role in the reaction. The calculated 27.64 kcal·mol⁻¹ free energy barrier of the distal proton shuttle mechanism is lower than that of eight-membered ring transition state. The distal proton shuttle mechanism studied in this work can provide new insights into the important biological peptide synthesis process.

Protonation of the Hydroperoxo Intermediate of Cytochrome P450 2B4 Is Slower in the Presence of Cytochrome P450 Reductase Than in the Presence of Cytochrome b5

Microsomal cytochromes P450 (P450) require two electrons and two protons for the oxidation of substrates. Although the two electrons can be provided by cytochrome P450 reductase, the second electron can also be donated by cytochrome b5 (b5). The steady-state activity of P450 2B4 is increased up to 10-fold by b5. To improve our understanding of the molecular basis of the stimulatory effect of b5 and to test the hypothesis that b5 stimulates catalysis by more rapid protonation of the anionic ferric hydroperoxo heme intermediate of P450 (Fe³⁺OOH)^- and subsequent formation of the active oxidizing species (Fe⁺⁴═O POR^•+), we have freeze-quenched the reaction mixture during a single turnover following reduction of oxyferrous P450 2B4 by each of its redox partners, b5 and P450 reductase. The electron paramagnetic resonance spectra of the freeze-quenched reaction mixtures lacked evidence of a hydroperoxo intermediate when b5 was the reductant presumably because hydroperoxo protonation and catalysis occurred within the dead time of the instrument. However, when P450 reductase was the reductant, a hydroperoxo P450 intermediate was observed. The effect of b5 on the enzymatic efficiency in D₂O and the kinetic solvent isotope effect under steady-state conditions are both consistent with the ability of b5 to promote rapid protonation of the hydroperoxo species and more efficient catalysis. In summary, by binding to the proximal surface of P450, b5 stimulates the activity of P450 2B4 by enhancing the rate of protonation of the hydroperoxo intermediate and formation of Compound I, the active oxidizing species, which allows less time for side product formation.

The Cytochrome P450-Catalyzed Oxidative Rearrangement in the Final Step of Pentalenolactone Biosynthesis: Substrate Structure Determines Mechanism

The final step in the biosynthesis of the sesquiterpenoid antibiotic pentalenolactone (1) is the highly unusual cytochrome P450-catalyzed, oxidative rearrangement of pentalenolactone F (2), involving the transient generation and rearrangement of a neopentyl cation. In Streptomyces arenae this reaction is catalyzed by CYP161C2 (PntM), with highly conserved orthologs being present in at least 10 other Actinomycetes. Crystal structures of substrate-free PntM, as well as PntM with bound substrate 2, product 1, and substrate analogue 6,7-dihydropentalenolactone F (7) revealed interactions of bound ligand with three residues, F232, M77, and M81 that are unique to PntM and its orthologs and absent from essentially all other P450s. Site-directed mutagenesis, ligand-binding measurements, steady-state kinetics, and reaction product profiles established there is no special stabilization of reactive cationic intermediates by these side chains. Reduced substrate analogue 7 did not undergo either oxidative rearrangement or simple hydroxylation, suggesting that the C1 carbocation is not anchimerically stabilized by the 6,7-double bond of 2. The crystal structures also revealed plausible proton relay networks likely involved in the generation of the key characteristic P450 oxidizing species, Compound I, and in mediating stereospecific deprotonation of H-3re of the substrate. We conclude that the unusual carbocation intermediate results from outer shell electron transfer from the transiently generated C1 radical to the tightly paired heme-•Fe(3+)-OH radical species. The oxidative electron transfer is kinetically dominant as a result of the unusually strong steric barrier to oxygen rebound to the neopentyl center C-1si, which is flanked on each neighboring carbon by syn-axial substituents.

Pore hydration states of KcsA potassium channels in membranes

Water-filled hydrophobic cavities in channel proteins serve as gateways for transfer of ions across membranes, but their properties are largely unknown. We determined water distributions along the conduction pores in two tetrameric channels embedded in lipid bilayers using neutron diffraction: potassium channel KcsA and the transmembrane domain of M2 protein of influenza A virus. For the KcsA channel in the closed state, the distribution of water is peaked in the middle of the membrane, showing water in the central cavity adjacent to the selectivity filter. This water is displaced by the channel blocker tetrabutyl-ammonium. The amount of water associated with the channel was quantified, using neutron diffraction and solid state NMR. In contrast, the M2 proton channel shows a V-shaped water profile across the membrane, with a narrow constriction at the center, like the hourglass shape of its internal surface. These two types of water distribution are therefore very different in their connectivity to the bulk water. The water and protein profiles determined here provide important evidence concerning conformation and hydration of channels in membranes and the potential role of pore hydration in channel gating.

Differences and comparisons of the properties and reactivities of iron(III)-hydroperoxo complexes with saturated coordination sphere

Heme and nonheme monoxygenases and dioxygenases catalyze important oxygen atom transfer reactions to substrates in the body. It is now well established that the cytochrome P450 enzymes react through the formation of a high-valent iron(IV)-oxo heme cation radical. Its precursor in the catalytic cycle, the iron(III)-hydroperoxo complex, was tested for catalytic activity and found to be a sluggish oxidant of hydroxylation, epoxidation and sulfoxidation reactions. In a recent twist of events, evidence has emerged of several nonheme iron(III)-hydroperoxo complexes that appear to react with substrates via oxygen atom transfer processes. Although it was not clear from these studies whether the iron(III)-hydroperoxo reacted directly with substrates or that an initial O-O bond cleavage preceded the reaction. Clearly, the catalytic activity of heme and nonheme iron(III)-hydroperoxo complexes is substantially different, but the origins of this are still poorly understood and warrant a detailed analysis. In this work, an extensive computational analysis of aromatic hydroxylation by biomimetic nonheme and heme iron systems is presented, starting from an iron(III)-hydroperoxo complex with pentadentate ligand system (L5(2)). Direct C-O bond formation by an iron(III)-hydroperoxo complex is investigated, as well as the initial heterolytic and homolytic bond cleavage of the hydroperoxo group. The calculations show that [(L5(2))Fe(III)(OOH)](2+) should be able to initiate an aromatic hydroxylation process, although a low-energy homolytic cleavage pathway is only slightly higher in energy. A detailed valence bond and thermochemical analysis rationalizes the differences in chemical reactivity of heme and nonheme iron(III)-hydroperoxo and show that the main reason for this particular nonheme complex to be reactive comes from the fact that they homolytically split the O-O bond, whereas a heterolytic O-O bond breaking in heme iron(III)-hydroperoxo is found.

Time-averaged distributions of solute and solvent motions: exploring proton wires of GFP and PfM2DH

Proton translocation pathways of selected variants of the green fluorescent protein (GFP) and Pseudomonas fluorescens mannitol 2-dehydrogenase (PfM2DH) were investigated via an explicit solvent molecular dynamics-based analysis protocol that allows for direct quantitative relationship between a crystal structure and its time-averaged solute-solvent structure obtained from simulation. Our study of GFP is in good agreement with previous research suggesting that the proton released from the chromophore upon photoexcitation can diffuse through an extended internal hydrogen bonding network that allows for the proton to exit to bulk or be recaptured by the anionic chromophore. Conversely for PfM2DH, we identified the most probable ionization states of key residues along the proton escape channel from the catalytic site to bulk solvent, wherein the solute and high-density solvent crystal structures of binary and ternary complexes were properly reproduced. Furthermore, we proposed a plausible mechanism for this proton translocation process that is consistent with the state-dependent structural shifts observed in our analysis. The time-averaged structures generated from our analyses facilitate validation of MD simulation results and provide a comprehensive profile of the dynamic all-occupancy solvation network within and around a flexible solute, from which detailed hydrogen-bonding networks can be inferred. In this way, potential drawbacks arising from the elucidation of these networks by examination of static crystal structures or via alternate rigid-protein solvation analysis procedures can be overcome. Complementary studies aimed at the effective use of our methodology for alternate implementations (e.g., ligand design) are currently underway.

An application of QM/MM simulation: the second protonation of cytochrome P450

The multiscale model strategy, hybrid quantum mechanics and molecular mechanics (QM/MM), has become more and more prevalent in the theoretical study of enzymatic reactions. It combines both the efficiency of the Newtonian molecular calculations and the accuracy of the quantum mechanical methods. Simulation using QM/MM multiscale model may be one of the most promising approaches that could further narrow the gap between the theoretical models and the real problems. It is capable of dealing with not only the conformational changes of biomacromolecules, but also the catalytic reactions. Herein, we reviewed some of our recent work to demonstrate the application of the QM/MM simulations in exploring the enzymatic reactions.

pH-dependent structural conformations of B-phycoerythrin from Porphyridium cruentum

B-phycoerythrin from the red alga Porphyridium cruentum was crystallized using the technique of capillary counter-diffusion. Crystals belonging to the space group R3 with almost identical unit cell constants and diffracting to 1.85 and 1.70 Å were obtained at pH values of 5 and 8, respectively. The most important difference between structures is the presence of the residue His88α in two different conformations at pH 8. This residue is placed next to the chromophore phycoerythrobilin PEB82α and the new conformation results in the relocation of the hydrogen-bond network and hydration around PEB82α, which probably contributes to the observed pH dependence of the optical spectrum associated with this chromophore. Comparison with the structures of B-phycoerythrin from other red algae shows differences in the conformation of the A-ring of the chromophore PEB139α. This conformational difference in B-phycoerythrin from P. cruentum enables the formation of several hydrogen bonds that connect PEB139α with the chromophore PEB158β at the (αβ)(3) hexamer association interface. The possible influence of these structural differences on the optical spectrum and the ability of the protein to perform energy transfer are discussed, with the two pH-dependent conformations of His88α and PEB82α being proposed as representing critical structural features that are correlated with the pH dependence of the optical spectrum and transient optical states during energy transfer.

Modeling the structure and proton transfer pathways of the mutant His-107-Tyr of human carbonic anhydrase II

We present molecular modeling of the structure and possible proton transfer pathways from the surface of the protein to the zinc-bound water molecule in the active site of the mutant His-107-Tyr of human carbonic anhydrase II (HCAII). No high-resolution structure or crystal structure is available till now for this particular mutant due to its lack of stability at physiological temperature. Our analysis utilizes as starting point a series of structures derived from high-resolution crystal structure of the wild type protein. While many of the structures investigated do not reveal a complete path between the zinc bound water and His-64, several others do indicate the presence of a transient connection even when His-64 is present in its outward conformation. Mutation at the residue 107 also reveals the formation of a new path into the active site. Competing contributions from His-64 sidechain rotation from its outward conformation are also evaluated in terms of optimal path analysis. No indication of a lower catalytic efficiency of the mutant is evident from our results under the condition of thermal stability of the mutant.

Energetics and dynamics of proton transfer reactions along short water wires

Proton transfer (pT) reactions in biochemical processes are often mediated by chains of hydrogen-bonded water molecules. We use hybrid density functional calculations to study pT along quasi one-dimensional water arrays that connect an imidazolium-imidazole proton donor-acceptor pair. We characterize the structures of intermediates and transition states, the energetics, and the dynamics of the pT reactions, including vibrational contributions to kinetic isotope effects. In molecular dynamics simulations of pT transition paths, we find that for short water chains with four water molecules, the pT reactions are semi-concerted. The formation of a high-energy hydronium intermediate next to the proton-donating group is avoided by a simultaneous transfer of a proton from the donor to the first water molecule, and from the first water molecule into the water chain. Lowering the dielectric constant of the environment and increasing the water chain length both reduce the barrier for pT. We study the effect of the driving force on the energetics of the pT reaction by changing the proton affinity of the donor and acceptor groups through halogen and methyl substitutions. We find that the barrier of the pT reaction depends linearly on the proton affinity of the donor but is nearly independent of the proton affinity of the acceptor, corresponding to Brønsted slopes of one and zero, respectively.

The rate-limiting step in P450 hydroxylation of hydrocarbons a direct comparison of the "somersault" versus the "consensus" mechanism involving compound I

Model theoretical quantum mechanical (QM) calculations are described for the P-450 hydroxylation of methane, isobutane, and camphor that compare the concerted somersault H-abstraction mechanism with the oxidation step involving Cpd I. Special emphasis has been placed on maintaining a balanced basis set in the oxidation step. QM calculations, employing the 6-311+G(d,p) basis set on the Fe atom and all of the key surrounding atoms involved in the C-H abstraction step, reaffirm a mechanism involving rearrangement of the iron hydroperoxide group (FeO-OH --> FeO...HO(*)) in concert with hydrogen abstraction from the C-H bond of the substrate by the incipient bound hydroxyl radical HO(*). The barrier for the somersault rearrangement of model Cpd 0 (FeO-OH) is calculated to be 21.4 kcal/mol in the absence of substrate. The overall activation energy for the oxidation of camphor involving the somersault motion of the FeO-OH group of P450 model porphyrin iron(III) hydroperoxide [Por(SH)Fe(III)-OOH(-)] --> [Por(SH)Fe(III)-O....HO(-)] in concert with hydrogen abstraction is DeltaE(++) = 12.4 kcal/mol. The corresponding abstraction of the hydrogen atom from the C-H bond of camphor by Cpd I has an activation barrier of 17.6 kcal/mol. Arguments are presented that the somersault rearrangement is induced by steric compression at the active site. Kinetic isotope effect data are discussed that provides compelling evidence for a rate-limiting step involving C-H bond cleavage.

Coupling and uncoupling mechanisms in the methoxythreonine mutant of cytochrome P450cam: a quantum mechanical/molecular mechanical study

The Thr252 residue plays a vital role in the catalytic cycle of cytochrome P450cam during the formation of the active species (Compound I) from its precursor (Compound 0). We investigate the effect of replacing Thr252 by methoxythreonine (MeO-Thr) on this protonation reaction (coupling) and on the competing formation of the ferric resting state and H₂O₂ (uncoupling) by combined quantum mechanical/molecular mechanical (QM/MM) methods. For each reaction, two possible mechanisms are studied, and for each of these the residues Asp251 and Glu366 are considered as proton sources. The computed QM/MM barriers indicate that uncoupling is unfavorable in the case of the Thr252MeO-Thr mutant, whereas there are two energetically feasible proton transfer pathways for coupling. The corresponding rate-limiting barriers for the formation of Compound I are higher in the mutant than in the wild-type enzyme. These findings are consistent with the experimental observations that the Thr252MeO-Thr mutant forms the alcohol product exclusively (via Compound I), but at lower reaction rates compared with the wild-type enzyme.

Formation of an Amine-Water Cyclic Pentamer: A New Type of Water Cluster in a Polyazacryptand

Structural analysis of an adduct of a thiophene-based cryptand with tosylic acid shows the formation of a hybrid amine-water cyclic pentamer composed of four water molecules and one protonated amine in the charged hydrophobic cavity. The bulky tosylate groups remain outside the cavity, making the ligand favorable for hosting water molecules. Ab initio calculations based on density functional theory (DFT) confirm that the hybrid amine-water pentamer is stabilized inside the hydrophobic cavity of the cryptand.

Resonance Raman characterization of the peroxo and hydroperoxo intermediates in cytochrome P450

Resonance Raman (RR) studies of intermediates generated by cryoreduction of the oxyferrous complex of the D251N mutant of cytochrome P450(cam) (CYP101) are reported. Owing to the fact that proton delivery to the active site is hindered in this mutant, the unprotonated peroxo-ferric intermediate is observed as the primary species after radiolytic reduction of the oxy-complex in frozen solutions at 77 K. In as much as previous EPR and ENDOR studies have shown that annealing of this species to approximately 180 K results in protonation of the distal oxygen atom to form the hydroperoxo intermediate, this system has been exploited to permit direct RR interrogation of the changes in the Fe-O and O-O bonds caused by the reduction and subsequent protonation. Our results show that the nu(O-O) mode decreases from a superoxo-like frequency near approximately 1130 cm(-1) to 792 cm(-1) upon reduction. The latter frequency, as well as its lack of sensitivity to H/D exchange, is consistent with heme-bound peroxide formulation. This species also exhibits a nu(Fe-O) mode, the 553 cm(-1) frequency of which is higher than that observed for the nonreduced oxy P450 precursor (537 cm(-1)), implying a strengthened Fe-O linkage upon reduction. Upon subsequent protonation, the resulting Fe-O-OH fragment exhibits a lowered nu(O-O) mode at 774 cm(-1), whereas the nu(Fe-O) increases to 564 cm(-1). Both modes exhibit a downshift upon H/D exchange, as expected for a hydroperoxo-ferric formulation. These experimental RR data are compared with those previously acquired for the wild-type protein, and the shifts observed upon reduction and subsequent protonation are discussed with reference to theoretical predictions.

Role of protein motions on proton transfer pathways in human carbonic anhydrase II

We report here a theoretical study on the formation of long-range proton transfer pathways in proteins due to side chain conformational fluctuations of amino acid residues and reorganization of interior hydration positions. The proton transfer pathways in such systems may be modeled as fluctuating hydrogen-bonded networks with both short- and long-lived connections between the networked nodes, the latter being formed by polar protein atoms and water molecules. It is known that these fluctuations may extend over several decades of time ranging from a few femtoseconds to a few milliseconds. We have shown in this article how the use of a variety of theoretical methods may be utilized to detect a generic set of pathways and assess the feasibility of forming one or more transient connections. We demonstrate the application of these methods to the enzyme human carbonic anhydrase II and its mutants. Our results reveal several alternative pathways in addition to the one mediated by His-64. We also probe at length the mechanism of key conformational fluctuations contributing to the formation of the detected pathways.

How is the reactivity of cytochrome P450cam affected by Thr252X mutation? A QM/MM study for X = serine, valine, alanine, glycine

Proton transfer reactions play a vital role in the catalytic cycle of cytochrome P450cam and are responsible for the formation of the iron-oxo species called Compound I (Cpd I) that is supposed to be the active oxidant. Depending on the course of the proton transfer, protonation of the last observable intermediate (ferric hydroperoxo complex, Cpd 0) can lead to either the formation of Cpd I (coupling reaction) or the ferric resting state (uncoupling reaction). The ratio of these two processes is drastically affected by mutation of the Thr252 residue. In this work, we study the effect of Thr252X (X = serine, valine, alanine, glycine) mutations on the formation of Cpd I by means of hybrid quantum mechanical/molecular mechanical (QM/MM) calculations and classical simulations. In the wild-type enzyme, the coupling reaction is favored since its rate-limiting barrier is 13 kcal/mol lower than that for uncoupling. This difference is reduced to 7 kcal/mol in the serine mutant. In the case of valine, alanine, and glycine mutants, an additional water molecule enters the active site and lowers the activation energy of the uncoupling reaction significantly. With the additional water molecule, coupling and uncoupling have similar barriers in the valine mutant, and the uncoupling reaction becomes favored in the alanine and glycine mutants. These findings agree very well with experimental results and thus confirm the assumption that uncontrolled proton delivery by solvent water networks is responsible for the uncoupling reaction. The present study provides a detailed mechanistic understanding of the role of the Thr252 residue.

A theoretical study on the detection of proton transfer pathways in some mutants of human carbonic anhydrase II

Structural and kinetic studies of mutants can give much insight into the function of an enzyme. We report the detection of possible proton transfer pathways into the active site of a number of mutants of the enzyme human carbonic anhydrase II (HCA II). Using a recently developed method of path search in the protein conformational space, we identify hydrogen-bonded networks (or proton paths) that can dynamically connect the protein surface to the active site through fluctuations in protein structure and hydration. The feasibility of establishing such dynamical connectivities is assessed by computing the change in free energy of conformational fluctuations and compared to those identified earlier in the wild type enzyme. It is found that the point mutation facilitates or suppresses one or more of the alternative pathways. Our results allow the use of a generic set of pathways to correlate qualitatively the residual activity in the mutants to the molecular mechanism of proton transfer in the absence of His at position 64. We also demonstrate how the detected pathways may be used to compare the efficiencies of the mutants His-64-Ala/Asn-62-His and His-64-Ala/Asn-67-His using the empirical valence bond theory.

Quantum and molecular mechanical study of the first proton transfer in the catalytic cycle of cytochrome P450cam and its mutant D251N

In the catalytic cycle of cytochrome P450cam, the hydroperoxo intermediate (Cpd 0) is formed by proton transfer from a reduced oxyheme complex (S5). This process is drastically slowed down when Asp251 is mutated to Asn (D251N). We report quantum mechanical/molecular mechanical (QM/MM) calculations that address this proton delivery in the doublet state through a hydrogen-bond network in the Asp251 channel, both for the wild-type enzyme and the D251N mutant, using four different active-site models. For the wild-type, we find a facile concerted mechanism for proton transfer from protonated Asp251 via Wat901 and Thr252 to the FeOO moiety, with a barrier of about 1 kcal/mol and a high exothermicity of more than 20 kcal/mol. In the D251N mutant with a neutral Asn251 residue, the proton transfer is almost thermoneutral or slightly exothermic in the three models considered. It is still very facile when the Asn251 residue adopts a conformation analogous to Asp251 in the wild-type enzyme, but the barrier increases significantly when the Asn251 side chain flips (as indicated by classical molecular dynamics simulations). This flip disrupts the hydrogen-bond network and hence the proton-transfer pathway, which causes a longer lifetime of S5 in the D251N mutant (consistent with experimental observations). The entry of an additional water molecule into the active site of D251N with flipped Asn251 regenerates the hydrogen-bond network and provides a viable mechanism for proton delivery in the mutant, with a moderate barrier of about 7 kcal/mol.

Alteration of P450 distal pocket solvent leads to impaired proton delivery and changes in heme geometry

Distal pocket water molecules have been widely implicated in the delivery of protons required in O-O bond heterolysis in the P450 reaction cycle. Targeted dehydration of the cytochrome P450cam (CYP101) distal pocket through mutagenesis of a distal pocket glycine to either valine or threonine results in the alteration of spin state equilibria, and has dramatic consequences on the catalytic rate, coupling efficiency, and kinetic solvent isotope effect parameters, highlighting an important role of the active-site hydration level on P450 catalysis. Cryoradiolysis of the mutant CYP101 oxyferrous complexes further indicates a specific perturbation of proton-transfer events required for the transformation of ferric-peroxo to ferric-hydroperoxo states. Finally, crystallography of the 248Val and 248Thr mutants in both the ferric camphor bound resting state and ferric-cyano adducts shows both the alteration of hydrogen-bonding networks and the alteration of heme geometry parameters. Taken together, these results indicate that the distal pocket microenvironment governs the transformation of reactive heme-oxygen intermediates in P450 cytochromes.

Identification of Proton-Transfer Pathways in Human Carbonic Anhydrase II

We investigate the probable proton-transfer pathways from the surface of human carbonic anhydrase II into the active site cavity through His-64 that has been widely implicated as a key residue along the proton-transfer path. A recursive analysis of hydrogen-bonded clusters in the static crystallographic structure shows that there is no complete path through His-64 in either of its experimentally detected conformations. Side chain conformational fluctuation of His-64 from its outward conformation toward the active site is found to provide a crucial dynamic connectivity needed to complete the path coupled to local reorganization of the protein structure and hydration. The energy and free energy barriers along the detected pathway have been estimated to derive the mechanism of His-64 rotation toward the active site. We also investigate a dynamical connectivity map that highlights networks of disordered water molecules that may promote a direct (and probably transient) access of the solvent to the active site. Our studies reveal how such solvent access channels may be related to the putative proton shuttle mediated by His-64. The paths thus identified can be potentially used as reaction coordinates for further studies on the molecular mechanism of enzyme action.

Energy barriers of proton transfer reactions between amino acid side chain analogs and water from ab initio calculations

Proton transfer reactions were studied in all titratable pairs of amino acid side chains where, under physiologically reasonable conditions, one amino acid may function as a donor and the other one as an acceptor. Energy barriers for shifting the proton from donor to acceptor atom were calculated by electronic structure methods at the MP2/6-31++G(d,p) level, and the well-known double-well potentials were characterized. The energy difference between both minima can be expressed by a parabola using as argument the donor-acceptor distance R(DA). In this work, the fit parameters of the quadratic expression are determined for each donor-acceptor pair. Moreover, it was found previously that the energy barriers of the reactions can be expressed by an analytical expression depending on the distance between donor and acceptor and the energy difference between donor and acceptor bound states. The validity of this approach is supported by the extensive new data set. This new parameterization of proton transfer barriers between titratable amino acid side chains allows us to very efficiently estimate proton transfer probabilities in molecular modelling studies or during classical molecular dynamics simulation of biomolecular systems.

Metastable water clusters in the nonpolar cavities of the thermostable protein tetrabrachion

Water expulsion from the protein core is a key step in protein folding. Nevertheless, unusually large water clusters confined into the nonpolar cavities have been observed in the X-ray crystal structures of tetrabrachion, a bacterial protein that is thermostable up to at least 403 K (130 degrees C). Here, we use molecular dynamics (MD) simulations to investigate the structure and thermodynamics of water filling the largest cavity of the right-handed coiled-coil stalk of tetrabrachion at 365 K (92 degrees C), the temperature of optimal bacterial growth, and at room temperature (298 K). Hydrogen-bonded water clusters of seven to nine water molecules are found to be thermodynamically stable in this cavity at both temperatures, confirming the X-ray studies. Stability, as measured by the transfer free energy of the optimal size cluster, decreases with increasing temperature. Water filling is thus driven by the energy of transfer and opposed by the transfer entropy, both depending only weakly on temperature. Our calculations suggest that cluster formation becomes unfavorable at approximately 384 K (110 degrees C), signaling the onset of drying just slightly above the temperature of optimal growth. "Drying" thus precedes protein denaturation. At room temperature, the second largest cavity in tetrabrachion accommodates a five water molecule cluster, as reported in the X-ray studies. However, the simulations show that at 365 K the cluster is unstable and breaks up. We suggest that the large hydrophobic cavities may act as binding sites for two proteases, possibly explaining the unusual thermostability of the resulting protease-stalk complexes (up to approximately 393 K, 120 degrees C).