A comparative study on carbon, boron-nitride, boron-phosphide and silicon-carbide nanotubes based on surface electrostatic potentials and average local ionization energies

A density functional theory study was carried out to predict the electrostatic potentials as well as average local ionization energies on both the outer and the inner surfaces of carbon, boron-nitride (BN), boron-phosphide (BP) and silicon-carbide (SiC) single-walled nanotubes. For each nanotube, the effect of tube radius on the surface potentials and calculated average local ionization energies was investigated. It is found that SiC and BN nanotubes have much stronger and more variable surface potentials than do carbon and BP nanotubes. For the SiC, BN and BP nanotubes, there are characteristic patterns of positive and negative sites on the outer lateral surfaces. On the other hand, a general feature of all of the systems studied is that stronger potentials are associated with regions of higher curvature. According to the evaluated surface electrostatic potentials, it is concluded that, for the narrowest tubes, the water solubility of BN tubes is slightly greater than that of SiC followed by carbon and BP nanotubes.

Exploration of various electronic properties along the reaction coordinate for hydration of Pt(II) and Ru(II) complexes; the CCSD, MPx, and DFT computational study

In the study behavior of molecular electrostatic potential, averaged local ionization energy, and reaction electronic flux along the reaction coordinate of hydration process of three representative Ru(II) and Pt(II) complexes were explored using both post-HF and DFT quantum chemical approximations. Previously determined reaction mechanisms were explored by more detailed insight into changes of electronic properties using ωB97XD functional and MP2 method with 6–311++G(2df,2pd) basis set and CCSD/6–31(+)G(d,p) approach. The dependences of all examined properties on reaction coordinate give more detailed understanding of the hydration process.

Modeling the effect of H-bonding interactions and molecular packing on the molecular structure of [Ag(ethylnicotinate)2]NO3 complex

The gas phase molecular structure of a single isolated molecule of [Ag(Etnic)₂NO₃];1 where Etnic = Ethylnicotinate was calculated using B3LYP method. The H-bonding interaction between 1 with one (complex 2) and two (complex 3) water molecules together with the dimeric formula [Ag(Etnic)₂NO₃]₂;4 and the tetrameric formula [Ag(Etnic)₂NO₃]₄;5 were calculated using the same level of theory to model the effect of intermolecular interactions and molecular packing on the molecular structure of the titled complex. The H-bond dissociation energies of complexes 2 and 3 were calculated to be in the range of 12.220–14.253 and 30.106–31.055 kcal?mol^?1, respectively, indicating the formation of relatively strong H-bonds between 1 and water molecules. The calculations predict bidentate nitrate ligand in the case of 1 and 2, leading to distorted tetrahedral geometry around the silver ion with longer Ag–O distances in case of 2 compared to 1, while 3 has a unidentate nitrate ligand leading to a distorted trigonal planar geometry. The packing of two [Ag(Etnic)₂NO₃] complex units; 4 does not affect the molecular geometry around Ag(I) ion compared to 1. In the case of 5, the two asymmetric units of the formula [Ag(Etnic)₂NO₃] differ in the bonding mode of the nitrate group, where the geometry around the silver ion is distorted tetrahedral in one unit and trigonal planar in the other. The calculations predicted almost no change in the charge densities at the different atomic sites except at the sites involved in the C–H?O interactions as well as at the coordinated nitrogen of the pyridine ring.

Density functional studies on photophysical properties and chemical reactivities of the triarylboranes: effect of the constraint of planarity

