Flephedrone is an analogue of cathinone - chemically similar to ephedrine, cathine and other amphetamines. Conformations of all isomers of flephedrone have been studied at the quantum-chemical level. Calculations have been performed using DFT and MP2 methods with two basis sets - 6-31G and 6-31G(d,p). Results show that there are low energy conformers for the ortho, meta, and para isomers that are connected by way of low-barrier transition states. Boltzmann distribution of population predicts the highest population for the 1-meta conformer with a 10 % increase in solution. The molecular electrostatic potential surface data for each molecule has been calculated revealing likely reaction sites.

Electronic supplementary material

The online version of this article (doi:10.1007/s00894-012-1673-z) contains supplementary material, which is available to authorized users.

Increase of the atmospheric concentration of halogenated organic compounds is partially responsible for a change of the global climate. In this work we have investigated the interaction between halogenated ether and water, which is one of the most important constituent of the atmosphere. The structures of the complexes formed by the two most stable conformers of enflurane (a volatile anaesthetic) with one and two water molecules were calculated by means of the counterpoise CP-corrected gradient optimization at the MP2/6–311++G(d,p) level. In these complexes the CH…Ow hydrogen bonds are formed, with the H…Ow distances varying between 2.23 and 2.32 Å. A small contraction of the CH bonds and the blue shifts of the ν(CH) stretching vibrations are predicted. There is also a weak interaction between one of the F atoms and the H atom of water, with the Hw…F distances between 2.41 and 2.87 Å. The CCSD(T)/CBS calculated stabilization energies in these complexes are between −5.89 and −4.66 kcal mol−1, while the enthalpies of formation are between −4.35 and −3.22 kcal mol−1. The Cl halogen bonding between enflurane and water has been found in two complexes. The intermolecular (Cl···O) distance is smaller than the sum of the corresponding van der Waals radii. The CCSD(T)/CBS stabilization energies for these complexes are about −2 kcal mol−1.

FigureComplex between enflurane and water molecules

One of the factors limiting the search of new compounds based on the structure of target proteins involved in diseases is the limited amount of target structural information. Great advances in the search for lead compounds could be achieved to find new cavities in protein structures that are generated using well established computational chemistry tools. In the case of dengue, the discovery of pockets in the crystallographic structure of the E protein has contributed to the search for lead compounds aimed at interfering in conformational transitions involved in the pH-dependent fusion process. This is a complex mechanism triggered by the acid pH of the endosomes that leads to the initial changes in the E protein assembly at the virus surface. In the present work, an arrangement of three ectodomain portions of the E protein present on the surface of the mature dengue virus was studied through long all-atom molecular dynamics simulations with explicit solvent. In order to identify new pockets and to evaluate the influence of the acid pH on these pockets, the physiological neutral pH conditions and the acid pH of the endosomes that trigger the fusion process were modeled. Several pockets presenting pH-dependent characteristics were found in the contact regions between the chains. Pockets at the protein-protein interfaces induced by a monomer in another monomer were also found. Some of the pockets are good candidates for the design of lead compounds that could interfere in the rearrangements in E proteins along the fusion process contributing to the development of specific inhibitors of the dengue disease.

Electronic supplementary material

The online version of this article (doi:10.1007/s00894-012-1687-6) contains supplementary material, which is available to authorized users.

Modern semiempirical methods are of sufficient accuracy when used in the modeling of molecules of the same type as used as reference data in the parameterization. Outside that subset, however, there is an abundance of evidence that these methods are of very limited utility. In an attempt to expand the range of applicability, a new method called PM7 has been developed. PM7 was parameterized using experimental and high-level ab initio reference data, augmented by a new type of reference data intended to better define the structure of parameter space. The resulting method was tested by modeling crystal structures and heats of formation of solids. Two changes were made to the set of approximations: a modification was made to improve the description of noncovalent interactions, and two minor errors in the NDDO formalism were rectified. Average unsigned errors (AUEs) in geometry and ΔHf for PM7 were reduced relative to PM6; for simple gas-phase organic systems, the AUE in bond lengths decreased by about 5 % and the AUE in ΔHf decreased by about 10 %; for organic solids, the AUE in ΔHf dropped by 60 % and the reduction was 33.3 % for geometries. A two-step process (PM7-TS) for calculating the heights of activation barriers has been developed. Using PM7-TS, the AUE in the barrier heights for simple organic reactions was decreased from values of 12.6 kcal/mol-1 in PM6 and 10.8 kcal/mol-1 in PM7 to 3.8 kcal/mol-1. The origins of the errors in NDDO methods have been examined, and were found to be attributable to inadequate and inaccurate reference data. This conclusion provides insight into how these methods can be improved.

