Hydropathic analysis of the free energy differences in anthracycline antibiotic binding to DNA

Molecular models of six anthracycline antibiotics and their complexes with 32 distinct DNA octamer sequences were created and analyzed using HINT (Hydropathic INTeractions) to describe binding. The averaged binding scores were then used to calculate the free energies of binding for comparison with experimentally determined values. In parsing our results based on specific functional groups of doxorubicin, our calculations predict a free energy contribution of –3.6 ± 1.1 kcal mol^–1 (experimental –2.5 ± 0.5 kcal mol^–1) from the groove binding daunosamine sugar. The net energetic contribution of removing the hydroxyl at position C9 is –0.7 ± 0.7 kcal mol^–1 (–1.1 ± 0.5 kcal mol^–1). The energetic contribution of the 3′ amino group in the daunosamine sugar (when replaced with a hydroxyl group) is –3.7 ± 1.1 kcal mol^–1 (–0.7 ± 0.5 kcal mol^–1). We propose that this large discrepancy may be due to uncertainty in the exact protonation state of the amine. The energetic contribution of the hydroxyl group at C14 is +0.4 ± 0.6 kcal mol^–1 (–0.9 ± 0.5 kcal mol^–1), largely due to unfavorable hydrophobic interactions between the hydroxyl oxygen and the methylene groups of the phosphate backbone of the DNA. Also, there appears to be considerable conformational uncertainty in this region. This computational procedure calibrates our methodology for future analyses where experimental data are unavailable.     

Energetics of echinomycin binding to DNA

Differential scanning calorimetry and UV thermal denaturation have been used to determine a complete thermodynamic profile for the bis-intercalative interaction of the peptide antibiotic echinomycin with DNA. The new calorimetric data are consistent with all previously published binding data, and afford the most rigorous and direct determination of the binding enthalpy possible. For the association of echinomycin with DNA, we found ΔG° = –7.6 kcal mol^–1, ΔH = +3.8 kcal mol^–1 and ΔS = +38.9 cal mol^–1 K^–1 at 20°C. The binding reaction is clearly entropically driven, a hallmark of a process that is predominantly stabilized by hydrophobic interactions, though a deeper analysis of the free energy contributions suggests that direct molecular recognition between echinomycin and DNA, mediated by hydrogen bonding and van der Waals contacts, also plays an important role in stabilizing the complex.     

Markov state models elucidate the stability of DNA influenced by the chiral 5S-Tg base

The static and dynamic structures of DNA duplexes affected by 5S-Tg (Tg, Thymine glycol) epimers were studied using MD simulations and Markov State Models (MSMs) analysis. The results show that the 5S,6S-Tg base caused little perturbation to the helix, and the base-flipping barrier was determined to be 4.4 kcal mol⁻¹ through the use of enhanced sampling meta-eABF calculations, comparable to 5.4 kcal mol⁻¹ of the corresponding thymine flipping. Two conformations with the different hydrogen bond structures between 5S,6R-Tg and A19 were identified in several independent MD trajectories. The 5S,6R-Tg:O6H_O6•••N1:A19 hydrogen bond is present in the high-energy conformation displaying a clear helical distortion, and near barrier-free Tg base flipping. The low-energy conformation always maintains Watson–Crick base pairing between 5S,6R-Tg and A19, and 5S-Tg base flipping is accompanied by a small barrier of ca. 2.0 K_BT (T = 298 K). The same conformations are observed in the MSMs analysis. Moreover, the transition path and metastable structures of the damaged base flipping are for the first time verified through MSMs analysis. The data clearly show that the epimers have completely different influence on the stability of the DNA duplex, thus implying different enzymatic mechanisms for DNA repair.     

Determination of base and backbone contributions to the thermodynamics of premelting and melting transitions in B DNA

