A delocalized proton-binding site within a membrane protein

The Functionality of the DC Pair in a Rhodopsin Guanylyl Cyclase from Catenaria anguillulae

Rhodopsin guanylyl cyclases (RGCs) belong to the class of enzymerhodopsins catalyzing the transition from GTP into the second messenger cGMP, whereas light-regulation of enzyme activity is mediated by a membrane-bound microbial rhodopsin domain, that holds the catalytic center inactive in the dark. Structural determinants for activation of the rhodopsin moiety eventually leading to catalytic activity are largely unknown. Here, we investigate the mechanistic role of the D283-C259 (DC) pair that is hydrogen bonded via a water molecule as a crucial functional motif in the homodimeric C. anguillulae RGC. Based on a structural model of the DC pair in the retinal binding pocket obtained by MD simulation, we analyzed formation and kinetics of early and late photocycle intermediates of the rhodopsin domain wild type and specific DC pair mutants by combined UV-Vis and FTIR spectroscopy at ambient and cryo-temperatures. By assigning specific infrared bands to S-H vibrations of C259 we are able to show that the DC pair residues are tightly coupled. We show that deprotonation of D283 occurs already in the inactive L state as a prerequisite for M state formation, whereas structural changes of C259 occur in the active M state and early cryo-trapped intermediates. We propose a comprehensive molecular model for formation of the M state that activates the catalytic moiety. It involves light induced changes in bond strength and hydrogen bonding of the DC pair residues from the early J state to the active M state and explains the retarding effect of C259 mutants.

Electrometric and Electron Paramagnetic Resonance Measurements of a Difference in the Transmembrane Electrochemical Potential: Photosynthetic Subcellular Structures and Isolated Pigment-Protein Complexes

A transmembrane difference in the electrochemical potentials of protons (ΔμH+) serves as a free energy intermediate in energy-transducing organelles of the living cell. The contributions of two components of the ΔμH+ (electrical, Δψ, and concentrational, ΔpH) to the overall ΔμH+ value depend on the nature and lipid composition of the energy-coupling membrane. In this review, we briefly consider several of the most common instrumental (electrometric and EPR) methods for numerical estimations of Δψ and ΔpH. In particular, the kinetics of the flash-induced electrometrical measurements of Δψ in bacterial chromatophores, isolated bacterial reaction centers, and Photosystems I and II of the oxygenic photosynthesis, as well as the use of pH-sensitive molecular indicators and kinetic data regarding pH-dependent electron transport in chloroplasts, have been reviewed. Further perspectives on the application of these methods to solve some fundamental and practical problems of membrane bioenergetics are discussed.

Protons in Gating the Kv1.2 Channel: A Calculated Set of Protonation States in Response to Polarization/Depolarization of the Channel, with the Complete Proposed Proton Path from Voltage Sensing Domain to Gate

We have in the past proposed that proton motion constitutes the gating current in the potassium channel K_v1.2 and is responsible for the gating mechanism. For this to happen, there must be a proton path between the voltage-sensing domain (VSD) and the channel gate, and here we present quantum calculations that lead to a specific pair of proton paths, defined at the molecular level, with well-defined water molecule linkages, and with hydrogen bonding between residues; there is also at least one interpath crossover, where protons can switch paths. Quantum calculations on the entire 563-atom system give the complete geometry, the energy, and atomic charges. Calculations show that three specific residues (in the pdb 3Lut numbering, H418, E327, R326), and the T1 intracellular moiety, all of which have been shown experimentally to be involved in gating, would necessarily be protonated or deprotonated in the path between the VSD and the gate. Hydroxyl reorientation of serine and threonine residues are shown to provide a means of adjusting proton directions of motion. In the deprotonated state for K312, a low energy state, our calculations come close to reproducing the X-ray structure. The demonstration of the existence of a double proton path between VSD and gate supports the proposed proton gating mechanism; when combined with our earlier demonstration of proton generation in the VSD, and comparison with other systems that are known to move protons, we are close to achieving the definition of a complete gating mechanism in molecular detail. The coupling of the paths to the VSD, and to the PVPV section that essentially forms the gate, can be easily seen from the results of the calculation. The gate itself remains for further computations.

