Targeted multicolor in vivo imaging over 1,000 nm enabled by nonamethine cyanines.

Recent progress has shown that using wavelengths between 1,000 and 2,000 nm, referred to as the shortwave-infrared or near-infrared (NIR)-II range, can enable high-resolution in vivo imaging at depths not possible with conventional optical wavelengths. However, few bioconjugatable probes of the type that have proven invaluable for multiplexed imaging in the visible and NIR range are available for imaging these wavelengths. Using rational design, we have generated persulfonated indocyanine dyes with absorbance maxima at 872 and 1,072 nm through catechol-ring and aryl-ring fusion, respectively, onto the nonamethine scaffold. Multiplexed two-color and three-color in vivo imaging using monoclonal antibody and dextran conjugates in several tumor models illustrate the benefits of concurrent labeling of the tumor and healthy surrounding tissue and lymphatics. These efforts are enabled by complementary advances in a custom-built NIR/shortwave-infrared imaging setup and software package for multicolor real-time imaging.

Red-Shifted Coumarin Luciferins for Improved Bioluminescence Imaging.

Multicomponent bioluminescence imaging in vivo requires an expanded collection of tissue-penetrant probes. Toward this end, we generated a new class of near-infrared (NIR) emitting coumarin luciferin analogues (CouLuc-3s). The scaffolds were easily accessed from commercially available dyes. Complementary mutant luciferases for the CouLuc-3 analogues were also identified. The brightest probes enabled sensitive imaging in vivo. The CouLuc-3 scaffolds are also orthogonal to popular bioluminescent reporters and can be used for multicomponent imaging applications. Collectively, this work showcases a new set of bioluminescent tools that can be readily implemented for multiplexed imaging in a variety of biological settings.

Covalent Inhibition by a Natural Product-Inspired Latent Electrophile.

Strategies to target specific protein cysteines are critical to covalent probe and drug discovery. 3-Bromo-4,5-dihydroisoxazole (BDHI) is a natural product-inspired, synthetically accessible electrophilic moiety that has previously been shown to react with nucleophilic cysteines in the active site of purified enzymes. Here, we define the global cysteine reactivity and selectivity of a set of BDHI-functionalized chemical fragments using competitive chemoproteomic profiling methods. Our study demonstrates that BDHIs capably engage reactive cysteine residues in the human proteome and the selectivity landscape of cysteines liganded by BDHI is distinct from that of haloacetamide electrophiles. Given its tempered reactivity, BDHIs showed restricted, selective engagement with proteins driven by interactions between a tunable binding element and the complementary protein sites. We validate that BDHI forms covalent conjugates with glutathione S-transferase Pi (GSTP1) and peptidyl-prolyl cis-trans isomerase NIMA-interacting 1 (PIN1), emerging anticancer targets. BDHI electrophile was further exploited in Bruton's tyrosine kinase (BTK) inhibitor design using a single-step late-stage installation of the warhead onto acrylamide-containing compounds. Together, this study expands the spectrum of optimizable chemical tools for covalent ligand discovery and highlights the utility of 3-bromo-4,5-dihydroisoxazole as a cysteine-reactive electrophile.

Expedient Synthesis and Characterization of π-Extended Luciferins.

Bioluminescence imaging enables the sensitive tracking of cell populations and the visualization of biological processes in living systems. Bioluminescent luciferase/luciferin pairs with far-red and near-infrared emission benefit from the reduced competitive absorption by blood and tissue while also facilitating multiplexing strategies. Luciferins with extended π-systems, such as AkaLumine and recently reported CouLuc-1 and -3, can be used for bioluminescence imaging in this long wavelength regime. Existing synthetic routes to AkaLumine and similar π-extended compounds require a multistep sequence to install the thiazoline heterocycle. Here we detail the development of a two-step strategy for accessing these molecules via a Horner-Wadsworth-Emmons reaction and cysteine condensation sequence from readily available aldehyde starting materials. We detail an improved synthesis of AkaLumine, as well as the corresponding two-carbon homologues, Tri- and Tetra-AkaLumine. We then extended this approach to prepare coumarin- and naphthalene-derived luciferins. These putative luciferins were tested against a panel of luciferases to identify capable emitters. Of these, an easily prepared naphthalene derivative exhibits photon emission on par with that of the broadly used Akaluc/AkaLumine pair with similar emission maxima. Overall, this chemistry provides efficient access to several bioluminescent probes for a variety of imaging applications.