In previous papers of this series the temperature-dependent Raman spectra of poly(dA)·poly(dT) and poly(dA–dT)·poly(dA–dT) were used to characterize structurally the melting and premelting transitions in DNAs containing consecutive A·T and alternating A·T/T·A base pairs. Here, we describe procedures for obtaining thermodynamic parameters from the Raman data. The method exploits base-specific and backbone-specific Raman markers to determine separate thermodynamic contributions of A, T and deoxyribosyl-phosphate moieties to premelting and melting transitions. Key findings include the following: (i) Both poly(dA)·poly(dT) and poly(dA–dT)· poly(dA–dT) exhibit robust premelting transitions, due predominantly to backbone conformational changes. (ii) The significant van’t Hoff premelting enthalpies of poly(dA)·poly(dT) [ΔH_vH^pm = 18.0 ± 1.6 kcal·mol^–1 (kilocalories per mole cooperative unit)] and poly(dA–dT)·poly(dA–dT) (ΔH_vH^pm = 13.4 ± 2.5 kcal·mol^–1) differ by an amount (~4.6 kcal·mol^–1) estimated as the contribution from three-centered inter-base hydrogen bonding in (dA)_n·(dT)_n tracts. (iii) The overall stacking free energy of poly(dA)· poly(dT) [–6.88 kcal·mol_bp^–1 (kilocalories per mole base pair)] is greater than that of poly(dA–dT)· poly(dA–dT) (–6.31 kcal·mol_bp^–1). (iv) The difference between stacking free energies of A and T is significant in poly(dA)·poly(dT) (ΔΔG_st = 0.8 ± 0.3 kcal· mol_bp^–1), but marginal in poly(dA–dT)·poly(dA–dT) (ΔΔG_st = 0.3 ± 0.3 kcal·mol_bp^–1). (v) In poly(dA)· poly(dT), the van’t Hoff parameters for melting of A (ΔH_vH^A = 407 ± 23 kcal·mol^–1, ΔS_vH^A = 1166 ± 67 cal·°K^–1·mol^–1, ΔG_vH(25°C)^A = 60.0 ± 3.2 kcal·mol^–1) are clearly distinguished from those of T (ΔH_vH^T = 185 ± 38 kcal·mol^–1, ΔS_vH^T = 516 ± 109 cal·°K^–1·mol^–1, ΔG_vH(25°C)^T = 27.1 ± 5.5 kcal·mol^–1). (vi) Similar relative differences are observed in poly(dA–dT)· poly(dA–dT) (ΔH_vH^A = 333 ± 54 kcal·mol^–1, ΔS_vH^A = 961 ± 157 cal·°K^–1·mol^–1, ΔG_vH(25°C)^A = 45.0 ± 7.6 kcal· mol^–1; ΔH_vH^T = 213 ± 30 kcal·mol^–1, ΔS_vH^T = 617 ± 86 cal·°K^–1·mol^–1, ΔG_vH(25°C)^T = 29.3 ± 4.9 kcal·mol^–1). The methodology employed here distinguishes thermodynamic contributions of base stacking, base pairing and backbone conformational ordering in the molecular mechanism of double-helical B DNA formation.     

Kinetics of formation of hypoxanthine containing base pairs by HIV-RT: RNA template effects on the base substitution frequencies

Hypoxanthine (H), the deamination product of adenine, has been implicated in the high frequency of A to G transitions observed in retroviral and other RNA genomes. Although H·C base pairs are thermodynamically more stable than other H·N pairs, polymerase selection may be determined in part by kinetic factors. Therefore, the hypoxanthine induced substitution pattern resulting from replication by viral polymerases may be more complex than that predicted from thermodynamics. We have examined the steady-state kinetics of formation of base pairs opposite template H in RNA by HIV-RT, and for the incorporation of dITP during first- and second-strand synthesis. Hypoxanthine in an RNA template enhances the k_2app for pairing with standard dNTPs by factors of 10–1000 relative to adenine at the same sequence position. The order of base pairing preferences for H in RNA was observed to be H·C >> H·T > H·A > H·G. Steady-state kinetics of insertion for all possible mispairs formed with dITP were examined on RNA and DNA templates of identical sequence. Insertion of dITP opposite all bases occurs 2–20 times more frequently on RNA templates. This bias for higher insertion frequencies on RNA relative to DNA templates is also observed for formation of mispairs at template A. This kinetic advantage afforded by RNA templates for mismatches and pairing involving H suggests a higher induction of mutations at adenines during first-strand synthesis by HIV-RT.     

Recombination R-triplex: H-bonds contribution to stability as revealed with minor base substitutions for adenine