Electronic supplementary material

The online version of this article (doi:10.1007/s00894-012-1667-x) contains supplementary material, which is available to authorized users.

The preliminary cytotoxic effect of 4-ethoxycarbonylmethyl-1-(piperidin-4-ylcarbonyl)-thiosemicarbazide hydrochloride (1)—a potent topoisomerase II inhibitor—was measured using a MTT assay. It was found that the compound decreased the number of viable cells in both estrogen receptor-positive MCF-7 and estrogen receptor-negative MDA-MB-231breast cancer cells, with IC50 values of 146 ± 2 and 132 ± 2 μM, respectively. To clarify the molecular basis of the inhibitory action of 1, molecular docking studies were carried out. The results suggest that 1 targets the ATP binding pocket.

Figure4-ethoxycarbonylmethyl-1-(piperidin-4-ylcarbonyl)-thiosemicarbazide hydrochloride

Electronic supplementary material

The online version of this article (doi:10.1007/s00894-012-1679-6) contains supplementary material, which is available to authorized users.

The performance of the reaction-field method of electrostatics is tested in molecular dynamics simulations of protein human interleukin-4 and a short DNA fragment in explicit solvent. Two truncation schemes are considered: one based on the position of atomic charges in water molecules and the other on the position of groups of charges. The group-based truncation leads to the melting of the DNA double helix. In contrast, the atom-based truncation maintains the helical structure intact. Similarly for the protein, the group-based truncation leads to an unfolding at pH 2 while the atom-based truncation produces stable trajectories at low and normal pH, in agreement with experiment. Artificial repulsion between charged residues associated with the group-based truncation is identified as the microscopic reason behind unfolding of the protein. Implications of different truncation schemes in reaction-field simulations of biomolecules are discussed.

The HIV-1 accessory protein Nef plays an active role in the pathogenesis of AIDS by its numerous cellular interactions that facilitate the release of virus particles. This 27 kDa protein is required for maintenance of the viral replication in HIV, and is also known to contribute to immune evasion, blocking of apoptosis in virus-infected cells and enhancement of virus infectivity. Nef has been shown to be secreted and is present on the surface of virus-infected cells. Recent studies from our laboratory have shown that the Nef protein is secreted from nef-transfected and HIV-1-infected cells in small exosome-like vesicles (40–100 nm diam.) that do not contain virions. We have identified three amino-terminal domains of Nef as necessary for secretion: (i) the four arginine residues (17,19,21, 22) comprising the basic region; (ii) the phosphofurin acidic cluster sequence (PACS) composed of four glutamic acid residues (61–64); (iii) a previously unknown motif spanning amino acid residues 65–69 (VGFPV) which we named the secretion modification region (SMR). In this study, we have used population-based phylogeny data and sequence analysis to characterize the conservation of the Nef SMR domain that regulates vesicle secretion. We have performed in silico computational chemistry analysis involving molecular dynamic structure modeling of mutations in the SMR motif. Sequence analysis of Nef from HIV-1-infected patients, including slow progressors (SP), long term progressors (LTP) and long term non-progressors (LTNP) demonstrated 99 % conservation of the Nef SMR motif. Computational analysis including modeling of wild-type HIV-1 Nef and V66A Nef SMR mutant using structural homology and molecular dynamics of ligand-associated interactions indicated significant structural changes in the Nef mutant, thus supporting the importance of the SMR domain for mediating Nef vesicle secretion.

With improvements in computer speed and algorithm efficiency, MD simulations are sampling larger amounts of molecular and biomolecular conformations. Being able to qualitatively and quantitatively sift these conformations into meaningful groups is a difficult and important task, especially when considering the structure-activity paradigm. Here we present a study that combines two popular techniques, principal component (PC) analysis and clustering, for revealing major conformational changes that occur in molecular dynamics (MD) simulations. Specifically, we explored how clustering different PC subspaces effects the resulting clusters versus clustering the complete trajectory data. As a case example, we used the trajectory data from an explicitly solvated simulation of a bacteria’s L11·23S ribosomal subdomain, which is a target of thiopeptide antibiotics. Clustering was performed, using K-means and average-linkage algorithms, on data involving the first two to the first five PC subspace dimensions. For the average-linkage algorithm we found that data-point membership, cluster shape, and cluster size depended on the selected PC subspace data. In contrast, K-means provided very consistent results regardless of the selected subspace. Since we present results on a single model system, generalization concerning the clustering of different PC subspaces of other molecular systems is currently premature. However, our hope is that this study illustrates a) the complexities in selecting the appropriate clustering algorithm, b) the complexities in interpreting and validating their results, and c) by combining PC analysis with subsequent clustering valuable dynamic and conformational information can be obtained.

