van der Waals potential barrier for cobaltocene encapsulation into single-walled carbon nanotubes: classical molecular dynamics and ab initio study

In this work, we carried out geometry optimisations and classical molecular dynamics for the problem of cobaltocene (CC) encapsulation into different carbon nanotubes (CNTs) ((7,7), (8,8), (13,0) and (14,0) tubes were used). CCs are molecules composed of two aromatic pentagonal rings (C5H5) sandwiching one cobalt atom. From our simulation results, we observed that CC was encapsulated into CNTs (8,8), (13,0) and (14,0). However, for CNT (7,7), the encapsulation could not occur, in disaggrement with some previous works in the literature. Our results show that the encapsulation process is mainly governed by van der Waals potential barriers.     

Acid-controlled release complexes of podophyllotoxin and etoposide with acyclic cucurbit[n]urils for low cytotoxicity

The targeted or responsive systems are appealing therapeutic platforms for the development of next-generation precision medications. So, we design and prepare acid-controlled release complexes of podophyllotoxin (POD) and etoposide (VP-16) with pH-labile acyclic cucurbit[n]urils, and their characteristics and inclusion complexation behaviors were investigated via fluorescence spectroscopy, nuclear magnetic resonance and X-ray power diffraction. Cells incubated with complexes have been analyzed by high-content analysis (HCA), and cytotoxicity tests have been completed by MTT assay. The results showed that complexes with different binding constants can release the drug substance in the physiological pH environment of cancer cells, maintain good anticancer activity, and have low cytotoxicity. This provides a strategy about targeted and responsive systems of POD and VP-16 for clinical application.     

Valence isomerization of 2-phospha-4-silabicyclo[1.1.0]butane: a high-level ab initio study

The rearrangements for 2-phospha-4-silabicyclo[1.1.0]butane, analogous to the valence isomerization of the hydrocarbons bicyclobutane, 1,3-butadiene, and cyclobutene, were studied at the (U)QCISD(T)/6-311+G**//(U)QCISD/6-31G* level of theory. The monocyclic 1,2-dihydro-1,2-phosphasiletes are shown to be the thermodynamically preferred product, in contrast to the isomerization of the hydrocarbons, which favors the 1,3-butadiene structure. Furthermore, an unprecedented direct isomerization pathway to the 1,2-dihydro-1,2-phosphasiletes was identified. This pathway is competitive with the isomerization via the open-chain butadienes and becomes favorable when electron-donating substituents are present on silicon. Figure 2-Phospha-4-silabicyclo[1.1.0]butane can isomerize directly into the more stable P,Si-cyclobutene via an unprecedented [sigma2s+sigma2a] process, which becomes favorable over the isomerization via the P,Si-butadiene when electron-donating substituents are present on silicon.     

van der Waals forces between cylinders: Arrays of thin cylinders and three body forces

We study the van der Waals energy of interaction of an array of parallel dielectric cylinders immersed in a dielectric medium. We consider only “thin” cylinders which have radius small compared to the separation of the cylinders. The energy is calculated as a sum of two body plus three body interactions. The case of hexagonally close packed cylinders is studied in more detail. Some biophysical applications are discussed and in particular the van der Waals energy of the myosin lattice in striated muscle is examined.     

Ab initio and DFT study on the electrophilic addition of bromine to endo-tricyclo[3.2.1.02,4]oct-6-ene

Full geometric optimization of endo-tricyclo[3.2.1.0^2,4]oct-6-ene (endo-TCO) by ab initio and DFT methods allowed us to investigate the structure of the molecule. The double bond is endo-pyramidalized and its two faces are no longer found to be equivalent. The exo face of the double bond has regions with far more electron density (q_i,HOMO) and more negative electrostatic potential. The endo-TCO-Br₂ system was investigated at the B3LYP/6-311+G** level and the endo-TCO···Br₂(exo) molecular complex was found to be relatively more stable than the endo-TCO···Br₂(endo) complex. The cationic intermediates of the reaction were studied by ab initio and DFT methods. The bridged exo-bromonium cation(I) is relatively more stable than the endo-bromonium cation(II). An absolute exo-facial selectivity should be observed in the addition reaction of Br₂ to endo-TCO, which is caused by steric and electronic factors. The nonclassical rearranged cation IV was found to be the most stable ion among the cationic intermediates and the ionic addition occurs via the formation of this cation. The mechanism of the addition reaction is also discussed.     