Several cellular processes involve alignment of three nucleic acids strands, in which the third strand (DNA or RNA) is identical and in a parallel orientation to one of the DNA duplex strands. Earlier, using 2-aminopurine as a fluorescent reporter base, we demonstrated that a self-folding oligonucleotide forms a recombination-like structure consistent with the R-triplex. Here, we extended this approach, placing the reporter 2-aminopurine either in the 5′- or 3′-strand. We obtained direct evidence that the 3′-strand forms a stable duplex with the complementary central strand, while the 5′-strand participates in non-Watson–Crick interactions. Substituting 2,6-diaminopurine or 7-deazaadenine for adenine, we tested and confirmed the proposed hydrogen bonding scheme of the A*(T·A) R-type triplet. The adenine substitutions expected to provide additional H-bonds led to triplex structures with increased stability, whereas the substitutions consistent with a decrease in the number of H-bonds destabilized the triplex. The triplex formation enthalpies and free energies exhibited linear dependences on the number of H-bonds predicted from the A*(T·A) triplet scheme. The enthalpy of the 10 nt long intramolecular triplex of −100 kJ·mol⁻¹ demonstrates that the R-triplex is relatively unstable and thus an ideal candidate for a transient intermediate in homologous recombination, t-loop formation at the mammalian telomere ends, and short RNA invasion into a duplex. On the other hand, the impact of a single H-bond, 18 kJ·mol⁻¹, is high compared with the overall triplex formation enthalpy. The observed energy advantage of a ‘correct’ base in the third strand opposite the Watson–Crick base pair may be a powerful mechanism for securing selectivity of recognition between the single strand and the duplex.     

Derivatization of DNAs with selenium at 6-position of guanine for function and crystal structure studies

To investigate nucleic acid base pairing and stacking via atom-specific mutagenesis and crystallography, we have synthesized for the first time the 6-Se-deoxyguanosine phosphoramidite and incorporated it into DNAs via solid-phase synthesis with a coupling yield over 97%. We found that the UV absorption of the Se-DNAs red-shifts over 100 nm to 360 nm (ε = 2.3 × 10⁴ M⁻¹ cm⁻¹), the Se-DNAs are yellow colored, and this Se modification is relatively stable in water and at elevated temperature. Moreover, we successfully crystallized a ternary complex of the Se-G-DNA, RNA and RNase H. The crystal structure determination and analysis reveal that the overall structures of the native and Se-modified nucleic acid duplexes are very similar, the selenium atom participates in a Se-mediated hydrogen bond (Se … H–N), and the ^SeG and C form a base pair similar to the natural G–C pair though the Se-modification causes the base-pair to shift (approximately 0.3 Å). Our biophysical and structural studies provide new insights into the nucleic acid flexibility, duplex recognition and stability. Furthermore, this novel selenium modification of nucleic acids can be used to investigate chemogenetics and structure of nucleic acids and their protein complexes.     

Determination of thermodynamic parameters for HIV DIS type loop-loop kissing complexes

The HIV-1 type dimerization initiation signal (DIS) loop was used as a starting point for the analysis of the stability of Watson–Crick (WC) base pairs in a tertiary structure context. We used ultraviolet melting to determine thermodynamic parameters for loop–loop tertiary interactions and compared them with regular secondary structure RNA helices of the same sequences. In 1 M Na⁺ the loop–loop interaction of a HIV-1 DIS type pairing is 4 kcal/mol more stable than its sequence in an equivalent regular and isolated RNA helix. This difference is constant and sequence independent, suggesting that the rules governing the stability of WC base pairs in the secondary structure context are also valid for WC base pairs in the tertiary structure context. Moreover, the effect of ion concentration on the stability of loop–loop tertiary interactions differs considerably from that of regular RNA helices. The stabilization by Na⁺ and Mg²⁺ is significantly greater if the base pairing occurs within the context of a loop–loop interaction. The dependence of the structural stability on salt concentration was defined via the slope of a T_m/log [ion] plot. The short base-paired helices are stabilized by 8°C/log [Mg²⁺] or 11°C/log [Na⁺], whereas base-paired helices forming tertiary loop–loop interactions are stabilized by 16°C/log [Mg²⁺] and 26°C/log [Na⁺]. The different dependence on ionic strength that is observed might reflect the contribution of specific divalent ion binding to the preformation of the hairpin loops poised for the tertiary kissing loop–loop contacts.     

ClaRNA: a classifier of contacts in RNA 3D structures based on a comparative analysis of various classification schemes

