Theoretical study of N<Subscript>4</Subscript>X (X = O,S, Se) systems

A series of N₄X (X = O, S, Se) compounds have been examined with ab initio and density functional theory (DFT) methods. To our knowledge, these compounds, except for the C_2v ring and the C_3v towerlike isomers of N₄O, are first reported here. The ring structures are the most energetically favored for N₄X (X = O and S) systems. For N₄Se, the cagelike structure is the most energetically favored. Several decomposition and isomerization pathways for the N₄X species have been investigated. The dissociation of C_2v ring N₄O and N₄S structures via ring breaking and the barrier height are only 1.1 and −0.2 kcal mol⁻¹ at the CCSD(T)/6-311+G*//MP2/6-311+G* level of theory. The dissociation of the cagelike N₄X species is at a cost of 12.1–16.2 kcal mol⁻¹. As for the towerlike and triangle bipyramidal isomers, their decomposition or isomerization barrier heights are all lower than 10.0 kcal mol⁻¹. Although the C_S cagelike N₄S isomer has a moderate isomerization barrier (18.3–29.1 kcal mol⁻¹), the low dissociation barrier (−1.0 kcal mol⁻¹) indicates that it will disappear when going to the higher CCSD(T) level. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

3D-QSAR studies on fluoropyrrolidine amides as dipeptidyl peptidase IV inhibitors by CoMFA and CoMSIA

Three-dimensional quantitative structure-activity relationship (3D-QSAR) analyses using CoMFA and CoMSIA methods were conducted on a series of fluoropyrrolidine amides as dipeptidyl peptidase IV (DP-IV) inhibitors. The selected ligands were docked into the binding site of the 3D model of DP-IV using the GOLD software, and the possible interaction models between DP-IV and the inhibitors were obtained. Based on the binding conformations of these fluoropyrrolidine amides and their alignment inside the binding pocket of DP-IV, predictive 3D-QSAR models were established by CoMFA and CoMSIA analyses, which had conventional r ² and cross-validated coefficient values () up to 0.982 and 0.555 for CoMFA and 0.953 and 0.613 for CoMSIA, respectively. The predictive ability of these models was validated by six compounds that were in the testing set. Structure-based investigations and the final 3D-QSAR results provide the guide for designing new potent inhibitors.     

The platinum-olefin binding energy in series of (PH<Subscript>3</Subscript>)<Subscript>2</Subscript>Pt(olefin) complexes - a theoretical study

Theoretical investigation of Pt(0)-olefin organometallic complexes containing tertiary phosphine ligands was focused on the strength of platinum-olefin electronic interaction. DFT theoretical study of electronic effects in a substantial number of ethylene derivatives was evaluated in terms of the Pt-olefin binding energy using MP2 correlation theory. Organometallics bearing coordinated olefins with general formula (R¹R²C = CR³R⁴)Pt(PH₃)₂ [R = various substituents] had been selected, including olefins containing both electron-donor substituents as well as electron-withdrawing groups. The stability of the corresponding complexes increases with a strengthening electron-withdrawal ability of the olefin substituents. Figure Representation of (CH₂ = CHR)Pt(PPh₃)₂ and the stability chart     

Density functional study of closed-shell attraction on X(ML)<Subscript>3</Subscript><Superscript>+</Superscript>(<Emphasis Type="Italic">X</Emphasis> = O,S, Se; <Emphasis Type="Italic">M</Emphasis> = Au,Ag, Cu) systems

With the help of quantum mechanical calculations, we have examined the series of central system X(ML)₃ ⁺(X = O, S, Se; M = Au, Ag, Cu). Using a scalar–relativistic density functional approach, we studied the geometry structures, Mulliken populations and charges of the systems. Structure parameters of the experimental systems are reproduced well with Xα method. The metallophilic interaction energy is analyzed and decomposed. For the systems with different central atoms and different metal atoms, the nature of the metallophilic attraction interaction is analyzed.     

Possible dynamic anchor points in a benzoxazinone derivative–human oxytocin receptor system — a molecular docking and dynamics calculation