Hydrogen bonding interactions in PN···HX complexes: DFT and ab initio studies of structure, properties and topology

Spin-restricted DFT (X3LYP and B3LYP) and ab initio (MP2(fc) and CCSD(fc)) calculations in conjunction with the Aug-CC-pVDZ and Aug-CC-pVTZ basis sets were performed on a series of hydrogen bonded complexes PN···HX (X = F, Cl, Br) to examine the variations of their equilibrium gas phase structures, energetic stabilities, electronic properties, and vibrational characteristics in their electronic ground states. In all cases the complexes were predicted to be stable with respect to the constituent monomers. The interaction energy (ΔE) calculated using a super-molecular model is found to be in this order: PN···HF > PN···HCl > PN···HBr in the series examined. Analysis of various physically meaningful contributions arising from the Kitaura-Morokuma (KM) and reduced variational space self-consistent-field (RVS-SCF) energy decomposition procedures shows that the electrostatic energy has significant contribution to the over-all interaction energy. Dipole moment enhancement (Δμ) was observed in these complexes expected of predominant dipole-dipole electrostatic interaction and was found to follow the trend PN···HF > PN···HCl > PN···HBr at the CCSD level. However, the DFT (X3LYP and B3LYP) and MP2 levels less accurately determined these values (in this order HF < HCl < HBr). Examination of the harmonic vibrational modes reveals that the PN and HX bands exhibit characteristic blue- and red shifts with concomitant bond contraction and elongation, respectively, on hydrogen bond formation. The topological or critical point (CP) analysis using the static quantum theory of atoms in molecules (QTAIM) of Bader was considered to classify and to gain further insight into the nature of interaction existing in the monomers PN and HX, and between them on H-bond formation. It is found from the analysis of the electron density ρ _c , the Laplacian of electron charge density ∇²ρ_c, and the total energy density (H _c ) at the critical points between the interatomic regions that the interaction N···H is indeed electrostatic in origin (ρ_c > 0, ∇²ρ_c > 0 and H_c > 0 at the BCP) whilst the bonds in PN (ρ_c > 0, ∇²ρ_c > 0 and H_c < 0) and HX ((ρ_c > 0, ∇²ρ_c < 0 and H_c < 0)) are predominantly covalent. A natural bond orbital (NBO) analysis of the second order perturbation energy lowering, E⁽²⁾, caused by charge transfer mechanism shows that the interaction N···H is n(N) → BD*(HX) delocalization.     

Linear free energy relationship for the anomeric effect: MP2, DFT and ab initio study of 2-substituted-1,4-dioxanes

The anomeric effect of 2-substituted 1,4-dioxane derivatives was calculated and compared with the values for substituted cyclohexane. The bond lengths, bond angles, torsion angles, and relative energies of axial and equatorial conformers of 2-substituted 1,4-dioxanes were calculated by the second-order Møller–Plesset (MP2), density functional theory (DFT/B3LYP), and Hartree–Fock (HF) methods using 6-31G^∗ basis set. The energy differences between the axial and equatorial conformers, endo and exo-anomeric effects, repulsive non-bond and H-bonding interactions were investigated. A linear free energy relationship (LFER) between calculated (MP2/6-31G^∗) anomeric effect and inductive substituent constants (σ_I) was obtained for 2-substituted-1,4-dioxanes (slope = 6.19 and r² = 0.967). The calculated energy differences indicate lower equatorial orientation for 2-substituted-1,4-dioxanes compared to the 2-substituted-tetrahydropyrans. The contribution of resonance, hyperconjugation, inductive, steric, hydrogen bonding, electrostatic interaction, and level of theory influences the anomeric effect.     

On the importance of van der Waals interaction in the groove binding of DNA with ligands: restrained molecular dynamics study

Comparison of interaction energy between an oligonucleotide and a DNA-binding ligand in the minor and major groove modes was made by use of restrained molecular dynamics. Distortion in DNA was found for the major groove mode whereas less significant changes for both ligand and DNA were detected for the minor groove binding after molecular dynamics simulation. The conformation of the ligand obtained from the major groove mode resembles that computed with the ligand soaked in water. The van der Waals contact energy was found to be as significant as electrostatic energy and more important for difference in binding energy between these two binding modes. The importance of van der Waals force in groove binding was supported by computations on the complex formed by the repressor peptide fragment from the bacteriophage 434 and its operator oligonucleotide.     