The geometric and electronic structures, absorption spectra, transporting properties, chemical reactivity indices and electrostatic potentials of the planar three-coordinate organoboron compounds 1-2 and twisted reference compound Mes ₃ B, have been investigated by employing density functional theory (DFT) and conceptual DFT methods to shed light on the planarity effects on the photophysical properties and the chemical reactivity. The results show that the planar compounds 1-2 exhibit significantly lower HOMO level than Mes ₃ B, owing to the stronger electronic induction effect of boron centers. This feature conspicuously induces a blue shifted absorption for 1, although 1 seemingly possesses more extended conjugation framework than Mes ₃ B. Importantly, the reactivity strength of the boron atoms in 1-2 is much lower than that in Mes ₃ B, despite the fact that the tri-coordinate boron centers of 1-2 are completely naked. The interesting and abnormal phenomenon is caused by the strong p-π electronic interactions, that is, the empty p-orbital of boron center is partly filled by π-electron of the neighbor carbon atoms in 1-2, which are confirmed by the analysis of Laplacian of the electron density and natural bond orbitals. Furthermore, the negative electrostatic potentials of the boron centers in 1-2 also interpret that they are not the most preferred sites for incoming nucleophiles. Moreover, it is also found that the planar compounds 1-2 can act as promising electron transporting materials since the internal reorganization energies for electron are really small.

Effects of protein binding on a lipid bilayer containing local anesthetic articaine, and the potential of mean force calculation: a molecular dynamics simulation approach

Articaine, as a local anesthetic drug has been simulated in neutral and charged forms, and its interaction with the dimyristoylphosphatidylcholine (DMPC) lipid bilayer membrane is investigated by molecular dynamics simulation using GROMACS software. In order to obtain the optimum location of the drug molecules, as they penetrate into the membrane, umbrella sampling is applied and the free energy is calculated. The effect of protein binding to DMPC membrane on the process of drug diffusion through the membrane is considered. Five simulation systems are designed and by applying the potential of mean force, the molecular dynamics simulation on the system is performed. In light of the obtained results, the electrostatic potential, variation of lipid bilayer’s order parameter and the diffusion coefficient of drug are discussed.

Perspectives on the reaction force constant

A synchronous, concerted chemical process is rigorously divided by the reaction force F(R), the negative gradient of V(R), into “reactant” and “product” regions which are dominated by structural changes and an intervening “transition” region which is electronically intensive. The reaction force constant κ(R), the second derivative of V(R), is negative throughout the transition region, not just at the nominal transition state, at which κ(R) has a minimum. This is consistent with experimental evidence that there is a transition region, not simply a specific point. We show graphically that significant nonsynchronicity in the process is associated with the development of a maximum of κ(R) in the transition region, which increases as the process becomes more nonsynchronous. (We speculate that for a nonconcerted process this maximum is actually positive.) Thus, κ(R) can serve as an indicator of the level of nonsynchronicity.

Quantitative characterization of protein tertiary motifs

Dimerization of 2-naphthalenecarbonitrile (2-NpCN) mediated by cucurbit[8]uril (CB[8]) has been investigated employing the density functional theory. Different structures of 2-NpCN dimers were generated by combining monomers in anti-head-to-head (A), anti-head-to-tail (B) and syn-head-to-tail (C) fashion. All these dimeric structures possess rigid cube-like architecture. On confinement within the CB[8] dimer A turns out to be the lowest energy structure. Calculated ¹H NMR spectra revealed that the 2-NpCN dimer exhibits large shielding for aromatic protons consistent with the experiment. The protons attached to cubane moiety on the other hand, led to down-field signals. Dimerization mediated with CB[8] cavitand is further accompanied by the frequency up-shift (blue shift) of methylene stretching vibration in its infrared spectra.

Cytotoxic effect and molecular docking of 4-ethoxycarbonylmethyl-1-(piperidin-4-ylcarbonyl)-thiosemicarbazide—a novel topoisomerase II inhibitor

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 IC₅₀ 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.

Electron correlation effects and density analysis of the first-order hyperpolarizability of neutral guanine tautomers

