Defining substrate interactions with calreticulin: an isothermal titration calorimetric study

Calreticulin (CRT) is a soluble, lectin chaperone found in the endoplasmic reticulum of eukaryotes. It binds the N-glycosylated polypeptides via the glycan intermediate Glc₁Man_5–9GlcNAc₂, present on the target glycoproteins. Earlier we have studied interactions of substrate with CRT by isothermal titration calorimetry (ITC) and molecular modeling, to establish that CRT recognizes the Glcα1–3 linkage and forms contacts with each saccharide moiety of the oligosaccharide Glcα1–3Manα1–2Manα1–2Man. We also delineated the amino acid residues in the sugar binding pocket of CRT that play a crucial role in sugar–CRT binding. Here, we have used mono-deoxy analogues of the trisaccharide unit Glcα1–3Manα1–2Man to determine the role of various hydroxyl groups of the sugar substrate in sugar–CRT interactions. Using the thermodynamic data obtained by ITC with these analogues we demonstrate that the 3-OH group of Glc1 plays an important role in sugar–CRT binding, whereas the 6-OH group does not. Also, the 4-OH, 6-OH of Man2 and 3-OH, 4-OH of Man3 in the trisaccharide are involved in binding, of which 6-OH of Man2 and 4-OH of Man3 have a more significant role to play. This study sheds light further on the interactions between the substrate sugar of glycoproteins and the lectin chaperone CRT.     

Thermodynamics of ion-induced RNA folding in the hammerhead ribozyme: an isothermal titration calorimetric study

The hammerhead ribozyme undergoes a well-defined two-stage conformational folding process, induced by the binding of magnesium ions. In this study, we have used isothermal titration calorimetry to analyze the thermodynamics of magnesium binding and magnesium ion-induced folding of the ribozyme. Binding to the natural sequence ribozyme is strongly exothermic and can be analyzed in terms of sequential interaction at two sites with association constants K(A) = 480 and 2840 M(-1). Sequence variants of the hammerhead RNA give very different isothermal titration curves. An A14G variant that cannot undergo ion-induced folding exhibits endothermic binding. By contrast, a deoxyribose G5 variant that can undergo only the first of the two folding transitions gives a complex titration curve. However, despite these differences the ITC data for all three species can be analyzed in terms of the sequential binding of magnesium ions at two sites. While the binding affinities are all in the region of 10(3) M(-1), corresponding to free energies of Delta G degrees = -3.5 to -4 kcal mol(-1), the enthalpic and entropic contributions show much greater variation. The ITC experiments are in good agreement with earlier conformational studies of the folding of the ion-induced folding of the hammerhead ribozyme.     

Insights into hydrophobicity and the chaperone-like function of alphaA- and alphaB-crystallins: an isothermal titration calorimetric study

Alpha-crystallin, composed of two subunits, alphaA and alphaB, has been shown to function as a molecular chaperone that prevents aggregation of other proteins under stress conditions. The exposed hydrophobic surfaces of alpha-crystallins have been implicated in this process, but their exact role has not been elucidated. In this study, we quantify the hydrophobic surfaces of alphaA- and alphaB-crystallins by isothermal titration calorimetry using 8-anilino-1-napthalenesulfonic acid (ANS) as a hydrophobic probe and analyze its correlation to the chaperone potential of alphaA- and alphaB-crystallins under various conditions. Two ANS binding sites, one with low and another with high affinity, were clearly detected, with alphaB showing a higher number of sites than alphaA at 30 degrees C. In agreement with the higher number of hydrophobic sites, alphaB-crystallin demonstrated higher chaperone activity than alphaA at this temperature. Thermodynamic analysis of ANS binding to alphaA- and alphaB-crystallins indicates that high affinity binding is driven by both enthalpy and entropy changes, with entropy dominating the low affinity binding. Interestingly, although the number of ANS binding sites was similar for alphaA and alphaB at 15 degrees C, alphaA was more potent than alphaB in preventing aggregation of the insulin B-chain. Although there was no change in the number of high affinity binding sites of alphaA and alphaB for ANS upon preheating, there was an increase in the number of low affinity sites of alphaA and alphaB. Preheated alphaA, in contrast to alphaB, exhibited remarkably enhanced chaperone activity. Our results indicate that although hydrophobicity appears to be a factor in determining the chaperone-like activity of alpha-crystallins, it does not quantitatively correlate with the chaperone function of alpha-crystallins.     

