Binding of cationic and neutral phenanthridine intercalators to a DNA oligomer is controlled by dispersion energy: quantum chemical calculations and molecular mechanics simulations

In silico and in vitro assessments of the mutagenicity of the azilsartan photoproduct

Photodegradation of azilsartan yields a phenanthridine derivative (APP). We suspected that APP could be a DNA-reactive substance, since many phenanthridine derivatives are mutagenic. In silico quantitative structure-activity relationship analysis indicated potential mutagenicity of APP, due to DNA reactivity at the 6-aminophenanthridine moiety. However, APP was not mutagenic in the Ames test. Density functional theory (DFT) calculations showed that APP cannot intercalate into DNA, due to its nonplanar structure, resulting from steric hindrance of its phenanthridine and benzimidazole moieties.

The Effective Synthesis of New Benzoquinoline Derivatives as Small Molecules with Anticancer Activity

In this study, some novel benzo[c]quinoline derivatives were synthesized, their structural characteristics were described, and their in vitro anticancer efficacy was investigated. The synthesis involves an initial quaternization of the nitrogen atom from benzo[c]quinoline and then a [3+2] dipolar cycloaddition reaction of the in situ formed ylide. The effectiveness of synthesis using traditional thermal heating (TH) compared to microwave (MW) and ultrasound (US) irradiation was investigated in detail. The setup of a reaction under MW or US irradiation offers a number of additional benefits: higher yields, a reduction in the amount of solvent used compared to TH, a reduction in the reaction time from hours to minutes, and a reduction in the amount of energy consumed. The structure of all the obtained compounds was proved by several spectral techniques (FTIR, HRMS, and NMR). All benzo[c]quinoline derivatives (quaternary salts and cycloadducts) along with ten other benzo[f]quinoline derivatives (quaternary salts and cycloadducts), previously obtained, were tested in an in vitro single-dose anticancer experiment. The results demonstrated that the cycloadducts 5a-c and 6a-c exhibit stronger anticancer activity than quaternary salts 3a-c. The most active compound is compound 5a, with anticancer activity on most of the cell lines studied, while the second most active compound is 6c, showing significant lethality for the SR leukemia cell line (17%). Structure-activity relationship (SAR) correlations are also included in the study.

Molecular photoswitches in aqueous environments

Molecular photoswitches enable dynamic control of processes with high spatiotemporal precision, using light as external stimulus, and hence are ideal tools for different research areas spanning from chemical biology to smart materials. Photoswitches are typically organic molecules that feature extended aromatic systems to make them responsive to (visible) light. However, this renders them inherently lipophilic, while water-solubility is of crucial importance to apply photoswitchable organic molecules in biological systems, like in the rapidly emerging field of photopharmacology. Several strategies for solubilizing organic molecules in water are known, but there are not yet clear rules for applying them to photoswitchable molecules. Importantly, rendering photoswitches water-soluble has a serious impact on both their photophysical and biological properties, which must be taken into consideration when designing new systems. Altogether, these aspects pose considerable challenges for successfully applying molecular photoswitches in aqueous systems, and in particular in biologically relevant media. In this review, we focus on fully water-soluble photoswitches, such as those used in biological environments, in both in vitro and in vivo studies. We discuss the design principles and prospects for water-soluble photoswitches to inspire and enable their future applications.

Ethidium bromide interactions with DNA: an exploration of a classic DNA-ligand complex with unbiased molecular dynamics simulations

Visualization of double stranded DNA in gels with the binding of the fluorescent dye ethidium bromide has been a basic experimental technique in any molecular biology laboratory for >40 years. The interaction between ethidium and double stranded DNA has been observed to be an intercalation between base pairs with strong experimental evidence. This presents a unique opportunity for computational chemistry and biomolecular simulation techniques to benchmark and assess their models in order to see if the theory can reproduce experiments and ultimately provide new insights. We present molecular dynamics simulations of the interaction of ethidium with two different double stranded DNA models. The first model system is the classic sequence d(CGCGAATTCGCG)₂ also known as the Drew–Dickerson dodecamer. We found that the ethidium ligand binds mainly stacked on, or intercalated between, the terminal base pairs of the DNA with little to no interaction with the inner base pairs. As the intercalation at the terminal CpG steps is relatively rapid, the resultant DNA unwinding, rigidification, and increased stability of the internal base pair steps inhibits further intercalation. In order to reduce these interactions and to provide a larger groove space, a second 18-mer DNA duplex system with the sequence d(GCATGAACGAACGAACGC) was tested. We computed molecular dynamics simulations for 20 independent replicas with this sequence, each with ∼27 μs of sampling time. Results show several spontaneous intercalation and base-pair eversion events that are consistent with experimental observations. The present work suggests that extended MD simulations with modern DNA force fields and optimized simulation codes are allowing the ability to reproduce unbiased intercalation events that we were not able to previously reach due to limits in computing power and the lack of extensively tested force fields and analysis tools.