Dipole moments (μ), charge distributions, and static electronic first-order hyperpolarizabilities (β _μ ) of the two lowest-energy keto tautomers of guanine (7H and 9H) were determined in the gas phase using Hartree–Fock, Møller–Plesset perturbation theory (MP2 and MP4), and DFT (PBE1PBE, B97-1, B3LYP, CAM-B3LYP) methods with Dunning’s correlation-consistent aug-cc-pVDZ and d-aug-cc-pVDZ basis sets. The most stable isomer 7H exhibits a μ value smaller than that of the 9H form by a factor of ca. 3.5. The β _μ value of the 9H tautomer is strongly dependent on the computational method employed, as it dramatically influences the β _μ (9H)/β _μ (7H) ratio, which at the highest correlated MP4/aug-cc-pVDZ level is predicted to be ca. 5. The Coulomb-attenuating hybrid exchange-correlation CAM-B3LYP method is superior to the conventional PBE1PBE, B3LYP, and B97-1 functionals in predicting the β _μ values. Differences between the largest diagonal hyperpolarizability components were clarified through hyperpolarizability density analyses. Dipole moment and first-order hyperpolarizability are molecular properties that are potentially useful for distinguishing the 7H from the 9H tautomer.

Theoretical evaluation of flotation performance of carboxyl hydroxamic acids with different number of polar groups on the surfaces of diaspore (010) and kaolinite (001)

The adsorption behaviors of three carboxyl hydroxamic acids on diaspore (010) and kaolinite (001) have been studied by density functional theory (DFT) and molecular dynamics (MD) method. The results indicated that carboxyl hydroxamic acids could adsorb on diaspore surface by ionic bonds and hydrogen bonds, and adsorb on kaolinite surface by hydrogen bonds. The models of carboxyl hydroxamic acids adsorbed on diaspore and kaolinite surfaces are proposed.

Computationally designed prodrugs of statins based on Kirby’s enzyme model

DFT calculations at B3LYP/6-31G(d,p) for intramolecular proton transfer in Kirby’s enzyme models 1–7 demonstrated that the reaction rate is dependent on the distance between the two reacting centers, r_GM, and the hydrogen bonding angle, α, and the rate of the reaction is linearly correlated with r_GM and α. Based on these calculation results three simvastatin prodrugs were designed with the potential to provide simvastatin with higher bioavailability. For example, based on the calculated log EM for the three proposed prodrugs, the interconversion of simvastatin prodrug ProD 3 to simvastatin is predicted to be about 10 times faster than that of either simvastatin prodrug ProD 1 or simvastatin ProD 2. Hence, the rate by which the prodrug releases the statin drug can be determined according to the structural features of the promoiety (Kirby’s enzyme model).

Prodrugs for masking bitter taste of antibacterial drugs—a computational approach

DFT calculations for the acid-catalyzed hydrolysis of several maleamic acid amide derivatives revealed that the reaction rate-limiting step is determined on the nature of the amine leaving group. Further, it was established that when the amine leaving group was a secondary amine, acyclovir or cefuroxime moiety the tetrahedral intermediate formation was the rate-limiting step such as in the cases of acyclovir ProD 1- ProD 4 and cefuroxime ProD 1- ProD 4. In addition, the linear correlation between the calculated and experimental rates provided a credible basis for designing prodrugs for masking bitter taste of the corresponding parental drugs which have the potential to release the parent drug in a sustained release fashion. For example, based on the DFT calculated rates the predicted t_1/2 (a time needed for 50 % of the reactant to be hydrolyzed to products) for cefuroxime prodrugs, cefuroxime ProD 1- ProD 4, were 12 min, 18 min, 200 min and 123 min, respectively.

Prodrugs of fumarate esters for the treatment of psoriasis and multiple sclerosis—a computational approach

