Search

Pol η belongs to the important Y family of DNA polymerases that can catalyze translesion synthesis across sites of damaged DNA. This activity involves the reduced fidelity of Pol η for 8-oxo-7,8-dhyedro-2'-deoxoguanosin(8-oxoG). The fundamental interest in Pol η has grown recently with the demonstration of the importance of a 3rd Mg2+ ion. The current work explores both the fidelity of Pol η and the role of the 3rd metal ion, by using empirical valence bond (EVB) simulations. The simulations reproduce the observed trend in fidelity and shed a new light on the role of the 3rd metal ion. It is found that this ion does not lead to a major catalytic effect, but most probably plays an important role in reducing the product release barrier. Furthermore, it is concluded, in contrast to some implications, that the effect of this metal does not violate transition state theory, and the evaluation of the catalytic effect must conserve the molecular composition upon moving from the reactant to the transition state. Proteins 2017; 85:1446-1453. © 2017 Wiley Periodicals, Inc., (© 2017 Wiley Periodicals, Inc.)

Elucidating the catalytic mechanism of DNA polymerase is crucial for a progress in the understanding of the control of replication fidelity. This work tries to advance the mechanistic understanding by analyzing the observed effect of mutations of the acidic groups in the active site of Polymerase β as well as the pH effect on the rate constant. The analysis involves both empirical valence bond (EVB) free energy calculations and considerations of the observed pH dependence of the reaction. The combined analysis indicates that the proton transfer (PT) from the nucleophilic O3' has two possible pathways, one to D256 and the second to the bulk. We concluded based on calculations and the experimental pH profile that the most likely path for the wild-type (WT) and the D256E and D256A mutants is a PT to the bulk, although the WT may also use a PT to Asp 256. Our analysis highlights the need for very extensive sampling in the calculations of the activation barrier and also clearly shows that ab initio QM/MM calculations that do not involve extensive sampling are unlikely to give a clear quantitative picture of the reaction mechanism. Proteins 2016; 84:1644-1657. © 2016 Wiley Periodicals, Inc., (© 2016 Wiley Periodicals, Inc.)

Gaining detailed understanding of the energetics of the proton-pumping process in cytochrome c oxidase (CcO) is one of the challenges of modern biophysics. Despite promising mechanistic proposals, most works have not related the activation barriers of the different assumed steps to the protein structure, and there has not been a physically consistent model that reproduced the barriers needed to create a working pump. This work reevaluates the activation barriers for the primary proton transfer (PT) steps by calculations that reflect all relevant free energy contributions, including the electrostatic energies of the generated charges, the energies of water insertion, and large structural rearrangements of the donor and acceptor. The calculations have reproduced barriers that account for the directionality and sequence of events in the primary PT in CcO. It has also been found that the PT from Glu-286 (E) to the propionate of heme a(3) (Prd) provides a gate for an initial back leakage from the high pH side of the membrane. Interestingly, the rotation of E that brings it closer to Prd appears to provide a way for blocking competing pathways in the primary PT. Our study elucidates and quantifies the nature of the control of the directionality in the primary PT in CcO and provides instructive insight into the role of the water molecules in biological PT, showing that "bridges" of several water molecules in hydrophobic regions present a problem (rather than a solution) that is minimized in the primary PT.

The light-induced proton transport in bacteriorhodopsin has been considered as a model for other light-induced proton pumps. However, the exact nature of this process is still unclear. For example, it is not entirely clear what the driving force of the initial proton transfer is and, in particular, whether it reflects electrostatic forces or other effects. The present work simulates the primary proton transfer (PT) by a specialized combination of the EVB and the QCFF/PI methods. This combination allows us to obtain sufficient sampling and a quantitative free energy profile for the PT at different protein configurations. The calculated profiles provide new insight about energetics of the primary PT and its coupling to the protein conformational changes. Our finding confirms the tentative analysis of an earlier work (A. Warshel, Conversion of light energy to electrostatic energy in the proton pump of Halobacterium halobium, Photochem. Photobiol. 30 (1979) 285-290) and determines that the overall PT process is driven by the energetics of the charge separation between the Schiff base and its counterion Asp85. Apparently, the light-induced relaxation of the steric energy of the chromophore leads to an increase in the ion-pair distance, and this drives the PT process. Our use of the linear response approximation allows us to estimate the change in the protein conformational energy and provides the first computational description of the coupling between the protein structural changes and the PT process. It is also found that the PT is not driven by twist-modulated changes of the Schiff base's pKa, changes in the hydrogen bond directionality, or other non-electrostatic effects. Overall, based on a consistent use of structural information as the starting point for converging free energy calculations, we conclude that the primary event should be described as a light-induced formation of an unstable ground state, whose relaxation leads to charge separation and to the destabilization of the ion-pair state. This provides the driving force for the subsequent PT steps.