Importance of model size in quantum mechanical studies of DNA intercalation

The convergence of DFT-computed interaction energies with increasing binding site model size was assessed. The data show that while accurate intercalator interaction energies can be derived from binding site models featuring only the flanking nucleotides for uncharged intercalators that bind parallel to the DNA base pairs, errors remain significant even when including distant nucleotides for intercalators that are charged, exhibit groove-binding tails that engage in noncovalent interactions with distant nucleotides, or that bind perpendicular to the DNA base pairs. Consequently, binding site models that include at least three adjacent nucleotides are required to consistently predict converged binding energies. The computationally inexpensive HF-3c method is shown to provide reliable interaction energies and can be routinely applied to such large models.

Spectroscopic analysis of tylosin adsorption on extracellular DNA reveals its interaction mechanism

Extracellular DNA (eDNA), which is commonly detected in aquatic and terrestrial environments, may be involved in gene transfer, increases in genetic diversity, and evolution. However, it has been reported that some small organic molecules or heavy metal ions can influence the transformation of DNA and even destroy its structure. We previously found that tylosin (TYL, a kind of antibiotic) is adsorbed onto salmon sperm DNA in a mixed solution. However, it is not clear whether this antibiotic affects the structure of DNA, and the mechanism of their interaction needs to be clarified. Therefore, we investigated the adsorption of TYL on different concentrations of salmon sperm DNA using agarose gel electrophoresis, ultraviolet-visible (UV-vis) spectroscopy, fluorescence spectroscopy, and surface enhanced Raman spectroscopy (SERS) to elucidate the interaction mechanism between TYL and DNA. The results showed that the adsorption of TYL decreased with increased concentrations of DNA. The electrophoresis band of pristine DNA was at 5000 bps. The brightness of the DNA band decreased with the TYL concentration and their incubation time. As the concentration of TYL increased, the fluorescence absorption intensity of DNA decreased significantly. Redshift and hyperchromicity were observed in the UV-vis adsorption spectrum with the presence of TYL in DNA solution, and they weakened as the DNA concentration increased. The Raman spectrum intensities of characteristic peaks in the mixed solution were weaker than that of pure TYL solution, and the peak intensity increased with increasing DNA concentration. Even a part of TYL characteristic peaks disappeared in the mixed solution. These results indicated that the pyran and macrolide of TYL might intercalate into the base pair plane of DNA. In addition, electrostatic attraction between TYL and DNA and interactions among TYL molecules may also play a role in the interaction mechanism. However, the double helix structure of DNA was not subject to the interaction of TYL.

Anticancer activity expressed by a library of 2,9-diazaperopyrenium dications

Polyaromatic compounds are well-known to intercalate DNA. Numerous anticancer chemotherapeutics have been developed upon the basis of this recognition motif. The compounds have been designed such that they interfere with the role of the topoisomerases, which control the topology of DNA during the cell-division cycle. Although many promising chemotherapeutics have been developed upon the basis of polyaromatic DNA intercalating systems, these candidates did not proceed past clinical trials on account of their dose-limiting toxicity. Herein, we discuss an alternative, water-soluble class of polyaromatic compounds, the 2,9-diazaperopyrenium dications, and report in vitro cell studies for a library of these dications. These investigations reveal that a number of 2,9-diazaperopyrenium dications show similar activities as doxorubicin toward a variety of cancer cell lines. Additionally, we report the solid-state structures of these dications, and we relate their tendency to aggregate in solution to their toxicity profiles. The addition of bulky substituents to these polyaromatic dications decreases their tendency to aggregate in solution. The derivative substituted with 2,6-diisopropylphenyl groups proved to be the most cytotoxic against the majority of the cell lines tested. In the solid state, the 2,6-diisopropylphenyl-functionalized derivative does not undergo π···π stacking, while in aqueous solution, dynamic light scattering reveals that this derivative forms very small (50-100 nm) aggregates, in contrast with the larger ones formed by dications with less bulky substituents. Alteration of the aromaticitiy in the terminal heterocycles of selected dications reveals a drastic change in the toxicity of these polyaromatic species toward specific cell lines.

