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.
In order to study the effects of N-oxide on structure and performance, six categories of energetic compounds were systemically investigated. The results indicated that the C–C bonds in the rings were shortened, and the C–N bonds close to the N?→?O bond were elongated when N atoms was oxidized to form N?→?O bonds. N?→?O bonds can increase the densities of most categories of compounds, and the increment will increase with the number of N?→?O bonds. As to their detonation performances, almost all categories of compounds had an increased trend, except for some NO2-, NHNO2- and ONO2-substituted compounds. The contribution of 1,2,3,4-tetrazine and 1,2,4,5-tetrazine to performances was better than that of pyrazine and [1,2,5] oxadiazolo [3,4-b] pyrazine on the whole, and the groups, especially energetic groups, made a huge contribution to performance. When R was a NH2 or ONO2 group, all compounds had lower impact sensitivities, and thus represent candidates for novel energetic compounds. However, other than the sixth category of compounds, all compounds had higher impact sensitivities when R was a NO2 or NHNO2 group, and have little significance in application.
Graphical abstract To study the effects of N-oxide on the structure and performance of energetic compounds, and to propose theoretical guidance for the design of novel compounds, the six categories (94 species) listed in the figure were investigated systemically by density functional theory methods and some empirical formulae
A molecular dynamics method was employed to study the binding energies associated with the cocrystallization (at selected crystal planes) of either 1,3,5-triamino-2,4,6-trinitro-benzene (TATB), 1,1-diamino-2,2-dinitroethylene, 3-nitro-1,2,4-triazol-5-one (TATB, FOX-7, and NTO, respectively, all of which are explosives), or N,N-dimethylformamide (DMF, a nonenergetic solvent) in various molar ratios with 1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane in its α and β conformations (α-HMX and β-HMX, respectively). The results showed that the cocrystals with low molar ratios (2:1, 1:1, 1:2, and 1:3) were the most stable. The binding energies of HMX/NTO and HMX/DMF were larger than those of HMX/TATB and HMX/FOX-7. According to the calculated stabilities, HMX prefers to adopt its α form in HMX/TATB and its β form in HMX/NTO, whereas the two forms coexist in HMX/FOX-7. For HMX/TATB, HMX/NTO, and α-HMX/FOX-7, increasing the proportion of the cocrystal component with the highest detonation heat (HMX in the first two cases, FOX-7 in the latter) increases the detonation heat, velocity, and pressure of the cocrystal. However, increasing the proportion of the component with the highest detonation heat in β-HMX/FOX-7 and γ-CL-20/FOX-7 increases the detonation heat of the cocrystal but decreases its detonation velocity. An investigation of the surface electrostatic potential revealed how the sensitivity changes upon cocrystal formation.
Beryllium telluride (BeTe) with cubic zinc-blende (ZB) structure was studied using ab initio constant pressure method under high pressure. The ab initio molecular dynamics (MD) approach for constant pressure was studied and it was found that the first order phase transition occurs from the ZB structure to the nickel arsenide (NiAs) structure. It has been shown that the MD simulation predicts the transition pressure PT more than the value obtained by the static enthalpy and experimental data. The structural pathway reveals MD simulation such as cubic → tetragonal → orthorhombic → monoclinic → orthorhombic → hexagonal, leading the ZB to NiAs phase. The phase transformation is accompanied by a 10% volume drop and at 80 GPa is likely to be around 35 GPa in the experiment. In the present study, our obtained values can be compared with the experimental and theoretical results.
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.
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
There has been considerable debate about the contribution of salt bridges to the stabilization of protein folds, in spite of their participation in crucial protein functions. Salt bridges appear to contribute to the activity–stability trade-off within proteins by bringing high-entropy charged amino acids into close contacts during the course of their functions. The current study analyzes the modes of association of salt bridges (in terms of networks) within globular proteins and at protein–protein interfaces. While the most common and trivial type of salt bridge is the isolated salt bridge, bifurcated salt bridge appears to be a distinct salt-bridge motif having a special topology and geometry. Bifurcated salt bridges are found ubiquitously in proteins and interprotein complexes. Interesting and attractive examples presenting different modes of interaction are highlighted. Bifurcated salt bridges appear to function as molecular clips that are used to stitch together large surface contours at interacting protein interfaces. The present work also emphasizes the key role of salt-bridge-mediated interactions in the partial folding of proteins containing long stretches of disordered regions. Salt-bridge-mediated interactions seem to be pivotal to the promotion of “disorder-to-order” transitions in small disordered protein fragments and their stabilization upon binding. The results obtained in this work should help to guide efforts to elucidate the modus operandi of these partially disordered proteins, and to conceptualize how these proteins manage to maintain the required amount of disorder even in their bound forms. This work could also potentially facilitate explorations of geometrically specific designable salt bridges through the characterization of composite salt-bridge networks.
The solubility advantage (SA) of meloxicam cocrystalized with mono- and dicarboxylic acids was expressed in terms of equilibrium constants involving active pharmaceutical ingredient and coformer in aqueous solutions. It is argued that SA can be quantified by concentration of pairs formed in water. The pH and concentration of dissolved components is included explicitly in the model. The alternative behavior of mono- and dicarboxylic acids was emphasized and addressed to different structural motifs. The structural and energetic properties of meloxicam and its complexes with carboxylic acids were characterized, including tautmerism and dissociation in aqueous media. In particular, performed in silico modeling confirmed experimental observation that meloxicam dissolved in water or modest acidic solutions is expected to be a mixture of anionic form in equilibrium with at least five neutral isomers. Tautomer-related diversity of pairs formation and the possibility of salt formation is also discussed.
Hyaluronan lyase from Streptococcus pneumoniae can degrade hyaluronic acid, which is one of the major components in the extracellular matrix. Hyaluronan can regulate water balance, osmotic pressure, and act as an ion exchange resin. Followed by our recent work on the catalytic reaction mechanism and substrate binding mode, we in this work further investigate the functional role of active site arginine residue, R462, in the degradation of hyaluronan. The site directed mutagenesis simulation of R462A and R462Q were modeled using a combined quantum mechanical and molecular mechanical method. The overall substrate binding features upon mutations do not have significant changes. The energetic profiles for the reaction processes are essentially the same as that in wild type enzyme, but significant activation barrier height changes can be observed. Both mutants were shown to accelerate the overall enzymatic activity, e.g., R462A can reduce the barrier height by about 2.8 kcal mol–1, while R462Q reduces the activation energy by about 2.9 kcal mol–1. Consistent with the active site model calculated using density functional theory, our results can support that the positive charge on R462 guanidino side chain group plays a negative role in the catalysis. Finally, the functional role of R462 was proposed to facilitate the formation of initial enzyme-substrate complex, but not in the subsequent catalytic degradation reaction.
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.
Recently, a series of xanthone analogues has been identified as α-glucosidase inhibitors. To provide deeper insight into the three-dimensional (3D) structural requirements for the activities of these molecules, CoMFA and CoMSIA approaches were employed on 54 xanthones to construct 3D-QSAR models. Their bioactive conformations were first investigated by docking studies and optimized by subsequent molecular dynamics (MD) simulations using the homology modeled structure of the target protein. Based on the docking/MD-determined conformers, 3D-QSAR studies generated several significant models in terms of 47 molecules as the training set. The best model (CoMSIA-SHA) yielded q2 of 0.713, r2 of 0.967 and F of 140.250. The robustness of the model was further externally confirmed by a test set of the remaining molecules (q2 = 0.793, r2 = 0.902, and k = 0.905). Contour maps provided much information for future design and optimization of new compounds with high inhibitory activities towards α-glucosidase.
Spin-polarized first-principles total-energy calculations have been performed to investigate the possible chain reaction of acetylene molecules mediated by hydrogen abstraction on hydrogenated hexagonal boron nitride monolayers. Calculations have been done within the periodic density functional theory (DFT), employing the PBE exchange correlation potential, with van der Waals corrections (vdW-DF). Reactions at two different sites have been considered: hydrogen vacancies on top of boron and on top of nitrogen atoms. As previously calculated, at the intermediate state of the reaction, when the acetylene molecule is attached to the surface, the adsorption energy is of the order of ?0.82 eV and ?0.20 eV (measured with respect to the energy of the non interacting molecule-substrate system) for adsorption on top of boron and nitrogen atoms, respectively. After the hydrogen abstraction takes place, the system gains additional energy, resulting in adsorption energies of ?1.52 eV and ?1.30 eV, respectively. These results suggest that the chain reaction is energetically favorable. The calculated minimum energy path (MEP) for hydrogen abstraction shows very small energy barriers of the order of 5 meV and 22 meV for the reaction on top of boron and nitrogen atoms, respectively. Finally, the density of states (DOS) evolution study helps to understand the chain reaction mechanism.
The main purpose of this study was to address the performance of virtual screening methods based on ligands and the protein structure of acetylcholinesterase (AChE) in order to retrieve novel human AChE (hAChE) inhibitors. In addition, a protocol was developed to identify novel hit compounds and propose new promising AChE inhibitors from the ZINC database with 10 million commercially available compounds. In this sense, 3D similarity searches using rapid overlay of chemical structures and similarity analysis through comparison of electrostatic overlay of docked hits were used to retrieve AChE inhibitors from collected databases. Molecular dynamics simulation of 100 ns was carried out to study the best docked compounds from similarity searches. Some key residues were identified as crucial for the dual binding mode of inhibitor with the interaction site. All results indicated the relevant use of EON and docking strategy for identifying novel hit compounds as promising potential anticholinesterase candidates, and seven new structures were selected as potential hAChE inhibitors.
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.
Mechanisms for the activation of water, ammonia, and other small molecules by the PCcarbeneP nickel pincer complex were studied computationally with the aid of density functional theory. The calculation results indicate that the strongly donating, nucleophilic carbene center can engage in a variety of heterolytic splitting of E?H (E=H, C, N, O) bonds, some of which are reversible. The cleavage of E?H bonds across the Ni=C bond represents a new mode of bond activation by ligand cooperativity in nickel pincer complex. On the basis of the calculations, we also demonstrate that reversible H2 activation across the Ir=C bond via the PCcarbeneP iridium pincer complex was observed in the experiments, while other E?H (E=C, N, O) bonds were not activated. Our calculations are in good agreement with experimental observations and could provide new insights into ligand cooperativity in nickel pincer complexes.
Our investigation is devoted to the theoretical study of the low-lying electronic structure of the LaCl molecule by using ab initio quantum methods. We are concerned with several methods such as the complete active space-self consistent field (CAS-SCF) and the multi reference of configuration interaction (MRCI + Q) methods. These methods are applied for the purpose of drawing the potential energy curves (PECs) and calculating the molecular spectroscopic constants for a given number of electronic states of singlet and triplet multiplicity. We count 26 2S+?1 Λ(±) electronic states located below 24,000 cm??1 neglecting the spin-orbit effects and 47 Ω(±) components taken into consideration these effects. Our calculations are performed via the quantum ab initio package MOLPRO (Werner and Knowles 2000).
Graphical Abstract A new set of low-lying electronic states on the theoretical energetic level diagram for the LaCl molecule among the first four lanthanum monhalides.
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.
