Molecular dynamics simulation directed rational design of inhibitors targeting drug-resistant mutants of influenza A virus M2

Influenza A virus M2 (A/M2) forms a homotetrameric proton selective channel in the viral membrane. It has been the drug target of antiviral drugs such as amantadine and rimantadine. However, most of the current virulent influenza A viruses carry drug-resistant mutations alongside the drug binding site, such as S31N, V27A, and L26F, etc., each of which might be dominant in a given flu season. Among these mutations, the V27A mutation was prevalent among transmissible viruses under drug selection pressure. Until now, V27A has not been successfully targeted by small molecule inhibitors, despite years of extensive medicinal chemistry research efforts and high throughput screening. Guided by molecular dynamics (MD) simulation of drug binding and the influence of drug binding on the dynamics of A/M2 from earlier experimental studies, we designed a series of potent spirane amine inhibitors targeting not only WT, but also both A/M2-27A and L26F mutants with IC(50)s similar to that seen for amantadine's inhibition of the WT channel. The potencies of these inhibitors were further demonstrated in experimental binding and plaque reduction assays. These results demonstrate the power of MD simulations to probe the mechanism of drug binding as well as the ability to guide design of inhibitors of targets that had previously appeared to be undruggable.     

Specific binding of adamantane drugs and direction of their polar amines in the pore of the influenza M2 transmembrane domain in lipid bilayers and dodecylphosphocholine micelles determined by NMR spectroscopy

The transmembrane domain of the influenza M2 protein (M2TM) forms a tetrameric proton channel important for the virus lifecycle. The proton-channel activity is inhibited by amine-containing adamantyl drugs amantadine and rimantadine, which have been shown to bind specifically to the pore of M2TM near Ser31. However, whether the polar amine points to the N- or C-terminus of the channel has not yet been determined. Elucidating the polar group direction will shed light on the mechanism by which drug binding inhibits this proton channel and will facilitate rational design of new inhibitors. In this study, we determine the polar amine direction using M2TM reconstituted in lipid bilayers as well as dodecylphosphocholine (DPC) micelles. (13)C-(2)H rotational-echo double-resonance NMR experiments of (13)C-labeled M2TM and methyl-deuterated rimantadine in lipid bilayers showed that the polar amine pointed to the C-terminus of the channel, with the methyl group close to Gly34. Solution NMR experiments of M2TM in DPC micelles indicate that drug binding causes significant chemical shift perturbations of the protein that are very similar to those seen for M2TM and M2(18-60) bound to lipid bilayers. Specific (2)H-labeling of the drugs permitted the assignment of drug-protein cross peaks, which indicate that amantadine and rimantadine bind to the pore in the same fashion as for bilayer-bound M2TM. These results strongly suggest that adamantyl inhibition of M2TM is achieved not only by direct physical occlusion of the channel, but also by perturbing the equilibrium constant of the proton-sensing residue His37. The reproduction of the pharmacologically relevant specific pore-binding site in DPC micelles, which was not observed with a different detergent, DHPC, underscores the significant influence of the detergent environment on the functional structure of this membrane protein.     

Photoinduced electron transfer and fluorophore motion as a probe of the conformational dynamics of membrane proteins: application to the influenza a M2 proton channel

The structure and function of the influenza A M2 proton channel have been the subject of intensive investigations in recent years because of their critical role in the life cycle of the influenza virus. Using a truncated version of the M2 proton channel (i.e., M2TM) as a model, here we show that fluctuations in the fluorescence intensity of a dye reporter that arise from both fluorescence quenching via the mechanism of photoinduced electron transfer (PET) by an adjacent tryptophan (Trp) residue and local motions of the dye molecule can be used to probe the conformational dynamics of membrane proteins. Specifically, we find that the dynamics of the conformational transition between the N-terminal open and C-terminal open states of the M2TM channel occur on a timescale of about 500 μs and that the binding of either amantadine or rimantadine does not inhibit the pH-induced structural equilibrium of the channel. These results are consistent with the direct occluding mechanism of inhibition which suggests that the antiviral drugs act by sterically occluding the channel pore.     

Gating Mechanism of the Voltage-Gated Proton Channel Studied by Molecular Dynamics Simulations