Sequence-dependent binding of flavonoids to duplex DNA

The whole family of structurally distinct flavonoids has been recognized as a valuable source of prospective anticancer agents. There is experimental evidence demonstrating that some flavonoids, like flavopiridol (FLP) and quercetin (QUE), bind to DNA influencing their key physiological function. FLP is involved in the combined mode of interaction (intercalation and minor groove binding), while QUE is viewed as a minor groove binder. From a physical standpoint, experimental and theoretical studies have not so far provided a sufficiently consistent picture of the nature of interaction with DNA. Herein the sequence-dependent binding of FLP and of QUE (two representative examples of the structurally different flavonoids) with duplex DNA, containing a variety of the sequences of eight nucleotides (I: GGGGCCCC, II: GGCCGGCC, III: AAAATTTT, IV: AAGCGCTT, V: GCGCGCGC) in the 5'-strand, is investigated using a sophisticated molecular dynamics (MD) approach. For various parts (helix, backbone, bases) of the DNA structure, the change of asymptotic (in terms of an infinite length of MD simulation) configurational entropy, being the thermodynamic consequence of DNA flexibility change due to ligand binding, is explored. As far as the sequence-dependent extent of DNA flexibility change upon QUE (or FLP) binding is concerned, for the entire double helix, increased flexibility is observed for I (or I ≈ II), while increased rigidity is found to be in the order of V > III > II > IV (or III > V > IV) for the rest of sequences. For the backbone, increased rigidity in the order of V > III > II > IV > I (or III > V > IV > I > II) is generally observed. For the nucleobases, increased flexibility is determined for I and II (I > II for both ligands), while increased rigidity in the order of V ≈ III > IV (or III > V > IV) is reported for the other sequences. Of the overall increased rigidity of the DNA structure upon ligand binding that is observed for the sequences III, IV, and V, about 50-70% comes from the sugar-phosphate backbone. Noteworthy is that the increased flexibility of the entire double helix and of the complete system of nucleobases upon ligand binding is only established for sequence I. The insights are further subtly substantiated by considering the configurational entropy contributions at the level of individual nucleobase pairs and of individual nucleo-base pair steps and by analyzing the sequence dependent estimates of intra-base pair entropy and inter-base pair entropy. The GGC triplet, which is part of the central tetramer (GGCC) of I, is concluded to be critical for binding of flavonoids, while the effect of the presence of ligand to the flexibility of nucleobases is localized through the intra-base pair motion of the intercalation site and its immediate vicinity. G-rich DNA sequences with consecutive Gs going before and/or after the critical GGC code (such as I: GGGGCCCC) are proposed to be uniquely specific for flavonoids. The configurational entropy contribution, as an upper bound of the true entropy contribution to the free energy in noncovalent binding, is demonstrated to influence the fundamental discrimination (intercalation vs groove binding) of DNA-flavonoid recognition modes. Some interesting implications for the structure-based design of optimal DNA binders are discussed.

Come-back of phenanthridine and phenanthridinium derivatives in the 21st century

Phenanthridine derivatives are one of the most intensively studied families of biologically active compounds with efficient DNA binding capability. Attracting attention since DNA structure discovery (1960s), they were early recognized as a symbol of DNA intercalative binding, for many decades applied as gold-standard DNA- and RNA-fluorescent markers (ethidium bromide), probes for cell viability (propidium iodide), but also “ill-famed” for various toxic (genotoxic) and mutagenic effects. After two decades of low interest, the discovery of phenanthridine alkaloids and new studies of antiparasitic/antitumor properties of phenanthridine derivatives resulted in the strong increase of the scientific interest about the turn of this century. Here are summarized phenanthridine-related advances in the 21st century (2000-present period) with emphasis on the supramolecular interactions and bioorganic chemistry, as well as novel or improved synthetic approaches.

Atomistic account of structural and dynamical changes induced by small binders in the double helix of a short DNA

