1
|
Juretić D, Bonačić Lošić Ž. Theoretical Improvements in Enzyme Efficiency Associated with Noisy Rate Constants and Increased Dissipation. ENTROPY (BASEL, SWITZERLAND) 2024; 26:151. [PMID: 38392406 PMCID: PMC10888251 DOI: 10.3390/e26020151] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/30/2023] [Revised: 01/18/2024] [Accepted: 02/05/2024] [Indexed: 02/24/2024]
Abstract
Previous studies have revealed the extraordinarily large catalytic efficiency of some enzymes. High catalytic proficiency is an essential accomplishment of biological evolution. Natural selection led to the increased turnover number, kcat, and enzyme efficiency, kcat/KM, of uni-uni enzymes, which convert a single substrate into a single product. We added or multiplied random noise with chosen rate constants to explore the correlation between dissipation and catalytic efficiency for ten enzymes: beta-galactosidase, glucose isomerase, β-lactamases from three bacterial strains, ketosteroid isomerase, triosephosphate isomerase, and carbonic anhydrase I, II, and T200H. Our results highlight the role of biological evolution in accelerating thermodynamic evolution. The catalytic performance of these enzymes is proportional to overall entropy production-the main parameter from irreversible thermodynamics. That parameter is also proportional to the evolutionary distance of β-lactamases PC1, RTEM, and Lac-1 when natural or artificial evolution produces the optimal or maximal possible catalytic efficiency. De novo enzyme design and attempts to speed up the rate-limiting catalytic steps may profit from the described connection between kinetics and thermodynamics.
Collapse
Affiliation(s)
- Davor Juretić
- Mediterranean Institute for Life Sciences, Šetalište Ivana Meštrovića 45, 21000 Split, Croatia
- Faculty of Science, University of Split, Ruđera Boškovića 33, 21000 Split, Croatia
| | | |
Collapse
|
2
|
Kim S, Hooker AC, Shi Y, Kim GHJ, Wong WK. Metaheuristics for pharmacometrics. CPT Pharmacometrics Syst Pharmacol 2021; 10:1297-1309. [PMID: 34562342 PMCID: PMC8592519 DOI: 10.1002/psp4.12714] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2020] [Revised: 08/06/2021] [Accepted: 09/07/2021] [Indexed: 12/22/2022] Open
Abstract
Metaheuristics is a powerful optimization tool that is increasingly used across disciplines to tackle general purpose optimization problems. Nature-inspired metaheuristic algorithms is a subclass of metaheuristic algorithms and have been shown to be particularly flexible and useful in solving complicated optimization problems in computer science and engineering. A common practice with metaheuristics is to hybridize it with another suitably chosen algorithm for enhanced performance. This paper reviews metaheuristic algorithms and demonstrates some of its utility in tackling pharmacometric problems. Specifically, we provide three applications using one of its most celebrated members, particle swarm optimization (PSO), and show that PSO can effectively estimate parameters in complicated nonlinear mixed-effects models and to gain insights into statistical identifiability issues in a complex compartment model. In the third application, we demonstrate how to hybridize PSO with sparse grid, which is an often-used technique to evaluate high dimensional integrals, to search for D -efficient designs for estimating parameters in nonlinear mixed-effects models with a count outcome. We also show the proposed hybrid algorithm outperforms its competitors when sparse grid is replaced by its competitor, adaptive gaussian quadrature to approximate the integral, or when PSO is replaced by three notable nature-inspired metaheuristic algorithms.
Collapse
Affiliation(s)
- Seongho Kim
- Department of OncologyWayne State UniversityDetroitMichiganUSA
| | | | - Yu Shi
- Department of BiostatisticsUniversity of California Los AngelesLos AngelesCaliforniaUSA
| | - Grace Hyun J. Kim
- Department of BiostatisticsUniversity of California Los AngelesLos AngelesCaliforniaUSA
| | - Weng Kee Wong
- Department of BiostatisticsUniversity of California Los AngelesLos AngelesCaliforniaUSA
| |
Collapse
|
3
|
Juretić D, Bonačić Lošić Ž. Comments on 'Flexibility of enzymatic transitions as a hallmark of optimized enzyme steady-state kinetics and thermodynamics'. Comput Biol Chem 2021; 95:107571. [PMID: 34531144 DOI: 10.1016/j.compbiolchem.2021.107571] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Davor Juretić
- Mediterranean Institute for Life Sciences, 21000 Split, Croatia; University of Split, Faculty of Science, 21000 Split, Croatia.
| | | |
Collapse
|
4
|
Gerlt JA. Evolution of Enzyme Function and the Development of Catalytic Efficiency: Triosephosphate Isomerase, Jeremy R. Knowles, and W. John Albery. Biochemistry 2021; 60:3529-3538. [PMID: 34015914 DOI: 10.1021/acs.biochem.1c00211] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Every reader knows that an enzyme accelerates a reaction by reducing the activation-energy barrier. However, understanding how this is achieved by the structure of the enzyme and its interactions with stable complexes and transition states and, then, using this to (re)design enzymes to catalyze novel reactions remain the "holy grail" of mechanistic enzymology. The necessary foundation is the free-energy profile that specifies the energies of the bound substate, product, and intervening intermediates as well as the transition states by which they are interconverted. When this free-energy profile is compared to that for the uncatalyzed reaction, strategies for establishing and enhancing catalysis can be identified. This Perspective reminds readers that the first free-energy profile determined for an enzyme-catalyzed reaction, that for triosephosphate isomerase, was published in Biochemistry in 1976 by Jeremy R. Knowles, W. John Albery, and co-workers. They used the profile to propose three steps of increasing "subtlety" that can be influenced by evolutionary pressure to increase the flux through the reaction coordinate: (1) "uniform binding" of the substrate, product, and intermediates; (2) "differential binding" of complexes so that these are isoenergetic (to minimize the energy of the intervening transition states); and (3) "catalysis of an elementary step" in which the transition state for the kinetically significant chemical step is stabilized so that flux can be determined by the rate of substrate binding or product dissociation. These papers continue to guide mechanistic studies of enzyme-catalyzed reactions and provide principles for the (re)design of novel enzymes.
Collapse
Affiliation(s)
- John A Gerlt
- Departments of Biochemistry and Chemistry and Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
5
|
Enzyme free energy profiles: Can substrate binding be nonspontaneous? Can ground state interactions enhance catalysis? Biophys Chem 2021; 274:106606. [PMID: 33945990 DOI: 10.1016/j.bpc.2021.106606] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2021] [Revised: 04/23/2021] [Accepted: 04/25/2021] [Indexed: 11/21/2022]
Abstract
Two influential enzymological theories were proposed in the late 1970s - that catalytic power stems only from transition state stabilization, while ground state interactions are either irrelevant or inhibitory; and enzyme substrate binding is nonspontaneous at low substrate concentrations ([S]0 << Km). I show here that ground state destabilization can be a very effective source of catalytic power, especially at high substrate concentrations, and enzyme-substrate binding thermodynamics are independent of initial substrate concentration. Binding free energy ranges from negative (spontaneous) under pre-steady state conditions up to a maximum of zero at steady state. Nonspontaneous binding can only occur under standard state conditions when c° is defined to be less than Km.
Collapse
|
6
|
Bearne SL. The role of Brønsted base basicity in estimating carbon acidity at enzyme active sites: a caveat. Org Biomol Chem 2019; 17:7161-7165. [PMID: 31317156 DOI: 10.1039/c9ob00863b] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Many enzymes catalyze the abstraction of a proton from a carbon acid substrate to initiate a variety of reactions; however, the development of a complete quantitative description of enzyme-catalyzed heterolytic cleavage of a C-H bond remains a challenge to enzymologists. To determine the pK value for such substrates bound at the active site, recent studies have estimated the equilibrium for formation of the deprotonated intermediate at the active site, however, accurate knowledge of the pK of the conjugate acid of the Brønsted base catalyst (BH+) is also required. Herein, it is shown that using the value of pK of the enzyme-substrate complex can underestimate the value of pK by an amount between zero and pδ, where pδ is the change in basicity of BH+ upon going from the enzyme-substrate complex to the enzyme-intermediate complex.
Collapse
Affiliation(s)
- Stephen L Bearne
- Department of Biochemistry & Molecular Biology, Dalhousie University, Halifax, NS B3H 4R2, Canada. and Department of Chemistry, Dalhousie University, Halifax, NS B3H 4R2, Canada
| |
Collapse
|
7
|
Abstract
The pKa values for substrates acting as carbon acids (i.e., C-H deprotonation reactions) in several enzyme active sites are presented. The information needed to calculate them includes the pKa of the active site acid/base catalyst and the equilibrium constant for the deprotonation step. Carbon acidity is obtained from the relation pKeq = pKar–pKap = ΔpKa for a proton transfer reaction. Five enzymatic free energy profiles (FEPs) were calculated to obtain the equilibrium constants for proton transfer from carbon in the active site, and six additional proton transfer equilibrium constants were extracted from data available in the literature, allowing substrate C-H pKas to be calculated for 11 enzymes. Active site-bound substrate C-H pKa values range from 5.6 for ketosteroid isomerase to 16 for proline racemase. Compared to values in water, enzymes lower substrate C-H pKas by up to 23 units, corresponding to 31 kcal/mol of carbanion stabilization energy. Calculation of Marcus intrinsic barriers (ΔG0‡) for pairs of non-enzymatic/enzymatic reactions shows significant reductions in ΔG0‡ for cofactor-independent enzymes, while pyridoxal phosphate dependent enzymes appear to increase ΔG0‡ to a small extent as a consequence of carbanion resonance stabilization. The large increases in carbon acidity found here are central to the large rate enhancements observed in enzymes that catalyze carbon deprotonation.
Collapse
Affiliation(s)
- Michael D Toney
- Department of Chemistry, University of California, Davis, Davis, CA, United States
| |
Collapse
|
8
|
Extracting Kinetic Isotope Effects From a Global Analysis of Reaction Progress Curves. Methods Enzymol 2017. [PMID: 28911785 DOI: 10.1016/bs.mie.2017.06.041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
Abstract
Enzyme reaction progress curves, or time course datasets, are often rich in information, yet their analysis typically reduces their information content to a single parameter, the initial velocity. An alternative approach is described here, where the time course is described by a model constructed from rate equations. In combination with global nonlinear regression, intrinsic rate and/or equilibrium constants can be directly obtained by fitting these data. This method can be greatly enhanced when combined with the measurement of (usually deuterium) isotope effects, which selectively perturb individual step(s) within the reaction, allowing better separation of fitted parameters and more robust model testing. This chapter focuses on practical considerations when using analytical and/or numerically integrated rate equations to model enzyme reactions. The emphasis is on the underlying methodology, which is demonstrated with specific examples alongside recommendations of suitable software.
Collapse
|
9
|
Bonačić Lošić Ž, Donđivić T, Juretić D. Is the catalytic activity of triosephosphate isomerase fully optimized? An investigation based on maximization of entropy production. J Biol Phys 2017; 43:69-86. [PMID: 28050739 PMCID: PMC5323346 DOI: 10.1007/s10867-016-9434-3] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2016] [Accepted: 11/21/2016] [Indexed: 11/29/2022] Open
Abstract
Triosephosphate isomerase (TIM) is often described as a fully evolved housekeeping enzyme with near-maximal possible reaction rate. The assumption that an enzyme is perfectly evolved has not been easy to confirm or refute. In this paper, we use maximization of entropy production within known constraints to examine this assumption by calculating steady-state cyclic flux, corresponding entropy production, and catalytic activity in a reversible four-state scheme of TIM functional states. The maximal entropy production (MaxEP) requirement for any of the first three transitions between TIM functional states leads to decreased total entropy production. Only the MaxEP requirement for the product (R-glyceraldehyde-3-phosphate) release step led to a 30% increase in enzyme activity, specificity constant kcat/KM, and overall entropy production. The product release step, due to the TIM molecular machine working in the physiological direction of glycolysis, has not been identified before as the rate-limiting step by using irreversible thermodynamics. Together with structural studies, our results open the possibility for finding amino acid substitutions leading to an increased frequency of loop six opening and product release.
Collapse
Affiliation(s)
| | - Tomislav Donđivić
- Medical High School, Šibenik, Ante Šupuka bb, 22000, Šibenik, Croatia
| | - Davor Juretić
- Mediterranean Institute for Life Sciences, Šetalište Ivana Meštrovića 45, 21000, Split, Croatia.
| |
Collapse
|
10
|
Effects of Glycosylation on the Enzymatic Activity and Mechanisms of Proteases. Int J Mol Sci 2016; 17:ijms17121969. [PMID: 27898009 PMCID: PMC5187769 DOI: 10.3390/ijms17121969] [Citation(s) in RCA: 73] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2016] [Revised: 11/07/2016] [Accepted: 11/10/2016] [Indexed: 02/06/2023] Open
Abstract
Posttranslational modifications are an important feature of most proteases in higher organisms, such as the conversion of inactive zymogens into active proteases. To date, little information is available on the role of glycosylation and functional implications for secreted proteases. Besides a stabilizing effect and protection against proteolysis, several proteases show a significant influence of glycosylation on the catalytic activity. Glycans can alter the substrate recognition, the specificity and binding affinity, as well as the turnover rates. However, there is currently no known general pattern, since glycosylation can have both stimulating and inhibiting effects on activity. Thus, a comparative analysis of individual cases with sufficient enzyme kinetic and structural data is a first approach to describe mechanistic principles that govern the effects of glycosylation on the function of proteases. The understanding of glycan functions becomes highly significant in proteomic and glycomic studies, which demonstrated that cancer-associated proteases, such as kallikrein-related peptidase 3, exhibit strongly altered glycosylation patterns in pathological cases. Such findings can contribute to a variety of future biomedical applications.
Collapse
|
11
|
Mak WS, Siegel JB. Computational enzyme design: Transitioning from catalytic proteins to enzymes. Curr Opin Struct Biol 2014; 27:87-94. [DOI: 10.1016/j.sbi.2014.05.010] [Citation(s) in RCA: 56] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2014] [Revised: 05/23/2014] [Accepted: 05/30/2014] [Indexed: 12/14/2022]
|