Characterization of the structures and dynamics of phosphoric acid doped benzimidazole mixtures: a molecular dynamics study

Benzimidazole-based polymer membranes like poly(2,5-benzimidazole) (ABPBI) doped with phosphoric acid (PA) are electrolytes that exhibit high proton conductivity in fuel cells at elevated temperatures. The benzimidazole (BI) moiety is an important constituent of these membranes, so the present work was performed in order to achieve a molecular understanding of the BI–PA interactions in the presence of varying levels of the PA dopant, using classical molecular dynamics (MD) simulations. The various hydrogen-bonding interactions, as characterized based on structural properties and hydrogen-bond lifetime calculations, show that both BI and PA molecules exhibit dual proton-acceptor/donor functionality. An examination of diffusion coefficients showed that the diffusion of BI decreases with increasing PA uptake, whereas the diffusion of PA slightly increases. The hydrogen-bond lifetime calculations pointed to the existence of competitive hydrogen bonding between various sites in BI and PA.

Coarse-grained simulations for organic molecular liquids based on Gay-Berne and electric multipole potentials

Coarse-grained studies of CH₃SH, CH₃CHO and CHCl₃ liquids, based on anisotropic Gay-Berne (GB) and electric multipole potentials (EMP), demonstrate that the coarse-grained model is able to qualitatively reproduce the results obtained from the atomistic model (AMOEBA polarizable force field) and allows for significant saving in computation time. It should be pointed out that the accuracy of the coarse-grained model is very sensitive to how well the anisotropic GB particle is defined and how satisfactorily the EMP sites are chosen.

Conformational entropy of a polymer chain grafted to rough surfaces

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.

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.

Computational investigation on redox-switchable nonlinear optical properties of a series of polycyclic p-quinodimethane molecules

The polycyclic p-quinodimethanes are proposed to be the novel candidates of the high-performance nonlinear optical (NLO) materials because of their large third order polarizabilities (γ). We investigate the switchable NLO responses of a series of polycyclic p-quinodimethanes with redox properties by employing the density functional theory (DFT). The polycyclic p-quinodimethanes are forecasted to exhibit obvious pure diradical characters because of their large y ₀ index (the y ₀ index is a value between 0 [closed-shell state] and 1 [pure biradical state]). The γ values of these polycyclic p-quinodimethanes and their corresponding one-electron and two-electron reduced/oxidized species are calculated by the (U)BHandHLYP method. The γ values of polycyclic p-quinodimethanes and their corresponding one-electron reduced species are all positive and significantly different. The large differences of the γ values are due to a change in the transition energy and are related to the different delocalization of the spin density, which demonstrates that the NLO switching is more effective on one-electron reduction reactions. Therefore, the study on these polycyclic p-quinodimethanes provides a guideline for a molecular design of highly efficient NLO switching.

In pursuit of negative Fukui functions: examples where the highest occupied molecular orbital fails to dominate the chemical reactivity

In our quest to explore molecules with chemically significant regions where the Fukui function is negative, we explored reactions where the frontier orbital that indicates the sites for electrophilic attack is not the highest occupied molecular orbital. The highest occupied molecular orbital (HOMO) controls the location of the regions where the Fukui function is negative, supporting the postulate that negative values of the Fukui function are associated with orbital relaxation effects and nodal surfaces of the frontier orbitals. Significant negative values for the condensed Fukui function, however, were not observed.

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.

Cubic Magnetically Guided Nanoaggregates for Inhalable Drug Delivery: In Vitro Magnetic Aerosol Deposition Study

The present work describes the in vitro aerosol deposition and enhanced deaggregation behavior of superparamagnetic iron oxide nanoaggregates (SPIONs). SPIONs were surface-coated with amine functionalized polyrotaxane and were proposed as a carrier for inhalation dry powders. Polyrotaxane is primarily composed of beta cyclodextrin rings which are spontaneously threaded on the block copolymer, poly(propylene glycol) bis(2-aminopropylether). Variable concentrations of surface coating polymers showed controlled manipulation of the crystal size and morphology. Magnetic nanoaggregates fabricated with low concentration of polyrotaxane showed cubic crystal morphology. However, these nanoaggregates exhibited rhombic dodecahedron crystal structure upon increasing the coating polymer concentration. In comparison to the spherical uncoated magnetic nanoparticles, cubic phase magnetic nanoaggregates demonstrated an enhanced in vitro aerosol deposition using magnetic field alignment. This enhancement can be accomplished at low inhalation flow rates (15 and 30 L/min). However, transformation to the cubic crystal structure was observed to be associated with a reduction in the powder geometric standard deviation. Using a mathematical modeling approach, we noted significant enhancement in the deaggregation behavior of inhalation dry powders; that can be achieved with small amounts of magnetic nanoaggregates. Aggregates of cubic nanoparticles showed promise for targeted pulmonary deposition of anticancer drugs.