Nucleic acids are flexible molecules and their dynamical properties play a key role in molecular recognition events. Small binders interacting with DNA fragments induce both structural and dynamical changes in the double helix. We study the dynamics of a DNA dodecamer and of its complexes with Hoechst 33258, which is a minor groove binder, and with the ethidium cation, which is an intercalator, by molecular dynamics simulation. The thermodynamics of DNA-drug interaction is evaluated in connection with the structure and the dynamics of the resulting complexes. We identify and characterize the relevant changes in the configurational distribution of the DNA helix and relate them to the corresponding entropic contributions to the binding free energy. The binder Hoechst locks the breathing motion of the minor groove inducing a reduction of the configurational entropy of the helix, which amounts to 20 kcal mol(-1). In contrast, intercalations with the ethidium cation enhance the flexibility of the double helix. We show that the balance between the energy required to deform the helix for the intercalation and the gain in configurational entropy is the origin of cooperativity in the binding of a second ethidium and of anti-cooperativity in the binding of a third one. The results of our study provide an understanding of the relation between structure, dynamics and energetics in the interaction between DNA fragments and small binders, highlighting the role of dynamical changes and consequent variation of the configurational entropy of the DNA double helix for both types of binders.

Effect of various intercalators on the fenton-type oxidative cleavage of double-stranded DNA

The intensity of the linear dichroism (LD) in the absorption region of DNA (about 260 nm) decreased with time in the presence of [Fe(EDTA)](2+) (EDTA = ethylenediaminetetraacetic acid), H2O2, and ascorbate. The decrease in the LD signal indicated either an increase in flexibility, a shortening of the DNA stem, or both, owing to oxidative cleavage, and was best described by the difference between the two single-exponential-decay curves, thereby suggesting the involvement of two sequential first-order reactions. The fast reaction was assigned to cleavage of one of two DNA strands, which increased the flexibility of the DNA. The slow reaction corresponded to cleavage at or near the first cleavage site, thereby shortening the DNA stem. The presence of an intercalator, including ethidium, propidium, 9-aminoacridine, and proflavine, inhibited the first step of the cleavage reaction. One of the possible reasons for the observed inhibition might be a change in the DNA conformation near the intercalation site. Intercalation caused an unwinding and elongation of the DNA and resulted in changes in the location of the H atoms of the sugar moiety, which is known to be the main site at which hydroxyl radicals react.

Expanding the palette of phenanthridinium cations

5,6-Disubstituted phenanthridinium cations have a range of redox, fluorescence and biological properties. Some properties rely on phenanthridiniums intercalating into DNA, but the use of these cations as exomarkers for the reactive oxygen species (ROS), superoxide, and as inhibitors of acetylcholine esterase (AChE) do not require intercalation. A versatile modular synthesis of 5,6-disubstituted phenanthridiniums that introduces diversity by Suzuki–Miyaura coupling, imine formation and microwave-assisted cyclisation is presented. Computational modelling at the density functional theory (DFT) level reveals that the novel displacement of the aryl halide by an acyclic N-alkylimine proceeds by an S(N)Ar mechanism rather than electrocyclisation. It is found that the displacement of halide is concerted and there is no stable Meisenheimer intermediate, provided the calculations consistently use a polarisable solvent model and a diffuse basis set.

Targeting abasic site-containing DNA with annelated quinolizinium derivatives: the influence of size, shape and substituents

A comparative analysis showed that the type and degree of annelation as well as methyl or chloro-substitution are relevant structural features that determine the interactions of quinolizinium derivatives with abasic site-containing DNA.

Spectroscopic probe on N-H⋯N, N-H⋯O and controversial C-H⋯O contact in A-T base pair: a DFT study

DNA base pair A-T has been investigated by IR and NMR spectroscopy using DFT methods. The results have been analyzed in terms of infrared vibrational frequencies and (1)H NMR chemical shifts. Different types of interactions N-H⋯N, N-H⋯O and C-H⋯O types have been investigated in DNA base pairs. Although, previous reports argued about the third C-H⋯O type interaction in A-T base pair, such typical interaction has been analyzed thoroughly by IR and NMR spectroscopy using DFT methods. Our results show that the CH⋯O interaction in the A-T base pair is a weak interaction compared to normal hydrogen bond interactions.

Drug-DNA intercalation: from discovery to the molecular mechanism

The ability of small molecules to perturb the natural structure and dynamics of nucleic acids is intriguing and has potential applications in cancer therapeutics. Intercalation is a special binding mode where the planar aromatic moiety of a small molecule is inserted between a pair of base pairs, causing structural changes in the DNA and leading to its functional arrest. Enormous progress has been made to understand the nature of the intercalation process since its idealistic conception five decades ago. However, the biological functions were detected even earlier. In this review, we focus mainly on the acridine and anthracycline types of drugs and provide a brief overview of the development in the field through various experimental methods that led to our present understanding of the subject. Subsequently, we discuss the molecular mechanism of the intercalation process, free-energy landscapes, and kinetics that was revealed recently through detailed and rigorous computational studies.

