NMR and NQR parameters of the SiC-doped on the (4,4) <Emphasis Type="Italic">armchair</Emphasis> single-walled BPNT: a computational study

The structural properties, NMR and NQR parameters in the pristine and silicon carbide (SiC) doped boron phosphide nanotubes (BPNTs) were calculated using DFT methods (BLYP, B3LYP/6-31G*) in order to evaluate the influence of SiC-doped on the (4,4) armchair BPNTs. Nuclear magnetic resonance (NMR) parameters including isotropic (CS^I) and anisotropic (CS^A) chemical shielding parameters for the sites of various ¹³C, ²⁹Si, ¹¹B, and ³¹P atoms and quadrupole coupling constant (C _Q ), and asymmetry parameter (η _Q ) at the sites of various ¹¹B nuclei were calculated in pristine and SiC- doped (4,4) armchair boron phosphide nanotubes models. The calculations indicated that doping of ¹¹B and ³¹P atoms by C and Si atoms had a more significant influence on the calculated NMR and NQR parameters than did doping of the B and P atoms by Si and C atoms. In comparison with the pristine model, the SiC- doping in Si_PC_B model of the (4,4) armchair BPNTs reduces the energy gaps of the nanotubes and increases their electrical conductance. The NMR results showed that the B and P atoms which are directly bonded to the C atoms in the SiC-doped BPNTs have significant changes in the NMR parameters with respect to the B and P atoms which are directly bonded to the Si atoms in the SiC-doped BPNTs. The NQR results showed that in BPNTs, the B atoms at the edges of nanotubes play dominant roles in determining the electronic behaviors of BPNTs. Also, the NMR and NQR results detect that the Fig. 1b (Si_PC_B) model is a more reactive material than the pristine and the Fig. 1a (Si_BC_p) models of the (4,4) armchair BPNTs.     

Empirical solvation models can be used to differentiate native from near-native conformations of bovine pancreatic trypsin inhibitor.

Several hydration models for peptides and proteins based on solvent accessible surface area have been proposed previously. We have evaluated some of these models as well as four new ones in the context of near-native conformations of a protein. In addition, we propose an empirical site-site distance-dependent correction that can be used in conjunction with any of these models. The set of near-native structures consisted of 39 conformations of bovine pancreatic trypsin inhibitor (BPTI) each of which was a local minimum of an empirical energy function (ECEPP) in the absence of solvent. Root-mean-square (rms) deviations from the crystallographically determined structure were in the following ranges: 1.06-1.94 A for all heavy atoms, 0.77-1.36 A for all backbone heavy atoms, 0.68-1.33 A for all alpha-carbon atoms, and 1.41-2.72 A for all side-chain heavy atoms. We have found that there is considerable variation among the solvent models when evaluated in terms of concordance between the solvation free energy and the rms deviations from the crystallographically determined conformation. The solvation model for which the best concordance (0.939) with the rms deviations of the C alpha atoms was found was derived from NMR coupling constants of peptides in water combined with an exponential site-site distance dependence of the potential of mean force. Our results indicate that solvation free energy parameters derived from nonpeptide free energies of hydration may not be transferrable to peptides. Parameters derived from peptide and protein data may be more applicable to conformational analysis of proteins. A general approach to derive parameters for free energy of hydration from ensemble-averaged properties of peptides in solution is described.     

High-pressure phase transition and thermoelastic properties of europium chalcogenides

The pressure-induced crystal properties of Eu chalcogenides were investigated using two different models: a modified charge-transfer potential model consisting of Coulomb screening due to the delocalization of the f electron of the rare earth atom, and modified by covalency and zero-point energy effects along with attractive and repulsive interactions; and a charge-transfer model that excluded the covalency and zero-point energy effects in the previous model. Both models were used to visualize the effect of covalency on the mechanism of interaction of the constituent atoms. Eu chalcogenides transform from the Fm3m to the Pm3m phase under the influence of sufficient pressure (P(T) = 39.52, 21.01, 14.31, and 10.58 GPa), and their equations of state indicated decreases in volume during this phase transition of 6.38, 12.32, 12.76, and 11.15%, respectively, for EuO, EuS, EuSe, and EuTe. The results obtained from the models were in good agreement with corresponding experimental data. The elastic constants and Debye temperatures were also computed at normal and high pressures. Both of the models were found to be capable of successfully explaining these properties.     