The understanding of folding and function of RNA molecules depends on the identification and classification of interactions between ribonucleotide residues. We developed a new method named ClaRNA for computational classification of contacts in RNA 3D structures. Unique features of the program are the ability to identify imperfect contacts and to process coarse-grained models. Each doublet of spatially close ribonucleotide residues in a query structure is compared to clusters of reference doublets obtained by analysis of a large number of experimentally determined RNA structures, and assigned a score that describes its similarity to one or more known types of contacts, including pairing, stacking, base–phosphate and base–ribose interactions. The accuracy of ClaRNA is 0.997 for canonical base pairs, 0.983 for non-canonical pairs and 0.961 for stacking interactions. The generalized squared correlation coefficient (GC²) for ClaRNA is 0.969 for canonical base pairs, 0.638 for non-canonical pairs and 0.824 for stacking interactions. The classifier can be easily extended to include new types of spatial relationships between pairs or larger assemblies of nucleotide residues. ClaRNA is freely available via a web server that includes an extensive set of tools for processing and visualizing structural information about RNA molecules.     

Amino acid–base interactions: a three-dimensional analysis of protein–DNA interactions at an atomic level

To assess whether there are universal rules that govern amino acid–base recognition, we investigate hydrogen bonds, van der Waals contacts and water-mediated bonds in 129 protein–DNA complex structures. DNA–backbone interactions are the most numerous, providing stability rather than specificity. For base interactions, there are significant base–amino acid type correlations, which can be rationalised by considering the stereochemistry of protein side chains and the base edges exposed in the DNA structure. Nearly two-thirds of the direct read-out of DNA sequences involves complex networks of hydrogen bonds, which enhance specificity. Two-thirds of all protein–DNA interactions comprise van der Waals contacts, compared to about one-sixth each of hydrogen and water-mediated bonds. This highlights the central importance of these contacts for complex formation, which have previously been relegated to a secondary role. Although common, water-mediated bonds are usually non-specific, acting as space-fillers at the protein–DNA interface. In conclusion, the majority of amino acid–base interactions observed follow general principles that apply across all protein–DNA complexes, although there are individual exceptions. Therefore, we distinguish between interactions whose specificities are ‘universal’ and ‘context-dependent’. An interactive Web-based atlas of side chain–base contacts provides access to the collected data, including analyses and visualisation of the three-dimensional geometry of the interactions.     

Stacking Interactions between Carbohydrate and Protein Quantified by Combination of Theoretical and Experimental Methods

Carbohydrate – receptor interactions are an integral part of biological events. They play an important role in many cellular processes, such as cell-cell adhesion, cell differentiation and in-cell signaling. Carbohydrates can interact with a receptor by using several types of intermolecular interactions. One of the most important is the interaction of a carbohydrate''s apolar part with aromatic amino acid residues, known as dispersion interaction or CH/π interaction. In the study presented here, we attempted for the first time to quantify how the CH/π interaction contributes to a more general carbohydrate - protein interaction. We used a combined experimental approach, creating single and double point mutants with high level computational methods, and applied both to Ralstonia solanacearum (RSL) lectin complexes with α-l-Me-fucoside. Experimentally measured binding affinities were compared with computed carbohydrate-aromatic amino acid residue interaction energies. Experimental binding affinities for the RSL wild type, phenylalanine and alanine mutants were −8.5, −7.1 and −4.1 kcal.mol⁻¹, respectively. These affinities agree with the computed dispersion interaction energy between carbohydrate and aromatic amino acid residues for RSL wild type and phenylalanine, with values −8.8, −7.9 kcal.mol⁻¹, excluding the alanine mutant where the interaction energy was −0.9 kcal.mol⁻¹. Molecular dynamics simulations show that discrepancy can be caused by creation of a new hydrogen bond between the α-l-Me-fucoside and RSL. Observed results suggest that in this and similar cases the carbohydrate-receptor interaction can be driven mainly by a dispersion interaction.     

A calorimetric investigation of single stranded base stacking in the ribo-oligonucleotide A7

Differential scanning calorimetry has been employed to determine the energy change associated with single stranded base stacking in the ribo-oligonucleotide A₇. A total enthaipy change of 20.3 kcal (mole of heptamer)^?1 was measured. This corresponds to 2.9 kcal (mole of adenine)^?1 or 3.4 kcal (mole of A-A stack)^?1 if one assumes that all six stacking interactions are energetically equivalent. These results represent the first direct determination of this important parameter for a ribo-oligonucleotide. It is noted that the calorimetrically determined value reported here is considerably lower than any of the previously published van 't Hoff enthalpies but is consistent with values that can be derived from other calorimetric data.     

