1
|
Zhu Q, Soudackov AV, Tommos C, Hammes-Schiffer S. Proton-Coupled Electron Transfer upon Oxidation of Tyrosine in a De Novo Protein: Analysis of Proton Acceptor Candidates. Biochemistry 2024; 63:1999-2008. [PMID: 39024184 PMCID: PMC11402026 DOI: 10.1021/acs.biochem.4c00211] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2024]
Abstract
Redox-active residues, such as tyrosine and tryptophan, play important roles in a wide range of biological processes. The α3Y de novo protein, which is composed of three α helices and a tyrosine residue Y32, provides a platform for investigating the redox properties of tyrosine in a well-defined protein environment. Herein, the proton-coupled electron transfer (PCET) reaction that occurs upon oxidation of tyrosine in this model protein by a ruthenium photosensitizer is studied by using a vibronically nonadiabatic PCET theory that includes hydrogen tunneling and excited vibronic states. The input quantities to the analytical nonadiabatic rate constant expression, such as the diabatic proton potential energy curves and associated proton vibrational wave functions, reorganization energy, and proton donor-acceptor distribution functions, are obtained from density functional theory calculations on model systems and molecular dynamics simulations of the solvated α3Y protein. Two possible proton acceptors, namely, water or a glutamate residue in the protein scaffold, are explored. The PCET rate constant is greater when glutamate is the proton acceptor, mainly due to the more favorable driving force and shorter equilibrium proton donor-acceptor distance, although contributions from excited vibronic states mitigate these effects. Nevertheless, water could be the dominant proton acceptor if its equilibrium constant associated with hydrogen bond formation is significantly greater than that for glutamate. Although these calculations do not definitively identify the proton acceptor for this PCET reaction, they elucidate the conditions under which each proton acceptor can be favored. These insights have implications for tyrosine-based PCET in a wide variety of biochemical processes.
Collapse
Affiliation(s)
- Qiwen Zhu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Alexander V Soudackov
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Cecilia Tommos
- Department of Biochemistry and Biophysics, Texas A&M University, College Station, Texas 77843, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
2
|
Weckel-Dahman H, Carlsen R, Swanson JM. Multiscale Responsive Kinetic Modeling: Quantifying Biomolecular Reaction Flux under Varying Electrochemical Conditions. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.08.01.606205. [PMID: 39131358 PMCID: PMC11312519 DOI: 10.1101/2024.08.01.606205] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Indexed: 08/13/2024]
Abstract
Attaining a complete thermodynamic and kinetic characterization for processes involving multiple interconnected rare-event transitions remains a central challenge in molecular biophysics. This challenge is amplified when the process must be understood under a range of reaction conditions. Herein, we present a condition-responsive kinetic modeling framework that can combine the strengths of bottom-up rate quantification from multiscale simulations with top-down solution refinement using experimental data. Although this framework can be applied to any process, we demonstrate its use for electrochemically driven transport through channels and transporters. Using the Cl- /H+ antiporter ClC-ec1 as a model system, we show how robust and predictive kinetic solutions can be obtained when the solution space is grounded by thermodynamic constraints, seeded through multiscale rate quantification, and further refined with experimental data, such as electrophysiology assays. Turning to the Shaker K+ channel, we demonstrate that robust solutions and biophysical insights can also be obtained with sufficient experimental data. This multi-pathway method proves capable of identifying single-pathway dominant mechanisms but also highlights that competing and off-pathway flux is still essential to replicate experimental findings and to describe concentration-dependent channel rectification.
Collapse
Affiliation(s)
- Hannah Weckel-Dahman
- Department of Chemistry, University of Utah, Salt Lake City, UT, 84112 – United States of America
| | - Ryan Carlsen
- Department of Chemistry, University of Utah, Salt Lake City, UT, 84112 – United States of America
| | - Jessica M.J. Swanson
- Department of Chemistry, University of Utah, Salt Lake City, UT, 84112 – United States of America
| |
Collapse
|
3
|
Liu Y, Li C, Gupta M, Stroud RM, Voth GA. Kinetic network modeling with molecular simulation inputs: A proton-coupled phosphate symporter. Biophys J 2024:S0006-3495(24)00216-9. [PMID: 38549372 DOI: 10.1016/j.bpj.2024.03.035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2024] [Revised: 03/14/2024] [Accepted: 03/26/2024] [Indexed: 04/13/2024] Open
Abstract
Phosphate, an essential metabolite involved in numerous cellular functions, is taken up by proton-coupled phosphate transporters of plants and fungi within the major facilitator family. Similar phosphate transporters have been identified across a diverse range of biological entities, including various protozoan parasites linked to human diseases, breast cancer cells with increased phosphate requirements, and osteoclast-like cells engaged in bone resorption. Prior studies have proposed an overview of the functional cycle of a proton-driven phosphate transporter (PiPT), yet a comprehensive understanding of the proposed reaction pathways necessitates a closer examination of each elementary reaction step within an overall kinetic framework. In this work, we leverage kinetic network modeling in conjunction with a "bottom-up" molecular dynamics approach to show how such an approach can characterize the proton-phosphate co-transport behavior of PiPT under different pH and phosphate concentration conditions. In turn, this allows us to reveal the prevailing reaction pathway within a high-affinity phosphate transporter under different experimental conditions and to uncover the molecular origin of the optimal pH condition of this transporter.
Collapse
Affiliation(s)
- Yu Liu
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois
| | - Chenghan Li
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois
| | - Meghna Gupta
- Department of Biochemistry and Biophysics, University of California, San Francisco, San Francisco, California
| | - Robert M Stroud
- Department of Biochemistry and Biophysics, University of California, San Francisco, San Francisco, California
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois.
| |
Collapse
|
4
|
George A, Zuckerman DM. From Average Transient Transporter Currents to Microscopic Mechanism─A Bayesian Analysis. J Phys Chem B 2024; 128:1830-1842. [PMID: 38373358 DOI: 10.1021/acs.jpcb.3c07025] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/21/2024]
Abstract
Electrophysiology studies of secondary active transporters have revealed quantitative mechanistic insights over many decades of research. However, the emergence of new experimental and analytical approaches calls for investigation of the capabilities and limitations of the newer methods. We examine the ability of solid-supported membrane electrophysiology (SSME) to characterize discrete-state kinetic models with >10 rate constants. We use a Bayesian framework applied to synthetic data for three tasks: to quantify and check (i) the precision of parameter estimates under different assumptions, (ii) the ability of computation to guide the selection of experimental conditions, and (iii) the ability of our approach to distinguish among mechanisms based on SSME data. When the general mechanism, i.e., event order, is known in advance, we show that a subset of kinetic parameters can be "practically identified" within ∼1 order of magnitude, based on SSME current traces that visually appear to exhibit simple exponential behavior. This remains true even when accounting for systematic measurement bias and realistic uncertainties in experimental inputs (concentrations) are incorporated into the analysis. When experimental conditions are optimized or different experiments are combined, the number of practically identifiable parameters can be increased substantially. Some parameters remain intrinsically difficult to estimate through SSME data alone, suggesting that additional experiments are required to fully characterize parameters. We also demonstrate the ability to perform model selection and determine the order of events when that is not known in advance, comparing Bayesian and maximum-likelihood approaches. Finally, our studies elucidate good practices for the increasingly popular but subtly challenging Bayesian calculations for structural and systems biology.
Collapse
Affiliation(s)
- August George
- Department of Biomedical Engineering, Oregon Health and Science University, Portland, Oregon 97239, United States
| | - Daniel M Zuckerman
- Department of Biomedical Engineering, Oregon Health and Science University, Portland, Oregon 97239, United States
| |
Collapse
|
5
|
Yue Z, Li C, Voth GA. The role of conformational change and key glutamic acid residues in the ClC-ec1 antiporter. Biophys J 2023; 122:1068-1085. [PMID: 36698313 PMCID: PMC10111279 DOI: 10.1016/j.bpj.2023.01.025] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2022] [Revised: 01/16/2023] [Accepted: 01/20/2023] [Indexed: 01/26/2023] Open
Abstract
The triple glutamine (Q) mutant (QQQ) structure of a Cl-/H+ antiporter from Escherichia coli (ClC-ec1) displaying a novel backbone arrangement has been used to challenge the long-held notion that Cl-/H+ antiporters do not operate through large conformational motions. The QQQ mutant substitutes the glutamine residue for an external glutamate E148, an internal glutamate E203, and a third glutamate E113 that hydrogen-bonds with E203. However, it is unknown if QQQ represents a physiologically relevant state, as well as how the protonation of the wild-type glutamates relates to the global dynamics. We herein apply continuous constant-pH molecular dynamics to investigate the H+-coupled dynamics of ClC-ec1. Although any large-scale conformational rearrangement upon acidification would be due to the accumulation of excess charge within the protein, protonation of the glutamates significantly impacts mainly the local structure and dynamics. Despite the fact that the extracellular pore enlarges at acidic pHs, an occluded ClC-ec1 within the active pH range of 3.5-7.5 requires a protonated E148 to facilitate extracellular Cl- release. E203 is also involved in the intracellular H+ transfer as an H+ acceptor. The water wire connection of E148 with the intracellular solution is regulated by the charge states of the E113/E203 dyad with coupled proton titration. However, the dynamics extracted from our simulations are not QQQ-like, indicating that the QQQ mutant does not represent the behavior of the wild-type ClC-ec1. These findings reinforce the necessity of having a protonatable residue at the E203 position in ClC-ec1 and suggest that a higher level of complexity exists for the intracellular H+ transfer in Cl-/H+ antiporters.
Collapse
Affiliation(s)
- Zhi Yue
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois
| | - Chenghan Li
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois.
| |
Collapse
|
6
|
Proton Coupling and the Multiscale Kinetic Mechanism of a Peptide Transporter. Biophys J 2022; 121:2266-2278. [PMID: 35614850 DOI: 10.1016/j.bpj.2022.05.029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2022] [Revised: 05/13/2022] [Accepted: 05/19/2022] [Indexed: 11/02/2022] Open
Abstract
Proton coupled peptide transporters (POTs) are crucial for the uptake of di- and tri-peptides as well as drug and pro-drug molecules in prokaryotes and eukaryotic cells. We illustrate from multiscale modeling how transmembrane proton flux couples within a POT protein to drive essential steps of the full functional cycle: 1) protonation of a glutamate on transmembrane helix (TM) 7 opens the extracellular gate, allowing ligand entry; 2) inward proton flow induces the cytosolic release of ligand by varying the protonation state of a second conserved glutamate on TM10; 3) proton movement between TM7 and TM10 is thermodynamically driven and kinetically permissible via water proton shuttling without the participation of ligand. Our results, for the first time, give direct computational confirmation for the alternating access model of POTs, and point to a quantitative multiscale kinetic picture of the functioning protein mechanism.
Collapse
|
7
|
Swanson JMJ. Multiscale kinetic analysis of proteins. Curr Opin Struct Biol 2022; 72:169-175. [PMID: 34923310 PMCID: PMC9288830 DOI: 10.1016/j.sbi.2021.11.005] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2021] [Revised: 10/30/2021] [Accepted: 11/01/2021] [Indexed: 02/03/2023]
Abstract
The stochasticity of molecular motion results in the existence of multiple kinetically relevant pathways in many biomolecular mechanisms. Because it is highly demanding to characterize them for complex systems, mechanisms are often described with a single-pathway perspective. However, kinetic network analysis and sub-ensemble experimental insight are increasingly demonstrating not only the existence of competing pathways but also the importance of kinetic selection in biology. This review focuses on advances in multiscale kinetic analysis of proteins, which connects molecular level information from simulations to macroscopic data to characterize mechanistic reaction networks and the reactive flux through them. We describe a range of methods used and highlight several examples where kinetic modeling has revealed functional importance of pathway heterogeneity.
Collapse
|