Theoretical study of the intermolecular potential energy and second virial coefficient in the mixtures of CH4 and Kr gases: a comparison with experimental data

The intermolecular potential energy surface (IPS) in the mixtures of CH₄–Kr gases from ab initio calculations has been explored. The ab initio calculation was performed at the second-order Møller–Plesset perturbation theory (MP₂), with the 6-311+G(2df,2pd) basis set, for three relative orientations of two CH₄–Kr molecules as a function of CH₄–Kr separation distance. In this work, the IPS, U(r), of the CH₄–Kr complex has been investigated, where the vertex (V), edge (E) and face (F) of CH₄ approaches to Kr have been considered. Then, adjustable parameters of the Lennard-Jones and Buckingham potential energy function are fitted to the ab initio MP2/6-311+G(2df,2pd) interaction energies for three different orientations. Assuming a given set of parameters, we theoretically obtained second virial coefficients for the CH₄–Kr system, and compared with the experimental data at different temperatures. Trivial differences can be observed between the experimental and computational results.     

QM/MM model study on properties and structure of some antibiotics in gas phase: Comparison of energy and NMR chemical shift

The combination of Quantum Mechanics (QM) and Molecular Mechanics (MM) methods has become an alternative tool for many applications for which pure QM and MM are not suitable. The QM/MM method has been used for different types of problems, for example: structural biology, surface phenomena, and liquid phase. In this paper, we have used these methods for antibiotics and then we compared results. The calculations were done by the full ab initio method (HF/3-21G) and the (HF/STO-3G) and QM/MM (ONIOM) method with HF (3-21G)/AM1/UFF and HF (STO-3G)/AM1/UFF. We found the geometry obtained by the QM/MM method to be very accurate, and we can use this rapid method in place of time consuming ab initio methods for large molecules. Comparison of energy values in the QM/MM and QM methods is given. In the present work, we compare chemical shifts and conclude that the QM/MM method is a perturbed full QM method. The work has been done on penicillin, streptomycin, benzyl penicillin, neomycin, kanamycin, gentamicin, and amoxicillin.     

Simulation of DNA bases in water: Comparison of the Monte Carlo algorithm with molecular mechanics force fields

The interaction between the nucleic acid bases and solvent molecules has an important effect in various biochemical processes. We have calculated total energy and free energy of the solvation of DNA bases in water by Monte Carlo simulation. Adenine, guanine, cytosine, and thymine were first optimized in the gas phase and then placed in a cubic box of water. We have used the TIP3 model for water and OPLS for the nucleic acid bases. The canonical (T, V, N) ensemble at 25°C and Metropolis sampling technique have been used. Good agreement with other available computational data was obtained. Radial distribution functions of water around each site of adenine, guanine, cytosine, and thymine have been computed and the results have shown the ability of the sites for hydrogen bonding and other interactions. The computations have shown that guanine has the highest value of solvation free energy and N7 and N6 in adenine and guanine, N3 in cytosine, and N3 and O4 in thymine have the largest radial distribution function. Monte Carlo simulation has also been performed using the CHARMM program under the same conditions, and the results of two procedures are compared.     

Thermodynamic study of interaction of TSPP,CoTsPc, and FeTsPc with calf thymus DNA

The interaction of two water soluble phthalocyanines, cobalt(II) 4,4′,4″,4‴-tetrasulfo-phthalocyanine (CoTsPc) and iron(II) 4,4′,4″,4‴-tetrasulfo-phthalocyanine (FeTsPc), and one water soluble porphyrin, tetra sodium mesotetrakis(p-sulfophenyl)porphyrin (TSPP), with calf thymus DNA has been studied by UV-Vis spectroscopy at five different temperatures (20, 25, 30, 35, and 40°C). The optical absorption spectra of these materials were analyzed to obtain binding constants and stoichiometries using SQUAD software. The results show that the best fitting corresponds to a 1: 1 complex model between a base pair of DNA and these materials. All of the studied porphyrin and phthalocyanines showed strong electrolyte effect, and increasing NaCl concentration induced self-aggregation of these materials. Published in Russian in Biokhimiya, 2007, Vol. 72, No. 6, pp. 799–804.     

Simulation of DNA bases in water: Comparison of the Monte Carlo algorithm with molecular mechanics force fields

The interaction between the nucleic acid bases and solvent molecules has an important effect in various biochemical processes. We have calculated total energy and free energy of the solvation of DNA bases in water by Monte Carlo simulation. Adenine, guanine, cytosine, and thymine were first optimized in the gas phase and then placed in a cubic box of water. We have used the TIP3 model for water and OPLS for the nucleic acid bases. The canonical (T, V, N) ensemble at 25 degrees C and Metropolis sampling technique have been used. Good agreement with other available computational data was obtained. Radial distribution functions of water around each site of adenine, guanine, cytosine, and thymine have been computed and the results have shown the ability of the sites for hydrogen bonding and other interactions. The computations have shown that guanine has the highest value of solvation free energy and N7 and N6 in adenine and guanine, N3 in cytosine, and N3 and O4 in thymine have the largest radial distribution function. Monte Carlo simulation has also been performed using the CHARMM program under the same conditions, and the results of two procedures are compared.     

QM/MM model study on properties and structure of some antibiotics in gas phase: Comparison of energy and NMR chemical shift

