Molecular dynamics study of the initial stages of catalyzed single-wall carbon nanotubes growth: force field development

Effective force fields for Ni-C interactions developed by Yamaguchi and Maruyama for the formation of metallofullerenes are modified to simulate the catalyzed growth of single-wall carbon nanotubes on Ni_n clusters with n >20, and the reactive empirical bond order Brenner potential for C-C interactions is also revised to include the effect of the metal atoms on such interactions. Figure Force field parameters for carbon-metal interactions obtained from DFT calculations in small clusters.     

Molecular modeling of the intercalation of porphyrins into α-zirconium phosphate

In this work, n-alkylamines (number of carbon atoms ranging from 3 to 10) were investigated in detail by molecular modeling as spacers for intercalating porphyrins into α-zirconium phosphate (α-ZrP). Pre-intercalated n-alkylamines can form either a flat monolayer or a canted bilayer in the gallery of α-ZrP. Based on the interlayer state and intercalative potential of the two modes in α-ZrP, it is suggested that the flat monolayer is a better spacer than the bilayer and that n-propylamine (PA) and n-butylamine (BA) in mobile monolayers are the best spacers among the n-alkylamines studied, as is also found experimentally. The intercalation behavior of TMPyP [5,10,15,20-tetrakis (1-methylpyridinium-4-yl) porphyrin] and several other porphyrins was investigated by calculating the intercalative potential. The calculated results showed that the porphyrins were densely packed in a canted monolayer model, and an increase of polarity of the substituent would facilitate the intercalation of the porphyrins. Figure Schematic representation of platform of intercalated spacers and guests taking n-butylamine and TMPyP as an example, respectively: a a flat monolayer of n-bultylamine in α-ZrP; b a canted monolayer of TMPyP in α-ZrP; c the top layer of the canted bilayer n-bultylamine in α-ZrP (the gray area indicates the amphiphilic distribution on the interface between α-ZrP layers and n-alkylamine/porphyrin).      

Theoretical study of [ XN5 ]− (X=O, S, Se, Te) systems

A series of [XN₅]⁻ (X=O, S, Se, Te) compounds has been examined with ab initio and Density Functional Theory (DFT) methods. The five-membered nitrogen ring series of structures are global minima and may exist or be characterized due to their significant dissociation barriers (29.7–32.7 kcal mol⁻¹). Nucleus-independent chemical shifts (NICS) criteria and the presence of (4n+2) π-electrons confirmed that the five-membered nitrogen ring in their structures exhibits characteristics of aromaticity. Thus, the strong stability of the five-membered nitrogen ring structures may be attributed partially to their aromaticity.      

A theoretical investigation of cytotoxic activity of celastroid triterpenoids

Quantum chemical calculations at the B3LYP/6-31G* level of theory have been carried out on 20 celastroid triterpenoids to obtain a set of molecular electronic properties and to correlate these with cytotoxic activities. The cytotoxic activities of these compounds can be roughly correlated with electronic effects related to nucleophilic addition to C(6) of the compounds: The energies of the frontier molecular orbitals (E _HOMO and E _LUMO), the HOMO-LUMO energy gap, the dipole moment, the charge on C(6), and the electrophilicity on C(6). Figure LUMO of Pristimerin.     

DFT investigation of the 1-octene metathesis reaction mechanism with the Phobcat precatalyst

The productive self-metathesis reaction of 1-octene in the presence of the Phobcat precatalyst [RuCl₂(Phoban-Cy)₂(=CHPh)] using density functional theory was investigated and compared to the Grubbs 1 precatalyst [RuCl₂(PCy₃)₂(=CHPh)]. At the GGA-PW91/DNP level, the geometry optimization of all the participating species and the PES scans of the various activation and catalytic cycles in the dissociative mechanism were performed. The formation of the catalytically active heptylidene species is kinetically and thermodynamically favored, while the formation of trans-tetradecene is thermodynamically favored.      

A density functional study towards substituent effects on anion sensing with urea receptors

Effects of substituents on anion binding in different urea based receptors have been examined using density functional (B3LYP/6-311+G**) level of theory. The complexes formed by a variety of substituted urea with a halide anion (fluoride) and an oxy-anion (acetate) have been calculated. The stronger complexes were predicted for receptors with fluoride ion than that of acetate ion, however, in water the preference was found to be reversed. The pK _a calculations showed the preferred sites of deprotonation for positional isomers, while interacting with anions. The position of the substituent in the receptor, however, could change the preferred sites of deprotonation compared to the site predicted with pK _a values.      

Alpha particle chemistry. On the formation of stable complexes between He2+ and other simple species: implications for atmospheric and interstellar chemistry