In this study, we performed a molecular docking and dynamics simulation for a benzoxazinone–human oxytocin receptor system to determine the possible hydrophobic and electrostatic interaction points in the dynamic complex. After the homology modeling, the ligand was docked into the putative active using AutoDock 3.05. After the application of energetic and structural filters, the complexes obtained were further refined with a simulated annealing protocol (AMBER8) to remove steric clashes. Three complexes were selected for subjection to the molecular dynamics simulation (5 ns), and the results on the occurrence of average anchor points showed a stable complex between the benzoxazinone derivative and the receptor. The complex could be used as a good starting point for further analysis with site-directed mutagenesis, or further computational research. Figure The location of the ligands (complex B – blue; complex E – red; and complex F – green) in the transmembrane regions (TM1 – red; TM2 – blue; TM3 – yellow; TM4 – purple; TM5 – orange; TM6 – cyan; TM7 – pink) of the hOTR. For clarity, the EC and IC loops are not shown Electronic Supplementary Material Supplementary material is available at     

Sequence analysis and rule development of predicting protein stability change upon mutation using decision tree model

Understanding the mechanism of the protein stability change is one of the most challenging tasks. Recently, the prediction of protein stability change affected by single point mutations has become an interesting topic in molecular biology. However, it is desirable to further acquire knowledge from large databases to provide new insights into the nature of them. This paper presents an interpretable prediction tree method (named iPTREE-2) that can accurately predict changes of protein stability upon mutations from sequence based information and analyze sequence characteristics from the viewpoint of composition and order. Therefore, iPTREE-2 based on a regression tree algorithm exhibits the ability of finding important factors and developing rules for the purpose of data mining. On a dataset of 1859 different single point mutations from thermodynamic database, ProTherm, iPTREE-2 yields a correlation coefficient of 0.70 between predicted and experimental values. In the task of data mining, detailed analysis of sequences reveals the possibility of the compositional specificity of residues in different ranges of stability change and implies the existence of certain patterns. As building rules, we found that the mutation residues in wild type and in mutant protein play an important role. The present study demonstrates that iPTREE-2 can serve the purpose of predicting protein stability change, especially when one requires more understandable knowledge.     

Long-time molecular dynamics simulations of botulinum biotoxin type-A at different pH values and temperatures

Botulinum neurotoxins type A (BoNT/A) are highly potent toxins, but are also useful in the treatment of illnesses. We studied the properties of BoNT/A at various temperatures and pH values in order to understand its toxicity and structure variations. The pH values of the environment of BoNT/A are obtained by changing the protonation states of certain titratable residue groups. Our results show that certain parts of the protein are active at acidic pH environments or at high temperatures. The protein is more stable in neutral environments at normal human body temperature, whereas, at high temperature, the protein is more stable in acidic environments. Also, the three domains of the protein tend to have relative motion rather than within individual domains.     

Controls for Phylogeny and Robust Analysis in Pareto Task Inference

Understanding the tradeoffs faced by organisms is a major goal of evolutionary biology. One of the main approaches for identifying these tradeoffs is Pareto task inference (ParTI). Two recent papers claim that results obtained in ParTI studies are spurious due to phylogenetic dependence (Mikami T, Iwasaki W. 2021. The flipping t-ratio test: phylogenetically informed assessment of the Pareto theory for phenotypic evolution. Methods Ecol Evol. 12(4):696–706) or hypothetical p-hacking and population-structure concerns (Sun M, Zhang J. 2021. Rampant false detection of adaptive phenotypic optimization by ParTI-based Pareto front inference. Mol Biol Evol. 38(4):1653–1664). Here, we show that these claims are baseless. We present a new method to control for phylogenetic dependence, called SibSwap, and show that published ParTI inference is robust to phylogenetic dependence. We show how researchers avoided p-hacking by testing for the robustness of preprocessing choices. We also provide new methods to control for population structure and detail the experimental tests of ParTI in systems ranging from ammonites to cancer gene expression. The methods presented here may help to improve future ParTI studies.     

Expression of <Emphasis Type="Italic">Ruminococcus albus</Emphasis> Xylanase Gene (<Emphasis Type="Italic">xyn</Emphasis>A) in <Emphasis Type="Italic">Streptococcus bovis</Emphasis> 12-U-1

The objective of this study was to ligate the xylanase gene A (xynA) isolated from Ruminococcus albus 7 into the promoter and signal-peptide region of the lichenase [β-(1,3-1,4)-glucanase] gene of Streptococcus bovis JB1. This fusion gene was inserted into the pSBE11 vector, and the resulting recombinant, plasmid pXA, was used to transform S. bovis 12-U-1 cells. The transformant, S. bovis 12UXA, secreted the xylanase, which was stable against freeze-thaw treatment and long-time incubation at 37°C. The introduction of pXA and production of xylanase did not affect cell growth, and the xylanase produced degraded xylan from oat-spelt and birchwood. Received: 24 June 2002 / Accepted: 7 October 2002