Aconitum and delphinium diterpenoid alkaloids of local anesthetic activity: comparative QSAR analysis based on GA-MLRA/PLS and optimal descriptors approach.

The duration of anesthesia (related to protein binding of a drug) and the onset time (determined by the pKa) are important characteristics in assessment of local anesthetic agents. They are known to be affected by a number of factors. Early studies of antiarrhythmic diterpenoid alkaloids from plants Aconitum and Delphinium suggested that they possess local anesthetic activity due to their ability to suppress sodium currents of excited membranes. In this study we utilized toxicity, duration, and onset of action as endpoints to construct Quantitative Structure-Activity Relationship (QSAR) models for the series of 34 diterpenoid alkaloids characterized by local anesthetic activity using genetic algorithm-based multiple linear regression analysis/partial least squares and simplified molecular input line entry system (SMILES)-based optimal descriptors approach. The developed QSAR models correctly reflected factors that determine three endpoints of interest. Toxicity correlates with descriptors describing partition and reactivity of compounds. The duration of anesthesia was encoded by the parameters defining the ability of a compound to bind at the receptor site. The size and number of H-bond acceptor atoms were found not to favor the speed of onset, while topographic electronic descriptor demonstrated strong positive effect on it. SMILES-based optimal descriptors approach resulted in overall improvement of models. This approach was shown to be more sensitive to structural peculiarities of molecules than regression methods. The results clearly indicate that obtained QSARs are able to provide distinct rationales for compounds optimization with respect to particular endpoint.

CORAL: QSPR modeling of rate constants of reactions between organic aromatic pollutants and hydroxyl radical.

The rate constants (K(OH)) of reactions between 78 organic aromatic pollutants and hydroxyl radical were examined. Simplified molecular input line entry system was used as representation of the molecular structure of the pollutants. Quantitative structure-property relationships was developed using CORAL software (http://www.insilico.eu/CORAL) for four random splits of the data into the subtraining, calibration, and test sets. The obtained results reveal good predictive potential of the applied approach: correlation coefficients (r(2)) for the test sets of the four random splits are 0.75, 0.91, 0.84, and 0.80. Using the Monte Carlo method CORAL software generated the optimal descriptors for one-variable models. The reproducibility of each model was tested performing three runs of the Monte Carlo optimization. The current data were compared to previous results and discussed.

Predictions of Gibbs free energies for the adsorption of polyaromatic and nitroaromatic environmental contaminants on carbonaceous materials: efficient computational approach.

The adsorption of benzene, polycyclic aromatic hydrocarbons (PAHs), and nitroaromatic compounds (NACs) on the carbonaceous surfaces from the gas phase and water solution was investigated. Several different levels of theory were applied, including DFT-, MP2-, and CCSD(T)-based methods, to find an approach that is computationally inexpensive and can provide accurate thermodynamic parameters for studied adsorption phenomena. The methods and techniques used (including cluster and periodic approximations) were evaluated on the basis of comparison with available experimental data. The optimized structures of calculated complexes are obtained, and the interaction energies and Gibbs free energies are predicted. Good agreement was revealed for the theoretical and experimental adsorption energies of benzene and PAHs adsorbed on the carbon surfaces. The adsorption of benzene, PAHs, and NACs on carbon is suggested to be effective from the gas phase for all studied compounds and for PAHs and NACs also from water solution at room temperature.

Identification of potential antivirals against 3CLpro enzyme for the treatment of SARS-CoV-2: A multi-step virtual screening study.

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) outbreak is posing a serious public health threat worldwide in the form of COVD-19. Herein, we have performed two-dimensional quantitative structure-activity relationship (2D-QSAR) and three-dimensional pharmacophore modelling analysis employing inhibitors of 3-chymotrypsin-like protease (3CLpro), the leading protease that is crucial for the replication of SARS-CoV-2. The investigation aims to identify the important structural features responsible for the enzyme inhibition and the search for novel 3CLpro enzyme inhibitors as effective therapeutics for treating SARS-CoV-2. Furthermore, we carried out molecular docking studies using the most and least active compounds in the dataset, aiming to validate the contributions of various features as appeared in the QSAR models. Later, the stringently validated 2D-QSAR model was used to estimate the 3CLpro inhibitory activity of compounds from five chemical databases. Compounds with the significant predicted activity were then subjected to pharmacophore-based virtual screening to screen the top-rated compounds, which were then further subjected to molecular docking analysis, absorption, distribution, metabolism, excretion - toxicity (ADMET) profiling, and molecular dynamics (MD) simulation. The multi-step virtual screening analyses suggested that compounds CASAntiV-865453-58-3, CASAntiV-865453-40-3, and CASAntiV-2043031-84-9 could be used as effective therapeutic agents for the treatment of SARS-CoV-2.