Protonation-Dependent Sequencing of 5-Formylcytidine in RNA.

Chemical modification of cytidine in noncoding RNAs plays a key role in regulating translation and disease. However, the distribution and dynamics of many of these modifications remain unknown due to a lack of sensitive site-specific sequencing technologies. Here, we report a protonation-dependent sequencing reaction for the detection of 5-formylcytidine (5fC) and 5-carboxycytidine (5caC) in RNA. First, we evaluate how protonation combined with electron-withdrawing substituents alters the molecular orbital energies and reduction of modified cytidine nucleosides, highlighting 5fC and 5caC as reactive species. Next, we apply this reaction to detect these modifications in synthetic oligonucleotides as well as endogenous human transfer RNA (tRNA). Finally, we demonstrate the utility of our method to characterize a patient-derived model of 5fC deficiency, where it enables facile monitoring of both pathogenic loss and exogenous rescue of NSUN3-dependent 5fC within the wobble base of human mitochondrial tRNAMet. These studies showcase the ability of protonation to enhance the reactivity and sensitive detection of 5fC in RNA and more broadly provide a molecular foundation for using optimized sequencing reactions to better understand the role of oxidized RNA cytidine residues in diseases.

Repurposing the Pummerer Rearrangement: Determination of Methionine Sulfoxides in Peptides.

The reversible oxidation of methionine residues in proteins has emerged as a biologically important post-translational modification. However, detection and quantitation of methionine sulfoxide in proteins is difficult. Our aim is to develop a method for specifically derivatizing methionine sulfoxide residues. We report a Pummerer rearrangement of methionine sulfoxide treated sequentially with trimethylsilyl chloride and then 2-mercaptoimidazole or pyridine-2-thiol to produce a dithioacetal product. This derivative is stable to standard mass spectrometry conditions, and its formation identified oxidized methionine residues. The scope and requirements of dithioacetal formation are reported for methionine sulfoxide and model substrates. The reaction intermediates have been investigated by computational techniques and by 13 C NMR spectroscopy. These provide evidence for an α-chlorinated intermediate. The derivatization allows for detection and quantitation of methionine sulfoxide in proteins by mass spectrometry and potentially by immunochemical methods.

Impact of Cyanine Conformational Restraint in the Near-Infrared Range.

Appending conformationally restraining ring systems to the cyanine chromophore creates exceptionally bright fluorophores in the visible range. Here, we report the application of this strategy in the near-infrared range through the preparation of the first restrained heptamethine indocyanine. Time-resolved absorption spectroscopy and fluorescence correlation spectroscopy verify that, unlike the corresponding parent unrestrained variant, the restrained molecule is not subject to photoisomerization. Notably, however, the room-temperature emission efficiency and the fluorescence lifetime of the restrained cyanine are not extended relative to the parent cyanine, even in viscous solvents. Thus, in contrast to prior reports, the photoisomerization of heptamethine cyanines does not contribute significantly to the excited-state chemistry of these molecules. We also find that the fluorescence lifetime of the restrained heptamethine cyanine is temperature-insensitive and significantly extended at moderately elevated temperatures relative to the parent cyanine. Finally, computational studies have been used to evaluate the impact of the conformational restraint on atomic and orbital structure across the cyanine series. These studies clarify the role of photoisomerization in the heptamethine cyanine scaffold and demonstrate the dramatic effect of restraint on the temperature sensitivity of these dyes.

Defining the conditional basis of silicon phthalocyanine near-IR ligand exchange.

Bond cleavage reactions initiated by long-wavelength light are needed to extend the scope of the caged-uncaged paradigm into complex physiological settings. Axially unsymmetrical silicon phthalocyanines (SiPcs) undergo efficient release of phenol ligands in a reaction contingent on three factors - near-IR light (690 nm), hypoxia, and a thiol reductant. These studies detail efforts to define the mechanistic basis for this unique conditionally-dependent bond cleavage reaction. Spectroscopic studies provide evidence for the formation of a key phthalocyanine radical anion intermediate formed from the triplet state in a reductant-dependent manner. Computational chemistry studies indicate that phenol ligand solvolysis proceeds through a heptacoordinate silicon transition state and that this solvolytic process is favored following SiPc radical anion formation. These results provide insight regarding the central role that radical anion intermediates formed through photoinduced electron transfer with biological reductants can play in long-wavelength uncaging reactions.