Esterification of fatty acids using nylon-immobilized lipase in n-hexane: kinetic parameters and chain-length effects.

The esterification of long-chain fatty acids in n-hexane catalyzed by nylon-immobilized lipase from Candida rugosa has been investigated. Butyl oleate (22 carbon atoms), oleyl butyrate (22 carbon atoms) and oleyl oleate (36 carbon atoms) were produced at maximum reaction rates of approximately equal to 60 mmol h(-1) g(-1) immobilized enzyme when the substrates were present in equimolar proportions at an initial concentration of 0.6 mol l(-1). The observed kinetic behavior of all the esterification reactions is found to follow a ping-pong bi-bi mechanism with competitive inhibition by both substrates. The effect of the chain-length of the fatty acids and the alcohols could be correlated to some mechanistic models, in accordance with the calculated kinetic parameters.     

Charge sensitivity approach to mutual polarization of reactants: molecular mechanics perspective

Charge sensitivity analysis (CSA) in force-field atoms resolution was applied to describe the mutual polarization of reactants as well as charge-transfer (CT) effects. An inclusion complex of β-cyclodextrin with salicylic acid was used as a model system. Three CSA models were taken into account and verified on a Born–Oppenheimer molecular dynamics (BOMD) trajectory. The models differed in terms of the equilibrium conditions imposed on the system. It was demonstrated that mutual polarization is an important source of stabilization, in contrast to the results obtained from static charge calculations. The energy lowering induced by CT was small and comparable to the CT stabilization that occurs in hydrogen-bonded systems. All models correctly described the main topological features of the BOMD energy surface. CSA in force-field atoms resolution qualitatively reproduced the charge reorganization accompanying hydrogen-bond formation. It was shown that CSA parameters are very sensitive to the bond formation process, which suggests that they could be applied in reactive force fields as detectors of newly formed chemical bonds.

QSAR study on (13)C NMR chemical shifts on carbinol carbon atoms

QSAR calculations of (13)C NMR chemical shifts (ppm, TMS=0) on carbinol carbon atoms have been attempted using a large set of distance based topological indices: Wiener (W)-, Szeged (Sz)-, PI (Padmakar-Ivan) and Connectivity ((m)chi, (m)chi(v)) indices. The regression analysis has shown that excellent results are obtained in multiparametric regression. The predictive power of the proposed models are discussed using cross-validated parameters.     

Free energy of helical transmembrane peptide dimerization in OPLS-AA/Berger force field simulations: inaccuracy and implications for partner-specific Lennard-Jones parameters between peptides and lipids

Interactions between transmembrane (TM) peptides are important in biophysical chemistry, but there are few studies assessing atomistic simulation parameters concerning the energetics of interactions of TM helical peptides. Our potential of mean force analysis using OPLS-AA protein/Berger lipid force fields (FFs) shows that the dimerisation energy of (AALALAA)₃ helical peptides in the dioleoylphosphatidylcholine bilayer is ?4.4 kJ/mol, which was much smaller than the reported experimental value (?12.7 kJ/mol), thus calling for improvement of parameters of the combined FFs. As each of the FFs has been independently developed, we then tested the effects of downscaling the Lennard-Jones (LJ) energy terms between the OPLS-AA atoms and Berger lipid atoms, preserving the parameters within each FF. A 0.9-fold rescaling of the LJ energies was found to render the dimerisation energy close to the experimental value. Solvation of backbone atoms as well as side chain atoms in lipids is crucial for the TM helix interaction. In similar analyses, GROMOS 53A6 FF exhibited as weak dimerisation propensity (~?5.2 kJ/mol) as OPLS-AA/Berger, but CHARMM36 showed relatively accurate propensity (~?9.9 kJ/mol). Challenges and strategies in rendering the TM interaction energy realistic within the framework of current FFs are discussed.     