Archaeal Lipids Regulating the Trimeric Structure Dynamics of Bacteriorhodopsin for Efficient Proton Release and Uptake

S-TGA-1 and PGP-Me are native archaeal lipids associated with the bacteriorhodopsin (bR) trimer and contribute to protein stabilization and native dynamics for proton transfer. However, little is known about the underlying molecular mechanism of how these lipids regulate bR trimerization and efficient photocycling. Here, we explored the specific binding of S-TGA-1 and PGP-Me with the bR trimer and elucidated how specific interactions modulate the bR trimeric structure and proton release and uptake using long-term atomistic molecular dynamic simulations. Our results showed that S-TGA-1 and PGP-Me are essential for stabilizing the bR trimer and maintaining the coherent conformational dynamics necessary for proton transfer. The specific binding of S-TGA-1 with W80 and K129 regulates proton release on the extracellular surface by forming a “Glu-shared” model. The interaction of PGP-Me with K40 ensures proton uptake by accommodating the conformation of the helices to recruit enough water molecules on the cytoplasmic side. The present study results could fill in the theoretical gaps of studies on the functional role of archaeal lipids and could provide a reference for other membrane proteins containing similar archaeal lipids.

True-atomic-resolution insights into the structure and functional role of linear chains and low-barrier hydrogen bonds in proteins

Hydrogen bonds are fundamental to the structure and function of biological macromolecules and have been explored in detail. The chains of hydrogen bonds (CHBs) and low-barrier hydrogen bonds (LBHBs) were proposed to play essential roles in enzyme catalysis and proton transport. However, high-resolution structural data from CHBs and LBHBs is limited. The challenge is that their 'visualization' requires ultrahigh-resolution structures of the ground and functionally important intermediate states to identify proton translocation events and perform their structural assignment. Our true-atomic-resolution structures of the light-driven proton pump bacteriorhodopsin, a model in studies of proton transport, show that CHBs and LBHBs not only serve as proton pathways, but also are indispensable for long-range communications, signaling and proton storage in proteins. The complete picture of CHBs and LBHBs discloses their multifunctional roles in providing protein functions and presents a consistent picture of proton transport and storage resolving long-standing debates and controversies.

Predicting Ion Channel Conductance via Dissipation-Corrected Targeted Molecular Dynamics and Langevin Equation Simulations

Ion channels are important proteins for physiological information transfer and functional control. To predict the microscopic origins of their voltage-conductance characteristics, here we applied dissipation-corrected targeted molecular dynamics in combination with Langevin equation simulations to potassium diffusion through the gramicidin A channel as a test system. Performing a nonequilibrium principal component analysis on backbone dihedral angles, we find coupled protein-ion dynamics to occur during ion transfer. The dissipation-corrected free energy profiles correspond well to predictions from other biased simulation methods. The incorporation of an external electric field in Langevin simulations enables the prediction of macroscopic observables in the form of I-V characteristics.

Quantum-Mechanical Molecular Dynamics Simulations on Secondary Proton Transfer in Bacteriorhodopsin Using Realistic Models

Bacteriorhodopsin (BR) transports a proton from intracellular to extracellular (EC) sites through five proton transfers. The second proton transfer is the release of an excess proton stored in BR into the EC medium, and an atomistic understanding of this whole process has remained unexplored due to its ubiquitous environment. Here, fully quantum mechanical (QM) molecular dynamics (MD) and metadynamics (MTD) simulations for this process were performed at the divide-and-conquer density-functional tight-binding level using realistic models (∼50000 and ∼20000 atoms) based on the time-resolved photointermediate structures from an X-ray free electron laser. Regarding the proton storage process, the QM-MD/MTD simulations confirmed the Glu-shared mechanism, in which an excess proton is stored between Glu194 and Glu204, and clarified that the activation occurs by localizing the proton at Glu204 in the photocycle. Furthermore, the QM-MD/MTD simulations elucidated a release pathway from Glu204 through Ser193 to the EC water molecules and clarified that the proton release starts at ∼250 μs. In the ubiquitous proton diffusion in the EC medium, the transient proton receptors predicted experimentally were assigned to carboxylates in Glu9 and Glu74. Large-scale QM-MD/MTD simulations beyond the conventional sizes, which provided the above findings and confirmations, were possible by adopting our Dcdftbmd program.