Hybrid Correlation Energy (HyCE): An Approach Based on Separate Evaluations of Internal and External Components.

A novel hybrid correlation energy (HyCE) approach is proposed that determines the total correlation energy via distinct computation of its internal and external components. This approach evolved from two related studies. First, rigorous assessment of the accuracies and size extensivities of a number of electron correlation methods, that include perturbation theory (PT2), coupled-cluster (CC), configuration interaction (CI), and coupled electron pair approximation (CEPA), shows that the CEPA(0) variant of the latter and triples-corrected CC methods consistently perform very similarly. These findings were obtained by comparison to near full CI results for four small molecules and by charting recovered correlation energies for six steadily growing chain systems. Second, by generating valence virtual orbitals (VVOs) and utilizing the CEPA(0) method, we were able to partition total correlation energies into internal (or nondynamic) and external (or dynamic) parts for the aforementioned six chain systems and a benchmark test bed of 36 molecules. When using triple-ζ basis sets it was found that per orbital internal correlation energies were appreciably larger than per orbital external energies and that the former showed far more chemical variation than the latter. Additionally, accumulations of external correlation energies were seen to proceed smoothly, and somewhat linearly, as the virtual space is gradually increased. Combination of these two studies led to development of the HyCE approach, whereby the internal and external correlation energies are determined separately by CEPA(0)/VVO and PT2/external calculations, respectively. When applied to the six chain systems and the 36-molecule benchmark test set it was found that HyCE energies followed closely those of triples-corrected CC and CEPA(0) while easily outperforming MP2 and CCSD. The success of the HyCE approach is more notable when considering that its cost is only slightly more than MP2 and significantly cheaper than the CC approaches.

Infrared Spectroscopy and Structure of (NO)(n) Clusters.

Nitrogen oxide clusters (NO)n have been studied in He droplets via infrared depletion spectroscopy and by quantum chemical calculations. The ν1 and ν5 bands of cis-ON-NO dimer have been observed at 1868.2 and 1786.5 cm(-1), respectively. Furthermore, spectral bands of the trimer and tetramer have been located in the vicinity of the corresponding dimer bands in accord with computed frequencies that place NO-stretch bands of dimer, trimer, and tetramer within a few wavenumbers of each other. In addition, a new line at 1878.1 cm(-1) close to the band origin of single molecules was assigned to van der Waals bound dimers of (NO)2, which are stabilized due to the rapid cooling in He droplets. Spectra of larger clusters (n > 5), have broad unresolved features in the vicinity of the dimer bands. Experiments and calculations indicate that trimers consist of a dimer and a loosely bound third molecule, whereas the tetramer consists of two weakly bound dimers.

The transition from the open minimum to the ring minimum on the ground state and on the lowest excited state of like symmetry in ozone: A configuration interaction study.

The metastable ring structure of the ozone 1(1)A1 ground state, which theoretical calculations have shown to exist, has so far eluded experimental detection. An accurate prediction for the energy difference between this isomer and the lower open structure is therefore of interest, as is a prediction for the isomerization barrier between them, which results from interactions between the lowest two (1)A1 states. In the present work, valence correlated energies of the 1(1)A1 state and the 2(1)A1 state were calculated at the 1(1)A1 open minimum, the 1(1)A1 ring minimum, the transition state between these two minima, the minimum of the 2(1)A1 state, and the conical intersection between the two states. The geometries were determined at the full-valence multi-configuration self-consistent-field level. Configuration interaction (CI) expansions up to quadruple excitations were calculated with triple-zeta atomic basis sets. The CI expansions based on eight different reference configuration spaces were explored. To obtain some of the quadruple excitation energies, the method of Correlation Energy Extrapolation by Intrinsic Scaling was generalized to the simultaneous extrapolation for two states. This extrapolation method was shown to be very accurate. On the other hand, none of the CI expansions were found to have converged to millihartree (mh) accuracy at the quadruple excitation level. The data suggest that convergence to mh accuracy is probably attained at the sextuple excitation level. On the 1(1)A1 state, the present calculations yield the estimates of (ring minimum-open minimum) ∼45-50 mh and (transition state-open minimum) ∼85-90 mh. For the (2(1)A1-(1)A1) excitation energy, the estimate of ∼130-170 mh is found at the open minimum and 270-310 mh at the ring minimum. At the transition state, the difference (2(1)A1-(1)A1) is found to be between 1 and 10 mh. The geometry of the transition state on the 1(1)A1 surface and that of the minimum on the 2(1)A1 surface nearly coincide. More accurate predictions of the energy differences also require CI expansions to at least sextuple excitations with respect to the valence space. For every wave function considered, the omission of the correlations of the 2s oxygen orbitals, which is a widely used approximation, was found to cause errors of about ±10 mh with respect to the energy differences.

