Herein we report a study of the switchable [3]rotaxane reported by Huang et al. (Appl Phys Lett 85(22):5391–5393, 1) that can be mounted to a surface to form a nanomechanical, linear, molecular motor. We demonstrate the application of semiempirical electronic structure theory to predict the average and instantaneous force generated by redox-induced ring shuttling. Detailed analysis of the geometric and electronic structure of the system reveals technical considerations essential to success of the approach. The force is found to be in the 100–200 pN range, consistent with published experimental estimates.
Development of new energetic salts is the key factor in replacing low performance compounds in conventional formulations of high explosives as well as propellants. Ten salts based on the nitroformate anion and various nitrogen-rich cations were designed and their geometric optimizations carried out using the density functional method. With reasonable oxygen balance (from ?36 % to 0 %), heats of formation (47–624 kJ mol?1) and high densities (1.81–1.89 g cm?3), the detonation velocity (D) and pressure (P) values of salts were calculated as 8.62–9.36 km s?1 and 33.10–40.01 GPa, respectively. Lastly, the nitroformate salts studied in this work are of prospective interest as high performance explosives.
The activation of human epidermal growth factor receptor (hEGFR) involves a large conformational change in its soluble extracellular domains (sECD, residues 1–620), from a tethered to an extended conformation upon binding of ligands, such as EGF. It has been reported that this dynamic process is pH-dependent, that is, hEGFR can be activated by EGF at high pH to form an extended dimer but remains as an inactive monomer at low pH. In this paper, we perform all-atom molecular dynamics (MD) simulations starting from the tethered conformation of sECD:EGF complex, at pH 5.0 and 8.5, respectively. Simulation results indicate that sECD:EGF shows different dynamic properties between the two pHs, and the complex may have a higher tendency of activation at pH 8.5. Twenty residues, including 13 histidines, in sECD:EGF have different protonation states between the two pHs (calculated by the H++ server). The charge distribution at pH 8.5 is more favorable for forming an extended conformation toward the active state of sECD than that at pH 5.0. Our study may shed light on the mechanism of pH dependence of hEGFR activation.
The chelating properties of diethanoldithiocarbamate (DEDC) and π-electron flow from the nitrogen atom to the sulfur atom via a plane-delocalized π-orbital system (quasi ring) was studied using a density functional theory method. The molecular structure of DEDC and its complexes with Zn(II), Cd(II), and Hg(II) were also considered. First, the geometries of this ligand and DEDC-Zn(II), DEDC-Cd(II), and DEDC-Hg(II) were optimized, and the formation energies of these complexes were then calculated based on the electronic energy, or sum of electronic energies, with the zero point energy of each species. Formation energies indicated the DEDC-Zn(II) complex as the most stable complex, and DEDC-Cd(II) as the least stable. Structural data showed that the N1–C2 π-bond was localized in the complexes rather than the ligand, and a delocalized π-bond over S7–C2–S8 was also present. The stability of DEDC-Zn(II), DEDC-Cd(II), and DEDC-Hg(II) complexes increased in the presence of the non-specific effects of the solvent (PCM model), and their relative stability did not change. There was π-electron flow or resonance along N1–C2–S7 and along S7–C2–S8 in the ligand. The π-electron flow or resonance along N1–C2–S7 was abolished when the metal interacted with sulfur atoms. Energy belonging to van der Waals interactions and non-covalent delocalization effects between the metal and sulfur atoms of the ligand was calculated for each complex. The results of nucleus-independent chemical shift (NICS) indicated a decreasing trend as Zn(II)?<?Cd(II)?<?Hg(II) for the aromaticity of the quasi-rings. Finally, by ignoring van der Waals interactions and non-covalent delocalization effects between the metal and sulfur atoms of the ligand, the relative stability of the complexes was changed as follows:
At high temperature, silicon oxycarbide (SiCO) exhibits excellent mechanical properties and thermal stability. The incorporation of boron in SiCO results in improved performance in creep temperatures. In this work, large-scale molecular dynamics calculations were applied to obtain amorphous SiCO structures containing boron. Phase separation of C–C, B–C and Si–O was achieved for three compositions, and silicon-centered mixed-bond tetrahedrons were reproduced successfully. As the boron content increases, the boron atoms tend to form B–C and B–Si bonds in the voids, which stretches the free carbon network in some instances, causing a increase in C–C distance. Young’s modulus remains stable at high temperature for the high-carbon case, which indicates that the free carbon network plays a critical role in the structural and thermal stability of SiBCO.
Graphical Abstract Three major typical structures in the cooling down process for silicon boron oxycarbide (Si5BC2O8). Bonds: red Si–O, blue Si–C, black C–C, green B–C, purple Si–B
A post-calculation correction is established for PM7 band gaps of transition-metal oxides. The correction is based on the charge on the metal cation of interest, as obtained from MOPAC PM7 calculations. Application of the correction reduces the average error in the PM7 band gap from ~3 eV to ~1 eV. The residual error after correction is shown to be uncorrelated to the Hartree–Fock method upon which PM7 is based.
Graphical Abstract Comparison between calculated band gaps and experimental band gaps for binary oxides. The orange crosses are for corrected PM7 band gaps. Blue squares are uncorrected values. The orange crosses fall closer to the diagonal dashed line, showing an overall improvement of the accuracy of calculated values
Density functional theory (DFT) was utilized to elucidate the reaction mechanisms of and the key factors that influence the Ni(0)-catalyzed cross-dimerization and -trimerization of trimethylsilylacetylene (R1) and diphenylacetylene (R2). Calculated results revealed that the electron-donating ability of the ligand plays a crucial role in determining the regionselectivity of this tandem reaction. The use of strongly electron-donating ligands favors the formation of cross-dimer intermediates, whereas cross-trimer products can easily be synthesized using weakly electron-donating ligands. A simple method of estimating the electron-donating abilities of different ligands based on the Mulliken charge distribution of the ligand–ligand pair was employed. The present theoretical results allow us to elucidate the reaction mechanisms for and to identify the factors that exert the greatest influence on the ligand-controlled cross-dimerization and -trimerization of trimethylsilylacetylene and diphenylacetylene. Guidelines for the design of novel ligand systems with Ni(0) catalysts are also proposed.
Vitamin C is one of the most abundant exogenous antioxidants in the cell, and it is of the utmost importance to elucidate its mechanism of action against radicals. In this study, the reactivity of vitamin C toward OH and \( {HO}_2/{O}_2^{-} \) radicals in aqueous medium was analyzed by ab initio molecular dynamics using CPMD code. The simulations led to results similar to those of static studies or experiments for the pair of \( {HO}_2/{O}_2^{-} \) radicals but bring new insights for the reactivity with hydroxyl radical: the reaction takes place before the formation of an adduct and consists of two steps: first an electron is transferred to hydroxyl radical and then the ascorbyl radical loses a proton.
Efficient design of ionic compounds requires a systematic understanding of cation–anion interactions. Weakening of electrostatic attraction is essential to increase the liquid range of the ionic compound and decrease its melting point. Here, we report simulations of the closest-approach cation–anion distances in a variety of ion pairs containing the tetrakis(pentafluorophenyl)borate (TFPB–) anion. Small alkali cations (Li+, Na+) penetrate the TFPB– core, whereas K+ and larger organic cations do not. In the latter case, the shortest possible distance from the cations to the boron atom of TFPB– ranges from 0.50 nm to 0.63 nm. TFPB– was shown to be substantially rigid, providing a steric hindrance to thermodynamically efficient cation–anion coordination. Our results prove that TFPB– is more efficient for electrostatic charge confinement than the tetraoctylammonium cation, whereas the perfluorophenyl group is more efficient than linear alkyl chains. These simulations will motivate development of TFPB–-based ionic liquids with low phase transition points.
The structure of a hairpin loop—in particular its large accessible surface area and its exposed hydrogen-bonding edges—facilitate an inherent possibility for interactions. Just like higher-order RNA macromolecules, pre-microRNAs possess a hairpin loop, and it plays a crucial role in miRNA biogenesis. Upon inspecting the crystal structures of RNAs with various functions, we noticed that, along with a fairly long double helix, the RNAs contained sequentially different hairpin loops comprising four residues. We therefore applied molecular dynamics simulation to analyze six of these previously unexplored tetraloops, along with GNRA (where N is any nucleotide and R is a purine nucleotide) tetraloops, to understand their structural and functional characteristics. A number of analyses quantifying loop stability by examining base–base stacking, base–sugar and base–phosphate hydrogen bonding, and backbone variability were performed. Importantly, we determined the different interbase stacking preferences of the single-stranded unpaired bases of the hairpin loops, which had not previously been quantified in any form. Furthermore, our study indicates that canonical GNRA structural properties are exhibited by some structures containing non-GNRA loop sequences.
Angiotensin-converting enzyme (ACE) converts angiotensin I to angiotensin II and degrades bradykinin and other vasoactive peptides. ACE inhibitors are used to treat diseases such as hypertension and heart failure. It is thus highly desirable to understand the catalytic mechanism of ACE, as this should facilitate the design of more powerful and selective ACE inhibitors. ACE exhibits two different active domains, the C-domain and the N-domain. In this work, we systematically investigated the inhibitor- and substrate-binding patterns in the N-domain of human ACE using a combined quantum mechanical and molecular mechanical approach. The hydrolysis of hippuryl–histidyl–leucine (HHL) as catalyzed by the N-domain of human somatic ACE was explored, and the effects of chloride ion on the overall reaction were also investigated. Two models, one with and one without a chloride ion at the first binding position, were then designed to examine the chloride dependence of inhibitor–substrate binding and the catalytic mechanism. Our calculations indicate that the hydrolysis reaction follows a stepwise general base/general acid catalysis path. The estimated mean free energy barrier height in the two models is about 15.6 kcal/mol, which agrees very well with the experimentally estimated value of 15.8 kcal/mol. Our simulations thus suggest that the N-domain is in a mixed form during ACE-catalyzed hydrolysis, with the single-chloride-ion and the double-chloride-ion forms existing simultaneously.
The characterization of the seleno-sulfide-bromo systems and the isomerization process on the [H, S, Se, Br] potential energy surface were investigated using state-of-the-art theoretical methods. The CCSD(T) and the MP2 levels of theory were employed along with the series of correlation consistent basis sets extrapolated to the complete basis set (CBS) limit in the optimization of the geometrical parameters and computation of electronic energies. The relative stability, in kcal mol?1, at the CCSD(T)/CBS follows the trend: HSSeBr (0) > HSeSBr (9.51) > SSeHBr (24.02) > SeSHBr (25.42). This order was observed in the previous study of the [H, S, Se, Cl] species. The structural parameters and vibrational frequencies of the [H, S, Se, Br] species are reported. This research work should be helpful to experimentalists in order to gain insights into these novel heteroatom molecules.
Graphical abstract Relative energy profile (in kcal mol?1) using the CCSD(T)/CBS and MP2/CBS (in parentheses) method of the stationary states on the [H, S, Se, Br] PES.
Ecdysone receptor (EcR) is a significant target in the identification of new environmentally friendly pesticides. There are two types of ecdysone agonists: steroidal ecdysone agonists and dibenzoylhydrazines (DBHs). In this study, various modeling methods (homology modeling, molecular docking, MD simulation, binding free energy calculation, and per-residue binding free energy decomposition) were utilized to study the different binding mechanisms of two types of ecdysone agonists. Our theoretical results indicated that the relative binding potencies of DBHs can be ranked sufficiently accurately using the MOE docking method. However, MM/PBSA calculations more accurately predicted the binding affinities between steroidal ecdysone agonists and EcR-LBD. To identify the key residues involved in ecdysone agonist binding, the binding free energy (ΔGBind) was decomposed into the energy contributions of individual residues. The results revealed that nine residues—Ile339, Thr343, Met380, Met381, Tyr403, Tyr408, Asp419, Gln503, and Asn504—determined the binding affinities of the DBHs. Glu309, Met342, Arg383, Arg387, and Leu396 were important influences on the binding affinities of the steroidal ecdysone agonists.
Designed multi-target ligand (DML) is an emerging strategy for the development of new drugs and involves the engagement of multiple targets with the same moiety. In the context of NSAIDs it has been suggested that targeting the thromboxane prostanoid (TP) receptor along with cyclooxygenase-2 (COX-2) may help to overcome cardiovascular (CVS) complications associated with COXIBs. In the present work, azaisoflavones were studied for their COX-2 and TP receptor binding activities using structure based drug design (SBDD) techniques. Flavonoids were selected as a starting point based on their known COX-2 inhibitory and TP receptor antagonist activity. Iterative design and docking studies resulted in the evolution of a new class scaffold replacing the benzopyran-4-one ring of flavonoids with quinolin-4-one. The docking and binding parameters of these new compounds are found to be promising in comparison to those of selective COX-2 inhibitors, such as SC-558 and celecoxib. Owing to the lack of structural information, a model for the TP receptor was generated using a threading base alignment method with loop optimization performed using an ab initio method. The model generated was validated against known antagonists for TP receptor using docking/MMGBSA. Finally, the molecules that were designed for selective COX-2 inhibition were docked into the active site of the TP receptor. Iterative structural modifications and docking on these molecules generated a series which displays optimum docking scores and binding interaction for both targets. Molecular dynamics studies on a known TP receptor antagonist and a designed molecule show that both molecules remain in contact with protein throughout the simulation and interact in similar binding modes.
Density functional theory and its time-dependent extension (DFT, TDDFT) were employed to establish the feasibility of using a series of 4,4-difluoro-4-bora-3a,4a-diaza-s-indacenes (BODIPYs) in photodynamic therapy. Their absorption electronic spectra, singlet–triplet energy gaps, and spin–orbit matrix elements were computed and are discussed here. The effects of bromine substitution on the photophysical properties of BODIPY were elucidated. The investigated compounds were found to possess different excited triplet states that lie below the energy of the bright excited singlet state (S1 or S2), depending on the positions occupied by the bromine atoms. The computed spin–orbit matrix elements for the radiationless intersystem crossing Sn?→ ?Tm and the relative singlet–triplet energy gaps allowed the prediction of plausible nonradiative decay pathways for the production of singlet excited molecular oxygen, the key cytotoxic agent in photodynamic therapy.
Graphical Abstract The photophysical properties affected by the presence of bromine atoms in different positions of a BODIPY core have been here elucidated. In particular it has been found that SOC values strongly depend on the position of heavy atoms into the BODIPY core, suggesting positions 1 and 7 as the best ones to enhance the ISC kinetics
In this study, the doped defects in nitromethane crystals were investigated using first-principles calculations for the first time. We introduce dopant atoms in the interstitial sites of the nitromethane lattice, aiming to study the effects of element-doping on the structural properties, electronic properties, and sensitivity characteristics. The obtained results show that doped defects obviously affect the neighboring nitromethane molecules. The modification of electronic properties shows that the band gaps are significantly influenced by doped defects. Partial density of states and population analysis further reveal the mechanism for sensitivity control of nitromethane. It is shown that the new electronic states were introduced in the forbidden bands and the doped defects resulted in charge redistributions in the systems.
A new compound based on the D-π-A concept, where D = dimethylamino-phenyl and A = naphthoic acid, separated by an imine motif, was designed, synthesized and characterized. The spectral, energetics, and structural characteristics of the compound were studied thoroughly theoretically by density functional theory (DFT) in the gas and aqueous phases and experimentally (steady-state absorption) in aqueous media with various degrees of polarity and hydrogen bonding ability. This compound shows high sensitivity to the polarity, basicity and proton affinity of the environment. Based on DFT, TD-DFT and NBO analysis, the compound exists in the ground-state with both intermolecular and intramolecular hydrogen bond conformations in association with the –COOH, with latter isomer calculated to be more stable. Furthermore, structural changes via intermolecular solute–solvent interactions, dictate electronic modifications and spectral changes.