Mg2+-induced triplex formation of an equimolar mixture of poly(rA) and poly(rU)

Magnesium ions strongly influence the structure and biochemical activity of RNA. The interaction of Mg²⁺ with an equimolar mixture of poly(rA) and poly(rU) has been investigated by UV spectroscopy, isothermal titration calorimetry, ultrasound velocimetry and densimetry. Measurements in dilute aqueous solutions at 20°C revealed two differ ent processes: (i) Mg²⁺ binding to unfolded poly(rA)·poly(rU) up to [Mg²⁺]/[phosphate] = 0.25; and (ii) poly(rA)·2poly(rU) triplex formation at [Mg²⁺]/[phosphate] between 0.25 and 0.5. The enthalpies of these two different processes are favorable and similar to each other, ~–1.6 kcal mol^–1 of base pairs. Volume and compressibility effects of the first process are positive, 8 cm³ mol^–1 and 24 × 10^–4 cm³ mol^–1 bar^–1, respectively, and correspond to the release of water molecules from the hydration shells of Mg²⁺ and the polynucleotides. The triplex formation is also accompanied by a positive change in compressibility, 14 × 10^–4 cm³ mol^–1 bar^–1, but only a small change in volume, 1 cm³ mol^–1. A phase diagram has been constructed from the melting experiments of poly(rA)·poly(rU) at a constant K⁺ concentration, 140 mM, and various amounts of Mg²⁺. Three discrete regions were observed, corresponding to single-, double- and triple-stranded complexes. The phase boundary corresponding to the transition between double and triple helical conformations lies near physiological salt concentrations and temperature.     

Incorporation of a cationic aminopropyl chain in DNA hairpins: thermodynamics and hydration

We report on the physicochemical effects resulting from incorporating a 5-(3-aminopropyl) side chain onto a 2′-deoxyuridine (dU) residue in a short DNA hairpin. A combination of spectroscopy, calorimetry, density and ultrasound techniques were used to investigate both the helix–coil transition of a set of hairpins with the following sequence: d(GCGACTTTTTGNCGC) [N = dU, deoxythymidine (dT) or 5-(3-aminopropyl)-2′-deoxyuridine (dU*)], and the interaction of each hairpin with Mg²⁺. All three molecules undergo two-state transitions with melting temperatures (T_M) independent of strand concentration that indicates their intramolecular hairpin formation. The unfolding of each hairpin takes place with similar T_M values of 64–66°C and similar thermodynamic profiles. The unfavorable unfolding free energies of 6.4–6.9 kcal/mol result from the typical compensation of unfavorable enthalpies, 36–39 kcal/mol, and favorable entropies of ~110 cal/mol. Furthermore, the stability of each hairpin increases as the salt concentration increases, the T_M-dependence on salt yielded slopes of 2.3–2.9°C, which correspond to counterion releases of 0.53 (dU and dT) and 0.44 (dU*) moles of Na⁺ per mole of hairpin. Absolute volumetric and compressibility measurements reveal that all three hairpins have similar hydration levels. The electrostatic interaction of Mg²⁺ with each hairpin yielded binding affinities in the order: dU > dT > dU*, and a similar release of 2–4 electrostricted water molecules. The main result is that the incorporation of the cationic 3-aminopropyl side chain in the major groove of the hairpin stem neutralizes some local negative charges yielding a hairpin molecule with lower charge density.     

Thermodynamic criteria for high hit rate antisense oligonucleotide design

Antisense oligonucleotides are used for therapeutic applications and in functional genomic studies. In practice, however, many of the oligonucleotides complementary to an mRNA have little or no antisense activity. Theoretical strategies to improve the ‘hit rate’ in antisense screens will reduce the cost of discovery and may lead to identification of antisense oligonucleotides with increased potency. Statistical analysis performed on data collected from more than 1000 experiments with phosphorothioate-modified oligonucleotides revealed that the oligo-probes, which form stable duplexes with RNA (ΔG^o₃₇ ≤ –30 kcal/mol) and have small self-interaction potential, are more frequently efficient than molecules that form less stable oligonucleotide–RNA hybrids or more stable self-structures. To achieve optimal statistical preference, the values for self-interaction should be (ΔG^o₃₇) ≥ –8 kcal/mol for inter-oligonucleotide pairing and (ΔG^o₃₇) ≥ –1.1 kcal/mol for intra-molecular pairing. Selection of oligonucleotides with these thermodynamic values in the analyzed experiments would have increased the ‘hit rate’ by as much as 6-fold.     