The possibility that stable complexes may be formed between alpha particles (He²⁺) and small molecules is investigated using QCISD quantum mechanical calculations. Implications for their presence in the terrestrial atmosphere and/or in interstellar space are discussed. Figure Optimized structure of a stable H₂OHe²⁺ complex     

Phenolic esters with potential anticancer activity - the structural variable

The conformational preferences of several potential anticancer dihydroxycinnamic esters with a variable length alkyl chain were studied by quantum-mechanical (DFT) calculations (both for the isolated molecule and for aqueous solutions). The orientation of the hydroxyl ring substituents and of the alkyl ester moiety relative to the carbonyl group showed these to be the most determinant factors for the overall stability of this type of phenolic systems, strongly dependent on an effective π-electron delocalization. Compared to the parent caffeic acid (dihydroxycinnamic acid), esterification was found to lead to a higher conformational freedom, and to affect mainly the energy barrier corresponding to the (O=)C-OR internal rotation. No particular differences were verified to occur upon lengthening of the ester alkyl chain, except when this is branched instead of linear. The vibrational spectra of the whole series of compounds were simulated, based on their calculated harmonic vibrational frequencies, and a preliminary assignment was performed. Figure Schematic representation of the dihydroxycinnamic esters studied in the present work and of the main internal rotations affecting the overall stability of the molecules. (R=(CH₂)_n, n = 0,1,2,3,7,11 for MC, EC, PC, BC, OC and DC, respectively; R=(CHCH₃) for IPC. The atom numbering is included, with the exception of the alkyl ester group)     

Interactions in pseudorotoxanes based on crown ether-secondary ammonium motifs. A theoretical study

A theoretical analysis of the nature of the interactions in dibenzo[24]crown-8 (DB24C8)-n-dibutylammonium (DBM)—pseudorotaxane complex at the MP2 and DFT levels shows that the main contribution to the binding energy is the electrostatic interaction with moderate (20–25%) correlation stabilization. The total binding energy in the DB24C8-DBM complex represents a sum of the binding energies of two NH–O and one CH–O hydrogen bonds and the latter constitutes about 25% of the total interaction energy, giving the total binding energy of −41.2 kcal mol⁻¹ at the BHandHLYP/6-311++G** level. Deprotonation of the DB24C8-DBM complex reduces the binding energy by some 50 kcal mol⁻¹, giving metastable complexes DB24C8-DBA-1 or DB24C8-DBA-2, which will dissociate to give free crown ether and n-dibutylamine because of the strong exchange repulsion that prevails in neutral complexes. Figure Formation of DB24C8-DBM pseudorotoxane complex     

Study of a ligand complexed with Cdk2/Cdk4 by computer simulation

Cyclin-dependent kinases (Cdks) play important roles in the regulation of the cell cycle. Their inhibitors have entered clinical trials to treat cancer. Very recently, Davis et al. (Nat Struct Biol 9:745–749, 2002) have found a ligand NU6102, which has a high affinity with cyclin-dependent kinase 2 (K _i =6 nM) but a low affinity with cyclin-dependent kinase 4 (K _i =1,600 nM). To understand the selectivity, we use homology modeling, molecular docking, molecular dynamics and free-energy calculations to analyze the interactions. A rational 3D model of the Cdk4–NU6102 complex is built. Asp86 is a key residue that recognizes NU6102 more effectively with Cdk2 rather than Cdk4. Good binding free energies are obtained. Energetic analysis reveals that van der Waals interaction and nonpolar contributions to solvent are favorable in the formation of complexes and the sulfonamide group of the ligand plays a crucial role for binding selectivity between Cdk2 and Cdk4. Figure Two-dimensional representative for the interacting model of NU6102 complexed with the Cdk4 from a predicted structure by LIGPLOT.      

A computational investigation of sulfur-containing heterocyclic components from the anal sac secretions of Mustela species