In-vitro assessment of the binding mechanism of oxyfluorfen (herbicide) with garlic phytocystatin: multi-spectroscopic and isothermal titration calorimetric study

Abstract

Oxyfluorfen (2-chloro-1-(3-ethoxy-4-nitrophenoxy)-4-(trifluoromethyl)benzene) is a nitrophenyl ether herbicide. Phytocystatins are crucial plant proteins which regulate various physiological processes and are also responsible for maintaining protease–antiprotease balance within plants. Thus, the present article deciphers the interaction of oxyfluorfen with garlic phytocystatin (GPC) through various spectroscopic and calorimetric techniques. The cysteine proteinase inhibitory assay was done to assess the inhibitory action of GPC in the presence of oxyfluorfen. The GPC loses its inhibitory activity in the presence of oxyfluorfen. The complex formation of GPC-oxyfluorfen was shown by UV absorption spectroscopy. The intrinsic fluorescence experiment affirmed the quenching of GPC in the presence of oxyfluorfen. The Stern–Volmer quenching constant and binding constant was obtained as 6.89?×?10³ M^?1 and 9.72?×?10³ M^?1, respectively. Synchronous fluorescence showed the alteration in the microenvironment around tyrosine residues. 3D fluorescence suggested the perturbation in the polarity around aromatic residues. The isothermal titration experiment suggests that the interaction of oxyfluorfen with GPC is a thermodynamically favorable reaction. Secondary structure alteration of GPC in the presence of oxyfluorfen was studied by circular dichroism (CD). The CD result showed a reduction in the α-helical content of GPC on interaction with oxyfluorfen. Consequently, all these outcomes affirmed the formation of GPC–oxyfluorfen complex along with the structural and conformational alteration. This study identifies and signifies that the exposure of oxyfluorfen induces stress within the plant system.

Communicated by Ramaswamy H. Sarma     

Energetics of carbohydrate binding to Momordica charantia (bitter gourd) lectin: an isothermal titration calorimetric study

Physico-chemical and carbohydrate binding studies have been carried out on the Momordica charantia (bitter gourd) seed lectin (MCL). The lectin activity is maximal in the pH range 7.4-11.0, but decreases steeply below pH 7.0. The lectin activity is mostly unaffected in the temperature range 4-50 degrees C, but a sharp decrease is seen between 50 and 60 degrees C, which could be correlated to changes in the structure of the protein as seen by circular dichroism and fluorescence spectroscopy. Isothermal titration calorimetric studies show that the tetrameric MCL binds two sugar molecules and the binding constants (Kb), determined at 288.15 K, for various saccharides were found to vary between 7.3 x 10(3) and 1.52 x 10(4)M(-1). The binding reactions for all the saccharides investigated were essentially enthalpy driven, with the binding enthalpies (DeltaHb) at 288.15 K being in the range of -50.99 and -43.39 kJ mol(-1), whereas the contribution to the binding reaction from the entropy of binding was negative, with values of binding entropy (DeltaSb) ranging between -99.2 and -72.0 J mol(-1)K(-1) at 288.15 K. Changes in heat capacity (DeltaCp) for the binding of disaccharides, lactose and lactulose, were significantly larger in magnitude than those obtained for the monosaccharides, methyl-beta-D-galactopyranoside, and methyl-alpha-D-galactopyranoside, and could be correlated reasonably well with the surface areas of these ligands. Enthalpy-entropy compensation was observed for all the sugars studied, suggesting that water structure plays an important role in the overall binding reaction. CD spectroscopy indicates that carbohydrate binding does not lead to significant changes in the secondary and tertiary structures of MCL, suggesting that the carbohydrate binding sites on this lectin are mostly preformed.     

The development of a continuous isothermal titration calorimetric method for equilibrium studies