Electronic supplementary material

The online version of this article (doi:10.1007/s00894-012-1563-4) contains supplementary material, which is available to authorized users.

A polymer molecule (represented by a statistical chain) end-grafted to a topologically rough surface was studied by static MC simulations. A modified self-avoiding walk on a cubic lattice was used to model the polymer in an athermal solution. Different statistical models of surface roughness were applied. Conformational entropies of chains attached to uncorrelated Gaussian, Brownian, and fractional Brownian surfaces were calculated. Results were compared with the predictions of a simple analytical model of a macromolecule end-grafted to a fractal surface.

FigureVisualization of SAW generated by the (023) algorithm on a 3D cubic lattice

Electronic supplementary material

The online version of this article (doi:10.1007/s00894-012-1546-5) contains supplementary material, which is available to authorized users.

Molecular dynamics simulations and free energy calculation have been performed to study how the single-chain variable fragment (scFv) binds Methamphetamine (METH) and Amphetamine (AMP). The structures of the scFv:METH and the scFv:AMP complexes are analyzed by examining the time-dependence of their RMSDs, by analyzing the distance between some key atoms of the selected residues, and by comparing the averaged structures with their corresponding crystallographic structures. It is observed that binding an AMP to the scFv does not cause significant changes to the binding pocket of the scFv:ligand complex. The binding free energy of scFv:AMP without introducing an extra water into the binding pocket is much stronger than scFv:METH. This is against the first of the two scenarios postulated in the experimental work of Celikel et al. Protein Science 18, 2336 (2009). However, adding a water to the AMP (at the position of the methyl group of METH), the binding free energy of the scFv:AMP-H2O complex, is found to be significantly weaker than scFv:METH. This is consistent with the second of the two scenarios given by Celikel et al. Decomposition of the binding energy into ligand-residue pair interactions shows that two residues (Tyr175 and Tyr177) have nearly-zero interactions with AMP in the scFv:AMP-H2O complex, whereas their interactions with METH in the scFv:METH complex are as large as -0.8 and -0.74 kcal/mol. The insights gained from this study may be helpful in designing more potent antibodies in treating METH abuse.

Protein surface roughness is a structural property associated with ligand-protein and protein-protein binding interfaces. In this work we apply for the first time the concept of surface roughness, expressed as the fractal dimension, to address structure and function of G protein-coupled receptors (GPCRs) which are an important group of drug targets. We calculate the exposure ratio and the fractal dimension for helix-forming residues of the β2 adrenergic receptor (β2AR), a model system in GPCR studies, in different conformational states: in complex with agonist, antagonist and partial inverse agonists. We show that both exposure ratio and roughness exhibit periodicity which results from the helical structure of GPCRs. The pattern of roughness and exposure ratio of a protein patch depends on its environment: the residues most exposed to membrane are in general most rough whereas parts of receptors mediating interhelical contacts in a monomer or protein complex are much smoother. We also find that intracellular ends (TM3, TM5, TM6 and TM7) which are relevant for G protein binding and thus receptor signaling, are exposed but smooth. Mapping the values of residual fractal dimension onto receptor 3D structures makes it possible to conclude that the binding sites of orthosteric ligands as well as of cholesterol are characterized with significantly higher roughness than the average for the whole protein. In summary, our study suggests that identification of specific patterns of roughness could be a novel approach to spot possible binding sites which could serve as original drug targets for GPCRs modulation.

FigureThe significance of surface roughness for protein structure and function.

Quantum-chemical calculations {DFT(B3LYP)/6-311+G(d,p)} were performed for all possible tautomers (aromatic and nonaromatic) of neutral 2- and 4-aminopyridines and their oxidized and reduced forms. One-electron oxidation has no important effect on the tautomeric preference for 2-aminopyridine. The amine tautomer is favored. However, oxidation increases the stability of the imine NH tautomer, and its contribution in the tautomeric mixture cannot be neglected. In the case of 4-aminopyridine, one-electron oxidation increases the stability of both the amine and imine NH tautomers. Consequently, they possess very close energies. As major tautomers, they dictate the composition of the tautomeric mixture. The CH tautomers may be considered as very rare forms for both neutral and oxidized aminopyridines. A reverse situation takes place for the reduced forms of aminopyridines. One-electron reduction favors the C3 atom for the labile proton for both aminopyridines. This may partially explain the origin of the CH tautomers for the anionic states of nucleobases containing the exo NH2 group.