The combination of Quantum Mechanics (QM) and Molecular Mechanics (MM) methods has become an alternative tool for many applications for which pure QM and MM are not suitable. The QM/MM method has been used for different types of problems, for example: structural biology, surface phenomena, and liquid phase. In this paper, we have used these methods for antibiotics and then we compare results. The calculations were done by the full ab initio method (HF/3-21G) and the (HF/STO-3G) and QM/MM (ONIOM) method with HF (3-21G)/AM1/UFF and HF (STO-3G)/AM1/UFF. We found the geometry that has obtained by the QM/MM method to be very accurate, and we can use this rapid method in place of time consuming ab initio methods for large molecules. Comparison of energy values in the QM/MM and QM methods is given. In the present work, we compare chemical shifts and conclude that the QM/MM method is a perturbed full QM method. The work has been done on penicillin, streptomycin, benzyl penicillin, neomycin, kanamycin, gentamicin, and amoxicillin.     

Solvent Dielectric Effect and Side Chain Mutation on the Structural Stability of Burkholderia cepacia Lipase Active Site: A Quantum Mechanical/Molecular Mechanics Study

Quantum mechanical and molecular dynamics methods were used to analyze the structure and stability of neutral and zwitterionic configurations of the extracted active site sequence from a Burkholderia cepacia lipase, histidyl-seryl-glutamin (His86-Ser87-Gln88) and its mutated form, histidyl-cysteyl-glutamin (His86-Cys87-Gln88) in vacuum and different solvents. The effects of solvent dielectric constant, explicit and implicit water molecules and side chain mutation on the structure and stability of this sequence in both neutral and zwitterionic forms are represented. The quantum mechanics computations represent that the relative stability of zwitterionic and neutral configurations depends on the solvent structure and its dielectric constant. Therefore, in vacuum and the considered non-polar solvents, the neutral form of the interested sequences is more stable than the zwitterionic form, while their zwitterionic form is more stable than the neutral form in the aqueous solution and the investigated polar solvents in most cases. However, on the potential energy surfaces calculated, there is a barrier to proton transfer from the positively charged ammonium group to the negatively charged carboxylat group or from the ammonium group to the adjacent carbonyl oxygen and or from side chain oxygen and sulfur to negatively charged carboxylat group. Molecular dynamics simulations (MD) were also performed by using periodic boundary conditions for the zwitterionic configuration of the hydrated molecules in a box of water molecules. The obtained results demonstrated that the presence of explicit water molecules provides the more compact structures of the studied molecules. These simulations also indicated that side chain mutation and replacement of sulfur with oxygen leads to reduction of molecular flexibility and packing.     

In silico study of effects of polymorphisms on biophysical chemical properties of oxidized N-terminal domain of X-ray cross-complementing group 1 protein

Base excision repair (BER) is the major pathway involved in removal of endogenous and mutagen-induced DNA damage. The X-ray cross-complementing group 1 protein (XRCC1), which participates in BER, is a scaffolding protein. The oxidized XRCC1 N-terminal domain (NTD) forms additional interactions with DNA polymerase β (Pol β). Any change in the residues of a protein (XRCC1, XRCC4, etc.) may alter its stability and function. Many coding regions of genes have single nucleotide polymorphisms (SNPs) that change the conformation of their products, and they are probably involved in some diseases. The R7L and R107H mutations are located in the XRCC1-NTD. In the present study, biophysical chemical properties of oxidized XRCC1-NTD (wild type or mutants) were investigated at different temperatures (290, 295, 298, 301, 304, 309, 310, 311, and 312 K) in water using in silico molecular mechanic computational methods. Comparison of the average calculated potential energies of oxidized XRCC1-NTD reveals that the R7L mutation increases stability, but the R107H and R7L&R107H mutations are destabilizing. Therefore, mutant types of this protein (R107H or R7L&R107H) may not function correctly. Furthermore, quantitative structure-activity relationship (QSAR) of oxidized XRCC1-NTD and docking assay showed that the R7L mutation is advantageous but the R107H and R7L&R107H mutations are disadvantageous for XRCC1-NTD, and in the latter cases it cannot interact with Pol β as well as the wild type does. Hence, DNA repair may be defective. Also, using the equation dE = ?E/(?T)_V·dT + ?E/(?V)_T·dV, it was determined that the best temperature for normal activity of oxidized XRCC1-NTD is exactly the natural body temperature (310 K).     

Metal-doped graphene layers composed with boron nitride–graphene as an insulator: a nano-capacitor

A new family of energetic azacubane N-oxides were designed by introducing N-oxides into azacubanes and investigated by using density functional theory. Introducing the N-oxides into the azacubanes could improve their detonation performance significantly due to the increase of the OB and ρ but would also increase the sensitivity to some extent. These effects would be further enhanced as the numbers of N-oxides increase. Among all the designed azacubane N-oxides, D6-4 (1,3,5,7-tetraazacubane-1,3,5,7-tetraoxides) has higher detonation performance than one famous high explosive HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocane) and lower sensitivity than one very insensitive explosive TNT (1-methyl-2,4,6-trinitrobenzene), suggesting that its overall performance is outstanding and may be considered as the potential candidate of insensitive high explosives. The internal small cage C-N skeleton of D6-4 is surrounded by the external big cage hydrogen bonds and this special double cage structure may be an important reason why it has low sensitivity.