Metal ions in sugar binding, sugar specificity and structural stability of Spatholobus parviflorus seed lectin

Spatholobus parviflorus seed lectin (SPL) is a heterotetrameric lectin, with two α and two β monomers. In the crystal structure of SPL α monomer, two residues at positions 240 and 241 are missing. This region was modeled based on the positional and sequence similarities. The role of metal ions in SPL structure was analyzed by 10 ns molecular dynamics simulation. MD simulations were performed in the presence and absence of metal ions to explain the loss of haemagglutinating property of the lectin due to demetallization. Demetallized structure was found to deviate drastically at the metal binding loop region. Affinity of different sugars like N-acetyl galactosamine (GalNAc), D-galactose and lactose towards the native and demetallized protein was calculated by molecular docking studies. It was found that the sugar binding site got severely distorted in demetallized lectin. Consequently, sugar binding ability of lectin might be decreasing in the demetallized condition. Isothermal titration calorimetric (ITC) analysis of the sugars in the presence of native and demetallized protein confirmed the in silico results. It was observed after molecular dynamics simulations, that significant structural deviations were not caused in the quaternary structure of demetallized lectin. It was confirmed that the structural changes modified the sugar binding ability, as well as sugar specificity of the present lectin. The role of metal ions in sugar binding is described based on the in silico studies and ITC analysis. A comprehensive analysis of the ITC data suggests that the sugar specificity of the metal bound lectin and the loss of sugar specificity due to metal chelation are not linear.

Electric field effect on the zigzag (6,0) single-wall BC2N nanotube for use in nano-electronic circuits

We have analyzed the effect of external electric field on the zigzag (6,0) single-wall BC₂N nanotube using density functional theory calculations. Analysis of the structural parameters indicates that the nanotube is resistant against the external electric field strengths. Analysis of the electronic structure of the nanotube indicates that the applied parallel electric field strengths have a much stronger interaction with the nanotube with respect to the transverse electric field strengths and the nanotube is easier to modulate by the applied parallel electric field. Our results show that the properties of the nanotube can be controlled by the proper external electric field for use in nano-electronic circuits.

Quantitative analysis of molecular surfaces: areas, volumes, electrostatic potentials and average local ionization energies

We describe a procedure for performing quantitative analyses of fields f(r) on molecular surfaces, including statistical quantities and locating and evaluating their local extrema. Our approach avoids the need for explicit mathematical representation of the surface and can be implemented easily in existing graphical software, as it is based on the very popular representation of a surface as collection of polygons. We discuss applications involving the volumes, surface areas and molecular surface electrostatic potentials, and local ionization energies of a group of 11 molecules.

Cubitane: a rare diterpenoid skeleton

There are only 15 diterpenoids known sharing the rare cubitane skeleton (1,3-diisopropyl-6,10-dimethylcyclododecane, Fig. 1) named after the natural product (+)-cubitene ((+)-1) from the termite Cubitermes umbratus. Cubitane-type diterpenoids (“cubitanoids”) have since then been isolated from gorgonian corals and there are reports that cubitene also occurs in plants. The twelve-membered ring of the cubitanoids is an interesting feature for which we have developed a novel synthesis involving bicyclic precursors. In this review, we discuss the state-of-research regarding isolation, structure elucidation, biological activity, biosynthesis, and total synthesis of cubitane-type diterpenoids.

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.

Cyclo-hexa-peptides at the water/cyclohexane interface: a molecular dynamics simulation

Molecular dynamic (MD) simulations have been performed to study the behaviors of ten kinds of cyclo-hexa-peptides (CHPs) composed of amino acids with the diverse hydrophilic/hydrophobic side chains at the water/cyclohexane interface. All the CHPs take the “horse-saddle” conformations at the interface and the hydrophilicity/hydrophobicity of the side chains influences the backbones’ structural deformations. The orientations and distributions of the CHPs at the interface and the differences of interaction energies (ΔΔE) between the CHPs and the two liquid phases have been determined. RDF analysis shows that the H-bonds were formed between the O_C atoms of the CHPs’ backbones and H_w atoms of water molecules. N atoms of the CHPs’ backbones formed the H-bonds or van der Waals interactions with the water solvent. It was found that there is a parallel relationship between ΔΔE and the lateral diffusion coefficients (D _xy ) of the CHPs at the interface. The movements of water molecules close to the interface are confined to some extent, indicating that the dynamics of the CHPs and interfacial water molecules are strongly coupled.