Reliable studies of proton transfer (PT) reactions in solution and in enzymes by combined quantum mechanical/molecular mechanics (QM/MM) approaches with an ab initio description of the quantum region present a major challenge to computational chemists. The main problem is the need for extensive computer time to evaluate the QM energy, which in turn makes it extremely challenging to perform proper configurational sampling. The present work presents a new effective way for performing such calculations by using the frozen density functional (FDFT) approach to generate diabatic surfaces that are used to generate a mapping potential that takes the system from the reactant to the product state. The resulting umbrella sampling/free energy perturbation (US/FEP) mapping is done in full analogy with the approach used in the empirical valence bond (EVB) treatment, moving from the diabatic mapping potential to the adiabatic ground state surface, while an ab initio Hamiltonian is used for the QM part. The present approach provides a particularly effective way for evaluating the free energy associated with both the substrate and the solvent motions. This allows us to obtain a free energy barrier that properly reflects the solute entropy. This advance allows one to obtain ab initio QM/MM (QM(ai)/MM) free energy surfaces for very challenging cases such as the autodissociation of water in water, proton transfer between methanol and water in water, and the effect of Mg2+ ion on such a reaction. We also consider as a benchmark the initial PT reaction in the catalytic cycle of triose phosphate isomerase and obtain excellent results without any adjustable parameters. Our results point out that the present implementation of the FDFT approach provides a very promising approach for evaluating QM(ai)/MM free energy surfaces.

The origin of the barrier for proton transport through the aquaporin channel is a problem of general interest. It is becoming increasingly clear that this barrier is not attributable to the orientation of the water molecules across the channel but rather to the electrostatic penalty for moving the proton charge to the center of the channel. However, the reason for the high electrostatic barrier is still rather controversial. It has been argued by some workers that the barrier is due to the so-called NPA motif and/or to the helix macrodipole or to other specific elements. However, our works indicated that the main reason for the high barrier is the loss of the generalized solvation upon moving the proton charge from the bulk to the center of the channel and that this does not reflect a specific repulsive electrostatic interaction but the absence of sufficient electrostatic stabilization. At this stage it seems that the elucidation and clarification of the origin of the electrostatic barrier can serve as an instructive test case for electrostatic models. Thus, we reexamine the free-energy surface for proton transport in aquaporins using the microscopic free-energy perturbation/umbrella sampling (FEP/US) and the empirical valence bond/umbrella sampling (EVB/US) methods as well as the semimacroscopic protein dipole Langevin dipole model in its linear response approximation version (the PDLD/S-LRA). These extensive studies help to clarify the nature of the barrier and to establish the "reduced solvation effect" as the primary source of this barrier. That is, it is found that the barrier is associated with the loss of the generalized solvation energy (which includes of course all electrostatic effects) upon moving the proton charge from the bulk solvent to the center of the channel. It is also demonstrated that the residues in the NPA region and the helix dipole cannot be considered as the main reasons for the electrostatic barrier. Furthermore, our microscopic and semimacroscopic studies clarify the problems with incomplete alternative calculations, illustrating that the effects of various electrostatic elements are drastically overestimated by macroscopic calculations that use a low dielectric constant and do not consider the protein reorganization. Similarly, it is pointed out that microscopic potential of mean force calculations that do not evaluate the electrostatic barrier relative to the bulk water cannot be used to establish the origin of the electrostatic barrier. The relationship between the present study and calculations of pK(a)s in protein interiors is clarified, pointing out that approaches that are applied to study the aquaporin barrier should be validated by pK(a)s calculations. Such calculations also help to clarify the crucial role of solvation energies in establishing the barrier in aquaporins., ((c) 2006 Wiley-Liss, Inc.)

