Journal Pre-proof On the electrical properties of collagen macromolecule solutions: Role of collagen-water interactions
M. Ashoorirad, M. Saviz, A. Fallah PII:
S0167-7322(19)35154-2
DOI:
https://doi.org/10.1016/j.molliq.2019.112344
Reference:
MOLLIQ 112344
To appear in:
Journal of Molecular Liquids
Received date:
14 September 2019
Revised date:
13 December 2019
Accepted date:
17 December 2019
Please cite this article as: M. Ashoorirad, M. Saviz and A. Fallah, On the electrical properties of collagen macromolecule solutions: Role of collagen-water interactions, Journal of Molecular Liquids(2018), https://doi.org/10.1016/j.molliq.2019.112344
This is a PDF file of an article that has undergone enhancements after acceptance, such as the addition of a cover page and metadata, and formatting for readability, but it is not yet the definitive version of record. This version will undergo additional copyediting, typesetting and review before it is published in its final form, but we are providing this version to give early visibility of the article. Please note that, during the production process, errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
© 2018 Published by Elsevier.
Journal Pre-proof
On the electrical properties of collagen macromolecule solutions: Role of collagen-water interactions M. Ashoorirad1, M. Saviz1 , A. Fallah1,* 1
Department of Biomedical Engineering, Amirkabir University of Technology (Tehran Polytechnic), Tehran, Iran.
Abstract
re
-p
ro
of
Collagen is the most abundant protein of the extracellular matrix in connective tissue in the mammalians with potential applications in the coverage of cardiovascular prostheses, regenerative medicine, pharmaceutical industry, controlled drug delivery, and proliferation, migration and cellular differentiation. In this research, the electrical properties (conductivity and permittivity) of type I collagen solution with different concentrations were measured in the frequency range 10KHz - 5MHz and electrical analysis of the measurement results was performed. Analysis shows that special interactions with water molecules have a significant effect on the electrical properties of collagen solutions beyond what is predicted considering the rules of physical composition, such as the Maxwell-Garnet model of mixing rules. Very small variations in the amount of collagen in water, significantly changed solution conductivity and dielectric permittivity. Finally, to simplify the representation of the impedance measurements and to better understanding of the electrical behavior of the solutions, an equivalent electrical circuit model was developed to describe the measured data. The results obtained from the fitted electrical model were in good agreement with the experimental results.
lP
Keywords: Collagen; Electrical Properties; Bioimpedance Spectroscopy; Electrical Circuit Model. 1. Introduction
Jo
ur
na
The superiority of collagen as a biomaterial essentially depends on the fact that collagen is a component of the extracellular matrix and therefore is considered as an intrinsic element of the body [1– 5]. Collagen, one of the closest approximation to extracellular matrix [6], is a complex supramolecular structure, with very diverse morphologies across different tissues which give them a wide range of biological functions in the body [3,4]. Collagen is the most important protein in the connective tissues and is the most abundant protein in mammals [7–10]. In the vertebrates, collagen is a major component of the specialized and non-specialized connective tissues, which accounts for about a quarter of the body's total protein, about three-quarters of the human dry skin, more than 90% of the tissues of the tendon and the cornea, and about 80% of organic matter of bones [5]. Up to now, 29 types of collagen have been identified [5]. The difference in collagens is due to the difference in the triple polypeptide chains, the difference in helical length, the different interruptions in the helix and the difference in the helical terminals [5]. Generally, collagen is divided into two types: fibrous and non-fibrous. Collagen types I, II, and III are the most common fibrillar collagens. Collagen type I, is the predominant collagen in most tissues of higher order animals [11] and accounts for 80-85% of collagen in the body [12]. Collagen biopolymer is the most abundant protein in the mammalian body with a molecular weight of 300 kg/mol [13], which consists of amino acid sequences. Collagen contains different amino acids. Amino acids such as Glysine, Prolyne, Hydroxyproline, Glutamine, Arginine, Aspartic acid and Hydroxylysine are present in collagen *
Corresponding Authors: Ali Fallah, Ph.D., Department of Biomedical Engineering, Amirkabir University of Technology, Tehran, Iran, Email:
[email protected]
1
Journal Pre-proof with more content [14]. Hydroxyproline and Hydroxylysine are not significantly present in proteins other than collagen [14].
of
In other words, amino acids with a hydroxyl group are specific to collagen and do not exist in a sufficiently great in other proteins. The presence of these hydroxyl groups in the collagen macromolecule leads to the formation of hydrogen bonds between collagen fibers with each other, and collagen fibers with the water molecule which results in increase of the conductivity of the collagen solution [15]. By linking these amino acids via the peptide bonds, a chain of amino acids is formed. The collagen macromolecule consists of three polypeptide chains that are weaved together with hydrogen bonds and forms a triple helical structure [4,16]. Hydrogen bonds hold amino acids firmly in a mesh of fibers. Each collagen helix consists of approximately 1.5 nm in diameter and 300 nm in length [17,18]. The amino acids in collagen macromolecule have a regular arrangement in each of the three chains. This sequence often follows the pattern of Gly-Pro-X or Gly-X-Hyp, in which Gly is glysine, Pro is hydroxyproline, Hyp is hydroxyproline, and X can be any of different amino acids.
re
-p
ro
The Gly-Pro-Hyp chain is the most stable amino acid sequence in collagen. Collagen triple helical structure stabilizes with hydrogen bonds between carbonyl and hydroxyl groups of hydroxyproline [19]. Collagen is a natural protein used as a dielectric gate for organic thin film transistors (OTFTs) in vacuum and air conditions. The high mobility value of the field effect in OTFTs is attributed to the interaction between water molecules and hydroxyl groups in collagen [19] which leads to conductivity enhancement. Hydroxyl groups are found in the hydroxyproline and hydroxylysine amino acids in the collagen structure.
Jo
ur
na
lP
Generally, polymers are classified as polar and nonpolar depending on their chemical structure. Nonpolar polymers are usually insulating materials with a pseudo-constant permittivity of 2.5 to 3 [20]. The structure of these proteins is made of carbon and hydrogen atoms that are symmetrically arranged along a carbon chain. Polar polymers are composed of molecules including electronegative atoms such as nitrogen, oxygen, chloride and fluoride, and create dipoles within the polymer. The higher dielectric coefficient is a characteristic of these polymers [20]. Proteins are polar biopolymers whose polar index can be weaker or stronger depending on the lateral chain of amino acids present in them [20]. Contrary to complex molecular chains, proteins provide some interesting points, especially their charge content, which results in electrical properties (conductivity and dielectric permittivity) of proteins. The electrical properties of proteins play an important role in evaluating their functions. Many studies have been performed to investigate the electrical, mechanical and structural properties of proteins [20–23]. Electrostatic interactions are important for the overall stability of the protein structure [14], and in total, are considered as the contributions of van der Waals forces, hydrogen bonds, covalent bonds, and hydrophobic interactions [24]. Therefore, these bonds not only stabilize the protein structure, but also contribute to the electrical properties of proteins. Hence, proteins are considered as polarizable materials that can be characterized by their electrical properties [15]. Polarizing the components of a material is due to a change in the redistribution of electrical charges, when the material is exposed to an electric field. Polarization can be due to various effects such as charge accumulation at surfaces of multilayer materials with different dielectric permittivity (interfacial polarization), bipolar orientation or electronic and atomic ( or ionic) polarization [20]. Each of these polarization effects contribute to overall electrical properties. The electrical properties of the materials depend on their chemical composition, in particular the stable dipole moments associated with water molecules. The higher the number of dipoles inside the solution, the higher the total charge transfer and the higher conductivity of the solution [25]. Water molecules enter the amino acids chains of protein and form hydrogen bonds through their hydroxyl group. 2
Journal Pre-proof With the penetration of water into the protein network, the structure and mobility of the protein peptide chains change [20]. Therefore, water molecules can effectively influence the electrical properties of proteins. It has been reported in [26] that the formation of hydrogen bonds between DNA bases enhance conductance. The shorter bond length conforms stronger hydrogen bond and results in more conductance.
of
Like other proteins, for collagen, amino acids and hydroxyl groups form hydrogen bonds with water molecules and increase the conductivity of collagen solution. In addition, quasi-ring structures have improved conductivity due to bonds [26]. Some of the amino acids in collagen macromolecules, such as phenylalanine, tyrosine and histidine, have ring structures with bonds, and therefore, can lead to increased conductivity due to free electrons. Since collagen protein consists of polar repeated units of amino acids (-NH-CR-CO) and other polar groups such as COO-, OSO3- and OH- [13], electrical methods can be used to analyze this biomaterial. These methods are used to study the effects of water molecules and frequency of electric field on the electrical properties of collagen. In this study, the electrical properties of type 1 collagen solutions are measured and investigated.
lP
re
-p
ro
In this paper, the electrical properties of type I collagen (the most abundant collagen type [5]) in soluble phase are investigated and the effect of increasing collagen concentrations on these properties is studied and analyzed. To theoretically confirm the role of water in the electrical properties of collagen protein, Maxwell-Garnet model for the mixing rules of physical composition of materials was implemented in MATLAB environment. Finally, an equivalent electrical circuit model was fitted to measured data at different concentrations. The advantage of this model is that by determining the various parameters of the model for each concentration, better intuition can be obtained of the electrical events occurring in collagen solution under the influence of the electric field. 2. Materials and methods 2.1. Preparation of collagen solutions
ur
na
Dry Type 1 Collagen, obtained from calf skin, was purchased from the Pasteur Institute (Tehran). A total of 20 mg of collagen was dissolved in 2 ml of acetic acid 2.5% (v/v in deionized water) using a magnetic stirrer and therefore, 2 ml of collagen solution at a concentration of 10 mg/ml was prepared. Using this solution, collagen solutions with the concentration of 250, 500, 750 and 1000 were prepared by diluting the original collagen solution with deionized water.
Jo
2.2. Biosensor fabrication and impedance measurements To measure the electrical properties of the collagen solutions, flat gold electrodes with size of 1.9 cm × 1 cm were constructed on a FR4 substrate. With the help of two other FR4 plates with dimensions of 1.7cm × 1cm, two other walls of the enclosure measuring the electrical properties of the collagen solution were prepared and carefully sealed. After preparing the compartment containing electrodes, existence of any leakage on it was tested. Then, it was washed well with detergent and rinsed with deionized water. After drying, it was sprayed with 70% ethanol and used for measurements after alcohol evaporation. In order to keep the sensor stable during different measurements, the sensor was fixed on a circular base of Plexiglas material (Figure 1. ). Impedance measurements of collagen solutions were performed at room temperature ( ) and frequency range of 10 KHz to 5 MHz at 200 linearly spaced frequencies using a Hioki IM3570 impedance analyzer and an average of four measurements per frequency. Before data collection the Hioki IM3570 was calibrated using the open/short/load procedure based on its user manual. Further details of this calibration are provided in [27]. Also, as far as possible, the interconnecting wires were twisted pairs to minimize the induction of electromagnetic noise. 3
Journal Pre-proof Using a sampler, 3.5 ml of collagen solution was slowly poured in the sensing enclosure in such a way to prevent bubbles formation inside the solution. Then, a sinusoidal stimulation current of 100 A amplitude was applied to the electrodes at a frequency range of 10KHz to 5MHz. These measurements were performed using 3.5 ml of collagen solution with four concentrations of 250, 500, 750 and 1000μg/mL to study the effects of increasing collagen fibers concentration on the conductivity and permittivity of the solutions. Measurement errors were less than 1%. This means that the difference between the measured impedances among the at least three samples of each concentration has been less than 1%. The measurements errors of impedance magnitude and phase for different concentrations have been indicated using error bars in Figure 2. 2.3. Statistical Analysis
of
All the experimental data are presented as the mean of three repetitions. At each repetition, three samples were measured for each concentration. Statistical analysis was done by student’s T-test. Differences were statistically significant at P-values less than 0.01.
ro
3. Results and discussion 3.1. Electrical properties of proteins
lP
re
-p
The response of a material to an electric field is determined with its electrical properties (conductivity and dielectric permittivity). In aqueous media, conductivity is mainly due to the presence of charged ions in the environment. The higher the number of these charged factors, the greater the expected conductivity. The dielectric permittivity is defined as the power of the material in energy storage (electrical capacitor). The greater the dielectric permittivity of a material, the greater the ability of the material to store energy when exposed to an electric field, in other words, its electrical capacity is increased.
ur
na
Since proteins are composed of polar repeating units of amino acids and neutral, polar, and charged side chains [13], they can be considered as polarizable materials. In a peptide bond in a protein molecule, the two amino acid residues are linked by removing a water molecule from the carboxyl group from an amino acid and an amino group from another acid through a covalent bond, and charged agents are generated in the protein molecule. Due to this natural charge, proteins are sensitive to electrical excitation and exhibit specific electrical properties under certain conditions. As the electric properties of proteins depend on their composition and chemical structure, environmental parameters and material processing also affect the charge and consequently the electrical properties of protein solutions [20].
Jo
3.2. Measuring electrical properties of type I collagen solution The Nyquist diagram of the measurements is shown in Figure 3. The experimental curves demonstrate a distinguishable electrical sensitivity to even small concentrations of collagen. An elevation of collagen concentration causes a significant change in the real part of the impedance which results in increasing conductivity. For a more accurate comparison of the electrical properties of collagen solution at different concentrations, the electrical properties (conductivity and dielectric permittivity) of collagen solution were obtained at different concentrations by calculating the geometric coefficient of the measuring electrode. Considering the dimensions of the electrode used in the sensor fabrication (section 2.2), the geometric factor of the measuring electrodes was calculated and the electrical complex conductivity was obtained. Specifically [28,29]: (1) (2)
4
Journal Pre-proof Collagen Protein Electrical Model
Cdl
Ccol
Cdl
Rs Rct
IM 3570 Impedance Analyzer
Rcol
Rct
Collagen fibers in solution
L2000 4-Terminal Probe
USB Port
𝑃𝑜𝑡
𝐻𝑃𝑜𝑡 𝑢𝑟𝑟𝑒𝑛𝑡
1 cm
1.9 cm
lP
re
-p
ro
of
𝐻
(b)
Jo
(a)
ur
na
Figure 1. Measurement setup for measuring electrical properties of type I collagen protein solution.
Figure 2. Measurements errors of (a) impedance magnitude and (b) impedance phase for different concentrations of collagen solution using error bars representation. Differences were statistically significant at P-values less than 0.01 and measurements error was less than 1%.
5
Journal Pre-proof Where Z is the measured impedance, A is the area of the electrodes, and L is the distance between the electrodes. The
ratio is the geometric factor of the measurement chamber. The
is the vacuum
lP
re
-p
ro
of
permittivity and is the relative permittivity of the collagen solution. is complex conductivity. The resulting conductivity is a complex value, its real part is the conductivity ( ) of collagen solution, and its imaginary part is a factor of the relative dielectric permittivity of collagen solution ( ). Conductivity and relative dielectric permittivity of collagen solutions are presented in Figure 4.
Jo
ur
na
Figure 3. Nyquist diagram of collagen protein solutions with different concentrations at the frequency range of 10KHz to 5MHz. Differences were statistically significant at P-values less than 0.01 and measurements error was less than 1%.
125KHz
(a)
(b)
Figure 4. (a) Conductivity and (b) dielectric permittivity of collagen solution with different concentrations at the frequency range of 10KHz to 5MHz. For frequencies more than 125KHZ, relationship between concentration and dielectric permittivity is inversed. Differences were statistically significant at P-values less than 0.01 and measurements error was less than 1%.
As shown in Figure 4, the conductivity of the collagen solutions significantly increases with increasing concentrations. Various factors contribute to the increase in conductivity by increasing the concentration of collagen fibers in the solution. Electrical conduction of collagen fibers in a soluble phase can be due to the electrical interaction between bipolar molecules of water and active chemical groups on 6
Journal Pre-proof
ro
of
collagen fibers. The intertwined structure of collagen is retained together by electrostatic and hydrophobic bonds trapping large quantities of liquids and allowing the exchange of ions and metabolites with surrounding tissues [5]. Besides, with increasing number of hydroxyl groups in the polymer, mobility increases [30]. Hydroxyl groups are significantly present in the amino acid chains of collagen helical macromolecules. Therefore, the existence of these functional groups increases mobility which leads to an increase in electrical conductivity of collagen solution. On the other hand, O- and H3O + (hydronium ion) groups are produced in the collagen solution using the following reaction:
-p
Figure 5. Interaction of the hydroxyl group of hydroxyproline amino acid in the peptide chains of collagen fibers with the water molecule and formation of hydronium ion.
lP
re
The residue of the hydroxyproline amino acid with the O- group in the collagen macromolecule is immobilized due to its large size [19]. In contrast, H3O + ions are able to move freely and easily inside the solution and therefore, when the electric field is applied to collagen solution, these ions move to the negative pole of the field, resulting in high conductivity. Since the number of these ions increases with increasing collagen fiber content (increasing the concentration of collagen solution), the measured conductivity would be proportional to the concentration of collagen solution and increases (Figure 4(a)).
Jo
ur
na
Except for hydroxyl groups of amino acids in collagen macromolecules that form hydrogen bonds with water molecules and increase the conductivity of collagen solution, another feature of the chemical structure of amino acids of the collagen macromolecule also improves conductivity. Ring structures lead to improved conductivity due to their -bonds [26]. Some of the amino acids in collagen macromolecule, such as phenylalanine, tyrosine and histidine, have a benzene ring with -bonds and, due to free electron, increase the conductivity of the collagen solution. The higher the concentration of collagen solution, the higher the number of these amino acids inside the solution, resulting in greater conductivity. However, in our measurements on the dry collage (see supplementary data) there was observed no difference between bare electrodes and electrodes with the same mass of collagen related to the different concentrations. Meanwhile, for the same amount of dry collagen dissolved in acetic acid and diluted with deionized water, very significant changes in electrical properties were observed. Thus, we think that this factor may have no important effect on the electrical properties of the collagen solution. On the other hand, the formation of bipolar moments in the solution, essentially contributes to the reduction of the total impedance of the solution [25]. Since the hydroxyproline residues with the O- group (Figure 5) and other charged groups in collagen fibers form local electric dipoles inside the solution, the collagen solution would present dielectric property when exposed to an electric field. Since the number of dipoles in the solution is proportional to the collagen concentration, changing the concentration of collagen changes the dielectric behavior of the solution. In the low-frequency range (f < 125KHz), the dielectric coefficient increases with increasing concentration of collagen solution. In the higher frequency range (f > 125KHz), dielectric permittivity of the collagen solution decreases with
7
Journal Pre-proof increasing the concentration of collagen solution. A more detailed study of the effect of frequency on the electrical properties of collagen solution is given in section 3.4. 3.3. Maxwell-Garnet mixing model To theoretically demonstrate that the measured electrical properties for collagen solutions are due to the presence of water and the formation of hydrogen bonds and other chemical interactions between the charged agents present in the collagen amino acid chains with the water molecules, the electrical properties of collagen with different fractional volumes were obtained using the Maxwell-Garnet mixing model. In the Maxwell-Garnet model, only physical interactions between the various components of a mixture are considered.
ro
of
Dielectric mixing rules are algebraic formulas that can be used to calculate the effective permittivity of a mixture as a function of constituent permittivities, their fractional volumes, and some other possible parameters that characterize the microstructure of the mixture. The mixture can be discrete, meaning that the homogeneous inclusions are embedded in another homogeneous environment. The size of the inclusions in the mixture needs to be much smaller than the wavelength [31].
re
-p
In practice, the inclusions can take many different shapes. The only shape for which simple analytical solutions are available is the ellipsoid shape that fortunately covers many practical situations [31]. In the case of collagen molecule, which is an inclusion in the 2.5% acetic acid medium, it is properly approximated as an ellipsoid shape.
(3)
na
lP
The properties of a dilute mixture with a fractional volume less than 0.1 can be correctly calculated using the Maxwell-Garnet mixing law [31]:
Jo
ur
Where is a dimensionless quantity, the volume fraction of the inclusions in the mixture and subscripts e and i refer to the environment and inclusion, respectively. This formula is widely used in various fields. If f → 0 then the inclusion phase is vanished and and if f → 1 then the environment phase is vanished and . To calculate the effective permittivity of the mixture of collagen inclusions in 2.5% acid acetic medium using mixing rules, given that the collagen fibers are 1.5 nm in diameter and 300 nm in length, the collagen molecules approximately form randomly oriented ellipsoid inclusions with axial dimensions of 0.75 nm, 0.75 nm and 150 nm. In this calculation, the permittivity of solid collagen fibers is assumed to be 2.5 [20] respectively at room temperature. The permittivity of 2.5% acetic acid as the base material (medium) was considered as 59 (see supplementary data). To calculate the fractional volume of collagen inclusions in acetic acid medium, assuming the protein density is 1.35 gr/mL [32], for the collagen solution at a concentration of 250 g/mL fractional volume would be: (4) To obtain the fractional volumes for concentrations of 500, 750 and 1000 , this value is multiplied by 2, 3, and 4, respectively. The corresponding permittivity was calculated and plotted in MATLAB environment using the mixing rules ( Figure 6 ). As can be seen from these graphs, for different fractional volumes corresponding to different concentrations of collagen solution, the 8
Journal Pre-proof permittivity variation is negligible respect to medium permittivity. Thus, in terms of mixing rules and only considering physical interactions without regard to chemical interactions, there is no significant difference between the electrical properties of dry collagen with different fractional volumes. However, in the measured data for collagen solutions, although the mass of collagen varies similarly, the electrical properties change significantly. This fact is due to the presence of water and the chemical interactions between the water molecules and the charged agents present in the structure of collagen fibers, which significantly increases the conductivity and permittivity.
na
lP
re
-p
ro
of
Furthermore, in our measurements on dry collagen, no significant relations were observed between recorded impedance characteristics of bare electrodes and electrodes coated by different amount of dry collagen corresponding with different concentrations of the collagen solutions (see supplementary data).
Jo
ur
Figure 6: Permittivity of solid state collagen fibers at different fractional volumes corresponding to the collagen solutions at different concentrations calculated by the Maxwell-Garnet mixing rule in MATLAB environment at frequency range of 10KHz to 5MHz.
3.4. The effect of frequency on the electrical properties of collagen solution Frequency is an important factor that influences the electrical behavior of different materials, including proteins. As the frequency increases, the dipoles in the protein are not able to follow faster changes in the electric field polarity, and as a result, the effect of these dipoles decreases with increasing frequency, and the dielectric permittivity of the protein decreases [20]. In performed measurements on collagen protein with different concentrations, this decreasing trend of the dielectric permittivity respect to the frequency is clearly seen (Figure 4). As shown in Figure 4 (b), the dielectric behavior of collagen solution in low and high frequencies is quite different. For frequencies below 125 KHz, with increasing collagen concentration, the relative dielectric permittivity is increased, while for frequencies higher than 125 KHz, with increasing the concentration of collagen solution, the dielectric permittivity decreases. At low frequencies (<125KHz), electrical dipole in collagen macromolecules have the opportunity to be aligned with the electric field before changing the direction of the field, and their overall results create a larger capacitor, but at higher frequencies (>125 KHz) because of the rapid changes in the electric field, the electrical dipoles in the 9
Journal Pre-proof structure of collagen fibers do not have the opportunity to rotate and align with the electric field and results in decreasing dielectric permittivity of collagen solution respect to its concentration. As the concentration of collagen fibers in the solution is increased, due to the greater number of fibers and their more intertwining, a smaller number of dipoles in the solution are capable of aligning to the electric field and exhibiting their capacitive (dielectric permittivity) properties. This leads to a decrease in the dielectric permittivity (capacitance) at higher frequencies (> 125KHz) respect to collagen concentration increasing. Due to the large size of charged residues in the collagen fibers, they immobilize in the solution at higher frequencies and result in a further reduction in the dielectric permittivity by increasing the concentration. By increasing the concentration of collagen solution, the number of large residues would increase and the dielectric permittivity of the solution further decreases.
ro
of
On the other hand, as shown in Figure 4(a), as the frequency increases, the conductivity increases with respect to collagen solutions concentration. This is due to the decrease in dielectric properties (capacitance) of collagen solution respect to frequency increment. By decreasing the capacitance, the overall impedance is reduced and the conductivity is increased.
-p
3.5. Fitting an electrical equivalent circuit to measured bio-impedance data
na
lP
re
Bioimpedance or electrical impedance of biological samples describes the electrical properties of these materials. The purpose of collecting these measurements is to provide details of electrochemical structures and electrical processes within a tissue or material being studied [33]. The bioimpedance spectrum is a safe, noninvasive, and inexpensive way to measure the electrical responses of living tissues and other biomaterials that can reflect the physiological and pathological status of the tissue or biomaterial being measured. To measure bioimpedance, a sinusoidal current (or voltage) is applied through the stimulation electrodes to the measured biomaterial sample, and its response is measured as voltage (or current) through the recording electrodes. By dividing the recorded response to the stimulation applied to the measured sample, a complex quantity is obtained called bioimpedance.
Jo
ur
To simplify the representation of the impedance spectrum measurements and to better understand the electrical behavior of the samples being measured, an equivalent electrical circuit schematic is desirable; electrical circuits models are regularly used to fit and interpret bioimpedance data [33–36]. In this method the measured spectrum is fitted on to an equivalent electrical circuit model which is constructed based on resemblances between the electrical circuit model and the measured physical phenomena (e.g. double layers). Thus, the spectral behavior is confined to several descriptive parameters of the model reducing a bioimpedance dataset from the potentially hundreds of data points to several parameters (dependent on the number of circuit components in the model). In this study, an equivalent circuit model is fitted to measured bioimpedance data and can be used to describe and analyze the electrical behavior of collagen protein solution. Due to the measured electrical properties of the collagen protein solution mentioned in the previous sections, this protein solution has both ohmic and capacitive components. Therefore, a conventional parallel RC circuit could model its electrical behavior. The analyzed collagen protein solutions are electrolytes, and they tend to induce a specific interfacial layer near to the sensor electrodes that is named a double layer ( ). The sensor electrode becomes polarized with regard to the collagen solution that attracts the ions of contrary charge. This effect can be illustrated by a double layer capacitor and a faradic resistance at the electrodeelectrolyte interface. The double layers of the electrode-solution interface are also conventionally modeled as parallel RC sections in series with the solution. As shown in Figure 1. the electrical equivalent 10
Journal Pre-proof circuit is symmetrically modeled on the basis of the electrodes structure and electrical behavior of collagen protein solution. The parameters , , ، , , and stand for double layer capacity, charge transfer resistance, bulk solution resistance, ionic conduction and dielectric behavior of collagen protein solution, respectively. Since the equivalent electrical circuit models for both electrodes (paired parallel ) are similar and in series, they can be considered equivalent to a single combined parallel circuit for modeling the electrode-solution behavior. Therefore, the equivalent circuit that was considered to fit the measured bioimpedance data is simplified as shown in Figure 7. The equivalent impedance of the circuit is described in Eq. (5).
of
Collagen Protein Electrical Model
Cdl
lP
Rct
re
-p
Rs
ro
Ccol
Rcol
na
Figure 7. Equivalent electrical circuit model fitted to measured bioimpedance data to analyze impedance characteristics of collagen protein solutions.
3.5.1 Fitting algorithm for determination of electrical circuit model parameters
Jo
ur
A complex nonlinear least square algorithm in MATLAB software was used to find model parameters that cause best agreement between measured impedance spectrum and model (estimated) spectrum. This algorithm attempts to minimize the distance between the measured bioimpedance spectrum and the estimated impedance with the model by solving an optimization problem: 𝑛∑
(6)
𝑡
Where is the vector of bioimpedance parameters ( ، ، ، ، ), is the estimated impedance (Figure 7) calculated using , 𝑡 is the measured bioimpedance, and 𝑡 are the calculated and measured dataset at frequency , and is the total number of data points in the measured dataset. The Eq.(6) is simultaneously applied on real and imaginary parts of impedance dataset. The optimization algorithm aims find bioimpedance parameters in a way to ideally reduce the least square error to zero or realistically to the level of error within the dataset [35]. Several factors such as an incorrect model for the dataset being fitted, poor estimation for the initial values, and noise may cause the algorithm to fail on converging to a useful fit. The measured and fitted impedance spectroscopy of type 1 collagen solution at concentration range of 250-1000 g/mL were plotted as a Nyquist curve in Figure 8. The values of equivalent electrical circuit 11
Journal Pre-proof model parameters for different concentrations are given in Table 1. Fitting quality is indicated by sum of squared errors (SSE) at the final row of the table. SSE is the sum of squared differences between real and imaginary parts of the measured and fitted impedances. It is a measure of the discrepancy between the measured data and an estimation model. SSE is defined as below:
SSE=∑
(7)
Where, , , , , and N are the real part of the measured impedance, real part of the estimated impedance using the proposed electrical model, imaginary part of the measured impedance, imaginary part of the estimated impedance using the proposed electrical model, and N is the number of frequencies in which the measurements have been performed, respectively.
lP
re
-p
ro
of
As it can be seen in Figure 3, the elevation of the collagen concentration decreases the semicircle diameter of the Nyquist diagram. From the point of the electrochemical impedance spectroscopy, this observation can be explained as a reduction of charge transfer resistance ( ) in electrical model proposed in Figure 7 and is consistent with the fitted values. From experimental data (Figure 4), the collagen solutions with various concentrations have specific electrical properties (conductivity and dielectric permittivity) which can be modeled by a RC circuit ( ) in series with bulk solution resistance ( ). As it can be seen from the values of the model parameters in Table 1, with increasing collagen protein concentration, the value of decreases which means that with elevation in collagen concentration the conductivity improves which is consistent with the measured results presented in Figure 4(a). The value of is increased respect to collagen concentration which is again consistent with our understanding of the measured data (Figure 4(b)). is bulk solution resistance which is affected by both ohmic and capacitive components.
750 g/mL
1000 g/mL
1227
1274
1193
1039
1785
752.3
294.2
94.13
22.46
47.32
142.5
733
98.25
86.35
56.64
42.12
𝑛
235.07
206.99
305.63
366.67
SSE
0.001
0.003
0.004
0.009
𝑝
250 g/mL
ur
500 g/mL
Jo
Parameters of Fitted Electrical Circuit Model
na
Table 1: Calculated parameters for the electrical circuit model fitted to the measured impedance data for different concentrations of collagen protein.
12
Journal Pre-proof
(b)
lP
re
-p
ro
of
(a)
(d)
ur
na
(c)
Jo
Figure 8. Nyquist diagram of measured and fitted data using electrical equivalent circuit presented in Figure 7 for collagen solution with concentrations of (a) 250μg/mL, (b) 500μg/mL, (c) 750μg/mL, and (d) 1000μg/mL.
4. Conclusion
In this paper, the electrical properties of collagen protein were investigated in soluble phase at the frequency range of 10KHz to 5MHz and an equivalent electrical circuit was fitted to the measured data using the least squares algorithm. Proteins are sensitive to electrostatic methods due to their polarization, and exhibit special electrical properties when subjected to an electric field. Because of this polarization property, impedance spectroscopy can be used as an effective method for protein analysis. Simulations results using the Maxwell-Garnet model of mixing rules for physical composition of materials, showed that only with respect to the rules of physical composition there is no difference between the different fractional volumes of collagen. However, the same weight of collagen in soluble state, significantly changed its conductivity and dielectric permittivity. Hence, the main cause of the variation in electrical properties of collagen solution appears to be the formation of hydronium ions due to the interaction of water molecules with the collagen macromolecule. The number of dipoles in the solution is proportional to the collagen concentration. Hence, changing the concentration of collagen changes the dielectric behavior of the solution. 13
Journal Pre-proof Conflict of interest statement We, Masoomeh Ashoorirad, Mehrdad Saviz, and Ali Fallah certify that there is no conflict of interest. Acknowledgment This work was supported by Iran National Science Foundation (INSF) through the financial agreement number 97015120.
References J. Štancl, J. Skočilas, A. Landfeld, R. Ţitný, M. Houška, ELECTRICAL AND THERMODYNAMIC PROPERTIES OF A COLLAGEN SOLUTION, Acta Polytech. 57 (2017) 229–234.
[2]
M. Gauza, L. Kubisz, Effect of water content on electrical conductivity of fish skin collagen, J. Non. Cryst. Solids. 357 (2011) 686–690.
[3]
D. Gopalakrishnan, V. Balachandar, A.C. Kumbharkhane, R. Sampathkumar, Dielectric relaxation studies of collagen–surfactant complexes in aqueous buffer solution, Int. J. Biol. Macromol. 138 (2019) 215–223.
[4]
R. Gendron, M.R. Kumar, H. Paradis, D. Martin, N. Ho, D. Gardiner, E.F. Merschrod S, K.M. Poduska, Controlled cell proliferation on an electrochemically engineered collagen scaffold, Macromol. Biosci. 12 (2012) 360–366.
[5]
A. Sorushanova, L.M. Delgado, Z. Wu, N. Shologu, A. Kshirsagar, R. Raghunath, A.M. Mullen, Y. Bayon, A. Pandit, M. Raghunath, The collagen suprafamily: from biosynthesis to advanced biomaterial development, Adv. Mater. 31 (2019) 1801651.
[6]
E. Nam, W.C. Lee, S. Takeuchi, Formation of highly aligned collagen nanofibers by continuous cyclic stretch of a collagen hydrogel sheet, Macromol. Biosci. 16 (2016) 995–1000.
[7]
M.D. Shoulders, I.A. Guzei, R.T. Raines, 4‐Chloroprolines: Synthesis, conformational analysis, and effect on the collagen triple helix, Biopolym. Orig. Res. Biomol. 89 (2008) 443–454.
[8]
K. Lewandowska, A. Sionkowska, S. Grabska, Chitosan blends containing hyaluronic acid and collagen. Compatibility behaviour, J. Mol. Liq. 212 (2015) 879–884.
[9]
A. Sionkowska, J. Skopinska-Wisniewska, M. Wisniewski, Collagen–synthetic polymer interactions in solution and in thin films, J. Mol. Liq. 145 (2009) 135–138.
[10]
H. Li, D. Wang, S. Li, B. Liu, L. Gao, Sustained Release of BSA from a Novel Drug Delivery Matrix–Bullfrog Skin Collagen Film, Macromol. Biosci. 4 (2004) 454–457.
[11]
C. Dong, Y.L.- Polymers, undefined 2016, Application of collagen scaffold in tissue engineering: recent advances and new perspectives, Mdpi.Com. (n.d.). https://www.mdpi.com/20734360/8/2/42 (accessed 15 December 2018).
[12]
F.F. Felician, C. Xia, W. Qi, H. Xu, Collagen from marine biological sources and medical applications, Chem. Biodivers. 15 (2018) e1700557.
[13]
S.D. Figueiro, A.A.M. Macedo, M.R.S. Melo, A.L.P. Freitas, R.A. Moreira, R.S. De Oliveira, J.C. Goes, A.S.B. Sombra, On the dielectric behaviour of collagen–algal sulfated polysaccharide blends: effect of glutaraldehyde crosslinking, Biophys. Chem. 120 (2006) 154–159.
Jo
ur
na
lP
re
-p
ro
of
[1]
14
Journal Pre-proof J.E. Eastoe, The amino acid composition of mammalian collagen and gelatin, Biochem. J. 61 (1955) 589.
[15]
A. Bonincontro, C. Mari, M. Mengoni, G. Risuleo, A study of the dielectric properties of E. coli ribosomal RNA and proteins in solution, Biophys. Chem. 67 (1997) 43–50.
[16]
J.A.J. Van Der Rijt, K.O. Van Der Werf, M.L. Bennink, P.J. Dijkstra, J. Feijen, Micromechanical testing of individual collagen fibrils, Macromol. Biosci. 6 (2006) 697–702.
[17]
J.C. Góes, S.D. Figueiró, J.A.C. De Paiva, I.F. De Vasconcelos, A.S.B. Sombra, On the piezoelectricity of anionic collagen films, J. Phys. Chem. Solids. 63 (2002) 465–470.
[18]
M. Beauvais, T. Degabriel, N. Aissaoui, V. Dupres, E. Colaço, K. El Kirat, R.F. Domingos, D. Brouri, C. Pradier, J. Spadavecchia, Supramolecular Self‐Assembly and Organization of Collagen at Solid/Liquid Interface: Effect of Spheroid‐and Rod‐Shaped TiO2 Nanocrystals, Adv. Mater. Interfaces. (2019) 1900195.
[19]
C.-Y. Hsieh, J.-C. Hwang, T.-H. Chang, J.-Y. Li, S.-H. Chen, L.-K. Mao, L.-S. Tsai, Y.-L. Chueh, P.-C. Lyu, S.S.H. Hsu, Enhanced mobility of organic thin film transistors by water absorption of collagen hydrolysate gate dielectric, Appl. Phys. Lett. 103 (2013) 118_1.
[20]
F. Bibi, M. Villain, C. Guillaume, B. Sorli, N. Gontard, A review: Origins of the dielectric properties of proteins and potential development as bio-sensors, Sensors. 16 (2016) 1232.
[21]
B.E. García-Moreno, J.J. Dwyer, A.G. Gittis, E.E. Lattman, D.S. Spencer, W.E. Stites, Experimental measurement of the effective dielectric in the hydrophobic core of a protein, Biophys. Chem. 64 (1997) 211–224.
[22]
W. Li, B.J. Dobraszczyk, J.D. Schofield, Stress relaxation behavior of wheat dough, gluten, and gluten protein fractions, Cereal Chem. 80 (2003) 333–338.
[23]
D. Van Belle, I. Couplet, M. Prevost, S.J. Wodak, D. Van Belle, I. Couplet, M. Prevost, S.J. Wodak, Calculations of electrostatic properties in proteins: Analysis of contributions from induced protein dipoles, J. Mol. Biol. 198 (1987) 721–735. doi:https://doi.org/10.1016/00222836(87)90213-0.
[24]
A. Warshel, S.T. Russell, Calculations of electrostatic interactions in biological systems and in solutions, Q. Rev. Biophys. 17 (1984) 283–422.
[25]
G. Bhatt, K. Mishra, G. Ramanathan, S. Bhattacharya, Dielectrophoresis assisted impedance spectroscopy for detection of gold-conjugated amplified DNA samples, Sensors Actuators B Chem. 288 (2019) 442–453.
[26]
X. Chen, Theoretical electrical conductivity of hydrogen-bonded benzamide-derived molecules and single DNA bases, J. Biol. Phys. 39 (2013) 607–624.
[27]
R. Baghbani, M.H. Moradi, M.B. Shadmehr, The Development of a Four-Electrode BioImpedance Sensor for Identification and Localization of Deep Pulmonary Nodules, Ann. Biomed. Eng. (2018) 1–12.
[28]
R. Baghbani, M.H. Moradi, M.B. Shadmehr, Z.M. Sanat, A New Bio-Impedance Forceps Sensor for Measuring Electrical Conductivity of the Biological Tissues, IEEE Sens. J. (2019).
[29]
B.S. Terrones, Broadband electrical impedance spectroscopy for dynamic electrical bio-impedance characterization, (2012).
[30]
S. Lee, B. Koo, J. Shin, E. Lee, H. Park, H. Kim, Effects of hydroxyl groups in polymeric
Jo
ur
na
lP
re
-p
ro
of
[14]
15
Journal Pre-proof dielectrics on organic transistor performance, Appl. Phys. Lett. 88 (2006) 162109. A. Sihvola, Mixing rules with complex dielectric coefficients, Subsurf. Sens. Technol. Appl. 1 (2000) 393–415.
[32]
H. Fischer, I. Polikarpov, A.F. Craievich, Average protein density is a molecular‐weight‐ dependent function, Protein Sci. 13 (2004) 2825–2828.
[33]
T. Freeborn, B. Fu, Fatigue-induced cole electrical impedance model changes of biceps tissue bioimpedance, Fractal Fract. 2 (2018) 27.
[34]
J. Yun, H.W. Kim, H.-I. Kim, J.-H. Lee, Electrical impedance spectroscopy on a needle for safer Veress needle insertion during laparoscopic surgery, Sensors Actuators B Chem. 250 (2017) 453– 460.
[35]
T.J. Freeborn, A survey of fractional-order circuit models for biology and biomedicine, IEEE J. Emerg. Sel. Top. Circuits Syst. 3 (2013) 416–424.
[36]
C.I. Trainito, O. Français, B. Le Pioufle, Analysis of pulsed electric field effects on cellular tissue with Cole–Cole model: Monitoring permeabilization under inhomogeneous electrical field with bioimpedance parameter variations, Innov. Food Sci. Emerg. Technol. 29 (2015) 193–200.
na
CRediT author statement
lP
re
-p
ro
of
[31]
Masoomeh Ashoorirad: Investigation, Conceptualization, Methodology, Writing-Original Draft Preparation, Visualization, Software, Validation, Formal Analysis.
ur
Mehrdad Saviz: Supervision, Writing-Reviewing and Editing, Formal Analysis, Conceptualization.
Jo
Ali Fallah: Supervision, Writing-Reviewing and Editing, Project Administration, Funding Acquisition.
Highlights: Measuring electrical properties of collagen solution with different concentrations Relations between electrical properties and the mass of collagen in the solution Water role on the electrical properties of the collagen solution Fitting an electrical model on to the measured data of the collagen solution 16
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8