1
|
Jiang W. Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics. J Chem Theory Comput 2024; 20:3085-3095. [PMID: 38568961 DOI: 10.1021/acs.jctc.4c00191] [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: 04/05/2024]
Abstract
Ligand binding free energy simulations (LB-FES) that involve sampling of protein functional conformations have been longstanding challenges in research on molecular recognition. Particularly, modeling of the conformational transition pathway and design of the heuristic biasing mechanism are severe bottlenecks for the existing enhanced configurational sampling (ECS) methods. Inspired by the key role of hydration in regulating conformational dynamics of macromolecules, this report proposes a novel ECS approach that facilitates binding-associated structural dynamics by accelerated hydration transitions in combination with the λ-exchange of free energy perturbation (FEP). Two challenging protein-ligand binding processes involving large configurational transitions of the receptor are studied, with hydration transitions at binding sites accelerated by Hamiltonian-simulated annealing of the hydration layer. Without the need for pathway analysis or ad hoc barrier flattening potential, LB-FES were performed with FEP/λ-exchange molecular dynamics simulation at a minor overhead for annealing of the hydration layer. The LB-FES studies showed that the accelerated rehydration significantly enhances the collective conformational transitions of the receptor, and convergence of binding affinity calculations is obtained at a sweet-spot simulation time scale. Alchemical LB-FES with the proposed ECS strategy is free from the effort of trial and error for the setup and realizes efficient on-the-fly sampling for the collective functional response of the receptor and bound water and therefore presents a practical approach to high-throughput screening in drug discovery.
Collapse
Affiliation(s)
- Wei Jiang
- Computational Science Division, Argonne National Laboratory, 9700 South Cass Avenue, Building 240, Argonne, Illinois 60439, United States
| |
Collapse
|
2
|
Ohler A, Taylor PE, Bledsoe JA, Iavarone AT, Gilbert NC, Offenbacher AR. Identification of the Thermal Activation Network in Human 15-Lipoxygenase-2: Divergence from Plant Orthologs and Its Relationship to Hydrogen Tunneling Activation Barriers. ACS Catal 2024; 14:5444-5457. [PMID: 38601784 PMCID: PMC11003420 DOI: 10.1021/acscatal.4c00439] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2024] [Revised: 03/05/2024] [Accepted: 03/13/2024] [Indexed: 04/12/2024]
Abstract
The oxidation of polyunsaturated fatty acids by lipoxygenases (LOXs) is initiated by a C-H cleavage step in which the hydrogen atom is transferred quantum mechanically (i.e., via tunneling). In these reactions, protein thermal motions facilitate the conversion of ground-state enzyme-substrate complexes to tunneling-ready configurations and are thus important for transferring energy from the solvent to the active site for the activation of catalysis. In this report, we employed temperature-dependent hydrogen-deuterium exchange mass spectrometry (TDHDX-MS) to identify catalytically linked, thermally activated peptides in a representative animal LOX, human epithelial 15-LOX-2. TDHDX-MS of wild-type 15-LOX-2 was compared to two active site mutations that retain structural stability but have increased activation energies (Ea) of catalysis. The Ea value of one variant, V427L, is implicated to arise from suboptimal substrate positioning by increased active-site side chain rotamer dynamics, as determined by X-ray crystallography and ensemble refinement. The resolved thermal network from the comparative Eas of TDHDX-MS between wild-type and V426A is localized along the front face of the 15-LOX-2 catalytic domain. The network contains a clustering of isoleucine, leucine, and valine side chains within the helical peptides. This thermal network of 15-LOX-2 is different in location, area, and backbone structure compared to a model plant lipoxygenase from soybean that exhibits a low Ea value of catalysis compared to the human ortholog. The presented data provide insights into the divergence of thermally activated protein motions in plant and animal LOXs and their relationships to the enthalpic barriers for facilitating hydrogen tunneling.
Collapse
Affiliation(s)
- Amanda Ohler
- Department
of Chemistry, East Carolina University, Greenville, North Carolina 27858, United States
| | - Paris E. Taylor
- Department
of Biological Sciences, Louisiana State
University, Baton
Rouge, Louisiana 70803, United States
| | - Jasmine A. Bledsoe
- Department
of Biological Sciences, Louisiana State
University, Baton
Rouge, Louisiana 70803, United States
| | - Anthony T. Iavarone
- QB3/Chemistry
Mass Spectrometry Facility, University of
California, Berkeley, Berkeley, California 94720, United States
| | - Nathaniel C. Gilbert
- Department
of Biological Sciences, Louisiana State
University, Baton
Rouge, Louisiana 70803, United States
| | - Adam R. Offenbacher
- Department
of Chemistry, East Carolina University, Greenville, North Carolina 27858, United States
| |
Collapse
|
3
|
Feid C, Luma L, Fischer T, Löffler JG, Grebenovsky N, Wachtveitl J, Heckel A, Bredenbeck J. Iminothioindoxyl Donors with Exceptionally High Cross Section for Protein Vibrational Energy Transfer. Angew Chem Int Ed Engl 2024; 63:e202317047. [PMID: 38103205 DOI: 10.1002/anie.202317047] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2023] [Revised: 12/10/2023] [Accepted: 12/15/2023] [Indexed: 12/18/2023]
Abstract
Various protein functions are related to vibrational energy transfer (VET) as an important mechanism. The underlying transfer pathways can be experimentally followed by ultrafast Vis-pump/IR-probe spectroscopy with a donor-sensor pair of non-canonical amino acids (ncAAs) incorporated in a protein. However, so far only one donor ncAA, azulenylalanine (AzAla), exists, which suffers from a comparably low Vis extinction coefficient. Here, we introduce two novel donor ncAAs based on an iminothioindoxyl (ITI) chromophore. The dimethylamino-ITI (DMA-ITI) and julolidine-ITI (J-ITI) moieties overcome the limitation of AzAla with a 50 times higher Vis extinction coefficient. While ITI moieties are known for ultrafast photoswitching, DMA-ITI and J-ITI exclusively form a hot ground state on the sub-ps timescale instead, which is essential for their usage as vibrational energy donor. In VET measurements of donor-sensor dipeptides we investigate the performance of the new donors. We observe 20 times larger signals compared to the established AzAla donor, which opens unprecedented possibilities for the study of VET in proteins.
Collapse
Affiliation(s)
- Carolin Feid
- Institute of Biophysics, Goethe University Frankfurt, Max-von-Laue-Straße 1, 60438, Frankfurt (Main), Germany
| | - Larita Luma
- Institute for Organic Chemistry and Chemical Biology, Goethe University Frankfurt, Max-von-Laue-Straße 7, 60438, Frankfurt (Main), Germany
| | - Tobias Fischer
- Institute for Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Straße 7, 60438, Frankfurt (Main), Germany
| | - Jan Gerrit Löffler
- Institute of Biophysics, Goethe University Frankfurt, Max-von-Laue-Straße 1, 60438, Frankfurt (Main), Germany
| | - Nikolai Grebenovsky
- Institute for Organic Chemistry and Chemical Biology, Goethe University Frankfurt, Max-von-Laue-Straße 7, 60438, Frankfurt (Main), Germany
| | - Josef Wachtveitl
- Institute for Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Straße 7, 60438, Frankfurt (Main), Germany
| | - Alexander Heckel
- Institute for Organic Chemistry and Chemical Biology, Goethe University Frankfurt, Max-von-Laue-Straße 7, 60438, Frankfurt (Main), Germany
| | - Jens Bredenbeck
- Institute of Biophysics, Goethe University Frankfurt, Max-von-Laue-Straße 1, 60438, Frankfurt (Main), Germany
| |
Collapse
|
4
|
Nam K, Shao Y, Major DT, Wolf-Watz M. Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development. ACS OMEGA 2024; 9:7393-7412. [PMID: 38405524 PMCID: PMC10883025 DOI: 10.1021/acsomega.3c09084] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/14/2023] [Revised: 01/15/2024] [Accepted: 01/19/2024] [Indexed: 02/27/2024]
Abstract
Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.
Collapse
Affiliation(s)
- Kwangho Nam
- Department
of Chemistry and Biochemistry, University
of Texas at Arlington, Arlington, Texas 76019, United States
| | - Yihan Shao
- Department
of Chemistry and Biochemistry, University
of Oklahoma, Norman, Oklahoma 73019-5251, United States
| | - Dan T. Major
- Department
of Chemistry and Institute for Nanotechnology & Advanced Materials, Bar-Ilan University, Ramat-Gan 52900, Israel
| | | |
Collapse
|
5
|
Hill TD, Basnet S, Lepird HH, Rightnowar BW, Moran SD. Anisotropic dynamics of an interfacial enzyme active site observed using tethered substrate analogs and ultrafast 2D IR spectroscopy. J Chem Phys 2023; 159:165101. [PMID: 37870142 PMCID: PMC10597647 DOI: 10.1063/5.0167991] [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: 07/16/2023] [Accepted: 09/29/2023] [Indexed: 10/24/2023] Open
Abstract
Enzymes accelerate the rates of biomolecular reactions by many orders of magnitude compared to bulk solution, and it is widely understood that this catalytic effect arises from a combination of polar pre-organization and electrostatic transition state stabilization. A number of recent reports have also implicated ultrafast (femtosecond-picosecond) timescale motions in enzymatic activity. However, complications arising from spatially-distributed disorder, the occurrence of multiple substrate binding modes, and the influence of hydration dynamics on solvent-exposed active sites still confound many experimental studies. Here we use ultrafast two-dimensional infrared (2D IR) spectroscopy and covalently-tethered substrate analogs to examine dynamical properties of the promiscuous Pyrococcus horikoshii ene-reductase (PhENR) active site in two binding configurations mimicking proposed "inactive" and "reactive" Michaelis complexes. Spectral diffusion measurements of aryl-nitrile substrate analogs reveal an end-to-end tradeoff between fast (sub-ps) and slow (>5 ps) motions. Fermi resonant aryl-azide analogs that sense interactions of coupled oscillators are described. Lineshape and quantum beat analyses of these probes reveal characteristics that correlate with aryl-nitrile frequency fluctuation correlation functions parameters, demonstrating that this anisotropy is an intrinsic property of the water-exposed active site, where countervailing gradients of fast dynamics and disorder in the reactant ground state are maintained near the hydration interface. Our results suggest several plausible factors leading to state-selective rate enhancement and promiscuity in PhENR. This study also highlights a strategy to detect perturbations to vibrational modes outside the transparent window of the mid-IR spectrum, which may be extended to other macromolecular systems.
Collapse
Affiliation(s)
| | - Sunil Basnet
- School of Chemical and Biomolecular Sciences, Southern Illinois University Carbondale, 1245 Lincoln Drive MC 4409, Carbondale, Illinois 62901, USA
| | - Hannah H. Lepird
- School of Chemical and Biomolecular Sciences, Southern Illinois University Carbondale, 1245 Lincoln Drive MC 4409, Carbondale, Illinois 62901, USA
| | - Blaze W. Rightnowar
- School of Chemical and Biomolecular Sciences, Southern Illinois University Carbondale, 1245 Lincoln Drive MC 4409, Carbondale, Illinois 62901, USA
| | - Sean D. Moran
- School of Chemical and Biomolecular Sciences, Southern Illinois University Carbondale, 1245 Lincoln Drive MC 4409, Carbondale, Illinois 62901, USA
| |
Collapse
|
6
|
Conti Nibali V, Sacchetti F, Paciaroni A, Petrillo C, Tarek M, D'Angelo G. Intra-protein interacting collective modes in the terahertz frequency region. J Chem Phys 2023; 159:161101. [PMID: 37870134 DOI: 10.1063/5.0142381] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2023] [Accepted: 09/25/2023] [Indexed: 10/24/2023] Open
Abstract
Understanding how proteins work requires a thorough understanding of their internal dynamics. Proteins support a wide range of motions, from the femtoseconds to seconds time scale, relevant to crucial biological functions. In this context, the term "protein collective dynamics" refers to the complex patterns of coordinated motions of numerous atoms throughout the protein in the sub-picosecond time scale (terahertz frequency region). It is hypothesized that these dynamics have a substantial impact on the regulation of functional dynamical mechanisms, including ligand binding and allosteric signalling, charge transport direction, and the regulation of thermodynamic and thermal transport properties. Using the theoretical framework of hydrodynamics, the collective dynamics of proteins had previously been described in a manner akin to that of simple liquids, i.e. in terms of a single acoustic-like excitation, related to intra-protein vibrational motions. Here, we employ an interacting-mode model to analyse the results from molecular dynamics simulations and we unveil that the vibrational landscape of proteins is populated by multiple acoustic-like and low-frequency optic-like modes, with mixed symmetry and interfering with each other. We propose an interpretation at the molecular level of the observed scenario that we relate to the side-chains and the hydrogen-bonded networks dynamics. The present insights provide a perspective for understanding the molecular mechanisms underlying the energy redistribution processes in the interior of proteins.
Collapse
Affiliation(s)
- Valeria Conti Nibali
- Department of Mathematical and Computational Sciences, Physical Science and Earth Science, Messina University, Viale Ferdinando Stagno D'Alcontres 31, 98166 Messina, Italy
| | - Francesco Sacchetti
- Department of Physics and Geology, Perugia University, Via Alessandro Pascoli, I-06123 Perugia, Italy
| | - Alessandro Paciaroni
- Department of Physics and Geology, Perugia University, Via Alessandro Pascoli, I-06123 Perugia, Italy
| | - Caterina Petrillo
- Department of Physics and Geology, Perugia University, Via Alessandro Pascoli, I-06123 Perugia, Italy
| | - Mounir Tarek
- Université de Lorraine, CNRS, LPCT, F-54000 Nancy, France
| | - Giovanna D'Angelo
- Department of Mathematical and Computational Sciences, Physical Science and Earth Science, Messina University, Viale Ferdinando Stagno D'Alcontres 31, 98166 Messina, Italy
| |
Collapse
|
7
|
Azizi K, Gori M, Morzan U, Hassanali A, Kurian P. Examining the origins of observed terahertz modes from an optically pumped atomistic model protein in aqueous solution. PNAS NEXUS 2023; 2:pgad257. [PMID: 37575674 PMCID: PMC10416812 DOI: 10.1093/pnasnexus/pgad257] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/28/2022] [Revised: 07/14/2023] [Accepted: 07/26/2023] [Indexed: 08/15/2023]
Abstract
The microscopic origins of terahertz (THz) vibrational modes in biological systems are an active and open area of current research. Recent experiments [Phys Rev X. 8, 031061 (2018)] have revealed the presence of a pronounced mode at ∼0.3 THz in fluorophore-decorated bovine serum albumin (BSA) protein in aqueous solution under nonequilibrium conditions induced by optical pumping. This result was heuristically interpreted as a collective elastic fluctuation originating from the activation of a low-frequency phonon mode. In this work, we show that the sub-THz spectroscopic response emerges in a statistically significant manner (> 2 σ ) from such collective behavior, illustrating how photoexcitation can alter specific THz vibrational modes. We revisit the theoretical analysis with proof-of-concept molecular dynamics that introduce optical excitations into the simulations. Using information theory techniques, we show that these excitations can give rise to a multiscale response involving two optically excited chromophores (tryptophans), other amino acids in the protein, ions, and water. Our results motivate new experiments and fully nonequilibrium simulations to probe these phenomena, as well as the refinement of atomistic models of Fröhlich condensates that are fundamentally determined by nonlinear interactions in biology.
Collapse
Affiliation(s)
- Khatereh Azizi
- The Abdus Salam International Centre for Theoretical Physics, Strada Costiera 11, 34151 Trieste, Italy
- Quantum Biology Laboratory, Howard University, Washington, DC 20060, USA
| | - Matteo Gori
- Quantum Biology Laboratory, Howard University, Washington, DC 20060, USA
| | - Uriel Morzan
- The Abdus Salam International Centre for Theoretical Physics, Strada Costiera 11, 34151 Trieste, Italy
| | - Ali Hassanali
- The Abdus Salam International Centre for Theoretical Physics, Strada Costiera 11, 34151 Trieste, Italy
| | - Philip Kurian
- Quantum Biology Laboratory, Howard University, Washington, DC 20060, USA
| |
Collapse
|
8
|
Abstract
Debate has been simmering for some years regarding the importance of internal thermal motions of enzymes to catalysis. Recent developments in protein design may bring resolution of the more contentious points a little closer.
Collapse
Affiliation(s)
- Jeremy R. H. Tame
- Protein Design Laboratory, Graduate School of Medical Life Science, Yokohama City University, Suehiro 1-7-29, Tsurumi, Yokohama, 230-0045 Japan
| |
Collapse
|
9
|
Yu CC, Raj N, Chu JW. Statistical Learning of Protein Elastic Network from Positional Covariance Matrix. Comput Struct Biotechnol J 2023; 21:2524-2535. [PMID: 37095762 PMCID: PMC10121796 DOI: 10.1016/j.csbj.2023.03.033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2022] [Revised: 03/20/2023] [Accepted: 03/20/2023] [Indexed: 03/30/2023] Open
Abstract
Positional fluctuation and covariance during protein dynamics are key observables for understanding the molecular origin of biological functions. A frequently employed potential energy function for describing protein structural variation at the coarse-gained level is elastic network model (ENM). A long-standing issue in biomolecular simulation is thus the parametrization of ENM spring constants from the components of positional covariance matrix (PCM). Based on sensitivity analysis of PCM, the direct-coupling statistics of each spring, which is a specific combination of position fluctuation and covariance, is found to exhibit prominent signal of parameter dependence. This finding provides the basis for devising the objective function and the scheme of running through the effective one-dimensional optimization of every spring by self-consistent iteration. Formal derivation of the positional covariance statistical learning (PCSL) method also motivates the necessary data regularization for stable calculations. Robust convergence of PCSL is achieved in taking an all-atom molecular dynamics trajectory or an ensemble of homologous structures as input data. The PCSL framework can also be generalized with mixed objective functions to capture specific property such as the residue flexibility profile. Such physical chemistry-based statistical learning thus provides a useful platform for integrating the mechanical information encoded in various experimental or computational data.
Collapse
|
10
|
Rapp C, Nidetzky B. Hydride Transfer Mechanism of Enzymatic Sugar Nucleotide C2 Epimerization Probed with a Loose-Fit CDP-Glucose Substrate. ACS Catal 2022; 12:6816-6830. [PMID: 35747200 PMCID: PMC9207888 DOI: 10.1021/acscatal.2c00257] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2022] [Revised: 05/12/2022] [Indexed: 11/29/2022]
Abstract
![]()
Transient oxidation–reduction
through hydride transfer with
tightly bound NAD coenzyme is used by a large class of sugar nucleotide
epimerases to promote configurational inversion of carbon stereocenters
in carbohydrate substrates. A requirement for the epimerases to coordinate
hydride abstraction and re-addition with substrate rotation in the
binding pocket poses a challenge for dynamical protein conformational
selection linked to enzyme catalysis. Here, we studied the thermophilic
C2 epimerase from Thermodesulfatator atlanticus (TaCPa2E) in combination with a slow CDP-glucose
substrate (kcat ≈ 1.0 min–1; 60 °C) to explore the sensitivity of the enzymatic hydride
transfer toward environmental fluctuations affected by temperature
(20–80 °C). We determined noncompetitive primary kinetic
isotope effects (KIE) due to 2H at the glucose C2 and showed
that a normal KIE on the kcat (Dkcat) reflects isotope sensitivity of
the hydrogen abstraction to enzyme-NAD+ in a rate-limiting
transient oxidation. The Dkcat peaked at 40 °C was 6.1 and decreased to 2.1 at low (20 °C)
and 3.3 at high temperature (80 °C). The temperature profiles
for kcat with the 1H and 2H substrate showed a decrease in the rate below a dynamically
important breakpoint (∼40 °C), suggesting an equilibrium
shift to an impaired conformational landscape relevant for catalysis
in the low-temperature region. Full Marcus-like model fits of the
rate and KIE profiles provided evidence for a high-temperature reaction
via low-frequency conformational sampling associated with a broad
distribution of hydride donor–acceptor distances (long-distance
population centered at 3.31 ± 0.02 Å), only poorly suitable
for quantum mechanical tunneling. Collectively, dynamical characteristics
of TaCPa2E-catalyzed hydride transfer during transient
oxidation of CDP-glucose reveal important analogies to mechanistically
simpler enzymes such as alcohol dehydrogenase and dihydrofolate reductase.
A loose-fit substrate (in TaCPa2E) resembles structural
variants of these enzymes by extensive dynamical sampling to balance
conformational flexibility and catalytic efficiency.
Collapse
Affiliation(s)
- Christian Rapp
- Institute of Biotechnology and Biochemical Engineering, Graz University of Technology, Petersgasse 10-12/1, 8010 Graz, Austria
| | - Bernd Nidetzky
- Institute of Biotechnology and Biochemical Engineering, Graz University of Technology, Petersgasse 10-12/1, 8010 Graz, Austria
- Austrian Centre of Industrial Biotechnology (ACIB), Petersgasse 14, 8010 Graz, Austria
| |
Collapse
|
11
|
Varvdekar B, Prabhakant A, Krishnan M. Response of Terahertz Protein Vibrations to Ligand Binding: Calmodulin-Peptide Complexes as a Case Study. J Chem Inf Model 2022; 62:1669-1679. [PMID: 35312312 DOI: 10.1021/acs.jcim.1c01344] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Terahertz vibrations are sensitive reporters of the structure and interactions of proteins. Ligand binding alters the nature and distribution of these collective vibrations. The ligand-induced changes in the terahertz protein vibrations contribute to the binding entropy and to the overall thermodynamic stability of the resultant protein-ligand complexes. Here, we have examined the response of the low-frequency (below 6 terahertz) collective vibrations of the calcium-loaded calmodulin (CaM) to binding to five different ligands, both in the presence and absence of water, using normal-mode analysis and molecular dynamics simulations. A comparison of the vibrational spectra of hydrated and dry systems reveals that protein-solvent interactions stiffen the terahertz protein vibrations and that these solvent-coupled collective vibrations contribute significantly to the hydration-sensitive variation in the vibrational entropy of CaM. In the absence of water, the low-frequency vibrations of CaM are stiffened by ligand binding. On the contrary, the number and the cumulative vibrational entropy of low-frequency vibrational modes (ω < 200 cm-1) of the hydrated CaM are increased noticeably after binding to the peptides, indicating binding-induced softening of collective vibrations of the protein. Although the calculated and experimental binding affinities of the chosen complexes correlated reasonably well, no systematic correlation was observed between the protein vibrational entropy and the binding affinity. The results underscored the importance of the interplay of protein-ligand and solvent interactions in modulating the low-frequency vibrations of proteins.
Collapse
Affiliation(s)
- Bhagyesh Varvdekar
- Center for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology, Gachibowli, Hyderabad 500032, India
| | - Akshay Prabhakant
- Center for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology, Gachibowli, Hyderabad 500032, India
| | - Marimuthu Krishnan
- Center for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology, Gachibowli, Hyderabad 500032, India
| |
Collapse
|
12
|
Weng W, Weberg AB, Gera R, Tomson NC, Anna JM. Probing Ligand Effects on the Ultrafast Dynamics of Copper Complexes via Midinfrared Pump-Probe and 2DIR Spectroscopies. J Phys Chem B 2021; 125:12228-12241. [PMID: 34723540 DOI: 10.1021/acs.jpcb.1c06370] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The effects of ligand structural variation on the ultrafast dynamics of a series of copper coordination complexes were investigated using polarization-dependent mid-IR pump-probe spectroscopy and two-dimensional infrared (2DIR) spectroscopy. The series consists of three copper complexes [(R3P3tren)CuIIN3]BAr4F (1PR3, R3P3tren = tris[2-(phosphiniminato)ethyl]amine, BAr4F = tetrakis(pentafluorophenyl)borate) where the number of methyl and phenyl groups in the PR3 ligand are systematically varied across the series (PR3 = PMe3, PMe2Ph, PMePh2). The asymmetric stretching mode of azide in the 1PR3 series is used as a vibrational probe of the small-molecule binding site. The results of the pump-probe measurements indicate that the vibrational energy of azide dissipates through intramolecular pathways and that the bulkier phenyl groups lead to an increase in the spatial restriction of the diffusive reorientation of bound azide. From 2DIR experiments, we characterize the spectral diffusion of the azide group and find that an increase in the number of phenyl groups maps to a broader inhomogeneous frequency distribution (Δ2). This indicates that an increase in the steric bulk of the secondary coordination sphere acts to create more distinct configurations in the local environment that are accessible to the azide group. This work demonstrates how ligand structural variation affects the ultrafast dynamics of a small molecular group bound to the metal center, which could provide insight into the structure-function relationship of the copper coordination complexes and transition-metal coordination complexes in general.
Collapse
Affiliation(s)
- Wei Weng
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Alexander B Weberg
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Rahul Gera
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Neil C Tomson
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Jessica M Anna
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
13
|
Ojeda-May P, Mushtaq AU, Rogne P, Verma A, Ovchinnikov V, Grundström C, Dulko-Smith B, Sauer UH, Wolf-Watz M, Nam K. Dynamic Connection between Enzymatic Catalysis and Collective Protein Motions. Biochemistry 2021; 60:2246-2258. [PMID: 34250801 PMCID: PMC8297476 DOI: 10.1021/acs.biochem.1c00221] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
![]()
Enzymes employ a wide range of protein motions to achieve efficient catalysis of
chemical reactions. While the role of collective protein motions in substrate binding,
product release, and regulation of enzymatic activity is generally understood, their
roles in catalytic steps per se remain uncertain. Here, molecular dynamics simulations,
enzyme kinetics, X-ray crystallography, and nuclear magnetic resonance spectroscopy are
combined to elucidate the catalytic mechanism of adenylate kinase and to delineate the
roles of catalytic residues in catalysis and the conformational change in the enzyme.
This study reveals that the motions in the active site, which occur on a time scale of
picoseconds to nanoseconds, link the catalytic reaction to the slow conformational
dynamics of the enzyme by modulating the free energy landscapes of subdomain motions. In
particular, substantial conformational rearrangement occurs in the active site following
the catalytic reaction. This rearrangement not only affects the reaction barrier but
also promotes a more open conformation of the enzyme after the reaction, which then
results in an accelerated opening of the enzyme compared to that of the reactant state.
The results illustrate a linkage between enzymatic catalysis and collective protein
motions, whereby the disparate time scales between the two processes are bridged by a
cascade of intermediate-scale motion of catalytic residues modulating the free energy
landscapes of the catalytic and conformational change processes.
Collapse
Affiliation(s)
- Pedro Ojeda-May
- Department of Chemistry, Umeå University, Umeå SE-90187, Sweden.,High Performance Computing Centre North (HPC2N), Umeå University, Umeå SE-90187, Sweden
| | | | - Per Rogne
- Department of Chemistry, Umeå University, Umeå SE-90187, Sweden
| | - Apoorv Verma
- Department of Chemistry, Umeå University, Umeå SE-90187, Sweden
| | - Victor Ovchinnikov
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| | | | - Beata Dulko-Smith
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, Texas 76019, United States
| | - Uwe H Sauer
- Department of Chemistry, Umeå University, Umeå SE-90187, Sweden
| | | | - Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, Texas 76019, United States
| |
Collapse
|
14
|
Through bonds or contacts? Mapping protein vibrational energy transfer using non-canonical amino acids. Nat Commun 2021; 12:3284. [PMID: 34078890 PMCID: PMC8172543 DOI: 10.1038/s41467-021-23591-1] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2021] [Accepted: 05/06/2021] [Indexed: 11/08/2022] Open
Abstract
Vibrational energy transfer (VET) is essential for protein function. It is responsible for efficient energy dissipation in reaction sites, and has been linked to pathways of allosteric communication. While it is understood that VET occurs via backbone as well as via non-covalent contacts, little is known about the competition of these two transport channels, which determines the VET pathways. To tackle this problem, we equipped the β-hairpin fold of a tryptophan zipper with pairs of non-canonical amino acids, one serving as a VET injector and one as a VET sensor in a femtosecond pump probe experiment. Accompanying extensive non-equilibrium molecular dynamics simulations combined with a master equation analysis unravel the VET pathways. Our joint experimental/computational endeavor reveals the efficiency of backbone vs. contact transport, showing that even if cutting short backbone stretches of only 3 to 4 amino acids in a protein, hydrogen bonds are the dominant VET pathway.
Collapse
|
15
|
Schirò G, Fichou Y, Brogan APS, Sessions R, Lohstroh W, Zamponi M, Schneider GJ, Gallat FX, Paciaroni A, Tobias DJ, Perriman A, Weik M. Diffusivelike Motions in a Solvent-Free Protein-Polymer Hybrid. PHYSICAL REVIEW LETTERS 2021; 126:088102. [PMID: 33709739 DOI: 10.1103/physrevlett.126.088102] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/04/2020] [Accepted: 01/11/2021] [Indexed: 05/16/2023]
Abstract
The interaction between proteins and hydration water stabilizes protein structure and promotes functional dynamics, with water translational motions enabling protein flexibility. Engineered solvent-free protein-polymer hybrids have been shown to preserve protein structure, function, and dynamics. Here, we used neutron scattering, protein and polymer perdeuteration, and molecular dynamics simulations to explore how a polymer dynamically replaces water. Even though relaxation rates and vibrational properties are strongly modified in polymer coated compared to hydrated proteins, liquidlike polymer dynamics appear to plasticize the conjugated protein in a qualitatively similar way as do hydration-water translational motions.
Collapse
Affiliation(s)
- Giorgio Schirò
- Université Grenoble Alpes, CEA, CNRS, Institut de Biologie Structurale, F-38000 Grenoble, France
| | - Yann Fichou
- CNRS, Chemistry and Biology of Membranes and Nanoobjects (CBMN) UMR 5348, Institut Europeen de Chimie et Biologie (IECB), University of Bordeaux, 33600 Pessac, France
| | - Alex P S Brogan
- School of Biochemistry, Medical Sciences Building, University Walk, Bristol BS8 1TD, United Kingdom
| | - Richard Sessions
- Department of Chemistry, King's College London, Britannia House, London SE1 1DB, United Kingdom
| | - Wiebke Lohstroh
- Heinz Maier-Leibnitz Zentrum, Technische Universität München, Garching, Germany
| | - Michaela Zamponi
- Forschungszentrum Jülich GmbH, Jülich Centre for Neutron Science at Heinz Maier-Leibnitz Zentrum, Garching, Germany
| | - Gerald J Schneider
- Department of Chemistry and Department of Physics & Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - François-Xavier Gallat
- Université Grenoble Alpes, CEA, CNRS, Institut de Biologie Structurale, F-38000 Grenoble, France
| | - Alessandro Paciaroni
- Dipartimento di Fisica e Geologia, Università degli Studi di Perugia, 06123 Perugia, Italy
| | - Douglas J Tobias
- Department of Chemistry, University of California, Irvine, California, CA 92697, USA
| | - Adam Perriman
- School of Cellular and Molecular Medicine, University of Bristol, Bristol BS8 1TL, United Kingdom
| | - Martin Weik
- Université Grenoble Alpes, CEA, CNRS, Institut de Biologie Structurale, F-38000 Grenoble, France
| |
Collapse
|
16
|
Schafer JW, Schwartz SD. Directed Evolution's Influence on Rapid Density Fluctuations Illustrates How Protein Dynamics Can Become Coupled to Chemistry. ACS Catal 2020; 10:8476-8484. [PMID: 33163256 DOI: 10.1021/acscatal.0c01618] [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] [Indexed: 11/28/2022]
Abstract
Protein engineering is a growing field with a variety of experimental techniques available for altering protein function. However, creating an enzyme de novo is still in its infancy, so far yielding enzymes of modest catalytic efficiency. In this study, a system of artificial retro-aldolase enzymes found to have chemistry coupled to protein dynamics was examined. The original design was created computationally, and this protein was then subjected to directed evolution to improve the initial low catalytic efficiency. We found that this re-engineering of the enzyme resulted in rapid density fluctuations throughout the enzyme being reshaped via alterations in the hydrogen bonding network. This work also led to the discovery of a second important motion which aids in the release of an intermediate product. These results provide compelling evidence that to engineer efficient protein catalysts, fast protein dynamics need to be considered in the design.
Collapse
Affiliation(s)
- Joseph W. Schafer
- Department of Chemistry and Biochemistry, University of Arizona, Tucson, Arizona 85721, United States
| | - Steven D. Schwartz
- Department of Chemistry and Biochemistry, University of Arizona, Tucson, Arizona 85721, United States
| |
Collapse
|