Anharmonic Vibrational Calculations Based on Group-Localized Coordinates: Applications to Internal Water Molecules in Bacteriorhodopsin

An efficient anharmonic vibrational method is developed exploiting the locality of molecular vibration. Vibrational coordinates localized to a group of atoms are employed to divide the potential energy surface (PES) of a system into intra- and inter-group contributions. Then, the vibrational Schrödinger equation is solved based on a PES, in which the inter-group coupling is truncated at the harmonic level while accounting for the intra-group anharmonicity. The method is applied to a pentagonal hydrogen bond network (HBN) composed of internal water molecules and charged residues in a membrane protein, bacteriorhodopsin. The PES is calculated by the quantum mechanics/molecular mechanics (QM/MM) calculation at the level of B3LYP-D3/aug-cc-pVDZ. The infrared (IR) spectrum is computed using a set of coordinates localized to each water molecule and amino acid residue by second-order vibrational quasi-degenerate perturbation theory (VQDPT2). Benchmark calculations show that the proposed method yields the N-D/O-D stretching frequencies with an error of 7 cm-1 at the cost reduced by more than five times. In contrast, the harmonic approximation results in a severe error of 150 cm-1. Furthermore, the size of QM regions is carefully assessed to find that the QM regions should include not only the pentagonal HBN itself but also its HB partners. VQDPT2 calculations starting from transient structures obtained by molecular dynamics simulations have shown that the structural sampling has a significant impact on the calculated IR spectrum. The incorporation of anharmonicity, sufficiently large QM regions, and structural samplings are of essential importance to reproduce the experimental IR spectrum. The computational spectrum paves the way for decoding the IR signal of strong HBNs and helps elucidate their functional roles in biomolecules.

Protein Motifs for Proton Transfers That Build the Transmembrane Proton Gradient

Biological membranes are barriers to polar molecules, so membrane embedded proteins control the transfers between cellular compartments. Protein controlled transport moves substrates and activates cellular signaling cascades. In addition, the electrochemical gradient across mitochondrial, bacterial and chloroplast membranes, is a key source of stored cellular energy. This is generated by electron, proton and ion transfers through proteins. The gradient is used to fuel ATP synthesis and to drive active transport. Here the mechanisms by which protons move into the buried active sites of Photosystem II (PSII), bacterial RCs (bRCs) and through the proton pumps, Bacteriorhodopsin (bR), Complex I and Cytochrome c oxidase (CcO), are reviewed. These proteins all use water filled proton transfer paths. The proton pumps, that move protons uphill from low to high concentration compartments, also utilize Proton Loading Sites (PLS), that transiently load and unload protons and gates, which block backflow of protons. PLS and gates should be synchronized so PLS proton affinity is high when the gate opens to the side with few protons and low when the path is open to the high concentration side. Proton transfer paths in the proteins we describe have different design features. Linear paths are seen with a unique entry and exit and a relatively straight path between them. Alternatively, paths can be complex with a tangle of possible routes. Likewise, PLS can be a single residue that changes protonation state or a cluster of residues with multiple charge and tautomer states.

A protet-based, protonic charge transfer model of energy coupling in oxidative and photosynthetic phosphorylation