Water-membrane partition thermodynamics of an amphiphilic lipopeptide: an enthalpy-driven hydrophobic effect

To shed light on the driving force for the hydrophobic effect that partitions amphiphilic lipoproteins between water and membrane, we carried out an atomically detailed thermodynamic analysis of a triply lipid modified H-ras heptapeptide anchor (ANCH) in water and in a DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine) bilayer. Combining molecular mechanical and continuum solvent approaches with an improved technique for solute entropy calculation, we obtained an overall transfer free energy of ~−13 kcal mol⁻¹. This value is in qualitative agreement with free energy changes derived from a potential of mean force calculation and indirect experimental observations. Changes in free energies of solvation and ANCH conformational reorganization are unfavorable, whereas ANCH-DMPC interactions—especially van der Waals—favor insertion. These results are consistent with an enthalpy-driven hydrophobic effect, in accord with earlier calorimetric data on the membrane partition of other amphiphiles. Furthermore, structural and entropic analysis of molecular dynamics-generated ensembles suggests that conformational selection may play a hitherto unappreciated role in membrane insertion of lipid-modified peptides and proteins.     

Enkephalin analogs containing 4,4-difluoro-2-aminobutyric acid: Synthesis and fluorine effect on the biological activity

Analogs of Met-enkephalin and [d -Pen², d -Pen⁵]enkephalin (DPDPE) containing the partially fluorinated amino acid 4,4-difluoro-2-aminobutyric acid (DFAB) in the 2- or 3-position of the peptide sequence were synthesized and their opioid activities and receptor selectivities were determined in vitro. The linear fluorinated [d -DFAB², Met⁵-NH₂]enkephalin showed μ and δ agonist potencies comparable to those of natural [Leu⁵]enkephalin. The partially fluorinated DPDPE analogs behaved differently as compared with their non-fluorinated correlates. While l -amino acid substitution in position 3 of DPDPE usually resulted in higher δ agonist potency than d -amino acid substitution, [d -DFAB³]DPDPE turned out to be a more potent δ agonist than [l -DFAB³]DPDPE. Furthermore, [d -DFAB³]DPDPE showed over 100-fold higher δ agonist potency than [d -Abu³]DPDPE (Abu=2-aminobutyric acid), indicating that the fluorine substituents interact favorably with a δ opioid receptor subsite. © 1998 European Peptide Society and John Wiley & Sons, Ltd.     

Identification of NH...N hydrogen bonds by magic angle spinning solid state NMR in a double-stranded RNA associated with myotonic dystrophy

RNA plays a central role in biological processes and exhibits a variety of secondary and tertiary structural features that are often stabilized via hydrogen bonds. The distance between the donor and acceptor nitrogen nuclei involved in NH…N hydrogen bonds in nucleic acid base pairs is typically in the range of 2.6–2.9 Å. Here, we show for the first time that such spatial proximity between ¹⁵N nitrogen nuclei can be conveniently monitored via magic angle spinning solid state NMR on a uniformly ¹⁵N-labelled RNA. The presence of NH…N hydrogen bonds is reflected as cross-peaks between the donor and acceptor nitrogen nuclei in 2D ¹⁵N dipolar chemical shift correlation spectra. The RNA selected for this experimental study was a CUG repeat expansion implicated in the neuromuscular disease myotonic dystrophy. The results presented provide direct evidence that the CUG repeat expansion adopts a double-stranded conformation.     

Water Transport in Intact Yeast Cells as Assessed by Fluorescence Self-Quenching

Intact yeast cells loaded with 5- and-6-carboxyfluorescein were used to assess water transport. The results were similar to those previously reported for protoplasts assessed by using either fluorescence or light scattering, and the activation energies were 8.0 and 15.1 kcal mol⁻¹ (33.4 and 63.2 kJ mol⁻¹) for a strain overexpressing AQY1 aquaporin and a parental strain, respectively.