The voltage-gated proton channel Hv1 has important roles in proton extrusion, pH homeostasis, sperm motility, and cancer progression. The Hv1 channel has also been found to be highly expressed in cell lines and tissue samples from patients with breast cancer. A high-resolution closed-state structure has been reported for the mouse Hv1 chimera channel (mHv1cc), solved by X-ray crystallography, but the open-state structure of Hv1 has not been solved. Since Hv1 is a promising drug target, various groups have proposed open conformations by molecular modeling and simulation studies. However, the gating mechanism and the open-state conformation under the membrane potential are still debate. Here, we present a molecular dynamics study considering membrane potential and pH conditions. The closed-state structure of mHv1cc was used to run molecular dynamics (MD) simulations with respect to electric field and pH conditions in order to investigate the mechanism of proton transfer. We observed a continuous hydrogen bond chain of water molecules called a water-wire to be formed through the channel pore in the channel opening, triggered by downward displacement of the S2 helix and upward movement of the S4 helix relative to other helices. Due to the movement of the S2 and S4 helices, the internal salt bridge network was rearranged, and the hydrophobic gating layers were destroyed. In line with previous experimental and simulation observations, our simulation results led us to propose a new gating mechanism for the Hv1 proton channel, and may provide valuable information for novel drug discovery.     

Membrane-dependent effects of a cytoplasmic helix on the structure and drug binding of the influenza virus M2 protein

The influenza A M2 protein forms a proton channel for virus infection and also mediates virus assembly and budding. The minimum protein length that encodes both functions contains the transmembrane (TM) domain (roughly residues 22-46) for the amantadine-sensitive proton-channel activity and an amphipathic cytoplasmic helix (roughly residues 45-62) for curvature induction and virus budding. However, structural studies involving the TM domain with or without the amphipathic helix differed on the drug-binding site. Here we use solid-state NMR spectroscopy to determine the amantadine binding site in the cytoplasmic-helix-containing M2(21-61). (13)C-(2)H distance measurements of (13)C-labeled protein and (2)H-labeled amantadine showed that in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayers, the first equivalent of drug bound S31 inside the M2(21-61) pore, similar to the behavior of M2 transmembrane peptide (M2TM) in DMPC bilayers. The nonspecific surface site of D44 observed in M2TM is disfavored in the longer peptide. Thus, the pharmacologically relevant drug-binding site in the fully functional M2(21-61) is S31 in the TM pore. Interestingly, when M2(21-61) was reconstituted into a virus-mimetic membrane containing 30% cholesterol, no chemical shift perturbation was observed for pore-lining residues, whereas M2TM in the same membrane exhibited drug-induced chemical shift changes. Reduction of the cholesterol level and the use of unsaturated phospholipids shifted the conformational equilibrium of M2TM fully to the bound state but did not rescue drug binding to M2(21-61). These results suggest that the amphipathic helix, together with cholesterol, modulates the ability of the TM helix to bind amantadine. Thus, the M2 protein interacts with the lipid membrane and small-molecule inhibitors in a complex fashion, and a careful examination of the environmental dependence of the protein conformation is required to fully understand the structure-function relation of this protein.     

Free energy calculations on the two drug binding sites in the M2 proton channel

Two alternative binding sites of adamantane-type drugs in the influenza A M2 channel have been suggested, one with the drug binding inside the channel pore and the other with four drug molecule S-binding to the C-terminal surface of the transmembrane domain. Recent computational and experimental studies have suggested that the pore binding site is more energetically favorable but the external surface binding site may also exist. Nonetheless, which drug binding site leads to channel inhibition in vivo and how drug-resistant mutations affect these sites are not completely understood. We applied molecular dynamics simulations and potential of mean force calculations to examine the structures and the free energies associated with these putative drug binding sites in an M2-lipid bilayer system. We found that, at biological pH (~7.4), the pore binding site is more thermodynamically favorable than the surface binding site by ~7 kcal/mol and, hence, would lead to more stable drug binding and channel inhibition. This result is in excellent agreement with several recent studies. More importantly, a novel finding of ours is that binding to the channel pore requires overcoming a much higher energy barrier of ~10 kcal/mol than binding to the C-terminal channel surface, indicating that the latter site is more kinetically favorable. Our study is the first computational work that provides both kinetic and thermodynamic energy information on these drug binding sites. Our results provide a theoretical framework to interpret and reconcile existing and often conflicting results regarding these two binding sites, thus helping to expand our understanding of M2-drug binding, and may help guide the design and screening of novel drugs to combat the virus.     

Molecular dynamics simulation of a hydrated diphytanol phosphatidylcholine lipid bilayer containing an alpha-helical bundle of four transmembrane domains of the influenza A virus M2 protein