Polycyclic azoniahetarenes: assessing the binding parameters of complexes between unsubstituted ligands and G-quadruplex DNA

Polycyclic azoniahetarenes were employed to determine the effect of the structure of unsubstituted polyaromatic ligands on their quadruplex-DNA binding properties. The interactions of three isomeric diazoniadibenzo[b,k]chrysenes (4 a-c), diazoniapentaphene (5), diazoniaanthra[1,2-a]anthracene (6), and tetraazoniapentapheno[6,7-h]pentaphene (3) with quadruplex DNA were examined by DNA melting studies (FRET melting) and fluorimetric titrations. In general, penta- and hexacyclic azoniahetarenes bind to quadruplex DNA (K(b) ≈10(6) M(-1)) even in the absence of additional functional side chains. The binding modes of 4 a-c and 3 were studied in more detail by ligand displacement experiments, isothermal titration calorimetry, and CD and NMR spectroscopy. All experimental data indicate that terminal π stacking of the diazoniachrysenes to the quadruplex is the major binding mode; however, because of different electron distributions of the π systems of each isomer, these ligands align differently in the binding site to achieve ideal binding interactions. It is proposed that tetraazonia ligand 3 binds to the quadruplex by terminal stacking with a small portion of its π system, whereas a significant part of the bulky ligand most likely points outside the quadruplex structure, and is thus partially placed in the grooves. Notably, 3 and the known tetracationic porphyrin TMPyP4 exhibit almost the same binding properties towards quadruplex DNA, with 3 being more selective for quadruplex than for duplex DNA. Overall, studies on azonia-type hetarenes enable understanding of some parameters that govern the quadruplex-binding properties of parent ligand systems. Since unsubstituted ligands were employed in this study, complementary and cooperative effects of additional substituents, which may interfere with the ligand properties, were eliminated.

Theoretical study of the stability of the DNA duplexes modified by a series of hydrophobic base analogues

The geometries of a 13 mer of a DNA double helix (5'-GCGTACACATGCG-3') were determined by molecular dynamics simulations using a Cornell et al. empirical force field. The bases in the central base pair (shown in bold) were replaced (one or both) by a series of hydrophobic base analogues (phenyl, biphenyl, phenylnaphathalene, phenylanthracene and phenylphenanthrene). Due to the large fluctuations of the systems, an average geometry could not be determined. The interaction energies of the Model A, which consisted of three central steps of a duplex without a sugar phosphate backbone, taken from molecular dynamics simulations (geometry sampled every 1 ps), were calculated by the self-consistent charge density functional based tight-binding (SCC-DFTB-D) method and were subsequently averaged. The higher the stability of the systems the higher the aromaticity of the base analogues. To estimate the desolvation energy of the duplex, the COSMO continuum solvent model was used and the calculations were provided on a larger model, Model B (the three central steps of the duplex with a sugar phosphate backbone neutralised by H atoms), taken from molecular dynamics simulations (geometry sampled every 200 ps) and subsequently averaged. The selectivity of the base analogue pairs was ascertained (Model B) by including the desolvation energy and the interaction energy of both strands, as determined by the SCC-DFTB-D method. The highest selectivity was found for a phenylphenanthrene. Replacing the nucleic acid bases with a base analogue leads to structural changes of the central pair. Only with the smallest base analogues (phenyl) does the central base pair stay planar. When passing to larger base analogues the central base pair is usually stacked.

Effects of the biological backbone on DNA-protein stacking interactions