New information content in RNA base pairing deduced from quantitative analysis of high-resolution structures

Non-canonical base pairs play important roles in organizing the complex three-dimensional folding of RNA. Here, we outline methodology developed both to analyze the spatial patterns of interacting base pairs in known RNA structures and to reconstruct models from the collective experimental information. We focus attention on the structural context and deformability of the seven pairing patterns found in greatest abundance in the helical segments in a set of well-resolved crystal structures, including (i–ii) the canonical A·U and G·C Watson–Crick base pairs, (iii) the G·U wobble pair, (iv) the sheared G·A pair, (v) the A·U Hoogsteen pair, (vi) the U·U wobble pair, and (vii) the G·A Watson–Crick-like pair. The non-canonical pairs stand out from the canonical associations in terms of apparent deformability, spanning a broader range of conformational states as measured by the six rigid-body parameters used to describe the spatial arrangements of the interacting bases, the root-mean-square deviations of the base-pair atoms, and the fluctuations in hydrogen-bonding geometry. The deformabilties, the modes of base-pair deformation, and the preferred sites of occurrence depend on sequence. We also characterize the positioning and overlap of the base pairs with respect to the base pairs that stack immediately above and below them in double-helical fragments. We incorporate the observed positions of the bases, base pairs, and intervening phosphorus atoms in models to predict the effects of the non-canonical interactions on overall helical structure.     

A versatile modelling approach to determine the hydrophobicity of peptides at the atomic level

This study describes a versatile computational method to determine the hydrophobicity of small peptides at the atomic level. Free energies of transfer for individual atoms in peptide structures were derived, utilising two specifically defined parameters: (i) the water-excluding distance to define the dynamic interface between a peptide solute and its surrounding solvent and (ii) the corresponding hydrophobicity index as a relative measure for water occlusion/repulsion. The method was tested on a range of small peptide models (Ac-X-NH₂, G-X-G, Ac-WL-X-LL and Ac-GG-X-GG-NH₂) and several derivatives of these structures, whereby X was any of the 20 most common amino acids that naturally occur in polypeptides or proteins. The advantage of this new method lies in its versatility, ease to implement and capability to provide information on the hydrophobicity characteristics at the atomic level. The approach also encapsulates the impact of factors that influence these properties, but which have hitherto been difficult to accurately quantify, e.g. steric hindrance or proximity effects due to nearby polarised atoms. The method is not conditional on the knowledge of hydrophobicity parameters from the literature and does not require a sophisticated computer software/hardware to enable the atomic solvent-accessible surface areas or other hydrophobicity parameters to be de novo obtained.     

Antitumor promoting activities of 3-O-acyl-(-)epigallocatechins

As an exploratory investigation of antitumor promoting compounds, 3-O-acyl-(-)-epigallocatechins possessing a straight-, branched-, phenyl-inserted- or 1,4-phenylene-inserted-acyl chain of varying length from C4 to C18 were synthesized and evaluated their inhibitory effects against the activation of the Epstein-Barr virus early antigen (EBV-EA). It was indicated that the epigallocatechin derivatives having the straight- or branched-acyl chain of C8 to C11 carbon atoms achieve marked effects.     

3D-QSAR study indicates an enhancing effect of membrane ions on psychiatric drugs targeting serotonin receptor 5-HT1A