An alpha-helical bundle composed of four transmembrane portions of the M2 protein from the Influenza A virus has been studied in a hydrated diphytanol phosphatidylcholine bilayer using molecular dynamics (MD) calculations. Experimentally, the sequence utilized is known to aggregate as a four-helix bundle and act as a pH-gated proton-selective ion channel, which is blocked by the drug amantadine hydrochloride. In the presented simulation, the ion channel was initially set up as a parallel four-helix bundle. The all-atom simulation consisted of almost 16,000 atoms, described classically, using a forcefield from the CHARMM22 database. Bilayers with and without the bundle were shown to be stable throughout the nanosecond timescale of the MD simulation. Structural and dynamical properties of the bilayer both with and without the transmembrane protein are reported.     

Multiple Proton Confinement in the M2 Channel from the Influenza A Virus

The tetrameric M2 protein bundle of the influenza A virus is the proton channel responsible for the acidification of the viral interior, a key step in the infection cycle. Selective proton transport is achieved by successive protonation of the conserved histidine amino acids at position 37. A recent X-ray structure of the tetrameric transmembrane (TM) domain of the protein (residues 22-46) resolved several water clusters in the channel lumen, which suggest possible proton pathways to the His37 residues. To explore this hypothesis, we have carried out molecular dynamics (MD) simulations of a proton traveling towards the His37 side chains using MD with classical and quantum force fields. Diffusion through the first half of the channel to the "entry" water cluster near His37 may be hampered by significant kinetic barriers due to electrostatic repulsion. However, once in the entry cluster, a proton can move to one of the acceptor His37 in a nearly barrierless fashion, as evidenced both by MD simulations and a scan of the potential energy surface (PES). Water molecules of the entry cluster, although confined in the M2 pore and restricted in their motions, can conduct protons with a rate very similar to that of bulk water.     

Side-chain conformation of the M2 transmembrane peptide proton channel of influenza a virus from 19F solid-state NMR

The M2 transmembrane peptide (M2TMP) of the influenza A virus forms a tetrameric helical bundle that acts as a proton-selective channel important in the viral life cycle. The side-chain conformation of the peptide is largely unknown and is important for elucidating the proton-conducting mechanism and the channel stability. Using a 19F spin diffusion NMR technique called CODEX, we have measured the oligomeric states and interhelical side chain-side chain 19F-19F distances at several residues using singly fluorinated M2TMP bound to DMPC bilayers. 19F CODEX data at a key residue of the proton channel, Trp41, confirm the tetrameric state of the peptide and yield a nearest-neighbor interhelical distance of approximately 11 A under both neutral and acidic pH. Since the helix orientation is precisely known from previous 15N NMR experiments and the backbone channel diameter has a narrow allowed range, this 19F distance constrains the Trp41 side-chain conformation to t90 (chi1 approximately 180 degrees , chi2 approximately 90 degrees ). This Trp41 rotamer, combined with a previously measured 15N-13C distance between His37 and Trp411, suggests that the His37 rotamer is t-160. The implication of the proposed (His37, Trp41) rotamers to the gating mechanism of the M2 proton channel is discussed. Binding of the antiviral drug amantadine to the peptide does not affect the F-F distance at Trp41. Interhelical 19F-19F distances are also measured at residues 27 and 38, each mutated to 4-19F-Phe. For V27F-M2TMP, the 19F-19F distances suggest a mixture of dimers and tetramers, whereas the L38F-M2TMP data indicate two tetramers of different sizes, suggesting side chain conformational heterogeneity at this lipid-facing residue. This work shows that 19F spin diffusion NMR is a valuable tool for determining long-range intermolecular distances that shed light on the mechanism of action and conformational heterogeneity of membrane protein oligomers.     

Kinetic analysis of the M2 proton conduction of the influenza virus

The M2 protein of the flu virus forms a proton selective channel that is necessary for viral replication. The channel has a slow rate of conduction but attains near perfect selectivity for protons. Many models have been proposed to explain the mechanism of proton conduction based on whole cell channel recordings and molecular dynamics simulations, but a detailed kinetic analysis of the channel activity has not yet been performed. We obtained detailed conduction vs pH measurements for M2 and a number of its variants using a sensitive and reproducible liposome proton flux assay. The proton transport follows Michaelis-Menten-like kinetics with two saturation steps: one pseudosaturation at pH ～5.5, and another full saturation at pH ～4. The heart of the mechanism is the pore-lining His37 and Trp41. NMR measurements suggest that histidine and tryptophan act in unison to transport protons down the concentration gradient. The log of apparent K(m) derived from the kinetics data matches closely to the histidine pK(a) and correlates with chemical shift perturbation of the Trp41 gate, indicating that histidine protonation and opening of the channel gate are synchronized events. Finally, mutagenesis and structural analysis identified key residues that affect the rate of conduction.     