Guanine holes are prominent targets for mutation in cancer and inherited disease.

Single base substitutions constitute the most frequent type of human gene mutation and are a leading cause of cancer and inherited disease. These alterations occur non-randomly in DNA, being strongly influenced by the local nucleotide sequence context. However, the molecular mechanisms underlying such sequence context-dependent mutagenesis are not fully understood. Using bioinformatics, computational and molecular modeling analyses, we have determined the frequencies of mutation at G • C bp in the context of all 64 5'-NGNN-3' motifs that contain the mutation at the second position. Twenty-four datasets were employed, comprising >530,000 somatic single base substitutions from 21 cancer genomes, >77,000 germline single-base substitutions causing or associated with human inherited disease and 16.7 million benign germline single-nucleotide variants. In several cancer types, the number of mutated motifs correlated both with the free energies of base stacking and the energies required for abstracting an electron from the target guanines (ionization potentials). Similar correlations were also evident for the pathological missense and nonsense germline mutations, but only when the target guanines were located on the non-transcribed DNA strand. Likewise, pathogenic splicing mutations predominantly affected positions in which a purine was located on the non-transcribed DNA strand. Novel candidate driver mutations and tissue-specific mutational patterns were also identified in the cancer datasets. We conclude that electron transfer reactions within the DNA molecule contribute to sequence context-dependent mutagenesis, involving both somatic driver and passenger mutations in cancer, as well as germline alterations causing or associated with inherited disease.

The somatic autosomal mutation matrix in cancer genomes.

DNA damage in somatic cells originates from both environmental and endogenous sources, giving rise to mutations through multiple mechanisms. When these mutations affect the function of critical genes, cancer may ensue. Although identifying genomic subsets of mutated genes may inform therapeutic options, a systematic survey of tumor mutational spectra is required to improve our understanding of the underlying mechanisms of mutagenesis involved in cancer etiology. Recent studies have presented genome-wide sets of somatic mutations as a 96-element vector, a procedure that only captures the immediate neighbors of the mutated nucleotide. Herein, we present a 32 × 12 mutation matrix that captures the nucleotide pattern two nucleotides upstream and downstream of the mutation. A somatic autosomal mutation matrix (SAMM) was constructed from tumor-specific mutations derived from each of 909 individual cancer genomes harboring a total of 10,681,843 single-base substitutions. In addition, mechanistic template mutation matrices (MTMMs) representing oxidative DNA damage, ultraviolet-induced DNA damage, (5m)CpG deamination, and APOBEC-mediated cytosine mutation, are presented. MTMMs were mapped to the individual tumor SAMMs to determine the maximum contribution of each mutational mechanism to the overall mutation pattern. A Manhattan distance across all SAMM elements between any two tumor genomes was used to determine their relative distance. Employing this metric, 89.5% of all tumor genomes were found to have a nearest neighbor from the same tissue of origin. When a distance-dependent 6-nearest neighbor classifier was used, 10.4% of the SAMMs had an Undetermined tissue of origin, and 92.2% of the remaining SAMMs were assigned to the correct tissue of origin. [corrected]. Thus, although tumors from different tissues may have similar mutation patterns, their SAMMs often display signatures that are characteristic of specific tissues.

Aminolysis of an N-diazeniumdiolated amidine as an approach to diazeniumdiolated ammonia.