Gaining a detailed understanding of the molecular nature of the redox coupled proton transfer in cytochrome c oxidase (COX) is one of the challenges of modern biophysics. The present work addresses this by integrating approaches for simulations of proton transport (PTR) and electron transfer (ET). The resulting method converts the electrostatic energies of different charge configurations and reorganization energies to free-energy profiles for different PTR and ET pathways. This approach provides for the first time a tool to study the actual activation barriers and kinetics of different feasible PTR processes in the cycle of COX. Using this tool, we explore the PTR through the bottleneck water molecules. It is found that a stepwise PTR along this commonly assumed path leads to far too high barriers and is, thus, inconsistent with the observed kinetics. Furthermore, the simulated free-energy profile does not provide a simple gating mechanism. Fortunately, we obtain reasonable kinetics when we consider a PTR that involves a concerted transfer of protons to and from E286. Finally, semi-qualitative considerations of the forward and backward barriers point toward open questions about the actual gating process and offer a feasible pumping mechanism. Although further studies are clearly needed, we believe that our approach offers a general and effective tool for correlating the structure of COX with its function.

The nature of proton transduction (PTR) through a file of water molecules, along the gramicidin A (gA) channel, has long been considered as being highly relevant to PTR in biological systems. Previous attempts to model this process implied that the so-called Grotthuss mechanism and the corresponding orientation of the water file plays a major role. The present work reexamines the PTR in gA by combining a fully microscopic empirical valence bond (EVB) model and a recently developed simplified EVB-based model with Langevin dynamics (LD) simulations. The full model is used first to evaluate the free energy profile for a stepwise PTR process. The corresponding results are then used to construct the effective potential of the simplified EVB. This later model is then used in Langevin dynamics simulations, taking into account the correct physics of possible concerted motions and the effect of the solvent reorganization. The simulations reproduce the observed experimental trend and lead to a picture that is quite different from that assumed previously. It is found that the PTR in gA is controlled by the change in solvation energy of the transferred proton along the channel axis. Although the time dependent electrostatic fluctuations of the channel and water dipoles play their usual role in modulating the proton-transfer process (Proc. Natl. Acad. Sci. U.S.A. 1984, 81, 444), the PTR rate is mainly determined by the free energy profile. Furthermore, the energetics of the reorientation of the unprotonated water file do not appear to provide a consistent way of assessing the activation barrier for the PTR process. It seems to us that in the case of gA, and probably other systems with significant electrostatic barriers for the transfer of the proton charge, the PTR rate is controlled by the electrostatic barrier. This finding has clear consequences with regards to PTR processes in biological systems.

Proton transport (PTR) processes play a major role in bioenergetics and thus it is important to gain a molecular understanding of these processes. At present the detailed description of PTR in proteins is somewhat unclear and it is important to examine different models by using well-defined experimental systems. One of the best benchmarks is provided by carbonic anhydrase III (CA III), because this is one of the few systems where we have a clear molecular knowledge of the rate constant of the PTR process and its variation upon mutations. Furthermore, this system transfers a proton between several water molecules, thus making it highly relevant to a careful examination of the "proton wire" concept. Obtaining a correlation between the structure of this protein and the rate of the PTR process should help to discriminate between alternative models and to give useful clues about PTR processes in other systems. Obviously, obtaining such a correlation requires a correct representation of the "chemistry" of PTR between different donors and acceptors, as well as the ability to evaluate the free energy barriers of charge transfer in proteins, and to simulate long-time kinetic processes. The microscopic empirical valence bond (Warshel, A., and R. M. Weiss. 1980. J. Am. Chem. Soc. 102:6218-6226; and Aqvist, J., and A. Warshel. 1993. Chem. Rev. 93:2523-2544) provides a powerful way for representing the chemistry and evaluating the free energy barriers, but it cannot be used with the currently available computer times in direct simulation of PTR with significant activation barriers. Alternatively, one can reduce the empirical valence bond (EVB) to the modified Marcus' relationship and use semimacroscopic electrostatic calculations plus a master equation to determine the PTR kinetics (Sham, Y., I. Muegge, and A. Warshel. 1999. Proteins. 36:484-500). However, such an approximation does not provide a rigorous multisite kinetic treatment. Here we combine the useful ingredients of both approaches and develop a simplified EVB effective potential that treats explicitly the chain of donors and acceptors while considering implicitly the rest of the protein/solvent system. This approach can be used in Langevin dynamics simulations of long-time PTR processes. The validity of our new simplified approach is demonstrated first by comparing its Langevin dynamics results for a PTR along a chain of water molecules in water to the corresponding molecular dynamics simulations of the fully microscopic EVB model. This study examines dynamics of both models in cases of low activation barriers and the dependence of the rate on the energetics for cases with moderate barriers. The study of the dependence on the activation barrier is next extended to the range of higher barriers, demonstrating a clear correlation between the barrier height and the rate constant. The simplified EVB model is then examined in studies of the PTR in carbonic anhydrase III, where it reproduces the relevant experimental results without the use of any parameter that is specifically adjusted to fit the energetics or dynamics of the reaction in the protein. We also validate the conclusions obtained previously from the EVB-based modified Marcus' relationship. It is concluded that this approach and the EVB-based model provide a reliable, effective, and general tool for studies of PTR in proteins. Finally in view of the behavior of the simulated result, in both water and the CA III, we conclude that the rate of PTR in proteins is determined by the electrostatic energy of the transferred proton as long as this energy is higher than a few kcal/mol., (Copyright 2004 Biophysical Society)