Exploring organosilane amines as potent inhibitors and structural probes of influenza a virus M2 proton channel

We describe the use of organosilanes as inhibitors and structural probes of a membrane protein, the M2 proton channel from influenza A virus. Organosilane amine inhibitors were found to be generally as potent as their carbon analogues in targeting WT A/M2 and more potent against the drug-resistant A/M2-V27A mutant. In addition, intermolecular NOESY spectra with dimethyl-substituted organosilane amine inhibitors clearly located the drug binding site at the N-terminal lumen of the A/M2 channel close to V27.     

In Situ Spectroscopic Detection of Large-Scale Reorientations of Transmembrane Helices During Influenza A M2 Channel Opening**

Viroporins are small ion channels in membranes of enveloped viruses that play key roles during viral life cycles. To use viroporins as drug targets against viral infection requires in-depth mechanistic understanding and, with that, methods that enable investigations under in situ conditions. Here, we apply surface-enhanced infrared absorption (SEIRA) spectroscopy to Influenza A M2 reconstituted within a solid-supported membrane, to shed light on the mechanics of its viroporin function. M2 is a paradigm of pH-activated proton channels and controls the proton flux into the viral interior during viral infection. We use SEIRA to track the large-scale reorientation of M2’s transmembrane α-helices in situ during pH-activated channel opening. We quantify this event as a helical tilt from 26° to 40° by correlating the experimental results with solid-state nuclear magnetic resonance-informed computational spectroscopy. This mechanical motion is impeded upon addition of the inhibitor rimantadine, giving a direct spectroscopic marker to test antiviral activity. The presented approach provides a spectroscopic tool to quantify large-scale structural changes and to track the function and inhibition of the growing number of viroporins from pathogenic viruses in future studies.     

Simulating the proton transfer in gramicidin A by a sequential dynamical Monte Carlo method

The large interest in long-range proton transfer in biomolecules is triggered by its importance for many biochemical processes such as biological energy transduction and drug detoxification. Since long-range proton transfer occurs on a microsecond time scale, simulating this process on a molecular level is still a challenging task and not possible with standard simulation methods. In general, the dynamics of a reactive system can be described by a master equation. A natural way to describe long-range charge transfer in biomolecules is to decompose the process into elementary steps which are transitions between microstates. Each microstate has a defined protonation pattern. Although such a master equation can in principle be solved analytically, it is often too demanding to solve this equation because of the large number of microstates. In this paper, we describe a new method which solves the master equation by a sequential dynamical Monte Carlo algorithm. Starting from one microstate, the evolution of the system is simulated as a stochastic process. The energetic parameters required for these simulations are determined by continuum electrostatic calculations. We apply this method to simulate the proton transfer through gramicidin A, a transmembrane proton channel, in dependence on the applied membrane potential and the pH value of the solution. As elementary steps in our reaction, we consider proton uptake and release, proton transfer along a hydrogen bond, and rotations of water molecules that constitute a proton wire through the channel. A simulation of 8 mus length took about 5 min on an Intel Pentium 4 CPU with 3.2 GHz. We obtained good agreement with experimental data for the proton flux through gramicidin A over a wide range of pH values and membrane potentials. We find that proton desolvation as well as water rotations are equally important for the proton transfer through gramicidin A at physiological membrane potentials. Our method allows to simulate long-range charge transfer in biological systems at time scales, which are not accessible by other methods.     

NMR detection of pH-dependent histidine-water proton exchange reveals the conduction mechanism of a transmembrane proton channel