Recent theoretical studies have suggested that the parent diazeniumdiolate ion, H2N-N(O)═NO(-) ("diazeniumdiolated ammonia"), might be stable enough to be isolated and that it could potentially serve as a uniquely advantageous prodrug form of bioactive nitroxyl (HNO). Here, we report on an attempt to isolate its O(2)-benzylated derivative by aminolysis of the C═N bond in PhC(NH2)═N-N(O)═NOBn. The reaction proved remarkably sluggish in comparison to aminolysis of unsubstituted benzamidine, and the desired product could not be isolated, apparently because of base sensitivity of the NH2 group. Consistent with this interpretation, O-benzylhydroxylamine and N2O were recovered from the reaction mixture in high yields, along with N,N'-dibutylbenzamidine. Theoretical calculations rationalize the observed slow aminolysis by demonstrating that the diazeniumdiolate group greatly suppresses the electrophilicity of the adjacent C═N carbon center, rendering attack at that position endothermic. The data provide significant insights into the challenges inherent to the pursuit of diazeniumdiolated ammonia.

Covalent bonds are created by the drive of electron waves to lower their kinetic energy through expansion.

An analysis based on the variation principle shows that in the molecules H2 (+), H2, B2, C2, N2, O2, F2, covalent bonding is driven by the attenuation of the kinetic energy that results from the delocalization of the electronic wave function. For molecular geometries around the equilibrium distance, two features of the wave function contribute to this delocalization: (i) Superposition of atomic orbitals extends the electronic wave function from one atom to two or more atoms; (ii) intra-atomic contraction of the atomic orbitals further increases the inter-atomic delocalization. The inter-atomic kinetic energy lowering that (perhaps counter-intuitively) is a consequence of the intra-atomic contractions drives these contractions (which per se would increase the energy). Since the contractions necessarily encompass both, the intra-atomic kinetic and potential energy changes (which add to a positive total), the fact that the intra-atomic potential energy change renders the total potential binding energy negative does not alter the fact that it is the kinetic delocalization energy that drives the bond formation.

Decoding nitric oxide release rates of amine-based diazeniumdiolates.

Amine-based diazeniumdiolates (NONOates) have garnered widespread use as nitric oxide (NO) donors, and their potential for nitroxyl (HNO) release has more recently been realized. While NO release rates can vary significantly with the type of amine, half-lives of seconds to days under physiological conditions, there is as yet no way to determine a priori the NO or HNO production rates of a given species, and no discernible trends have manifested other than that secondary amines produce only NO (i.e., no HNO). As a step to understanding these complex systems, here we describe a procedure for modeling amine-based NONOates in water solvent that provides an excellent correlation (R(2) = 0.94) between experimentally measured dissociation rates of seven secondary amine species and their computed NO release activation energies. The significant difference in behavior of NONOates in the gas and solvent phases is also rigorously demonstrated via explicit additions of quantum mechanical water molecules. The presented results suggest that the as-yet unsynthesized simplest amine-based NONOate, the diazeniumdiolated ammonia anion [H2N-N(O)═NO(-)], could serve as an unperturbed HNO donor. These results provide a step forward toward the accurate modeling of general NO and/or HNO donors as well as for the identification of tailored prodrug candidates.

Categorizing biases in high-confidence high-throughput protein-protein interaction data sets.

We characterized and evaluated the functional attributes of three yeast high-confidence protein-protein interaction data sets derived from affinity purification/mass spectrometry, protein-fragment complementation assay, and yeast two-hybrid experiments. The interacting proteins retrieved from these data sets formed distinct, partially overlapping sets with different protein-protein interaction characteristics. These differences were primarily a function of the deployed experimental technologies used to recover these interactions. This affected the total coverage of interactions and was especially evident in the recovery of interactions among different functional classes of proteins. We found that the interaction data obtained by the yeast two-hybrid method was the least biased toward any particular functional characterization. In contrast, interacting proteins in the affinity purification/mass spectrometry and protein-fragment complementation assay data sets were over- and under-represented among distinct and different functional categories. We delineated how these differences affected protein complex organization in the network of interactions, in particular for strongly interacting complexes (e.g. RNA and protein synthesis) versus weak and transient interacting complexes (e.g. protein transport). We quantified methodological differences in detecting protein interactions from larger protein complexes, in the correlation of protein abundance among interacting proteins, and in their connectivity of essential proteins. In the latter case, we showed that minimizing inherent methodology biases removed many of the ambiguous conclusions about protein essentiality and protein connectivity. We used these findings to rationalize how biological insights obtained by analyzing data sets originating from different sources sometimes do not agree or may even contradict each other. An important corollary of this work was that discrepancies in biological insights did not necessarily imply that one detection methodology was better or worse, but rather that, to a large extent, the insights reflected the methodological biases themselves. Consequently, interpreting the protein interaction data within their experimental or cellular context provided the best avenue for overcoming biases and inferring biological knowledge.