Antidepressants and antipsychotics are psychiatric agents used for the treatment of various types of psychiatric diseases. Although currently among the most commonly prescribed drugs, their effectiveness and adverse effects are the topic of many studies and controversial claims. Here we generate QSAR models based on compounds series including 20 drugs recommended for two critical psychiatric diseases: depression and schizophrenia and we use these QSAR models to predict the biological activity of these 20 antidepressants and antipsychotics. We establish the membrane ions' contributions (sodium, potassium, calcium and iron) mediated by water to the antagonism of these drugs at the 5-HT1A receptor. The reliability of our QSAR models in predicting compounds activity is indicated by significant values for cross-validated correlation q2 (0.60-0.76) and fitted correlation r2 (0.96-0.98) coefficients. Our results indicate that potassium, calcium and iron play a key role for the antagonistic activity of drugs at the 5-HT1A receptor. Moreover, based on the established QSAR equations, we analysed 24 new escitalopram derivatives as possibly improved antidepressants targeting the 5-HT1A receptor. We identified that the presence of methyl groups and hydrogen atoms improves antidepressant activity while the simultaneous presence of ethyl, propyl or halogens decreased drastically antidepressant activity at the 5-HT1A site.     

Anharmonic Normal Mode Analysis of Elastic Network Model Improves the Modeling of Atomic Fluctuations in Protein Crystal Structures

Protein conformational dynamics, despite its significant anharmonicity, has been widely explored by normal mode analysis (NMA) based on atomic or coarse-grained potential functions. To account for the anharmonic aspects of protein dynamics, this study proposes, and has performed, an anharmonic NMA (ANMA) based on the C_α-only elastic network models, which assume elastic interactions between pairs of residues whose C_α atoms or heavy atoms are within a cutoff distance. The key step of ANMA is to sample an anharmonic potential function along the directions of eigenvectors of the lowest normal modes to determine the mean-squared fluctuations along these directions. ANMA was evaluated based on the modeling of anisotropic displacement parameters (ADPs) from a list of 83 high-resolution protein crystal structures. Significant improvement was found in the modeling of ADPs by ANMA compared with standard NMA. Further improvement in the modeling of ADPs is attained if the interactions between a protein and its crystalline environment are taken into account. In addition, this study has determined the optimal cutoff distances for ADP modeling based on elastic network models, and these agree well with the peaks of the statistical distributions of distances between C_α atoms or heavy atoms derived from a large set of protein crystal structures.     

1H and 31P Fourier transform magnetic resonance studies of the conformation of enzyme-bound propionyl coenzyme A on the transcarboxylase.

