Mechanism of a proton pump analyzed with computer simulations

Thermal isomerization rates in retinal analogues using Ab-Initio molecular dynamics

For a detailed understanding of chemical processes in nature and industry, we need accurate models of chemical reactions in complex environments. While Eyring transition state theory is commonly used for modeling chemical reactions, it is most accurate for small molecules in the gas phase. A wide range of alternative rate theories exist that can better capture reactions involving complex molecules and environmental effects. However, they require that the chemical reaction is sampled by molecular dynamics simulations. This is a formidable challenge since the accessible simulation timescales are many orders of magnitude smaller than typical timescales of chemical reactions. To overcome these limitations, rare event methods involving enhanced molecular dynamics sampling are employed. In this work, thermal isomerization of retinal is studied using tight-binding density functional theory. Results from transition state theory are compared to those obtained from enhanced sampling. Rates obtained from dynamical reweighting using infrequent metadynamics simulations were in close agreement with those from transition state theory. Meanwhile, rates obtained from application of Kramers' rate equation to a sampled free energy profile along a torsional dihedral reaction coordinate were found to be up to three orders of magnitude higher. This discrepancy raises concerns about applying rate methods to one-dimensional reaction coordinates in chemical reactions.

Tools for analyzing protonation states and for tracing proton transfer pathways with examples from the Rb. sphaeroides photosynthetic reaction centers

Protons participate in many reactions. In proteins, protons need paths to move in and out of buried active sites. The vectorial movement of protons coupled to electron transfer reactions establishes the transmembrane electrochemical gradient used for many reactions, including ATP synthesis. Protons move through hydrogen bonded chains of waters and hydroxy side chains via the Grotthuss mechanism and by proton binding and release from acidic and basic residues. MCCE analysis shows that proteins exist in a large number of protonation states. Knowledge of the equilibrium ensemble can provide a rational basis for setting protonation states in simulations that fix them, such as molecular dynamics (MD). The proton path into the Q_B site in the bacterial reaction centers (RCs) of Rb. sphaeroides is analyzed by MD to provide an example of the benefits of using protonation states found by the MCCE program. A tangled web of side chains and waters link the cytoplasm to Q_B. MCCE analysis of snapshots from multiple trajectories shows that changing the input protonation state of a residue in MD biases the trajectory shifting the proton affinity of that residue. However, the proton affinity of some residues is more sensitive to the input structure. The proton transfer networks derived from different trajectories are quite robust. There are some changes in connectivity that are largely restricted to the specific residues whose protonation state is changed. Trajectories with Q_B^•- are compared with earlier results obtained with Q_B [Wei et. al Photosynthesis Research volume 152, pages153-165 (2022)] showing only modest changes. While introducing new methods the study highlights the difficulty of establishing the connections between protein conformation.

Collective exchange processes reveal an active site proton cage in bacteriorhodopsin

Proton translocation across membranes is vital to all kingdoms of life. Mechanistically, it relies on characteristic proton flows and modifications of hydrogen bonding patterns, termed protonation dynamics, which can be directly observed by fast magic angle spinning (MAS) NMR. Here, we demonstrate that reversible proton displacement in the active site of bacteriorhodopsin already takes place in its equilibrated dark-state, providing new information on the underlying hydrogen exchange processes. In particular, MAS NMR reveals proton exchange at D85 and the retinal Schiff base, suggesting a tautomeric equilibrium and thus partial ionization of D85. We provide evidence for a proton cage and detect a preformed proton path between D85 and the proton shuttle R82. The protons at D96 and D85 exchange with water, in line with ab initio molecular dynamics simulations. We propose that retinal isomerization makes the observed proton exchange processes irreversible and delivers a proton towards the extracellular release site.

Daniel Friedrich et al. show that reversible proton translocation occurs in the dark–state of bacteriorhodopsin, involving the retinal Schiff base and D85 exchanging protons with H₂O. They find evidence of an active site proton cage and possible proton transfer via R82.

Primary Transfer Step in the Light-Driven Ion Pump Bacteriorhodopsin: An Irreversible U-Turn Revealed by Dynamic Nuclear Polarization-Enhanced Magic Angle Spinning NMR