A continuous isothermal titration calorimetry (cITC) method for microcalorimeters has been developed. The method is based on continuous slow injection of a titrant into the calorimetric vessel. The experimental time for a cITC binding experiment is 12-20 min and the number of data points obtained is on the order of 1000. This gives an advantage over classical isothermal titration calorimetry (ITC) binding experiments that need 60-180 min to generate 20-30 data points. The method was validated using two types of calorimeters, which differ in calorimetric principle, geometry, stirring, and way of delivering the titrant into the calorimetric vessel. Two different experimental systems were used to validate the method: the binding of Ba(2+) to 18-crown-6 and the binding of cytidine 2'-monophosphate to RNAse A. Both systems are used as standard test systems for titration calorimetry. Computer simulations show that the dynamic range for determination of equilibrium constants can be increased by three orders of magnitude compared to that of classical ITC, making it possible to determine high affinities. Simulations also show an improved possibility to elucidate the actual binding model from cITC data. The simulated data demonstrate that cITC makes it easier to discriminate between different thermodynamic binding models due to the higher density of data points obtained from one experiment.     

Entropy-driven folding of an RNA helical junction: an isothermal titration calorimetric analysis of the hammerhead ribozyme

Helical junctions are extremely common motifs in naturally occurring RNAs, but little is known about the thermodynamics that drive their folding. Studies of junction folding face several challenges: non-two-state folding behavior, superposition of secondary and tertiary structural energetics, and drastically opposing enthalpic and entropic contributions to folding. Here we describe a thermodynamic dissection of the folding of the hammerhead ribozyme, a three-way RNA helical junction, by using isothermal titration calorimetry of bimolecular RNA constructs. By using this method, we show that tertiary folding of the hammerhead core occurs with a highly unfavorable enthalpy change, and is therefore entropically driven. Furthermore, the enthalpies and heat capacities of core folding are the same whether supported by monovalent or divalent ions. These properties appear to be general to the core sequence of bimolecular hammerhead constructs. We present a model for the ion-induced folding of the hammerhead core that is similar to those advanced for the folding of much larger RNAs, involving ion-induced collapse to a structured, non-native state accompanied by rearrangement of core residues to produce the native fold. In agreement with previous enzymological and structural studies, our thermodynamic data suggest that the hammerhead structure is stabilized in vitro predominantly by diffusely bound ions. Our approach addresses several significant challenges that accompany the study of junction folding, and should prove useful in defining the thermodynamic determinants of stability in these important RNA motifs.     

An easy-to-use tool for planning and modeling a calorimetric titration

Isothermal titration calorimetry (ITC) is widely employed to measure thermodynamic properties of binding interactions between two macromolecules or a macromolecule and a small ligand. No labeling of interacting species is required for ITC, but this advantage is offset by potentially material-consuming experimental optimization complicated by an indirect readout of an ITC titration. Here we present a simple, practical, and portable spreadsheet-based tool for planning and modeling an ITC titration experiment accompanied by basic guidelines.     

Interaction of superoxide dismutase with phospholipid liposomes. An uptake, spin label and calorimetric study

The effect of superoxide dismutases from five species upon phospholipid bilayers has been investigated. The uptake by egg phosphatidylcholine bilayers of the holo and apo forms of bovine superoxide dismutase increases with enzyme concentration and only a fraction of each is removed by treatment with trypsin. These uptake data indicate that both forms of the enzyme associate with and are embedded within lipid bilayers. From the spectrum of the spin label 2-(3-carboxypropyl)-4,4-dimethyl-2-tridecyl-3-oxazolidinyloxyl, the binding of superoxide dismutase to egg phosphatidylcholine bilayers can be shown to disorder the lipid packing. The disordering by the bovine holoenzyme is small but increases with increasing enzyme concentration and period of incubation. The disordering effects of the apoenzyme are much larger and are reversible by Cu2+, Zn2+ reconstitution of the apoenzyme. The disordering effect of the apoenzyme is further confirmed by differential scanning calorimetry. The gel to liquid crystalline phase transition of egg phosphatidylcholine is lowered 7 degrees C by 25% by weight apo-superoxide dismutase to lipid. Human, dog, swordfish and yeast superoxide dismutases also disorder, and to a greater extent than the bovine enzyme. The greatest perturbation is produced by yeast superoxide dismutase; a 20% decrease in the order parameter by 50% by weight enzyme to lipid.     

Interaction of cationic dodecyl‐trimethyl‐ammonium bromide with oxy‐HbGp by isothermal titration and differential scanning calorimetric studies: Effect of proximity of isoelectric point