Tunable differential conductance of single wall C/BN nanotube heterostructure

The transport properties and differential conductance of the heterostructures constructed by (5,5) single wall carbon nanotube (SWCNT) and (5,5) single wall boron nitride nanotube (SWBNNT) are investigated using density functional theory in combination with non-equilibrium Green’s functions. We find that the transmission conductance of (5,5) BN/C nanotube heterostructure is not only continually depressed as the BNNT region increases but also the drop of the conductance is uniform in the energy window (?1.43 eV, 1 eV), which leads to linear I–V dependence for the systems when the bias is within this energy range. Moreover, the differential conductance linearly decreases when n?≤?3 but exponentially decreases when n?≥?3 for (5,5)(BN) _n /C heterostructure. Such tunable differential conductance of (5,5) BN/C nanotube heterostructure mainly derives from the blockage of the transport channels induced by the semiconductive BN segment.

Effect of pristine and functionalized single- and multi-walled carbon nanotubes on CO<Subscript>2</Subscript> separation of mixed matrix membranes based on polymers of intrinsic microporosity (PIM-1): a molecular dynamics simulation study

Molecular dynamics (MD) and grand canonical Monte Carlo (GCMC) simulations were conducted to investigate the transport properties of carbon dioxide, methane, nitrogen, and oxygen through pure and mixed matrix membranes (MMMs) based on polymers of intrinsic microporosity (PIM-1). For this purpose, first, 0.5 to 3 wt% of pristine single-walled carbon nanotube (p-SWCNT) and multi-walled carbon nanotube (p-MWCNT) were embedded into the pure PIM-1, and then for better dispersion of CNT particles into the polymer matrix and to improve the performance of the resulting MMMs, polyethylene glycol (PEG) functionalized SWCNT and MWCNT (f-SWCNT and f-MWCNT, respectively) were loaded. The characterization of the obtained MMMs was carried out by using density, glass transition temperature, X-ray pattern, and fractional free volume calculations. Comparing the obtained results with the available reported experimental data, indicate the authenticity of the applied simulation approach. The simulation results exhibit that the pristine and PEG-functionalized CNT particles improve the transport properties such as diffusivity, solubility, and permeability of the PIM-1 membranes, without sacrificing their selectivity. Also, the MMMs incorporated with 2 wt% of the functionalized CNT particles indicate better performance for the CO₂ separation from other gases. According to the calculated results, the highest permeability and diffusivity for CO₂ are observed in the [PIM-1/f-SWCNT] MMM among the other membranes which represent that the loading of the f-SWCNTs can enhance the CO₂ separation performance of PIM-1 more than other CNTs studied in this work.     

Simulations on the possibility of formation of complexes between fluorouracil drug and cucurbit[n]urils: ab initio van der Waals DFT study

The binding geometry of fluorouracil/cucurbit[n]urils (CB[n]s) complexes with n?=?5–8 is investigated using the first-principles van der Waals density functional (vdW-DF) method, involving full geometry optimization. Such host-guest complexes are typically calculated using conventional DFT method, which significantly underestimates non-local dispersion forces (or vdW contributions) and therefore affects interactions between respected entities. We address here the role of vdW forces for the fluorouracil and CB[n]s molecules which can form directional hydrogen bonds with each other. It was found that the inclusion of dispersion interactions significantly affects the host-guest binding properties and the hydrogen bonding between the molecules provides the main binding mechanism, while results in the same geometries for the considered complexes. The 0.84 eV binding energy, for the thermodynamically favorable state, reveals that the interaction of fluorouracil with CB[n]s is an exothermic interaction and typical for strong hydrogen bonding. Furthermore, we have investigated the binding nature of these host-guest systems in aqueous solution with ab initio MD simulations adopting vdW-DF method. These findings afford evidence for the applicability of the vdW-DF approach and provide a realistic benchmark for the investigation of the host-guest complexes.

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.

Entropy versus aromaticity in the conformational dynamics of aromatic rings

Comparison of the results of Car-Parrinello molecular dynamics simulations of isolated benzene, pyrimidine and 1,2,4-triazine molecules reveals that the unusually low population of planar geometry of the benzene ring is caused by entropy effects despite its high aromaticity. The decrease in symmetry of the molecule results in smaller changes in entropy and Gibbs free energy due to out-of-plane deformations of the ring, leading to an increase in the population of planar geometry of the ring. This leads to differences in the topology of potential energy and Gibbs free energy surfaces.