CORAL: quantitative structure-activity relationship models for estimating toxicity of organic compounds in rats.

For six random splits, one-variable models of rat toxicity (minus decimal logarithm of the 50% lethal dose [pLD50], oral exposure) have been calculated with CORAL software (http://www.insilico.eu/coral/). The total number of considered compounds is 689. New additional global attributes of the simplified molecular input line entry system (SMILES) have been examined for improvement of the optimal SMILES-based descriptors. These global SMILES attributes are representing the presence of some chemical elements and different kinds of chemical bonds (double, triple, and stereochemical). The "classic" scheme of building up quantitative structure-property/activity relationships and the balance of correlations (BC) with the ideal slopes were compared. For all six random splits, best prediction takes place if the aforementioned BC along with the global SMILES attributes are included in the modeling process. The average statistical characteristics for the external test set are the following: n = 119 ± 6.4, R(2) = 0.7371 ± 0.013, and root mean square error = 0.360 ± 0.037.

Adsorption of thymine and uracil on 1:1 clay mineral surfaces: comprehensive ab initio study on influence of sodium cation and water.

This computational study performed using the density functional theory shows that hydrated and non-hydrated tetrahedral and octahedral kaolinite mineral surfaces in the presence of a cation adsorb the nucleic acid bases thymine and uracil well. Differences in the structure and chemistry of specific clay mineral surfaces led to a variety of DNA bases adsorption mechanisms. The energetically most predisposed positions for an adsorbate molecule on the mineral surface were revealed. The target molecule binding with the surface can be characterized as physisorption, which occurs mainly due to a cation-molecular oxygen interaction, with hydrogen bonds providing an additional stabilization. The adsorption strength is proportional to the number of intermolecular interactions formed between the target molecule and the surface. From the Atoms in Molecules analysis and comparison of binding energy values of studied systems it is concluded that the sorption activity of kaolinite minerals for thymine and uracil depends on various factors, among which are the structure and accessibility of the organic compounds. The adsorption is governed mostly by the surface type, its properties and presence of cation, which cause a selective binding of the nucleobase. Adsorbate stabilization on the mineral surface increases only slightly with explicit addition of water. Comparison of activity of different studied kaolinite mineral models reveals the following order for stabilization: octahedral-Na-water > octahedral-Na > tetrahedral-Na > tetrahedral-Na-water. Further investigation of the electrostatic potentials helps understanding of the adsorption process and confirmation of the active sites on the kaolinite mineral surfaces. Based on the conclusions that clay mineral affinity for DNA and RNA bases can vary due to different structural and chemical properties of the surface, a hypothesis on possible role of clays in the origin of life was made.

Can the Gibbs free energy of adsorption be predicted efficiently and accurately: an M05-2X DFT study.

This study presents new insight into the prediction of partitioning of organic compounds between a carbon surface (soot) and water, and it also sheds light on the sluggish desorption of interacting molecules from activated and nonactivated carbon surfaces. This paper provides details about the structure and interactions of benzene, polycyclic aromatic hydrocarbons, and aromatic nitrocompounds with a carbon surface modeled by coronene using a density functional theory approach along with the M05-2X functional. The adsorption was studied in vacuum and from water solution. The molecules studied are physisorbed on the carbon surface. While the intermolecular interactions of benzene and hydrocarbons are governed by dispersion forces, nitrocompounds are adsorbed also due to quite strong electrostatic interactions with all types of carbon surfaces. On the basis of these results, we conclude that the method of prediction presented in this study allows one to approach the experimental level of accuracy in predicting thermodynamic parameters of adsorption on a carbon surface from the gas phase. The empirical modification of the polarized continuum model leads also to a quantitative agreement with the experimental data for the Gibbs free energy values of the adsorption from water solution.

Abdominal aortic Doppler waveform in patients with aorto-iliac disease.

OBJECTIVES: The mid-systolic deceleration (notch) in the proximal descending aortic Doppler waveform was reported to be common in patients with aorto-iliac disease. However, evaluation of the descending aorta is limited to echocardiography and may be technically difficult. Therefore, we decided to check whether similar Doppler flow disturbance can be found in abdominal aorta, which is easily evaluated in wider range of patients undergoing general abdominal and vascular ultrasound, as well as echocardiography. METHODS: We evaluated 115 consecutive symptomatic patients with severe peripheral artery disease admitted for vascular surgery, and 60 controls. The presence or absence of the mid-systolic deceleration in the Doppler waveform was evaluated retrospectively, by the single echocardiographer blinded to the localisation of the arterial occlusion or stenosis. RESULTS: The mid-systolic notch in the proximal abdominal aorta was present in 58 of 71 patients (82%) with significant aorto-iliac disease, seven of 44 (16%) patients with occlusion or significant stenosis distally to the external iliac artery (P < 0.001) and in none of the patients from the control group. Sensitivity, specificity and positive predictive value of the mid-systolic notch in the abdominal aortic Doppler waveform in the detection of aorto-iliac disease in patients with peripheral artery disease were 82%, 84% and 89%, respectively. CONCLUSION: The mid-systolic deceleration (notch) in the proximal abdominal Doppler waveform is a simple ultrasonographic marker of significant aorto-iliac disease.

Loops in globular proteins: a novel category of secondary structure.

The protein loop, a novel category of nonregular secondary structure, is a segment of contiguous polypeptide chain that traces a "loop-shaped" path in three-dimensional space; the main chain of an idealized loop resembles a Greek omega (omega). A systematic study was made of 67 proteins of known structure revealing 270 omega loops. Although such loops are typically regarded as "random coil," they are, in fact, highly compact substructures and may also be independent folding units. Loops are almost invariably situated at the protein surface where they are poised to assume important roles in molecular function and biological recognition. They are often observed to be modules of evolutionary exchange and are also natural candidates for bioengineering studies.

Thermodynamics and kinetics of intramolecular water assisted proton transfer in Na+ -1-methylcytosine water complexes.

High-level ab initio predictions of the tautomerization equilibrium and rate constants of water-assisted proton transfer of 1-methyl-cytosine (MeC) to its MeC* imino tautomeric form in the presence of up to two water molecules (W) and the Na(+) cation were carried out. The calculated energy values were used to obtain the thermodynamic parameters and equilibrium concentration of MeC, its rare tautomer, and their complexes with up to two water molecules and the Na (+) cation. The rate constants for the tautomerization were obtained by using the instanton approach (a quasiclassical method based on the least-action principle). Hydration of MeC by one water molecule leads to an increase of the concentration of the MeC* tautomer in the equilibrium mixture and a decrease of the barrier of the MeC* formation (to 15.6 kcal/mol). If the Na(+) cation is present, the tautomeric form is much less favored, and the tautomerization barrier increases to 25.2 kcal/mol. It was found that MeC monohydrate has both the highest equilibrium (2.9 x 10(-2)) and rate (7.9 x 10(5) s(-1)) constants of tautomerization in comparison to the MeC*NaW and MeC*Na2W complexes containing the Na(+) cation. Moreover, this study also allows one to estimate the concentration of MeC present in the cell during DNA synthesis as the unwanted tautomer, which in forming a mismatched base pair can cause spontaneous point mutations. Kinetic simulations have demonstrated that the low values of equilibrium (10(-14)-10(-13)) and rate constants (10(-17)-10(-16) s(-1)) of tautomerization make contribution of the MeC*Na(+)W and MeC*Na(+)2W complexes to the point mutations in DNA unlikely. In contrast to these Na(+) complexes, MeC*W can reach an equilibrium concentration of 2.9 x 10(-2) within 10(-7) s; thus, it can increase the probability of the point mutations.

Influence of the central metal ion on nonlinear optical and two-photon absorption properties of push-pull transition metal porphyrins.

We present a quantum-chemical analysis of the central metal ion's effect on first hyperpolarizabilities and two-photon absorption (TPA) cross sections at the infrared region of a series of push-pull porphyrins whose synthesis and NLO properties have been reported earlier (J. Am. Chem. Soc. 2005, 127, 9710). The molecular geometries are obtained via the B3LYP/6-31G(d,p) level optimization including SCRF/PCM approach, and the NLO and TPA properties are calculated with the ZINDO/CV method including solvent effects. It is found that the CT transition between the metal ion's d orbital and the macrocycle pi orbitals plays an important role on NLO and TPA properties of metal porphyrins. Our data suggest a new approach to enhance TPA properties of porphyrin materials. We also present a quantum-chemical analysis on porphyrin dimers and trimers to understand the relationship between structural and collective NLO properties. It has been observed that beta values can be improved about an order of magnitude and TPA properties can be enhanced by 2 orders of magnitude by the formation of a trimer. The importance of our results with respect to the design of photonic and photodynamic therapy materials have been discussed.

The effect of nitroaromatics' composition on their toxicity in vivo: novel, efficient non-additive 1D QSAR analysis.

Novel 1D QSAR approach that allows analysis of non-additive effects of molecular fragments on toxicity has been proposed. Twenty-eight nitroaromatic compounds including some well-known explosives have been chosen for this study. The 50% lethal dose concentration for rats (LD50) was used as the estimation of toxicity in vivo to develop 1D QSAR models on the framework of Simplex representation of molecular structure. The results of 1D QSAR analysis show that even the information about the composition of molecules provides the main trends of toxicity changes. The necessity of consideration of substituents' mutual impacts for the development of adequate QSAR models of nitroaromatics' toxicity was demonstrated. Statistic characteristics for all the developed partial least squares QSAR models, except the additive ones are quite satisfactory (R2=0.81-0.92; Q2=0.64-0.83; R2 test=0.84-0.87). A successful performance of such models is due to their non-additivity i.e. possibility of taking into account the mutual influence of substituents in benzene ring which plays the governing role for toxicity change and could be mediated through the different C-H fragments of the ring. The correspondence between observed and predicted by these models toxicity values is good. This allowing combine advantages of such approaches and develop adequate consensus model that can be used as a toxicity virtual screening tool.

Molecular modelling and QSAR analysis of the estrogenic activity of terpenoids isolated from Ferula plants.

The relationship between chemical structure and estrogenic activity in a series of terpenoid esters with aromatic and aliphatic acid substituents isolated from Ferula plants, was studied. The fragments of the terpenoid structure that are potentially responsible for estrogenic activity were revealed. A quantitative structure-estrogenic activity study has been carried out using the QSAR approach with use of data derived from quantum-chemical calculations as well as data generated from three-dimensional structures of terpenoids. A number of molecular descriptors was obtained from the density functional theory (DFT) at the B3LYP/6-31G(d, p) level of calculation. Comparative analysis of the quantum-chemical computational data was also performed to confirm hypothesis concerning importance of the distance between the oxygen of alcohol hydroxyl group and the functional group in the para-position of the benzene ring (the hydroxyl or methoxy group). Use of the Genetic Algorithm in the QSAR analysis allowed the structural and physicochemical parameters of the terpenoids responsible for estrogenic activity to be determined. A significant QSAR model was obtained with an r(2) value of 0.892. The resulting model showed a reliable dependence of estrogenic activity of the terpenoids on such parameters as molecular shape, number of phenolic groups, surface polarity and the energy of the highest occupied molecular orbital.

A quantum chemical based toxicity study of estimated reduction potential and hydrophobicity in series of nitroaromatic compounds.

Nitroaromatic compounds and the products of their degradation are toxic to bacteria, cells and animals. Various studies have been carried out to better understand the mechanism of toxicity of aromatic nitrocompounds and their relationship to humans and the environment. Recent data relate cytotoxicity of nitroaromatic compounds to their single- or two-electron enzymatic reduction. However, mechanisms of animal toxicity could be more complex. This work investigates the estimated reduction and oxidation potentials of 34 nitroaromatic compounds using quantum chemical approaches. All geometries were optimized with density functional theory (DFT) using the solvation model based on density (SMD) and polarizable continuum model (PCM) solvent model protocols. Quantitative structure-activity/property (QSAR/QSPR) models were developed using descriptors obtained from quantum chemical optimizations as well as the DRAGON software program. The QSAR/QSPR equations developed consist of two to four descriptors. Correlations have been identified between electron affinity (ELUMO) and hydrophobicity (log P).

In silico kinetics of alkaline hydrolysis of 1,3,5-trinitro-1,3,5-triazinane (RDX): M06-2X investigation.

Alkaline hydrolysis of RDX (hexahydro-1,3,5-trinitro-1,3,5-triazine), as one of the most promising methods for nitrocompound remediation, was investigated computationally at the PCM(Pauling)/M06-2X/6-311++G(d,p) level of theory. Computational simulation shows that RDX hydrolysis is a highly exothermic multistep process involving initial deprotonation and nitrite elimination, cycle cleavage, further transformation of cycle-opened intermediates to end products caused by a series of C-N bond ruptures, hydroxide attachments, and proton transfers. Computationally predicted products of RDX hydrolysis such as nitrite, nitrous oxide, formaldehyde, formate, and ammonia correspond to experimentally observed ones. Accounting of specific hydration of hydroxide is critical to create an accurate kinetic model for alkaline hydrolysis. Simulated kinetics of the hydrolysis are in good agreement with available experimental data. A period of one month is necessary for 99% RDX decomposition at pH 10. Computations predict significant increases of the reaction rate of hydrolysis at pH 11, pH 12, and pH 13.

Theoretical study of adsorption of sarin and soman on tetrahedral edge clay mineral fragments.

This study provides details of the structure and interactions of Sarin and Soman with edge tetrahedral fragments of clay minerals. The adsorption mechanism of Sarin and Soman on these mineral fragments containing the Si(4+) and Al(3+) central cations was investigated. The calculations were performed using the B3LYP and MP2 levels of theory in conjunction with the 6-31G(d) basis set. The studied systems were fully optimized. Optimized geometries, adsorption energies, and Gibbs free energies of Sarin and Soman adsorption complexes were computed. The number and strength of formed intermolecular interactions have been analyzed using the AIM theory. The charge of the systems and a termination of the mineral fragment are the main contributing factors on the formation of intermolecular interactions in the studied systems. In the neutral complexes, Sarin and Soman is physisorbed on these mineral fragments due to the formation of C-H...O, and O-H...O hydrogen bonds. The chemical bond is formed between a phosphorus atom of Sarin and Soman and an oxygen atom of the -2 charged clusters containing an Al(3+) central cation and -1 charged complex containing a Si(4+) central cation (chemisorption). Sarin and Soman interact mostly in the same way with the same terminated edge mineral fragments containing different central cations. However, the interaction energies of the complexes with an Al(3+) central cation are larger than these values for the Si(4+) complexes. The interaction enthalpies of all studied systems corrected for the basis set superposition error were found to be negative. However, on the basis of the Gibbs free energy values, only strongly interacting complexes containing a charged edge mineral fragment with an Al(3+) central cation are stable at room temperature. We can conclude that Sarin and Soman will be adsorbed preferably on this type of edge mineral surfaces. Moreover, on the basis of the character of these edge surfaces, a tetrahedral edge mineral fragment can provide effective centers for the dissociation.

Prediction of CL-20 chemical degradation pathways, theoretical and experimental evidence for dependence on competing modes of reaction.

Highest occupied and lowest unoccupied molecular orbital energies, formation energies, bond lengths and FTIR spectra all suggest competing CL-20 degradation mechanisms. This second of two studies investigates recalcitrant, toxic, aromatic CL-20 intermediates that absorb from 370 to 430 nm. Our earlier study (Struct. Chem., 15, 2004) revealed that these intermediates were formed at high OH(-) concentrations via the chemically preferred pathway of breaking the C-C bond between the two cyclopentanes, thereby eliminating nitro groups, forming conjugated pi bonds, and resulting in a pyrazine three-ring aromatic intermediate. In attempting to find and make dominant a more benign CL-20 transformation pathway, this current research validates hydroxylation results from both studies and examines CL-20 transformations via photo-induced free radical reactions. This article discusses CL-20 competing modes of degradation revealed through: computational calculation; UV/VIS and SF spectroscopy following alkaline hydrolysis; and photochemical irradiation to degrade CL-20 and its byproducts at their respective wavelengths of maximum absorption.

Reliable but Timesaving: In Search of an Efficient Quantum-chemical Method for the Description of Functional Fullerenes.

Fullerene and its derivatives are currently one of the most intensively investigated species in the area of nanomedicine and nanochemistry. Various unique properties of fullerenes are responsible for their wide range applications in industry, biology and medicine. A large pool of functionalized C60 and C70 fullerenes is investigated theoretically at different levels of quantum-mechanical theory. The semiempirial PM6 method, density functional theory with the B3LYP functional, and correlated ab initio MP2 method are employed to compute the optimized structures, and an array of properties for the considered species. In addition to the calculations for isolated molecules, the results of solution calculations are also reported at the DFT level, using the polarizable continuum model (PCM). Ionization potentials (IPs) and electron affinities (EAs) are computed by means of Koopmans' theorem as well as with the more accurate but computationally expensive ΔSCF method. Both procedures yield comparable values, while comparison of IPs and EAs computed with different quantum-mechanical methods shows surprisingly large differences. Harmonic vibrational frequencies are computed at the PM6 and B3LYP levels of theory and compared with each other. A possible application of the frequencies as 3D descriptors in the EVA (EigenVAlues) method is shown. All the computed data are made available, and may be used to replace experimental data in routine applications where large amounts of data are required, e.g. in structure-activity relationship studies of the toxicity of fullerene derivatives.

Ab Initio Studies of Anatase TiO2 (101) Surface-supported Au8 Clusters.

Supported transition metals on TiO2 surfaces have shown exceptional catalytic properties in many important process such as CO oxidation, selective propane oxidation, hydrogenation, water adsorption and other catalytic and photocatalytic oxidation reaction at low-temperature. Among the three polymorphs of TiO2, the anatase crystal is the more photoactive. The anatase (101) surface attracts more attention since it has lower surface energy relative to (001) and (100) surfaces and it is observed to adsorb small molecules on its surface. Using density-functional theory (DFT) with on-site Coulomb interactions corrections, we have computed the structural and electronic properties of selected Au8 clusters interacting with clean and reduced anatase TiO2(101) surfaces. The computed adsorption energies are suggesting that the considered Au8 clusters are only physisorbed onto pristine TiO2(101) surface. Oxygen vacancies are found to enhance the absorption of Au8 on the Ti2(101) surface. Accurate simulations required spin polarized DFT since the ground state of Au8 interacting with defective TiO2(101) shows magnetic solutions. The results show that Au8 clusters are chemically bonded to the surface around the locality of the oxygen vacancy. The surface oxygen vacancy is found to be energetically more favourable than sub-surface oxygen vacancy configuration. These vacancy sites may act as nucleation sites for small Au clusters or Au atoms. Finally, the computed electronic structure of all the Au8/TiO2(101) configurations considered in this work are analysed in the light of available experimental data.

The molecular structure and electronic spectrum of the C@C60 endohedral complex: an ab initio study.

We performed ab initio quantum-chemical studies at the HF and B3LYP levels of theory to determine the geometry and electronic spectra of the C@C60 complex. The STO-3G and 6-31G(d) basis sets were employed. Two different types of stable conformations for the endohedral atom related to the centers of five- and six-member rings were found. The estimated potential barrier between those conformations is small; under certain conditions, the endohedral atom can transfer from one location in the molecular cage to another. The influence of the location of endohedral carbon atom on the electronic spectrum of C60 is discussed.