High-level theoretical study of the NO dimer and tetramer: has the tetramer been observed?

The ground-state properties of (NO)(2) and (NO)(4) have been investigated using multireference second-order perturbation theory (MRMP2) and include a two-tier extrapolation to the complete basis set (CBS) limit. For the NO dimer the MRMP2(18,14)/CBS predicted structure, binding energy (with respect to 2NO; D(e) = 3.46 kcal/mol), and dipole moment (u(e) = 0.169 D) are in excellent agreement with experimental measurements (D(e) = 2.8-3.8 kcal/mol; u(e) = 0.171 D). Additionally, three of four intermolecular anharmonic MRMP2(18,14)/CBS-estimated frequencies (143 cm(-1), 238 cm(-1), 421 cm(-1)) are in excellent agreement with recent gas-phase experimental measurements (135 cm(-1), 239 cm(-1), 429/428 cm(-1)); however, the predicted value of 151 cm(-1) for the out-of-plane torsion (A(2)) is elevated compared to recent experimental estimates of 97-117 cm(-1). Our finding that this infrared-forbidden vibration is also predicted to have an extremely low Raman activity (0.04 Å/amu at the MP2/QZ level of theory) conflicts with Raman measurements of a strong intensity for a low frequency band; however, these studies were performed for low temperature solid and liquid phases. Probing the possibility of the presence of higher order clusters we investigated the stability of (NO)(4) and discovered three isomers, each resembling pairs of dimers, that were stable to dissociation to 2(NO)(2), with the lowest-energy isomer (C(i) structure) having a predicted binding energy almost identical to that of the dimer. Computed vibrational frequencies of the C(i) isomer indicate that the 12 highest-frequency modes resemble barely shifted NO dimer-combined bands while the 13th highest-frequency mode at ~100 cm(-1) is exclusive to (NO)(4). Moreover, this tetramer-unique vibration is infrared inactive but has a very high predicted Raman activity of some 24 Å/amu. Guided by the theoretical results, we reexamined and reassigned experimental Raman and infrared data going back to 1951 and determined that our best predictions of vibrational frequencies of (NO)(2) and (NO)(4) are consistent with experimental observations. We thus postulate the existence and observation of (NO)(4).

Ketone Incorporation Extends the Emission Properties of the Xanthene Scaffold Beyond 1000 nm.

Imaging in the shortwave-infrared region (SWIR, λ = 1000-2500 nm) has the potential to enable deep tissue imaging with high resolution. Critical to the development of these methods is the identification of low molecular weight, biologically compatible fluorescent probes that emit beyond 1000 nm. Exchanging the bridging oxygen atom on the xanthene scaffold (C10' position) with electron withdrawing groups has been shown to lead to significant redshifts in absorbance and emission. Guided by quantum chemistry computational modeling studies, we investigated the installation of a ketone bridge at the C10' position. This simple modification extends the absorbance maxima to 860 nm and the emission beyond 1000 nm, albeit with reduced photon output. Overall, these studies demonstrate that broadly applied xanthene dyes can be extended into the SWIR range.

Modified norcyanines enable ratiometric pH imaging beyond 1000 nm.

Activatable fluorophores with emission beyond 1000 nm have the potential to enable high contrast imaging in complex in vivo settings. However, there are few scaffolds that can be applied to this challenge. Here we detail the synthesis and evaluation of benzo[c,d]indole-substituted norcyanines that enable pH responsive fluorescence imaging in the long wavelength (>1150 nm) range. A key component of our molecular design is the installation of a hydrophilic substituted quaternary amine in the central dihydropyridine ring system. A compound with a C4'-phenyl substituent, but not the C4'-protio homologue, exhibits absorbance maxima of 740 nm and 1130 nm in basic and acidic media, respectively, with evidence of J-aggregate-like properties. These two distinct absorbances enabled ratiometric imaging of probe internalization in a tumor model. Overall, these studies provide a new class of activatable long-wavelength responsive fluorophores with promising photophysical properties.