In this work, isothermal titration and differential scanning calorimetric methods, in combination with pyrene fluorescence emission and dynamic light scattering have been used to investigate the interaction of dodecyltrimethylammonium bromide (DTAB) with the giant extracellular Glossoscolex paulistus hemoglobin (HbGp) in the oxy‐form, at pH values around the isoelectric point (pI ≈ 5.5). Our ITC results have shown that the interaction of DTAB with the hemoglobin is more intense at pH 7.0, with a smaller cac (critical aggregation concentration) value. The increase of protein concentration does not influence the cac value of the interaction, at both pH values. Therefore, the beginning of the DTAB‐oxy‐HbGp premicellar aggregates formation, in the cac region, is not affected by the increase of protein concentration. HSDSC studies show higher T_m values at pH 5.0, in the absence and presence of DTAB, when compared with pH 7.0. Furthermore, at pH 7.0, an aggregation process is observed with DTAB in the range from 0.75 to 1.5 mmol/L, noticed by the exothermic peak, and similar to that observed for pure oxy‐HbGp, at pH 5.0, and in the presence of DTAB. DLS melting curves show a decrease on the hemoglobin thermal stability for the oxy‐HbGp‐DTAB mixtures and formation of larger aggregates, at pH 7.0. Our present data, together with previous results, support the observation that the protein structural changes, at pH 7.0, occur at smaller DTAB concentrations, as compared with pH 5.0, due to the acidic pI of protein that favors the oxy‐HbGp‐cationic surfactant interaction at neutral pH. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 199–211, 2016.     

Elucidating the binding thermodynamics of 8-anilino-1-naphthalene sulfonic acid with the A-state of alpha-lactalbumin: an isothermal titration calorimetric investigation

Isothermal titration calorimetry has been used to demonstrate that the heat profile associated with the binding of 8-anilino-1-naphthalene sulfonic acid (ANS) with the acid induced molten globule state (A-state) of alpha-lactalbumin (alpha-LA) is different from that with the native and denatured states of the protein. The results corroborate the spectroscopic observations that ANS binds more strongly to the partially folded states of the protein compared to that with the native and denatured states. ANS binds to the A-state of alpha-LA at two independent binding sites that remain nearly the same in the temperature range of 10-35 degrees C. The number of moles of ANS binding at site 1 at 10 degrees C is 14.0+/-0.2 and remains nearly the same with rise in temperature. However, the number of moles of ANS molecules binding at site 2 show an increase from 1.6+/-0.2 at 10 degrees C to 4.1+/-0.1 at 35 degrees C. The deviation of the slope of enthalpy-entropy compensation plot from unity and nonadherence to van't Hoff dictates implies that the binding sites on the A-state of alpha-LA for ANS are not well defined and specific; rather, these binding sites are formed due to greater exposure of hydrophobic clusters in the A-state of the protein. The results for the first time demonstrate the use of isothermal titration calorimetry in characterizing the A-state of alpha-LA both qualitatively and quantitatively.     

Synthesis of chiral chloroquine and its analogues as antimalarial agents

In this investigation, we describe a new approach to chiral synthesis of chloroquine and its analogues. All tested compounds displayed potent activity against chloroquine sensitive as well as chloroquine resistant strains of Plasmodium falciparum in vitro and Plasmodium yoelii in vivo. Compounds S-13b, S-13c, S-13d and S-13i displayed excellent in vitro antimalarial activity with an IC₅₀ value of 56.82, 60.41, 21.82 and 7.94 nM, respectively, in the case of resistant strain. Furthermore, compounds S-13a, S-13c and S-13d showed in vivo suppression of 100% parasitaemia on day 4 in the mouse model against Plasmodium yoelii when administered orally. These results underscore the application of synthetic methodology and need for further lead optimization.     

An isothermal titration calorimetric method to determine the kinetic parameters of enzyme catalytic reaction by employing the product inhibition as probe.

An isothermal titration calorimetric (ITC) method was developed to measure the kinetic parameters of ribonuclease A catalytic hydrolysis of cytidine 2',3'-cyclic monophosphate. Employing the inhibition of product as a probe, the K(m), K(i), k(c), and DeltaH(m) can be determined by two simple calorimetric measurements. First, the substrate was titrated into the cell containing high concentration of enzyme. The molar reaction heat was calculated from the titration peak area divided by substrate moles per titration, and the initial catalytic reaction rate in the presence of various concentrations of product can be calculated from the peak height and the molar reaction heat. From Michaelis-Menten function in the presence of inhibitors, the relationship between K(m) and K(i) can be obtained. Then, the dissociation constant, which is equal to K(i), was measured by a regular ITC experiment. Thus, K(m) and k(c) can be calculated. The method developed here can be applied in other enzyme catalytic systems with inhibitive products.     