The acid-activated proton channel formed by the influenza M2 protein is important for the life cycle of the virus. A single histidine, His37, in the M2 transmembrane domain (M2TM) is responsible for pH activation and proton selectivity of the channel. Recent studies suggested three models for how His37 mediates proton transport: a shuttle mechanism involving His37 protonation and deprotonation, a H-bonded imidazole-imidazolium dimer model, and a transporter model involving large protein conformational changes in synchrony with proton conduction. Using magic-angle-spinning (MAS) solid-state NMR spectroscopy, we examined the proton exchange and backbone conformational dynamics of M2TM in a virus-envelope-mimetic membrane. At physiological temperature and pH, (15)N NMR spectra show fast exchange of the imidazole (15)N between protonated and unprotonated states. To quantify the proton exchange rates, we measured the (15)N T(2) relaxation times and simulated them for chemical-shift exchange and fluctuating N-H dipolar fields under (1)H decoupling and MAS. The exchange rate is 4.5 × 10(5) s(-1) for Nδ1 and 1.0 × 10(5) s(-1) for Nε2, which are approximately synchronized with the recently reported imidazole reorientation. Binding of the antiviral drug amantadine suppressed both proton exchange and ring motion, thus interfering with the proton transfer mechanism. By measuring the relative concentrations of neutral and cationic His as a function of pH, we determined the four pK(a) values of the His37 tetrad in the viral membrane. Fitting the proton current curve using the charge-state populations from these pK(a)'s, we obtained the relative conductance of the five charge states, which showed that the +3 channel has the highest time-averaged unitary conductance. At physiologically relevant pH, 2D correlation spectra indicated that the neutral and cationic histidines do not have close contacts, ruling out the H-bonded dimer model. Moreover, a narrowly distributed nonideal helical structure coexists with a broadly distributed ideal helical conformation without interchange on the sub-10 ms time scale, thus excluding the transporter model in the viral membrane. These data support the shuttle mechanism of proton conduction, whose essential steps involve His-water proton exchange facilitated by imidazole ring reorientations.     

Recent progress in computational approaches to studying the M2 proton channel and its implication to drug design against influenza viruses

For quite a long period of time in history, many intense efforts have been made to determine the 3D (three-dimensional) structure of the M2 proton channel. The reason why the M2 proton channel has attracted so many attentions is because (1) it is the key for really understanding the life cycle of influenza viruses, and (2) it is indispensable for conducting rational drug design against the flu viruses. Recently, the long-sough 3D structures of the M2 proton channels for both influenza A and B viruses were consecutively successfully determined by the high-resolution NMR spectroscopy (Schnell J.R. and Chou, J.J., Nature, 2008, 451: 591-595; Wang, J., Pielak, R.M., McClintock, M.A., and Chou, J.J., Nature Structural & Molecular Biology, 2009,16: 1267-1271). Such a milestone work has provided a solid structural basis for in-depth understanding the action mechanism of the M2 channel and rationally designing effective drugs against influenza viruses. This review is devoted to, with the focus on the M2 proton channel of influenza A, addressing a series of relevant problems, such as how to correctly understand the novel allosteric inhibition mechanism inferred from the NMR structure that is completely different from the traditional view, what the possible impacts are to the previous functional studies in this area, and what kind of new strategy can be stimulated for drug development against influenza.     

Simulation of proton transport in the gramicidin A channel

The free energy profiles for proton transfer along the oriented water file inside the gramicidin A channel were calculated. An original implementation of the rigid-body molecular dynamics method was used for describing the peptide groups of the channel and outer water molecules. The inner water wire was simulated using the PM6 force field parameters, which adequately describe the formation and cleavage of chemical and hydrogen bonds in water molecules. Different mechanisms of proton transfer through the gramicidin A channel were considered, namely, proton H⁺ translocation, transfer of the anion defect OH^?, and reorientation of the water file inside the channel. To facilitate parallel calculations of trajectories, the reaction coordinate was divided into segments, and the results were combined by the weighted histogram analysis method. The first two processes, H⁺ and OH^? transfers, were shown to be barrierless. Only the stage of reorientation of the water file inside the channel has an energy barrier.     

Parallel proton transfer pathways in aqueous acid-base reactions

We study the mechanism of proton transfer (PT) between the photoacid 8-hydroxy-1,3, 6-pyrenetrisulfonic acid (HPTS) and the base chloroacetate in aqueous solution. We investigate both proton and deuteron transfer reactions in solutions with base concentrations ranging from 0.25 M to 4 M. Using femtosecond midinfrared spectroscopy, we probe the vibrational responses of HPTS, its conjugate photobase, the hydrated proton/deuteron, and chloroacetate. The measurement of these four resonances allows us to follow the sequence of proton departure from the acid, its uptake by the water solvent, and its arrival at the base. In recent studies it was shown that proton transfer to carboxylate bases proceeds via Grotthuss conduction through a water wire connecting the acid and the base [Mohammed et al., Science 310, 83 (2005);Agnew. Chem. Int. Ed. 46, 1458 (2007);Siwick and Bakker, J. Am. Chem. Soc. 129, 13412 (2007); J. Phys. Chem. B 112, 378 (2008)]. Here we show that, for the weaker base chloroacetate, an alternative channel for proton transfer arises. In this channel the proton is first transferred to the water solvent and only later taken up from the water by the base. We study the base concentration dependence of the two competing channels.     