Textbooks of biochemistry will explain that the otherwise endergonic reactions of ATP synthesis can be driven by the exergonic reactions of respiratory electron transport, and that these two half-reactions are catalyzed by protein complexes embedded in the same, closed membrane. These views are correct. The textbooks also state that, according to the chemiosmotic coupling hypothesis, a (or the) kinetically and thermodynamically competent intermediate linking the two half-reactions is the electrochemical difference of protons that is in equilibrium with that between the two bulk phases that the coupling membrane serves to separate. This gradient consists of a membrane potential term Δψ and a pH gradient term ΔpH, and is known colloquially as the protonmotive force or pmf. Artificial imposition of a pmf can drive phosphorylation, but only if the pmf exceeds some 150-170mV; to achieve in vivo rates the imposed pmf must reach 200mV. The key question then is 'does the pmf generated by electron transport exceed 200mV, or even 170mV?' The possibly surprising answer, from a great many kinds of experiment and sources of evidence, including direct measurements with microelectrodes, indicates it that it does not. Observable pH changes driven by electron transport are real, and they control various processes; however, compensating ion movements restrict the Δψ component to low values. A protet-based model, that I outline here, can account for all the necessary observations, including all of those inconsistent with chemiosmotic coupling, and provides for a variety of testable hypotheses by which it might be refined.

Density-functional tight-binding: basic concepts and applications to molecules and clusters

The scope of this article is to present an overview of the Density Functional based Tight Binding (DFTB) method and its applications. The paper introduces the basics of DFTB and its standard formulation up to second order. It also addresses methodological developments such as third order expansion, inclusion of non-covalent interactions, schemes to solve the self-interaction error, implementation of long-range short-range separation, treatment of excited states via the time-dependent DFTB scheme, inclusion of DFTB in hybrid high-level/low level schemes (DFT/DFTB or DFTB/MM), fragment decomposition of large systems, large scale potential energy landscape exploration with molecular dynamics in ground or excited states, non-adiabatic dynamics. A number of applications are reviewed, focusing on -(i)- the variety of systems that have been studied such as small molecules, large molecules and biomolecules, bare orfunctionalized clusters, supported or embedded systems, and -(ii)- properties and processes, such as vibrational spectroscopy, collisions, fragmentation, thermodynamics or non-adiabatic dynamics. Finally outlines and perspectives are given.

Estimation of Protein-Ligand Unbinding Kinetics Using Non-Equilibrium Targeted Molecular Dynamics Simulations

We here report on nonequilibrium targeted molecular dynamics simulations as a tool for the estimation of protein-ligand unbinding kinetics. Correlating simulations with experimental data from SPR kinetics measurements and X-ray crystallography on two small molecule compound libraries bound to the N-terminal domain of the chaperone Hsp90, we show that the mean nonequilibrium work computed in an ensemble of trajectories of enforced ligand unbinding is a promising predictor for ligand unbinding rates. We furthermore investigate the molecular basis determining unbinding rates within the compound libraries. We propose ligand conformational changes and protein-ligand nonbonded interactions to impact on unbinding rates. Ligands may remain longer at the protein if they exhibit strong electrostatic and/or van der Waals interactions with the target. In the case of ligands with a rigid chemical scaffold that exhibit longer residence times, transient electrostatic interactions with the protein appear to facilitate unbinding. Our results imply that understanding the unbinding pathway and the protein-ligand interactions along this path is crucial for the prediction of small molecule ligands with defined unbinding kinetics.

X-ray structure analysis of bacteriorhodopsin at 1.3 Å resolution

Bacteriorhodopsin (bR) of Halobacterium salinarum is a membrane protein that acts as a light-driven proton pump. bR and its homologues have recently been utilized in optogenetics and other applications. Although the structures of those have been reported so far, the resolutions are not sufficient for elucidation of the intrinsic structural features critical to the color tuning and ion pumping properties. Here we report the accurate crystallographic analysis of bR in the ground state. The influence of X-rays was suppressed by collecting the data under a low irradiation dose at 15 K. Consequently, individual atoms could be separately observed in the electron density map at better than 1.3 Å resolution. Residues from Thr5 to Ala233 were continuously constructed in the model. The twist of the retinal polyene was determined to be different from those in the previous models. Two conformations were observed for the proton release region. We discuss the meaning of these fine structural features.

pH-sensitive vibrational probe reveals a cytoplasmic protonated cluster in bacteriorhodopsin