Nucleotide binding to creatine kinase: an isothermal titration microcalorimetry study

We investigated the binding of ATP in the presence and absence of Mg(2+) to dimeric muscle creatine kinase (CK) by isothermal titration microcalorimetry as a function of pH and temperature. The thermodynamic parameters for these events show that (1) binding of nucleotide to the CK active site does not involve proton exchange with the buffer and (2) the active sites are the only nucleotide binding sites on CK. Interdependence of the active sites in the dimer could not be demonstrated. As CK undergoes major structural changes upon Mg-nucleotide binding, a thermodynamic cycle was employed to calculate the contributions of domain movements to the observed enthalpies.     

A study of statistical error in isothermal titration calorimetry

In isothermal titration calorimetry (ITC), the two main sources of random (statistical) error are associated with the extraction of the heat q from the measured temperature changes and with the delivery of metered volumes of titrant. The former leads to uncertainty that is approximately constant and the latter to uncertainty that is proportional to q. The role of these errors in the analysis of ITC data by nonlinear least squares is examined for the case of 1:1 binding, M+X right arrow over left arrow MX. The standard errors in the key parameters-the equilibrium constant Ko and the enthalpy DeltaHo-are assessed from the variance-covariance matrix computed for exactly fitting data. Monte Carlo calculations confirm that these "exact" estimates will normally suffice and show further that neglect of weights in the nonlinear fitting can result in significant loss of efficiency. The effects of the titrant volume error are strongly dependent on assumptions about the nature of this error: If it is random in the integral volume instead of the differential volume, correlated least-squares is required for proper analysis, and the parameter standard errors decrease with increasing number of titration steps rather than increase.     

Heat capacity changes associated with guanine quadruplex formation: an isothermal titration calorimetry study