Despite much attention, the path of the highly consequential primary proton transfer in the light-driven ion pump bacteriorhodopsin (bR) remains mysterious. Here we use DNP-enhanced magic angle spinning (MAS) NMR to study critical elements of the active site just before the Schiff base (SB) deprotonates (in the L intermediate), immediately after the SB has deprotonated and Asp85 has become protonated (in the M_o intermediate), and just after the SB has reprotonated and Asp96 has deprotonated (in the N intermediate). An essential feature that made these experiments possible is the 75-fold signal enhancement through DNP. ¹⁵N(SB)-¹H correlations reveal that the newly deprotonated SB is accepting a hydrogen bond from an alcohol and ¹³C-¹³C correlations show that Asp85 draws close to Thr89 before the primary proton transfer. Concurrently, ¹⁵N-¹³C correlations between the SB and Asp85 show that helices C and G draw closer together just prior to the proton transfer and relax thereafter. Together, these results indicate that Thr89 serves to relay the SB proton to Asp85 and that creating this pathway involves rapprochement between the C and G helices as well as chromophore torsion.

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.

Density functional tight binding: values of semi-empirical methods in an ab initio era

Semi-empirical (SE) methods are derived from Hartree-Fock (HF) or Density Functional Theory (DFT) by neglect and approximation of electronic integrals. Thereby, parameters are introduced which have to be determined from reference calculations and/or by fitting to available experimental data. This leads to computational methods that are about 2-3 orders of magnitude faster than the standard HF/DFT methods using medium sized basis sets while being about 3 orders of magnitude slower than empirical force field methods (Molecular Mechanics: MM). Therefore, SE methods are most appropriate for a specific range of applications. These include the study of systems that contain a large number of atoms and therefore being too large for ab initio or DFT methods and also problems where dynamic or entropic effects are particularly important. In the latter case, the errors made by considering a very limited number of molecular structures or neglecting entropic contributions can be much larger than the accuracy lost due to the use of SE methods. Another area where SE methods are attractive concerns the analysis of systems for which reliable MM models are not readily available. Therefore, even in an era when rapid progress is being made in ab initio methods, there is considerable interest in further developing SE methods. We illustrate this point by focusing on the discussion of recent development and application of the Density Functional Tight Binding method.

A modified QM/MM Hamiltonian with the Self-Consistent-Charge Density-Functional-Tight-Binding Theory for highly charged QM regions

To improve the description of electrostatic interaction between QM and MM atoms when the QM is SCC-DFTB, we adopt a Klopman-Ohno (KO) functional form which considers the finite size of the QM and MM charge distributions. Compared to the original implementation that used a simple Coulombic interaction between QM Mulliken and MM point charges, the KO based QM/MM scheme takes charge penetration effect into consideration and therefore significantly improves the description of QM/MM interaction at short range, especially when the QM region is highly charged. To be consistent with the third-order formulation of SCC-DFTB, the Hubbard parameter in the KO functional is dependent on the QM charge. As a result, the effective size of the QM charge distribution naturally adjusts as the QM region undergoes chemical transformations, making the KO based QM/MM scheme particularly attractive for describing chemical reactions in the condensed phase. Together with the van der Waals parameters for the QM atom, the KO based QM/MM model introduces four parameters for each element type. They are fitted here based on microsolvation models of small solutes, focusing on negatively charged molecular ions, for elements O, C, H and P with a specific version of SCC-DFTB (SCC-DFTBPR). Test calculations confirm that the KO based QM/MM scheme significantly improves the interactions between QM and MM atoms over the original point charge based model and it is transferable due to the small number of parameters. The new form of QM/MM Hamiltonian will greatly improve the applicability of SCC-DFTB based QM/MM methods to problems that involve highly charged QM regions, such as enzyme catalyzed phosphoryl transfers.

Coupling of retinal, protein, and water dynamics in squid rhodopsin

The light-induced isomerization of the retinal from 11-cis to all-trans triggers changes in the conformation of visual rhodopsins that lead to the formation of the activated state, which is ready to interact with the G protein. To begin to understand how changes in the structure and dynamics of the retinal are transmitted to the protein, we performed molecular dynamics simulations of squid rhodopsin with 11-cis and all-trans retinal, and with two different force fields for describing the retinal molecule. The results indicate that structural rearrangements in the binding pocket, albeit small, propagate toward the cytoplasmic side of the protein, and affect the dynamics of internal water molecules. The sensitivity of the active-site interactions on the retinal force-field parameters highlights the coupling between the retinal molecule and its immediate protein environment.