A computational investigation of the sulfur-containing heterocyclic components (substituted thietanes and 1,2-dithiolanes) of Mustela anal sac secretions has been carried out. A cluster analysis of the chemical compositions of Mustela anal sac volatiles reveals little similarity with established phylogenetic relationships between members of the genus. Ab initio calculations [MP2/6–311++G(2df,2p)//B3LYP/6–311++G**] show the lowest-energy C₅H₁₀S isomeric thietane to be 2,2-dimethylthietane, which is also the most abundant of the Mustela thietanes. Similarly, 3,3-dimethyl-1,2-dithiolane is the lowest-energy C₅H₁₀S₂ compound. 2-n-Propylthietane is the highest-energy C₆H₁₂S compound, but the most abundant Mustela C₆H₁₂S compound produced, whereas cis-2-ethyl-4-methylthietane, the lowest-energy C₆H₁₂S thietane, has never been observed in Mustela anal sac secretions. A molecular docking analysis of the Mustela sulfur-containing heterocycles into both porcine and bovine odorant binding proteins reveals the interactions of the docked ligands with the proteins to be largely hydrophobic, and have binding energies generally lower than typical odorant molecules such as linalool or eugenol. Figure Mustela anal sac volatile components, 2,2-dimethylthietane and cis-3,4-dimethyl–1,2-dithiolane.     

Pattern recognition based on color-coded quantum mechanical surfaces for molecular alignment

A pattern recognition algorithm for the alignment of drug-like molecules has been implemented. The method is based on the calculation of quantum mechanical derived local properties defined on a molecular surface. This approach has been shown to be very useful in attempting to derive generalized, non-atom based representations of molecular structure. The visualization of these surfaces is described together with details of the methodology developed for their use in molecular overlay and similarity calculations. In addition, this paper also introduces an additional local property, the local curvature (C _L), which can be used together with the quantum mechanical properties to describe the local shape. The method is exemplified using some problems representing common tasks encountered in molecular similarity. Figure Molecular surfaces for Lorazepam (left) and Diazepam (right)     

TCNE-aniline charge transfer complex: ab initio and TDDFT investigations in gas phase

The geometric and electronic structure of tetracyanoethylene (TCNE)-aniline (donor-acceptor type) complex has been investigated in gas phase using ab initio and time dependent density functional theory calculations. Both the above calculations predict a composed structure for the complex, in which the interacting site is a C≡N and C=C bond center in the TCNE and, –NH₂ and π-electrons of aniline. The N atom of aniline is oriented toward the TCNE molecule. The charge transfer transition energy, estimated by calculating the ground-to-excited state transition electric dipole moments of the complex, agree well with the reported experimental value in chloroform medium. TCNE-aniline at ground state. TCNE-aniline at excited state     

Molecular models of protein kinase 6 from Plasmodium falciparum

Cyclin-dependent kinases (CDKs) have been identified as potential targets for development of drugs, mainly against cancer. These studies generated a vast library of chemical inhibitors of CDKs, and some of these molecules can also inhibit kinases identified in the Plasmodium falciparum genome. Here we describe structural models for Protein Kinase 6 from P. falciparum (PfPK6) complexed with Roscovitine and Olomoucine. These models show clear structural evidence for differences observed in the inhibition, and may help designing inhibitors for PfPK6 generating new potential drugs against malaria. Figure Ribbon diagram of PfPK6 complexed with a roscovitine and b olomoucine     

The platinum-olefin binding energy in series of (PH<Subscript>3</Subscript>)<Subscript>2</Subscript>Pt(olefin) complexes - a theoretical study

Theoretical investigation of Pt(0)-olefin organometallic complexes containing tertiary phosphine ligands was focused on the strength of platinum-olefin electronic interaction. DFT theoretical study of electronic effects in a substantial number of ethylene derivatives was evaluated in terms of the Pt-olefin binding energy using MP2 correlation theory. Organometallics bearing coordinated olefins with general formula (R¹R²C = CR³R⁴)Pt(PH₃)₂ [R = various substituents] had been selected, including olefins containing both electron-donor substituents as well as electron-withdrawing groups. The stability of the corresponding complexes increases with a strengthening electron-withdrawal ability of the olefin substituents. Figure Representation of (CH₂ = CHR)Pt(PPh₃)₂ and the stability chart     

G2, G3, and complete basis set calculations on the thermodynamic properties of triazane

As a follow-up study to our study on tetrazane (N₄H₆), we present computed thermodynamic properties of triazane (N₃H₅). Calculated properties include optimized geometries, infrared vibrations, enthalpy of formation, enthalpy of combustion, and proton affinities. We have also mapped the potential energy surface as the molecule is rotated about the N-N bond. We have predicted a specific enthalpy of combustion for triazane of about -20 kJ g⁻¹. Figure Schematic diagram of the dielectric barrier discharge (left) and typical temporal profiles of voltage and current, as obtained from the simulations (right)     

Experimental and theoretical investigations of spectroscopic properties of azobenzene derivatives in solution

The UV-Vis spectra of series of polymethylmethacrylate (PMMA) copolymers with attached trans-azobenzene derivatives were measured in 1,1,2-trichloroethane. In order to gain some insight into the recorded spectra, the quantum chemical calculations were performed for the substituted azobenzenes using both configuration interaction with single excitations method (CIS) as well as density functional theory (DFT) with B3LYP and PBE0 functionals. The calculations were performed in solvent. In particular, we found that the PBE0 excitation energies are in very good agreement with the experimental data. Figure The plots of orbital contour surfaces for molecule II. The molecular orbitals were calculated at the PBE0/6-311++G(d,p) level of theory. The upper plot presents contour surface of HOMO and the lower presents contour surface of LUMO. Shown are the contour surfaces of orbital amplitude 0.04 (red) and -0.04 (blue)