The relaxation rates of the carbon-bound protons and of the three assigned phosphorus resonances of propionyl-CoA were measured in solutions of free propionyl-CoA and of the transcarboxylase-propionyl-CoA complex. In free propionyl-CoA, analysis of the 1/T1 values of 15 protons at 100 and 220 MHz and of 1/T1 and 1/T2 of the three phosphorus atoms at 40.5 MHz indicated free rotation of the propionyl region (taur approximately 3 x 10(-11) sec) but hindered motion of the remainder of the molecule with correlation times of 1-3. 5 x 10(-10) sec, approaching the tumbling time of the entire molecule (taur - 6 x 10(-10) sec. The correlation times of the three phosphorus atoms were indistinguishable from those of their nearest neighbor protons. The effects of three homogeneous enzyme preparations with varying contents of Zn(II), Co(II), and Cu(II) on 1/T1 of 12 protons and 3 phosphorus atoms of prionyl-CoA were analyzed with the help of simultaneous equations to yield the individual contributions at the three metal sites. Only diamagnetic effects were detected on the relaxation rates of the three phosphorus atoms. From the diamagnetic effects it was calculated that the motions of the prionyl side chain and of the terminal pantetheine methylene protons were hindered on the enzyme by an order of magnitude (taur approximately 6 x 10(-10) sec) and that the phosphorus atoms were hindered by two orders of magnitude (taur approximately 1 x 10(-8) sec) over the taur values found in free propionyl-CoA, but that these taur values remained well below that of the entire protein molecule (taur =6 x 10(-7) sec)...     

Selecting an appropriate genetic evaluation model for selection in a developing dairy sector

This study aimed to identify genetic evaluation models (GEM) to accurately select cattle for milk production when only limited data are available. It is based on a data set from the Pakistani Sahiwal progeny testing programme which includes records from five government herds, each consisting of 100 to 350 animals, with lactation records dating back to 1968. Different types of GEM were compared, namely: (1) multivariate v. repeatability model when using the first three lactations, (2) an animal v. a sire model, (3) different fixed effects models to account for effects such as herd, year and season; and (4) fitting a model with genetic parameters fixed v. estimating the genetic parameters as part of the model fitting process. Two methods were used for the comparison of models. The first method used simulated data based on the Pakistani progeny testing system and compared estimated breeding values with true breeding values. The second method used cross-validation to determine the best model in subsets of actual Australian herd-recorded data. Subsets were chosen to reflect the Pakistani data in terms of herd size and number of herds. Based on the simulation and the cross-validation method, the multivariate animal model using fixed genetic parameters was generally the superior GEM, but problems arise in determining suitable values for fixing the parameters. Using mean square error of prediction, the best fixed effects structure could not be conclusively determined. The simulation method indicated the simplest fixed effects structure to be superior whereas in contrast, the cross-validation method on actual data concluded that the most complex one was the best. In conclusion it is difficult to propose a universally best GEM that can be used in any data set of this size. However, some general recommendations are that it is more appropriate to estimate the genetic parameters when evaluating for selection purposes, the animal model was superior to the sire model and that in the Pakistani situation the repeatability model is more suitable than a multivariate.     

Isotopic effects on retention times of caffeine and its metabolites 1,3,7-trimethyluric acid, theophylline, theobromine and paraxanthine

Physicochemical parameters that influence gas chromatographic separation are numerous. Consequently, isotope labelling, because it modifies physicochemical properties, can induce isotopic effects on retention time. Caffeine has been chosen to study this influence because as itself and its metabolites, it allows the preparation of different methylxanthine isotopomers and thus is one of the best models to study isotopic effects induced by stable isotope labelling. Using a caffeine molecule labelled with deuterium at different positions and rat hepatocytes to obtain metabolites, it was possible to study the influence of labelling on retention time [(14% cyanopropylphenyl)methylpolysiloxane] and to point out the role of each labelled site. It appears that isotopic effects induced by the labelling depend not only on the number of labelling atoms but also on whether this labelling is at position 1, 3 or 7 and, consequently, on the role of the labelled site on the function of the molecule.     

Harmonic force field for nitro compounds

Molecular simulations leading to sensors for the detection of explosive compounds require force field parameters that can reproduce the mechanical and vibrational properties of energetic materials. We developed precise harmonic force fields for alanine polypeptides and glycine oligopeptides using the FUERZA procedure that uses the Hessian tensor (obtained from ab initio calculations) to calculate precise parameters. In this work, we used the same procedure to calculate generalized force field parameters of several nitro compounds. We found a linear relationship between force constant and bond distance. The average angle in the nitro compounds was 116°, excluding the 90° angle of the carbon atoms in the octanitrocubane. The calculated parameters permitted the accurate molecular modeling of nitro compounds containing many functional groups. Results were acceptable when compared with others obtained using methods that are specific for one type of molecule, and much better than others obtained using methods that are too general (these ignore the chemical effects of surrounding atoms on the bonding and therefore the bond strength, which affects the mechanical and vibrational properties of the whole molecule).     

Carbon doped boron phosphide nanotubes: A computational study

A computational study based on density functional theory (DFT) calculations has been performed to investigate the properties of the electronic structure of carbon doped boron phosphide nanotube (C-doped BPNT). Pristine and the C-doped structures of two representative (6,0) zigzag and (4,4) armchair BPNTs have been investigated. At first, the geometries of the structures have been allowed to relax by optimization. Subsequently, NMR parameters have been calculated in the optimized structures. The results indicated that the influence of C-doping was more significant on the geometries of the zigzag model than the armchair one. The difference of band gap energies between the pristine and C-doped armchair BPNT was larger than the zigzag model. Significant differences of NMR parameters of those nuclei directly contributed to the C-doping atoms have been observed.