Nontargeted parallel cascade selection molecular dynamics (nt-PaCS-MD) is a rare event sampling method of proteins, which does not rely on knowledge of the target structure. nt-PaCS-MD is an extension of targeted PaCS-MD (t-PaCS-MD). In nt-PaCS-MD, it makes use of cyclic resampling from some relevant initial structures to expand the searched conformational subspace. Reliable identification of these initial structures is the key to using nt-PaCS-MD. In the present study, we introduce the moving root-mean-square deviation (mRMSD) as a metric for identification of these statistical conformation outliers. mRMSD can be calculated for any ith geometry in the trajectory generated by short MD runs. The reference to which the mRMSD relates is the close surrounding of the ith conformation, often the (i-1) st one. Based on mRMSD, we show that it increases its effectiveness compared to the conventional MD.
To investigate the effects of forsterite structure on processes of H2O decomposition and proton diffusion, we performed molecular dynamics calculations on forsterites in glassy and crystalline states. The result shows that the decomposition rate of H2云南1快乐十分走势一定牛O on forsterite in the glassy state is higher than that in crystalline state in a temperature range of 200 to 400 K. Furthermore, the decomposed proton permeates into the internal part of glassy forsterite through hopping, whereas no permeation was observed for the crystalline state. These results suggest that the glassy structure of forsterite is an important factor for chemical evolution processes in interstellar spaces.
we demonstrate that accurate prediction of the crystal structures of anthracene and its derivatives is possible with the van der waals density functional. based on the calculated crystal structures, we investigate their electronic structures and discuss the correlation between crystal structure/molecular configuration and the electronic structure, in particular, the dispersion of the bands composed of the highest occupied molecular orbitals, which are important to discuss the charge carrier transport property.
An RNA aptamer that binds with high affinity and specificity to human immunoglobulin G (IgG) is a 24-mer single-stranded oligonucleotide containing 2′-fluoro pyrimidine nucleotides. Using X-ray crystallographic analysis, a Ca2+ ion was reported as being located near the G7 phosphate of the RNA aptamer. Surface plasmon resonance analysis showed that the RNA aptamer could not bind to IgG in a buffer without Ca2+ ions. To elucidate the role of Ca2+ ions in the binding of the RNA aptamer to IgG, we performed molecular dynamics simulation for the RNA aptamer/IgG complex with and without Ca2+ ions in a solvent. In the presence of Ca2+ ions, the RMSD of the RNA aptamer backbone remained below approximately 3.0 Å from the crystal structure during a 1,500-ns simulation. The distance between the centroids of the RNA aptamer and IgG was maintained between the centroids in the crystal structure. However, in the absence of Ca2+ ions, the RMSD increased and the structure of the RNA aptamer changed from the initial structure. These results indicate that Ca2+ ions play a role in maintaining the conformation of the RNA aptamer to the binding form.
Ferulic acid is known to have strong antioxidant properties. In the present study, we investigate the electronic structures of ferulic acid and its radical species extracting the hydrogen atom from its phenolic hydroxyl group. The relation of the results by several machine learning models using R/caret package, such as partial least squares, random forest, radial basis function kernel regularized least squares, and baysian regularized neural network, with the radical scavenging activity with the DPPH reagent, IC50, measured by Sakamoto et al. is discussed. We found all four methods gave reasonable correlation coefficients which means the possible prediction of the IC50 values with the results of the molecular orbital calculations only.
in this study, we attempted to quantify the shape of various metastable morphologies of microphase separation in block copolymers that can be observed in a computer simulation under the same condition as that for double gyroid (stable state). for the quantification, we used multifractal analysis, which is a tool to evaluate the shape having heterogeneous density inside it. our analysis showed that the indicators relating to the multifractal analysis of metastable structures are distributed far from that of double gyroid.
Ionic liquids (ILs) are non-volatile and highly stable solvents, which consist of organic/inorganic cations and anions. ILs are attractive since they can be utilized as gas-separating materials. Using a conductor-like screening model for realistic solvation (COSMO-RS) method, in this study, Henry's law constants of 15 gases (HCN, SO2, H2O, NH3, H2S, NO2, NO, CO, N2O, CO2, C2H6, C2H4, CH4, N2, H2) in 68,775 ILs were calculated to construct an ab initio database to search for ILs which are suitable for gas separation. Based on the database, it was numerically confirmed that not only the ions' geometric structures but also electronic structures are important to predict ILs' gas absorption abilities with high accuracy.
umbrella sampling is a method to calculate free energy by using molecular dynamics simulation. in the previous study, free energy reaction route mapping method (ferrmap) [5,6] is proposed, which is the method to calculate free energy reaction networks by using the umbrella integration method [1,2] and scaled hypersphere searching method [3,4]. in this study, we calculated ferns of alanine octapeptide (ala8) in water by using the ferrmap method on the dihedral of the ramachandran plot, which is 12 dimensional. we found 613 equation structures and 835 transition structures on the fern of ala8. this fern is too complicated to explain the folding of ala8, we propose an effective flow analysis method. by using this method, we found the representative paths connecting the beta-strand and alpha-helix structure.
A band gap prediction model of π-conjugated polymers was constructed using aromatic/quinoid, donor/acceptor, and torsion properties to predict quantitatively the band gap of π-conjugated polymers from properties of monomers. Quinoid stabilization energy (QSE), energy difference between HOMO of donor and LUMO of acceptor, torsion angle in homo-dimer of monomers were used as descriptors of aromatic/quinoid, donor/acceptor, and torsion properties. The neural network, which was constructed by 2 hidden layers with 5 neurons per layer, quantitatively predicts (RMSD = 0.207 eV, R2 = 0.885) the band gap of the π-conjugated polymers from descriptors of monomers.
According to the current superposition law in nanoscale, the conductance gets more than doubled when the number of paths is doubled. The superposition law can be applied to the molecular graph, which is a discrete mathematical expression for a structure consisting of nodes and edges. We previously revealed that the aromatic parallel circuit i.e. benzene molecule breaks the superposition law. In this paper, we show that [4n云南1快乐十分走势一定牛+2] annulenes also break the current superposition law from the viewpoint of orbital interaction between two fragments in the parallel circuit.
Reaction mechanisms for the isomerization of nitrile N-oxide to isocyanate catalyzed by hydroxide ion were investigated by density functional theory calculations. In this mechanism, nitrile N-oxide first dimerizes to a six-membered ring intermediate, and then undergoes two rearrangement reactions to isomerize into two isocyanate molecules. The process with the highest activation energy as an elementary reaction was the second rearrangement reaction, with an activation free energy of 31.1 kcal/mol, which is reasonable for a reaction proceeding slowly at room temperature. This mechanism can also explain the experimental results that the oxygen atom exchanges between nitrile N云南1快乐十分走势一定牛-oxide molecules during the isomerization reaction.
云南1快乐十分走势一定牛among ester synthesis methods, there is a method using a pyridine derivative as a catalyst, and the esterification reaction rate may be increased through a highly active acylpyridinium cation intermediate. it has been reported that the stability of acylpyridinium cations correlates with the reaction rate of esterification. in this study, we try to clarify the reason for the correlation between the reaction rate and the stability, then we investigate the relationship between the stability and structures of acylpyridinium cations. as a result, we found that the stability largely depends on the electronic donor property of the substituents, and acylpyridinium cation may act as an effective catalyst for esterification when the substituents have strong donor property at positions 3, 4, 5 of the pyridine ring.
云南1快乐十分走势一定牛all-atom molecular dynamics were performed on 4-heptyl-4’-cyanobiphenyl (7cb) to study the mechanism of heat conduction in the nematic liquid crystal. the non-equilibrium molecular dynamics simulation using our amended force field for 7cb has reproduced the anisotropy of the thermal conductivity in the nematic phase of 7cb. furthermore, the calculated phonon density of state shows that the higher thermal conductivity parallel to the director than that perpendicular to the director is due to the heat conduction derived from the stretching vibration of the well-oriented covalent bonds.