Enantiomer-specific binding of ruthenium(II) molecular wires by the amine oxidase of Arthrobacter globiformis

The copper amine oxidase from Arthrobacter globiformis (AGAO) is reversibly inhibited by molecular wires comprising a Ru(II) complex head group and an aromatic tail group joined by an alkane linker. The crystal structures of a series of Ru(II)-wire-AGAO complexes differing with respect to the length of the alkane linker have been determined. All wires lie in the AGAO active-site channel, with their aromatic tail group in contact with the trihydroxyphenylalanine quinone (TPQ) cofactor of the enzyme. The TPQ cofactor is consistently in its active ("off-Cu") conformation, and the side chain of the so-called "gate" residue Tyr296 is consistently in the "gate-open" conformation. Among the wires tested, the most stable complex is produced when the wire has a -(CH2)4- linker. In this complex, the Ru(II)(phen)(bpy)2 head group is level with the protein molecular surface. Crystal structures of AGAO in complex with optically pure forms of the C4 wire show that the linker and head group in the two enantiomers occupy slightly different positions in the active-site channel. Both the Lambda and Delta isomers are effective competitive inhibitors of amine oxidation. Remarkably, inhibition by the C4 wire shows a high degree of selectivity for AGAO in comparison with other copper-containing amine oxidases.     

Unusual "amphiphilic" association of hydrated protons in strong acid solution

The solvation structure of the hydrated excess proton in concentrated aqueous HCl solution is studied using the self-consistent iterative multi-state empirical valence bond method. At 0.43-0.85 M concentrations, hydronium cations are found to form unusual cation pairs. This behavior is consistent with our earlier finding that hydronium cations can have an "amphiphilic" character due in part to the asymmetric nature of their hydrogen bonding to nearby water molecules. The existence of these hydronium amphiphilic pairs is further supported by a Car-Parrinello ab initio molecular dynamics simulation at 1.0 M HCl concentration. It is also found that the hydronium cation pairs are stabilized by a delocalization of the hydrated excess proton charge defects involving additional water molecules. At the higher concentrations of 1.68 and 3.26 M, the abundance of such hydronium pairs decreases, and the analysis of the radial distribution functions indicates the possible formation of an aggregate structure with longer-ranged order.     

Intricate role of water in proton transport through cytochrome c oxidase

Cytochrome c oxidase (CytcO), the final electron acceptor in the respiratory chain, catalyzes the reduction of O(2) to H(2)O while simultaneously pumping protons across the inner mitochondrial or bacterial membrane to maintain a transmembrane electrochemical gradient that drives, for example, ATP synthesis. In this work mutations that were predicted to alter proton translocation and enzyme activity in preliminary computational studies are characterized with extensive experimental and computational analysis. The mutations were introduced in the D pathway, one of two proton-uptake pathways, in CytcO from Rhodobacter sphaeroides . Serine residues 200 and 201, which are hydrogen-bonded to crystallographically resolved water molecules halfway up the D pathway, were replaced by more bulky hydrophobic residues (Ser200Ile, Ser200Val/Ser201Val, and Ser200Val/Ser201Tyr) to query the effects of changing the local structure on enzyme activity as well as proton uptake, release, and intermediate transitions. In addition, the effects of these mutations on internal proton transfer were investigated by blocking proton uptake at the pathway entrance (Asp132Asn replacement in addition to the above-mentioned mutations). Even though the overall activities of all mutant CytcO's were lowered, both the Ser200Ile and Ser200Val/Ser201Val variants maintained the ability to pump protons. The lowered activities were shown to be due to slowed oxidation kinetics during the P(R) → F and F → O transitions (P(R) is the "peroxy" intermediate formed at the catalytic site upon reaction of the four-electron-reduced CytcO with O(2), F is the oxoferryl intermediate, and O is the fully oxidized CytcO). Furthermore, the P(R) → F transition is shown to be essentially pH independent up to pH 12 (i.e., the apparent pK(a) of Glu286 is increased from 9.4 by at least 3 pK(a) units) in the Ser200Val/Ser201Val mutant. Explicit simulations of proton transport in the mutated enzymes revealed that the solvation dynamics can cause intriguing energetic consequences and hence provide mechanistic insights that would never be detected in static structures or simulations of the system with fixed protonation states (i.e., lacking explicit proton transport). The results are discussed in terms of the proton-pumping mechanism of CytcO.