The pi-pi stacking (face-to-face) interactions between the five natural DNA or RNA nucleobases and the four aromatic amino acids were compared using three different types of dimers: (1) a truncated nucleoside (nucleobase) stacked with a truncated amino acid; (2) a truncated nucleoside (nucleobase) stacked with an extended amino acid; and (3) a nucleoside (extended nucleobase) stacked with a truncated amino acid. Systematic (MP2/6-31G*(0.25)) potential energy surface scans reveal important information about the effects of the deoxyribose sugar and protein backbone on the structure and binding energy between truncated nucleobase and amino acid models that are typically implemented in the literature. Most notably, electrostatic and steric interactions arising from the bulkiness of the biological backbones can change the preferred relative orientations of DNA and protein pi-systems. More importantly, the protein backbone can strengthen the stacking energy (by up to 10 kJ mol(-1)), while the deoxyribose moiety can strengthen or weaken the stacking interaction depending on the positioning of the amino acid relative to the sugar residue. These effects are likely due to additional interactions between the amino acid or nucleobase ring and the backbone in the extended monomer rather than significant changes in the properties of the biological pi-systems upon model extension. Since the present work reveals that all calculated DNA-protein stacking interactions are significant and approach the strength of other noncovalent interactions between biomolecules, both pi-pi and backbone-pi interactions must be considered when attempting to gain a complete picture of DNA-protein binding.

Recent developments in the chemistry of deoxyribonucleic acid (DNA) intercalators: principles, design, synthesis, applications and trends

In the present overview, we describe the bases of intercalation of small molecules (cationic and polar neutral compounds) in DNA. We briefly describe the importance of DNA structure and principles of intercalation. Selected syntheses, possibilities and applications are shown to exemplify the importance, drawbacks and challenges in this pertinent, new, and exciting research area. Additionally, some clinical applications (molecular processes, cancer therapy and others) and trends are described.

Electrostatic contribution to the energy of binding of aromatic ligands with DNA

The method of solution of the nonlinear Poisson-Boltzmann equation was used to calculate electrostatic energy of binding of various aromatic ligands with DNA oligomers of different length. Analysis of the electrostatic contribution was made in terms of a two-step DNA binding process: formation of the intercalation cavity and insertion of the ligand. The total electrostatic energy was also partitioned into components: the energy of atom-atom coulombic interactions and the energy of interaction with surrounding water. The results indicate that electrostatic interactions are, as a whole, unfavorable to the intercalation process and that a correct analysis of structure-energy interrelation for Ligand-DNA interactions should only be accomplished at the level of the components rather than at the level of total electrostatic energy.

Interaction of acetone with single wall carbon nanotubes at cryogenic temperatures: a combined temperature programmed desorption and theoretical study

The interaction of acetone with single wall carbon nanotubes (SWCNTs) at low temperatures was studied by a combination of temperature programmed desorption (TPD) and dispersion-augmented density-functional-based tight binding (DFTB-D) theoretical simulations. On the basis of the results of the TPD study and theoretical simulations, the desorption peaks of acetone can be assigned to the following adsorption sites: (i) sites with energy of approximately 75 kJ mol (-1) ( T des approximately 300 K)endohedral sites of small diameter nanotubes ( approximately 7.7 A); (ii) sites with energy 40-68 kJ mol (-1) ( T des approximately 240 K)acetone adsorption on accessible interstitial, groove sites, and endohedral sites of larger nanotubes ( approximately 14 A); (iii) sites with energy 25-42 kJ mol (-1) ( T des approximately 140 K)acetone adsorption on external walls of SWCNTs and multilayer adsorption. Oxidatively purified SWCNTs have limited access to endohedral sites due to the presence of oxygen functionalities. Oxygen functionalities can be removed by annealing to elevated temperature (900 K) opening access to endohedral sites of nanotubes. Nonpurified, as-received SWCNTs are characterized by limited access for acetone to endohedral sites even after annealing to elevated temperatures (900 K). Annealing of both purified and as-produced SWCNTs to high temperatures (1400 K) leads to reduction of access for acetone molecules to endohedral sites of small nanotubes, probably due to defect self-healing and cap formation at the ends of SWCNTs. No chemical interaction between acetone and SWCNTs was detected for low temperature adsorption experiments. Theoretical simulations of acetone adsorption on finite pristine SWCNTs of different diameters suggest a clear relationship of the adsorption energy with tube sidewall curvature. Adsorption of acetone is due to dispersion forces, with its C-O bond either parallel to the surface or O pointing away from it. No significant charge transfer or polarization was found. Carbon black was used to model amorphous carbonaceous impurities present in as-produced SWCNTs. Desorption of acetone from carbon black revealed two peaks at approximately 140 and approximately 180-230 K, similar to two acetone desorption peaks from SWCNTs. The characteristic feature of acetone desorption from SWCNTs was peak at approximately 300 K that was not observed for carbon black. Care should be taken when assigning TPD peaks for molecules desorbing from carbon nanotubes as amorphous carbon can interfere.