The in vitro antifungal potency of six series of 4-arylthiosemicarbazides was evaluated. Two isoquinoline derivatives with an ortho-methoxy or ortho-methyl group at the phenyl ring were the most potent antifungal agents. Molecular modeling studies and docking of all 4-arylthiosemicarbazides into the active sites of sterol 14α-demethylase (CYP51), topoisomerase II (topo II), l-glutamine: d-fructose-6-phosphate amidotransferase (GlcN-6-P), secreted aspartic proteinase (SAP), N-myristoyltransferase (NMT), and UDP-N-acetylmuramoyl-l-alanine:d-glutamate ligase (MurD) indicated the importance of both structural and electronic factors in ligand recognition and thus for the antifungal effectiveness of 4-arylthiosemicarbazides. A possible antifungal target was identified (NMT) and isoquinoline-thiosemicarbazides showed more favorable affinity than the native ligand.

FigureElectrostatic potential surface of isoquiniline derivative compound 6o with antifungal activity

Electronic supplementary material

The online version of this article (doi:10.1007/s00894-012-1420-5) contains supplementary material, which is available to authorized users.

A comprehensive ab initio investigation using coupled cluster theory with the aug-cc-pVnZ, n = D,T basis sets is carried out to identify distinct structures of the Al4H14— cluster anion and to evaluate its fragmentation stability. Both thermodynamic and mechanistic aspects of the fragmentation reactions are studied. The observation of this so far the most hydrogenated aluminum tetramer was reported in the recent mass spectrometry study of Li et al. (2010) J Chem Phys 132:241103–241104. The four Al4H14— anion structures found are chain-like with the multiple-coordinate Al center and can be viewed approximately as comprising Al2H7— and Al2H7 moieties. Locating computationally some of the Al4H14— minima on the correlated ab initio potential energy surfaces required the triple-zeta quality basis set to describe adequately the Al multi-coordinate bonding. For the two most stable Al4H14— isomers, the mechanism of their low-barrier interconversion is described. The dissociation of Al4H14— into the Al2H7— and Al2H7 units is predicted to require 20-22 (10-13) kcal mol-1 in terms of ΔH (ΔG) estimated at T = 298.15 K and p = 1 atm. However, Al4H14— is found to be a metastable species in the gas phase: the H2 loss from the radical moiety of its most favorable isomer is exothermic by 18 kcal mol-1 in terms of ΔH (298.15 K) and by 25 kcal mol-1 in terms of ΔG(298.15 K), with the enthalpic/free energy barrier involved being less than 1 kcal mol-1. By contrast with alane Al4H14—, only a weakly bound complex between Ga4H12— and H2 has been identified for the gallium analogue using the relativistic effective core potential.

COPI, a 600 kD heptameric complex (consisting of subunits α, β, γ, δ, ε, ζ, and β′) “coatomer,” assembles non-clathrin-coated vesicles and is responsible for intra-Golgi and Golgi-to-ER protein trafficking. Here, we report the three-dimensional structures of the entire sequences of yeast Sec21 (γ-COPI mammalian ortholog), yeast Ret3 (ζ-COPI mammalian ortholog), and the results of successive molecular dynamics investigations of the subunits and assembly based on a protein–protein docking experiment. The three-dimensional structures of the subunits in their complexes indicate the residues of the two subunits that impact on assembly, the conformations of Ret3 and Sec21, and their binding orientations in the complexed state. The structure of the appendage domain of Sec21, with its two subdomains—the platform and the β-sandwich, was investigated to explore its capacity to bind to accessory protein recruitment motifs. Our study shows that a binding site on the platform is capable of binding the Eps15 DPF and epsin DPW2 peptides, whereas the second site on the platform and the site on the β-sandwich subdomain were found to selectively bind to the amphiphysin FXDXF and epsin DPW1 peptides, respectively. Identifying the regions of both the platform and sandwich subdomains involved in binding each peptide motif clarifies the mechanism through which the appendage domain of Sec21 engages with the accessory proteins during the trafficking process of non-clathrin-coated vesicles.

Electronic supplementary material

The online version of this article (doi:10.1007/s00894-011-1324-9) contains supplementary material, which is available to authorized users.