Density functional theory (DFT) calculations at B3LYP/6-31 G (d,p) and B3LYP/6-311?+?G(d,p) levels for the substituted pyridine-catalyzed isomerization of monomethyl maleate revealed that isomerization proceeds via four steps, with the rate-limiting step being proton transfer from the substituted pyridinium ion to the C=C double bond in INT1. In addition, it was found that the isomerization rate (maleate to fumarate) is solvent dependent. Polar solvents, such as water, tend to accelerate the isomerization rate, whereas apolar solvents, such as chloroform, act to slow down the reaction. A linear correlation was obtained between the isomerization activation energy and the dielectric constant of the solvent. Furthermore, linearity was achieved when the activation energy was plotted against the pK _a value of the catalyst. Substituted-pyridine derivatives with high pK _a values were able to catalyze isomerization more efficiently than those with low pK _a values. The calculated relative rates for prodrugs 1–6 were: 1 (406.7), 2 (7.6?×?10⁶), 3 (1.0), 4 (20.7), 5 (13.5) and 6 (2.2?×?10³). This result indicates that isomerizations of prodrugs 1 and 3–5 are expected to be slow and that of prodrugs 2 and 6 are expected to be relatively fast. Hence, prodrugs 2 and 3–5 have the potential to be utilized as prodrugs for the slow release of monomethylfumarate in the treatment of psoriasis and multiple sclerosis.

Halogen bond tunability II: the varying roles of electrostatic and dispersion contributions to attraction in halogen bonds

In a previous study we investigated the effects of aromatic fluorine substitution on the strengths of the halogen bonds in halobenzene…acetone complexes (halo?=?chloro, bromo, and iodo). In this work, we have examined the origins of these halogen bonds (excluding the iodo systems), more specifically, the relative contributions of electrostatic and dispersion forces in these interactions and how these contributions change when halogen σ-holes are modified. These studies have been carried out using density functional symmetry adapted perturbation theory (DFT-SAPT) and through analyses of intermolecular correlation energies and molecular electrostatic potentials. It is found that electrostatic and dispersion contributions to attraction in halogen bonds vary from complex to complex, but are generally quite similar in magnitude. Not surprisingly, increasing the size and positive nature of a halogen’s σ-hole dramatically enhances the strength of the electrostatic component of the halogen bonding interaction. Not so obviously, halogens with larger, more positive σ-holes tend to exhibit weaker dispersion interactions, which is attributable to the lower local polarizabilities of the larger σ-holes.

Molecular simulation study of PAMAM dendrimer composite membranes

Pure polysulfone (PSF) and its composites with chitosan (CST), hyaluronic acid (HA), conventional poly(amidoamine), and hydroxyl poly(amidoamine) dendrimers as the membranes for separation of the gases, methane, carbon dioxide, hydrogen sulfide, nitrogen, and oxygen have been studied by molecular dynamics (MD) and grand canonical Monte Carlo (GCMC) simulations. The transport properties (solubility, diffusivity, and permeability) of pure and gas mixtures in the membranes were calculated and the results of the simulations were compared with the available experimental data. The simulated structural properties of the pure and composite PSF membranes including occupied volume, free volume, surface area, fractional free volume (FFV), and radius of gyration (R _g ) were evaluated and their effects on the separability of the gases by the membranes were analyzed and interpreted by the obtained results.

Investigation of the binding network of IGF-I on the cavity surface of IGFBP4

Insulin-like growth factor-binding proteins (IGFBPs) control bioactivity and distribution of insulin-like growth factors (IGFs) through high-affinity complex of IGFBP and IGF. To get more insight into the binding interaction of IGF system, the site-directed mutagenesis and force-driving desorption methods were employed to study the interaction mechanism of IGFBP4 and IGF-I by molecular dynamics (MD) simulation. In IGF-I, residues Gly7 to Asp12 were found to be the hot spots and they mainly anchored on the N-domain of IGFBP4. The contact area, the shape and size of protein, the surroundings of the binding site, the hydrophobic and electrostatic interaction between the two proteins worked as a complex network to regulate the protein-protein interaction. It was also found that the unfolding of the helix was not inevitable in the mutant, and it could be regulated by careful selection of the substituted amino acid.