Importance of van der Waals contact between Glu 35 and Trp 109 to the catalytic action of human lysozyme.

The importance of van der Waals contact between Glu 35 and Trp 109 to the active-site structure and the catalytic properties of human lysozyme (HL) has been investigated by site-directed mutagenesis. The X-ray analysis of mutant HLs revealed that both the replacement of Glu 35 by Asp or Ala, and the replacement of Trp 109 by Phe or Ala resulted in a significant but localized change in the active-site cleft geometry. A prominent movement of the backbone structure was detected in the region of residues 110 to 120 and in the region of residues 100 to 115 for the mutations concerning Glu 35 and Trp 109, respectively. Accompanied by the displacement of the main-chain atoms with a maximal deviation of C alpha atom position ranging from 0.7 A to 1.0 A, the mutant HLs showed a remarkable change in the catalytic properties against Micrococcus luteus cell substrate as compared with native HL. Although the replacement of Glu 35 by Ala completely abolished the lytic activity, HL-Asp 35 mutant retained a weak but a certain lytic activity, showing the possible involvement of the side-chain carboxylate group of Asp 35 in the catalytic action. The kinetic consequence derived from the replacement of Trp 109 by Phe or Ala together with the result of the structural change suggested that the structural detail of the cleft lobe composed of the residues 100 to 115 centered at Ala 108 was responsible for the turnover in the reaction of HL against the bacterial cell wall substrate. The results revealed that the van der Waals contact between Glu 35 and Trp 109 was an essential determinant in the catalytic action of HL.     

H-bonded complexes between acetylacetone and two molecules of methanol: HF and DFT level study

Five stable H-bonded complexes (supersystems) between acetylacetone and two methanol molecules were investigated at the B3LYP and HF levels of theory using the 6-311G** and 6-11++G** basis sets. The most stable complex was found as the one with the highest relative bonding and interaction energies. All vibrational frequencies resulting from calculations with the 6-311++G** basis set were compared with the recorded IR spectrum of acetylacetone/methanol mixture in a molar ratio 1:2.     

Ab initio and DFT study of the inner mechanism and dynamic stereochemistry of electrophilic addition reaction of bromine to bisbenzotetracyclo[6.2.2.23,6.02,7]tetradeca-4,9,11,13-tetraene

The inner mechanism and dynamic stereochemistry of electrophilic addition of bromine to bisbenzotetracyclo[6.2.2.2(3,6).0(2,7)]tetradeca-4,9,11,13-tetraene(BBTT) molecule have been investigated by the methods of quantum chemistry. The structure of the BBTT molecule has been studied by ab initio and DFT/B3LYP methods using the 6-31G(d) and 6-311G(d) basis sets. The double bonds of BBTT molecule are endo-pyramidalized. The structure and stability of the cationic intermediates and products of the addition reaction have been investigated by HF/6-311G(d), HF/6-311G(d,p), B3LYP/6-311G(d) and B3LYP/6-311++G(2d,p)//B3LYP/6-311G(d) methods. The bridged bromonium cation isomerized into the more stable nonclassical delocalized N- and U-type cations and the difference between the stability of these cations is small. For the determination of the direction of addition reaction and the stereochemistry of the products, the stability of nonclassical delocalized N- and U-type ions and the structure of their cationic centres play a vital role. Since the cationic centre of the N-type ion is in interaction with the benzene ring from the exo face, the nucleofilic attackof the bromide anion to this centre occurs from the endo face and the exo,endo-isomer of the N-type product is obtained. The attack of bromide anion, towards the cationic centre of U-type ion from the endo face is sterically hindered by the hydrogen atom therefore the attack occurs from the exo face, which interacts with the benzene ring and the more stable exo,exo-isomer of U-type product is formed. Although, the U-type cation was 2.232 kcal mol(-1) more stable than the N-type cation, the U-type product was 0.587 kcal mol(-1) less stable than the N-type product.     

Effect of microsolvation on zwitterionic glycine: an ab initio and density functional theory study

