1
|
Helmer N, Wolf S, Stock G. Energy Transport and Its Function in Heptahelical Transmembrane Proteins. J Phys Chem B 2022; 126:8735-8746. [PMID: 36261792 DOI: 10.1021/acs.jpcb.2c05892] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Photoproteins such as bacteriorhodopsin (bR) and rhodopsin (Rho) need to effectively dissipate photoinduced excess energy to prevent themselves from damage. Another well-studied seven transmembrane (TM) helices protein is the β2 adrenergic receptor (β2AR), a G protein-coupled receptor for which energy dissipation paths have been linked with allosteric communication. To study the vibrational energy transport in the active and inactive states of these proteins, a master equation approach [J. Chem. Phys.2020, 152, 045103] is employed, which uses scaling rules that allow us to calculate energy transport rates solely based on the protein structure. Despite their overall structural similarity, the three 7TM proteins reveal quite different strategies to redistribute excess energy. While bR quickly removes the energy using the TM7 helix as a "lightning rod", Rho exhibits a rather poor energy dissipation, which might eventually require the hydrolysis of the Schiff base between the protein and the retinal chromophore to prevent overheating. Heating the ligand adrenaline of β2AR, the resulting energy transport network of the protein is found to change significantly upon switching from the active state to the inactive state. While the energy flow may highlight aspects of the inter-residue couplings of β2AR, it seems not particularly suited to explain allosteric phenomena.
Collapse
Affiliation(s)
- Nadja Helmer
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104Freiburg, Germany
| | - Steffen Wolf
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104Freiburg, Germany
| | - Gerhard Stock
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104Freiburg, Germany
| |
Collapse
|
2
|
Domino effect in allosteric signaling of peptide binding. J Mol Biol 2022; 434:167661. [DOI: 10.1016/j.jmb.2022.167661] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2022] [Revised: 05/24/2022] [Accepted: 05/24/2022] [Indexed: 11/22/2022]
|
3
|
Yamato T, Wang T, Sugiura W, Laprévote O, Katagiri T. Computational Study on the Thermal Conductivity of a Protein. J Phys Chem B 2022; 126:3029-3036. [PMID: 35416670 DOI: 10.1021/acs.jpcb.2c00958] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Protein molecules are thermally fluctuating and tightly packed amino acid residues strongly interact with each other. Such interactions are characterized in terms of heat current at the atomic level. We calculated the thermal conductivity of a small globular protein, villin headpiece subdomain, based on the linear response theory using equilibrium molecular dynamics simulation. The value of its thermal conductivity was 0.3 ± 0.01 [W m-1 K-1], which is in good agreement with experimental and computational studies on the other proteins in the literature. Heat current along the main chain was dominated by local vibrations in the polypeptide bonds, with amide I, II, III, and A bands on the Fourier transform of the heat current autocorrelation function.
Collapse
Affiliation(s)
- Takahisa Yamato
- Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan
| | - Tingting Wang
- Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan
| | - Wataru Sugiura
- Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan
| | - Olivier Laprévote
- Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan
| | - Takahiro Katagiri
- Information Technology Center, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan
| |
Collapse
|
4
|
Gong W, Liu Y, Zhao Y, Wang S, Han Z, Li C. Equally Weighted Multiscale Elastic Network Model and Its Comparison with Traditional and Parameter-Free Models. J Chem Inf Model 2021; 61:921-937. [PMID: 33496590 DOI: 10.1021/acs.jcim.0c01178] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Dynamical properties of proteins play an essential role in their function exertion. The elastic network model (ENM) is an effective and efficient tool in characterizing the intrinsic dynamical properties encoded in biomacromolecule structures. The Gaussian network model (GNM) and anisotropic network model (ANM) are the two often-used ENM models. Here, we introduce an equally weighted multiscale ENM (equally weighted mENM) based on the original mENM (denoted as mENM), in which fitting weights of Kirchhoff/Hessian matrixes in mENM are removed since they neglect the details of pairwise interactions. Then, we perform its comparison with the mENM, traditional ENM, and parameter-free ENM (pfENM) in reproducing dynamical properties for the six representative proteins whose molecular dynamics (MD) trajectories are available in http://mmb.pcb.ub.es/MoDEL/. In the results, for B-factor prediction, mENM performs best, while the equally weighted mENM performs also well, better than the traditional ENM and pfENM models. As to the dynamical cross-correlation map calculation, mENM performs worst, while the results produced from the equally weighted mENM and pfENM models are close to those from MD trajectories with the latter a little better than the former. Furthermore, encouragingly, the equally weighted mANM displays the best performance in capturing the functional motional modes, followed by pfANM and traditional ANM models, while the mANM fails in all the cases. This work is helpful for strengthening the understanding of the elastic network model and provides a valuable guide for researchers to utilize the model to explore protein dynamics.
Collapse
Affiliation(s)
- Weikang Gong
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Yang Liu
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Yanpeng Zhao
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Shihao Wang
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Zhongjie Han
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| | - Chunhua Li
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing 100124, China.,Beijing International Science and Technology Cooperation Base for Intelligent Physiological Measurement and Clinical Transformation, Beijing University of Technology, Beijing 100124, China
| |
Collapse
|
5
|
Poudel H, Reid KM, Yamato T, Leitner DM. Energy Transfer across Nonpolar and Polar Contacts in Proteins: Role of Contact Fluctuations. J Phys Chem B 2020; 124:9852-9861. [PMID: 33107736 DOI: 10.1021/acs.jpcb.0c08091] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Molecular dynamics simulations of the villin headpiece subdomain HP36 have been carried out to examine relations between rates of vibrational energy transfer across non-covalently bonded contacts and equilibrium structural fluctuations, with focus on van der Waals contacts. Rates of energy transfer across van der Waals contacts vary inversely with the variance of the contact length, with the same constant of proportionality for all nonpolar contacts of HP36. A similar relation is observed for hydrogen bonds, but the proportionality depends on contact pairs, with hydrogen bonds stabilizing the α-helices all exhibiting the same constant of proportionality, one that is distinct from those computed for other polar contacts. Rates of energy transfer across van der Waals contacts are found to be up to 2 orders of magnitude smaller than rates of energy transfer across polar contacts.
Collapse
Affiliation(s)
- Humanath Poudel
- Department of Chemistry, University of Nevada, Reno, Nevada 89557, United States
| | - Korey M Reid
- Department of Chemistry, University of Nevada, Reno, Nevada 89557, United States
| | - Takahisa Yamato
- Graduate School of Science, Division of Material Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan
| | - David M Leitner
- Department of Chemistry, University of Nevada, Reno, Nevada 89557, United States
| |
Collapse
|
6
|
Kadantsev VN, Goltsov A. Collective excitations in α-helical protein structures interacting with the water environment. Electromagn Biol Med 2020; 39:419-432. [PMID: 33023315 DOI: 10.1080/15368378.2020.1826961] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Abstract
Low-frequency vibrational excitations of protein macromolecules in the terahertz frequency region are suggested to contribute to many biological processes such as enzymatic catalysis, intra-protein energy/charge transport, recognition, and allostery. To explain high effectiveness of these processes, two possible mechanisms of the long-lived excitation were proposed by H. Fröhlich and A.S. Davydov, which relate to either vibrational modes or solitary waves, respectively. In this paper, we developed a quantum dynamic model of vibrational excitation in α-helical proteins interacting with the aqueous environment. In the model, we distinguished three coupled subsystems, i.e., (i) a chain of hydrogen-bonded peptide groups (PGs), interacting with (ii) the subsystem of the side-chain residuals which in turn interact with (iii) the environment, surrounding water responsible for dissipation and fluctuation in the system. It was shown that the equation of motion for phonon variables of the PG chain can be transformed to nonlinear Schrodinger equation which admits bifurcation into the solution corresponding to the weak-damped vibrational modes (Fröhlich-type regime) and Davydov solitons. A bifurcation parameter is derived through the strength of phonon-phonon interaction between the side-chains and hydration-shell water molecules. As shown, the energy of these excited states is pumped through the interaction of the side-chains with fluctuating water environment of the proteins. The suggested mechanism of the collective vibrational mode excitation is discussed in connection with the recent experiments on the long-lived collective protein excitations in the terahertz frequency region and vibrational energy transport pathways in proteins.
Collapse
Affiliation(s)
| | - Alexey Goltsov
- Russian Technological University (MIREA) , Moscow, Russia
| |
Collapse
|
7
|
Drug design by machine-trained elastic networks: predicting Ser/Thr-protein kinase inhibitors' activities. Mol Divers 2020; 25:899-909. [PMID: 32222890 DOI: 10.1007/s11030-020-10074-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2020] [Accepted: 03/11/2020] [Indexed: 12/23/2022]
Abstract
An elastic network model (ENM) represents a molecule as a matrix of pairwise atomic interactions. Rich in coded information, ENMs are hereby proposed as a novel tool for the prediction of the activity of series of molecules, with widely different chemical structures, but a common biological activity. The new approach is developed and tested using a set of 183 inhibitors of serine/threonine-protein kinase enzyme (Plk3) which is an enzyme implicated in the regulation of cell cycle and tumorigenesis. The elastic network (EN) predictive model is found to exhibit high accuracy and speed compared to descriptor-based machine-trained modeling. EN modeling appears to be a highly promising new tool for the high demands of industrial applications such as drug and material design.
Collapse
|
8
|
Valiño Borau L, Gulzar A, Stock G. Master equation model to predict energy transport pathways in proteins. J Chem Phys 2020; 152:045103. [DOI: 10.1063/1.5140070] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Luis Valiño Borau
- Biomolecular Dynamics, Institute of Physics, Albert Ludwigs University, 79104 Freiburg, Germany
| | - Adnan Gulzar
- Biomolecular Dynamics, Institute of Physics, Albert Ludwigs University, 79104 Freiburg, Germany
| | - Gerhard Stock
- Biomolecular Dynamics, Institute of Physics, Albert Ludwigs University, 79104 Freiburg, Germany
| |
Collapse
|
9
|
Leitner DM, Pandey HD, Reid KM. Energy Transport across Interfaces in Biomolecular Systems. J Phys Chem B 2019; 123:9507-9524. [DOI: 10.1021/acs.jpcb.9b07086] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- David M. Leitner
- Department of Chemistry, University of Nevada, Reno, Nevada 89557, United States
| | - Hari Datt Pandey
- Department of Chemistry, University of Nevada, Reno, Nevada 89557, United States
| | - Korey M. Reid
- Department of Chemistry, University of Nevada, Reno, Nevada 89557, United States
| |
Collapse
|
10
|
Gulzar A, Valiño Borau L, Buchenberg S, Wolf S, Stock G. Energy Transport Pathways in Proteins: A Non-equilibrium Molecular Dynamics Simulation Study. J Chem Theory Comput 2019; 15:5750-5757. [PMID: 31433644 DOI: 10.1021/acs.jctc.9b00598] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
To facilitate the observation of biomolecular energy transport in real time and with single-residue resolution, recent experiments by Baumann et al. ( Angew. Chem. Int. Ed. 2019 , 58 , 2899 , DOI: 10.1002/anie.201812995 ) have used unnatural amino acids β-(1-azulenyl)alanine (Azu) and azidohomoalanine (Aha) to site-specifically inject and probe vibrational energy in proteins. To aid the interpretation of such experiments, non-equilibrium molecular dynamics simulations of the anisotropic energy flow in proteins TrpZip2 and PDZ3 domains are presented. On this account, an efficient simulation protocol is established that accurately mimics the excitation and probing steps of Azu and Aha. The simulations quantitatively reproduce the experimentally found cooling times of the solvated proteins at room temperature and predict that the cooling slows by a factor 2 below the glass temperature of water. In PDZ3, vibrational energy is shown to travel from the initially excited peptide ligand via a complex network of inter-residue contacts and backbone transport to distal regions of the protein. The supposed connection of these energy transport pathways with pathways of allosteric communication is discussed.
Collapse
Affiliation(s)
- Adnan Gulzar
- Biomolecular Dynamics, Institute of Physics , Albert Ludwigs University , 79104 Freiburg , Germany
| | - Luis Valiño Borau
- Biomolecular Dynamics, Institute of Physics , Albert Ludwigs University , 79104 Freiburg , Germany
| | - Sebastian Buchenberg
- Biomolecular Dynamics, Institute of Physics , Albert Ludwigs University , 79104 Freiburg , Germany
| | - Steffen Wolf
- Biomolecular Dynamics, Institute of Physics , Albert Ludwigs University , 79104 Freiburg , Germany
| | - Gerhard Stock
- Biomolecular Dynamics, Institute of Physics , Albert Ludwigs University , 79104 Freiburg , Germany
| |
Collapse
|
11
|
Reid KM, Yamato T, Leitner DM. Scaling of Rates of Vibrational Energy Transfer in Proteins with Equilibrium Dynamics and Entropy. J Phys Chem B 2018; 122:9331-9339. [DOI: 10.1021/acs.jpcb.8b07552] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Korey M. Reid
- Department of Chemistry, University of Nevada, Reno, Nevada 89557, United States
| | - Takahisa Yamato
- Graduate School of Science, Division of Material Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan
| | - David M. Leitner
- Department of Chemistry, University of Nevada, Reno, Nevada 89557, United States
| |
Collapse
|