A detailed analysis of the selected DFT functionals for the calculations of interaction-induced dipole moment, polarizability and first-order hyperpolarizability has been carried out. The hydrogen-bonded model chains consisting of HF, H2CO and H3N molecules have been chosen as a case study. The calculations of the components of the static electric properties using the diffuse Dunning’s basis set (aug-cc-pVDZ) have been performed employing different types of density functionals (B3LYP, LC-BLYP, PBE0, M06-2X and CAM-B3LYP). Obtained results have been compared with those gained at the CCSD(T) level of theory. The counterpoise correction scheme, namely site-site function counterpoise, has been applied in order to eliminate basis set superposition error. The performed tests allow to conclude that the DFT functionals can provide a useful tool for prediction of the interaction-induced electric properties, however a caution has to be urged to their decomposition to the two- and many-body terms.

Figure Hydrogen-bonded model chains consisting of HF, H2CO and H3N molecules

Quantum-chemical computations were used to investigate the structure–antioxidant parameter relationships of α-lipoic acid and its natural metabolites bisnorlipoic acid and tetranorlipoic acid in their oxidized and reduced forms. The enantiomers of lipoic and dihydrolipoic acid were optimized using the B3LYP/6-311+G(3df,2p), B3LYP/aug-cc-pVDZ and MP2(full)/6-31+G(d,p) levels of theory as isolated molecules and in the presence of water. The geometries of the metabolites and the values of their antioxidant parameters (proton affinity, bond dissociation enthalpy, adiabatic ionization potential, spin density, and the highest occupied molecular orbital energy) were calculated at the B3LYP/6-311+G(3df,2p) level of theory. The results obtained reveal similarities between these structures: a pentatomic, nonaromatic ring is present in the oxidized forms, while an unbranched aliphatic chain (as found in saturated fatty acids) is present in both the oxidized and the reduced forms. Analysis of the spin density and the highest occupied molecular orbital energy revealed that the SH groups exhibited the greatest electron-donating activities. The values obtained for the proton affinity, bond dissociation enthalpy and adiabatic ionization potential indicate that the preferred antioxidant mechanisms for α-lipoic acid and its metabolites are sequential proton loss electron transfer in polar media and hydrogen atom transfer in vacuum.

Electronic supplementary material

The online version of this article (doi:10.1007/s00894-011-1306-y) contains supplementary material, which is available to authorized users.

Lincosamides are a class of antibiotics used both in clinical and veterinary practice for a wide range of pathogens. This group of drugs inhibits the activity of the bacterial ribosome by binding to the 23S RNA of the large ribosomal subunit and blocking protein synthesis. Currently, three X-ray structures of the ribosome in complex with clindamycin are available in the Protein Data Bank, which reveal that there are two distinct conformations of the pyrrolidinyl propyl group of the bound clindamycin. In this work, we used quantum mechanical methods to investigate the probable conformations of clindamycin in order to explain the two binding modes in the ribosomal 23S RNA. We studied three lincosamide antibiotics: clindamycin, lincomycin, and pirlimycin at the B3LYP level with the 6-31G** basis set. The focus of our work was to connect the conformational landscape and electron densities of the two clindamycin conformers found experimentally with their physicochemical properties. For both functional conformers, we applied natural bond orbital (NBO) analysis and the atoms in molecules (AIM) theory, and calculated the NMR parameters. Based on the results obtained, we were able to show that the structure with the intramolecular hydrogen bond C=O…H–O is the most stable conformer of clindamycin. The charge transfer between the pyrrolidine-derivative ring and the six-atom sugar (methylthiolincosamide), which are linked via an amide bond, was found to be the dominant factor influencing the high stability of this conformer.

FigureMolecular graph of more stable conformer of clindamycin.

In the absence of experimental structures, comparative modeling continues to be the chosen method for retrieving structural information of target proteins. However, models lack the accuracy of experimental structures. Alignment error and structural divergence (between target and template) influence model accuracy the most. Here, we examine the potential additional impact of backbone geometry, as our previous studies have suggested that the structural class (all-α, αβ, all-β) of proteins may influence the accuracy of their models. In the twilight zone (sequence identity ≤ 30%) and at a similar level of target-template divergence, model accuracy of proteins does indeed follow the trend all-α > αβ > all-β. This is mainly a result of the alignment accuracy following the same trend (all-α > αβ > all-β) with backbone geometry playing only a minor role. Differences in the diversity of sequences belonging to different structural classes leads to the observed accuracy differences thus enabling a priori accuracy estimates of alignments/models in a class-dependent manner. This study provides a systematic description and quantification of structural class-dependent effect in comparative modeling. The study also suggests datasets for large-scale sequence/structure analyses should have equal representation of different structural classes to avoid class-dependent bias.