Infrared spectroscopy has been used in the past to probe the dynamics of internal proton transfer reactions taking place during the functional mechanism of proteins but has remained mostly silent to protonation changes in the aqueous medium. Here, by selectively monitoring vibrational changes of buffer molecules with a temporal resolution of 6 µs, we have traced proton release and uptake events in the light-driven proton-pump bacteriorhodopsin and correlate these to other molecular processes within the protein. We demonstrate that two distinct chemical entities contribute to the temporal evolution and spectral shape of the continuum band, an unusually broad band extending from 2,300 to well below 1,700 cm^-1 The first contribution corresponds to deprotonation of the proton release complex (PRC), a complex in the extracellular domain of bacteriorhodopsin where an excess proton is shared by a cluster of internal water molecules and/or ionic E194/E204 carboxylic groups. We assign the second component of the continuum band to the proton uptake complex, a cluster with an excess proton reminiscent to the PRC but located in the cytoplasmic domain and possibly stabilized by D38. Our findings refine the current interpretation of the continuum band and call for a reevaluation of the last proton transfer steps in bacteriorhodopsin.

Perspective: Quantum mechanical methods in biochemistry and biophysics

In this perspective article, I discuss several research topics relevant to quantum mechanical (QM) methods in biophysical and biochemical applications. Due to the immense complexity of biological problems, the key is to develop methods that are able to strike the proper balance of computational efficiency and accuracy for the problem of interest. Therefore, in addition to the development of novel ab initio and density functional theory based QM methods for the study of reactive events that involve complex motifs such as transition metal clusters in metalloenzymes, it is equally important to develop inexpensive QM methods and advanced classical or quantal force fields to describe different physicochemical properties of biomolecules and their behaviors in complex environments. Maintaining a solid connection of these more approximate methods with rigorous QM methods is essential to their transferability and robustness. Comparison to diverse experimental observables helps validate computational models and mechanistic hypotheses as well as driving further development of computational methodologies.

Local Mode Analysis: Decoding IR Spectra by Visualizing Molecular Details

Integration of experimental and computational approaches to investigate chemical reactions in proteins has proven to be very successful. Experimentally, time-resolved FTIR difference-spectroscopy monitors chemical reactions at atomic detail. To decode detailed structural information encoded in IR spectra, QM/MM calculations are performed. Here, we present a novel method which we call local mode analysis (LMA) for calculating IR spectra and assigning spectral IR-bands on the basis of movements of nuclei and partial charges from just a single QM/MM trajectory. Through LMA the decoding of IR spectra no longer requires several simulations or optimizations. The novel approach correlates the motions of atoms of a single simulation with the corresponding IR bands and provides direct access to the structural information encoded in IR spectra. Either the contributions of a particular atom or atom group to the complete IR spectrum of the molecule are visualized, or an IR-band is selected to visualize the corresponding structural motions. Thus, LMA decodes the detailed information contained in IR spectra and provides an intuitive approach for structural biologists and biochemists. The unique feature of LMA is the bidirectional analysis connecting structural details to spectral features and vice versa spectral features to molecular motions.

QM/MM Analysis of Transition States and Transition State Analogues in Metalloenzymes

Enzymology is approaching an era where many problems can benefit from computational studies. While ample challenges remain in quantitatively predicting behavior for many enzyme systems, the insights that often come from computations are an important asset for the enzymology community. Here we provide a primer for enzymologists on the types of calculations that are most useful for mechanistic problems in enzymology. In particular, we emphasize the integration of models that range from small active-site motifs to fully solvated enzyme systems for cross-validation and dissection of specific contributions from the enzyme environment. We then use a case study of the enzyme alkaline phosphatase to illustrate specific application of the methods. The case study involves examination of the binding modes of putative transition state analogues (tungstate and vanadate) to the enzyme. The computations predict covalent binding of these ions to the enzymatic nucleophile and that they adopt the trigonal bipyramidal geometry of the expected transition state. By comparing these structures with transition states found through free energy simulations, we assess the degree to which the transition state analogues mimic the true transition states. Technical issues worth treating with care as well as several remaining challenges to quantitative analysis of metalloenzymes are also highlighted during the discussion.

