1
|
Barton A, Sesin P, Diambra L. Simplifications and approximations in a single-gene circuit modeling. Sci Rep 2024; 14:12498. [PMID: 38822072 PMCID: PMC11143231 DOI: 10.1038/s41598-024-63265-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2024] [Accepted: 05/27/2024] [Indexed: 06/02/2024] Open
Abstract
The absence of detailed knowledge about regulatory interactions makes the use of phenomenological assumptions mandatory in cell biology modeling. Furthermore, the challenges associated with the analysis of these models compel the implementation of mathematical approximations. However, the constraints these methods introduce to biological interpretation are sometimes neglected. Consequently, understanding these restrictions is a very important task for systems biology modeling. In this article, we examine the impact of such simplifications, taking the case of a single-gene autoinhibitory circuit; however, our conclusions are not limited solely to this instance. We demonstrate that models grounded in the same biological assumptions but described at varying levels of detail can lead to different outcomes, that is, different and contradictory phenotypes or behaviors. Indeed, incorporating specific molecular processes like translation and elongation into the model can introduce instabilities and oscillations not seen when these processes are assumed to be instantaneous. Furthermore, incorporating a detailed description of promoter dynamics, usually described by a phenomenological regulatory function, can lead to instability, depending on the cooperative binding mechanism that is acting. Consequently, although the use of a regulating function facilitates model analysis, it may mask relevant aspects of the system's behavior. In particular, we observe that the two cooperative binding mechanisms, both compatible with the same sigmoidal function, can lead to different phenotypes, such as transcriptional oscillations with different oscillation frequencies.
Collapse
Affiliation(s)
- Alejandro Barton
- Centro Regional de Estudios Genómicos, Universidad Nacional de La Plata, La Plata, Argentina
- Consejo Nacional de Investigaciones Científicas y Técnicas, Buenos Aires, Argentina
| | - Pablo Sesin
- Departamento de Física Teórica, GAIDI, Comisión Nacional de Energía Atómica, 1429, Buenos Aires, Argentina
| | - Luis Diambra
- Centro Regional de Estudios Genómicos, Universidad Nacional de La Plata, La Plata, Argentina.
- Consejo Nacional de Investigaciones Científicas y Técnicas, Buenos Aires, Argentina.
| |
Collapse
|
2
|
Cooperativity, absolute interaction, and algebraic optimization. J Math Biol 2020; 81:1169-1191. [PMID: 32968839 PMCID: PMC7557646 DOI: 10.1007/s00285-020-01540-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2019] [Revised: 05/04/2020] [Indexed: 11/03/2022]
Abstract
We consider a measure of cooperativity based on the minimal interaction required to generate an observed titration behavior. We describe the corresponding algebraic optimization problem and show how it can be solved using the nonlinear algebra tool SCIP. Moreover, we compute the minimal interactions and minimal molecules for several binding polynomials that describe the oxygen binding of various hemoglobins under different conditions. We compare their minimal interaction with the maximal slope of the Hill plot, and discuss similarities and discrepancies with a view towards the shapes of the binding curves.
Collapse
|
3
|
Sevlever F, Di Bella JP, Ventura AC. Discriminating between negative cooperativity and ligand binding to independent sites using pre-equilibrium properties of binding curves. PLoS Comput Biol 2020; 16:e1007929. [PMID: 32497065 PMCID: PMC7297384 DOI: 10.1371/journal.pcbi.1007929] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2019] [Revised: 06/16/2020] [Accepted: 05/06/2020] [Indexed: 12/02/2022] Open
Abstract
Negative cooperativity is a phenomenon in which the binding of a first ligand or substrate molecule decreases the rate of subsequent binding. This definition is not exclusive to ligand-receptor binding, it holds whenever two or more molecules undergo two successive binding events. Negative cooperativity turns the binding curve more graded and cannot be distinguished from two independent and different binding events based on equilibrium measurements only. The need of kinetic data for this purpose was already reported. Here, we study the binding response as a function of the amount of ligand, at different times, from very early times since ligand is added and until equilibrium is reached. Over those binding curves measured at different times, we compute the dynamic range: the fold change required in input to elicit a change from 10 to 90% of maximum output, finding that it evolves in time differently and controlled by different parameters in the two situations that are identical in equilibrium. Deciphering which is the microscopic model that leads to a given binding curve adds understanding on the molecular mechanisms at play, and thus, is a valuable tool. The methods developed in this article were tested both with simulated and experimental data, showing to be robust to noise and experimental constraints.
Collapse
Affiliation(s)
- Federico Sevlever
- Department of Physiology, Molecular and Cellular Biology, University of Buenos Aires, Buenos Aires, Argentina
- Institute of Physiology, Molecular Biology and Neurosciences, National Research Council (CONICET), Buenos Aires, Argentina
| | - Juan Pablo Di Bella
- Department of Physiology, Molecular and Cellular Biology, University of Buenos Aires, Buenos Aires, Argentina
- Institute of Physiology, Molecular Biology and Neurosciences, National Research Council (CONICET), Buenos Aires, Argentina
| | - Alejandra C. Ventura
- Department of Physiology, Molecular and Cellular Biology, University of Buenos Aires, Buenos Aires, Argentina
- Institute of Physiology, Molecular Biology and Neurosciences, National Research Council (CONICET), Buenos Aires, Argentina
| |
Collapse
|
4
|
Rong C, Zhao D, Zhou T, Liu S, Yu D, Liu S. Homogeneous Molecular Systems are Positively Cooperative, but Charged Molecular Systems are Negatively Cooperative. J Phys Chem Lett 2019; 10:1716-1721. [PMID: 30916572 DOI: 10.1021/acs.jpclett.9b00639] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Molecular systems bound together through noncovalent interactions are ubiquitous in nature, many of which are involved in essential life processes, yet little is known about the principles governing their structure, stability, and function. Cooperativity as one of the intrinsic properties in these systems plays a key role. In this work, on the basis of our recent quantification scheme of the cooperativity effect, we present a general pattern to identify which systems are positively cooperative and which are negatively cooperative. We show that cooperativity in homogeneous molecular systems is positive, but cooperativity in charged molecular systems is negative. We also employ analytical tools from energetics and information perspectives to appreciate the origin of the cooperativity effect. We find that positive cooperativity is dominated by the exchange-correlation interaction and steric effect, whereas negative cooperativity is governed by the electrostatic interaction. Our results should have strong implications for better understanding molecular recognition, protein folding, signal transduction, allosteric regulation, and other processes.
Collapse
Affiliation(s)
- Chunying Rong
- College of Chemistry and Chemical Engineering , Hunan Normal University , Changsha Hunan 410081 , People's Republic of China
| | - Dongbo Zhao
- School of Chemistry and Chemical Engineering , Nanjing University , Nanjing 210093 , People's Republic of China
| | - Tianjing Zhou
- College of Chemistry and Chemical Engineering , Hunan Normal University , Changsha Hunan 410081 , People's Republic of China
| | - Siyuan Liu
- College of Chemistry and Chemical Engineering , Hunan Normal University , Changsha Hunan 410081 , People's Republic of China
| | - Donghai Yu
- College of Chemistry and Chemical Engineering , Hunan Normal University , Changsha Hunan 410081 , People's Republic of China
| | - Shubin Liu
- Research Computing Center , University of North Carolina , Chapel Hill , North Carolina 27599-3420 , United States
| |
Collapse
|
5
|
Zhou T, Liu S, Yu D, Zhao D, Rong C, Liu S. On the negative cooperativity of argon clusters containing one lithium cation or fluorine anion. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2018.12.019] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
|
6
|
Ren Y, Martini JWR, Torres J. Decoupled molecules with binding polynomials of bidegree (n, 2). J Math Biol 2018; 78:879-898. [PMID: 30283984 PMCID: PMC6510845 DOI: 10.1007/s00285-018-1295-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2017] [Revised: 09/09/2018] [Indexed: 01/15/2023]
Abstract
We present a result on the number of decoupled molecules for systems binding two different types of ligands. In the case of n and 2 binding sites respectively, we show that there are \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$2(n!)^{2}$$\end{document}2(n!)2 decoupled molecules to a generic binding polynomial. For molecules with more binding sites for the second ligand, we provide computational results.
Collapse
Affiliation(s)
- Yue Ren
- Max-Planck-Institut für Mathematik in den Naturwissenschaften, Inselstrasse 22, 04103, Leipzig, Germany.
| | | | - Jacinta Torres
- Max-Planck-Institut für Mathematik in den Naturwissenschaften, Inselstrasse 22, 04103, Leipzig, Germany
| |
Collapse
|
7
|
Rong C, Zhao D, Yu D, Liu S. Quantification and origin of cooperativity: insights from density functional reactivity theory. Phys Chem Chem Phys 2018; 20:17990-17998. [PMID: 29927447 DOI: 10.1039/c8cp03092h] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
Cooperativity is a widely used chemical concept whose existence is ubiquitous in chemical and biological systems but whose quantification is still controversial and origin much less appreciated. In this work, using the interaction energy of a molecular system, which is composed of multiple copies of a building block, we propose a quantitative measurement to evaluate the cooperativity effect. This quantification approach is then applied to six molecular systems, i.e., water cluster, argon cluster, protonated water cluster, zinc atom cluster, water cluster on top of a graphene sheet, and alpha helix of glycine amino acids, each with up to 20 copies of the building block. Cooperativity is seen in all these systems. Both positive and negative cooperativity effects are observed. Employing the two energy partition schemes in density functional theory and the information-theoretic quantities such as Shannon entropy, Fisher information, information gain, etc., we then examine the origin of the cooperativity effect for these systems. Strong linear correlations between the cooperativity measure and some of these theoretical quantities have been unveiled. With these correlations, we are able to quantitatively account for their origin of cooperativity. Our results show that the interactions governing the existence and validity of the cooperativity effect are complicated. An opposite mechanism in enthalpy-entropy compensation for positive and negative cooperativity has been unveiled. These results should provide new insights and understandings from a different viewpoint about the nature and origin of cooperativity to appreciate this vastly important chemical concept.
Collapse
Affiliation(s)
- Chunying Rong
- College of Chemistry and Chemical Engineering, Hunan Normal University, Changsha Hunan 410081, P. R. China
| | | | | | | |
Collapse
|
8
|
Martini JWR. A measure to quantify the degree of cooperativity in overall titration curves. J Theor Biol 2017; 432:33-37. [PMID: 28803910 DOI: 10.1016/j.jtbi.2017.08.010] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2017] [Revised: 07/31/2017] [Accepted: 08/09/2017] [Indexed: 11/17/2022]
Abstract
In the framework of the grand canonical ensemble, different definitions of cooperativity commonly used in the context of ligand binding are not equivalent. A unifying definition is the existence of non-real roots of the binding polynomial. Using this qualitative criterion, an open question is how to quantify the degree of cooperativity. In this work, we introduce a theoretical measure to quantify the degree of cooperativity of a titriation curve. Its definition is based on a minimal energy approach mapping a given binding polynomial to the minimal interaction energy which is required to generate it. We show that the degree of cooperativity can be calculated easily, if the molecule under consideration is assumed to consist of energetically identical binding sites. Moreover, the property of sub-multiplicativity allows us to determine upper bounds for the degree of cooperativity in asymmetric systems. The approach is consistent with the qualitative definition of cooperativity based on the existence of non-real roots of the binding polynomial, and thus helps to put the concept of cooperativity on a solid theoretical ground. It connects macro- and microstates, but takes here also into account that an infinite number of different molecules can cause the same macroscopic ligand binding behavior, which means that the underlying microsystem cannot be uniquely identified based on the titration curve only.
Collapse
|
9
|
Liu C, Liu Z, Wang J. Uncovering the molecular and physiological processes of anticancer leads binding human serum albumin: A physical insight into drug efficacy. PLoS One 2017; 12:e0176208. [PMID: 28426740 PMCID: PMC5398698 DOI: 10.1371/journal.pone.0176208] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2017] [Accepted: 04/06/2017] [Indexed: 12/24/2022] Open
Abstract
Human serum albumin (HSA) has its ability to bind drug molecules and influence their efficacies. Although anticancer leads NSC48693 and NSC290956 functioned at the same mechanism, the drug efficacies were obviously distinct. To gain insight into the distinct drug efficacy, the molecular and physiological processes of anticancer leads binding HSA have been investigated via a combined experimental and theoretical approach. The binding site, as characterized by fluorescence quenching and molecular modeling, is found to be located at site II in subdomain III A for NSC48693 with tight binding and at site FA1 in subdomain I B for NSC290956 with negatively cooperative binding, respectively. As indicated by the thermodynamic analysis, NSC48693 binds to HSA with an enthalpy driven mechanism, while NSC290956 binding with HSA is entropically driven. The further kinetic analysis indicates that the association rates appear to be similar to these two anticancer leads, however, the dissociation rate of NSC48693 is approximately 5-fold slower than that of NSC290956. For NSC48693, the pharmacodynamic efficacy is less than that of NSC290956, while its pharmacokinetic behavior is better than that of NSC290956. These parameters influence the pharmacodynamic efficacy and pharmacokinetic behavior, which will give further impacts on drug efficacy in vivo.
Collapse
Affiliation(s)
- Chuanbo Liu
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin, P.R. China
- University of Chinese Academy of Sciences, Beijing, P.R. China
| | - Zuojia Liu
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin, P.R. China
- * E-mail: (ZL); (JW)
| | - Jin Wang
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin, P.R. China
- Department of Chemistry and Physics, State University of New York, Stony Brook, New York, United States of America
- * E-mail: (ZL); (JW)
| |
Collapse
|
10
|
Abstract
Cooperativity as a concept is easy to grasp intuitively, but surprisingly hard to define. Two recent papers shed light on the issue and continue the debate on how best to define cooperative binding.
Collapse
Affiliation(s)
- Melanie I Stefan
- Centre for Integrative Physiology, Edinburgh Medical School: Biomedical Sciences, University of Edinburgh, Edinburgh, UK.
| |
Collapse
|
11
|
Martini JWR, Schlather M, Schütz S. A Model for Carrier-Mediated Biological Signal Transduction Based on Equilibrium Ligand Binding Theory. Bull Math Biol 2016; 78:1039-57. [PMID: 27230608 DOI: 10.1007/s11538-016-0173-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2015] [Accepted: 05/03/2016] [Indexed: 10/21/2022]
Abstract
Different variants of a mathematical model for carrier-mediated signal transduction are introduced with focus on the odor dose-electrophysiological response curve of insect olfaction. The latter offers a unique opportunity to observe experimentally the effect of an alteration in the carrier molecule composition on the signal molecule-dependent response curve. Our work highlights the role of involved carrier molecules, which have largely been ignored in mathematical models for response curves in the past. The resulting model explains how the involvement of more than one carrier molecule in signal molecule transport can cause dose-response curves as observed in experiments, without the need of more than one receptor per neuron. In particular, the model has the following features: (1) An extended sensitivity range of neuronal response is implemented by a system consisting of only one receptor but several carrier molecules with different affinities for the signal molecule. (2) Given that the sensitivity range is extended by the involvement of different carrier molecules, the model implies that a strong difference in the expression levels of the carrier molecules is absolutely essential for wide range responses. (3) Complex changes in dose-response curves which can be observed when the expression levels of carrier molecules are altered experimentally can be explained by interactions between different carrier molecules. The principles we demonstrate here for electrophysiological responses can also be applied to any other carrier-mediated biological signal transduction process. The presented concept provides a framework for modeling and statistical analysis of signal transduction processes if sufficient information on the underlying biology is available.
Collapse
Affiliation(s)
- Johannes W R Martini
- Animal Breeding and Genetics Group, Department of Animal Sciences, Georg-August University, Göttingen, Germany.
| | - Martin Schlather
- Fakultät für Wirtschaftsinformatik und Wirtschaftsmathematik, Universität Mannhein, Mannheim, Germany
| | - Stefan Schütz
- Büsgen-Institut, Abteilung Forstzoologie und Waldschutz, Georg-August-Universität, Göttingen, Germany
| |
Collapse
|
12
|
On Hill coefficients and subunit interaction energies. J Math Biol 2016; 73:1399-1411. [PMID: 27038164 DOI: 10.1007/s00285-016-1001-9] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2015] [Revised: 02/29/2016] [Indexed: 10/22/2022]
Abstract
The study of cooperative ligand binding to multimeric proteins aims to explain complex cooperative binding phenomena using concepts derived from ideal binding isotherms. The purpose of such efforts is the dissection of the cooperative binding isotherm into its interacting components, a result with a clear mechanistic value. Historically, cooperative binding is usually quantified using the Hill coefficient, [Formula: see text], defined as the slope of the Hill plot at 50 % saturation. It was previously shown that the slope of the Hill plot throughout the titration is equal to the ratio of the binding variance in the system under study, to the binding variance of a reference non-interacting system. In the present contribution, this leads to a broader approach towards quantifying cooperativity, which empirically links cooperativity to the ensemble average of the subunit interaction energy. The resulting equations can be used to derive average differential subunit interaction energies directly from experimental binding isotherms. Combined with recent experimental advances in assessing binding distributions in multimeric proteins, these equations can also be used to calculate individual subunit interaction energies for specific n-ligated protein species.
Collapse
|