A molecular modeling strategy is proposed to describe the temperature (T) dependence of solubility parameter (δ) for the amorphous polymers which exhibit glass-rubber transition behavior. The commercial forcefield “COMPASS” is used to support the atomistic simulations of the polymer. The temperature dependence behavior of δ for the polymer is modeled by running molecular dynamics (MD) simulation at temperatures ranging from 250 up to 650 K. Comparing the MD predicted δ value at 298 K and the glass transition temperature (Tg) of the polymer determined from δ–T curve with the experimental value confirm the accuracy of our method. The MD modeled relationship between δ and T agrees well with the previous theoretical works. We also observe the specific volume (v), cohesive energy (Ucoh), cohesive energy density (ECED) and δ shows a similar temperature dependence characteristics and a drastic change around the Tg. Meanwhile, the applications of δ and its temperature dependence property are addressed and discussed.

Structure and energy calculations of pristine and COOH-modified model single wall carbon nanotubes (SWCNTs) of different length were performed at B3LYP/6-31G* level of theory. From 1 to 9 COOH groups were added at the end of the nanotube. The differences in structure and energetics of partially and fully functionalized SWCNTs at one end of the nanotube are observed. Up to nine COOH groups could be added at one end of (9,0) zigzag SWCNT in case of full functionalization. However, for (5,5) armchair SWCNT, the full functionalization was impossible due to steric crowding and rim deformation. The dependence of substituent attachment energy on the number of substituents at the carbon nanotube rim was observed.

FigureStructure and energy calculations of pristine and COOH-modified model single wall carbon nanotubes (SWCNTs) of different length were performed at B3LYP/6-31G* level of theory and from 1 to 9 COOH groups was added at the end of the nanotube

Enzootic pneumonia caused by Mycoplasma hyopneumoniae is a major constraint to efficient pork production throughout the world. This pathogen has a small genome with 716 coding sequences, of which 418 are homologous to proteins with known functions. However, almost 42% of the 716 coding sequences are annotated as hypothetical proteins. Alternative methodologies such as threading and comparative modeling can be used to predict structures and functions of such hypothetical proteins. Often, these alternative methods can answer questions about the properties of a model system faster than experiments. In this study, we predicted the structures of seven proteins annotated as hypothetical in M. hyopneumoniae, using the structure-based approaches mentioned above. Three proteins were predicted to be involved in metabolic processes, two proteins in transcription and two proteins where no function could be assigned. However, the modeled structures of the last two proteins suggested experimental designs to identify their functions. Our findings are important in diminishing the gap between the lack of annotation of important metabolic pathways and the great number of hypothetical proteins in the M. hyopneumoniae genome.

The structures of ideal armchair (5,5) single-wall carbon nanotubes (SWCNTs) of different lengths (3.7, 8.8, and 16.0 Å for C40H20, C80H20, and C140H20) and with 1–10 hydroxyl groups at the end of the nanotube were fully optimized at the B3LYP/3-21G level, and in some cases at the B3LYP/6-31G* level, and the energy associated with the attachment of the OH substituent was determined. The OH-group attachment energy was compared with the OH functionalization of phenanthrene and picene models and with previous results for zigzag (9.0) SWCNT systems. In comparison to zigzag SWCNTs, the armchair form is more (by about 5 to 10 kcal mol−1) reactive toward hydroxylation.

FigureThe structures of ideal armchair (5,5) single-wall carbon nanotubes (SWCNTs) of different lengths (3.7, 8.8, and 16.0 Å for C40H20, C80H20, and C140H20) and with 1–10 hydroxyl groups at the end of the nanotube were fully optimized at the B3LYP/3-21 G level, and in some cases at the B3LYP/6-31 G* level, and the energy associated with the attachment of the OH substituent was determined.

The present study analyzed binding of Cu2+ to tetrapeptides in water solution at several levels of theoretical approximation. The methods used to study the energetic and structural properties of the complexes in question include semiempirical hamiltonians, density functional theory as well as ab initio approaches including electron correlation effects. In order to shed light on the character of interactions between Cu2+ and peptides, which are expected to be mainly electrostatic in nature, decomposition of interaction energy into physically meaningful components was applied.