The nature of the control of water/proton selectivity in biological channels is a problem of a fundamental importance. Most studies of this issue have proposed that an interference with the orientational requirements of the so-called proton wire is the source of selectivity. The elucidation of the structures of aquaporins, which have evolved to prevent proton transfer (PT), provided a clear benchmark for exploring the selectivity problem. Previous simulations of this system have not examined, however, the actual issue of PT, but only considered the much simpler task of the transfer of water molecules. Here we take aquaporin as a benchmark and quantify the origin of the water/proton selectivity in this and related systems. This is done by evaluating in a consistent way the free energy profile for transferring a proton along the channel and relating this profile to the relevant PT rate constants. It is found that the water/proton selectivity is controlled by the change in solvation free energy upon moving the charged proton from water to the channel. The reason for the focus on the elegant concept of the proton wire and the related Grotthuss-type mechanism is also considered. It is concluded that these mechanisms are clearly important in cases with flat free energy surfaces (e.g., in bulk water, in gas phase water chains, and in infinitely long channels). However, in cases of biological channels, the actual PT mechanism is much less important than the energetics of transferring the proton charge from water to different regions in the channels.

Zinc is a vital trace element crucial for the proper function of some 3,000 cellular proteins. Specifically, zinc is essential for key physiological processes including nucleic acid metabolism, regulation of gene expression, signal transduction, cell division, immune- and nervous system functions, wound healing, and apoptosis. Consequently, impairment of zinc homeostasis disrupts key cellular functions resulting in various human pathologies. Mammalian zinc transport proceeds via two transporter families ZnT and ZIP. However, the detailed mechanism of action of ZnT2, which is responsible for vesicular zinc accumulation and zinc secretion into breast milk during lactation, is currently unknown. Moreover, although the putative coupling of zinc transport to the proton gradient in acidic vesicles has been suggested, it has not been conclusively established. Herein we modeled the mechanism of action of ZnT2 and demonstrated both computationally and experimentally, using functional zinc transport assays, that ZnT2 is indeed a proton-coupled zinc antiporter. Bafilomycin A1, a specific inhibitor of vacuolar-type proton ATPase (V-ATPase) which alkalizes acidic vesicles, abolished ZnT2-dependent zinc transport into intracellular vesicles. Moreover, using LysoTracker Red and Lyso-pHluorin, we further showed that upon transient ZnT2 overexpression in intracellular vesicles and addition of exogenous zinc, the vesicular pH underwent alkalization, presumably due to a proton-zinc antiport; this phenomenon was reversed in the presence of TPEN, a specific zinc chelator. Finally, based on computational energy calculations, we propose that ZnT2 functions as an antiporter with a stoichiometry of 2H+/Zn2+ ion. Hence, ZnT2 is a proton motive force-driven, electroneutral vesicular zinc exchanger, concentrating zinc in acidic vesicles on the expense of proton extrusion to the cytoplasm. [ABSTRACT FROM AUTHOR]

The molecular basis of the function of transporters is a problem of significant importance, and the emerging structural information has not yet been converted to a full understanding of the corresponding function. This work explores the molecular origin of the function of the bacterial Na+/H+ antiporter NhaA by evaluating the energetics of the Na+ and H+ movement and then using the resulting landscape in Monte Carlo simulations that examine two transport models and explore which model can reproduce the relevant experimental results. The simulations reproduce the observed transport features by a relatively simple model that relates the protein structure to its transporting function. Focusing on the two key aspartic acid residues of NhaA, D163 and D164, shows that the fully charged state acts as an Na+ trap and that the fully protonated one poses an energetic barrier that blocks the transport of Na+. By alternating between the former and latter states, mediated by the partially protonated protein, protons, and Na+ can be exchanged across the membrane at 2:1 stoichiometry. Our study provides a numerical validation of the need of large conformational changes for effective transport. Furthermore, we also yield a reasonable explanation for the observation that some mammalian transporters have 1:1 stoichiometry. The present coarse-grained model can provide a general way for exploring the function of transporters on a molecular level. [ABSTRACT FROM AUTHOR]