The effect of microsolvation on zwitterionic glycine, considering both (-NH3(+)) as proton donor and (-COO(-)) as proton acceptor at correlated ab initio (MP2) level and density functional methods (B3LYP, PW91, MPW1PW91 and PBE) using 6-311++G** basis set has been reported. DFT methods have been employed so as to compare the performance/quality of different gradient-corrected correlation functionals (PW91, PBE), hybrid functionals (B3LYP, MPW1PW91) and to predict the near quantitative structural and vibrational properties, at reduced computational cost. B3LYP method outperforms among the different DFT methods for the computed hydrogen bond distances and found closer to the value obtained by correlated MP2 level, whereas MPW1PW91 and PBE methods shows very similar values but approximately 0.03 A less, compared to B3LYP method. MP2 calculation and single point CCSD(T)//MP2 calculation have been considered to decompose the interaction energy, including corrections for basis set superposition error (BSSE). Moreover, charge distribution analysis has also been carried out to understand the long raised questions, how and why the two body energies have significant contribution to the total binding energy.     

The interactions of square platinum(II) complexes with guanine and adenine: a quantum-chemical ab initio study of metalated tautomeric forms

The influence of binding of square planar platinum complexes on tautomeric equilibria of the DNA bases guanine and adenine was investigated using the density functional B3LYP method. Neutral trans-dichloro(amine)-, +1 charged chloro(diamine)-, and +2 charged triamine-platinum(II) species were chosen for coordination to bases. Only the N₇ interaction site of the bases was considered. The calculations demonstrate that the neutral platinum adduct does not change the tautomeric equilibria of the bases. Furthermore, N₇ binding of the neutral Pt adduct moderately reduces the probability of protonation of the N₁ position of adenine. Larger effects can be observed for +1 and mainly +2 adducts, but these can be rationalized by electrostatic effects. Since the electrostatic effects are expected to be efficiently compensated for by a charged backbone of DNA and counterions in a polar solvent, no dramatic increase in mispair formation is predicted for Pt(II) adducts, which is in agreement with experiment. The interaction energies between Pt adducts and the nucleobases were also evaluated. These interaction energies range from ca. 210 kJ/mol for neutral adducts, interacting with both bases and their tautomers, up to 500 kJ/mol for the +2 charged adducts, interacting with the keto-guanine tautomer and the anti-imino-adenine tautomer. The surprisingly large interaction energy for the latter structure is due to the strong H-bond between the NH₃ ligand group of the metal adduct and the N₆ nitrogen atom of the base. Received: 6 July 1999 / Accepted: 7 December 1999     

Theoretical spectroscopy of model-nonheme [Fe(IV)OL5]2+ complexes in their lowest triplet and quintet states using multireference ab initio and density functional theory methods

The structure, energies and spectroscopic properties of a simple [FeO(NH(3))(5)](2+) model with ground states (3)A(2g) and (5)A(1g) (in approximate C(4v) symmetry) have been studied in some detail using density functional (DFT) and simplified correlated multireference ab initio methods. The results reveal similarities as well as some pronounced differences in the properties of the molecule in the two alternative spin states.     

Actions between neonicotinoids and key residues of insect nAChR based on an ab initio quantum chemistry study: hydrogen bonding and cooperative pi-pi interaction

Neonicotinoid insecticides show selective actions on insect nicotinic acetylcholine receptor (nAChR). Two key residues (Trp and Arg/Lys) have been identified as contributing to the neonicotinois binding. To investigate the selective mechanism, a computational model was set up to simulate the interaction between residues (Trp and Arg) of insect nAChR and neonicotinoids by quantum chemistry method. Three analogues of neonicotinoid derivatives without the chloropyridinyl moiety and 3-methyl-indole (3MI), guanidinium (Gua) were used to mimic the neonicotinoids and the side chain of key residues Trp and Arg accordingly. Interaction features of 3MI-analogues, analogues-Gua and 3MI-analogues -Gua complexes were analyzed comparatively. Hydrogen bonding between the nitro group of analogues and Gua was found to be the most important for binding. Moreover, the cooperative pi-pi interaction between analogues and the indole ring, which is strengthened by the existence of Gua, also contributes to the binding. The alternative binding model of neonicotinoids proposed here, although slightly different from others, might be close to the actual.     

Theoretical ab initio study of the effects of methylation on structure and stability of G:C Watson-Crick base pair