This study addresses the temperature dependence of the enthalpy of formation for several unimolecular quadruplexes in the presence of excess monovalent salt. We examined a series of biologically significant guanine-rich DNA sequences: thrombin binding aptamer (TBA) (d(G(2)T(2)G(2)TGTG(2)T(2)G(2)), PS2.M, a catalytically active aptamer (d(GTG(3)TAG(3)CG(3)T(2)G(2))), and the human telomere repeat (HT) (d(AG(3)(T(2)AG(3))(3))). Using CD spectra and UV melting, we confirmed the presence of quadruplex structures and established the temperature range in which quadruplex conformation is stable. We then performed ITC experiments, adding DNA to a solution containing excess NaCl or KCl. In this approach, only several additions are made, and only the enthalpy of quadruplex formation is measured. This measurement was repeated at different temperatures to determine the temperature dependence of the enthalpy change accompanying quadruplex formation. To control for the effect of nonspecific salt interactions during DNA folding, we repeated the experiment by replacing the quadruplex-forming sequences with a similar but nonfolding sequence. Dilution enthalpies were also subtracted to obtain the final enthalpy value involving only the quadruplex folding process. For all sequences studied, quadruplex formation was exothermic but with an increasing magnitude with increasing temperature. These results are discussed in terms of the change in heat capacity associated with quadruplex formation.     

Ferrous ion binding to recombinant human H-chain ferritin. An isothermal titration calorimetry study

Iron deposition within the iron storage protein ferritin involves a complex series of events consisting of Fe(2+) binding, transport, and oxidation at ferroxidase sites and mineralization of a hydrous ferric oxide core, the storage form of iron. In the present study, we have examined the thermodynamic properties of Fe(2+) binding to recombinant human H-chain apoferritin (HuHF) by isothermal titration calorimetry (ITC) in order to determine the location of the primary ferrous ion binding sites on the protein and the principal pathways by which the Fe(2+) travels to the dinuclear ferroxidase center prior to its oxidation to Fe(3+). Calorimetric titrations show that the ferroxidase center is the principal locus for Fe(2+) binding with weaker binding sites elsewhere on the protein and that one site of the ferroxidase center, likely the His65 containing A-site, preferentially binds Fe(2+). That only one site of the ferroxidase center is occupied by Fe(2+) implies that Fe(2+) oxidation to form diFe(III) species might occur in a stepwise fashion. In dilute anaerobic protein solution (3-5 microM), only 12 Fe(2+)/protein bind at pH 6.51 increasing to 24 Fe(2+)/protein at pH 7.04 and 7.5. Mutation of ferroxidase center residues (E62K+H65G) eliminates the binding of Fe(2+) to the center, a result confirming the importance of one or both Glu62 and His65 residues in Fe(2+) binding. The total Fe(2+) binding capacity of the protein is reduced in the 3-fold hydrophilic channel variant S14 (D131I+E134F), indicating that the primary avenue by which Fe(2+) gains access to the interior of ferritin is through these eight channels. The binding stoichiometry of the channel variant is one-third that of the recombinant wild-type H-chain ferritin whereas the enthalpy and association constant for Fe(2+) binding are similar for the two with an average values (DeltaH degrees = 7.82 kJ/mol, binding constant K = 1.48 x 10(5) M(-)(1) at pH 7.04). Since channel mutations do not completely prevent Fe(2+) binding to the ferroxidase center, iron gains access to the center in approximately one-third of the channel variant molecules by other pathways.     

Interaction of gymnemic acid with cyclodextrins analyzed by isothermal titration calorimetry, NMR and dynamic light scattering

The physiological phenomenon that the antisweet taste effect of gymnemic acid (GA) is diminished by application of gamma-cyclodextrin (gamma-CD) to the mouth was evaluated at the molecular level using isothermal titration calorimetry, NMR and dynamic light scattering. These analyses showed that GA specifically binds to gamma-CD. Thermodynamic analysis using isothermal titration calorimetry revealed that the association constant of GA and gamma-CD is 10(5)-10(6) m(-1) with favorable enthalpy and entropy changes. The heat capacity change was negative and large, despite the change in accessible surface area upon binding being small. These thermodynamics indicate that the binding is dominated by hydrophobic interactions, which is in agreement with inclusion complex formation of gamma-CD. In addition, NMR measurements showed that in solution the spectra of GA are broad and sharpened by the addition of gamma-CD, indicating that unbound GA is in a water-soluble aggregate that is dispersed when it forms a complex with gamma-CD. Dynamic light scattering showed that the average diameter of unbound GA is > 30 nm and that of GA and gamma-CD complex is 2.2 nm, similar to unbound gamma-CD, supporting the aggregate property of GA and the inclusion complexation of GA by gamma-CD.     

Interaction of a designed interleukin-10 epitope mimic with an antibody studied by isothermal titration microcalorimetry

The mechanism of recognition of proteins and peptides by antibodies and the factors determining binding affinity and specificity are mediated by essentially the same features. However, additional effects of the usually unfolded and flexible solution structure of peptide ligands have to be considered. In an earlier study we designed and optimized six peptides (pepI to pepVI) mimicking the discontinuous binding site of interleukin-10 for the anti-interleukin-10 monoclonal antibody (mab) CB/RS/1. Three of them were selected for analysis of their solution conformation by circular dichroism measurements. The peptides differ in the content of alpha-helices and in the inducibility of helical secondary structures by trifluoroethanol. These properties, however, do not correlate with the binding affinity. PepVI, a 32-mer cyclic epitope mimic, has the highest affinity to mab CB/RS/1 identified to date. CD difference spectroscopy suggests an increase of the alpha-helix content of pepVI with complex formation. Binding of pepVI to mab CB/RS/1 is characterized by a large negative, favorable binding enthalpy and a smaller unfavorable loss of entropy (DeltaH degrees = -16.4 kcal x mol(-1), TDeltaS degrees = -6.9 kcal x mol(-1)) resulting in DeltaG degrees = -9.5 kcal x mol(-1) at 25 degrees C as determined by isothermal titration calorimetry. Binding of pepVI is enthalpically driven over the entire temperature range studied (10-35 degrees C). Complex formation is not accompanied by proton uptake or release. A negative heat capacity change DeltaC(p) of -0.354 kcal x mol(-1) x K(-1) was determined from the temperature dependence of DeltaH degrees. The selection of protein mimics with the observed thermodynamic properties is promoted by the applied identification and iterative optimization procedure.