Original Research ARTICLE
Computational Model for Membrane Transporters. Potential Implications for Cancer
- 1Instituto de Ciencias, Universidad Nacional de General Sarmiento, Buenos Aires, Argentina
- 2National Scientific and Technical Research Council, Buenos Aires, Argentina
- 3Departament de Física de la Matèria Condensada, Universitat de Barcelona, Barcelona, Spain
To explain the increased transport of nutrients and metabolites and to control the movement of drug molecules through the transporters to the cancer cells, it is important to understand the exact mechanism of their structure and activity, as well as their biological and physical characteristics. We propose a computational model that reproduces the functionality of membrane transporters by quantifying the flow of substrates through the cell membrane. The model identifies the force induced by conformational changes of the transporter due to hydrolysis of ATP, in ABC transporters, or by an electrochemical gradient of ions, in secondary transporters. The transport rate is computed by averaging the velocity generated by the force along the paths followed by the substrates. The results obtained are in accordance with the experiments. The model provides an overall framework for analyzing the membrane transport proteins that regulate the flows of ions, nutrients and other molecules across the cell membranes, and their activities.
Cancer cells synthesize increased amount of fatty acids and protein building blocks to successfully divide and metastasize. To support their metabolism cancer cells require increased supply of metabolic substrates and nutrients. Plasma membrane transporters secure import of a wide range of substrates into the cytoplasm. Consistently, increased expression of several transporting proteins has been correlated to the increased metabolic activity of cancer cells and poor disease prognosis (Natecz, 2020; Sampedro-Núñez et al., 2020; Xu et al., 2020; Yamada et al., 2020).
There are two types of transporters, ATP biding cassette (ABC) and secondary active transporters. Secondary active transporters carry substrates across the plasma membrane using electrochemical gradient of ions. They bind the substrate on one side of the membrane followed by a conformational change allowing to release the substrate on the other side of the membrane (Boudker and Verdon, 2010). Several of these transporters have been shown to play an important role in cancer cells by increasing uptake of such substrates as glucose, glutamine, lactate, etc., and supporting cancer cell metabolism (Payen et al., 2017; Reckzeh et al., 2019; Scalise et al., 2020).
ABC transporters belong to a large family of transporters carrying several different substrates across the plasma membrane. To translocate substrates ABC transporters use the energy from ATP hydrolysis. They play an important role in keeping cellular homeostasis by regulating the level of several molecules including peptides, lipids, amino acids, or drugs (El-Awady et al., 2017; Neumann et al., 2017). Importantly, ABC transporters have also been shown to contribute to anti-cancer drug resistance presenting a major problem in finding an efficient anti-cancer therapy (Robey et al., 2018; Bock et al., 2019; Asif et al., 2020). These transporters carry out anti-cancer drugs efflux from cancer cells significantly decreasing the efficiency of anti-cancer therapies.
To be able to successfully overcome the problem of the increased transport of nutrients and metabolites and control drugs influx/efflux through transporters in cancer cells it is important to understand the exact mechanism of their structure and activity as well as their biological and physical features. This has been the objective of this article.
We have proposed a model that describes the trans-membrane transport and calculates the flow of substrates through the membrane. The model assumes that changes in transporter conformations induce forces that contribute to the translocation of substrates. Since entropy is a measure of conformations, the forces are referred to as entropic forces. These forces thus encode the information of the change in the available space of the substrates due to conformation changes (Zwanzig, 1992; Reguera and Rubí, 2001; Vazquez et al., 2008; Carusela and Rubi, 2017, 2018; Rubi, 2019). The computed translocation rates scale with the ratio between the time it takes for conformations to change and the time in which substrates diffuse through the transporter. This scaling behavior makes our model a general tool in the study of membrane transport.
2. The Model
We model a transporter as a small funnel-shaped motor whose structure changes over time by alternately opening and closing to the intra/extracellular medium, as represented in Figure 1. This form has been observed for example in P-glycoprotein whose structure is narrow at the cytoplasmic side, of about 9–25 Å in the middle, and wider at the extracellular surface (Loo and Clarke, 2001). The conical form has also been observed in pumps (Rubi et al., 2017) by means of crystallization experiments (Olesen et al., 2007).
Figure 1. Conformational changes of a transporter modeled by an oscillating conical-shaped channel which closes/opens to the intra/extra cellular medium. stands for the cross-sectional area at the narrowest part whereas denotes that area at positions inside the transporter and time. Figures at the bottom represent the initial and final states and an intermediate state of the cone oscillations in the model.
This periodic movement makes it possible for substrates, such as aminoacids, ions, neurotransmitters, nutrients and different drugs, to overcome the potential barriers generated by interactions allowing them to pass to the other side of the plasma membrane.
Alternating gating increases transport efficiency with respect to that of diffusion. Transport of substrates through the membrane is the result of changes in its conformation. Primary active transporters couple substrate movements to a source of chemical energy, such as ATP hydrolysis. Secondary active transporters are driven by electrochemical gradients of ions. Transporters differ from ion channels in that their turnover rate is much slower than that of channels which is typically of the order of 106s−1. The rate of ABC transporters such as LeuT, MsbA, and of MFS secondary transporters frequently falls within the range (10−1 − 103)s−1 (Ashcroft et al., 2009; Liu et al., 2018; Fitzgerald et al., 2019).
2.1. Force Induced by Transporter Conformation Changes
Our model considers that the changes in the conformation of both types of transporters that allow the passage of the substrates entail a variation in the space they have to move. This fact affects the entropy of the substrates as this quantity measures the degree of disorder of a system which in our case is less in the narrow area of the transporter, where the substrates have fewer positions to occupy, and more in the wider area where the space available is greater. This difference of entropies between the narrow and the wide part of the transporter gives rise to a gradient of free energy and consequently to a force on the substrates that we will call entropic force (see Figure 2).
Figure 2. Forces acting on the substrates. Changes in the shape of the transporter induce an entropic force Fent, larger than the force Fμ generated by a difference of the concentrations c1 and c2, which helps to expel the substrates when the transporter is opened, and prevents the passage of them when it is closed. Substrates are also affected by the random motion of the molecules of the solvent whose average kinetic energy is proportional to kBT. The associated force has been denoted by Fr in Equation (1). The radius of the transporter h(x, t) = (hmax − hmin)(x/L − 1/2) sin ωt + (hmax + hmin)/2 varies with position and time reaching maximum and minimum values, hmax and hmin, respectively. x = 0 is located at one extreme of the transporter.
The force arising from the uneven shape of the transporter must therefore be proportional to ΔA = A2 − A1, with A1 and A2 the cross-sectional areas at the entrance and at the exit of the transporter (see Figure 2) which depend on time. Studies on ion translocation in Ca2+-ATPase and in Na+/K+-ATPase have revealed that structural changes in a protein channel and their induced entropic forces contribute significantly to the transport of the ions (Rubi et al., 2017).
Changes in the concentration, or equivalently, in the chemical potential of the substrates on both sides of the plasma membrane generate a mass flow from high to low concentrations. The force associated with this effect, proportional to Δc = c2 − c1, with c1 and c2 the concentrations at the entrance and at the exit of the transporter (see Figure 2), is lower than the entropic force and is directed in the opposite direction, so the net effect is a flow of substrates against the gradient, an active transport due to the catalyst effect of the transporter.
Substrates are also affected by the thermal motion of the solvent which exert a random force on them. The thermal energy is kBT and the thermal random force is given by , with η(t) a Gaussian random quantity of mean zero and correlation < η(t)η(t′)>=δ(t−t′), T the temperature, kB the Boltzmann constant and the friction coefficient which is the inverse of the diffusion coefficient in kBT units (Gardiner, 2004).
The previous forces acting on the substrates capture the essential factors involved in the translocation process. In the model, we assume that the resulting velocity v is given through the Langevin equation (Gardiner, 2004)
It has been shown that entropic forces are given by (de Groot et al., 1963; Zwanzig, 1992; Reguera and Rubí, 2001; Kalinay and Percus, 2006; Vazquez et al., 2008; Rubi, 2019) and . In Equation (1), these forces act on the substrates at their time-dependent positions. The entropic force depends on the local radius h(x,t) of the channel (Reguera and Rubí, 2001) and has the direction of the cross-sectional area gradient, i.e., it contributes to expel the substrates, whereas the diffusive force has the opposite sign of the concentration gradient. A similar Langevin equation was proposed to model particle translocation through microfluidic channels promoted by an oscillating external potential (Tan et al., 2017).
To study how changes in conformation of the transporter affect the velocity of the substrate, we must model its geometry. A simple yet representative form is that of a conical shaped region of length L that oscillates periodically in time with a frequency ω (Carusela and Rubi, 2017, 2018), opening and closing to the intra/extra cellular environment, as sketched in Figure 2. The radius h(x, t) of the channel changes from a maximum value hmax to a minimum value hmin evolving in time as
The value x = 0 is located at one extreme of the transporter.
2.2. Computation Protocol
Measurements of translocation rates are performed over a time interval long enough to comprise many time periods of the conformational change cycle of transporters of the same kind. To obtain a representative value of the velocity of the substrates, we must thus average the instantaneous velocity v(t), given in Equation (1), in time and over an ensemble of identical transporters. The average velocity is thus obtained as
where < ... > means average over an ensemble of realizations or initial states of the system.
To calculate < v(t) >, we integrate Equation (1) using a stochastic Velocity-Verlet algorithm (Gränbech-Jensen and Farago, 2013). The computation method runs as follows. We consider a particle at the entrance of the channel, on the extracellular environment, and set its initial velocity according to a Boltzmann distribution, at T = 300K. We then calculate its local velocity according to Equation (1). Once the particle reaches the exit of the transporter into the intracellular medium, we calculate the time average , assuming that when the particle leaves the channel it cannot re-enter. This protocol is repeated for a large set of initial particle conditions at the transporter entrance. The average (< .. >) of the obtained results gives us the value of . In the model, we have considered that the force due to the concentration gradient in Equation (1) is practically constant along the transporter which means that . The ratio and therefore fμ typically takes values in the range (10−1 − 1) (Sperelakis, 2000; Tashiro et al., 2005; Chu et al., 2013). The value of the entropic forces falls in the interval (1 − 101), therefore which means that changes in conformation is the main mechanism that regulates transport.
Following the protocol described in the previous section, we compute the transport rate which is plotted in Figure 3 vs. . This quantity represents the ratio between the time in which the conformation of the transporter changes and the time that substrates take to diffuse through the transporter. The presence of resonant peaks at , practically independent of the values of fμ, reveals the occurrence of an amplification of the velocity of the substrates at a certain value of the oscillation frequency and therefore shows that transporters work under optimal transport conditions. The peaks become more pronounced when fμ increases or equivalently the rates are higher when the concentration of substrates increases, as observed for example in experiments for LeuT transporters (Fitzgerald et al., 2019). The existence of an optimal frequency for transport of particles in microfluidic devices subjected to an external oscillatory potential was also found in Tan et al. (2017).
Figure 3. Transport rate Γ vs. ratio between the time conformations take to change and the diffusion time: . The quantities , fμ, L, and ω which are defined in the text, have been measured in units of , kBT/L0 (T = 300 K), L0 = 10nm and , respectively. An enhancement of the velocity is observed at regardless of the values of fμ, thus showing that transporters work in an optimal scenario dominated by entropic forces induced by conformational changes.
Considering transporters with length of the order of L0 = 10nm and typical diffusion coefficients for membrane proteins of Eukaryotic cells, (Kaňa, 2013), we obtain values for Γ in the range of (10−1 − 103)seg−1. Our model also provides values of the rates in ion channels. In this case, one would expect that since particles are lighter they move faster. Typical values of for ions are in the range (103 − 104)μm2/seg, therefore Γ takes values between (105 − 106)seg−1. These model predictions are in good agreement with data found in the literature for transport rates of slow and fast protein channels in membranes (Ashcroft et al., 2009). At the maximum opening configuration of the transporter, the entropic force changes from a value of the order of fμ at the widest part, to about 20fμ at the narrowness. Our results then show that transporters work in a regime dominated by entropic forces induced by conformational changes.
Membrane proteins transporters regulate the efflux/influx of substrates across plasma membrane of cancer cells and play a paramount role in the efficiency of cancer therapies as they regulate the retention of anticancer drugs in the cells. It is believed that the effectiveness of chemotherapy may be largely dependent on the activity of transporters. Knowing which is the intimate mechanism that makes possible the passage of substrates through the cell membrane is therefore a matter of vital importance.
In this article, we have proposed a computational model that analyses the activity of transporters. The model shows that changes in the transport configuration produced by the energy of ATP hydrolysis (ABC transporters) or by a electrochemical gradient of ions (secondary transporters) give rise to a force that makes the substrates able to overcome the potential barrier generated by the constrictions in order to pass to the other side of the membrane. Such a force results from the variation of the space that the substrates have to move within the transporter when its configuration changes, which generates an entropy gradient.
Transport rates are computed by following a simulation protocol in which we first obtain the time average of the velocity of the substrates which is subsequently averaged over a set of initial conditions of the position of a substrate. Our model reproduces experimental values of the rates for different ABC and secondary transporters and shows that they depend on substrate concentration, in accordance with data reported in recent experiments (Fitzgerald et al., 2019). It can also be used to compute translocation rates in ion channels showing that they are greater than for transporters, as observed in the experiments (Ashcroft et al., 2009). We have shown that the entropic force is greater than that produced by a mere concentration gradient and is directed in the opposite direction which shows that transport is active.
How readily substrates cross the membrane depends on the frequency of oscillation ω, which regulates the entropic force and therefore the changes in the conformation of the transporter, on the size of the substrates which is involved in and on the size of the transporter L. The rates found scale with the combination of these quantities: which represents the ratio between the time in which conformations change and the time it takes for the substrates to diffuse. Our results could therefore be applied in general to membrane transport, to transport in microfluidic devices and peristaltic channels and to capture cellular heterogeneity beyond cancer cells (Pu et al., 2016; Łapińska et al., 2019; Rhia et al., 2020). Monitoring the entropic force by means of drugs could help to regulate transport activity leading to the design of better transport modulators that can be used in anti-cancer therapies.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
All authors have made a substantial contribution to the work and approved it for publication.
MFC thanks PIO-CONICET and PICT 18/2036, JMR wants to thank financial support of MICIU of the Spanish Government, under Grant No. PGC2018-098373-B-I00.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Ashcroft, F., Gadsby, D., and Miller, C. (2009). Introduction. The blurred boundary between channels and transporters. Philos. Trans. R. Soc. Lond. B Biol. Sci. 364, 145–147. doi: 10.1098/rstb.2008.0245
Asif, M., Usman, M., Ayub, S., Farhat, S., Huma, Z., Ahmed, J., et al. (2020). Role of ATP-binding cassette transporter proteins in CNS tumors: resistance-based perspectives and clinical updates. Curr. Pharm. Design 26, 4747–4763. doi: 10.2174/1381612826666200224112141
Bock, C., Zollmann, T., Lindt, K., Tampé, R., and Abele, R. (2019). Peptide translocation by the lysosomal abc transporter tapl is regulated by coupling efficiency and activation energy. Sci. Rep. 15:11884. doi: 10.1038/s41598-019-48343-6
Chu, X., Korzekwa, K., Elsby, R., Fenner, K., Galetin, A., Lai, Y., et al. (2013). Intracellular drug concentrations and transporters: measurement, modeling, and implications for the liver. Clin. Pharmacol. Ther. 94, 126–141. doi: 10.1038/clpt.2013.78
El-Awady, R., Saleh, E., Hashim, A., Soliman, N., Dallah, A., Elrasheed, A., et al. (2017). The role of eukaryotic and prokaryotic abc transporter family in failure of chemotherapy. Front. Pharmacol. 7:535. doi: 10.3389/fphar.2016.00535
Łapińska, U., Glover, G., Capilla-Lasheras, P., Young Andrew, J., and Pagliara, S. (2019). Bacterial ageing in the absence of external stressors. Philos. Trans. R. Soc. B 374:1786. doi: 10.1098/rstb.2018.0442
Liu, Y., Liu, Y., He, L., Zhao, Y., and Zhang, X. (2018). Single-molecule fluorescence studies on the conformational change of the abc transporter msba. Biophys. Rep. 4, 153–165. doi: 10.1007/s41048-018-0057-z
Loo, T. W., and Clarke, D. M. (2001). Determining the dimensions of the drug-binding domain of human p-glycoprotein using thiol cross-linking compounds as molecular rulers. J. Biol. Chem. 276:36877. doi: 10.1074/jbc.C100467200
Neumann, J., Rose-Sperling, D., and Hellmich, U. A. (2017). Diverse relations between abc transporters and lipids: an overview. Biochim. Biophys. Acta Biomemb. 1859, 605–618. doi: 10.1016/j.bbamem.2016.09.023
Olesen, C., Picard, M., Winther, A.-M. L., Gyrup, C., Morth, J. P., Oxvig, C., et al. (2007). The structural basis of calcium transport by the calcium pump. Nature 450, 1036–42. doi: 10.1038/nature06418
Payen, V. L., Hsu, M. Y., Rädecke, K. S., Wyart, E., Vazeille, T., Bouzin, C., et al. (2017). Monocarboxylate transporter mct1 promotes tumor metastasis independently of its activity as a lactate transporter. Cancer Res. 77, 5591–5601. doi: 10.1158/0008-5472.CAN-17-0764
Pu, Y., Zhao, Z., Li, Y., Zou, J., Ma, Q., Zhao, Y., et al. (2016). Enhanced efflux activity facilitates drug tolerance in dormant bacterial cells. Mol. Cell 62, 284–294. doi: 10.1016/j.molcel.2016.03.035
Reckzeh, E. S., Karageorgis, G., Schwalfenberg, M., Ceballos, J., Nowacki, J., Stroet, M. C., et al. (2019). Inhibition of glucose transporters and glutaminase synergistically impairs tumor cell growth. Cell Chem. Biol. 26, 1214–1228.e25. doi: 10.1016/j.chembiol.2019.06.005
Rhia, M., Stone, L., Łapińska, U., Pagliara, S., Masi, M., Blanchfield, J. T., et al. (2020). Fluorescent macrolide probes-synthesis and use in evaluation of bacterial resistance. RSC Chem. Biol. 1, 395–404. doi: 10.1039/D0CB00118J
Robey, R., Pluchino, K., Hall, M., Fojo, A., Bates, S., and Gottesman, M. (2018). Revisiting the role of abc transporters in multidrug-resistant cancer. Nat. Rev. Cancer. 7, 452–464. doi: 10.1038/s41568-018-0005-8
Sampedro-Núñez, M., Bouthelier, A., Serrano-Somavilla, A., Martínez-Hernández, R., Adrados, M., Martín-Pérez, E., et al. (2020). LAT-1 and GLUT-1 carrier expression and its prognostic value in gastroenteropancreatic neuroendocrine tumors. Cancers 12, 1–17. doi: 10.3390/cancers12102968
Scalise, M., Pochini, L., Galluccio, M., Console, L., and Indiveri, C. (2020). Glutamine transporters as pharmacological targets: from function to drug design. Asian J. Pharm. Sci. 15, 207–219. doi: 10.1016/j.ajps.2020.02.005
Tashiro, M., Tursun, P., and Konishi, M. (2005). Intracellular and extracellular concentrations of Na+ modulate Mg2+ transport in rat ventricular myocytes. Biophys. J. 89, 3235–3247. doi: 10.1529/biophysj.105.068890
Vazquez, M.-V., Berezhkovskii, A. M., and Dagdug, L. (2008). Diffusion in linear porous media with periodic entropy barriers: a tube formed by contacting spheres. J. Chem. Phys. 129:046101. doi: 10.1063/1.2955447
Xu, L., Chen, J., Jia, L., Chen, X., Awaleh Moumin, F., and Cai, J. (2020). SLC1A3 promotes gastric cancer progression via the PI3K/AKT signalling pathway. J. Cell. Mol. Med. 24, 14392–14404. doi: 10.1111/jcmm.16060
Yamada, Y., Yoshimatsu, K., Yokomizo, H., Okayama, S., and Shiozawa, S. (2020). Expression of ATP-binding cassette transporter 11 (ABCC11) protein in colon cancer. Anticancer Res. 40, 5405–5409. doi: 10.21873/anticanres.14549
Keywords: cancer transporters, membrane transport mechanisms, entropic forces, Langevin equation, cancer therapies
Citation: Carusela MF and Miguel Rubi J (2021) Computational Model for Membrane Transporters. Potential Implications for Cancer. Front. Cell Dev. Biol. 9:642665. doi: 10.3389/fcell.2021.642665
Received: 16 December 2020; Accepted: 04 February 2021;
Published: 22 February 2021.
Edited by:Maria E. Mycielska, University Medical Center Regensburg, Germany
Reviewed by:Leonardo Dagdug, Metropolitan Autonomous University, Mexico
Stefano Pagliara, University of Exeter, United Kingdom
Copyright © 2021 Carusela and Miguel Rubi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: J. Miguel Rubi, email@example.com