Methylation of DNA occurs most readily at N(3), N(7), and O(6) of purine bases and N(3) and O(2) of pyrimidines. Methylated bases are continuously formed through endogenous and exogenous mechanisms. The results of a theoretical ab initio study on the methylation of G:C base pair components are reported. The geometries of the local minima were optimized without symmetry restrictions by the gradient procedure at DFT level of theory and were verified by energy second derivative calculations. The standard 6-31G(d) basis set was used. The single-point calculations have been performed at the MP2/6-31G(d,p), MP2/6-31++G(d,p), and MP2/6-311++G(2d,2p) levels of theory. The geometrical parameters, relative stability and counterpoise corrected interaction energies are reported. Also, using a variation-perturbation energy decomposition scheme we have found the vital contributions to the total interaction energy.     

High-level ab initio calculations on the energetics of low-lying spin states of biologically relevant transition metal complexes: a first progress report

Although DFT is the unrivaled method of choice for quantum chemical studies of bioinorganic problems, little is known about its ability to predict the energetics of the low-lying electronic states of transition metal complexes. The first high-level ab initio calculations aimed at calibrating DFT vis-a-vis this issue indicate that, despite its many successes, DFT is far from infallible. In the short term, additional calibration of DFT against more elaborate ab initio methods remains an important goal for computational bioinorganic researchers. In the longer term, we are optimistic that high-level ab initio methods such as CASPT2 and CCSD(T) will be regularly used to study realistic molecules of genuine biochemical interest.     

Spectrophotometric study of the interaction between rhenium complexes and phosphatidylcholine during liposome formation

The electron absorption spectra of halogenides and halogencarboxylate complex compounds of rhenium (III) having cluster structure with phosphatydilcholine and their lyposome forms were investigated. Some results which evidence for the interaction of these compounds with phosphatydilcholine were obtained. The possible mechanism of this interaction is discussed.     

DFT and MP2 study on the electrophilic addition reaction of bromine to exo-tricyclo[3.2.1.02.4]oct-6-ene

The geometry and electronic structure of exo-tricyclo[3.2.1.0^2,4]oct-6-ene (exo-TCO) was investigated using DFT methods. The two faces of the endo-pyramidalised double bond of the molecule are not equivalent. The exo face of the double bond has regions with high electron density (q _i,HOMO) and greater negative potential. Molecular complexes of exo-TCO with bromine were investigated using the B3LYP/6-311++G(d,p) method; the exo-TCO . . . Br₂(exo) molecular complex was found to be relatively more stable than the exo-TCO . . . Br₂(endo) complex. The cationic intermediates of the reaction were studied by DFT and MP2 methods. The solvent effect was evaluated using the self-consistent isodensity polarised continuum model (SCI-PCM). The exo-bromonium cation was found to be more stable than the endo-bromonium cation. Exo-facial selectivity due to electronic and steric factors was observed upon addition of bromine to exo-TCO. The multicentre nonclassical delocalised bromocarbonium cation IV and the exo-bridged-bromonium cation I are more stable than the rearrangement cation V. The reaction products are formed via exo-bridged-bromonium I and nonclassical IV cations, which are the most stable intermediates and whose stabilities barely differ. The mechanism of the addition reaction is also discussed.     

Conformational flexibility of the catalytic Asp dyad in HIV-1 protease: An ab initio study on the free enzyme

The enzyme protease from the human immunodeficiency virus type 1 (HIV-1 PR) is one of the main targets for therapeutic intervention in AIDS. Computer modeling is useful for probing the binding of novel ligands, yet empirical force field-based methods have encountered problems in adequately describing interactions of the catalytic aspartyl pair. In this work we use ab initio dynamic methods to study the molecular interactions and the conformational flexibility of the Asp dyad in the free enzyme. Calculations are performed on model complexes that include, besides the Asp dyad, the conserved Thr26 and Gly27 residues and water molecules present in the active site channel. Our calculations provide proton location and binding mode of the active-site water molecule, which turn out to be different from those of the eukariotic isoenzyme. Furthermore, the calculations reproduce well the structural features of the aspartyl dyad in the protein. Finally, they allow the identification of both dipole/charge interactions and a low-barrier hydrogen bond as important stabilizing factors for the peculiar conformation of the active site. These findings are consistent with site-directed mutagenesis experiments on the 27, 27; positions (Bagossi et al., Protein Eng 1996;9:997-1003). The electric field of the protein frame (included in some of the calculations) does not affect significantly the chemical bonding at the cleavage site. Proteins 2000;39:26-36.