Proteorhodopsin Activation Is Modulated by Dynamic Changes in Internal Hydration

Proteorhodopsin, a member of the microbial rhodopsin family, is a seven-transmembrane α-helical protein that functions as a light-driven proton pump. Understanding the proton-pumping mechanism of proteorhodopsin requires intimate knowledge of the proton transfer pathway via complex hydrogen-bonding networks formed by amino acid residues and internal water molecules. Here we conducted a series of microsecond time scale molecular dynamics simulations on both the dark state and the initial photoactivated state of blue proteorhodopsin to reveal the structural basis for proton transfer with respect to protein internal hydration. A complex series of dynamic hydrogen-bonding networks involving water molecules exists, facilitated by water channels and hydration sites within proteorhodopsin. High levels of hydration were discovered at each proton transfer site-the retinal binding pocket and proton uptake and release sites-underscoring the critical participation of water molecules in the proton-pumping mechanism. Water-bridged interactions and local water channels were also observed and can potentially mediate long-distance proton transfer between each site. The most significant phenomenon is after isomerization of retinal, an increase in water flux occurs that connects the proton release group, a conserved arginine residue, and the retinal binding pocket. Our results provide a detailed description of the internal hydration of the early photointermediates in the proteorhodopsin photocycle under alkaline pH conditions. These results lay the fundamental groundwork for understanding the intimate role that hydration plays in the structure-function relationship underlying the proteorhodopsin proton-pumping mechanism, as well as providing context for the relationship of hydration in proteorhodopsin to other microbial retinal proteins.

Quantum effects in a simple ring with hydrogen bonds

Complexes containing multiple arginines are common in proteins. The arginines are typically salt-bridged or hydrogen-bonded, so that their charges do not repel. Here we present a quantum calculation of a ring in which the components of a salt bridge composed of a guanidinium, the arginine side chain, and a carboxylic acid are separated by water molecules. When one water molecule is displaced from the ring, atomic charges of the other water molecule, as well as other properties, are significantly affected. The exchange and correlation energy differences between optimized and displaced rings are larger than thermal energy at room temperature, and larger than the sum of other energy differences. This suggests that calculations on proteins and other systems where such a ring may occur must take quantum effects into account; charges on certain atoms shift as substituents are added to the system: another water molecule, an -OH, or -CN bonded to either moiety. Also, charge shifts accompany proton shifts from the acid to guanidinium to ionize the salt bridge. The consequences of moving one water out of the ring give evidence for electron delocalization. Bond order and atomic charges are determined using natural bond orbital calculations. The geometry of the complex changes with ionization as well as the -OH and -CN additions but not in a simple manner. These results help in understanding the role of groups of arginines in salt-bridged clusters in proteins.

Molecular simulation of water and hydration effects in different environments: challenges and developments for DFTB based models

We discuss the description of water and hydration effects that employs an approximate density functional theory, DFTB3, in either a full QM or QM/MM framework. The goal is to explore, with the current formulation of DFTB3, the performance of this method for treating water in different chemical environments, the magnitude and nature of changes required to improve its performance, and factors that dictate its applicability to reactions in the condensed phase in a QM/MM framework. A relatively minor change (on the scale of kBT) in the O-H repulsive potential is observed to substantially improve the structural properties of bulk water under ambient conditions; modest improvements are also seen in dynamic properties of bulk water. This simple change also improves the description of protonated water clusters, a solvated proton, and to a more limited degree, a solvated hydroxide. By comparing results from DFTB3 models that differ in the description of water, we confirm that proton transfer energetics are adequately described by the standard DFTB3/3OB model for meaningful mechanistic analyses. For QM/MM applications, a robust parametrization of QM-MM interactions requires an explicit consideration of condensed phase properties, for which an efficient sampling technique was developed recently and is reviewed here. The discussions help make clear the value and limitations of DFTB3 based simulations, as well as the developments needed to further improve the accuracy and transferability of the methodology.