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.
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.
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
Density functional theory (B3LYP, B3LYP-D2 and wB97XD functionals) was used in finite models of zigzag carbon nanotubes (CNT), (n,0)×k with n?=?6–9 and k?=?2–4, to systematically investigate the effects of size on their structural and electronic properties. We found that the ratio between the length (Lt) and the diameter (dt) of the pristine CNT has to be larger than 2, i.e., Lt/dt?>?2, in order to provide the observed experimental trends of C=C bond distances, as well as to maintain the atomic charges nearly constant and zero around the center of the tube. Therefore, the concepts of useful length and volume were developed and tested for the encapsulation process of HCN and C2H2 into CNTs. The energies involved in these processes, as well as the changes in molecular structure and electronic properties of the dopants and the CNTs are discussed and rationalized by the amount of charge transferred between dopant and CNT.
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.
In this work, through a docking analysis of compounds from the ZINC chemical library on human β-tubulin using high performance computer cluster, we report new polycyclic aromatic compounds that bind with high energy on the colchicine binding site of β-tubulin, suggesting three new key amino acids. However, molecular dynamic analysis showed low stability in the interaction between ligand and receptor. Results were confirmed experimentally in in vitro and in vivo models that suggest that molecular dynamics simulation is the best option to find new potential β-tubulin inhibitors.
In this work, we address the effects of molecular doping on the electronic properties of fluorinated and chlorinated silicon nanowires (SiNWs), in comparison with those corresponding to hydrogen-passivated SiNWs. Adsorption of n-type dopant molecules on hydrogenated and halogenated SiNWs and their chemisorption energies, formation energies, and electronic band gap are studied by using density functional theory calculations. The results show that there are considerable charge transfers and strong covalent interactions between the dopant molecules and the SiNWs. Moreover, the results show that the energy band gap of SiNWs changes due to chemical surface doping and it can be further tuned by surface passivation. We conclude that a molecular based ex-situ doping, where molecules are adsorbed on the surface of the SiNW, can be an alternative path to conventional doping.
A series of three star-shaped compounds containing both donor (carbazole) and acceptor (2,4,6-triphenyl-1,3,5-triazine) moieties linked through various linking bridges was studied theoretically at the linear response TD-DFT level of theory to describe their absorption and fluorescence spectra. The concept of a localized charge-transfer excited state has been applied successfully to explain the observed strong solvatochromic effect in the emission spectra of the studied molecules, which can be utilized for the fabrication of color tunable solution-processable OLEDs. The concept is in particularly applicable to donor–acceptor species with a C3 symmetry point group where the static dipole moment changes dramatically upon electronic excitation. An important peculiarity of the studied molecules is that they are characterized by non-zero values of the HOMO and LUMO orbitals in the same common part of molecular space that provides a large electric dipole transition moment for both light absorption and emission.
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.
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.
The anti-hypertensive drugs amlodipine, atenolol and lisinopril, in ordinary and PEGylated forms, with different combined-ratios, were studied by molecular dynamics simulations using GROMACS software. Twenty simulation systems were designed to evaluate the interactions of drug mixtures with a dimyristoylphosphatidylcholine (DMPC) lipid bilayer membrane, in the presence of water molecules. In the course of simulations, various properties of the systems were investigated, including drug location, diffusion and mass distribution in the membrane; drug orientation; the lipid chain disorder as a result of drug penetration into the DMPC membrane; the number of hydrogen bonds; and drug surface area. According to the results obtained, combined drugs penetrate deeper into the DMPC lipid bilayer membrane, and the lipid chains remain ordered. Also, the combined PEGylated drugs, at a combination ratio of 1:1:1, enhance drug penetration into the DMPC membrane, reduce drug agglomeration, orient the drug in a proper angle for easy penetration into the membrane, and decrease undesirable lipotoxicity due to distorted membrane self-assembly and thickness.
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.
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.
It is known that inhibiting α-amylase, an important enzyme in digestion of starch and glycogen, is a useful strategy for treating disorders in carbohydrate uptake. Two natural components distributed in many fruits and plants, oleanolic acid and ursolic acid, are endowed with important pharmacological activities and wide therapeutic possibilities. Until now, only a tiny fraction of their applications have been identified and exploited. Our in vitro inhibition studies demonstrated that oleanolic acid and ursolic acid non-competitively inhibit the activity and function of human salivary α-amylase. The molecular simulations revealed that oleanolic acid and ursolic acid interact with amino acid residues within the binding pocket of human salivary α-amylase, among which the side chain of Arg195 and Asp 197 was supposed to be important in imparting the inhibitory activity of triterpenoids. The present work will provide meaningful information for future development of functional drugs for the treatment of disorders in carbohydrate metabolism.
Graphical abstract This work is valuable for providing a deeper insight into the interaction mechanism of oleanolic acid and ursolic acid with α-amylase
The activation mechanism of dopamine receptors is unknown. The amino acids S5.42, S5.43, and S5.46 located in helix 5 appear to be crucial, but their specific roles in receptor activation have not been studied. We modeled the D1 dopamine receptor using the crystal structures of the D3 dopamine and β2 adrenergic receptors. Molecular dynamics simulations show that the interaction of dopamine with the D1 receptor leads to the formation of a hydrogen-bond network with its catechol group and helices 3, 5, and 6, including water molecules. The para hydroxyl group of dopamine binds directly to S5.42 and N6.55, the latter also interacting with S5.43. Unexpectedly, S5.46 does not interact directly with the catechol; instead, it interacts through a water molecule with S5.42 and directly with T3.37. The formation of this hydrogen-bond network, part of which was previously observed in docking studies with dopamine agonists, triggers the opening of the E6.30–R3.60 ionic lock associated with the activation of GPCRs. These changes do not occur in the unbonded (apo) receptor or when it is in a complex with the antagonist 3-methoxy-5,6,7,8,9,14-hexahydrodibenz[d,g]azecine. Our results provide valuable insight into the T3.37–S5.46–water–S5.43–ligand interaction, which may be crucial to the activation of the D1 dopamine receptor and should be considered during the design of novel agonists.
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.
